微分方程建模案例.doc

第五章 微分方程建模案例 微分方程作为数学科学的中心学科，已经有三百多年的发展历史，其解法和理论已日臻完善，可以为分析和求得方程的解（或数值解）提供足够的方法，使得微分方程模型具有极大的普遍性、有效性和非常丰富的数学内涵。微分方程建模包括常微分方程建模、偏微分方程建模、差分方程建模及其各种类型的方程组建模。微分方程建模对于许多实际问题的解决是一种极有效的数学手段，对于现实世界的变化，人们关注的往往是其变化速度、加速度以及所处位置随时间的发展规律，其规律一般可以用微分方程或方程组表示，微分方程建模适用的领域比较广，涉及到生活中的诸多行业，其中的连续模型适用于常微分方程和偏微分方程及其方程组建模，离散模型适用于差分方程及其方程组建模。本章主要介绍几个简单的用微分方程建立的模型，让读者一窥方程的应用。下面简要介绍利用方程知识建立数学模型的几种方法：1利用题目本身给出的或隐含的等量关系建立微分方程模型这就需要我们仔细分析题目，明确题意，找出其中的等量关系，建立数学模型。 例如在光学里面，旋转抛物面能将放在焦点处的光源经镜面反射后成为平行光线，为了证明具有这一性质的曲线只有抛物线，我们就是利用了题目中隐含的条件入射角等于反射角来建立微分方程模型的。2从一些已知的基本定律或基本公式出发建立微分方程模型我们要熟悉一些常用的基本定律、基本公式。例如从几何观点看，曲线上某点的切线斜率即函数在该点的导数；力学中的牛顿第二运动定律：，其中加速度就是位移对时间的二阶导数，也是速度对时间的一阶导数等等。从这些知识出发我们可以建立相应的微分方程模型。 例如在动力学中，如何保证高空跳伞者的安全问题。对于高空下落的物体，我们可以利用牛顿第二运动定律建立其微分方程模型，设物体质量为，空气阻力系数为，在速度不太大的情况下，空气阻力近似与速度的平方成正比；设时刻时物体的下落速度为，初始条件：. 由牛顿第二运动定律建立其微分方程模型：求解模型可得：由上式可知，当时，物体具有极限速度：，其中，阻力系数，为与物体形状有关的常数，为介质密度，s为物体在地面上的投影面积。根据极限速度求解式子，在一定时，要求落地速度不是很大时，我们可以确定出来，从而设计出保证跳伞者安全的降落伞的直径大小来。3利用导数的定义建立微分方程模型导数是微积分中的一个重要概念，其定义为，商式表示单位自变量的改变量对应的函数改变量，就是函数的瞬时平均变化率，因而其极限值就是函数的变化率。函数在某点的导数，就是函数在该点的变化率。由于一切事物都在不停地发展变化，变化就必然有变化率，也就是变化率是普遍存在的，因而导数也是普遍存在的。这就很容易将导数与实际联系起来，建立描述研究对象变化规律的微分方程模型。 例如在考古学中，为了测定某种文物的绝对年龄，我们可以考察其中的放射性物质（如镭、铀等），已经证明其裂变速度（单位时间裂变的质量，即其变化率）与其存余量成正比。我们假设时刻时该放射性物质的存余量是的函数，由裂变规律，我们可以建立微分方程模型:期中是一正的比例常数，与放射性物质本身有关。求解该模型，我们解得：，其中c是由初始条件确定的常数。从这个关系式出发，我们就可以测定某文物的绝对年龄。（参考碳定年代法）另外，在经济学领域中，导数概念有着广泛的应用，将各种函数的导函数（即函数变化率）称为该函数的边际函数，从而得到经济学中的边际分析理论。4利用微元法建立微分方程模型一般的，如果某一实际问题中所求的变量符合下列条件：是与一个变量的变化区间有关的量；对于区间具有可加性；部分量的近似值可表示为。那么就可以考虑利用微元法来建立微分方程模型，其步骤是：首先根据问题的具体情况，选取一个变量例如为自变量，并确定其变化区间；在区间中随便选取一个任意小的区间并记作，求出相应于这个区间的部分量的近似值。如果能近似的标示为上的一个连续函数在处的值与的乘积，我们就把称为量的微元且记作.这样，我们就可以建立起该问题的微分方程模型：.对于比较简单的模型，两边积分就可以求解该模型。 例如在几何上求曲线的弧长、平面图形的面积、旋转曲面的面积、旋转体体积、空间立体体积；代数方面求近似值以及流体混合问题；物理上求变力做功、压力、平均值、静力矩与重心；这些问题都可以先建立他们的微分方程模型，然后求解其模型。5熟悉一些经典的微分方程模型，对一些类似的问题，经过稍加改进或直接套用这些模型。 多年来，在各种领域里，人们已经建立起了一些经典的微分方程模型，熟悉这些模型对我们是大有裨益的。案例1 设警方对司机饮酒后驾车时血液中酒精含量的规定为不超过.现有一起交通事故,在事故发生个小时后,测得司机血液中酒精含量是,又过两个小时后, 测得其酒精含量降为,试判断: 事故发生时,司机是否违反了酒精含量的规定? 解 模型建立设为时刻的血液中酒精的浓度, 则在时间间隔内, 酒精浓度的改变量,即其中为比例常数, 式前负号表示浓度随时间的推移是递减的, 两边除以, 并令, 则得到且满足以及.模型求解 容易求得通解为, 代入,得到.则为所求. 又由代入可得将代入得 80.故事故发生时,司机血液中的酒精浓度已超出规定.案例2 在凌晨1时警察发现一具尸体, 测得尸体温度是, 当时环境温度是.一小时后尸体温度下降到,若人的正常体温是,估计死者的死亡时间.解 运用牛顿冷却定律,得到它的通解为,这里是当时尸体的温度,也就是所求的死亡时间时尸体的温度,将题目提供的参数代入解得 和 则，进一步得.这时求得的是死者从死亡起到尸体被发现所经历的时间, 因此反推回去可推测死者的死亡时间大约是前一天的夜晚.案例3 建立铅球掷远模型不考虑阻力，设铅球初速度为，出手高度为，出手角度为（与地面夹角），建立投掷距离与的关系式，并求一定的条件下求最佳出手角度ORxyha 解 在图5-1坐标下铅球运动方程为， ，解出，后，可以得铅球掷远为 . 图5-1这个关系还可表为 由此计算，得最佳出手角度和最佳成绩分别为：， 设，则，案例4 在一种溶液中，化学物质分解而形成，其速度与未转换的的浓度成比例转换的一半用了分钟，把的浓度表示为时间的函数，并作出图象 解 记的浓度为时间的函数，的浓度为 一、假设 1分解后产生 2容体的体积在反应过程中不变 二、建立模型，求解 有假设知，的消耗速度与的浓度成比例，故有下列方程成立，其中为比例系数 设反应开始时，的浓度为，由题中条件知当（分）时，的浓度为解初值问题得 ,它应满足.解得 ,所以得 .由于的浓度为浓度减少量的倍，故有Otynx0.三、作图（如图5-2） 图5-2案例5 车间空气清洁问题某生产车间内有一台机器不断排出，为了清洁车间里的空气，用一台鼓风机通入新鲜空气来降低车间空气中的含量，那么，上述做法的清洁效果如何呢？这一问题是利用平衡原理来建模，即建立其微分方程模型请注意，平衡原理在建立微分方程模型时常表现为区间上的微元形式：某个量在该区间上的增加量等于该区间段内进入量与迁出量的差解 1问题分析与假设 上述清洁空气的原理是通过鼓风机通入新鲜的空气，其含量尽管也有但较低新鲜空气与车间内空气混合后再由鼓风机排出室外，从而降低含量为讨论问题方便，假设通入的新鲜空气能与原空气迅速均匀混合，并以相同风量排出车间.此问题中的主要变量及参数设为：车间体积：（单位：立方米），时间：（单位：分钟），机器产生速度：（单位：立方米/分钟），鼓风机风量：（单位：立方米/分钟）新鲜空气中含量：，开始时刻车间空气中含量：，t时刻车间空气中含量：.2模型建立 考虑时间区间，并利用质量守恒定律：内车间空气中含量的“增加”等于时间内，通入的新鲜空气中的量加上机器产生的的量减去鼓风机排出的的量，即增加量=新鲜空气中含有+机器产生的排出的数学上表示出来就是.其中. 于是令，取极限便得其中3模型求解与分析 此问题是一阶线性非齐次常微分方程的初值问题. 解之得这就是时刻车间空气中含的百分比.显然，否则含量只能增加. 令则有这说明了，车间空气中的含量最多只能降到.由此可见，鼓风机风量越大（越大），新鲜空气中含量越低（越小），净化效果越好.4模型的优缺点分析及改进方向：优点：模型简洁，易于分析和理解，并体现了建立微分方程模型的基本思想，而且所得到的结果与常识基本一致.缺点：建立数学模型时所作出的假设过于简单.改进方向：（1） 考虑新鲜空气和车间内的空气的混合扩散过程重新建模；（2）若要使得车间空气中的含量达到一定的指标，确定最优的实施方案.案例6 某人的食量是（焦/天），其中（焦/天）用于基本的新陈代谢（即自动消耗）。在健身训练中，他所消耗的热量大约是（焦/公斤天）乘以他的体重（公斤）。假设以脂肪形式贮藏的热量地有效，而公斤脂肪含热量（焦）。试研究此人的体重随时间变化的规律。解 模型分析 在问题中并未出现“变化率”、“导数”这样的关键词，但要寻找的是体重（记为）关于时间的函数。如果我们把体重看作是时间的连续可微函数，我们就能找到一个含有的微分方程.模型假设 1.以表示时刻某人的体重，并设一天开始时人的体重为.2体重的变化是一个渐变的过程。因此可认为是关于连续而且充分光滑的.3体重的变化等于输入与输出之差，其中输入是指扣除了基本新陈代谢之后的净食量吸收；输出就是进行健身训练时的消耗.模型建立 问题中所涉及的时间仅仅是“每天”，由此，对于“每天”体重的变化=输入-输出.由于考虑的是体重随时间的变化情况，因此，可得体重的变化/天=输入/天输出/天.代入具体的数值，得 输入/天 = （焦/天）（焦/天）= （焦/天）， 输出/天 = （焦/公斤天）（公斤）= （焦/天）.体重的变化/天=（公斤/天）.考虑单位的匹配，利用 “公斤/天=”,可建立如下微分方程模型.模型求解 用变量分离法求解，模型方程等价于，积分得，从而求得模型解就描述了此人的体重随时间变化的规律.模型讨论 现在我们再来考虑一下：此人的体重会达到平衡吗?显然由的表达式，当时，体重有稳定值.我们也可以直接由模型方程来回答这个问题。在平衡状态下，是不发生变化的，所以.这就非常直接地给出了.至此，问题已基本上得以解决.案例7、人口预测模型 由于资源的有限性,当今世界各国都注意有计划地控制人口的增长,为了得到人口预测模型,必须首先搞清影响人口增长的因素,而影响人口增长的因素很多,如人口的自然出生率、人口的自然死亡率、人口的迁移、自然灾害、战争等诸多因素,如果一开始就把所有因素都考虑进去，则无从下手.因此,先把问题简化,建立比较粗糙的模型,再逐步修改,得到较完善的模型. 模型1（ 马尔萨斯 （Malthus） 模型） 英国人口统计学家马尔萨斯（17661834）在担任牧师期间,查看了教堂100多年人口出生统计资料,发现人口出生率是一个常数,于1789年在人口原理一书中提出了闻名于世的马尔萨斯人口模型,他的基本假设是：在人口自然增长过程中,净相对增长（出生率与死亡率之差）是常数,即单位时间内人口的增长量与人口成正比,比例系数设为,在此假设下,推导并求解人口随时间变化的数学模型.解 设时刻的人口为,把当作连续、可微函数处理（因人口总数很大,可近似地这样处理,此乃离散变量连续化处理）,据马尔萨斯的假设,在到时间段内,人口的增长量为,并设时刻的人口为,于是 这就是马尔萨斯人口模型,用分离变量法易求出其解为,此式表明人口以指数规律随时间无限增长. 模型检验：据估计1961年地球上的人口总数为,而在以后7年中,人口总数以每年2%的速度增长,这样, ,于是 .这个公式非常准确地反映了在17001961年间世界人口总数.因为,这期间地球上的人口大约每35年翻一番,而上式断定34.6年增加一倍（请读者证明这一点）但是,后来人们以美国人口为例,用马尔萨斯模型计算结果与人口资料比较,却发现有很大的差异,尤其是在用此模型预测较遥远的未来地球人口总数时,发现更令人不可思议的问题,如按此模型计算,到2670年,地球上将有36 000亿人口.如果地球表面全是陆地（事实上,地球表面还有80%被水覆盖）,我们也只得互相踩着肩膀站成两层了,这是非常荒谬的,因此,这一模型应该修改. 模型2（逻辑Logistic模型） 马尔萨斯模型为什么不能预测未来的人口呢？这主要是地球上的各种资源只能供一定数量的人生活,随着人口的增加,自然资源环境条件等因素对人口增长的限制作用越来越显著,如果当人口较少时,人口的自然增长率可以看作常数的话,那么当人口增加到一定数量以后,这个增长率就要随人口的增加而减小.因此,应对马尔萨斯模型中关于净增长率为常数的假设进行修改.1838年,荷兰生物数学家韦尔侯斯特(Verhulst)引入常数,用来表示自然环境条件所能容许的最大人口数（一般说来,一个国家工业化程度越高,它的生活空间就越大,食物就越多,从而就越大）,并假设将增长率等于,即净增长率随着的增加而减小,当时,净增长率趋于零,按此假定建立人口预测模型.解 由韦尔侯斯特假定,马尔萨斯模型应改为上式就是逻辑模型,该方程可分离变量,其解为,.下面,我们对模型作一简要分析.（1）当,即无论人口的初值如何,人口总数趋向于极限值；（2）当时,这说明是时间的单调递增函数；（3）由于,所以当时,单增；当时,单减,即人口增长率由增变减,在处最大,也就是说在人口总数达到极限值一半以前是加速生长期,过这一点后