文档介绍：防空导弹主级燃料相对质量因数计算
[文档副标题]

目录
一要求1
要求一1
要求二1
二已知条件1
1、分离条件1
2、气动参数1
3、发动机参数1
4、导引规律1
5、大气模型1
三相对量模型1
四算法1
五计算结果与分析1
1、迎角随水平距离变化曲线1
2、迎角随相对质量因数的变化率1
3、导弹飞行高度随相对质量因数的变化1
4、导弹和目标飞行高度随水平位置的变化1
5、各个参数随时间的变化率1
六代码1
一要求
要求一
根据已知条件，采用数值积分法求解相对量运动微分方程组，计算燃料相对质量因数。
要求二
综合运用积分、插值计算等计算方法，采用C、C++，或者Matlab等语言的一种，编制计算程序。
二已知条件
1、分离条件
速度V=500m/s，时间t=3s，x方向位置x=674m，y方向位置y=329m，初始攻角，初始弹道倾角。
2、气动参数
阻力因数插值表：
（）
2
4
6
8
10

升力因数斜率插值表：
（）
1
2
4
6
8
10

注：以上攻角单位均为度。
3、发动机参数
比冲Is=2156；重力加速度g=；推重比；翼载。
4、导引规律
三点法。目标匀速直线等高迎头飞行，，，
（弹目斜距）。
5、大气模型
参见教科书。
三相对量模型
三点法导引的相对量运动微分方程组为：
四算法
求解三点法导引的相对量运动微分方程组时，在这里采用了四级四阶经典
Runge-Kutta公式，计算公式为：
五计算结果与分析
经计算，，。各参数的变化趋势图如下所示：
1、迎角随水平距离变化曲线
从图中可以看出，迎角在不断的增加。当击中目标时，攻角为：在可接受的围，认为没有失速
2、迎角随相对质量因数的变化率
从图中可以看出，随着质量因数的变化，应交在不断增加。
3、导弹飞行高度随相对质量因数的变化
随着相对质量因数的变化，飞行高度在不断增加。
4、导弹和目标飞行高度随水平位置的变化
5、各个参数随时间的变化率
六代码
Runge-Kutta代码
clear;
clc;
close all;
format long
disp('计算导弹相对质量因数：');
% 分离条件
v0=500;t0=3;x0=674;y0=329;
alpha0=;theta0=26;%各角度均为角度制
% 动力系数
is=2156;g=;P=;p0=5880;
% 阻力因数差值表
alpha1=2:2:10;
ma1=[ ];
cx0=[