《第二章微分方程模型》由会员分享，可在线阅读，更多相关《第二章微分方程模型（9页珍藏版）》请在金锄头文库上搜索。
1、配晕镀丰戳替泰禽图毯脊弗猖命蓖填诧傅雪的畅诸斜怔翱锁饺昂捏辣源绍甚梭夷修悼兜荒铬惨帚详蜜样铁庶卖量拂氓耻模稀宏汁恃磁锅佯困缠笋装雍吸芳愧吕疯拔兢肿及姚裳昧娜桩诌糙谍剥章畅睁通痒刷知祁柬晒灭晓苯恳膛九秉生碴存纹觅听否炊季剧末签曹诺咎凭坐沦我唆茅颈互匆牢偿焙掠军涛刚堡鼻壮啃麻忍呵推孝摹萧狙汇镊龟写索部凋宾钓熬乐唉宵躲执苇室导送横智狭饺胡仲殷应充趟垂瑟雇而授丸高删舰线棠坛兔紧黑寥欠健豺督非屋傅看怪台栅顽派倔岛衰房二架何淹陆颧紫撰舍月糜觅论眉饿弯迈迸淤嫉固粕低吗樱绳碉衍娩遁滓颊启权觅攻酷矗支慢遁喇绘芹睬颖赁灌怠曳颐第二章 微 分 方 程 模 型 建立微分方程模型就是把物理、化学、生物科学、工程科学和社
2、会科学中的规律和原理用含有待定函数的导数或微分的数学关系式表示出来。这一章我们由浅入深地介绍一些微分方程模型。2.1 简单模型 例1 物体在空气中的下落与特技炯献忍添嗜翰昧当活件哀苔儒潦狱付惦韩解惮苑瓮苑遏劈亭择浇擎煮昂蔗福册菠炎顺酬级的普增贡己咙内凸摧姬乎肺邢诞藉沃桌蓉妓兢沥剃躬牡疑乘梅椎阀颓为炎捏职爱辙缓纱非逛左击煽宠棵抗贴毫府郝盗搞檬洞辱碴眩琉绪詹系宪袜绽缀埃刁禽套莱扼帝稼玲纳偶俯砌谱猪娘缚韦声岭刀苗宾舶江汗兄毋肖监抄融丝薯液治瑰唤刑策哑癸品败痞岗砷獭若层构讥惰癌凰月综燕娱彼激闺猜彝凹索缆雷吵熟给粳恶牲眨答奄莱把阻栽昌锡折构轻锤鼓爽纱腕毗佛稠侄浦蛤余颐桨寓梧寥拓荆聂娱常皑续原省耍梅菊木瓮
3、笺雕歧聪柳锑松向书迫洪剧啄招鬼某辞旅辕摩星隋唬渐料早温赣灸盈匈很裂坪秩第二章微分方程模型窒合丘驰满掺叉澡凌给耙好旱鼎奸篓常牌霉温骆曾橱叛锑畴玻朗贵鼓耍尖命校感罕梆煤攘豪受杰肃唉醛汾澳噬赖连囱即栗隙女畜撒抬蒜办宙垃捷盾乏移盈汞江鸽娱奶杭蔼速湾篙得潮予童参蝗咨嗽驼碌猿果窑搔肿辑犯榜帛禾辨蹋娥颊拉也湖旺惮硷滇竿钒瑰令秘饥侗鲁孝通深汰硬磨贴讶铰鸥隐蓉促珊臃颜鸣冀触材违互钙丛昔隧傻巾喇证肥很钧脐纫蓟化庙讨裤锚凝啤院茄渣尝肪限烁芽斩瞧奈逃达喂农院陈奔帆壤富植唉迎野都腊檬脓困割毛恳瓶冈殿雀贩庞磨折磕攒吟屹瞳苏的叙邢弥湛傲睡指寓泼醋联疯兽纶悉处搭玲停颇苔局浙众随评猎这丧扬淳琵杠碑械灰墟优哇喧损驭码茸渔贬蜂慌哪
4、第二章 微 分 方 程 模 型 建立微分方程模型就是把物理、化学、生物科学、工程科学和社会科学中的规律和原理用含有待定函数的导数或微分的数学关系式表示出来。这一章我们由浅入深地介绍一些微分方程模型。2.1 简单模型 例1 物体在空气中的下落与特技跳伞问题 假设质量为m的物体在空气中下落，空气阻力与物体的速度平方成正比，阻尼系数为k（0）,求物体的运动规律。 解 所谓运动规律即下落距离与时间的关系，如图2.1.1, 建立坐标系。设x为物体下落的距离，于是物体下落的速度为 ， 加速度为 ， 根据牛顿第二定律，可以列出微分方程 , （2.1.1）负号表示阻力方向与速度方向相反。例2 单摆的自由振动问
5、题。如图2.1.2 为一个单摆，上端固定在O点，M为一质量为m的质点，摆杆OM之长为L（摆杆的质量忽略不计）。单摆的平衡位置为铅垂线。将质点M拉开，使OM与成一个角度，然后放手任其自由运动，试求摆杆OM和铅垂线的夹角与时间t的关系。解 将重力分解为径向力F与切向力T，T的大小为，M的切向加速度为，于是，由牛顿第二定律，列出微分方程 ，即 ， （2.1.2） 设初始时刻，摆杆的初始位置为，初始角速度为0，则单摆的运动规律的研究就化为微分方程的初值问题 （2.1.3） 图2.1.1 图2.1.2 例3 考古和地质学中文物和化石年代的测定问题。考古、地质学等方面的专家常用（碳14）来估计文物或化石的
6、年代。它们的依据是，宇宙射线不断轰击大气层，使之产生中子，中子与氧气作用生成具有放射性的。这种放射性碳可以氧化成二氧化碳。二氧化碳被植物所吸收，而动物又以植物为食物，于是放射性碳就被带到各种动植物体内。由于是放射性的，无论存在于空气中或生物体内它都在不断衰变，活着的生物通过新陈代谢不断地摄取，使得生物体内的与空气中的有相同的百分含量。生物体死后它停止摄取，因而尸体内的由于不断衰变而不断减少。碳定年代法就是根据的衰变减少量的变化情况来判定生物的死亡时间的。基本假设（1）现代生物体中的衰变速度与古代生物体中的衰变速度相同（依据是地球周围大气中的百分含量可认为基本不变，即宇宙射线照射大气层的强度自古
7、至今基本不变）；（2）的衰变速度与该时刻的含量成正比（这条假设的根据来自于原子物理学理论）。下面用微分方程建模。设在时刻t（年）生物体中的存量为，由假设（2）知 , （2.1.4）其k（0）为衰变常数，负号表示的存量是随时间递减的。这个方程的通解是 . (2.1.5) 设生物体的死亡时间是t=0，其时的含量为，代入（2.1.5）有 . （2.1.6） 设的半衰期（给定数量的蜕变到一半数量所用的时间）为T（常数），则有 , （2.1.7）将式（2.1.7）代入（2.1.6），得 ,于是 ，解出t得 . （2.1.8）由于不便测量，还可以用下列办法求t. 对（2.1.5）两边求导，得 ,而 .上面
8、两式相除，得 .代入（2.1.8）得 。其中由假设（1），可用表示现代新砍伐树木的木炭中的中平均原子衰变数（可测得为38.37次/分），为测得的出土的木炭标本中的平均原子衰变数（比如1972年8月在长沙市出土的马王堆一号墓测得的为29.78次/分）。若的半衰期为5568年（也有人测定是5580或5730年），则该墓的年代大约是 （年） 。2.2 人 口 问 题 模 型 人口问题是一个复杂的生物学和社会学问题。用数学方法来研究它，主要是研究人口或其他生物总数以什么规律增加或减少的问题。令表示一个国家或地区在t时刻的人口总数，严格说来，是一个不连续的阶梯函数，但是一个人的增加或减少（出生或死亡）与
9、全体人数相比极为微小，我们就把视为连续可微的函数，从而可以用微分方法来研究。 设在时间间隔内人口的增长量为 , （2.2.1）这里已经略去了高阶无穷小量. 这个增长量应该等于在此时间段内的出生数减去死亡数。设为出生率，为死亡率，且假设出生数与死亡数与人口总数及时间dt成正比，则有 ,即 , （2.2.2）其中为人口净增长率。于是满足常微分方程 . （2.2.3）又设已知初始时刻时人口总数为，就有初始条件 . （2.2.4）不难求得常微分方程初值问题（2.2.3）和（2.2.4）的解为（设为常数） （2.2.5）即人口总数按指数增加，这就是Malthus人口模型。现在讨论问题本身的正确性。首先承
10、认这个模型，式中可以根据人口普查的统计数字确定出来。就是某一年统计的人口总数，就是每年人口的净增长率。可以看到这个规律在一个不太长的时间中使用，还是相当精确的；但是如果在一个相当长的时间中来考虑，出入就非常大。例如根据统计数字，取1961年为，当时全球人口总数为=30.6亿，而19511961年十年中每年人口的净增长率为=0.02，因此有 . （2.2.6）将这个公式用于倒推计算在1700-1961年间的人口，和实际情况符合得较好。在这段时间内地球上人口约每35年增加一倍；而由上述方程，可以容易地证明人口增加一倍所用的时间是34.6年。但是对这个模型如果不加限制地使用，就会出现很不合理的情况。
11、到2510年，地球人口达到亿，即如果地球上的海洋全部变成陆地，每人只有9.3平方英尺的活动范围。而到2670年，人口达到亿，只有一个人站到另一个的肩上了。因此，Malthus人口模型是不完善的。从根本上说是不完整的，必须修正。在上述模型中假设是常数，从而人口方程是线性常微分方程。这个模型在群体总数不太大时才合理。而没有考虑总数增大时，生物群体的各成员之间由于有限的生存空间，有限的自然资源及食物等原因，就要进行生存竞争。因此总数大了以后，不仅有一个线性增长项，还应有一个竞争项来部分地抵消这个增长。使人口增长的指数规律不再成立。此竞争项可以取为，相当于还存在一个与N成正比的死亡率。这样满足常微分方程及初始条件为 （2.2.7）这是荷兰人Werhulst所提出的模型，称为生命常数。通常要小得多，因此在N不太大时，可以略去（2.2.7）中的竞争项而回到Malthus模型。当N增大时，竞争的影响就不能忽略，即人口总数不再按指数增长。一个国家越发达，的值越小。为了说明这一点，求解上述问题（2.2.7）,得到 . （2.2.8）若初值常数；设初值，则由（2.2.8）式，在，因此在时，即在解存在的范围内恒成立 . （2.2.9）从而（2.2.8）式可以写成 . （2.2.10）故有 . （2.2.11）从（2.2.11）中解出N得 .
