function?dydt=fun（ty）
dydt=zeros（41）;
d=0.229;
s=pi*d*d/4;
m=225;
Ma=[0.60.80.91.01.1];
Cxo=[0.149420.157540.171400.288320.34312];
Cxa2=[0.001860.001900.001920.002000.00209];
Cya=[0.069770.069480.069330.071140.07427];
ad=[-1.16-1.17-1.20-1.25-1.17];
Cydz=[0.033110.032440.031770.036690.03711];
Mzdz=[-0.1323-0.1323-0.1322-0.1531-0.1525];
Mzwz=[-13.0898-12.9252-13.8327-14.5677-15.3393];
rou=1.225*（1.0-2.0323*0.00001*y（4））^4.83;
T=288.34-5.86*10^（-3）.*y（4）;
sonic=20.046*T^0.5;
g=9.806*（1.0-2*y（4）/6371000）;
ma=y（1）./sonic;
Cx_o=interp1（MaCxoma‘spline‘）;
Cx_a2=interp1（MaCxa2ma‘spline‘）;
Cy_a=interp1（MaCyama‘spline‘）;
a_d=interp1（Maadma‘spline‘）;
Cy_dz=interp1（MaCydzma‘spline‘）;
if?t>3
????a=5;
else
????a=0;
end
Cx=Cx_o+Cx_a2*a*a;
Cy=Cy_a*a+Cy_dz*a/a_d;
X=0.5*Cx*rou*s*y（1）^2;
Y=0.5*Cy*rou*s*y（1）^2;
dydt（1）=（-X-m*g*sin（y（2）））/m;
dydt（2）=（Y-m*g*cos（y（2）））/（m*y（1））;
dydt（3）=y（1）*cos（y（2））;
dydt（4）=y（1）*sin（y（2））;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2016-01-05?19:21??王放?912000720220\
?????目錄???????????0??2016-01-05?19:21??王放?912000720220\導彈飛行理學作業--王放912000720220\
?????目錄???????????0??2016-01-05?19:21??王放?912000720220\導彈飛行理學作業--王放912000720220\ti1\
?????文件???????40448??2016-01-05?00:33??王放?912000720220\導彈飛行理學作業--王放912000720220\ti1\T1.doc
?????目錄???????????0??2016-01-05?19:21??王放?912000720220\導彈飛行理學作業--王放912000720220\ti2\
?????文件???????45590??2016-01-04?20:05??王放?912000720220\導彈飛行理學作業--王放912000720220\ti2\2.jpg
?????文件?????????860??2016-01-04?20:02??王放?912000720220\導彈飛行理學作業--王放912000720220\ti2\TI2.m
?????文件????????1000??2015-11-20?09:40??王放?912000720220\導彈飛行理學作業--王放912000720220\ti2\fun.m
?????文件?????1146505??2015-11-13?10:25??王放?912000720220\導彈飛行理學作業--王放912000720220\ti2\ti2答案.fig
?????目錄???????????0??2016-01-05?19:21??王放?912000720220\導彈飛行理學作業--王放912000720220\ti3\
?????文件???????62460??2016-01-04?20:32??王放?912000720220\導彈飛行理學作業--王放912000720220\ti3\3.jpg
?????文件????????1262??2015-11-20?10:35??王放?912000720220\導彈飛行理學作業--王放912000720220\ti3\fun.m
?????文件????????1028??2015-11-24?22:09??王放?912000720220\導彈飛行理學作業--王放912000720220\ti3\ti3.m
?????文件??????185673??2015-11-24?22:09??王放?912000720220\導彈飛行理學作業--王放912000720220\ti3\結果.fig
?????目錄???????????0??2016-01-05?19:21??王放?912000720220\導彈飛行理學作業--王放912000720220\ti4\
?????文件??????416666??2015-11-12?08:16??王放?912000720220\導彈飛行理學作業--王放912000720220\ti4\1.fig
?????文件??????357188??2015-11-12?08:16??王放?912000720220\導彈飛行理學作業--王放912000720220\ti4\2.fig
?????文件???????24031??2016-01-04?20:50??王放?912000720220\導彈飛行理學作業--王放912000720220\ti4\41.jpg
?????文件???????23610??2016-01-04?21:01??王放?912000720220\導彈飛行理學作業--王放912000720220\ti4\42.jpg
?????文件??????????28??2015-11-25?21:53??王放?912000720220\導彈飛行理學作業--王放912000720220\ti4\Untitled.m
?????文件?????????125??2015-11-25?22:24??王放?912000720220\導彈飛行理學作業--王放912000720220\ti4\fun1.m
?????文件?????????180??2015-11-26?16:52??王放?912000720220\導彈飛行理學作業--王放912000720220\ti4\fun2.m
?????文件?????????465??2016-01-04?20:50??王放?912000720220\導彈飛行理學作業--王放912000720220\ti4\t4.m
?????文件?????????599??2015-11-12?08:15??王放?912000720220\導彈飛行理學作業--王放912000720220\ti4\t42.m
?????目錄???????????0??2016-01-05?19:21??王放?912000720220\導彈飛行理學作業--王放912000720220\ti5\
?????文件?????????521??2015-11-26?21:38??王放?912000720220\導彈飛行理學作業--王放912000720220\ti5\t5.m
?????目錄???????????0??2016-01-05?19:21??王放?912000720220\導彈飛行理學作業--王放912000720220\ti6\
?????文件?????????131??2015-11-26?21:45??王放?912000720220\導彈飛行理學作業--王放912000720220\ti6\Ti6.m
?????目錄???????????0??2016-01-05?19:21??王放?912000720220\導彈飛行理學作業--王放912000720220\ti7\
?????文件?????????353??2015-11-26?22:01??王放?912000720220\導彈飛行理學作業--王放912000720220\ti7\ti7.m
?????目錄???????????0??2016-01-05?19:21??王放?912000720220\導彈飛行理學作業--王放912000720220\ti8\
............此處省略7個文件信息