【气动学】三级火箭发射弹道主动段仿真（三次点火达到目标轨道）【含Matlab源码 4711期】

本文链接：https://blog.csdn.net/weixin_63266434/article/details/139744552

在这里插入图片描述

⛄一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【气动学】基于matlab三级火箭发射弹道主动段仿真（三次点火达到目标轨道）【含Matlab源码 4711期】
点击上面蓝色字体，直接付费下载，即可。

备注：
点击上面蓝色字体付费专栏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)