常微分方程在数学建模中地应用91806

1、实用标准文案常微分方程在数学建模中的应用这里介绍几个典型的用微分方程建立数学模型的例子一、人口预测模型由于资源的有限性，当今世界各国都注意有计划地控制人口的增长，为了得到人口预测模型，必须首先搞清影响人口增长的因素，而影响人口增长的因素很多，如人口的自然出生率、人口的自然死亡率、人口的迁移、自然灾害、战争等诸多因素，如果一开始就把所有因素都考虑进去，则无从下手.因此，先把问题简化，建立比较粗糙的模型，再逐步修改，得到较完善的 模型.例1 (马尔萨斯 (Malthus ) 模型) 英国人口统计学家马尔萨斯(1766 1834 ) 在担任牧师期间，查看了教堂100多年人口出生统计资料 ，发现人口出

2、生率是一个常数，于1789年在人口原理一书中提出了闻名于世的马尔萨斯人口模型，他的基本假设是：在人口自然增长过程中，净相对增长(出生率与死亡率之差)是常数，即单位时间内人口的增长量与人口成正比，比例系数设为r,在此假设下，推导并求解人口随时间变化的数学模型.解 设时刻t的人口为N(t)，把N(t)当作连续、可微函数处理(因人口总数很大，可近似地这样处理，此乃离散变量连续化处理)，据马尔萨斯的假设，在t到t t时间段内，人口的文档大全增长量为N(tt) N(t) rN(t) t,并设tto时刻的人口为No,于dNdtN(to)No.这就是马尔萨斯人口模型，用分离变量法易求出其解为N(t)Noer

3、(tto)此式表明人口以指数规律随时间无限增长.9 模型检验：据估计1961年地球上的人口总数为3.06 10，而在以后7年中，人口总数以每年2%的速度增长，这样心1961,$ 3.06 109 ,r 0.02，于是9 0.02(t 1961) N(t) 3.06 10 e ().这个公式非常准确地反映了在1700 1961年间世界人口总数.因为，这期间地球上的人口大约每35年翻一番，而上式断定34.6年增加一倍(请读者证明这一点).但是，后来人们以美国人口为例，用马尔萨斯模型计算结果与人口资料比较，却发现有很大的差异，尤其是在用此模型预测较遥远的未来地球人口总数时，发现更令人不可思议的问题如

4、按此模型计算，到2670年，地球上将有36 000亿人口 .如果地球表面全是陆地(事实上，地球表面还有 80%被水覆盖)，我们也只得互相踩着肩膀站成两层了，这是非常荒谬的，因此，这一模型应该修改.例2 (逻辑Logistic模型)马尔萨斯模型为什么不能预测未来的人口呢？这主要是地球上的各种资源只能供一定数量的人生活，随着人口的增加，自然资源环境条件等因素对人口增长的限制作用越来越显著，如果当人口较少时，人口的自然增长率可以看作常数的话，那么当人口增加到一定数量以后，这个增长率就要随人口的增加而减小.因此，应对马尔萨斯模型中关于净增长率为常数的假设进行修改.1838年,荷兰生物数学家韦尔侯斯特

5、(Verhulst)引入常数Nm用来表示自然环境条件所 能容许的最大人口数(一般说来，一个国家工业化程度越高，它的生活空间就越大，食物就越多， 从而Nm就越大)，并假设将增长率等于r 1 旦©，即净增长率随着 N(t)的增加而减小，Nm当N(t) Nm时，净增长率趋于零，按此假定建立人口预测模型 .解由韦尔侯斯特假定，马尔萨斯模型应改为dN 一 rdtN(to)1 -NN( No-N, o上式就是逻辑模型该方程可分离变量淇解为,N(t)1Nm瓦Nm1 e r(t t0)卜面，我们对模型作一简要分析,N(t) Nm,即无论人口的初值如何，人口总数趋向于极限值 Nm；(2)当 0 N2

6、2 dNNm时,前N0,这说明N(t)是时间t的单调递增函Nmr2 12NmN,所以当N5空c dN0 ,单增;dt2Nm 口 d N一时，-2-2 dt2c dN0 , dt单减，即人口增长率 dN由增变减，在Nm处最大dt，也就是说在人口总数达到极限值一半以前是加速生长期，过这一点后，生长的速率逐渐变小，并且迟早会达到零这是减速生长期;(4)用该模型检验美国从1790年到1950年的人口，发现模型计算的结果与实际人口在1930年以前都非常吻合 启从1930年以后，误差愈来愈大，一个明显的原因是在 20世纪60年代美国的实际人口数已经突破了20世纪初所设的极限人口.由此可见该模型的缺点之一是

7、Nm不易确定，事实上，随着一个国家经济的腾飞，它所拥有的食物就越丰富，Nm的值也就越大；(5)用逻辑模型来预测世界未来人口总数.某生物学家估计,r 0.029，又当人口总数为93.06 10时，人口每年以2%的速率增长，由逻辑模型得1 dN 彳 N r 1 -N dtNm即从而得0.02 0.029 193.06 109NT-_ _ _ 一 9Nm 9.86 10 ,即世界人口总数极限值近100亿.值得说明的是：人也是一种生物，因此，上面关于人口模型的讨论 ，原则上也可以用于在自 然环境下单一物种生存着的其他生物 ，如森林中的树木、池塘中的鱼等，逻辑模型有着广泛的 应用.二、市场价格模型对于纯

8、粹的市场经济来说，商品市场价格取决于市场供需之间的关系，市场价格能促使商品的供给与需求相等(这样的价格称为(静态)均衡价格).也就是说，如果不考虑商品价格形成 的动态过程，那么商品的市场价格应能保证市场的供需平衡，但是，实际的市场价格不会恰好等于均衡价格，而且价格也不会是静态的，应是随时间不断变化的动态过程.例3试建立描述市场价格形成的动态过程的数学模型解假设在某一时刻t,商品的价格为p，它与该商品的均衡价格间有差别，此时，存在供需差，此供需差促使价格变动.对新的价格，又有新的供需差，如此不断调节，就构成市场价格 形成的动态过程，假设彳格p(t)的变化率dp与需求和供给之差成正比，并记f(p,

9、r)为需求dt函数，g( p)为供给函数(r为参数)，于是dpf-f P,r g P ，dtP(0) P0 ，ap b , g( p) cp d ,则上式变为其中P0为商品在t 0时亥U的价格,为正常数.若设f (p,r)其中a,b,c,d均为正常数，其解为P(t)Po卜面对所得结果进行讨论(1)设p为静态均衡价格于是得pb d ，从而价格函数a c，取极限得这说明，市场价格逐步趋于均衡价格上，整个动态过程就化为静态过程;(2)由于所以，当p°p时，电dtdp dt(a c) p (b d),p(0) po ，(a c)t，则其应满足f(p,r) g(p)a p b cpp(t)可写

10、为p(t) (Po p)etlim p(t)0,(a c)t p ，.又若初始价格pop,则动态价格就维持在均衡价格pdpdt(p p°) (a c)e (a c)t0,p(t)单调下降向p靠拢；当p0p时,型 0,p(t)单调增加dt，且逐步靠近均衡价格；，而需求与向p靠拢这说明：初始价格高于均衡价格时，动态价格就要逐步降低 否则，动态价格就要逐步升高.因此，式在一定程度上反映了价格影响需求与供给 供给反过来又影响价格的动态过程，并指出了动态价格逐步向均衡价格靠拢的变化趋势三、混合溶液的数学模型例4 设一容器内原有 100L盐，内含有盐10kg,现以3L/min的速度注入质量浓度为

11、0.01kg/L的淡盐水，同时以2L/min的速度抽出混合均匀的盐水，求容器内盐量变化的数学模型.解 设t时刻容器内的盐量为x(t) kg,考虑t到t dt时间内容器中盐的变化情况，在dt时间内容器中盐的改变量注入的盐水中所含盐量一抽出的盐水中所含盐量容器内盐的改变量为dx,注入的盐水中所含盐量为0.01 3dt,t时刻容器内溶液的质量浓度为x，假设t至iJ t dt时间内容器内溶液的质量浓度不变(事实上，容器内的100 (3 2)t溶液质量浓度时刻在变，由于dt时间很短，可以这样看).于是抽出的盐水中所含盐量为100x2dt,这样即可列出方程(3 2)tdx 0.03dtdhdt，dx dt

12、0.032x100 t又因为t 0时,容器内有盐10kg,于是得该问题的数学模型为dx 2xdt 100 t0.03,x(0) 10,这是一阶非齐次线性方程的初值问题，其解为x(t) 0.01(100 t)9 1042(100 t)2卜面对该问题进行一下简单的讨论，由上式不难发现:t时刻容器内溶液的质量浓度为小 x(t)9 104p0.013,100 t(100 t)3且当t时，p（t）0.01，即长时间地进行上述稀释过程，容器内盐水的质量浓度将趋于注入溶液的质量浓度.溶液混合问题的更一般的提法是：设有一容器装有某种质量浓度的溶液，以流量Vi注入质量浓度为Ci的溶液（指同一种类溶液，只是质量浓

13、度不同），假定溶液立即被搅匀，并以V2的流量流出这种混合溶液，试建立容器中质量浓度与时间的数学模型首先设容器中溶质的质量为x（t）,原来的初始质量为Xo ,t =0时溶液的体积为V2,在d t时间内，容器内溶质的改变量等于流入溶质的数量减去流出溶质的数量，即dx C1V1dt C2V2dt ,其中C1是流入溶液的质量浓度,C2为t时刻容器中溶液的质量浓度,C2Vo(ViV2)t'于是，有混合溶液的数学模型dx 一 一C1V1 C2V2 , dtx(0)x。.该模型不仅适用于液体的混合，而且还适用于讨论气体的混合四、振动模型.在自然界中，许多振动是生活与工程中的常见现象.研究振动规律有着

14、极其重要的意义 振动现象都可以抽象为下述振动问题例5设有一个弹簧，它的上端固定，下端挂一个质量为 m的物体，试研究其振动规律.解 假设（1 ）物体的平衡位置位于坐标原点 ，并取X轴的正向铅直向下（见图 4）.物体的平衡位置指物体处于静止状态时的位置.此时，作用在物体上的重力与弹性力大小相等方向相反；（2）在一定的初始位移 xo及初始速度v0下，物体离开平衡位置，并在平衡位置附近作没有摇摆的上下振动；(3)物体在t时刻的位置坐标为 x x(t),即t时刻物体偏离平衡位置的位移；(4)在振动过程中，受阻力作用.阻力的大小与物体速度成正比，阻力的方向总是与速度方向相反，因此阻力为 h0,h为阻尼系数

15、；(5)当质点有位移x(t)时，假设所受的 dt弹簧恢复力是与位移成正比的，而恢复力的方向总是指向平衡位置，也就是总与偏离平衡位置的位移方向相反，因此所受弹簧恢复力为kx,其中k为劲度系数；(6)在振动过程中受外力f(t)的作用.在上述假设下，根据牛顿第二定律得d2xdxm2hkx f(x),dtdt这就是该物体的强迫振动方程.由于方程中，f(t)的具体形式没有给出，所以，不能对式直接求解.下面我们分四种情形对其进行讨论.1 .无阻尼自由振动在这种，f#况下，假定物体在振动过程中，既无阻力、又不受外力作用.此时方程变为d2xm 2 kx 0 , dt2A k2令一，方程变为md2x 22x 0

16、,dt特征方程为220,特征根为1,2i ,通解为xC1sin tC2 cos t或将其写为x . Ci2 C；_Ci：c；一sinClC2 cos t ,C； C|A cos sinsincosA sin( t )22其中 A . CiC2 ,sinC2C2,cosCl,C12C2这就是说，无阻尼自由振动的振幅C c22濒率JK均为常数.m2 .有阻尼自由振动在该种，f#况下，考虑物体所受到的阻力,不考虑物体所受的外力.此时方程变为d2x m -dt2hdXdtkx0,k 2 h令一 ，2 ，方程变为m md2x 小 dx22 -dt dt特征方程为 2 22 0,特征根i20,22 .根据

17、与 的关系，又分为如下三种情形:(1)大阻尼情形>.特征根为二不等实根，通解为xC1e(22)tC2e(22)t(2)临界阻尼情形.特征根为重根，通解为x (C1 C2t)e这两种情形，由于阻尼比较大，都不发生振动.当有一初始扰动以后，质点慢慢回到平衡位置，位移随时间t的变化规律分别如图5和图6所示.图5Xoo t(3)小阻尼情形，.特征根为共轲复根，通解为e t(C1sin . 22tC2 sin 22t)将其简化为Ae t sin( . 22t其中 A C12 C22 ,sinC222CiC2Ci ,c0s2Ci一，振幅A2C 2t随时间t的增加而减小.因此，这是一种衰减振动.位移随时间t的变化规律见图7.3.无阻尼强迫振动在这种，f#形下，设物体不受阻力作用，其所受外力为x简谐力f(t) msin pt,此时方程化为d2xm-rdtkxmsin pt,d2x dt2sin pt ,O根据ip是否等于特征根i淇通解分为如下两种情形:(1)当p时，其通解为x2 sin pt C1sinpt C 2 cos t ,此时，特解的振幅象;时，其通解为122为常数，但当p接近于时，将会导致振幅增大，发生类似共振的现p1八八t cos pt C 1sint C 2 cos t ,2p此时，特解的振幅1t随时间t的增加而增大，这种现象称