%龍格庫塔右端函數，dy:導數值；y:積分值

function?dy=grkt2f（y）

global?AJZ;
global?AJZT;
global?AMa;
global?Aa;
global?ACX;
global?ACY;
global?AP;
global?APT;
global?AXG;
global?AXGT;
global?AMZWZ1;
global?AMZWZ2;
global?AMZAF;

RHO=1.225;
SONIC=343.13;

s=0.0227;?????????????%特征面積
q=0.5*RHO*y（1）*y（1）;??????????%動壓頭
Ma=y（1）/SONIC;??????????%計算馬赫數
L=1.8;??????????????????%特征長度
B=[2.1126?44.0832];??????%兩級發動機結束工作時間??????
g=9.81;


alpha=y（4）-y（2）;
ABS_alpha=abs（alpha）;
%插值部分開始
P=interp1（APTAPy（8））*g;
CX=twointerp（AaAMaACXABS_alphaMa）;
CY=twointerp（AaAMaACXABS_alphaMa）;
if?alpha<0
????CY=-CY;
end
XF=CX*q*s;
YF=CY*q*s;
wzt=y（3）*L/y（1）;

if?y（8）<=B（1）?
????MC=2.362;
else??if?y（8）>B（2）?
?????????MC=0;
?????????P=0;
??????else?
?????????MC=0.21059;
??????end
??end
??
if?y（8）XG=interp1（AXGTAXGy（8））;??%插值重心
end

if?y（8）??aa=twointerp（AaAMaAMZWZ1ABS_alphaMa）;
??bb=twointerp（AaAMaAMZWZ2ABS_alphaMa）;
??ab=[aabb];
??a=[.9831.8896];
??MZWZ=interp1（aabXG）;
else
????MZWZ=twointerp（AaAMaAMZWZ2ABS_alphaMa）;
end

MZAF=twointerp（AaAMaAMZAFABS_alphaMa）;?
if?ABS_alpha>0
????MZAF=MZAF*57.3+CY*（XG-0.9381）/（L*alpha）;?%由于重心的改變，對MZAF進行修正
end
%插值部分結束

JZ=interp1（AJZTAJZy（8））;
dy（1）=（P*cos（alpha）-XF-9.81*y（7）*sin（y（2）））/y（7）;
dy（2）=（P*sin（alpha）+YF-9.81*y（7）*cos（y（2）））/（y（7）*y（1））;
dy（3）=（MZAF*alpha+MZWZ*wzt）*q*s*L/JZ;
dy（4）=y（3）;
dy（5）=y（1）*cos（y（2））;
dy（6）=y（1）*sin（y（2））;
dy（7）=-MC;
dy（8）=1;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????????18??2004-11-17?05:17??對導彈進行彈道分析\Aa.txt

?????文件????????385??2004-11-16?09:05??對導彈進行彈道分析\ACX.txt

?????文件????????376??2004-11-16?09:10??對導彈進行彈道分析\ACY.txt

?????文件?????????82??2004-11-16?09:13??對導彈進行彈道分析\AJZ.txt

?????文件?????????75??2004-11-16?09:15??對導彈進行彈道分析\AJZT.txt

?????文件?????????43??2004-11-16?09:22??對導彈進行彈道分析\AMa.txt

?????文件????????385??2004-11-16?09:27??對導彈進行彈道分析\AMZAF.txt

?????文件????????378??2004-11-16?09:32??對導彈進行彈道分析\AMZWZ1.txt

?????文件????????376??2004-11-16?09:39??對導彈進行彈道分析\AMZWZ2.txt

?????文件?????????74??2004-11-18?11:07??對導彈進行彈道分析\AP.txt

?????文件?????????64??2004-11-18?11:27??對導彈進行彈道分析\APT.txt

?????文件?????????68??2004-11-16?09:18??對導彈進行彈道分析\AXG.txt

?????文件?????????54??2004-11-16?09:20??對導彈進行彈道分析\AXGT.txt

?????文件???????1466??2004-11-20?07:22??對導彈進行彈道分析\grkt2f.asv

?????文件???????1570??2004-11-29?11:07??對導彈進行彈道分析\grkt2f.m

?????文件????????413??2004-06-24?10:58??對導彈進行彈道分析\oneinterp1.m

?????文件????????440??2004-06-13?11:18??對導彈進行彈道分析\oneinterp2.m

?????文件??????35806??2009-05-07?13:34??對導彈進行彈道分析\result.dat

?????文件????????353??2004-11-16?07:38??對導彈進行彈道分析\rk.m

?????文件???????1462??2004-11-20?07:17??對導彈進行彈道分析\trajectory.asv

?????文件???????1462??2004-11-20?11:51??對導彈進行彈道分析\trajectory.m

?????文件????????562??2004-06-24?10:05??對導彈進行彈道分析\twointerp.m

?????文件???????1909??2009-05-07?13:33??對導彈進行彈道分析\untitled.fig

?????文件???????7824??2009-05-07?13:35??對導彈進行彈道分析\untitled1.fig

?????文件???????2118??2009-05-07?13:35??對導彈進行彈道分析\untitled2.fig

?????文件???????1838??2009-05-07?13:36??對導彈進行彈道分析\untitled3.fig

?????目錄??????????0??2009-05-07?13:36??對導彈進行彈道分析

?????文件???????3543??2009-02-01?16:19??使用幫助：新手必看.htm

?????文件????????183??2009-02-02?09:47??Matlab中文論壇--助努力的人完成畢業設計.url

-----------?---------??----------?-----??----

............此處省略2個文件信息