数学建模微分方程的应用举例.doc

1、第八节 数学建模微分方程的应用举例微分方程在物理学、力学、经济学和管理科学等实际问题中具有广泛的应用，本节我们将集中讨论微分方程的实际应用，读者可从中感受到应用数学建模的理论和方法解决实际问题的魅力.内容分布衰变问题逻辑斯谛方程价格调整问题人才分配问题模型追迹问题内容要点：一、衰变问题镭、铀等放射性元素因不断放射出各种射线而逐渐减少其质量, 这种现象称为放射性物质的衰变. 根据实验得知, 衰变速度与现存物质的质量成正比, 求放射性元素在时刻 t 的质量.用 x 表示该放射性物质在时刻 t 的质量, 则 表示 x 在时刻 t 的衰变速度, 于是“衰dt变速度与现存的质量成正比”可表示为(8.1)

2、.dxkt这是一个以 x 为未知函数的一阶方程, 它就是放射性元素 衰变的数学模型, 其中是比例常数, 称为衰变常数, 因元素的不同而异. 方程右端的负号表示当时间 t 增加0k时, 质量 x 减少 .解方程(8.1)得通解 若已知当 时, 代入通解 中可得.ktCe0t0,xktxCe则可得到方程(8.1)特解0,ktCe 0(),ktxe它反映了某种放射性元素衰变的规律.注: 物理学中 , 我们称放射性物质从最初的质量到衰变为该质量自身的一半所花费的时间为半衰期, 不同物质的半衰期差别极大. 如铀的普通同位素( )的半衰期约为 50 亿年;238U通常的镭( )的半衰期是 1600 年.半

3、衰期是上述放射性物质的特征 , 然而半衰期却不依26Ra赖于该物质的初始量, 一克 衰变成半克所需要的时间与一吨 衰变成半吨所需要26Ra26Ra的时间同样都是 1600 年, 正是这种事实才构成了确定考古发现日期时使用的著名的碳-14测验的基础.二、 逻辑斯谛（Logistic）方程:逻辑斯谛方程是一种在许多领域有着广泛应用的数学模型, 下面我们借助树的增长来建立该模型.一棵小树刚栽下去的时候长得比较慢, 渐渐地, 小树长高了而且长得越来越快, 几年不见, 绿荫底下已经可乘凉了; 但长到某一高度后, 它的生长速度趋于稳定, 然后再慢慢降下来. 这一现象很具有普遍性. 现在我们来建立这种现象的

4、数学模型.如果假设树的生长速度与它目前的高度成正比, 则显然不符合两头尤其是后期的生长情形, 因为树不可能越长越快 ; 但如果假设树的生长速度正比于最大高度与目前高度的差, 则又明显不符合中间一段的生长过程. 折衷一下, 我们假定它的生长速度既与目前的高度,又与最大高度与目前高度之差成正比.设树生长的最大高度为 H(m), 在 t(年) 时的高度为 h(t), 则有(8.2)(dhtkt其中 是比例常数. 这个方程为 Logistic 方程. 它是可分离变量的一阶常数微分方程.0k下面来求解方程(8.2). 分离变量得 ,()dhktH两边积分 ,()t得 11ln,hktC或 12kHtHt

5、ee故所求通解为 2(),1ktHkHthtCe其中的 是正常数.120CHe函数 的图象称为 Logistic 曲线. 图 8-8-1 所示的是一条典型的 Logistic 曲线, 由于()ht它的形状, 一般也称为 S 曲线. 可以看到, 它基本符合我们描述的树的生长情形. 另外还可以算得 lim().thH这说明树的生长有一个限制, 因此也称为限制性增长模式（阻滞增长模型）.注: Logistic 的中文音译名是“逻辑斯谛”. “逻辑”在字典中的解释是“客观事物发展的规律性”, 因此许多现象本质上都符合这种 S 规律. 除了生物种群的繁殖外, 还有信息的传播、新技术的推广、传染病的扩散以

6、及某些商品的销售等. 例如流感的传染、在任其自然发展( 例如初期未引起人们注意 )的阶段, 可以设想它的速度既正比于得病的人数又正比于未传染到的人数. 开始时患病的人不多因而传染速度较慢; 但随着健康人与患者接触, 受传染的人越来越多, 传染的速度也越来越快; 最后, 传染速度自然而然地渐渐降低, 因为已经没有多少人可被传染了.下面举两个例子说明逻辑斯谛的应用.人口阻滞增长模型 1837 年, 荷兰生物学家 Verhulst 提出一个人口模型(8.3)0(),()dykbytt其中 的称为生命系数.,kb我们不详细讨论这个模型, 只提应用它预测世界人口数的两个有趣的结果.有生态学家估计 k 的

7、自然值是 0.029. 利用本世纪 60 年代世界人口年平均增长率为 2%以及 1965 年人口总数 33.4 亿这两个数据, 计算得 从而估计得:2,b(1)世界人口总数将趋于极限 107.6 亿.(2)到 2000 年时世界人口总数为 59.6 亿.后一个数字很接近 2000 年时的实际人口数, 世界人口在 1999 年刚进入 60 亿.新产品的推广模型 设有某种新产品要推向市场, t 时刻的销量为 由于产品性能良(),xt好, 每个产品都是一个宣传品, 因此, t 时刻产品销售的增长率 与 成正比, 同时, 考dt虑到产品销售存在一定的市场容量 N, 统计表明 与尚未购买该产品的潜在顾客

8、的数量xt也成正比 , 于是有()Nxt(8.4)()dxkt其中 k 为比例系数. 分离变量积分, 可以解得(8.5)()1kNttCe由 的图像知，当 时, 则有 即销量 单调增加. 当 时, dxt:*x0,dxt()xt*()2Nxt当 时, 当 时, 即当销量达到最大需求量20;t*()2t20;dt*()22;tN 的一半时, 产品最为畅销 , 当销量不足 N 一半时, 销售速度不断增大, 当销量超过一半时, 销售速度逐渐减少.国内外许多经济学家调查表明. 许多产品的销售曲线与公式(8.5)的曲线(逻辑斯谛曲线)十分接近. 根据对曲线性状的分析, 许多分析家认为, 在新产品推出的初

9、期, 应采用小批量生产并加强广告宣传, 而在产品用户达到 20%到 80%期间 , 产品应大批量生产; 在产品用户超过 80%时, 应适时转产, 可以达到最大的经济效益 .三、价格调整模型假设某种商品的价格变化主要服从市场供求关系. 一般情况下,商品供给量 S 是价格 P的单调递增函数, 商品需求量 Q 是价格 P 的单调递减函数, 为简单起见, 分别设该商品的供给函数与需求函数分别为(8.6)(),()SPabQP其中 均为常数, 且,ab0.当供给量与需求量相等时, 由(8.6)可得供求平衡时的价格 eaPb并称 为均衡价格.eP一般地说, 当某种商品供不应求, 即 时, 该商品价格要涨

10、, 当供大于求, 即SQ时, 该商品价格要落. 因此, 假设 t 时刻的价格 的变化率与超额需求量 成SQ ()Pt QS正比, 于是有方程 ()dPkSt其中 用来反映价格的调整速度.0,k将(8.6)代入方程, 可得(8.7)()edPt其中常数 方程(8.7)的通解为()0,bk()tetC假设初始价格 代入上式, 得 于是上述价格调整模型的解为0()P0,eP()()tet由于 知, 时, 说明随着时间不断推延, 实际价格 将逐渐趋近0t. ()Pt均衡价格 .eP四、人才分配问题模型每年大学毕业生中都要有一定比例的人员留在学校充实教师队伍, 其余人员将分配到国民经济其他部门从事经济和

11、管理工作. 设 t 年教师人数为 科学技术和管理人员数目1(),xt为 又设 1 个教员每年平均培养 个毕业生, 每年从教育、科技和经济管理岗位退休、2(),xt 死亡或调出人员的比率为 表示每年大学生毕业生中从事教师职业所占比率(01)于是有方程(01),(8.8)11dxxt(8.9)212()dt方程(8.8)有通解(8.10)()1txCe若设 则 于是得特解110(),x0,(8.11)1()txe将(8.11)代入(8.9)方程变为(8.12)1()20(tdxet 求解方程(8.12)得通解(8.13)1()02()t txCee若设 则 于是得特解220(),x2100,x(8

12、.14)211()000t txexe(8.11)式和(8.14)式分别表示在初始人数分别为 情况, 对应于 的取值, 在12()t 年教师队伍的人数和科技经济管理人员人数 . 从结果看出 , 如果取 即毕业生全部留1在教育界, 则当 时, 由于 必有 而 说明教师队伍将迅t,1()xt2()0,xt速增加. 而科技和经济管理队伍不断萎缩, 势必要影响经济发展, 反过来也会影响教育的发展. 如果将 接近于零. 则 同时也导致 说明如果不保证适当比例的1()0xt2(),t毕业生充实教师选择好比率 , 将关系到两支队伍的建设, 以及整个国民经济建设的大局.五、追迹问题设开始时甲、乙水平距离为 1

13、 单位, 乙从 A 点沿垂直于 OA 的直线以等速 向正北行0v走; 甲从乙的左侧 O 点出发, 始终对准乙以 的速度追赶. 求追迹曲线方程, 并问0(1)nv乙行多远时, 被甲追到.建立如图 8-8-2 所示的坐标系, 设所求追迹曲线方程为 经过时刻 t, 甲在追迹().yx曲线上的点为 乙在点 于是有()Pxy0(1,).Bvt(8.15)tan,yx由题设, 曲线的弧长 OP 为 2001,xydnvt解出 代入(8.15), 得0vt 20().xxyyn两边对 x 求导, 整理得 21().这就是追迹问题的数学模型.这是一个不显含 y 的可降阶的方程, 设 , 代入方程得(),ypx或 21()xpn2,(1)dxn两边积分, 得 2 11l()ln|l|,xC即 21.np将初始条件 代入上式, 得 于是0|xxy 1.C(8.16)2,nyx两边同乘 并化简得21,y(8.17)21,nyx(8.16)与(8.17)式相加, 得11,2nnyx两边积分, 得 112()().nnyxxC 代入初始条件 得 故所求追迹曲线方程为0|x2,C112()()(1),nnxy n甲追到乙时, 即曲线上点 P 的横坐标 此时 即乙行走至离 A 点 个单1,2.y21n位距离时被甲追到.