Exercise_4.16

1.课堂回顾

第四章第一节主要是地球绕太阳做圆周运动，不考虑其他因素的影响。可画出圆轨道。
而第二节中考虑偏心率之后，在远日点速度最小，在近日点速度最大。可画出图像如下：

捕获.PNG

水星进动2.PNG

第三节，水星进动，考虑广义相对论修正后为

F_G\approx\frac{GM_SM_M}{r^2}(1+\frac{\alpha}{r^2})

利用欧拉迭代法可画出图像如下：

水星进动.PNG

第四节，我们考虑三体问题，在地球绕太阳转动时，加入木星对其的万有引力的影响。

F_{EJ,x}=-\frac{GM_JM_E(x_e-x_j)}{r_{EJ}^2}

F_{EJ,y}=-\frac{GM_JM_E(y_e-y_j)}{r_{EJ}^2}

可画出实际质量时，三体运动的轨迹：

实际质量.PNG

考虑太阳，地球，木星的三体问题，假设太阳不是静止不动的在原点，而是由三者的质心代替，给太阳一个初速度，使总动量为零，探究地球不同初值时的运动，和木星质量为实际质量的10,100,1000倍时的运动情况。

F_x=-G\frac{M_jM_s(x_j-x_s)}{r_sj^3}-G\frac{M_jM_e(x_j-x_e)}{r_sj^3}

F_y=-G\frac{M_jM_s(y_j-y_s)}{r_sj^3}-G\frac{M_jM_e(y_j-y_e)}{r_sj^3}

同理可得地球和太阳的受力。
确定质心：
设三个行星在一条直线上，太阳在左，地球和木星在右，距离为原始距离a(1+e)。利用质心公式，设太阳坐标为（-r_s），有：

M_s*r_s=M_e(a_e-r_s)+M_j(a_j-r_s)

可解得

r_s=\frac{M_ea_e+M_ja_j}{M_s+M_e+M_j}

因而可得三者的初始坐标，而，地球和木星的初始速度为v_min，再利用总动量为零可以求得太阳的初始速度：

v_{sy}=-\frac{M_j*v_{jy}+M_e*v_{ey}}{M_s}

3-body.PNG

对于改变地球的初值：

改变初始速度：
当地球速度变为2倍时，已经脱离了太阳的束缚，当木星与地球没有距离十分近的时候，木星仍然继续绕太阳运动。速度更大时，也将脱离太阳系。

v_x_e=2times.PNG
改变木星速度
当木星速度变为原来两倍时，木星将脱离太阳束缚而飞出太阳系。

V_x_j=2times.PNG
改变地球质量：
10倍

100times.PNG

1000倍

M_e=1000times.PNG
100000倍
M_e=100000times.PNG
1000000倍
M_e=1000000times.PNG

当地球的质量比太阳数量级小一时，使太阳的轨道产生较大影响，地球数量大于等于太阳的时，它们将远离质心。
改变木星质量
10倍

m_x_j=10times.PNG

100倍

m_x_j=100times.PNG
1000倍
m_x_j=1000times.PNG
从上面三幅图可以看出，当木星质量数量级远小于太阳的时，质心接近于在太阳处，当木星质量数量级比太阳小一时，质心将偏离太阳，三者都做圆周运动。当木星质量数量级大于等于太阳的时，三者将偏离质心，发散出去。