数模(微分方程模型).pptVIP

微分方程模型 重庆邮电大学 数理学院 沈世云 三 放射性核废料处理问题 此模型曾用于分析越南战争，乙方为美国. 根据类似于以上的计算及上一世纪四五十年代发生在马来亚、菲律宾、老挝等地的混合战争的实际情况估计出，正规部队一方要取胜至少要投入8倍于游击部队的兵力，而美国最多只能派出6倍于越南的兵力，越南战争最后以美国失败告终. 例 硫磺岛战役 硫磺岛是太平洋上一座由火山熔岩冷却后形成的火山岛，地形起伏，沟壑纵横，熔洞密布，悬崖峭壁临海高耸。从高空俯瞰，20余平方公里大小的硫磺岛像一只砍去双腿，又被拔光毛的火鸡，火鸡头位于岛的西南端，高度168米的摺钵山雄踞鸡头，是全岛制高点，一个鸡嘴状的岬角一直伸进滚滚波涛中。北部从鸡背一直到东北部鸡尾部分，是一片错落起伏的高地，由一系列小山岗和陡峭的峡谷构成。小山岗高程大多百米左右，地形复杂，可伏重兵。南部鸡脖子和鸡胸部位，地势低平，有一小片被梯状台地逼住的海滩，勉强可作登陆场。除此以外，全岛没有任何可供船舶停靠的锚地或港湾。 硫磺岛虽是弹丸小岛，却处在战略要津。它正当东京与美军新占领的塞班岛之间，距二地各约1200公里。美军占领塞班岛以后，一直以塞班岛为基地空袭东京。但因硫磺岛的报警作用，美军对东京的空袭一直效果不佳。驻硫磺岛的日军战斗机还不时升空拦截，冲散美国机群。为总攻日本，美军势必要夺占硫磺岛。而为东京安全，日军也势必要死守硫磺岛。结果，这座杳无人迹的小小火山岛，在太平洋战争后期，就成为日美必争之地。 美军在控制菲律宾后，于1945年1月3日开始，对硫磺岛实施轰炸。 2月16日开始硫磺岛战役，美军出动舰艇200多艘次，飞机400多架次，对硫磺岛日军3个机场和滩头阵地进行炮击和轰炸。 19日，美海军陆战第4师和第5师在600多架次飞机和舰炮火力掩护下，以250多艘登陆艇和500多辆水陆坦克、装甲车组成5个登陆波，在硫磺岛东南部登陆。经过两周激战，残余的3000多日军退进山洞死守。 3月8日，日军师团长栗原中将与800余名残兵在山洞内集体自杀。 此战役日军被击毙2.2万人，被俘1000人；美军为攻占该岛阵亡7000人，负伤1.9万人。 硫黄岛战役 J. H. Engel用二次大战末期美日硫黄岛战役中美军的战地记录，对正规战争模型进行了验证，发现模型结果与实际数据吻合很好. 设美军人数为A(t)，日军人数为J(t)，日方无援军，数学模型为： dA/dt = ?aJ(t) +u(t), dJ/dt = ?bA(t), (19) A(0) = 0, J(0) = 21500, 美军的战地记录给出增援率u(t)为 u(t) = 54000, 0 ? t 1, u(t) = 6000, 2 ? t 3, (20) u(t) = 13000, 5 ? t 6, u(t) = 0, 其它. 用A(k+1)?A(k)来代替dA/dt，有A(k+1) = A(k)?aJ(k) +u(k)，可得 A(t) = A(0)? a[J(0) +…+ J(t ? 1)] + [u(0) +…+u(t?1)], (21) 类似有 J(t) = J(0)?b[A(0) +…+A(t ?1)]. (22) 战役进行了36天，t = 36. J(36) = 0，各A(k)可由战地记录得到，这样由(22)式可得参数b，再利用求出的b代入(22)可得各J(k)，最后可由(21)得出a. 把求出的a代回(21)可得各A(t)的理论值， 微分方程的解的数值解法 Euler 折线法 初值问题的Euler折线法 Euler 折线法举例 Euler 折线法源程序 Euler折线法举例（续） Runge-Kutta 方法 Runge-Kutta 方法 四阶 R-K 方法源程序 Runge-Kutta 方法 Euler 法与 R-K法误差比较 Matlab 解初值问题 dsolve 求解析解 dsolve 的使用 dsolve 举例 dsolve 举例 dsolve 举例 Matlab函数数值求解 Matlab提供的ODE求解器 参数说明 数值求解举例 数值求解举例 Matlab 求解微分方程小结 例 2：求微分方程 在初值条件 下的特解，并画出解函数的图形。 y=dsolve(x*Dy+y-exp(x)=0,y(1)=2*exp(1),x) ezp