弹道导弹最大射程弹道设计

项目介绍

这是本人的课设。由于我的同学都用的祖传代码，而本人想自己写一份，做祖传代码的始祖。老实说并不完美，不过也是我自己努力的证明（嘻嘻~）

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