资源简介
垂直发射弹道仿真程序(采用四元数法);matlab语言源代码,可以满足弹道仿真类毕业设计的计算程序
代码片段和文件信息
%龙格库塔右端函数,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
评论
共有 条评论