1、 高等数学是现代各科知识的理论基础,在数高等数学是现代各科知识的理论基础,在数 学建模中有广泛的应用,极限、连续和积分学建模中有广泛的应用,极限、连续和积分等数学思想是建立数学模型的基本思想,抽等数学思想是建立数学模型的基本思想,抽象思维和逻辑思维能力是数学建模必备的能力。象思维和逻辑思维能力是数学建模必备的能力。在教学中,融入数学建模思想和方法,让学生在教学中,融入数学建模思想和方法,让学生养成数学建模的习惯。养成数学建模的习惯。培养他们建立数学模型和解决数学模型的能力。培养他们建立数学模型和解决数学模型的能力。高等数学在数学建模中的应用高等数学在数学建模中的应用举例举例某航空母舰派其护卫舰
2、去搜寻其跳伞的飞某航空母舰派其护卫舰去搜寻其跳伞的飞 行行员,护卫舰找到飞行员后,航母通知它尽快员,护卫舰找到飞行员后,航母通知它尽快 返回与其汇合并通报了航母当前的航速与方返回与其汇合并通报了航母当前的航速与方 向,问护卫舰应怎样航行,才能与航母汇合。向,问护卫舰应怎样航行,才能与航母汇合。例例1 舰艇的会舰艇的会合合12,11222aabrbaah令:令:则上式可简记成则上式可简记成:222rh-yx)(A(0,b)XYB(0,-b)P(x,y)O航母航母 护卫舰护卫舰 1 2)()(22222b-yx a byx即:即:22222222)1(411ababaayx可化为:可化为:记记v2
3、/v1=a通常通常a1 222|AP|a|BP|则则汇合点汇合点 p必位于此圆上。必位于此圆上。bxy)(tan1(护卫舰的路线方程)(护卫舰的路线方程)bxy)(tan2(航母的路线方程(航母的路线方程)即可求出即可求出P点的坐标和点的坐标和2 的值。的值。本模型虽简单,但分析极清晰且易于本模型虽简单,但分析极清晰且易于实际应用实际应用 例例2 双层玻璃的功效双层玻璃的功效在寒冷的北方,在寒冷的北方,许多住房的许多住房的 玻璃窗都是双层玻璃窗都是双层玻璃的,现在我们来建立一个简单玻璃的,现在我们来建立一个简单 的数学模的数学模型,研究一下双层玻璃到底有多型,研究一下双层玻璃到底有多 大的功效
4、。大的功效。比较两座其他条件完全相同的房屋,它们比较两座其他条件完全相同的房屋,它们 的的差异仅仅在窗户不同。差异仅仅在窗户不同。不妨可以提出以下不妨可以提出以下 假设:假设:1、设室内热量的流失是热传导引起的,不存在户、设室内热量的流失是热传导引起的,不存在户内外的空气对流。内外的空气对流。2、室内温、室内温 度度T1与户外温与户外温 度度T2均为常数。均为常数。3、玻璃是均匀的,热传导系数为常数。、玻璃是均匀的,热传导系数为常数。设玻璃的热传导系数设玻璃的热传导系数 为为k1,空气的热传导系数,空气的热传导系数 为为k2,单位时间通过单位面积由温度高的一侧流向温度低的一单位时间通过单位面积
5、由温度高的一侧流向温度低的一侧的热量为侧的热量为 ddl室室外外T2室室内内T1TaTb由热传导公式由热传导公式 =kT/d dTTklTTkdTTkbbaa21211)/()(21212121dklkTTdklkTa解得:解得:dklkdTTkddklkTTdklkTk212112121211122)1(此函数的图形为此函数的图形为dd室室外外T2室室内内T1dTTk2211)/()(2221dklk类似有类似有 321621kk一般一般dl/811故故记记h=l/d并令并令f(h)=181h01234567891000.10.20.30.40.50.60.70.80.91hf(h)考虑到美
6、观和使用上考虑到美观和使用上 的方便,的方便,h不必取得过大,例如,可不必取得过大,例如,可 取取h=3,即,即l=3d,此时房屋热量,此时房屋热量的损失不超过单层玻璃窗时的的损失不超过单层玻璃窗时的 3%。例例3 崖高的估算崖高的估算假如你站在崖顶且身上带着一只具有跑表功假如你站在崖顶且身上带着一只具有跑表功 能的计算器,你也许会出于好奇心想用扔下能的计算器,你也许会出于好奇心想用扔下 一块石头听回声的方法来估计山崖的高度,一块石头听回声的方法来估计山崖的高度,假定你能准确地测定时间,你又怎样来推算假定你能准确地测定时间,你又怎样来推算 山崖的高度呢,请你分析一下这一问题。山崖的高度呢,请你
7、分析一下这一问题。我有一只具有跑我有一只具有跑 表功能的计表功能的计算器。算器。方法一方法一假定空气阻力不计,可以直接利用自由落体运动的公式假定空气阻力不计,可以直接利用自由落体运动的公式来计算。例如,来计算。例如,设设t=4秒,秒,g=9.81米米/秒秒2,则可求得,则可求得h78.5米。米。221gth 我学过微积分,我可以做我学过微积分,我可以做 得更好,呵呵。得更好,呵呵。vKmgdtdvmF除去地球吸引力外,对石块下落影响最大的当除去地球吸引力外,对石块下落影响最大的当 属空气阻力。根据流体力学知识,此时可属空气阻力。根据流体力学知识,此时可设空气阻力正比于石块下落的速度,阻力系设空
8、气阻力正比于石块下落的速度,阻力系 数数K为常数,因而,由牛顿第二定律可得:为常数,因而,由牛顿第二定律可得:kgcevkt令令k=K/m,解得,解得 代入初始条件代入初始条件 v(0)=0,得,得c=g/k,故有,故有 ktekgkgv再积分一次,得:再积分一次,得:cekgtkghkt2若设若设k=0.05并仍设并仍设 t=4秒,则可求秒,则可求 得得h73.6米。米。听到回声再按跑表,计算得到的时间中包含了听到回声再按跑表,计算得到的时间中包含了 反应时间反应时间 不妨设平均反应时间不妨设平均反应时间 为为0.1秒秒,假如仍,假如仍 设设t=4秒,扣除反应时间后应秒,扣除反应时间后应 为
9、为3.9秒,代入秒,代入 式式,求得,求得h69.9米。米。222)1(kgektkgkgekgtkghktkt多测几次,取平均值多测几次,取平均值代入初始条代入初始条 件件h(0)=0,得到计算山崖高度的公式:,得到计算山崖高度的公式:将将e-kt用泰勒公式展开并用泰勒公式展开并 令令k 0+,即可得出前面不考虑空气阻,即可得出前面不考虑空气阻力时的结果。力时的结果。还应考虑回声传回来所需要的时间。为此,令石块下落还应考虑回声传回来所需要的时间。为此,令石块下落 的真正时间的真正时间 为为t1,声音传回来的,声音传回来的时间记时间记 为为t2,还得解一个方程组:,还得解一个方程组:93340
10、1212211.ttthkg)ekt(kghkt这一方程组是非线性这一方程组是非线性的,求解不太容易,的,求解不太容易,为了估算崖高竟要去为了估算崖高竟要去解一个非线性主程组解一个非线性主程组似乎不合情理似乎不合情理 相对于石块速度,声音速度要快得多,我们可相对于石块速度,声音速度要快得多,我们可 用方法二先求一次用方法二先求一次 h,令,令t2=h/340,校正,校正t,求石块下落时间,求石块下落时间 t1t-t2将将t1代入式代入式再算一次,得出崖再算一次,得出崖高的近似值。例如,高的近似值。例如,若若h=69.9米,则米,则 t20.21秒,故秒,故 t13.69秒,求得秒,求得 h62
11、.3米。米。例例4 录像带还能录多长时间录像带还能录多长时间录像机上有一个四位计数器,一盘录像机上有一个四位计数器,一盘 180分钟分钟的录像带在开始计数时为的录像带在开始计数时为 0000,到结束时计,到结束时计数为数为1849,实际走时为,实际走时为185分分20秒。我们从秒。我们从0084观察到观察到0147共用时间共用时间3分分21秒。若录像秒。若录像机目前的计数为机目前的计数为1428,问是否还能录下一个,问是否还能录下一个 60分钟的节目?分钟的节目?rRl由由 vt)r(R22得到得到212rvtR又又 因和因和 得得 Rl tvl tRv 积分得到积分得到tdt)rvtv(d0
12、2120rrvt2rvt2t2212021)()(即即从而有从而有rrvtn212)(12我们希望建立一个录像带已录像时我们希望建立一个录像带已录像时 间间t与计数器计与计数器计 数数n之间的函数关系。为建立一之间的函数关系。为建立一个正确的模型,首个正确的模型,首 先必须搞清哪些量是常量,哪些量是变量。首先,录像先必须搞清哪些量是常量,哪些量是变量。首先,录像 带的磁带的磁带的厚带的厚 度是度是 常量,它被绕在一个半径常量,它被绕在一个半径 为为r的园盘上,见图。磁带转动中的线速的园盘上,见图。磁带转动中的线速 度度v显然也是常数,否则图象声音必然会失真。此外,计数器的读显然也是常数,否则图
13、象声音必然会失真。此外,计数器的读 数数n与转过的圈与转过的圈数有关,从而与转过的角数有关,从而与转过的角 度度成正比。成正比。rRlrrvtn212)(12 此式中的三个参数此式中的三个参数、v和和r均不易精确测得,虽然我们可以从上式均不易精确测得,虽然我们可以从上式解出解出t与与n的函数关系,但效果不佳,故令的函数关系,但效果不佳,故令 则可将上式简化为:则可将上式简化为:v v/r2tn故故nnnt21222令令21a b2上式又可化简记成上式又可化简记成 t=an2+bn t=an2+bn rRl上式以上式以a、b为参数显然是一个十分明智的做法,它为公式的最终确为参数显然是一个十分明智
14、的做法,它为公式的最终确立即参数求解提供了方便。将已知条件代入,得方程组:立即参数求解提供了方便。将已知条件代入,得方程组:3.351471478484185.331849(1849)12122tbatbaba从后两式中消从后两式中消 去去t1,解得,解得a=0.0000291,b=0.04646,故故t=0.0000291 n2+0.04646n,令,令n=1428,得到,得到t=125.69(分)由于一盒录像带实际可录像时间为(分)由于一盒录像带实际可录像时间为185.33分,故尚可分,故尚可录像时间录像时间 为为59.64分,已不能再录下一个分,已不能再录下一个60分钟的节目了。分钟的节
15、目了。例例5 5 将形状质量相同的砖块一一向右往外叠放,欲尽可能地延伸将形状质量相同的砖块一一向右往外叠放,欲尽可能地延伸到远方,问最远可以延伸多大距离。到远方,问最远可以延伸多大距离。设砖块是均质的,长度与重量均设砖块是均质的,长度与重量均 为为1 1,其,其 重心在中点重心在中点1/21/2砖长处,现砖长处,现用归纳法推导。用归纳法推导。Zn(n1)n(n1)由第由第 n块砖受到的两个力的力矩相等,有:块砖受到的两个力的力矩相等,有:1/2-Zn=(n1)Zn故故Zn=1/(2n),从而上面,从而上面 n块砖向右推出的总距离为块砖向右推出的总距离为 ,nkk121112121,时nnknk
16、n 例例6 6 某人住在某公交线附近,该公交线路为在某人住在某公交线附近,该公交线路为在A A、B B两地间运行,两地间运行,每隔每隔 1010分钟分钟A A、B B两地各发出一班车,此人常在离家最近的两地各发出一班车,此人常在离家最近的 C C点等点等车,他发现了一个令他感到奇怪的现象:在绝大多数情况下,先车,他发现了一个令他感到奇怪的现象:在绝大多数情况下,先到站的总是由到站的总是由 B B去去A A的车,难道由的车,难道由 B B去去A A的车次多些吗?请你帮助的车次多些吗?请你帮助他找一下原因他找一下原因由于距离不同,设由于距离不同,设 A A到到C C行驶行驶3131分钟,分钟,B
17、B到到C C要行驶要行驶 3030分钟,考察一个时间长度分钟,考察一个时间长度 为为1010分钟的区间,例如,可分钟的区间,例如,可以从以从 A A方向来的车驶方向来的车驶 离离C C站时开始,在其后的站时开始,在其后的 9 9分钟分钟内到达的乘客见到先来的车均为内到达的乘客见到先来的车均为 B B开往开往A A的,仅有最的,仅有最 后后1 1分钟到达的乘客才见到分钟到达的乘客才见到 由由A A来的车先到。由此可见,来的车先到。由此可见,如果此人如果此人 到到C C站等车的时间是随机的,则他先遇站等车的时间是随机的,则他先遇 上上B B方向来的车的概率为方向来的车的概率为 90%90%。例例7
18、 方桌问题方桌问题将一张四条腿的方桌放在不平的地面上,不将一张四条腿的方桌放在不平的地面上,不 允许将桌子移到别处,但允许其绕中心旋转允许将桌子移到别处,但允许其绕中心旋转 ,是否总能设法使其四条腿同时落地?,是否总能设法使其四条腿同时落地?不附加任何条件,答案不附加任何条件,答案 显然显然 是否定的,是否定的,因此我们假设因此我们假设 (1)地面为连续曲面)地面为连续曲面(2)方桌的四条腿)方桌的四条腿长度相同长度相同 (3)相对于)相对于地面的弯曲程度而言,方桌的腿是足够长地面的弯曲程度而言,方桌的腿是足够长的的 (4)方桌的腿只要)方桌的腿只要有一点接触地面就算着地。有一点接触地面就算着
19、地。总可以使三条腿同时着地。总可以使三条腿同时着地。现在,我们来证明:如果上述假设条件成立,那么答案是肯定的。以方桌的中心为坐标原点作直现在,我们来证明:如果上述假设条件成立,那么答案是肯定的。以方桌的中心为坐标原点作直角坐标系如角坐标系如 图所示,方桌的四条腿分别在图所示,方桌的四条腿分别在A、B、C、D处,处,A、C的初始位置在的初始位置在x轴上,而轴上,而B、D则在则在y轴上,当方桌绕中轴上,当方桌绕中 心心0旋转时,对角线旋转时,对角线 AC与与x轴的夹角记为轴的夹角记为。容易看出,当四条腿尚未全部着地时,腿到地面的距离是不确定的。为消除这一不确定性,令容易看出,当四条腿尚未全部着地时
20、,腿到地面的距离是不确定的。为消除这一不确定性,令 f()为为A、C离地距离之和,离地距离之和,g()为为B、D离地距离之和,它们的值离地距离之和,它们的值 由由唯一确定。由假设(唯一确定。由假设(1),),f()、g()均为均为的连续函数。又的连续函数。又 由假设(由假设(3),三条腿总能同时着地,),三条腿总能同时着地,故故f()g()=0必成立(必成立()。)。不妨设不妨设f(0)=0,g(0)0(若(若g(0)也为也为0,则初始时刻已四条腿着地,不必再旋转),于是问题归结为:,则初始时刻已四条腿着地,不必再旋转),于是问题归结为:yxCDABo已知已知f()、g()均为均为的连续函数,
21、的连续函数,f(0)=0,g(0)0且对任意且对任意有有f()g()=0,求证存在某一,求证存在某一0,使,使f(0)=g(0)=0。(证法一)当(证法一)当=/2时,时,AC与与BD互换位置,故互换位置,故f(/2)0,g(/2)=0。令。令h()=f()-g(),显然,显然,h()也是也是的连续函数,的连续函数,h(0)=f(0)-g(0)0,由连续函数的取零值定理,由连续函数的取零值定理,存在存在 o,0o 0,g(/2)=0。令。令o=sup|f()=0,0,显然显然0 0,总有,总有0且且0。因为。因为f(0+)g(o+)=0,故必有,故必有g(0+)=0,由,由可任意小且可任意小且
22、g连续,可知必连续,可知必 有有 g(0)=0,证毕。证法二除用证毕。证法二除用 到到f、g的连续性外,还用到了上确界的性质。的连续性外,还用到了上确界的性质。圆周率是人类获得的最古老的数学概念之一,早在大约圆周率是人类获得的最古老的数学概念之一,早在大约37003700年年前(即公元前前(即公元前17001700年左右)的古埃及人就已经在年左右)的古埃及人就已经在 用用256/81256/81(约(约3.16053.1605)作为)作为的近似值了。几千年来,人们一直没有停止过求的近似值了。几千年来,人们一直没有停止过求的努力。的努力。例例8 的计算的计算 古古 典典 方方 法法 分分 析析
23、方方 法法 其其 它它 方方 法法 概率方法概率方法 数值积分方法数值积分方法 古典方法古典方法用什么方法来计用什么方法来计 算算的近似值呢?显然,不可能仅根据圆周率的定义,用圆的周长的近似值呢?显然,不可能仅根据圆周率的定义,用圆的周长去除以直径。起先,人们采用的都是用圆内接正多边形和圆外切正多边形来逼近的古去除以直径。起先,人们采用的都是用圆内接正多边形和圆外切正多边形来逼近的古典方法。典方法。6边形边形12边形边形24边形边形圆圆 阿基米德曾用圆内接阿基米德曾用圆内接 9696边形和圆外切边形和圆外切9696边形夹逼的方法证明边形夹逼的方法证明了了72271223 由由和和 导出导出 t
24、ansin 96 公元公元5世纪,祖冲之指出世纪,祖冲之指出3.14159273.1415926 比西方得到同样结果几乎早比西方得到同样结果几乎早了了1000年年 十五世纪中叶,阿尔十五世纪中叶,阿尔卡西给出卡西给出的的16位小数,打破了祖冲之位小数,打破了祖冲之的纪录的纪录 1579年年,韦达证明韦达证明373.14159265 353.14159265 1630年年,最后一位用古典方法求最后一位用古典方法求的人格林伯格也只求到了的人格林伯格也只求到了的第的第39位小数位小数 分析方法分析方法从十七世纪中叶起,人们开始用更先进的分析方法来求从十七世纪中叶起,人们开始用更先进的分析方法来求的近
25、似的近似值,其中应用的主要工具是收敛的无穷乘积和无穷级数,在本节值,其中应用的主要工具是收敛的无穷乘积和无穷级数,在本节中我们将介绍一些用此类方法求中我们将介绍一些用此类方法求近似值的实例。近似值的实例。067702.321201920543432122103516.341403940543432122 取取20 k取取10 k 1656年,沃里斯年,沃里斯(Wallis)证明证明 112212227656543432122kkkkk 在微积分中我们学过泰勒级数,其中有在微积分中我们学过泰勒级数,其中有12)1(53arctan12053 kxxxxxkkk),(x当当1 x121)1(513
26、1140 kkk 取取20 k189184.3411391513114 取取10 k232316.3211191513114 在中学数学中证明过下面的等式在中学数学中证明过下面的等式31arctan21arctan1arctan4 左边三个正方形组成的矩形左边三个正方形组成的矩形中,中,由由 和和 可得可得CBA DC 和和 的展开式的收的展开式的收敛速度都比敛速度都比 快得多快得多21arctan1arctan31arctanACBD 麦琴麦琴(Machin)给出给出2391arctan51arctan44 (Machin公式公式)51arctan 44 2391tan 记记 ,得,得此式求
27、得了此式求得了的第的第100位小数且全部正确位小数且全部正确 其它方法其它方法除用古典方法与分析方法求除用古典方法与分析方法求的近似值以外,还有人用其他方法的近似值以外,还有人用其他方法来求来求的近似值。这里我们将介绍两种方法:的近似值。这里我们将介绍两种方法:概率方法概率方法 数值积分方法数值积分方法 概率方法概率方法取一个二维数组(取一个二维数组(x,yx,y),取一个充分大的正整),取一个充分大的正整 数数n n,重复,重复n n次,每次,每次独立地从次独立地从 (0 0,1 1)中随机地取一对)中随机地取一对 数数x x和和y y,分别检验,分别检验x2+y21x2+y21是否成立。是否成立。设设n n次试验中等式成立的共有次试验中等式成立的共有m m次,令次,令4m/n4m/n。但这种方法很难得到但这种方法很难得到的较好的近似值。的较好的近似值。数值积分方法数值积分方法 10214dxx 102114dxx 还可用其它数值积分公式来求还可用其它数值积分公式来求,但用此类方法,但用此类方法效果也很难做效果也很难做得比用幂级数展开更好得比用幂级数展开更好