16世纪之前,人们认为地球处于中心,而太阳和其他行星围绕着它转。后来,哥白尼提出了“日心说”,把太阳放在中心,地球和其他行星围绕着太阳转。虽然哥白尼的日心说是一项壮举,但该理论中行星是圆轨道,为了与观测数据相符,需要添加一个又一个的本轮。
16世纪末期,开普勒分析了他的老师第谷留下来的观测数据,整理出了著名的开普勒三大定律,其中第一定律又称为椭圆定律:所有行星轨道都是椭圆,太阳处于其焦点。我们想,开普勒是从实践中得出椭圆轨道,那么如何才能从理论上推导出椭圆轨道。
椭圆的数学表达
为了做到这一点,我们需要知道椭圆的数学表达式。在高中,我们都学过直角坐标系下的椭圆标准方程:
如上图,参数a、b、c所代表的含义都已标明,我们主要还有:偏心率e=c/a,半正焦弦l=a(1-e^2)。
事实上,在分析行星运动时,坐标原点取太阳中心比取椭圆中心更为自然。因此,我们将用极坐标系(r, θ)代替直角坐标系(x, y)。我们选取焦点F2为原点,于是有下面的关系:
把它代入上述直角坐标系方程中,得到:
求解这个方程,并剔除掉使r成负数的解,我们就可以得到椭圆在极坐标系下的方程:
为了等下推导行星轨道更容易看出椭圆方程,我们用半正焦弦l代替,并重新排列方程:
行星轨道方程
在这里,我们要写出极坐标下的行星能量方程和角动量方程。
请注意,在极坐标下,速度包含了两部分:平行r方向的速度和垂直r方向的速度。因此,我们可以列出能量方程=动能+势能:
其中,k=GMm。而角动量方程为:
联立两个方程,我们可以得到:
因为行星轨道方程与时间无关,所以我们可以利用链式法则消去dt,如下:
代入上面的式子就得到:
注意,这里还要用到一个常用的技巧:令u=1/r,这会使数学计算更容易一点。这里我们就不给出繁琐的计算过程,直接给出最终的轨道方程:
是不是和极坐标下的椭圆方程形式相同,其中偏心率e: