⛄一、获取代码方式
获取代码方式1：
完整代码已上传我的资源：【气动学】基于matlab三级火箭发射弹道主动段仿真（三次点火达到目标轨道）【含Matlab源码 4711期】
点击上面蓝色字体，直接付费下载，即可。
获取代码方式2：
付费专栏Matlab物理应用（初级版）
备注：
点击上面蓝色字体付费专栏Matlab物理应用（初级版），扫描上面二维码，付费29.9元订阅海神之光博客付费专栏Matlab物理应用（初级版），凭支付凭证，私信博主，可免费获得1份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：1份本博客上传CSDN资源代码
⛄二、部分源代码
clear all
%u=[v,sita,x,y,h];
R0=6371000;
miu=3.986e14;
u1=[1;pi/2;0;1;1];
tspan=0:0.1:187.5;
[t,u]=ode45(@fun,tspan,u1);
data=[t,u];
N=length(t);
v=zeros(N,1);
sita=zeros(N,1);
x=zeros(N,1);
y=zeros(N,1);
h=zeros(N,1);
m=zeros(N,1);
Pe=zeros(N,1);
fai=zeros(N,1);
rho=zeros(N,1);
alfa=zeros(N,1);
Ma=zeros(N,1);
vsonic=zeros(N,1);
ny=zeros(N,1);
shecheng0=zeros(N,1);
dasita=zeros(N,1);
beta=zeros(N,1);
sudu=zeros(N,1);
for i=1:N
t(i)=data(i,1);
v(i)=data(i,2);
sita(i)=data(i,3);
x(i)=data(i,4);
y(i)=data(i,5);
h(i)=data(i,6);
Pe(i)=Pulse(t(i));
m(i)=mass(t(i));
vsonic(i)=Vsonic(h(i));
Ma(i)=v(i)/Vsonic(h(i));
fai(i)=faipr(t(i));
alfa(i)=(fai(i)-sita(i))/pi*180;
beta(i)=atan(x(i)/(R0+y(i)));
dasita(i)=(sita(i)+beta(i))/pi*180;
sudu(i)=sqrt(miu/(R0+h(i)));
end
subplot(2,3,1),plot(t,fai),title(‘程序俯仰角’)
subplot(2,3,2),plot(t,sita),title(‘sita’)
subplot(2,3,3),plot(t,alfa),title(‘攻角’)
subplot(2,3,4),plot(t,v),title(‘速度’)
subplot(2,3,5),plot(t,h),title(‘高度’)
subplot(2,3,6),plot(t,dasita),title(‘当地速度倾角’)
plot(x,y)
plot(t,sita)
⛄三、运行结果
⛄四、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1] 门云阁.MATLAB物理计算与可视化[M].清华大学出版社，2013.
3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除