• 大小: 15KB
    文件类型: .zip
    金币: 2
    下载: 1 次
    发布日期: 2021-05-12
  • 语言: Matlab
  • 标签: 弹道仿真  

资源简介

垂直发射弹道仿真程序(采用四元数法);matlab语言源代码,可以满足弹道仿真类毕业设计的计算程序

代码片段和文件信息

%龙格库塔右端函数,dy:导数值;y:积分值

function dy=grkt2fy

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(7XG=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

评论

共有 条评论