实验二 导弹跟踪问题.pptVIP

* * 实验二 导弹跟踪问题 一、实验目的 本实验主要涉及常微分方程。通过实验复习微分方程的建模和求解；介绍两种微分方程的数值方法：Euler法和改进的Euler法；还介绍了仿真方法。 二、实际问题 某军的一导弹基地发现正北方向120km处海面上有敌舰一艘以 90km/h 的速度向正东方向行驶。该基地立即发射导弹跟踪追击敌艇，导弹速度为 450km/h，自动导航系统使导弹在任一时刻都能对准敌艇。试问导弹在何时何处击中敌舰？ 三、数学模型 微分方程建模的方法主要是依据守恒律来建立等量关系式。对于这个问题，寻求等量关系是比较简单的。 设坐标系如图3.1所示，取导弹基地为原点 O(0,0)，x 轴指向正东方,y 轴指向正北方。 当 t=0 时，导弹位于点O，敌艇位于点 A(0,H)， 其中H=120（km）。设导弹在t时刻的位置为 P(x(t),y(t))，由题意, 其中 另外在 t 时刻, 敌艇位置为 , 其中 。由于导弹轨迹的切线方向必须指向敌艇，即直线 PM 的方向就是导弹轨迹上点 P 的切线方向, 故有 或写为 方程(3.1),(3.3)连同初值条件 构成了一个关于时间变量t的一阶微分方程组的初值问题。 为了寻求x与y的关系,要设法消去变量t, 由式(3.2)得 两边对t 求导 即有 把式(3.1)写为 代入上式,就得到轨迹方程.这是一个二阶非 线性微分方程,加上初值条件,则初值问题 就是导弹轨迹的数学模型。值得注意的是，前面导出的一阶微分方程组(3.1)，(3.3)和(3.4)实际上已经是一个数学模型了，不过多一个变量（或说参数）而已。 四、解析方法 方程（3.5）可以降阶. 令 则（3.5）化为一阶可分离变量 方程 即 易得 由初值条件(3.7)即 注意到上式可改写为 于是有 这样我们又得到一个可分离变量方程 积分得 利用 于是导弹轨迹方程为 设导弹击中敌艇于B(L,H), 以 Y=H 代入(3.9)式 , 得 而导弹击中敌艇的时刻 将数据 代入(3.10) , (3.11)式 , 得 五、数值方法 将初值问题(3.5) ~ (3.7) 化为一阶微分方程组 取自变量 y 的步长为 ,于是得分点 相应点上的 x 的值和 p 的值记为 显然 , 有初值 条件 我们将介绍两种近似算法来进行数值处理. Ⅰ Euler 方法 Euler 方法十分简单 , 就是用差商代替微商, 即以 代之以 而 代之以 这样有 设导弹到达 处的时刻为 那么得到计算的迭代格式 于是 表 3.1 是取 n=4 时的计算结果 , 读者可以用来检验程序或应用软件的正确性. 表 3.1 k yk xk pk 0 0 0 0 1 30 0 0.05 2 60