%num是循环次数，y_j是击中时的y坐标（因为会有误差，所以击中时的导弹y坐标会与120km有偏差），L是击中时的x坐标，T是导弹击中时一共用的时间。x,y为导弹的位置，t是多长时间模拟一次，eps是最小误差
function [ num,y_j,L,T] = simu1( x,y,t,eps )
k=0;
while k<1000 %最多循环1000次
p=90*k*t-x; %kt时刻敌舰与导弹的横向距离
q=120-y; %kt时刻敌舰与导弹的纵向距离
d1=p/(p^2+q^2)^0.5;%计算cos(αk)
d2=q/(p^2+q^2)^0.5;%计算sin(αk)
%t时刻导弹的位置
x=x+450*t*d1;
y=y+450*t*d2;
if (abs(120-y)<eps)
break;
end
k=k+1;
end
num=k;
L=x;
y_j=y;
T=L/90;
end
%运行代码
%t=0.01; x=0;y=0;eps=1.5
%[num,y_j,L,T]=simu1(x,y,t,eps)
