matlab ode45仿真,ode45仿真地空导弹，找不出错误！SOS

导弹弹道仿真，已纠结多日，尚无进展。求帮助！！

目标以300m/s的速度进入我方防空阵地，飞行高度为15km，作水平匀速直线运动。

我地空弹对其实施拦截。发射时目标视线角为30°，推力与初速在程序中有给出。

导引方式为比例导引法。

function ydot=bili(t,y)

% y(1)=theta  y(2)=V  y(3)=x  y(4)=y y(5)=xm  y(6)=ym  y(7)=m

T=288.15-0.0065*y(4);%温度与高度的关系，用于算当地音速

gamma=1.4;

r=287;

a=sqrt(T*r*gamma);%当地音速

Re=6370000;%地球半径

g=9.8*(Re/(Re+y(4)))^2;%当地重力加速度

rho=1.225*exp(-0.00015*y(4));%当地空气密度

s=1.94;       %机翼面积

M=y(2)/a;%马赫数

Vm=-300;%目标速度

K=3;%比例系数

e=180/pi;

q=atan((y(6)-y(4))/(y(5)-y(3)));

eta=q*e-y(1);

R=(y(6)-y(4))/sin(q);

%导弹推力，质量与时间的关系

t1=2.79;t2=t1+40.82;

dotm=219.45*(t<=t1)+7.42*(t>