%龍格庫塔右端函數，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=asin（2*（y（9）*y（10）+y（8）*y（11）））-y（2）;
ABS_alpha=abs（alpha）;
%插值部分開始
P=interp1（APTAPy（7））*g*10;
CX=twointerp（AaAMaACXABS_alphaMa）;

CY=twointerp（AaAMaACYABS_alphaMa）;

if?alpha<0
????CY=-CY;
end
XF=CX*q*s;
YF=CY*q*s;
wzt=y（3）*L/y（1）;

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

if?y（7）??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（7））;
dy（1）=（P*cos（alpha）-XF-9.81*y（6）*sin（y（2）））/y（6）;
dy（2）=（P*sin（alpha）+YF-9.81*y（6）*cos（y（2）））/（y（6）*y（1））;
dy（3）=（MZAF*alpha+MZWZ*wzt）*q*s*L/JZ;
dy（4）=y（1）*cos（y（2））;
dy（5）=y（1）*sin（y（2））;
dy（6）=-MC;
dy（7）=1;
dy（8）=-0.5*y（11）*y（3）;
dy（9）=0.5*y（10）*y（3）;
dy（10）=-0.5*y（9）*y（3）;
dy（11）=0.5*y（8）*y（3）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-05-12?10:30??垂直發射彈道仿真程序（采用四元數法）\
?????文件?????????385??2004-11-16?09:05??垂直發射彈道仿真程序（采用四元數法）\ACX.txt
?????文件?????????421??2007-04-18?14:39??垂直發射彈道仿真程序（采用四元數法）\ACY.txt
?????文件??????????82??2004-11-16?09:13??垂直發射彈道仿真程序（采用四元數法）\AJZ.txt
?????文件??????????75??2004-11-16?09:15??垂直發射彈道仿真程序（采用四元數法）\AJZT.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
?????文件??????????43??2004-11-16?09:22??垂直發射彈道仿真程序（采用四元數法）\AMa.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
?????文件??????????18??2004-11-17?05:17??垂直發射彈道仿真程序（采用四元數法）\Aa.txt
?????文件????????1686??2007-04-19?21:43??垂直發射彈道仿真程序（采用四元數法）\grkt2f.asv
?????文件????????1689??2007-04-19?21:48??垂直發射彈道仿真程序（采用四元數法）\grkt2f.m
?????文件?????????413??2004-06-24?10:58??垂直發射彈道仿真程序（采用四元數法）\oneinterp1.m
?????文件?????????440??2004-06-13?11:18??垂直發射彈道仿真程序（采用四元數法）\oneinterp2.m
?????文件???????13084??2008-11-05?16:18??垂直發射彈道仿真程序（采用四元數法）\result.dat
?????文件?????????353??2004-11-16?07:38??垂直發射彈道仿真程序（采用四元數法）\rk.m
?????文件????????2034??2007-04-19?21:48??垂直發射彈道仿真程序（采用四元數法）\trajectory.asv
?????文件????????2106??2008-11-05?16:18??垂直發射彈道仿真程序（采用四元數法）\trajectory.m
?????文件?????????562??2004-06-24?10:05??垂直發射彈道仿真程序（采用四元數法）\twointerp.m