常微分方程第五章微分方程建模案例

《常微分方程第五章微分方程建模案例》由会员分享，可在线阅读，更多相关《常微分方程第五章微分方程建模案例（14页珍藏版）》请在金锄头文库上搜索。

1、209第五章 微分方程建模案例微分方程作为数学科学的中心学科，已经有三百多年的发展历史，其解法和理论已日臻完善，可以为分析和求得方程的解（或数值解）提供足够的方法，使得微分方程模型具有极大的普遍性、有效性和非常丰富的数学内涵。微分方程建模包括常微分方程建模、偏微分方程建模、差分方程建模及其各种类型的方程组建模。微分方程建模对于许多实际问题的解决是一种极有效的数学手段，对于现实世界的变化，人们关注的往往是其变化速度、加速度以及所处位置随时间的发展规律，其规律一般可以用微分方程或方程组表示，微分方程建模适用的领域比较广，涉及到生活中的诸多行业，其中的连续模型适用于常微分方程和偏微分方程及其方程组建

2、模，离散模型适用于差分方程及其方程组建模。本章主要介绍几个简单的用微分方程建立的模型，让读者一窥方程的应用。下面简要介绍利用方程知识建立数学模型的几种方法：1利用题目本身给出的或隐含的等量关系建立微分方程模型这就需要我们仔细分析题目，明确题意，找出其中的等量关系，建立数学模型。例如在光学里面，旋转抛物面能将放在焦点处的光源经镜面反射后成为平行光线，为了证明具有这一性质的曲线只有抛物线，我们就是利用了题目中隐含的条件入射角等于反射角来建立微分方程模型的。2从一些已知的基本定律或基本公式出发建立微分方程模型我们要熟悉一些常用的基本定律、基本公式。例如从几何观点看，曲线上某点的切线斜率即函数 在该点

3、的导数；力学中的牛顿第二运)(xy)(xy动定律： ，其中加速度 就是位移对时间的二阶导数，也是速度对时间maFa的一阶导数等等。从这些知识出发我们可以建立相应的微分方程模型。例如在动力学中，如何保证高空跳伞者的安全问题。对于高空下落的物体，我们可以利用牛顿第二运动定律建立其微分方程模型，设物体质量为 ，空气m阻力系数为 ，在速度不太大的情况下，空气阻力近似与速度的平方成正比；k210设时刻 时物体的下落速度为 ，初始条件： . 由牛顿第二运动定律建立t v0)(v其微分方程模型： 2kvmgdt求解模型可得： )12(expmkgtv由上式可知，当 时，物体具有极限速度：t，kgvtli1其

4、中，阻力系数 ， 为与物体形状有关的常数， 为介质密度，s 为物sk体在地面上的投影面积。根据极限速度求解式子，在 一定时，要求落地,m速度 不是很大时，我们可以确定出 来，从而设计出保证跳伞者安全的降落1vs伞的直径大小来。3利用导数的定义建立微分方程模型导数是微积分中的一个重要概念，其定义为，xyxffxf 00lim)(lim)(商式 表示单位自变量的改变量对应的函数改变量，就是函数的瞬时平均变化xy率，因而其极限值就是函数的变化率。函数在某点的导数，就是函数在该点的变化率。由于一切事物都在不停地发展变化，变化就必然有变化率，也就是变化率是普遍存在的，因而导数也是普遍存在的。这就很容易将

5、导数与实际联系起来，建立描述研究对象变化规律的微分方程模型。例如在考古学中，为了测定某种文物的绝对年龄，我们可以考察其中的放射性物质（如镭、铀等） ，已经证明其裂变速度（单位时间裂变的质量，即其变化率）与其存余量成正比。我们假设时刻 时该放射性物质的存余量 是 的函t Rt数，由裂变规律，我们可以建立微分方程模型:211kRdt期中 是一正的比例常数，与放射性物质本身有关。求解该模型，我们解得：k，其中 c 是由初始条件确定的常数。从这个关系式出发，我们就可以tCeR测定某文物的绝对年龄。 （参考碳定年代法）另外，在经济学领域中，导数概念有着广泛的应用，将各种函数的导函数（即函数变化率）称为该

6、函数的边际函数，从而得到经济学中的边际分析理论。4利用微元法建立微分方程模型一般的，如果某一实际问题中所求的变量 符合下列条件： 是与一个变pp量 的变化区间 有关的量； 对于区间 具有可加性；部分量 的近t,bap,bai似值可表示为 。那么就可以考虑利用微元法来建立微分方程模型，其itf)(步骤是：首先根据问题的具体情况，选取一个变量例如 为自变量，并确定其变t化区间 ；在区间 中随便选取一个任意小的区间并记作 ，求出,ba,ba dt,相应于这个区间的部分量 的近似值。如果 能近似的标示为 上的一个ppba连续函数在 处的值 与 的乘积，我们就把 称为量 的微元且记作t)(tfddtf)

7、(p.这样，我们就可以建立起该问题的微分方程模型：dp.tfp)(对于比较简单的模型，两边积分就可以求解该模型。例如在几何上求曲线的弧长、平面图形的面积、旋转曲面的面积、旋转体体积、空间立体体积；代数方面求近似值以及流体混合问题；物理上求变力做功、压力、平均值、静力矩与重心；这些问题都可以先建立他们的微分方程模型，然后求解其模型。5熟悉一些经典的微分方程模型，对一些类似的问题，经过稍加改进或直接套用这些模型。 多年来，在各种领域里，人们已经建立起了一些经典的微分方程模型，熟212悉这些模型对我们是大有裨益的。案例 1 设警方对司机饮酒后驾车时血液中酒精含量的规定为不超过 %80.现有一起交通事

8、故,在事故发生 个小时后,测得司机血液中酒精含量)/(mlg 3是 ,又过两个小时后, 测得其酒精含量降为 ,试判断: %56 )/(40mlg事故发生时,司机是否违反了酒精含量的规定?解 模型建立设 为时刻 的血液中酒精的浓度, 则在时间间隔 内, 酒精浓度)(txt ,t的改变量 ,即tkxttx)()(其中 为比例常数, 式前负号表示浓度随时间的推移是递减的, 两边除以 , 0k t并令 , 则得到t ,dkxt且满足 以及 .40)5(,6)3(x0)(x模型求解容易求得通解为 , 代入 ,得到ktcte)(0)(x.ktxe则 为所求. 又由 代入 可得)0(x,4)5(,6)3(

9、0)(17.e0e250 kxkk将 代入得 80.17.0k 25.935617.0317.3 x故事故发生时,司机血液中的酒精浓度已超出规定.案例 2 在凌晨 1 时警察发现一具尸体, 测得尸体温度是 , 当时环境Co温度是 .一小时后尸体温度下降到 ,若人的正常体温是 ,估计死Co1Co2737者的死亡时间.解 运用牛顿冷却定律 ,得到它的通解为T)(out,0outtTte这里 是当 时尸体的温度,也就是所求的死亡时间时尸体的温度,将题目提0Tt213供的参数代入 27e)137(29)(t解得和 则 ，68et6)1(t34进一步得.)(09.2)(,27.0hLnt这时求得的 是死

10、者从死亡起到尸体被发现所经历的时间, 因此反推回去可t推测死者的死亡时间大约是前一天的夜晚 .35:1案例 3 建立铅球掷远模型不考虑阻力，设铅球初速度为 ，出手高度为v，出手角度为 （与地面夹角） ，建立投掷距离与 的关系式，并求 一h,hvh,定的条件下求最佳出手角度解 在图 5-1 坐标下铅球运动方程为， ， ， ，0xgy0)(xhy)(， cos)(vsinv解出 ， 后，可以得铅球掷远为t. 图 5-1cos)2si(sin122 vghgR这个关系还可表为 tanco2Rv由此计算 ，得最佳出手角度和最佳成绩分别为：0d， )(2sin1ghv ghv2设 ，则 ， mvh/0,

11、5.14.m4.1R案例 4 在一种溶液中，化学物质 分解而形成 ，其速度与未转换的 的ABA浓度成比例转换 的一半用了 分钟，把 的浓度 表示为时间的函数，并A20y作出图象解 记 的浓度为时间 的函数 ， 的浓度为 Bt)(ty)(tx一、假设O R xyh 2141 分解后产生 molAnmolB2容体的体积在反应过程中不变二、建立模型，求解有假设知， 的消耗速度与 的浓度成比例，故有下列方程成立A，kxtd其中 为比例系数k设反应开始时 ， 的浓度为 ，由题中条件知当 （分）时，0tA0 20t的浓度为 解初值问题A21)(x0)(dxkt得 ,ktte)(0它应满足.0201)(xx

12、k解得 ,2ln0k所以得 .)2ln0e)(（txt由于 的浓度为 浓度减少量的 倍，故有Bx.)e1(e)( 2ln002ln00 ttxxnty三、作图（如图 5-2）O tynx0215图 5-2案例 5 车间空气清洁问题某生产车间内有一台机器不断排出 ，为了清洁车间里的空气，用一台2CO鼓风机通入新鲜空气来降低车间空气中的 含量，那么，上述做法的清洁效果如何呢？这一问题是利用平衡原理来建模，即建立其微分方程模型请注意，平衡原理在建立微分方程模型时常表现为区间 上的微元形式：某个量在该,x区间上的增加量等于该区间段内进入量与迁出量的差解 1问题分析与假设 上述清洁空气的原理是通过鼓风机

13、通入新鲜的空气，其 含量尽管也有2CO但较低新鲜空气与车间内空气混合后再由鼓风机排出室外，从而降低 含2量为讨论问题方便，假设通入的新鲜空气能与原空气迅速均匀混合，并以相同风量排出车间.此问题中的主要变量及参数设为：车间体积： （单位：立方米） ，V时间： （单位：分钟） ，t机器产生 速度： （单位：立方米/分钟） ，2COr鼓风机风量： （单位：立方米/分钟）K新鲜空气中 含量： ，2%m开始时刻车间空气中 含量： ，2xt 时刻车间空气中 含量： .CO)(t2模型建立 考虑时间区间 ，并利用质量守恒定律： 内车间空气中,t,t含量的“增加”等于 时间内，通入的新鲜空气中 的量加上机器2

14、CO,t 2CO产生的 的量减去鼓风机排出的 的量，即2CO增加量=新鲜空气中含有 +机器产生的 排出的2 量 2量 2量数学上表示出来就是216.tdsKxrtKmtxtxV)%(%)(其中 . 于是令 ，取极限便得0t0.)(,0xtbadt其中 .,10VKbrma3模型求解与分析 此问题是一阶线性非齐次常微分方程的初值问题. 解之得 ,exp)10(10exp)( tVKrKmrKbtabtx 这就是 时刻车间空气中含 的百分比.显然， 否则 含量只t 2CO,02CO能增加. 令 则有,10)(limKrmrKtxt 这说明了，车间空气中 的含量最多只能降到 .由此可见，鼓风2CO%

15、0机风量越大（ 越大） ，新鲜空气中 含量越低（ 越小） ，净化效果越好.K24模型的优缺点分析及改进方向：优点：模型简洁，易于分析和理解，并体现了建立微分方程模型的基本思想，而且所得到的结果与常识基本一致.缺点：建立数学模型时所作出的假设过于简单.改进方向：（1） 考虑新鲜空气和车间内的空气的混合扩散过程重新建模；（2）若要使得车间空气中的 含量达到一定的指标，确定最优的实施方2CO案.案例 6 某人的食量是 （焦/天），其中 （焦/天）用于基本的新104675038陈代谢（即自动消耗）。在健身训练中，他所消耗的热量大约是 （焦/公斤69217天）乘以他的体重（公斤）。假设以脂肪形式贮藏的热量 地有效，而 公%101斤脂肪含热量 （焦）。试研究此人的体重随时间变化的规律。4186解 模型分析 在问题中并未出现“变化率”、“导数”这样的关键词，但要寻找的是体重（记为 ）关于时间 的函数。如果我们把体重 看作是时间 的连续可微函Wt Wt数，我们就能找到一个含有的 微分方程.dt模型假设 1.以 表示