参考
http://zhidao.baidu.com/question/479644084.html?oldq=1 回答的一样的问题,还有下面的程序
function dy=odefun(t,y)
%y(1) x
%y(2) x'
%y(3) y
%y(4) y'
%y(5) z
%y(6) z'
u=9.8*1;
dy=zeros(6,1);
dy(1)=y(2);
dy(2)=-u*y(1)/sqrt(y(1)^2+y(3)^2+y(5)^2)^3;
dy(3)=y(4);
dy(4)=-u*y(3)/sqrt(y(1)^2+y(3)^2+y(5)^2)^3;
dy(5)=y(6);
dy(6)=-u*y(5)/sqrt(y(1)^2+y(3)^2+y(5)^2)^3;
end
命令行中
tspan=[0,50];
y0=[909;-156;173;445;173;445];
[t,y]=ode45(@odefun,[0,10],y0);
注意:u=GM需要自己修改
tspan为你的时间区间
y(end,:),就为你50s的解。