daodandafeiji.m

来自「此段程序描述了一个导弹打击目标飞机的过程的导弹的飞行路线问题。导弹迅速飞行」· M 代码 · 共 19 行

clear
ts=0:0.1:12;
y0=[0,0]; 
[t,y]=ode45('fun2',ts,y0); 
%fun2为定义的函数，表示微分方程组： 
%           dx/dt=800(-500+500t-x)/sqrt((-500+500t-x)^2+(5000-y)^2); 
%           dy/dt=800(5000-y)/sqrt((-500+500t-x)^2+(5000-y)^2)];
grid on 
plot(y(:,1),y(:,2))   %导弹轨迹曲线 
title('导弹轨迹曲线') 
grid on 
n=1;
while n<121 
    if(5000-y(n,2)<5)
        fprintf('t=%f(s)时导弹击中目标，此时导弹高度y=%f(m)\n',ts(n),y(n,2))
        break;
    else n=n+1;
    end
end