数学电工数模学习教案

1、数学数学(shxu)电工数模电工数模第一页，共120页。第1页/共120页第二页，共120页。第2页/共120页第三页，共120页。第3页/共120页第四页，共120页。第4页/共120页第五页，共120页。第5页/共120页第六页，共120页。第6页/共120页第七页，共120页。ModelingModeling）应用知识从实际课题中抽象、应用知识从实际课题中抽象、提炼出数学模型的过程。提炼出数学模型的过程。1.1 数学数学(shxu)模型与数学模型与数学(shxu)建模建模 第7页/共120页第八页，共120页。 1.了解问题的实际背景，明确建模目的，收集掌握必要的数据资料。 2. 通过(

2、tnggu)对资料的分析计 算， 找出起主要作用的因素，经必要的精炼、简化，提出若干符合客观实际的假设。 3.在所作假设的基础上，利用适当的数学工具去刻划各变量之间的关系，建立相应的数学结构 即建立数学模型。 4.模型求解。 5.模型的分析与检验。 在难以得出解析解时，也应当借助在难以得出解析解时，也应当借助 计算机计算机 求出数值解。求出数值解。 实体信实体信息息(数据数据)假设假设建模建模求解求解验证验证应用应用第8页/共120页第九页，共120页。对某个实际问题对某个实际问题了解的深入程度了解的深入程度白箱模型、灰箱模型、黑箱模型白箱模型、灰箱模型、黑箱模型模型中变量的特模型中变量的特征

3、征连续型模型、离散型模型或确定性连续型模型、离散型模型或确定性模型、随机型模型等模型、随机型模型等建模中所用的数建模中所用的数学方法学方法初等模型、微分方程模型、差分方初等模型、微分方程模型、差分方程模型、优化模型等程模型、优化模型等研究课题的实际研究课题的实际范畴范畴人口模型、生态系统模型人口模型、生态系统模型 、交通、交通流模型、经济模型、流模型、经济模型、 基因模型等基因模型等第9页/共120页第十页，共120页。数学建模实践的数学建模实践的 每一步中都每一步中都 蕴含着能力上的蕴含着能力上的 锻炼，在调锻炼，在调查研究阶段，需查研究阶段，需 要用到观察能力、分析能力和数据处理能力要用到

4、观察能力、分析能力和数据处理能力等。在提出假设等。在提出假设 时，又需要用到时，又需要用到 想象力和归纳想象力和归纳 简化能力。简化能力。 在真正开始自己的研究之前，还应当尽可能先了解一下前在真正开始自己的研究之前，还应当尽可能先了解一下前人或别人的工作，使自己的工人或别人的工作，使自己的工 作成为作成为(chngwi)(chngwi)别人研究别人研究工作工作 的继续而不是别人工作的重复，你可以把某些已知的研的继续而不是别人工作的重复，你可以把某些已知的研究结果用作你的假设，去探索新的奥秘。因此我们还应当学究结果用作你的假设，去探索新的奥秘。因此我们还应当学会在尽可能短的时间会在尽可能短的时间

5、 内查到并学会我想应用的知识的本领。内查到并学会我想应用的知识的本领。 还需要你多少要有点还需要你多少要有点 创新的能力。这种能力不是生来就有创新的能力。这种能力不是生来就有的，建模实践就为你提供了一个培养创新能力的机会。的，建模实践就为你提供了一个培养创新能力的机会。 开设数学建模课的主要目的为了提高学开设数学建模课的主要目的为了提高学 生的生的综合素质综合素质，增强，增强 应用数学知识应用数学知识 解决实际问解决实际问 题的本领。题的本领。第10页/共120页第十一页，共120页。例例1 某人平时下班总是按预定时间某人平时下班总是按预定时间(shjin)到达某处，然到达某处，然然后他妻子开

6、车接他回家。有一天，他比平时提早然后他妻子开车接他回家。有一天，他比平时提早了三十分钟到达该处，于是此人就沿着妻子来接他了三十分钟到达该处，于是此人就沿着妻子来接他的方向步行回去并在途中遇到了妻子，这一天，他的方向步行回去并在途中遇到了妻子，这一天，他比平时提前了十分钟到家，问此人共步行了多长时比平时提前了十分钟到家，问此人共步行了多长时间？间？ 似乎条件似乎条件(tiojin)不够哦不够哦 。 换一种想法，问题就迎刃而解了。换一种想法，问题就迎刃而解了。假如他的妻子遇到他后仍载着他假如他的妻子遇到他后仍载着他开往会合地点，那么这一天他就开往会合地点，那么这一天他就不会提前回家了。提前的十分钟

7、不会提前回家了。提前的十分钟时间从何而来？时间从何而来？ 显然是由于节省了从相遇点到会合点，又从会合点返回相遇点这一段路的缘故，故由相遇点到会合点需开显然是由于节省了从相遇点到会合点，又从会合点返回相遇点这一段路的缘故，故由相遇点到会合点需开5分钟。而此人提前了三十分钟到达会合点，故相遇时他已步行了二十五分钟。分钟。而此人提前了三十分钟到达会合点，故相遇时他已步行了二十五分钟。 第11页/共120页第十二页，共120页。例例2 2 交通灯在绿灯转换成红灯时，有一交通灯在绿灯转换成红灯时，有一个过渡状态个过渡状态亮一段时间的黄灯。请分亮一段时间的黄灯。请分析黄灯应当亮多久。析黄灯应当亮多久。设想

8、一下黄灯的作用是什么，不难看设想一下黄灯的作用是什么，不难看出，黄灯起的是警告的作用，意思是出，黄灯起的是警告的作用，意思是马上要转红灯了，假如你能停住，请马上要转红灯了，假如你能停住，请立即停车。停车是需要时间的，在这立即停车。停车是需要时间的，在这段时间内，车辆仍将向前行驶一段距段时间内，车辆仍将向前行驶一段距离离 L。这就是说，在离街口距离为。这就是说，在离街口距离为 L处存在着一条停车线（尽管它没被画处存在着一条停车线（尽管它没被画在地上），见图。对于那些黄灯亮时在地上），见图。对于那些黄灯亮时已过线的车辆，则应当保证它们仍能已过线的车辆，则应当保证它们仍能穿过马路。穿过马路。 马路的

9、宽度马路的宽度 D是容易测得是容易测得 的，问题的关键在的，问题的关键在 于于L的确定。为确定的确定。为确定 L，还应当将，还应当将 L划分为两段：划分为两段：L1和和L2，其中其中 L1是司机在发现黄灯亮及判断应是司机在发现黄灯亮及判断应当刹车的反应时间内驶过的路程当刹车的反应时间内驶过的路程 ，L2为刹车制为刹车制动后车辆驶过的路程。动后车辆驶过的路程。L1较容易计算，交通部门较容易计算，交通部门对司机的平均反应时间对司机的平均反应时间 t1早有测算，反应时间早有测算，反应时间过长将考不出驾照），而此街道的行驶速度过长将考不出驾照），而此街道的行驶速度 v 也是交管部门早已定好的，目的是使

10、交通流量最也是交管部门早已定好的，目的是使交通流量最大，可另建模型研究，从而大，可另建模型研究，从而 L1=v*t1。刹车距。刹车距离离 L2既可用曲线拟合方法得出，也可利用牛顿既可用曲线拟合方法得出，也可利用牛顿第二定律计算出来第二定律计算出来 。黄灯究竟应当亮多久现在已经变得清楚多了。第黄灯究竟应当亮多久现在已经变得清楚多了。第一步，先计算出一步，先计算出 L应多大才能使看见黄灯的司机应多大才能使看见黄灯的司机停得住车。第二步，黄灯亮的时间应当让已过线停得住车。第二步，黄灯亮的时间应当让已过线的车顺利穿过马路，即的车顺利穿过马路，即T 至少应当达到至少应当达到 （L+D）/v。 DL第12

11、页/共120页第十三页，共120页。第13页/共120页第十四页，共120页。某航空母舰派其护卫舰去搜寻其跳伞的飞行某航空母舰派其护卫舰去搜寻其跳伞的飞行员，护卫舰找到飞行员后，航母通知员，护卫舰找到飞行员后，航母通知(tngzh)它尽快它尽快 返回与其汇合并通报了航母当前的航速与方返回与其汇合并通报了航母当前的航速与方 向，问护卫舰应怎样航行，才能与航母汇合。向，问护卫舰应怎样航行，才能与航母汇合。2.1 舰舰艇艇(jintng)的会合的会合第14页/共120页第十五页，共120页。12,11222aabrbaah令：令：则上式可简记成则上式可简记成 ： 222rh-yx)(A(0,b)XY

12、B(0,-b)P(x,y)O航母航母 护卫舰护卫舰 1 2 )()(22222b-yx a byx即：即：22222222) 1(411ababaayx可化为：可化为：记记v2/ v1=a通常通常a1 222|AP|a|BP|则则汇合点汇合点 p必位于此圆上。必位于此圆上。 bxy)(tan1（护卫舰的路线方程）（护卫舰的路线方程）bxy)(tan2（航母的路线方程（航母的路线方程 ）即可求出即可求出P点的坐标点的坐标(zubio)和和2 的值。的值。本模型虽简单，但分析本模型虽简单，但分析极清晰且易于实际应用极清晰且易于实际应用 第15页/共120页第十六页，共120页。在寒冷的北方，在寒冷

13、的北方， 许多住房的玻璃窗都是双层许多住房的玻璃窗都是双层玻璃的，现在我们来建立一个简单玻璃的，现在我们来建立一个简单 的数学模的数学模型，研究一下双层玻璃到底有多大的功效。型，研究一下双层玻璃到底有多大的功效。比较比较(bjio)两座其他条件完全相同的房屋，它们的两座其他条件完全相同的房屋，它们的差异仅仅在窗户不同。差异仅仅在窗户不同。 不妨可以提出以下不妨可以提出以下 假设假设：1、设室内热量的流失是热传导、设室内热量的流失是热传导引起的，不存在户内外的空气对引起的，不存在户内外的空气对流。流。2、室内温、室内温 度度T1与户外温与户外温 度度T2均均为常数。为常数。3、玻璃是均匀的，热传

14、导系数、玻璃是均匀的，热传导系数为常数。为常数。第16页/共120页第十七页，共120页。设玻璃的热传导系数设玻璃的热传导系数 为为k1，空气的热，空气的热传导系数传导系数 为为k2，单位时间通过，单位时间通过(tnggu)单位面积由温度高的一侧流单位面积由温度高的一侧流向温度低的一侧的热量为向温度低的一侧的热量为 ddl室室外外T2室室内内T1TaTb由热传导公式由热传导公式 =kT/d 12121aabbTTTTTTkkkdld)/()(21212121dklkTTdklkTa解得：解得：dklkdTTkddklkTTdklkTk212112121211122)1 (第17页/共120页第

15、十八页，共120页。此函数此函数(hnsh)的图形的图形为为dd室室外外T2室室内内T1dTTk2211)/()(2221dklk类似类似(li s)有有 321621kk一般一般(ybn)dl /811故故记记h=l/d并令并令f(h)= 181h01234567891000.10.20.30.40.50.60.70.80.91hf(h)考虑到考虑到美观美观和使用上和使用上 的的方便方便，h不必取得过大，例如，可不必取得过大，例如，可 取取h=3，即，即l=3d，此时房屋热量的损失不超过单层玻璃窗时，此时房屋热量的损失不超过单层玻璃窗时的的 3% 。 第18页/共120页第十九页，共120页

16、。假如你站在崖顶且身上带着一只具有假如你站在崖顶且身上带着一只具有(jyu)跑表功跑表功 能的计算器，你也许会出于好奇心想用扔下能的计算器，你也许会出于好奇心想用扔下 一块石头听回声的方法来估计山崖的高度，一块石头听回声的方法来估计山崖的高度， 假定你能准确地测定时间，你又怎样来推算假定你能准确地测定时间，你又怎样来推算 山崖的高度呢，请你分析一下这一问题。山崖的高度呢，请你分析一下这一问题。我有一只具有跑我有一只具有跑 表功能的计算器。表功能的计算器。第19页/共120页第二十页，共120页。假定空气阻力不计，可以直接利用自由落体运动的公式假定空气阻力不计，可以直接利用自由落体运动的公式来计

17、算。例如，来计算。例如， 设设t=4秒，秒，g=9.81米米/秒秒2，则可求得，则可求得h78.5米。米。221gth 我学过微积分，我可以我学过微积分，我可以(ky)做做 得更好，呵呵。得更好，呵呵。 第20页/共120页第二十一页，共120页。vKmgdtdvmF除去地球吸引力外，对石块下落影响最大的当除去地球吸引力外，对石块下落影响最大的当 属属空气阻空气阻力力。根据流体力学知识，此时可设空气阻力正比于石块下。根据流体力学知识，此时可设空气阻力正比于石块下落的速度，阻力系落的速度，阻力系 数数K为常数，因而，由牛顿第二定律可为常数，因而，由牛顿第二定律可得：得： kgcevkt令令k=K

18、/m，解得解得 代入初始条件代入初始条件 v(0)=0，得，得c=g/k，故有，故有 ktekgkgv再积分再积分(jfn)一次，得：一次，得： cekgtkghkt2第21页/共120页第二十二页，共120页。若设若设k=0.05并仍设并仍设 t=4秒，则可求秒，则可求 得得h73.6米。米。 听到回声再按跑表听到回声再按跑表(pobio)，计算得到的时间中包含了，计算得到的时间中包含了 反应时间反应时间 不妨设平均反应时间不妨设平均反应时间 为为0.1秒秒 ，假如，假如(jir)仍仍 设设t=4秒，扣除反应秒，扣除反应时间后应时间后应 为为3.9秒，代入秒，代入 式，求得式，求得h69.9

19、米。米。 222)1(kgektkgkgekgtkghktkt多测几次，取平均多测几次，取平均值值代入初始条代入初始条 件件h(0)=0，得到计算山崖高度的公式：，得到计算山崖高度的公式： 将将e-kt用泰勒公式展开并用泰勒公式展开并 令令k 0+ ，即可，即可得出前面不考虑空气阻力时的结果。得出前面不考虑空气阻力时的结果。第22页/共120页第二十三页，共120页。还应考虑还应考虑回声回声传回来所需要的时间。为此，令石块下落传回来所需要的时间。为此，令石块下落 的真正时间的真正时间 为为t1，声音传回来的时间记，声音传回来的时间记 为为t2，还得解一个，还得解一个方程组：方程组： 93340

20、1212211.ttthkg)ekt (kghkt这一方程组是这一方程组是非线性非线性的，求的，求解不太容易，解不太容易，为了估算崖高为了估算崖高竟要去解一个竟要去解一个非线性主程组非线性主程组似乎不合情理似乎不合情理 相对于石块速度，声音速度要快得多，我们可相对于石块速度，声音速度要快得多，我们可 用方法二先求一次用方法二先求一次 h，令，令t2=h/340，校正，校正t，求石，求石块下落时间块下落时间 t1t-t2将将t1代入式代入式再算一次，得出再算一次，得出崖高的近似值。例如，崖高的近似值。例如， 若若h=69.9米，则米，则 t20.21秒，故秒，故 t13.69秒，求得秒，求得 h

21、62.3米。米。 第23页/共120页第二十四页，共120页。 最小二乘法最小二乘法(chngf) 插值方法插值方法 当问题的机理非常不清楚难以直接利用其他当问题的机理非常不清楚难以直接利用其他知识来建模时，一个较为自然的方法是利用知识来建模时，一个较为自然的方法是利用数据进行曲线拟合，找出变量之间的近似依数据进行曲线拟合，找出变量之间的近似依赖关系即函数关系。赖关系即函数关系。2.4 经验模型经验模型第24页/共120页第二十五页，共120页。设经实际测量已得设经实际测量已得 到到n组数据（组数据（xi , yi），），i=1, n。将数。将数据画在平面直角坐标系中，见据画在平面直角坐标系中

22、，见 图。如果建模者判断图。如果建模者判断 这这n个个点很象是分布在某条直线附近，令点很象是分布在某条直线附近，令 该直线方程该直线方程 为为y=ax+b，进而利用数据来求参进而利用数据来求参 数数a和和b。由于该直线只是数据近似满。由于该直线只是数据近似满足足(mnz)的关系式，故的关系式，故 yi-(axi+b)=0一般不成立，但我们一般不成立，但我们希望希望 niiibaxy12)(最小最小此式对此式对a和和b的偏导数的偏导数(do sh)均均 为为0，解相应方程组，求得：，解相应方程组，求得： xaybxxyyxxaniiniii121)()(y=ax+byO(xi ,yi)x其中其中

23、 和和 分别为分别为xi和和yi的平均值的平均值 xy如果建模者判断变量间的关系并非线性关系而是其他类型的函如果建模者判断变量间的关系并非线性关系而是其他类型的函数，则可作数，则可作 变量替换变量替换使之转化为线性关系或用类似方法使之转化为线性关系或用类似方法拟合拟合。第25页/共120页第二十六页，共120页。显然，运动员体重越大，他能举起的重量也越大，但举显然，运动员体重越大，他能举起的重量也越大，但举重成绩和运动员体重到底是怎样关系的，不同量级运动重成绩和运动员体重到底是怎样关系的，不同量级运动员的成绩又如何比较优劣呢？运动成绩是包括生理条件、员的成绩又如何比较优劣呢？运动成绩是包括生理

24、条件、心理因素等等众多相关因素共同作用的结果，要建立精心理因素等等众多相关因素共同作用的结果，要建立精确的模型至少现在还无法办到。但我们拥有大量的比赛确的模型至少现在还无法办到。但我们拥有大量的比赛成绩纪录，根据这些数据不妨可以建立一些经验模型。成绩纪录，根据这些数据不妨可以建立一些经验模型。为简单起见，我们不妨取表中的数据为例。为简单起见，我们不妨取表中的数据为例。例例1（举重成绩的比较）（举重成绩的比较）举重举重是一种一般人都能看懂的运动，它共分是一种一般人都能看懂的运动，它共分九个重量级，有两种主要的比赛方法：抓举九个重量级，有两种主要的比赛方法：抓举和挺举。和挺举。 表中给出了到表中给

25、出了到1977年底为止九个年底为止九个重量级的世界纪录。重量级的世界纪录。255200110以上以上237.518511022118090207.517082.5195157.575180141.567.5161.513060151120.55614110952挺举（公斤）挺举（公斤）抓举（公斤）抓举（公斤）成绩成绩重量级（上限体重量级（上限体重）重）第26页/共120页第二十七页，共120页。模型模型(mxng)1（线（线性模型性模型(mxng)） 将数据画在直角坐标系中可以发现，运动成绩与体将数据画在直角坐标系中可以发现，运动成绩与体量近似满足线性关系，只有量近似满足线性关系，只有110公

26、斤级有点例外，两公斤级有点例外，两项成绩都显得较低。应用前面叙述的方法可求出近项成绩都显得较低。应用前面叙述的方法可求出近似关似关 系式系式L=kB+C，其中，其中B为体重，为体重，L为举重成绩。为举重成绩。你在作图你在作图 时时L轴可以放轴可以放 在在50公斤或公斤或52公斤处，因为公斤处，因为没有更轻级别的比赛，具体计算留给同学自己去完没有更轻级别的比赛，具体计算留给同学自己去完成。成。 第27页/共120页第二十八页，共120页。模型模型(mxng)2（幂函数（幂函数模型模型(mxng)） 线性模型并未得到广泛的接受，要改进结果，能够线性模型并未得到广泛的接受，要改进结果，能够想到的自然

27、首先是幂函数模型，即令想到的自然首先是幂函数模型，即令L=kBa，对此式，对此式取对数，得取对数，得 到到lnL=lnk+a lnB。将原始数据也取对数，。将原始数据也取对数，问题即转化了线性模型，可用最小二乘法求出参数。问题即转化了线性模型，可用最小二乘法求出参数。几十年前英国和爱尔兰采用的比较举重成绩优劣几十年前英国和爱尔兰采用的比较举重成绩优劣 的的Austin公式公式：L=L/B3/4就是用这一方法求得的。就是用这一方法求得的。 第28页/共120页第二十九页，共120页。模型模型(mxng)3（经典模（经典模型型(mxng)） 经典模型是根据生理学中的已知结果和比例关系推导出来的经典

28、模型是根据生理学中的已知结果和比例关系推导出来的公式，应当说，它并不属于经验公式。为建立数学模型，先公式，应当说，它并不属于经验公式。为建立数学模型，先提出如下一些假设：提出如下一些假设： (1)举重成绩举重成绩(chngj)正比于选手肌肉的平均横截正比于选手肌肉的平均横截 面积面积A，即，即L=k1A(2)A正比于身高正比于身高 L的平方，即的平方，即 A=k2L2(3)体重正比于身高体重正比于身高 L的三次方，的三次方， 即即B=k3L3根据上述假设根据上述假设(jish)，可得可得 3232321)(KBkBkkL显然，显然，K越大则成绩越好，故可用越大则成绩越好，故可用 来比较选手来比

29、较选手比赛成绩的优劣。比赛成绩的优劣。 32LBL32321kkkK第29页/共120页第三十页，共120页。模型模型(mxng)4（O Carroll公式）公式） 经验公式的主要依据是比例关系，其假设条件非常粗糙，可经验公式的主要依据是比例关系，其假设条件非常粗糙，可信度不大，因而大多数人认为它不能令人信服。信度不大，因而大多数人认为它不能令人信服。1967年，年，O Carroll基于动物学和统计分析得出了一个现在被广泛使用的基于动物学和统计分析得出了一个现在被广泛使用的公式。公式。O Carroll模型的假设条件是：模型的假设条件是： (1) L=k1Aa, a1 (2) A=k2Lb,

30、 bm）。由）。由于公式量纲齐次当且仅当它可用无量纲的量表示，故方于公式量纲齐次当且仅当它可用无量纲的量表示，故方程当且仅当可写程当且仅当可写 成成f(1，, m)=0时才是量纲齐次的，时才是量纲齐次的，定理证毕。定理证毕。 证证 设设x1,xk为方程中出现的变量与常数为方程中出现的变量与常数, ,对这些变量与对这些变量与常数的任一乘积常数的任一乘积 ,令令 函数函数g建立了建立了xi(i=1,k)的乘积所组成的空间的乘积所组成的空间 与与k维欧氏维欧氏空间之间的一个一一对应。现设涉及到的基本量纲有空间之间的一个一一对应。现设涉及到的基本量纲有n个个,它们它们 为为y1,yn.用这些基本量纲来

31、表达用这些基本量纲来表达 该该xi的乘幂的乘幂,设此设此乘幂的量纲为乘幂的量纲为 令令易见易见dg-1是是k维欧氏空间维欧氏空间 到到n维欧氏空间的一个变换，这维欧氏空间的一个变换，这里的里的g-1为为g的逆变换。的逆变换。 k1aka1xx )a,(a)xg(xk1aka1k1n1bnb1yy )b ,(b)xd(xn1aka1k1第44页/共120页第四十五页，共120页。考察质量集中于距支点为考察质量集中于距支点为 l 的质点上的无阻的质点上的无阻尼尼 单摆，（如图），其运动为某周单摆，（如图），其运动为某周 期期 t 的的左右摆动，现希望左右摆动，现希望(xwng)得到周期得到周期 t

32、 与其他量之间与其他量之间的的 关系。关系。lmg考考察察 ， 的的量量纲纲(lin n)为为MaLb+dTc-2b若若 无无量量纲纲(lin n)，则则有有ab c dem g t l 000e = 0abdc2b第45页/共120页第四十六页，共120页。量纲分析法虽然简单，但使用时在技巧方面的要求较高，稍量纲分析法虽然简单，但使用时在技巧方面的要求较高，稍一疏忽就会导出荒谬的结果或根本得不出任何有用的结果。一疏忽就会导出荒谬的结果或根本得不出任何有用的结果。首先，它要求建模者对研究的问题有正确而充分的了解，能首先，它要求建模者对研究的问题有正确而充分的了解，能正确列出与该问题相关的量及相

33、关的基本量纲，容易看出，正确列出与该问题相关的量及相关的基本量纲，容易看出，其后的分析正是通过对这些量的量纲研究而得出的，列多或其后的分析正是通过对这些量的量纲研究而得出的，列多或列少均不可能得出有用的结果。其次，在为寻找无量纲量而列少均不可能得出有用的结果。其次，在为寻找无量纲量而求解齐次线性方程组时，基向量组有无穷多种取法，如何选求解齐次线性方程组时，基向量组有无穷多种取法，如何选取也很重要，此时需依靠经验，并非任取一组基都能得出有取也很重要，此时需依靠经验，并非任取一组基都能得出有用的结果。此外，建模者在使用量纲分析法时对结果也不应用的结果。此外，建模者在使用量纲分析法时对结果也不应抱有

34、不切实际的过高要求，量纲分析法的基础是公式的量纲抱有不切实际的过高要求，量纲分析法的基础是公式的量纲齐次性，仅凭这一点又怎么可能得出十分深刻的结果，例如，齐次性，仅凭这一点又怎么可能得出十分深刻的结果，例如，公式可能包含某些无量纲常数或无量纲变量，对它们之间的公式可能包含某些无量纲常数或无量纲变量，对它们之间的关系，量纲分析法根本无法加以研究。关系，量纲分析法根本无法加以研究。第46页/共120页第四十七页，共120页。八人赛艇比赛和举重比赛一样，分成八人赛艇比赛和举重比赛一样，分成(fn chn)86公斤公斤的重量级和的重量级和 73公斤的轻量级。公斤的轻量级。1971年，年，T.A.McM

35、ahon比较了比较了1964-1970年期间两次年期间两次奥运会和两次世锦赛成绩，发现奥运会和两次世锦赛成绩，发现 86公斤级比公斤级比73公斤级的成绩大约好公斤级的成绩大约好5%，产生这一差异的，产生这一差异的原因何在呢？原因何在呢？ 我们将以我们将以L表示轻量级、以表示轻量级、以H表示重表示重量级，用量级，用S表示赛艇的浸水面积，表示赛艇的浸水面积，v表示赛艇速度，表示赛艇速度，W表示选手体重，表示选手体重，P表示选手的输出功率，表示选手的输出功率，I表示赛程，表示赛程，T表示比赛成绩（时间）。表示比赛成绩（时间）。 第47页/共120页第四十八页，共120页。考察优秀赛艇选手在比赛中的实

36、际表现可以发现考察优秀赛艇选手在比赛中的实际表现可以发现(fxin)，整个赛程大，整个赛程大致可以分三个阶段，致可以分三个阶段， 即初始时刻的加速阶段、中途的匀速阶段和到达终即初始时刻的加速阶段、中途的匀速阶段和到达终点的冲刺阶段点的冲刺阶段 。由于赛程较长，可以略去前后两段而只考虑中间一段。由于赛程较长，可以略去前后两段而只考虑中间一段 ，为此，提出以下建模假设。为此，提出以下建模假设。（1）设赛艇浸水部分的摩擦力是唯一阻力，摩擦力）设赛艇浸水部分的摩擦力是唯一阻力，摩擦力f正比正比 于于Sv2,（见流体力学），空气阻力等其他因素不计。（见流体力学），空气阻力等其他因素不计。（2）同一量级的

37、选手有相同的体重）同一量级的选手有相同的体重W，选手的输出功，选手的输出功 率率P正比于正比于W，且效率大体相同。，且效率大体相同。由由假设假设1，3svfvp，故，故13pvs 竞赛成绩竞赛成绩13IsTvp记比例系数记比例系数 为为k，则有，则有：31HHH31LLLPSkT ,PSkT第48页/共120页第四十九页，共120页。31HL31LHHLSSPPTT故故由由假设假设2， LHLHWWPP31HL31LHHLSSWWTT故故令令WH=86,WL=73,则有则有由于由于SL略小于略小于SH，故轻量级所化时间，故轻量级所化时间(shjin)比重量级比重量级所化时间所化时间(shjin

38、)约约 多多5%左右。左右。31HLHLSS TT1.056第49页/共120页第五十页，共120页。将一张四条腿的方桌放在不平的地面上，不将一张四条腿的方桌放在不平的地面上，不 允许将桌子移到别处，但允许其绕中心允许将桌子移到别处，但允许其绕中心(zhngxn)旋转旋转 ，是否总能设法使其四条腿同时落地？，是否总能设法使其四条腿同时落地？ 不附加任何条件不附加任何条件(tiojin)，答案，答案 显然显然 是否定的，是否定的， 因此我们因此我们假设假设 （1）地面为连续曲面地面为连续曲面 （2）方桌的四条腿长度相同方桌的四条腿长度相同 （3）相对于地面的弯曲程相对于地面的弯曲程度而言，方桌的

39、腿是足够长度而言，方桌的腿是足够长的的 （4）方桌的腿只要有一点接触地方桌的腿只要有一点接触地面就算着地。面就算着地。总可以使三条腿总可以使三条腿同时着地。同时着地。 第50页/共120页第五十一页，共120页。现在，我们来证明：如果上述假设条件成立，那么答案是肯定现在，我们来证明：如果上述假设条件成立，那么答案是肯定的。以方桌的中心为坐标原点作直角坐标系如的。以方桌的中心为坐标原点作直角坐标系如 图所示，方桌图所示，方桌的四条腿分别在的四条腿分别在A、B、C、D处，处，A、C的初始位置在的初始位置在x轴上，轴上，而而B、D则在则在y轴上，当方桌绕中轴上，当方桌绕中 心心0旋转时，对角线旋转时

40、，对角线 AC与与x轴轴的夹角记为的夹角记为。容易看出，当四条腿尚未全部着地时，腿到地面的距离是不确容易看出，当四条腿尚未全部着地时，腿到地面的距离是不确定的。为消除这一不确定性，令定的。为消除这一不确定性，令 f()为为A、C离地距离之和，离地距离之和，g()为为B、D离地距离之和，它们的值离地距离之和，它们的值 由由唯一确定。由唯一确定。由假设假设（1），），f()、g()均为均为的连续函数。又的连续函数。又 由由假设（假设（3），），三条腿三条腿总能同时着地，总能同时着地， 故故f()g()=0必成立（必成立（ ）。不妨设）。不妨设f(0)=0,g(0)0（若（若g(0)也为也为0，则初

41、始时刻已四条腿着地，不必，则初始时刻已四条腿着地，不必再旋转），于是问题归结为：再旋转），于是问题归结为：yxCDABo已知已知f()、g()均为均为的连续函数，的连续函数，f(0)=0,g(0)0且对任意且对任意有有f()g()=0，求证存在某一，求证存在某一0，使，使f(0)=g(0)=0。第51页/共120页第五十二页，共120页。 （证法一）（证法一）当当=/2时，时，AC与与BD互换位置，故互换位置，故f(/2)0 , g(/2)=0。作。作h()=f()-g()，显然，显然，h()也是也是的连续函数，的连续函数，h(0)=f(0)-g(0)0，由连续函数的取，由连续函数的取零值定理

42、，存在零值定理，存在 o，0o 0,g(/2)=0。令。令o =sup |f ()=0，0,显然显然0 0，总有，总有0且且0。因为。因为f(0+)g (o+)=0，故必有，故必有g (0+)=0，由，由可任意小且可任意小且g连续，可知必连续，可知必 有有 g (0)=0，证毕。证法，证毕。证法二除用二除用 到到f、g的连续性外，还用到了上确界的性质。的连续性外，还用到了上确界的性质。 第52页/共120页第五十三页，共120页。在解决实际问题时，注意观察和善于想象是十分重要的，在解决实际问题时，注意观察和善于想象是十分重要的，观察与想象不仅能发现问题隐含的某些属性，有时还能顺观察与想象不仅能

43、发现问题隐含的某些属性，有时还能顺理成章地找到解决实际问题的钥匙。本节的几个例子说明，理成章地找到解决实际问题的钥匙。本节的几个例子说明，猜测也是一种想象力。没有合理而又大胆的猜测，很难做猜测也是一种想象力。没有合理而又大胆的猜测，很难做出具有创新性的结果。开普勒的三大定律（尤其是后两条）出具有创新性的结果。开普勒的三大定律（尤其是后两条）并非一眼就能看出的，它们隐含在行星运动的轨迹之中，并非一眼就能看出的，它们隐含在行星运动的轨迹之中，隐含在第谷记录下来的一大堆数据之中。历史上这样的例隐含在第谷记录下来的一大堆数据之中。历史上这样的例子实在太多了。在获得了一定数量的资料数据后，人们常子实在太

44、多了。在获得了一定数量的资料数据后，人们常常会先去猜测某些结果，然后试图去证明它。猜测一经证常会先去猜测某些结果，然后试图去证明它。猜测一经证明就成了定理，而定理一旦插上想象的翅膀，又常常会被明就成了定理，而定理一旦插上想象的翅膀，又常常会被推广出许多更为广泛的结果。即使猜测被证明是错误的，推广出许多更为广泛的结果。即使猜测被证明是错误的，结果也决不是一无所获的失败而常常是对问题的更为深入结果也决不是一无所获的失败而常常是对问题的更为深入的了解。的了解。 2.9最短路径与最速方案问题最短路径与最速方案问题 第53页/共120页第五十四页，共120页。设有一个半径为设有一个半径为 r 的圆形湖，

45、圆心的圆形湖，圆心(yunxn)为为 O。A、B 位于湖的两侧，位于湖的两侧，AB连线过连线过O，见图。，见图。现拟从现拟从A点步行到点步行到B点，在不得进入湖中的限点，在不得进入湖中的限 制下，问怎样的路径最近。制下，问怎样的路径最近。 ABOr将湖想象成凸出地面的木桩，将湖想象成凸出地面的木桩， 在在AB间拉一根软线，当间拉一根软线，当线被拉紧时将得到最短路径。根据这样的想象，猜测线被拉紧时将得到最短路径。根据这样的想象，猜测 可以如下得到最短路径：可以如下得到最短路径： 过过A作圆的切线切圆于作圆的切线切圆于E，过，过B作圆的切线切圆作圆的切线切圆 于于F。最短路径为由线。最短路径为由线

46、 段段AE、弧、弧EF和线段和线段FB连接而成的连续曲线（根据对称性，连接而成的连续曲线（根据对称性，AE，弧弧EF，FB连接而成的连续曲线也是）。连接而成的连续曲线也是）。EFEF第54页/共120页第五十五页，共120页。以上只是一种猜测以上只是一种猜测(cic)，现在来证明这一猜测，现在来证明这一猜测(cic)是正确是正确的。为此，先介绍一下凸集与凸集的性质。的。为此，先介绍一下凸集与凸集的性质。定义定义2.1（凸集凸集）称集合）称集合 R为凸集，若为凸集，若x1、x2R及及0，1，总有总有x1+（1-）x2R。即若。即若x1、x2R，则，则x1、x2的的连线必整个地落连线必整个地落 在

47、在R中。中。定理定理2.2（分离定理分离定理）对平面中的凸）对平面中的凸 集集R与与R外的一点外的一点K，存在直线存在直线 l , l 分离分离R与与K，即，即R与与K分别位于分别位于 l 的两侧（注：的两侧（注：对一般的凸对一般的凸 集集R与与R外的一点外的一点K，则存在超平面分，则存在超平面分 离离R与与K），见图。），见图。klR下面证明下面证明(zhngmng)猜猜想想第55页/共120页第五十六页，共120页。猜测证明猜测证明(zhngmng)如下：如下：（方法一）（方法一）显然，显然， 由由AE、EF、FB及及AE，EF，FB围成围成的区域的区域 R是一凸集。利用是一凸集。利用分离

48、定理分离定理易证最短径不可能经过易证最短径不可能经过R外的点，若不然，设外的点，若不然，设 为最短路径，为最短路径，过过R外的一点外的一点M，则必存在直则必存在直 线线l分离分离M与与R，由于路径，由于路径是连续曲线，由是连续曲线，由A沿沿到到M，必交，必交l于于M1，由，由M沿沿到到B又必交又必交l于于M2。这样，直。这样，直线线 段段M1M2的长度必小于路的长度必小于路 径径M1MM2的长度，与的长度，与是是A到到B的最短路径矛盾，至此，我们已证明最短路径必在凸集的最短路径矛盾，至此，我们已证明最短路径必在凸集R内。不妨设路径经湖的上方到达内。不妨设路径经湖的上方到达B点，则弧点，则弧EF

49、必在路径必在路径F上，又直线段上，又直线段AE是由是由A至至E的最短路径，直线的最短路径，直线FB是由是由F到到B的最短路径，猜测得证。的最短路径，猜测得证。ABOrEFEFM1M2Ml第56页/共120页第五十七页，共120页。还可用微积分方法求弧长，根据计算证明满足还可用微积分方法求弧长，根据计算证明满足限止条件的其他连续曲线必具有更大的长度；限止条件的其他连续曲线必具有更大的长度；此外，本猜测也可用平面几何此外，本猜测也可用平面几何(pngminjh)知识加以证明等。知识加以证明等。 根据猜测不难看出，根据猜测不难看出， 例例5中的条件可以大大中的条件可以大大放松，可以不必放松，可以不必

50、 设设AB过圆心，甚至可不必过圆心，甚至可不必设湖是圆形的。例如对设湖是圆形的。例如对 下图，我们可断定由下图，我们可断定由A至至B的最短路径必的最短路径必 为为l1与与l2之一，其证明也之一，其证明也不难类似给出。不难类似给出。 ABl1l2D到此为止，我们的研讨还只局限于平面之中，其到此为止，我们的研讨还只局限于平面之中，其实上述猜测可十分自然地推广到一般空间实上述猜测可十分自然地推广到一般空间(kngjin)中去。中去。1973年，年，J.W.Craggs证明了以证明了以上结果：上结果：若可行区域的边界是光滑曲面。则最短路径必由下列弧组若可行区域的边界是光滑曲面。则最短路径必由下列弧组成

51、，它们或者是空间中的自然最短曲线，或者是可行区域成，它们或者是空间中的自然最短曲线，或者是可行区域的边界弧。而且，组成最短路径的各段弧在连接点处必定的边界弧。而且，组成最短路径的各段弧在连接点处必定相切。相切。第57页/共120页第五十八页，共120页。例例6 6 一辆汽车一辆汽车(qch)(qch)停于停于 A A处并垂直于处并垂直于ABAB方向，此方向，此汽车汽车(qch)(qch)可转的最小圆半径为可转的最小圆半径为 R R，求不倒车而由，求不倒车而由 A A到到B B的最短路径。的最短路径。解解（情况（情况1）若若|AB|2R，最短路径由，最短路径由 弧弧AC与切线与切线BC组组成（见

52、成（见图图 ）。）。（情况（情况2）若若|AB|0为推力，为推力，fS，故由连续函数的性质存在，故由连续函数的性质存在 某某TT，S（T）=S但这但这一结果与一结果与=(t)是最优方案下的车速的假设矛盾，因为用我是最优方案下的车速的假设矛盾，因为用我们猜测的推车方法推车，只们猜测的推车方法推车，只 需需T时间即可将车推到修车处，时间即可将车推到修车处， 而而TT。T T0 0S S( (t t) )d dt ta ad dt td dv vb botATTASy=aty=-b(t-T)第61页/共120页第六十二页，共120页。 圆周率是人类获得的最古老的数学概念圆周率是人类获得的最古老的数学

53、概念之一，早在大约之一，早在大约37003700年前（即公元前年前（即公元前17001700年左右）的古埃及人就已经在年左右）的古埃及人就已经在 用用256/81256/81（约约3.16053.1605）作为）作为的近似值了。几千年来的近似值了。几千年来，人们，人们(rn men)(rn men)一直没有停止过求一直没有停止过求的的努力。努力。2.10 的计算的计算第62页/共120页第六十三页，共120页。 古古 典典 方方 法法 分分 析析 方方 法法 其其 它它 方方 法法 概率概率(gil)(gil)方法方法 数值积分方数值积分方法法第63页/共120页第六十四页，共120页。 古典

54、古典(gdin)(gdin)方法方法用什么方法来计用什么方法来计 算算的近似值呢？显然的近似值呢？显然(xinrn)(xinrn)，不可能，不可能仅根据圆周率的定义，用圆的周长去除以直径。起先，人们仅根据圆周率的定义，用圆的周长去除以直径。起先，人们采用的都是用圆内接正多边形和圆外切正多边形来逼近的古采用的都是用圆内接正多边形和圆外切正多边形来逼近的古典方法。典方法。6边形边形12边形边形24边形边形圆圆第64页/共120页第六十五页，共120页。 阿基米德曾用圆内接阿基米德曾用圆内接 96 96边形和圆外切边形和圆外切(wi qi)96(wi qi)96边形夹逼的方法证明了边形夹逼的方法证明

55、了72271223 由由和和 导出导出 tansin 96 公元公元5 5世纪世纪(shj)(shj)，祖冲之指出，祖冲之指出3.14159273.1415926 比西方得到同样结比西方得到同样结果果(ji gu)(ji gu)几乎几乎早了早了10001000年年第65页/共120页第六十六页，共120页。 十五世纪中叶，阿尔十五世纪中叶，阿尔卡西给出卡西给出的的1616位小位小数，打破数，打破(d p)(d p)了祖冲之的纪录了祖冲之的纪录 1579 1579年年, ,韦达证明韦达证明(zhngmng)(zhngmng)373.14159265 353.14159265 1630 1630年

56、年, ,最后最后(zuhu)(zuhu)一位用古典方一位用古典方法求法求的人格林伯格也只求到了的人格林伯格也只求到了的第的第3939位小数位小数第66页/共120页第六十七页，共120页。 分析方法分析方法从十七世纪中叶起，人们开始用更先进的分析从十七世纪中叶起，人们开始用更先进的分析方法来求方法来求的近似值，其中应用的主要工具的近似值，其中应用的主要工具(gngj)(gngj)是收敛的无穷乘积和无穷级数，在是收敛的无穷乘积和无穷级数，在本节中我们将介绍一些用此类方法求本节中我们将介绍一些用此类方法求近似值近似值的实例。的实例。第67页/共120页第六十八页，共120页。067702.3212

57、01920543432122103516.341403940543432122 取取20 k取取10 k 1656 1656年，沃里斯年，沃里斯(Wallis)(Wallis)证明证明(zhngmng)(zhngmng) 112212227656543432122kkkkk 第68页/共120页第六十九页，共120页。 在微积分中我们在微积分中我们(w men)(w men)学过泰勒级数，其中有学过泰勒级数，其中有12)1(53arctan12053 kxxxxxkkk),( x当当1 x121)1(5131140 kkk 第69页/共120页第七十页，共120页。取取20 k189184.

58、3411391513114 取取10 k232316. 3211191513114 第70页/共120页第七十一页，共120页。 在中学数学中证明过下面在中学数学中证明过下面(xi (xi mian)mian)的等式的等式31arctan21arctan1arctan4 左边三个正方形左边三个正方形组成的矩形中，组成的矩形中， 由由 和和 可得可得CBA DC 和和 的展开式的收敛速度的展开式的收敛速度都比都比 快得多快得多21arctan1arctan31arctanACBD第71页/共120页第七十二页，共120页。 麦琴麦琴( (Machin) )给出给出2391arctan51arct

59、an44 (Machin公式公式) )51arctan 44 2391tan 记记 ， ，得，得此式求得了此式求得了的第的第100位小数且全部正确位小数且全部正确第72页/共120页第七十三页，共120页。 其它其它(qt)(qt)方法方法除用古典除用古典(gdin)(gdin)方法与分析方法求方法与分析方法求的的近似值以外，还有人用其他方法来求近似值以外，还有人用其他方法来求的近的近似值。这里我们将介绍两种方法：似值。这里我们将介绍两种方法： 概率概率(gil)(gil)方法方法 数值积分方法数值积分方法第73页/共120页第七十四页，共120页。 概率概率(gil)(gil)方法方法取一个

60、二维数组（取一个二维数组（x,yx,y），取一个充分），取一个充分(chngfn)(chngfn)大的正整大的正整 数数n n，重复，重复n n次，每次独次，每次独立地从立地从 （0 0，1 1）中随机地取一对）中随机地取一对 数数x x和和y y ，分，分别检验别检验x2+y21x2+y21是否成立。是否成立。 设设n n次试验中等式次试验中等式成立的共有成立的共有m m次，令次，令4m/n4m/n。但这种方法很难得到但这种方法很难得到的较好的近似值。的较好的近似值。第74页/共120页第七十五页，共120页。 数值积分方法数值积分方法(fngf) 10214dxx 102114dxx 还可