追击曲线问题追击曲线问题.doc

追击曲线问题 一敌舰在某海域内沿正北方向航行时，我方战舰恰位于敌舰的正西方1n mile（海里）处。我舰向敌舰发射制导鱼雷，敌舰速度为0.42n mile/min，鱼雷的速度为敌舰速度的两倍。试问敌舰航行多远时将被击中？ 一、微分方程模型建立和求解 设敌舰的速度为常数v0，曲线为y=y(x)，即在时刻t，鱼雷的位置在点P(x,y)处，这时敌舰的位置在点Q(1,v0t)处（如图6.1所示）。由于鱼雷在追击过程中始终指向敌舰，而鱼雷的运动的方向是沿曲线的切线方向，所以有 两边对x求导，得 即 由已知鱼雷的速度为2v0 因为，所以 代入式(6.11)，得线y=y(x)满足的微分方程模型 令，方程化为 分离变量，积分并代入初始条件计算得 即 而 上两式相减 直接积分并代入初始条件得 这就是鱼雷追击曲线的方程。因鱼雷击中敌舰时，它的横坐标x=1，代入曲线方程得，即敌舰行至处时将被击中，这段航程所需要的时间为95.2381s。 三、用Mathematica求解的计算过程 为了利用Mathematica求解常微分方程功能求解方程式（6.12），输入命令 Dsolve[{y″[x]= =Sqrt[1+y′[x]^2]/(2*(1-x),y[0]= =0，y′[0]= =0)，y[x],x} 得计算结果如下 数学软件的局限性，上面用Mathematica求解常微分方程初值时出现的两个问解，第一个是错误的，只有第二个解正确，它与用解析解求出的答案式（6.16）一致，即 为了绘出解曲线图形并计算微分方程解在自变量x=1时的函数值，输入下面命令 y[x]:=2/3-(1+(x-1)/3)*Sqrt[1-x] Polt[y[x],{x,0,1}] Y[1] 得计算结果如下（包括如图6.3所示图形） Out[3]=-Graphics- Out[4]=2/3 图6.3 数据结果说明，当敌舰行至处时将被击中。  三、用计算机模拟鱼雷追击敌舰的过程 现在我们先介绍什么是计算机模拟，用计算机模拟实物系统，对实物系统的结构和行为进行动态演示，评价或预测系统的行为效果，为决策提供信息。这一实验技术称为计算机模拟（又称仿真）。在对真实系统做实验时，可能时间太长、费用太高、危险太大，甚至很难进行。采用计算机模拟技术，常常能获得满意的结果。根据模拟对象的不同特点，计算机模拟分为确定性模拟和随机性模拟两大类，求解本问题的模拟属于确定性模拟。问题的解虽然可以用解析的方法获得，但分析和计算的过程都较复杂。计算机模拟的方法简单可行，能在短时间内观察到鱼雷追击敌舰的全过程，并且可以估计鱼雷（或敌舰）的速度变化对过程的影响，我们采用时间步长法，按照时间流逝的顺序一步一步对敌舰和鱼雷的活动进行动态模拟。在整个模拟过程中，时间步长是固定不变的。 建立直角坐标系（如图6.1），设敌舰为动点Q，鱼雷为动点P。Q点的初始位置为Q0（1，y轴方向运动且速度为常数v0，故Q点在时刻t=tk的位置为Qk（1，v0tk）。 x 图6.1 由于P点运动方向始终指向，设在t=tk时刻，P的位置是Pk（xkyk），则向量=(1-xk，v0tk-yk)，此时P点的运动方向可由下面的单位向量（方向余弦）表示 P点运动速度为常数v1=2v0，取时间步长△t=1秒（s），设在t=tk+1时刻，P点的位置为，于P点位置变化规律为 现对追击过程进行模拟，当两个动点的距离小于0.02n mile时，则认为P点已经追上Q点。模拟过程实际上是产生平面上两个点列：的过程。 94 0.99966 0.65164 可写PASCAL程序或用EXCEL将动点P在各个时刻的坐标构成了两个二维元素的表。可用几何画板或者Mathematica得到追击曲线的图形（如图6.2）所示。运行的结果可以看出，以后鱼雷的位置为：而敌舰的纵坐标为0.658。由于此时鱼雷和敌舰的距离不超过0.02 n mile，可以认为当敌舰行至0.672n mile处时将被击中。  0 0.00000 0.00000 40 0.54794 0.09757 80 0.95518 0.46195 1 0.01400 0.00010 41 0.56085 0.10298 81 0.96067 0.47482 2 0.02800 0.00030 42 0.57369 0.10857 82 0.96583 0.48784 3 0.04200 0.00059 43 0.58645 0.11433 83 0.97066 0.50098 4 0.05599 0.00099 44 0.59913 0.12027 84 0.97513 0.51425