小作业1--导弹

4. 导弹追踪问题：设位于坐标原点的甲舰向位于x轴上点A(1, 0)处的乙舰发射导弹，导弹头始终对准乙舰.如果乙舰以最大的速度(是常数)沿平行于y轴的直线行驶，导弹的速度是5，模拟导弹运行的轨迹.又乙舰行驶多远时，导弹将它击中？

轨迹

clear

x=0:0.01:1;

y=-5*(1-x).^(4/5)/8+5*(1-x).^(6/5)/12+5/24;

plot(x,y,'*')

解题

1．建立M文件eq1．m

function dy=eq1(x,y)

dy=zeros(2,1);

dy(1)=y(2);

dy(2)=1/5*sqrt(1+y(1)^2)/(1-x);

2．取x0=0，xf=0．9999，建立主程序ff6．m如下：

x0=0，xf=0．9999

[x,y]=ode15s('eq1',[x0 xf],[0 0]);

plot(x,y(:,1),’b．')

hold on

y=0:0．01:2;

plot(1,y,’b）