导弹跟踪问题
摘要
本文研究的是导弹跟踪敌艇问题。通过建立合理的数学模型,利用matlab软件模拟导弹跟踪敌艇的过程,并进行分析计算最终求得导弹击中敌艇的时间和位置。
首先我们根据问题中的条件粗略描绘出导弹和敌艇的运动轨迹并对其进行分析,利用在任何时刻导弹的飞行方向指向敌机的位置得出导弹的运动轨迹满足微分方程组的结论。然后结合初值条件,并经过严格的数学公式推导和合理的假设,求解出导弹运动轨迹的方程。在matlab中分别采用Euler法、改进的Euler法和仿真方法近似计算导弹轨迹上的一系列点,接着对这些点进行精确度分析,绘制出导弹运动轨迹的曲线,并比较这三种数值方法的精确度。
最后我们对该模型进行分析评价,指出该模型的优点和不足。
关键词 导弹跟踪 运动轨迹 matlab 常微分方程 Euler法
正文
1. 模型问题详述
某军的一导弹基地发现正北方向120km处海面上有敌艇一艘以90km/h的速度向正东方向行驶。该基地立即发射导弹跟踪追击敌艇,导弹的速度为450km/h,自动导航系统使导弹在任一时刻都能对准敌艇。试问导弹在何时何处击中敌艇?并绘制导弹轨迹曲线图。 2. 模型假设
假设一:导弹与敌艇的速率恒定。
假设二:导弹飞行的轨迹切线方向始终指向敌艇。
假设三:导弹飞行的轨迹和敌艇行驶的高度始终在同一平面内。 假设四:导弹与敌艇的长度可以忽略,均可看成物理质点。 假设五:外界对导弹和敌艇的运动没有影响
3. 数据字典
(x,y) 坐标位置 t 经历的时间 ve 敌艇的速度 vw 导弹飞行的速度
H 敌艇最初点与导弹的距离
4. 问题分类及分析
问题一 导弹在何时何处击中敌艇?
分析:由于敌艇的运行轨迹和导弹的运行轨迹是处于同一平面内,故可建立起平面直角坐标
系。又由于导弹飞行方向始终指向敌机,即导弹飞行方向随时间的改变而改变,故可建立起微分方程并求解。
问题二 绘制导弹轨迹曲线图。
分析:由于导弹运动轨迹满足微分方程,所以我们可以利用Euler法、改进Euler法、仿真方法分别对所建立的数学模型进行数值求解,并绘制出导弹运动轨迹曲线图像。 5. 建模
设坐标系如下,取导弹基地为原点0(0,0)。x轴指向正东方,y轴指向正北方。
当t=0时,导弹位于O,敌艇位于点(0,H)(H=120,(km))设导弹t时刻的位置为P(x(t),y(t)),由题意,
其中
(1)
vw?450(km/h)。
。由于导弹轨迹的切vt,H),其中v=90(km/h)
ee另外在t时刻,敌艇位置应该为M(线方向必须指向敌艇,即直线PM的方向就是导弹轨迹上点P的切线方向,故有
dyH?y? (2) dxvet?xdydxH?y?() (3) dtdtve?x
方程(3)初值条件想
x(0)=0,y(0)=0 (4) 构成了一个关于时间变量t的一阶微分方程组的初值问题。 由(2)得
dx(H?y)?vet?x dy两边对t求导得
xdydx?dy?dx
(H?y)??????2dy?dt?vedtdydt即有
d2xdy?H?y??ve 2dydt把(1)写为dydt?d2v?dx???dy????e代入上式,就得到轨迹方程。这是一个二阶非线性微分方程,
?12加上初值条件,则初值问题
?2x?H?y??ve?d22?dyvw?dx???1????? dy??????x?1?y?0??dx??dyy?0???上式分别为(5),(6),(7)。
就是导弹的轨迹的数学模型。 方程(5)可以降阶。令p?dx,记??ve,则式(5)化为一介可分离变量方程 dyvwdp?dy?p2?1H?y
易得
?H?y?由式(7)得C?H??????C?p??p2??1? ?,从而
p?p2?1???H??H?y?????
于是有
??????HH?y??? (8) 1?
p????????2?H?y?????H????于是积分又可以得到
?1?x?2???H?y???1?H???1?H?H?y???1??1??? ??C1??利用式(6)得
C1??H1??2,于是导弹轨迹方程为
? x?1?2???H?y???1?H???1?H???H?y?1??1??????H (9) ?1?2??设导弹击中敌艇于B(L,H),以y=H代入(9)得
L?而导弹击中敌艇的时刻
T??H1??Le2?Hvevwvw?ve2222 (10)
vvw?ve?Hvw (11)
将数据代入(10),(11)式,得
L=25(km), T?0.2778(h)
接下来用数值方法求解 (1)Euler法
由方程(1),(3)解出时导弹轨迹上点的坐标为
?x,y?,则Euler格式为
kkdxdy和的表达式,取时间步长?t?t,对应tk?kt dtdt?vw??xk?1?xk?2??H?yk??1????????vetk?xk? ???vw???yk?1yk2??vetk?xk???1?????H?y?k????x?0,y?00?0?上式分别为(12)(13)(14) 当计算到
yk?H,yk?1?H即停止,于是,
相关推荐: