Question

微分方程组如何用matlab解角度为3度,g为重力加速度,取9.8,如何得到v和t的曲线

Answer 1

clc

f=inline('[y(2);-1.204*(y(2)^2+y(4)^2)*cos(pi/60);y(4);-1.204*(y(2)^2+y(4)^2)*sin(pi/60)-9.8]','t','y');

[t y]=ode45(f,[0 1],[初值]);

v=sqrt(y(:,2).^2+y(:,4).^2);

plot(t,v)

带入初值,运行即可,还有角度真的是3度?