综合程序设计案例一:导弹追踪问题

设位于坐标原点的甲舰向位于x轴上点A(10, 20)处的乙舰发射导弹,导弹头始终对准乙舰。如果乙舰以最大的速度v0(是常数)行驶,行驶轨迹满足曲线方程y=-4x^2+80x+20,导弹的速度是20v0,绘图表示导弹和乙舰行驶轨迹的曲线方程,并标注图形说明。

思路:将模块离散化 每一小段近似为直线,通过 斜率 按比例预计算

示意绘图

结果:

 

结论: 上图为结果 以乙船速度V 为 单位1 计算   经过 1.1713 后被击中 被击中时的坐标点为 (10.0147 ,21.1712)

实现代码:

 %思路: 将模型离散化,每一小段近似为直线
 function main
 x=10:0.0000001:10.05;%乙船 x方向运动
 y=-4*(x-10).^2+80*(x-10)+20;%乙船 y方向运动
 [x1,y1,sum_L,goal]=myf(x,y,20);%自定义函数预测导弹轨迹
 plot(x,y,'r.')%打印乙船曲线
 axis([0,15,0,24]);
 hold on;
 plot(x1,y1,'b-')%打印导弹曲线
 t=sum_L/1.0%击中时间
 x(goal)%击中x点
 y(goal)%击中y点
 plot(x(goal),y(goal),'g*')%打印标记击中点
   legend('乙船轨迹','导弹轨迹','击中点');
 end
 function [x1,y1,sum_L,goal]=myf(x,y,Vb)%已知乙船的运动轨迹,速度倍速Vb
 x1=[];%导弹 x方向运动
 y1=[];%导弹 y方向运动
 x1(1)=0;%导弹 x 初始坐标
 y1(1)=0;%导弹 y 初始坐标
 sum_L=0;%乙船总路程
 goal=0;%击中目标点索引
     for i=2:length(x)% 通过已知轨迹计算导弹轨迹
        L=sqrt((x(i)-x(i-1))^2+(y(i)-y(i-1))^2);%乙船运动距离通过公式 根号(x^2+y^2)
        sum_L=sum_L+L;%更新乙船总路程     通过斜率
        x1(i)=x1(i-1)+Vb*L*(x(i)-x1(i-1)) /sqrt((x(i)-x1(i-1))^2+(y(i)-y1(i-1))^2);%求导弹下一个位子x1值
        y1(i)=y1(i-1)+Vb*L*(y(i)-y1(i-1))/sqrt((x(i)-x1(i-1))^2+(y(i)-y1(i-1))^2);%求导弹下一个位子y1值
       if(x1(i)>=x(i)) goal=i; break;end
    end
end