实验六 - 上海交通大学数学科学学院.pptVIP

改进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的计算结果