查看“︁开普勒定律”︁的源代码

{{refimprove|time=2015-07-26T16:00:58+00:00}}
{{noteTA
|T=zh-hans:开普勒定律; zh-hant:開普勒定律; zh-tw:克卜勒定律;
|G1=物理學
|1=zh-hans:开普勒; zh-hant:開普勒; zh-tw:克卜勒;
}}
[[File:Kepler laws diagram.svg|thumb|200px|图示遵守开普勒行星运动定律的两个行星轨道。（1）行星轨道是椭圆軌道。第一个行星的轨道焦点是f1与f2，第二个行星的轨道焦点是f1与f3。太阳的位置是在点f1。（2）A1与A2是两个面积相等的阴影区域。太陽与第一个行星的連線，扫过这两个阴影区域，所需的時間相等。（3）各个行星绕太阳公转周期的比率为a1<sup>3/2</sup>:a2<sup>3/2</sup> ；这里，a1与a2分别为第一个行星与第二个行星的半长轴长度]]

'''开普勒定律'''（{{lang-en|Kepler's law}}）是由[[德國]][[天文學家]]、[[數學家]][[約翰尼斯·克卜勒]]所发现的、关于[[行星]]运动的[[定律]]。开普勒於1609年在他出版的《'''[[新天文学]]'''》科學雜誌上发表了关于行星运动的两条定律，又於1618年，发现了第三条定律。

开普勒幸运地得到了著名丹麦天文学家[[第谷·布拉赫]]所观察与收集、且非常精确的天文資料。大约于1605年，根据布拉赫的行星位置資料，开普勒发现行星的移动遵守著三条相当简单的定律。同年年底，他撰寫完成了發表文稿。但是，直到1609年，才在《新天文学》科學雜誌發表，這是因為布拉赫的觀察數據屬於他的繼承人，不能隨便讓別人使用，因此產生的一些法律糾紛造成了延遲。

在天文学与物理学上、开普勒的定律给予[[亚里士多德]]派与[[托勒密]]派极大的挑战。他主张[[地球]]是不斷地移动的；[[行星轨道]]不是[[圓形]]（{{lang|en|epicycle}}）的，而是[[椭圆]]形的；行星[[公转]]的速度不等恒。这些论点，大大地动摇了当时的天文学与物理学。经过几乎一個世纪的研究，物理学家终于能够運用物理理论解释原理。[[艾萨克·牛顿|牛顿]]應用他的[[牛顿第二定律|第二定律]]和[[万有引力定律]]，在数学上严格地証明了开普勒定律，也让人了解了其中的[[物理]]意义。

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

=== 开普勒第一定律 ===
[[File:kepler-first-law.svg|thumb|根據开普勒第一定律，太陽位於橢圓軌道的一個焦點]] 
'''开普勒的第一定律'''，亦称為'''椭圆定律、軌道定律'''：每一个[[行星]]都沿著各自的[[椭圆]]轨道环绕著[[太阳]]，而太阳则处在椭圆的一个[[焦点]]上。<ref name="克卜勒的行星運動三定律">{{Cite web|url=http://teacher2.kyu.edu.tw/nstr/phx/phx06.pdf|title=克卜勒的行星運動三定律|accessdate=2015-07-26|publisher=[[高苑科技大學]]|archive-url=https://web.archive.org/web/20160305020238/http://teacher2.kyu.edu.tw/nstr/phx/phx06.pdf|archive-date=2016-03-05|dead-url=yes}}</ref>{{clear}}

=== 开普勒第二定律 ===
[[File:Kepler-second-law.svg|right|thumb|200px|根據开普勒第二定律，在同樣時間間隔內，行星繞著太陽公轉所掃過的面積相等]]
'''开普勒第二定律'''，也称為'''等面积定律'''：在相等[[时间]]内，[[太阳]]和[[運動 (物理學)|运动]]着的[[行星]]的连线所扫过的[[面积]]都是相等的。<ref name="克卜勒的行星運動三定律"/>

这一定律实际揭示了行星绕太阳公转的[[角动量]]守恒。用公式表示为
:<math>S_{AB}=S_{CD}=S_{EK}</math> 。
由此定律可以推知，行星离太阳较近时，运行速度较大，离太阳较远时，运行速度较小。{{clear}}

=== 开普勒第三定律 ===
[[File:Kepler third law.svg|right|thumb|200px|根據开普勒第三定律，行星繞著太陽公轉的周期平方和它们的橢圆轨道的[[半长轴]]立方成[[正比]]]]
'''开普勒第三定律'''，也稱為'''週期定律'''：各个[[行星]]绕[[太阳]][[公转]][[周期]]的[[平方]]及其[[椭圆]][[軌道 (力學)|轨道]]的[[半长轴]]的立方成[[正比]]。<ref name="克卜勒的行星運動三定律"/>
用公式表示为
:<math>\frac{a^3}{\tau^2}=k</math> ；
其中<math>k</math>是常量，它只与中心天体有关。

由这一定律导出：行星与太阳之间的引力与半径的平方成反比。这是[[艾萨克·牛顿]]的[[万有引力定律]]的一个重要基础。

用公式表示为
:<math>\frac{\tau^2}{a^3}=K</math> ；

这里，<math>a</math> 是行星公转轨道半长轴，<math>\tau</math> 是行星公转周期，<math>K</math> 是常数。

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

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

=== 开普勒第二定律推导 ===
[[牛顿万有引力定律]]表明，任意两个粒子由通过连线方向的力相互吸引。该引力的大小与它们的[[质量]]乘积成正比，与它们距离的平方成反比。由于太阳超重于行星，可以假设太阳是固定的。用方程式表示，
:<math> \boldsymbol{F}= - G \frac{mM}{r^2}\ \hat{\boldsymbol{r}}</math> ；

这裏，<math>\boldsymbol{F}</math> 是太阳作用於行星的万有引力，<math>m</math> 是行星的质量，<math>M</math> 是太阳的质量，<math>\boldsymbol{r}</math> 是行星相对于太阳的[[位移]]向量，<math>\hat{\boldsymbol{r}}</math> 是 <math>\boldsymbol{r}</math> 的[[单位向量]]。

[[牛顿第二定律]]表明，物體受力後所产生的[[加速度]] <math>\ddot{\boldsymbol{r}}</math> ，和其所受的淨力 <math>\boldsymbol{F}</math> 成正比，和其質量 <math>m</math> 成反比，以方程式表示，
:<math>\boldsymbol{F}=m\ddot{\boldsymbol{r}}</math> 。

合并这两个方程式，
:<math>\ddot{\boldsymbol{r}}= - G \frac{M}{r^2}\ \hat{\boldsymbol{r}}</math> 。 <span style="position:absolute;right:15%">(1)</span>

思考位置向量 <math>\boldsymbol{r}= r \hat{\boldsymbol{r}}</math> ，對於时间 <math>t</math> 微分一次可得到速度向量，再微分一次则可得到加速度向量：
:<math>\dot{\boldsymbol{r}}=\dot{r}\hat{\boldsymbol{r}} +r\dot{\theta}\hat{\boldsymbol{\theta}}</math> ，
:<math>\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}}</math> 。<span style="position:absolute;right:15%">(2)</span>

在这里，用到了单位向量微分方程式：
:<math>\frac{\mathrm{d}\hat{\boldsymbol{r}}}{\mathrm{d}t} = \dot{\theta}  \hat{\boldsymbol{\theta}}</math> ，
:<math>\frac{\mathrm{d}\hat{\boldsymbol{\theta}}}{\mathrm{d}t} = - \dot{\theta} \hat{\boldsymbol{r}}</math> 。

其中，<math>\theta </math> 是向量 <math>\hat{\boldsymbol{r}}</math> 的極角，<math>\hat{\boldsymbol{\theta}}</math> 是該向量逆時針旋轉 90 度的向量。

合并方程式 (1) 与 (2) ，可以得到向量运动方程式：
:<math> (\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}}</math>

取各个分量，可以得到两个[[常微分方程|常微分方程式]]，一个是关于径向加速度，另一个是关于切向加速度：
:<math>\ddot{r} - r\dot{\theta}^2 = - \frac{GM}{r^2}</math> ，<span style="position:absolute;right:15%">(3)</span>
:<math>r\ddot{\theta} + 2\dot{r}\dot{\theta} = 0</math> 。<span style="position:absolute;right:15%">(4)</span>

导引开普勒第二定律只需切向加速度方程式。

试想行星的[[角动量]] <math>\ell=m r^2\dot{\theta}</math> 。由于行星的质量是常数，角动量對於时间的导数为
:<math>\dot{\ell}=mr (2\dot{r}\dot{\theta}+r\ddot{\theta})=0</math> 。

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

在很小的時間 <math>dt</math> 中，恆星與行星連線掃過的面積 <math>dA</math> 近似於一個以 <math>r</math> 為底，
以 <math>r \tan d\theta \approx r d\theta</math> 為高的三角形，故

从时间 <math>t_1</math> 到时间 <math>t_2</math> 扫过的区域 <math>\Delta A</math> ，
:<math>\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) </math> 。

行星太阳连线扫过的区域面积相依于间隔时间 <math>t_2 - t_1</math> 。所以，开普勒第二定律是正确的。

=== 开普勒第一定律推导 ===
设定 <math>u =\frac{1}{r}</math> 。这样，角速度是
:<math>\dot{\theta} =\frac{\ell}{mr^2}=\frac{\ell u^2}{m}</math> 。

对时间微分和对角度微分有如下关系：
:<math>\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}</math> 。

根据上述关系，径向距离 <math>r=\frac{1}{u}</math> 对时间的导数为：
:<math>\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} </math> 。

再求一次导数：
:<math>\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}</math> 。

代入径向运动方程式 (3) ，<math>\ddot{r} - r\dot{\theta}^2 = - \frac{GM}{r^2}</math> ，
:<math> - \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</math> 。

将此方程式除以 <math> - \frac{\ell^2 u^2}{m^2}</math> ，则可得到一个简单的[[微分方程#常係数非齐次线性全微分方程|常係数非齐次线性全微分方程式]]来描述行星轨道：
:<math>\frac{\mathrm{d}^2u}{\mathrm{d}\theta^2} + u = \frac{GMm^2}{\ell^2}</math> 。

为了解这个微分方程，先列出一个特解
:<math>u=\frac{GMm^2}{\ell^2}</math> 。

再求解剩余的[[微分方程#常系数齐次线性全微分方程|常系数齐次线性全微分方程式]]，
:<math>\frac{\mathrm{d}^2u}{\mathrm{d}\theta^2} + u = 0</math> 。

它的解为
:<math>u=C\ \cos(\theta - \theta_0) </math> ；

这里，<math>C</math> 与 <math>\theta_0</math> 是常数。合并特解和与齐次方程解，可以得到通解
:<math>u=\frac{GMm^2}{\ell^2}+C\ \cos(\theta-\theta_0)</math> 。

选择坐标轴，让 <math>\theta_0=0</math> 。代回 <math>u=\frac{1}{r}</math> ，
:<math>\frac{1}{r}=\frac{GMm^2}{\ell^2}(1+ e\ \cos{\theta})</math> ；

其中，<math>e=C\ell^2/GMm^2</math> 是[[离心率]]。

这是[[圆锥曲线]]的[[圆锥曲线#极坐标|极坐标方程式]]，坐标系的原点是圆锥曲线的焦点之一。假若 <math>0<e< 1</math> ，则 <math>r</math> 所描述的是椭圆轨道。這证明了开普勒第一定律。

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

[[File:Elipse.svg|thumb|right|200px|行星环绕太阳（焦点 F1 ）的椭圆轨道。]]
开普勒第一定律阐明，行星环绕太阳的轨道是椭圆形的。[[椭圆#面积和周长|椭圆]]的面积是 <math>\pi ab</math> ；这里，<math>a</math> 与 <math>b</math> 分别为椭圆的[[椭圆|半长轴]]与[[椭圆|半短轴]]。在开普勒第二定律推导里，行星－太阳连线扫过区域速度 <math>\frac{\mathrm{d}A}{\mathrm{d}t}</math> 为
:<math>\frac{\mathrm{d}A}{\mathrm{d}t}=\frac{\ell}{2m}</math> 。

所以，行星公转周期 <math>\tau</math> 为
:<math>\tau=\frac{2m\pi ab}{\ell}</math> 。<span style="position:absolute;right:60%">(5)</span>

关于此行星环绕太阳，椭圆的半長軸 <math>a</math> ，半短軸 <math>b</math> 与'''近拱距''' <math>r_A</math> （[[拱点|近拱点]] A 与引力中心之间的距离），'''远拱距''' <math>r_B</math> （[[拱点|远拱点]] B 与引力中心之间的距离）的关系分别为
:<math>a=\frac{r_A+r_B}{2}</math> ，<span style="position:absolute;right:15%">(6)</span>
:<math>b=\sqrt{r_A\ r_B}</math> 。<span style="position:absolute;right:15%">(7)</span>

如果想要知道半長軸与半短軸，必须先求得近拱距与远拱距。依据[[能量守恒定律]]，
:<math>E=\frac{1}{2}m\dot{r}^2+\frac{1}{2}mr^2\dot{\theta}^2 - G\frac{mM}{r}</math> 。

在近拱点 A 与远拱点 B，径向速度都等于零：
:<math>\dot{r}=0</math> 。

所以，
:<math>E=\frac{1}{2}mr^2\dot{\theta}^2 - G\frac{mM}{r}= \frac{\ell^2}{2mr^2} - G\frac{mM}{r}</math> 。

稍微加以编排，可以得到 <math>r</math> 的一元二次方程式：
:<math>r^2+\frac{GmM}{E}r -  \frac{\ell^2}{2mE}=0</math> 。

其兩個根分别为椭圆轨道的近拱距 <math>r_A</math> 与远拱距 <math>r_B</math> 。
:<math>r_A=-\frac{\frac{GmM}{E} - \sqrt{\left(\frac{GmM}{E}\right)^2+\frac{2\ell^2}{mE}}}{2}</math> ；
:<math>r_B=-\frac{\frac{GmM}{E} + \sqrt{\left(\frac{GmM}{E}\right)^2+\frac{2\ell^2}{mE}}}{2}</math> 。

代入方程式 (6) 与 (7) ，
:<math>a= - \frac{GmM}{2E}</math> ，
:<math>b=\frac{\ell}{\sqrt{ - 2mE}}=\frac{\ell}{m}\frac{\sqrt{a}}{\sqrt{GM}}</math> 。

代入方程式 (5) ，周期的方程式为
:<math>\tau=\frac{2\pi a^{3/2}}{\sqrt{GM}}</math> 。

[[Q.E.D.|{{Unicode|■}}]]

== 參見 ==
* [[约翰内斯·开普勒]]
* [[牛頓第二運動定律]]
* [[万有引力定律]]
* 《[[新天文学]]》
* 《[[世界的和諧]]》
* [[牛頓旋轉軌道定理]]

==參考資料==
<references />

{{約翰尼斯·克卜勒}}
{{Authority control}}
[[Category:经典力学|K]]
[[Category:天体力学|K]]
[[Category:物理定律|K]]
[[Category:约翰内斯·开普勒]]
[[Category:哥白尼革命]]
[[Category:天文學方程式]]