克卜勒定律

維基百科,自由的百科全書
前往: 導覽搜尋
圖示遵守克卜勒行星運動定律的兩個行星軌道。 (1)行星軌道是橢圓軌道。第一個行星的軌道焦點是 f1f2 ,第二個行星的軌道焦點是 f1f3 。太陽的位置是在點 f1 。 (2)A1與A2是兩個面積相等的陰影區域。太陽與第一個行星的連線,掃過這兩個陰影區域,所需的時間相等。(3)各個行星繞太陽公轉周期的比率為a1^{3/2}:a2^{3/2} ;這裡,a1a2 分別為第一個行星與第二個行星的半長軸長度。

克卜勒定律克卜勒所發現、關於行星運動的定律。他於1609年在他出版的《新天文學》科學雜誌上發表了關於行星運動的兩條定律,又於1618年,發現了第三條定律。

克卜勒幸運地得到了著名丹麥天文學家第谷·布拉赫所觀察與收集、且非常精確的天文資料。大約於1605年,根據布拉赫的行星位置資料,克卜勒發現行星的移動遵守著三條相當簡單的定律。同年年底,他撰寫完成了發表文稿。但是,直到1609年,才在《新天文學》科學雜誌發表,這是因為布拉赫的觀察數據屬於他的繼承人,不能隨便讓別人使用,因此產生的一些法律糾紛造成了延遲。

在天文學與物理學上、克卜勒的定律給予亞里士多德派與托勒密派極大的挑戰。他主張地球是不斷地移動的;行星軌道不是周轉圓epicycle)的,而是橢圓形的;行星公轉的速度不等恆。這些論點,大大地動搖了當時的天文學與物理學。經過了幾乎一個世紀披星戴月,廢寢忘食的研究,物理學家終於能夠運用物理理論解釋其中的奧秘。牛頓應用他的第二定律萬有引力定律,在數學上嚴格地証明了克卜勒定律,也讓人們了解了其中的物理意義。

克卜勒定律[編輯]

克卜勒的三條行星運動定律改變了整個天文學,徹底摧毀了托勒密複雜的宇宙體系,完善並簡化了哥白尼日心說

克卜勒第一定律[編輯]

克卜勒第一定律

克卜勒第一定律,也稱橢圓定律、軌道定律:每一個行星都沿各自的橢圓軌道環繞太陽,而太陽則處在橢圓的一個焦點中。

克卜勒第二定律[編輯]

克卜勒第二定律

克卜勒第二定律,也稱等面積定律:在相等時間內,太陽運動著的行星的連線所掃過的面積都是相等的。

這一定律實際揭示了行星繞太陽公轉的角動量守恆。用公式表示為

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

克卜勒第三定律[編輯]

克卜勒第三定律

克卜勒第三定律,也稱週期定律:各個行星太陽公轉周期平方和它們的橢圓軌道半長軸立方正比

由這一定律不難導出:行星與太陽之間的重力與半徑的平方成反比。這是牛頓萬有引力定律的一個重要基礎。

用公式表示為

\frac{a^3}{\tau^2}=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}}(1)

思考位置向量 \boldsymbol{r}= r \hat{\boldsymbol{r}} ,對於時間 t 微分一次可得到速度向量,再微分一次則可得到加速度向量:

\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}}(2)

在這裡,用到了單位向量微分方程式:

\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}}

合併方程式 (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}}

取各個分量,可以得到兩個常微分方程式,一個是關於徑向加速度,另一個是關於切向加速度:

\ddot{r} - r\dot{\theta}^2 = - \frac{GM}{r^2}(3)
r\ddot{\theta} + 2\dot{r}\dot{\theta} = 0(4)

導引克卜勒第二定律只需切向加速度方程式。試想行星的角動量 \ell=m r^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}^2u}{\mathrm{d}\theta^2} + u = \frac{GMm^2}{\ell^2}

為了解這個微分方程式,先列出一個特解

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

再求解剩餘的常係數齊次線性全微分方程式

\frac{\mathrm{d}^2u}{\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} 。猜想這大概是牛頓發現萬有引力定律的思路,但這個猜想無法被證實,因為在他的計算本里,並沒有找到任何關於這方面的證據。

行星環繞太陽(焦點 F1 )的橢圓軌道。

克卜勒第一定律闡明,行星環繞太陽的軌道是橢圓形的。橢圓的面積是 \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}(5)

關於此行星環繞太陽,橢圓的半長軸 a ,半短軸 b近拱距 r_A近拱點 A 與重力中心之間的距離),遠拱距 r_B遠拱點 B 與重力中心之間的距離)的關係分別為

a=(r_A+r_B) / 2(6)
b=\sqrt{r_A\ r_B}(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=\left( - \frac{GmM}{E} - \sqrt{\left(\frac{GmM}{E}\right)^2+\frac{2\ell^2}{mE}}\right) / 2
r_B=\left( - \frac{GmM}{E} + \sqrt{\left(\frac{GmM}{E}\right)^2+\frac{2\ell^2}{mE}}\right) / 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}}

Q.E.D.

參見[編輯]