实验六-上海交通大学数学系.ppt

改进Euler法的几何解释 若 的解为x(t) , 代入方程从tk→tk+1 积分导出 O tk f t f =f (x(t),t) tk+1 Euler法 xk+1- xk=τf (xk,tk) 改进Euler法 xk+1- xk=τ[f (xk,tk)+ 舰的实际过程 仿真方法 θ1 =arctan[(H - vwτ)/veτ] P0 M0 模仿真实事件行为和过程 在这个问题上，就是一步步地模拟导弹追踪敌 在t = 0 时,敌艇在M0(0,H), 导弹在原点P0指向M0 M1 在t =τ时,敌艇的位置为 M1(veτ, H)， 导弹的位置为 P1 (x1,y1), 而 x1=0, y1= vwτ, P1 θ1 导弹飞行方向的倾角为 x2= x1+ vwτcosθ1 y2= y1+ vwτsinθ1 P1 P0 M0 M1 θ2= arctan[(H-y2)/(2veτ-x2)] θ2 在t =2τ时,敌艇的位置 为M2(2τve, H)， 位置为P2(x2, y2), 由θ1表达式可写出cosθ1和 sinθ1的表达式. 此时导弹飞行方向的倾角为 导弹的 M2 P2 θ1 依此类推… 仿真迭代格式 在t=(k+1)τ时，导弹位于Pk+1(xk+1, yk+1) 其中 未用微分方程得到同样结果! 例如τ = 0.05 使用Matlab的计算结果 （m4_4） k tk xk yk 1 0.05 0.000 00 22.500 00 2 0.10 1.037 36 44.976 07 3 0.15 3.412 05 67.350 41 4 0.20 7.646 15 89.448 43 5 0.25 14.867 90 110.757 96 6 0.30 29.194 80 128.107 02 仍可能不能改善结果，因此可以采取变步长法： 2.在对由模型直接得到的微分方程组用Euler法计 算时，我们计算到ykH,yk+1≥H即停止，但这样 的做法可能会有不小误差，即使步长减小而结果 实验任务 参见教材 当计算到ykHyk+1，两端距离H较远时，以xk, yk 为起点,缩小步长继续计算，将改善结果. 试用这个方法改进授课中计算的结果. 何地击中敌舰? 3.如果当基地发射导弹的同时,敌艇立即由仪器发 觉. 假定敌艇为一高速快艇，它即刻以135 km/h 的速度与导弹方向垂直的方向逃逸,问导弹何时 实验任务 参见教材 任务5是成其它角度 补充任务 若猎犬位于一个半径300m的圆形场地中央， 追踪问题有时也以猎犬追兔形式出现， 逃逸，问猎犬需要多久才能捕获兔子 速度追逐，同时兔子以 650m/min 的速度 发现兔子在场地边缘，即以 800m/min （首先尝试建立数学模型，再用数值方法 求解;其次考虑仿真方法）. 谢谢各位！ * 学习Euler,是认识数学的最好途径。 ― 高斯 纯粹数学与应用数学是理解世界及其发展的一把主要钥匙。 ― 里约热内卢宣言 数学实验 导弹跟踪问题 上海交通大学数学科学学院 实验目的 微分方程数学建模 回顾微分方程解析解法 介绍Euler方法 介绍改进Euler方法 介绍仿真方法 实际问题 某军的一导弹基地发现正北方向120 km 处海面上有敌艇一艘以90 km/h的速度向正东方 向行驶. 该基地立即发射导弹跟踪追击敌艇, 导弹速度为450 km/h，自动导航系统使导弹在 任一时刻都能对准敌艇.试问导弹在何时何处 击中敌艇？ （追踪问题） 120 P(x, y) 当t =0时,导弹位于原点O, 敌艇位于(0,H)点, H =120 (km). 当时刻t ,导弹位于P(x(t),y(t)) 敌艇位于(90t, H)点 O 建立坐标系 x 时刻t 导弹的位置为P(x(t),y(t)),其 Vw = 450 (km/h) Ve = 90 (km/h) 数学模型 速度可由水平分速度与垂直分速度合成 导弹方向指向敌艇,故有 导出一阶微分方程组 解析方法 技巧性较强: 消去 t ,化为二阶方程 连同初始条件 设 降阶得到一阶可分离变量方程 导弹击中艇时y=H, 得到此时其位置和时刻 进而解出 Euler 方法十分简单, 就是用差商代替微商 数值方法 (Euler方法) 第k+1步的值之间的关系式 即采用 通常取Δt为常数τ,就得到由第 k 步的值到 伟大的欧拉（Leonhard?Euler，1707～1783 ） ? 著作辉煌：涉及几乎每一数学领域 ? 数学天才：约翰·伯努利指引 13岁进巴塞尔大学 ? 毅力惊人：失明17年，完成书和约400篇论文