第七讲-微分方程模型(Ⅱ)课件.ppt

1、第三章第三章 微分方程建模（微分方程建模（） 3.4 战争模型战争模型 3.5 饿狼追兔问题饿狼追兔问题 3.6 放射性废物的处理问题放射性废物的处理问题 战争分类：正规战争，游击战争，混合战争战争分类：正规战争，游击战争，混合战争兵力兵力因因战斗战斗及及非战斗减员非战斗减员而减少，因而减少，因增援增援而增加而增加战斗力与射击次数及命中率有关战斗力与射击次数及命中率有关建模思路和方法为用数学模型讨论社会建模思路和方法为用数学模型讨论社会领域的实际问题提供了可借鉴的示例领域的实际问题提供了可借鉴的示例3.4战争模型战争模型第一次世界大战第一次世界大战Lanchester提出预测战役结局的模型提出

2、预测战役结局的模型只考虑双方兵力多少和战斗力强弱只考虑双方兵力多少和战斗力强弱0),(),()(0),(),()(tvyyxgtytuxyxftx3.4.1一般战争模型一般战争模型1）每方战斗减员率取决于双方的兵力和战斗力，）每方战斗减员率取决于双方的兵力和战斗力， f(x,y)甲方战斗减员率，甲方战斗减员率，g(x,y)乙方战斗减员率乙方战斗减员率2）每方非战斗减员率与本方兵力成正比，比）每方非战斗减员率与本方兵力成正比，比例系数分别为例系数分别为3）甲乙双方的增援率分别为）甲乙双方的增援率分别为u(t), v(t)f, g 取决于战争类型取决于战争类型x(t) 甲方兵力，甲方兵力，y(t)

3、 乙方兵力乙方兵力模模型型假假设设模型模型建立建立（1）,)2()(-)(-tvybxytuxayx3.4.2 正规战争模型正规战争模型甲甲(乙乙)方战斗减员率只取决于乙方战斗减员率只取决于乙(甲甲)方的兵力和战斗力方的兵力和战斗力双方均以正规部队作战双方均以正规部队作战xxprbbxg,忽略非战斗减员忽略非战斗减员 假设没有增援假设没有增援)3()0(,)0(00yyxxbxyayxf(x, y)=ay, a 乙方每个士兵的杀伤率乙方每个士兵的杀伤率a=ry py, ry 射击率，射击率， py 命中率命中率)(ty)(tx0ak0k0kbk0k正规战争模型正规战争模型为判断战争的结局，不求

4、为判断战争的结局，不求x(t), y(t)而在相平面上讨论而在相平面上讨论 x 与与 y 的关系的关系00)0(,)0(yyxxbxyayxaybxdxdy2020bxaykkbxay22000yxk时平方律平方律 模型模型甲方胜 0k平局0kyyxxprprabxy200乙方胜乙方胜 上式说明双方初始兵力之比上式说明双方初始兵力之比y0/x0以平方关系影响以平方关系影响着战争的结局，如乙方兵力增加到原来的着战争的结局，如乙方兵力增加到原来的2倍（甲倍（甲方不变），则影响战争结局的能力增加到方不变），则影响战争结局的能力增加到4倍，或倍，或者若甲方的战斗力如射击率者若甲方的战斗力如射击率ry增

5、加到原来的增加到原来的4倍（倍（px, rx, py不变），乙方只要将初始兵力不变），乙方只要将初始兵力y0增加到原来的增加到原来的2倍就可抗衡倍就可抗衡平方律平方律 模型模型yyxxprprabxy200结果结果分析分析乙方胜乙方胜3.4.3 游击战争模型游击战争模型双方都用游击部队作战双方都用游击部队作战甲方战斗减员率还随着甲方兵力的增加而增加甲方战斗减员率还随着甲方兵力的增加而增加忽略非战斗减员忽略非战斗减员假设没有增援假设没有增援yrxxxxssrprddxyyxg/,),(00)0(,)0(yyxxdxyycxyxf(x, y)=cxy, c 乙方每个士兵的杀伤率乙方每个士兵的杀伤率

6、c = ry pyry射击率射击率py 命中率命中率py=sry /sxsx 甲方活动面积甲方活动面积sry 乙方射击有效面积乙方射击有效面积)(tycm0dm)(tx0m0m0m游击战争模型游击战争模型00)0(,)0(yyxxdxyycxyx00dxcymmdxcy乙方胜时000yxmyryyxrxxssrssrcdxy00线性律线性律 模型模型甲方胜 0m平局 0mcddxdyyryyxrxxssrssrcdxy00线性律线性律 模型模型结果结果分析分析乙方胜乙方胜即初始兵力之比即初始兵力之比y0/x0以线性关系影响战争结局，并以线性关系影响战争结局，并且当射击率和射击有效面积一定时，增

7、加活动面且当射击率和射击有效面积一定时，增加活动面积积Sx与增加初始兵力与增加初始兵力y0起着同样的作用起着同样的作用)(ty)(tx0乙方胜, 0n平局, 0n甲方胜, 0n00)0(,)0(yyxxbxycxyx3.4.4混合战争模型混合战争模型甲方为游击部队，乙方为正规部队甲方为游击部队，乙方为正规部队020222bxcynnbxcy02002cxbxy乙方胜0n100)/(200 xy02002xsrsprxyryyxxx乙方必须乙方必须10倍于甲方的兵力方可取胜倍于甲方的兵力方可取胜.设设 x0=100, rx/ry=1/2, px=0.1, sx=1(km2), sry=1(m2)

8、正规战模正规战模型检验型检验 Engel J.H. 利用二次大战末期美日硫黄岛利用二次大战末期美日硫黄岛战役时美军的战地记录验证了正规战争模战役时美军的战地记录验证了正规战争模型型 美军美军1945年年2月月19日开始进攻硫黄岛，战斗日开始进攻硫黄岛，战斗36天，日军天，日军21500人全部阵亡或被俘，战地记录全部人全部阵亡或被俘，战地记录全部遗失，美军投入兵力遗失，美军投入兵力73000人，伤亡人，伤亡20265人人 美国曾用这个模型分析越南战争美国正规部队美国曾用这个模型分析越南战争美国正规部队（乙方）要取胜越南游击部队（甲方），至少要投入（乙方）要取胜越南游击部队（甲方），至少要投入8倍

9、于游击部队一方的兵力，而美国最多只能派出倍于游击部队一方的兵力，而美国最多只能派出6倍于倍于越南的兵力，美国得出不可取胜的结论，最终撤出了越南的兵力，美国得出不可取胜的结论，最终撤出了越南越南 美军有增援，日军没有增援美军有增援，日军没有增援.,0=（忽略非战斗减员） , 0tv),(,tJayxf),(,tbAyxg 其它065130003260001054000ttttu用用A(t)和和J(t)表示美军和日军第表示美军和日军第t t天的人数，天的人数，在模型（在模型（1 1）中取）中取且且0),(),()(0),(),()(tvyyxgtytuxyxftx（1） )13(215000, 0

10、0)()(JAtbAdtdJtutaJdtdA由由u(t)及每天的伤亡及每天的伤亡人数可算出人数可算出A(t)，t=1,2, ,36（见图见图3.4.4 中虚线中虚线）于是有于是有图图3.4.4 )14(011ttuJaAtA )15(01tAbJtJ 3612037000A由由A(t)的实际数据可得的实际数据可得 )13(215000, 00)()(JAtbAdtdJtutaJdtdA对式（对式（13）用求和代替积）用求和代替积分得分得,0106. 0203700021500b因因J(36)=0, J(0)=21500, 式式（15）中令中令t t=36求得求得再将再将b=0.0106代入代

11、入(15)式，式，又可算出又可算出J(t), t=1,2, ,36在式（在式（1414）中令）中令t=36，得，得 36136136JAua（16） )17(0544. 011ttuJtA代入式（代入式（14）得）得,0544. 037250020265a其中分子是美军的总伤亡数，为其中分子是美军的总伤亡数，为20265人，分母可人，分母可由已算出的由已算出的J(t)得到，为得到，为372500人，于是从（人，于是从（16）式有式有图图3.4.4由式（由式（17）就可算出）就可算出A(t)的理论值，图的理论值，图3.5.4中用中用实线画出，与实际值（虚线）相比，可以看出吻实线画出，与实际值（虚

12、线）相比，可以看出吻合的情况合的情况3.5 饿狼追兔问题饿狼追兔问题现有一只兔子，一匹狼，兔子位于狼的正西100米处假设兔子与狼同时发现对方并一起起跑，兔子往正北60米处的巢穴跑，而狼在追兔子已知兔子、狼是匀速跑而且狼的速度是兔子的两倍，问题是兔子能否安全回到巢穴问题问题)(xfy xyB),(yxChO)0 ,100(A兔子在O处，狼在A处由于狼要盯着兔子追，所以狼行走的是一条曲线，且在同一时刻，曲线上狼的位置与兔子的位置的连线为曲线上该点处的切线设狼的行走轨迹为y=f(x)， 首先建立坐标系如图3.5.160图3.5.10|0|100100 xxyydxxfhxfxyhx1002)(12)

13、(0则有又因狼的速度是兔子的两倍，所以在相同时间内，狼走的距离为兔子走的距离的两倍假设在某时刻，兔子跑到(0,h)处,而狼在C(x,y)处，则有)(xfy xyB),(yxChO)0 ,100(A60dxxfhxfxyhx1002)(12)(01002)(1)( 2xdxxfxxfy两边求导并整理) 1 (0| )(0| )()(1)(21001002 xxxfxfxfxf x)(xfy xB),(yxCO)0 ,100(A60 xyh320010301)(2123xxxf,603200)0(f解方程可得狼的行走轨迹为因所以狼追不上兔子若狼和兔子的行走路线不在坐标轴上，如何考虑？ 某些类型的导

14、弹对目标追击的数学模型与此数学模型相似 环境污染是人类面临的一大公害，放射性污染对人类生命安全和地球上的生物存在严重的威胁，所以特别为人们所关注和平利用原子能可为人类造福，但是核废物处置不好，将对人类是一大危害核废料如何处置为好，必须进行科学论证 过去一段时间，美国原子能委员会为了处理浓缩的放射性废物，他们把废物装入密封的圆桶，然后扔到水深为300英尺的海里 一些生态学家和科学家为此表示担心，圆桶是否会在运输过程中破裂而造成放射性污染？经过实验证明不会破裂. 圆桶扔到海洋中时是否因与海底碰撞而发生破裂？3.6 放射性废物的处理问题放射性废物的处理问题问题问题几位工程师进行了大量的实验以后发现：

15、当圆桶的速度超过40英尺/秒时，就会因为碰撞而破裂mg436.5272/2 .32秒英尺g如图3.6.1选取坐标系，圆桶重量，使圆桶向下，m圆桶质量， g重力加速度，建模与求解建模与求解B水作用在圆桶上浮力，B=470.327磅 D 水作用在圆桶上的阻力.图3.6.1为此我们需计算圆桶同海底碰撞时的速度，是否会超过40英尺/秒？D=0.08v磅秒/英尺（通过大量实验得出如下结论：圆桶方位对于阻力影响甚小，可以忽略不计）cvBF则作用在圆桶上的力为,maF 22dtyda 由牛顿第二定律：而所以) 1 (0)0()(122vcgcvBmdtyd（1）是一个二阶常微分方程 将 )2(1cgtect

16、其解为dtdvdtyd22代入（1）得0)0(vgvcgdtdv) 1 (0)0()(122vcgcvBmdtyd由(2)式可知，圆桶的速度为时间t的函数，要确定圆桶同海底的碰撞速度，就必须算出桶碰到海底所需的时间t但作为y的显函数求出t是困难的，所以不能用(2)来求圆桶同海底的碰撞速度 )2(1cgtect但由（2）我们可以得到圆桶的极限速度vTcecBvcgttT1 limcecBvcgttT1 lim Tvtv显然有如果极限速度小于40英尺/秒，那么圆桶就不可能因同海底碰撞而破裂然而秒英尺/86.713cBvT这个数值太大了，还不能断定v(t)究竟是否能超过40英尺/秒 .下面将速度v

17、作为位置y的函数v(y) 来考虑.由v(t)=vy(t),利用复合函数微分法 ) 3(00vcvdydvvg代入（1）中，得dydvvdtdydydvdtdv对于（3），利用微分方程数值解法，借助于计算机很容易求出v(300).下面用一个简便的方法得到v(300)的一个很好的近似值.) 1 (0)0()(122vcgcvBmdtyd秒）（英尺/7 .45436.527300109.572 .322300u在(3)中令c =0（即不考虑水的阻力），用u代替v得 ) 3(00vcvdydvvg )4(00udyduug积分，得)(2)(22Byguyug或由此，得)300(300u,300y ,3

18、00y磅7 . 330008. 0u1.无阻力时，圆桶的速度总会大一些，因此2. 当y增加时，速度v增加，所以对于有由此可以得出水作用在桶上的阻力D总是小于然而使圆桶向下的合力-B57.1磅.比D大得多，因而可忽略D 所以可以认为u(y)是的一个很好的近似值实际上，用数值解法可算出v(300)=45.1英尺/秒，与u(300)=45.7英尺/秒是比较接近的其理由是：u(300)就是v(300)的一个很好的近似值, 综上得出结论，圆桶能够因与海底碰撞而破裂工程师们的说法是正确的这一模型科学地论证了美国原子能委员会过去处理核废料的方法是错误的，从而改变了美国政府过去的错误做法现在美国原子能委员会条例明确禁止把低浓度的放射性废物抛到海里，改为在一些废弃的煤矿中修建放置核废料的深井这一模型为全世界其他国家处理核废料提供了经验教训 作业：第三章课后习题3 要求： 1、用白色A4纸单面打印。 2、论文题目用三号黑体字、一级标题用五号黑体字，并居中。论文中其他汉字一律采用五号宋体字。1.5倍行距。 3、建立微分方程并用Matlab软件求解。 再再 见见