Question

用matlab求解

哈雷彗星在1986年2月9日到达了近日点（最接近太阳的点，取太阳为原点），那时它的位置和速度分别为

位置单位为AU(天文单位，取地球轨道的长半轴为单位距离)，时间单位为年。彗星的三维运动方程为

其中参数 ， 。求微分方程的数值解，作出彗星三维轨道和彗星轨道在yz平面的射影。由r与t的关系，计算彗星的远日点距太阳的距离，预测下一次彗星到达近日点的时间。

Answer 1

这是从哪弄来的题目？你确定 运动方程没错？我做出的结果怎么不对啊，在近日点是圆弧，后面好像是一条直线，离太阳越来越远，不再回来
dt=0.01;u=4*pi^2;
t(1)=0;
x(1)=0.325514;y(1)=-0.459460;z(1)=0.166229;
vx(1)=-9.096111;vy(1)=-6.916686;vz(1)=-1.305721;
x(2)=vx(1)*dt+x(1);y(2)=vy(1)*dt+y(1);z(2)=vz(1)*dt+z(1);
r(1)=sqrt(x(1)^2+y(1)^2+z(1)^2);
i=1;
while t<=76
t(i+1)=t(i)+dt;
r(i+1)=sqrt(x(i+1)^2+y(i+1)^2+z(i+1)^2);
vx(i+1)=-u*x(i)/r(i)^3*dt+vx(i);
vy(i+1)=-u*y(i)/r(i)^3*dt+vy(i);
vz(i+1)=-u*z(i)/r(i)^3*dt+vz(i);
x(i+2)=vx(i+1)*dt+x(i+1);
y(i+2)=vy(i+1)*dt+y(i+1);
z(i+2)=vz(i+1)*dt+z(i+1);
i=i+1;
end