7常微分方程在数学建模中的应用这里介绍几个典型的用微分方程建立数学模型的例子.一、人口预测模型由于资源的有限性,当今世界各国都注意有计划地控制人口的增长,为了得到人口预测模型,必须首先搞清影响人口增长的因素,而影响人口增长的因素很多,如人口的自然出生率、人口的自然死亡率、人口的迁移、自然灾害、战争等诸多因素,如果一开始就把所有因素都考虑进去,则无从下手.因此,先把问题简化,建立比较粗糙的模型,再逐步修改,得到较完善的模型.例1(马尔萨斯(Malthus)模型)英国人口统计学家马尔萨斯(1766—1834)在担任牧师期间,查看了教堂100多年人口出生统计资料,发现人口出生率是一个常数,于1789年在《人口原理》一书中提出了闻名于世的马尔萨斯人口模型,他的基本假设是:在人口自然增长过程中,净相对增长(出生率与死亡率之差)是常数,即单位时间内人口的增长量与人口成正比,比例系数设为r,在此假设下,推导并求解人口随时间变化的数学模型.解设时刻t的人口为)(tN,把)(tN当作连续、可微函数处理(因人口总数很大,可近似地这样处理,此乃离散变量连续化处理),据马尔萨斯的假设,在t到tt时间段内,人口的增长量为ttrNtNttN)()()(,并设0tt时刻的人口为0N,于是.,00)(ddNtNrNtN这就是马尔萨斯人口模型,用分离变量法易求出其解为)(00e)(ttrNtN,此式表明人口以指数规律随时间无限增长.模型检验:据估计1961年地球上的人口总数为91006.3,而在以后7年中,人口总数以每年2%的速度增长,这样19610t,901006.3N,02.0r,于是)1961(02.09e1006.3)(ttN.这个公式非常准确地反映了在1700—1961年间世界人口总数.因为,这期间地球上的人口大约每35年翻一番,而上式断定34.6年增加一倍(请读者证明这一点).但是,后来人们以美国人口为例,用马尔萨斯模型计算结果与人口资料比较,却发现有很大的差异,尤其是在用此模型预测较遥远的未来地球人口总数时,发现更令人不可思议的问题,如按此模型计算,到2670年,地球上将有36000亿人口.如果地球表面全是陆地(事实上,地球表面还有80%被水覆盖),我们也只得互相踩着肩膀站成两层了,这是非常荒谬的,因此,这一模型应该修改.例2(逻辑Logistic模型)马尔萨斯模型为什么不能预测未来的人口呢?这主要是地8球上的各种资源只能供一定数量的人生活,随着人口的增加,自然资源环境条件等因素对人口增长的限制作用越来越显著,如果当人口较少时,人口的自然增长率可以看作常数的话,那么当人口增加到一定数量以后,这个增长率就要随人口的增加而减小.因此,应对马尔萨斯模型中关于净增长率为常数的假设进行修改.1838年,荷兰生物数学家韦尔侯斯特(Verhulst)引入常数mN,用来表示自然环境条件所能容许的最大人口数(一般说来,一个国家工业化程度越高,它的生活空间就越大,食物就越多,从而mN就越大),并假设将增长率等于mNtNr)(1,即净增长率随着)(tN的增加而减小,当mNtN)(时,净增长率趋于零,按此假定建立人口预测模型.解由韦尔侯斯特假定,马尔萨斯模型应改为,,000)(1ddNtNNNNrtN上式就是逻辑模型,该方程可分离变量,其解为,)(00e11)(ttrmmNNNtN.下面,我们对模型作一简要分析.(1)当t,mNtN)(,即无论人口的初值如何,人口总数趋向于极限值mN;(2)当mNN0时,01ddNNNrtNm,这说明)(tN是时间t的单调递增函数;(3)由于NNNNNrtNmm211dd222,所以当2mNN时,0dd22tN,tNdd单增;当2mNN时,0dd22tN,tNdd单减,即人口增长率tNdd由增变减,在2mN处最大,也就是说在人口总数达到极限值一半以前是加速生长期,过这一点后,生长的速率逐渐变小,并且迟早会达到零,这是减速生长期;(4)用该模型检验美国从1790年到1950年的人口,发现模型计算的结果与实际人口在1930年以前都非常吻合,自从1930年以后,误差愈来愈大,一个明显的原因是在20世纪60年代美国的实际人口数已经突破了20世纪初所设的极限人口.由此可见该模型的缺点之一是mN不易确定,事实上,随着一个国家经济的腾飞,它所拥有的食物就越丰富,mN的值也就越大;(5)用逻辑模型来预测世界未来人口总数.某生物学家估计,029.0r,又当人口总数9为91006.3时,人口每年以2%的速率增长,由逻辑模型得mNNrtNN1dd1,即mN91006.31029.002.0,从而得91086.9mN,即世界人口总数极限值近100亿.值得说明的是:人也是一种生物,因此,上面关于人口模型的讨论,原则上也可以用于在自然环境下单一物种生存着的其他生物,如森林中的树木、池塘中的鱼等,逻辑模型有着广泛的应用.二、市场价格模型对于纯粹的市场经济来说,商品市场价格取决于市场供需之间的关系,市场价格能促使商品的供给与需求相等(这样的价格称为(静态)均衡价格).也就是说,如果不考虑商品价格形成的动态过程,那么商品的市场价格应能保证市场的供需平衡,但是,实际的市场价格不会恰好等于均衡价格,而且价格也不会是静态的,应是随时间不断变化的动态过程.例3试建立描述市场价格形成的动态过程的数学模型解假设在某一时刻t,商品的价格为)(tp,它与该商品的均衡价格间有差别,此时,存在供需差,此供需差促使价格变动.对新的价格,又有新的供需差,如此不断调节,就构成市场价格形成的动态过程,假设价格)(tp的变化率tpdd与需求和供给之差成正比,并记),(rpf为需求函数,)(pg为供给函数(r为参数),于是,,0)0(,ddpppgrpftp其中0p为商品在0t时刻的价格,为正常数.若设baprpf),(,dcppg)(,则上式变为,,0)0()()(ddppdbpcatp①其中dcba,,,均为正常数,其解为cadbcadbptptca)(0e)(.10下面对所得结果进行讨论:(1)设p为静态均衡价格,则其应满足0)(),(pgrpf,即dpcbpa,于是得cadbp,从而价格函数)(tp可写为ppptptca)(0e)()(,令t,取极限得ptpt)(lim这说明,市场价格逐步趋于均衡价格.又若初始价格pp0,则动态价格就维持在均衡价格p上,整个动态过程就化为静态过程;(2)由于tcacapptp)(0e)()(dd,所以,当pp0时,0ddtp,)(tp单调下降向p靠拢;当pp0时,0ddtp,)(tp单调增加向p靠拢.这说明:初始价格高于均衡价格时,动态价格就要逐步降低,且逐步靠近均衡价格;否则,动态价格就要逐步升高.因此,式①在一定程度上反映了价格影响需求与供给,而需求与供给反过来又影响价格的动态过程,并指出了动态价格逐步向均衡价格靠拢的变化趋势.三、混合溶液的数学模型例4设一容器内原有100L盐,内含有盐10kg,现以3L/min的速度注入质量浓度为0.01kg/L的淡盐水,同时以2L/min的速度抽出混合均匀的盐水,求容器内盐量变化的数学模型.解设t时刻容器内的盐量为)(txkg,考虑t到ttd时间内容器中盐的变化情况,在dt时间内容器中盐的改变量注入的盐水中所含盐量-抽出的盐水中所含盐量容器内盐的改变量为xd,注入的盐水中所含盐量为td301.0,t时刻容器内溶液的质量浓度为ttx)23(100)(,假设t到ttd时间内容器内溶液的质量浓度不变(事实上,容器内的溶液质量浓度时刻在变,由于td时间很短,可以这样看).于是抽出的盐水中所含盐量为tttxd2)23(100)(,这样即可列出方程11ttxtxd1002d03.0d,即txtx100203.0dd.又因为0t时,容器内有盐10kg,于是得该问题的数学模型为d20.03d100(0)10xxttx,,这是一阶非齐次线性方程的初值问题,其解为24)100(109)100(01.0)(tttx.下面对该问题进行一下简单的讨论,由上式不难发现:t时刻容器内溶液的质量浓度为34)100(10901.0100)()(tttxtp,且当t时,01.0)(tp,即长时间地进行上述稀释过程,容器内盐水的质量浓度将趋于注入溶液的质量浓度.溶液混合问题的更一般的提法是:设有一容器装有某种质量浓度的溶液,以流量1V注入质量浓度为1C的溶液(指同一种类溶液,只是质量浓度不同),假定溶液立即被搅匀,并以2V的流量流出这种混合溶液,试建立容器中质量浓度与时间的数学模型.首先设容器中溶质的质量为)(tx,原来的初始质量为0x,t=0时溶液的体积为2V,在dt时间内,容器内溶质的改变量等于流入溶质的数量减去流出溶质的数量,即tVCtVCxddd2211,其中1C是流入溶液的质量浓度,2C为t时刻容器中溶液的质量浓度,,tVVVxC)(2102于是,有混合溶液的数学模型11220dd(0)xCVCVtxx,.该模型不仅适用于液体的混合,而且还适用于讨论气体的混合.四、振动模型12振动是生活与工程中的常见现象.研究振动规律有着极其重要的意义.在自然界中,许多振动现象都可以抽象为下述振动问题.例5设有一个弹簧,它的上端固定,下端挂一个质量为m的物体,试研究其振动规律.解假设(1)物体的平衡位置位于坐标原点,并取x轴的正向铅直向下(见图4).物体的平衡位置指物体处于静止状态时的位置.此时,作用在物体上的重力与弹性力大小相等,方向相反;(2)在一定的初始位移0x及初始速度0v下,物体离开平衡位置,并在平衡位置附近作没有摇摆的上下振动;(3)物体在t时刻的位置坐标为)(txx,即t时刻物体偏离平衡位置的位移;(4)在振动过程中,受阻力作用.阻力的大小与物体速度成正比,阻力的方向总是与速度方向相反,因此阻力为txhdd,h为阻尼系数;(5)当质点有位移)(tx时,假设所受的弹簧恢复力是与位移成正比的,而恢复力的方向总是指向平衡位置,也就是总与偏离平衡位置的位移方向相反,因此所受弹簧恢复力为kx,其中k为劲度系数;(6)在振动过程中受外力)(tf的作用.在上述假设下,根据牛顿第二定律得)(dddd22xfkxtxhtxm,①这就是该物体的强迫振动方程.由于方程①中,)(tf的具体形式没有给出,所以,不能对式①直接求解.下面我们分四种情形对其进行讨论.1.无阻尼自由振动在这种情况下,假定物体在振动过程中,既无阻力、又不受外力作用.此时方程①变为0dd22kxtxm,令2mk,方程变为0dd222xtx,特征方程为022,特征根为i2,1,通解为tCtCxcossin21,或将其写为tCCCtCCCCCxcossin22212222112221xO图413ttAcossinsincos,)sin(tA其中2221CCA,22212sinCCC,22211cosCCC.这就是说,无阻尼自由振动的振幅2221CCA,频率mk均为常数.2.有阻尼自由振动在该种情况下,考虑物体所受到的阻力,不考虑物体所受的外力.此时,方程①变为0dddd22kxtxhtxm,令2mk,2mh,方程变为0dd2dd222xtxtx,特征方程为0222,特征根222,1.