项目介绍
这是本人的课设。由于我的同学都用的祖传代码,而本人想自己写一份,做祖传代码的始祖。老实说并不完美,不过也是我自己努力的证明(嘻嘻~)
DF-X号是一款陆基单级液体燃料弹道导弹,本课程设计以该导弹为对象,在所给假设条件下,要求设计其最大射程飞行弹道,分析不同参数对导弹射程的影响规律,并给出主要参数的变化曲线。考虑①地球(视为标准椭球体)的引力②由导弹气动参数推算出的气动力③俯仰角控制程序④推力模型等,最终改变终端俯仰程序角 和起飞时的燃料装载质量mf的参数,以求得地球自转的最大射程解,输出最优轨迹。
求解过程
导弹控制方程
算了,太多了,不想打了
代码
主程序shechengjisuan.m
ae=6378145;
ee=1/298.257;
be=6356760.2713934;
we=7.292115e-5;
miu=3.986005e14;
J=1.5*1.08263e-3;
Cx=0.30;
Cy=0.05;
Sm=2.2;
long0=112.6;
long0=deg2rad(long0);
B0=37.5;
B0=deg2rad(B0);
A0=3*pi/2;
fai0=atan(tan(B0)*(1-ee)^2);
miu0=B0-fai0;
R0=ae*be/sqrt((ae*sin(fai0))^2 + (be*cos(fai0))^2);
R0x=-1*R0*sin(miu0)*cos(A0);
R0y=R0*cos(miu0);
R0z=R0*sin(miu0)*sin(A0);
R0shi=[R0x;R0y;R0z];
t=0;
vshi=[0.00001;0.00001;0.00001];
x=0.000001;
y=1500;
z=0.000001;
step=0.01;
h=1500;
m0=30000;
faiprf=24.7*pi/180;
m1=3000;
m2=1200;
P=4e5;
Is=320;
mdian=P/Is/9.8;
tk=(m0-m1)/mdian;
% tjilu=[];
% hjilu=[];
% shechengjilu=[];
% fiprjilu=[];
% alpjilu=[];
% sitajilu=[];
% rjilu=[];
% vjilu=[];
% mjilu=[];
% latijilu=[];
% longjilu=[];
while h>10 && t<=1800