function dy=Rocket2(t,y)
dy=[y(2);-9.8-0.4*y(2).^2/320];
end
ts=0;60;
x0=[0,0];
[t,x]=ode45(@Rocket2,ts,x0);
[t,x];
for n=1:2000
T=80-0.01*n;
tss=60:0.02:T;
y0=[x(61,1),x(61,2)];
option=odeset('reltol',1e-3,'abstol',1e-6);
[t2,y]=ode45(@Rocket2,tss,y0,option);
[t2,y];
if y(:,2)>=0
break
end
end
plot(t,x(:,1),'b't2,y(:,1),'r'),grid,
title('图1高度-时间')
xlabel('t/s')
ylabel('h/m')
pause
plot(t,x(:,2),'b't2,y(:,2),'r'),grid,
title('图2速度-时间')
xlabel('t/s')
ylabel('v/(m/s)')
pause
a=(32000-0.4*x(:,2).^2)./(1400-18*t)-9.8;
a2=-9.8-0.4*y(;,2).^2/320;
plot(t,a,'b',t2,a2'r'),grid,
title('图3加速度-时间')
xlabel('t/s')
ylabel('a/(m/s^2)')
你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答
本次提问扣除的有问必答次数，将会以问答VIP体验卡（1次有问必答机会、商城购买实体图书享受95折优惠）的形式为您补发到账户。
因为有问必答VIP体验卡有效期仅有1天，您在需要使用的时候【私信】联系我，我会为您补发。