克卜勒與行星運動定律

一個被釘在火星上的謎題

1600 年，年輕的德國數學家克卜勒(Johannes Kepler)抵達布拉格，成為丹麥天文學家第谷(Tycho Brahe)的助手。第谷交給他一個看似簡單的任務：算出火星的軌道。克卜勒當時誇口，八天就能搞定。結果，他足足耗費了將近八年。

火星之所以難纏，是因為它的軌道偏離正圓的程度，在當時可觀測的行星中相當顯著。第谷留下的觀測資料精度達到約 $2$ 角分(arcminute)——這是肉眼觀測的極限，望遠鏡尚未發明。正是這份近乎苛刻的精度，讓克卜勒無法用任何「圓形軌道加上小修正」的方案蒙混過關。當他發現自己最好的圓形模型仍與觀測差了 $8$ 角分時，他沒有把這個微小誤差掃到地毯底下，而是說：「就憑這 $8$ 角分，我們將重建整個天文學。」這句話，最終真的成立了。

第谷的遺產：一座沒有理論的資料金山

要理解克卜勒的成就，得先理解他站在誰的肩膀上。

第谷是望遠鏡發明前最偉大的觀測者。他在丹麥的烏拉尼堡(Uraniborg)天文台用巨大的象限儀(quadrant)與六分儀，數十年如一日地記錄行星位置，累積了人類史上第一批達到角分級精度的連續天文資料。然而第谷本人並不相信日心說，他抱持一種折衷的「第谷體系」，認為太陽繞地球、其他行星繞太陽。他擁有金山般的資料，卻沒有能解開它的理論框架。

克卜勒則相反：他是堅定的哥白尼(Copernicus)日心說支持者，也是一位帶有神祕主義色彩的數學家，深信宇宙的結構必定遵循簡潔的數學和諧。當這兩種特質——第谷的資料精度與克卜勒的數學執著——在布拉格相遇，近代天文學的轉捩點就此誕生。這是一個「資料驅動科學」的早期典範：理論必須臣服於觀測，而非觀測屈就於理論。

第一定律：行星走的是橢圓，不是圓

千百年來，從亞里斯多德(Aristotle)到托勒密(Ptolemy)，西方天文學有一條近乎信仰的鐵則：天體運動必須是「完美」的均速圓周運動。為了讓圓形模型符合觀測，前人疊加了本輪(epicycle)、均輪(deferent)、偏心(eccentric)等層層補丁，模型愈來愈臃腫。

克卜勒做出了驚人的一躍：他放棄了圓。

克卜勒第一定律（橢圓軌道定律）：每顆行星都沿著橢圓(ellipse)軌道運行，太陽位於橢圓的其中一個焦點(focus)上。

橢圓是一條封閉曲線，其上任意一點到兩個焦點的距離之和為定值。它的形狀由半長軸(semi-major axis) $a$ 與離心率(eccentricity) $e$ 描述。離心率衡量橢圓「壓扁」的程度：$e=0$ 時退化為正圓，$e$ 愈接近 $1$ 則愈狹長。

在極座標下，以焦點（太陽）為原點，行星到太陽的距離 $r$ 隨真近點角(true anomaly)$\theta$ 變化：

$$r = \frac{a(1 - e^2)}{1 + e\cos\theta}$$

行星離太陽最近的點稱為近日點(perihelion)，距離為 $a(1-e)$；最遠的點稱為遠日點(aphelion)，距離為 $a(1+e)$。

值得澄清一個常見迷思：地球軌道的離心率僅約 $e \approx 0.0167$，非常接近正圓。如果把地球軌道畫在紙上，肉眼幾乎分辨不出它與正圓的差別。因此四季的成因並非地球與太陽的遠近——事實上，北半球的冬天（每年一月初）正是地球通過近日點、離太陽最近的時候。四季真正的成因是地球自轉軸約 $23.5^\circ$ 的傾角。火星的離心率則高達 $e \approx 0.0934$，這正是它當年讓克卜勒焦頭爛額的原因。

第二定律：行星在近日點跑得快

確立了軌道形狀後，克卜勒接著問：行星沿橢圓運行的速度是均勻的嗎？答案是否定的。

克卜勒第二定律（面積速度定律）：連接行星與太陽的假想線段，在相等的時間內掃過相等的面積。

這意味著行星並非等速運動。當它靠近太陽（近日點附近）時，因為線段較短，必須跑得更快才能掃出相同面積；當它遠離太陽（遠日點附近）時則跑得較慢。地球在一月的近日點公轉速度約 $30.3\ \mathrm{km/s}$，七月的遠日點則減慢到約 $29.3\ \mathrm{km/s}$。

第二定律的物理本質，是角動量守恆(conservation of angular momentum)。由於太陽對行星的重力始終指向太陽（即所謂的中心力，central force），它對太陽不產生力矩，因此行星相對太陽的角動量保持不變。面積掃過的速率正比於角動量：

$$\frac{dA}{dt} = \frac{1}{2} r^2 \frac{d\theta}{dt} = \frac{L}{2m} = \text{常數}$$

其中 $L$ 是角動量、$m$ 是行星質量。克卜勒當然不知道「角動量」這個概念——它要再等大半個世紀才被牛頓(Newton)的力學系統化。克卜勒只是從第谷的數字中,赤手空拳地把這個規律「看」了出來。

第三定律：軌道愈大，一年愈長

前兩條定律各自描述單一行星，而第三定律是克卜勒花了十年才在 1619 年發現的「行星之間」的關係，它把整個太陽系串成了一個有機的整體。

克卜勒第三定律（週期定律）：行星公轉週期(orbital period)$T$ 的平方，正比於其軌道半長軸 $a$ 的立方。

$$T^2 \propto a^3$$

若以年(year)為週期單位、以天文單位(astronomical unit, AU；地球到太陽的平均距離)為長度單位，這個比例式對太陽系行星而言可以漂亮地簡化為：

$$\left(\frac{T}{\text{年}}\right)^2 = \left(\frac{a}{\text{AU}}\right)^3$$

動手算一下：火星的一年有多長？

火星的軌道半長軸約為 $a = 1.524\ \mathrm{AU}$。用第三定律來預測它的公轉週期：

$$T = a^{3/2} = (1.524)^{3/2}\ \text{年}$$

先算 $1.524^3 \approx 3.54$，再開平方根：

$$T = \sqrt{3.54} \approx 1.88\ \text{年}$$

也就是說，火星上的「一年」約等於地球的 $1.88$ 年，約 $687$ 個地球日。這與天文觀測值完全吻合。

我們再用木星檢驗一次。木星 $a \approx 5.20\ \mathrm{AU}$：

$$T = (5.20)^{3/2} = \sqrt{5.20^3} = \sqrt{140.6} \approx 11.86\ \text{年}$$

木星繞太陽一圈約需 $11.86$ 年，同樣與實測相符。第三定律的威力在於：只要量出一顆天體到太陽的距離，就能立刻推算它的「年」有多長，反之亦然。這成了後世測繪整個太陽系尺度的基石。

軌道要素：如何用六個數字鎖定一顆行星

克卜勒定律告訴我們軌道的形狀與運行規律，但要在三維空間中完整描述一條軌道、並指出行星此刻的位置，天文學家使用一組稱為軌道要素(orbital elements) 的參數。標準的六個要素為：

• 半長軸 $a$：決定軌道大小（連帶決定週期）。
• 離心率 $e$：決定軌道形狀（圓、橢圓、拋物線或雙曲線）。
• 軌道傾角(inclination) $i$：軌道平面相對於參考平面（太陽系常用黃道面，ecliptic）的傾斜角度。
• 升交點黃經(longitude of the ascending node) $\Omega$：軌道由南往北穿越參考平面之處的方位。
• 近日點幅角(argument of perihelion) $\omega$：近日點在軌道平面內的方位。
• 真近點角(true anomaly) $\theta$ 或平近點角(mean anomaly) $M$：指定行星在某時刻位於軌道上的哪個位置。

前兩個 $(a, e)$ 描述軌道的「大小與形狀」，中間三個 $(i, \Omega, \omega)$ 描述軌道在空間中的「方位」，最後一個則是「時間定位」。這六個數字，就是現代航太工程發射探測器、天文學家追蹤小行星與彗星時所依賴的座標語言。

重點回顧

• 第一定律（橢圓）：行星沿橢圓軌道運行，太陽位於其中一個焦點，而非軌道中心。
• 第二定律（面積速度）：行星與太陽的連線在等時間內掃過等面積，意味著近日點快、遠日點慢，其本質是角動量守恆。
• 第三定律（週期）：$T^2 \propto a^3$，把太陽系所有行星的軌道大小與週期統一在同一條定律下。
• 資料勝過信仰：克卜勒因不肯放過 $8$ 角分的誤差而推翻了兩千年的圓形軌道教條，是觀測精度驅動理論革命的經典案例。
• 軌道要素：六個參數 $(a, e, i, \Omega, \omega, \theta)$ 在三維空間中完整鎖定一條軌道，是現代航太與天文追蹤的通用語言。

深入探討（研究所視角）

克卜勒的三大定律是純粹從觀測歸納出的經驗律(empirical laws)——它們描述了行星「怎麼動」，卻沒回答「為什麼這樣動」。這個更深的問題，要等到牛頓在 1687 年的《自然哲學的數學原理》(Principia)中，以萬有引力與運動定律一舉解答。牛頓證明：克卜勒三大定律全都是平方反比重力(inverse-square gravity) 的數學後果。下面我們勾勒這條推導路徑。

二體問題與化約質量

考慮太陽（質量 $M$）與行星（質量 $m$）在彼此重力下運動，這就是二體問題(two-body problem)。兩者都繞著共同質心(barycenter)運動，但我們可以透過引入化約質量(reduced mass) $\mu = \dfrac{mM}{m+M}$，把問題轉化為「一個質量為 $\mu$ 的粒子，在固定中心的重力場中運動」的等效單體問題。相對位置向量 $\mathbf{r}$ 滿足：

$$\mu \ddot{\mathbf{r}} = -\frac{GMm}{r^2}\hat{\mathbf{r}}$$

其中 $G$ 是萬有引力常數。由於重力是中心力且沿連線方向，對質心的力矩為零，因此角動量 $\mathbf{L} = \mu\, \mathbf{r}\times\dot{\mathbf{r}}$ 守恆。角動量守恆有兩個直接推論：其一，運動被限制在一個固定平面內（軌道是平面曲線）；其二，面積速度 $\dfrac{dA}{dt}=\dfrac{L}{2\mu}$ 為常數——這正是克卜勒第二定律。

軌道方程與第一定律

利用能量守恆與角動量守恆，可將徑向運動寫成有效位能(effective potential)的問題：

$$E = \frac{1}{2}\mu \dot{r}^2 + \frac{L^2}{2\mu r^2} - \frac{GMm}{r}$$

式中 $\dfrac{L^2}{2\mu r^2}$ 是離心位能(centrifugal potential)，$-\dfrac{GMm}{r}$ 是重力位能。透過變數代換 $u = 1/r$ 並對軌道角度求解，可得到所謂的比涅方程(Binet equation)，其解為圓錐曲線(conic section)：

$$r = \frac{p}{1 + e\cos\theta}, \qquad p = \frac{L^2}{\mu\, GMm}$$

這個解的離心率為

$$e = \sqrt{1 + \frac{2EL^2}{\mu (GMm)^2}}$$

軌道的性質由總能量 $E$ 決定：$E<0$ 時 $e<1$，軌道是封閉的橢圓（束縛軌道，bound orbit）——這就是克卜勒第一定律；$E=0$ 時 $e=1$ 為拋物線（逃逸臨界）；$E>0$ 時 $e>1$ 為雙曲線（如某些掠日彗星或星際天體）。克卜勒只觀測到行星，所以只看到橢圓這一個特例；牛頓的理論則把它推廣成了一整族圓錐曲線。

推導第三定律 $T^2 \propto a^3$

最後，我們從重力導出第三定律。對橢圓軌道積分面積速度，一個完整週期 $T$ 掃過的總面積等於橢圓面積 $\pi a b$（$b$ 為半短軸）：

$$\frac{L}{2\mu}\, T = \pi a b$$

利用橢圓的幾何關係 $b = a\sqrt{1-e^2}$，以及半正焦弦 $p = a(1-e^2) = \dfrac{L^2}{\mu GMm}$，可解出 $L$ 並代回上式。經過整理（這是一段值得讀者自行動手完成的代數），得到：

$$T^2 = \frac{4\pi^2}{G(M+m)}\, a^3$$

這就是牛頓修正後的克卜勒第三定律。它揭示了一件克卜勒原版定律所沒有的深刻訊息：比例「常數」其實正比於 $\dfrac{1}{G(M+m)}$，它並非對所有系統都相同，而是取決於中心天體（加上繞行天體）的總質量。

這正是天文學最強大的「秤」：只要觀測一個天體繞另一個天體的軌道半長軸 $a$ 與週期 $T$，就能反推出系統的總質量 $M+m$。天文學家正是用這把秤，量出了太陽的質量、各行星的質量（透過它們的衛星）、雙星系統的質量，乃至銀河系中心那個質量達數百萬倍太陽的超大質量黑洞。在多數情況下行星質量 $m \ll M$，可近似 $T^2 \approx \dfrac{4\pi^2}{GM}a^3$，於是不同行星共享同一比例常數，回到克卜勒當年看到的簡潔形式。

從第谷一筆一筆抄錄的角分級數字，到克卜勒歸納的三條經驗律，再到牛頓以單一重力定律將它們統一——這條跨越近一個世紀的路徑，是科學如何從「描述現象」走向「解釋機制」的最佳縮影。今天我們發射的每一具行星探測器、追蹤的每一顆系外行星，都仍走在克卜勒與牛頓鋪下的這條軌道上。