开普勒定律

图示遵守开普勒行星运动定律的两个行星轨道。(1) 行星轨道是椭圆轨道。第一个行星的轨道焦点是f1与f2,第二个行星的轨道焦点是f1与f3。太阳的位置是在点f1。(2) A1与A2是两个面积相等的阴影区域。太阳与第一个行星的连线,扫过这两个阴影区域,所需的时间相等。(3) 各个行星绕太阳公转周期的比率为a1^{3/2}:a2^{3/2};这里,a1与a2分别为第一个行星与第二个行星的半长轴长度。

开普勒定律是由德国天文、数学家约翰尼斯·开普勒所发现、关于行星运动的定律。他于1609年在他出版的《新天文学》科学杂志上发表了关于行星运动的两条定律,又于1618年,发现了第三条定律。

开普勒幸运地得到了著名丹麦天文学家第谷·布拉赫所观察与收集、且非常精确的天文资料。大约于1605年,根据布拉赫的行星位置资料,开普勒发现行星的移动遵守着三条相当简单的定律。同年年底,他撰写完成了发表文稿。但是,直到1609年,才在《新天文学》科学杂志发表,这是因为布拉赫的观察数据属于他的继承人,不能随便让别人使用,因此产生的一些法律纠纷造成了延迟。

在天文学与物理学上、开普勒的定律给予亚里士多德派与托勒密派极大的挑战。他主张地球是不断地移动的;行星轨道不是圆形(epicycle)的,而是椭圆形的;行星公转的速度不等恒。这些论点,大大地动摇了当时的天文学与物理学。经过几乎一个世纪的研究,物理学家终于能够运用物理理论解释原理。牛顿应用他的第二定律和万有引力定律,在数学上严格地证明了开普勒定律,也让人了解了其中的物理意义。

一、开普勒定律

开普勒的三条行星运动定律改变了整个天文学,彻底摧毁了托勒密复杂的宇宙体系,完善并简化了哥白尼的日心说。

开普勒第一定律: 根据开普勒第一定律,太阳位于椭圆轨道的一个焦点。开普勒的第一定律,也称为椭圆定律、轨道定律:每一个行星都沿各自的椭圆轨道环绕太阳,而太阳则处在椭圆的一个焦点中。

开普勒第二定律:根据开普勒第二定律,在同样时间间隔内,行星绕着太阳公转所扫过的面积相等。开普勒第二定律,也称为等面积定律:在相等时间内,太阳和运动着的行星的连线所扫过的面积都是相等的。

这一定律实际揭示了行星绕太阳公转的角动量守恒。用公式表示为

S_{{AB}}=S_{{CD}}=S_{{EK}}

开普勒第三定律: 根据开普勒第三定律,行星绕着太阳公转的周期平方和它们的椭圆轨道的半长轴立方成正比。开普勒第三定律,也称为周期定律:各个行星绕太阳公转周期的平方和它们的椭圆轨道的半长轴的立方成正比。

由这一定律不难导出:行星与太阳之间的引力与半径的平方成反比。这是艾萨克·牛顿的万有引力定律的一个重要基础。用公式表示为

{\frac {\tau ^{2}}{a^{3}}}=K

这里,a是行星公转轨道半长轴,\tau是行星公转周期,K是常数。

二、数学推导

开普勒定律描述的是行星围绕太阳的运动,牛顿定律可以更广义地描述几个粒子因万有引力相互吸引而形成的运动。假设只有两个粒子,其中一个粒子超轻于另外一个粒子,则轻的粒子会绕着重的粒子运动,就好似行星根据开普勒定律绕著太阳运动。另外,牛顿定律还可计算出关于其它方面的解答,行星轨道可以呈抛物线运动或双曲线运动。这是开普勒定律所无法预测到的结果。在一个粒子并不超轻于另外一个粒子的状况下,依照广义二体问题的解答,每一个粒子会绕着它们的共同质心运动。这也是开普勒定律无法预测到的。

开普勒定律使用几何语言将行星的坐标及时间跟轨道参数相连结。牛顿第二定律是一个微分方程。开普勒定律的推导涉及一些解析微分方程的技巧。在推导开普勒第一定律之前,必须先推导出开普勒第二定律,因为开普勒第一定律需要用到开普勒第二定律里的一些计算结果。

牛顿万有引力定律表明,任意两个粒子由通过连线方向的力相互吸引。该引力的的大小与它们的质量乘积成正比,与它们距离的平方成反比。由于太阳超重于行星,可以假设太阳是固定的。用方程表示

{\boldsymbol {F}}=-G{\frac {mM}{r^{2}}}\ {\hat {{\boldsymbol {r}}}}

这里,{\boldsymbol {F}}是太阳作用于行星的万有引力,m是行星的质量,M是太阳的质量,{\boldsymbol {r}}是行星相对于太阳的位移矢量,\hat{\boldsymbol{r}}{\boldsymbol {r}}的单位矢量。

牛顿第二定律表明,物体受力后所产生的加速度{\ddot {{\boldsymbol {r}}}},和其所受的合力{\boldsymbol {F}} 成正比,和其质量m成反比,以方程表示:

{\boldsymbol {F}}=m{\ddot {{\boldsymbol {r}}}}

合并这两个方程,

{\ddot {{\boldsymbol {r}}}}=-G{\frac {M}{r^{2}}}\ {\hat {{\boldsymbol {r}}}} \tag{1}

思考位置矢量{\boldsymbol {r}}=r{\hat {{\boldsymbol {r}}}},对于时间t微分一次可得到速度矢量,再微分一次则可得到加速度矢量:

\begin{aligned} & {\dot {{\boldsymbol {r}}}}={\dot {r}}{\hat {{\boldsymbol {r}}}}+r{\dot {\theta }}{\hat {{\boldsymbol {\theta }}}} \\ & {\ddot {{\boldsymbol {r}}}}=\left({\ddot {r}}{\hat {{\boldsymbol {r}}}}+{\dot {r}}{\frac {{\mathrm {d}}{\hat {{\boldsymbol {r}}}}}{{\mathrm {d}}t}}\right)+\left({\dot {r}}{\dot {\theta }}{\hat {{\boldsymbol {\theta }}}}+r{\ddot {\theta }}{\hat {{\boldsymbol {\theta }}}}+r{\dot {\theta }}{\frac {{\mathrm {d}}{\hat {{\boldsymbol {\theta }}}}}{{\mathrm {d}}t}}\right)=({\ddot {r}}-r{\dot {\theta }}^{2}){\hat {{\boldsymbol {r}}}}+(r{\ddot {\theta }}+2{\dot {r}}{\dot {\theta }}){\hat {{\boldsymbol {\theta }}}} \end{aligned} \tag{2}

在这里,用到了单位矢量微分方程:

\begin{aligned} & {\frac {{\mathrm {d}}{\hat {{\boldsymbol {r}}}}}{{\mathrm {d}}t}}={\dot {\theta }}{\hat {{\boldsymbol {\theta }}}} , \\ & {\frac {{\mathrm {d}}{\hat {{\boldsymbol {\theta }}}}}{{\mathrm {d}}t}}=-{\dot {\theta }}{\hat {{\boldsymbol {r}}}} 。 \end{aligned}

合并方程(1)与(2),可以得到矢量运动方程:

({\ddot {r}}-r{\dot {\theta }}^{2}){\hat {{\boldsymbol {r}}}}+(r{\ddot {\theta }}+2{\dot {r}}{\dot {\theta }}){\hat {{\boldsymbol {\theta }}}}=-{\frac {GM}{r^{2}}}{\hat {{\boldsymbol {r}}}}

取各个分量,可以得到两个常微分方程,一个是关于径向加速度,另一个是关于切向加速度:

\begin{aligned} & {\ddot {r}}-r{\dot {\theta }}^{2}=-{\frac {GM}{r^{2}}} \\ & r{\ddot {\theta }}+2{\dot {r}}{\dot {\theta }}=0 \end{aligned}

导引开普勒第二定律只需切向加速度方程。试想行星的角动量\ell =mr^{2}{\dot {\theta }}。由于行星的质量是常数,角动量对于时间的导数为

{\dot {\ell }}=mr(2{\dot {r}}{\dot {\theta }}+r{\ddot {\theta }})=0

角动量\ell也是一个运动常数,即使距离r与角速度{\dot {\theta }}都可能会随时间变化。

从时间t_{1}到时间t_{2}扫过的区域\Delta A,

\Delta A=\int _{{t_{1}}}^{{t_{2}}}{\frac {1}{2}}\cdot r\cdot r{\dot \theta }\cdot {\mathrm {d}}t=\int _{{t_{1}}}^{{t_{2}}}{\frac {\ell }{2m}}{\mathrm {d}}t={\frac {\ell }{2m}}\cdot (t_{2}-t_{1}) 。

行星太阳连线扫过的区域面积相依于间隔时间t_{2}-t_{1}。所以,开普勒第二定律是正确的。

设定u={\frac {1}{r}}。这样,角速度是

{\dot {\theta }}={\frac {\ell }{mr^{2}}}={\frac {\ell u^{2}}{m}} 。

对时间微分和对角度微分有如下关系:

{\frac {{\mathrm {d}}}{{\mathrm {d}}t}}={\dot {\theta }}{\frac {{\mathrm {d}}}{{\mathrm {d}}\theta }}={\frac {\ell u^{2}}{m}}{\frac {{\mathrm {d}}}{{\mathrm {d}}\theta }} 。

根据上述关系,径向距离r={\frac {1}{u}}对时间的导数为:

{\dot {r}}={\frac {\ell u^{2}}{m}}{\frac {{\mathrm {d}}}{{\mathrm {d}}\theta }}{\frac {1}{u}}=-{\frac {\ell u^{2}}{m}}{\frac {1}{u^{2}}}{\frac {{\mathrm {d}}u}{{\mathrm {d}}\theta }}=-{\frac {\ell }{m}}{\frac {{\mathrm {d}}u}{{\mathrm {d}}\theta }} 。

再求一次导数:

{\ddot {r}}={\frac {\ell u^{2}}{m}}{\frac {{\mathrm {d}}{\dot {r}}}{{\mathrm {d}}\theta }}={\frac {\ell u^{2}}{m}}{\frac {{\mathrm {d}}}{{\mathrm {d}}\theta }}(-{\frac {\ell }{m}}{\frac {{\mathrm {d}}u}{{\mathrm {d}}\theta }})=-{\frac {\ell ^{2}u^{2}}{m^{2}}}{\frac {{\mathrm {d}}^{2}u}{{\mathrm {d}}\theta ^{2}}} 。

代入径向运动方程(3),{\ddot {r}}-r{\dot {\theta }}^{2}=-{\frac {GM}{r^{2}}},

-{\frac {\ell ^{2}u^{2}}{m^{2}}}{\frac {{\mathrm {d}}^{2}u}{{\mathrm {d}}\theta ^{2}}}-{\frac {\ell ^{2}u^{3}}{m^{2}}}=-GMu^{2} 。

将此方程除以-{\frac {\ell ^{2}u^{2}}{m^{2}}},则可得到一个简单的常系数非齐次线性全微分方程来描述行星轨道:

{\frac {{\mathrm {d}}^{2}u}{{\mathrm {d}}\theta ^{2}}}+u={\frac {GMm^{2}}{\ell ^{2}}} 。

为了解这个微分方程,先列出一个特解

u={\frac {GMm^{2}}{\ell ^{2}}} 。

再求解剩余的常系数齐次线性全微分方程,

{\frac {{\mathrm {d}}^{2}u}{{\mathrm {d}}\theta ^{2}}}+u=0 。

它的解为

u=C\ \cos(\theta -\theta _{0}) ;

这里,C\theta _{0}是常数。合并特解和与齐次方程解,可以得到通解

u={\frac {GMm^{2}}{\ell ^{2}}}+C\ \cos(\theta -\theta _{0}) 。

选择坐标轴,让\theta _{0}=0。代回u={\frac {1}{r}},

{\frac {1}{r}}={\frac {GMm^{2}}{\ell ^{2}}}(1+e\ \cos {\theta }) ;

其中,e=C\ell ^{2}/GMm^{2}是离心率。

这是圆锥曲线的极坐标方程,坐标系的原点是圆锥曲线的焦点之一。假若0<e<1,则r所描述的是椭圆轨道。这证明了开普勒第一定律。

在建立牛顿万有引力定律的概念与数学架构上,开普勒第三定律是牛顿依据的重要线索之一。假若接受牛顿运动定律。试想一个虚拟行星环绕着太阳公转,行星的移动轨道恰巧呈圆形,轨道半径为r。那么,太阳作用于行星的万有引力为F={\frac {mv^{2}}{r}}。行星移动速度为v={\frac {2\pi r}{\tau }}。依照开普勒第三定律,这速度v与半径的平方根{\sqrt {r}}成反比。所以,万有引力F\propto {\frac {1}{r^{2}}}。猜想这大概是牛顿发现万有引力定律的思路,但这个猜想无法被证实,因为在他的计算本里,并没有找到任何关于这方面的证据。

开普勒第一定律阐明,行星环绕太阳的轨道是椭圆形的。椭圆的面积是\pi ab;这里,ab分别为椭圆的半长轴与半短轴。在开普勒第二定律推导里,行星-太阳连线扫过区域速度{\frac {{\mathrm {d}}A}{{\mathrm {d}}t}}

{\frac {{\mathrm {d}}A}{{\mathrm {d}}t}}={\frac {\ell }{2m}} 。

所以,行星公转周期\tau

\tau ={\frac {2m\pi ab}{\ell }} \tag{5}

关于此行星环绕太阳,椭圆的半长轴a,半短轴b与近拱距r_{A}(近拱点A与引力中心之间的距离),远拱距r_{B}(远拱点 B 与引力中心之间的距离)的关系分别为

a={\frac {r_{A}+r_{B}}{2}} \tag{6}
b={\sqrt {r_{A}\ r_{B}}} \tag{7}

如果想要知道半长轴与半短轴,必须先求得近拱距与远拱距。依据能量守恒定律,

E={\frac {1}{2}}m{\dot {r}}^{2}+{\frac {1}{2}}mr^{2}{\dot {\theta }}^{2}-G{\frac {mM}{r}}

在近拱点A与远拱点B,径向速度都等于零:

{\dot {r}}=0

所以,

E={\frac {1}{2}}mr^{2}{\dot {\theta }}^{2}-G{\frac {mM}{r}}={\frac {\ell ^{2}}{2mr^{2}}}-G{\frac {mM}{r}}

稍为加以编排,可以得到r的一元二次方程:

r^{2}+{\frac {GmM}{E}}r-{\frac {\ell ^{2}}{2mE}}=0

其两个根分别为椭圆轨道的近拱距r_{A}与远拱距r_{B}

r_{A}=-{\frac {{\frac {GmM}{E}}-{\sqrt {\left({\frac {GmM}{E}}\right)^{2}+{\frac {2\ell ^{2}}{mE}}}}}{2}}
r_{B}=-{\frac {{\frac {GmM}{E}}+{\sqrt {\left({\frac {GmM}{E}}\right)^{2}+{\frac {2\ell ^{2}}{mE}}}}}{2}}

代入方程(6)与(7),

a=-{\frac {GmM}{2E}}
b={\frac {\ell }{{\sqrt {-2mE}}}}={\frac {\ell }{m}}{\frac {{\sqrt {a}}}{{\sqrt {GM}}}}

代入方程(5),周期的方程为

\tau ={\frac {2\pi a^{{3/2}}}{{\sqrt {GM}}}}