数学建模真题2.ppt

Clear[a,b,t,x,y,k] h=0.05 For[x[0]=0;y[0]=0;t[0]=0;k=0,k=100,k++, Print[k]; x[k+1]=x[k]+450*h/Sqrt[1+((120-y[k])/(90*t[k]- x[k]))^2];Print[x[k]//N]; y[k+1]=y[k]+450*h/Sqrt[1+((90*t[k]-x[k])/(120- y[k]))^2];Print[y[k]//N]; t[k]=h*k; If[y[k-1]120y[k]=120,Break[]]] * * 导弹跟踪问题 数学实验 北京交大数学系 Euler,他是我们大家的 读读Euler,读读 老师. ─ P.S.Laplace 对自然界的深刻研 究是数学最丰富的源泉. ─ Joseph Fourier 问导弹在何时何处 实际问题 某军的一导弹基地发现正北方向 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 时刻t 导弹的位于为P(x(t),y(t)),其 Vw =450 (km/h) Ve =90 (km/h) 数学模型 速度可由水平分速度与垂直分速度合成 导弹方向指向敌艇,故有 导出一阶微分方程组 解析方法 有较强技巧性: 消去 t ,化为二阶方程 连同初始条件 设 降阶得到一阶可分离变量方程 导弹击中艇时y=H, 得到此时其位置和时刻 Euler方法十分简单, 就是用差商代替 Euler方法(数值方法) 的值到第k+1步的值之间的关系式 微商,即 通常取Δt为常数τ,就得到由第k步 t = tk+1时导弹位置为(xk+1 ，yk+1) 计算到yk H, yk+1≥H 停止，取L≈xk Euler 迭代格式 利用Mathematica 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 L≈x6=29.195 T ≈0.3244 τ = 0.005 τ = 0.05 k tk xk yk L≈x6=25.667 T ≈0.2852 舰的实际过程. 仿真方法 θ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