克卜勒定律與軌道力學：從二體問題到活力公式

從經驗定律到力學必然

克卜勒在 17 世紀初僅憑第谷的觀測資料，歸納出三條描述行星運動的經驗定律；但這些定律「為什麼成立」要等到牛頓的萬有引力與微積分才獲得解答。本篇不停留在「橢圓、面積、週期」的口訣，而是直接從二體問題的運動方程出發，推導出三大定律如何成為平方反比力場下的數學必然，並把它接到現代軌道力學的語言：軌道要素、活力公式（vis-viva）與軌道轉移。

我們的出發點是兩個質點在萬有引力下的運動。設太陽質量 $M$、行星質量 $m$，相對位置向量為 $\vec{r}$，則牛頓第二定律給出

$$\mu \ddot{\vec{r}} = -\frac{GMm}{r^2}\hat{r}, \qquad \mu = \frac{Mm}{M+m}$$

其中 $\mu$ 為約化質量。將兩質點運動化為單一相對座標的運動後，等效於一個質量 $\mu$ 的粒子在中心力場 $-\dfrac{k}{r^2}\hat{r}$（$k=GMm$）中運動。整個軌道力學的骨架，都長在這條方程上。

角動量守恆與第二定律

中心力沿著 $\hat{r}$，因此對力心的力矩為零：$\vec{\tau}=\vec{r}\times\vec{F}=0$，於是角動量

$$\vec{L} = \mu\,\vec{r}\times\dot{\vec{r}} = \text{常數}$$

守恆。這一步立即蘊含兩個結論。第一，$\vec{r}$ 與 $\dot{\vec{r}}$ 永遠垂直於固定的 $\vec{L}$，所以運動限制在一個平面內——軌道是平面曲線。第二，考慮掃過的面積微元 $dA=\tfrac12|\vec{r}\times d\vec{r}|$，於是面速度

$$\frac{dA}{dt} = \frac{1}{2}\left|\vec{r}\times\dot{\vec{r}}\right| = \frac{L}{2\mu} = \text{常數}$$

這正是克卜勒第二定律（等面積定律）。值得強調：第二定律只依賴「力為中心力」這個性質，與力是否為平方反比無關。換言之，等面積定律比橢圓軌道更普遍。

軌道方程與第一定律

要得到軌道形狀，採用極座標 $(r,\theta)$。在極座標下徑向運動方程為

$$\mu\left(\ddot{r}-r\dot{\theta}^2\right) = -\frac{k}{r^2}$$

利用 $L=\mu r^2\dot\theta$ 消去時間，並作經典變數代換 $u=1/r$，可將方程化為線性形式。由 $\dot r = -\dfrac{L}{\mu}\dfrac{du}{d\theta}$、$\ddot r = -\dfrac{L^2 u^2}{\mu^2}\dfrac{d^2u}{d\theta^2}$ 代入，整理得到著名的 Binet 方程：

$$\frac{d^2u}{d\theta^2} + u = \frac{\mu k}{L^2}$$

這是一個常係數非齊次線性微分方程，通解為

$$u(\theta) = \frac{\mu k}{L^2}\left[1 + e\cos(\theta-\theta_0)\right]$$

取 $\theta_0=0$ 並令 $p=L^2/(\mu k)$，得到圓錐曲線的標準極式：

$$r(\theta) = \frac{p}{1+e\cos\theta}$$

其中 $p$ 為半通徑、$e$ 為離心率。當 $0\le e<1$ 時軌道為橢圓，力心位於其中一個焦點——這就是克卜勒第一定律。$e=0$ 為圓、$e=1$ 為拋物線、$e>1$ 為雙曲線（如某些彗星或星際天體）。第一定律之所以是橢圓而非其他曲線，本質上來自「平方反比」這個特殊指數。

能量、離心率與活力公式

軌道的能量積分為

$$E = \frac{1}{2}\mu\dot{\vec{r}}^2 - \frac{k}{r} = \text{常數}$$

把總能量、角動量與軌道幾何聯繫起來，可推得離心率與能量的關係

$$e = \sqrt{1 + \frac{2EL^2}{\mu k^2}}$$

由此看出能量的符號決定軌道類型：$E<0$ 為束縛橢圓，$E=0$ 為拋物線，$E>0$ 為雙曲線。對橢圓而言，半長軸 $a$ 僅由能量決定：

$$E = -\frac{k}{2a}$$

結合能量守恆，可整理出軌道力學中最實用的工具——活力公式（vis-viva equation）：

$$v^2 = GM\left(\frac{2}{r}-\frac{1}{a}\right)$$

它直接給出軌道上任一點的速率，是計算近地點、遠地點速度與設計軌道轉移的基礎。

週期與第三定律

對橢圓，整個軌道掃過的面積為 $A=\pi a b$，其中 $b=a\sqrt{1-e^2}$。由等面積定律 $dA/dt=L/(2\mu)$ 為常數，週期

$$T = \frac{A}{dA/dt} = \frac{2\mu\,\pi a b}{L}$$

代入 $b=\sqrt{ap}$、$p=L^2/(\mu k)$ 與 $k=GMm$、$\mu\approx m$（當 $m\ll M$），化簡得

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

這就是克卜勒第三定律的嚴格形式。注意牛頓的版本比克卜勒原始版本多了 $(M+m)$ 修正項——這個修正讓我們能藉由觀測週期與軌道反推天體質量，是測量恆星、行星乃至系外行星質量的核心方法。

定量小範例：低地球軌道的環繞速度與週期

考慮一顆衛星在地表上方 $400\text{ km}$ 的近圓軌道。取地球質量 $M=5.972\times10^{24}\text{ kg}$、半徑 $R=6.371\times10^6\text{ m}$，故軌道半徑

$$r = a = R + h = 6.371\times10^6 + 4.00\times10^5 = 6.771\times10^6\text{ m}$$

圓軌道下 $e=0$、$r=a$，活力公式退化為 $v^2=GM/r$：

$$v = \sqrt{\frac{GM}{r}} = \sqrt{\frac{6.674\times10^{-11}\times5.972\times10^{24}}{6.771\times10^6}} \approx 7.67\times10^{3}\ \text{m/s}$$

約 $7.67\text{ km/s}$。再由第三定律求週期：

$$T = 2\pi\sqrt{\frac{a^3}{GM}} = 2\pi\sqrt{\frac{(6.771\times10^6)^3}{6.674\times10^{-11}\times5.972\times10^{24}}} \approx 5.55\times10^{3}\ \text{s}$$

即約 $92.4$ 分鐘——與國際太空站約 $93$ 分鐘的繞行週期一致，差異主要來自其略高的實際軌道高度。這個小計算同時用到了活力公式與克卜勒第三定律，展示了從幾何到動力學的完整鏈條。

深入探討（研究所視角）

把克卜勒問題放進分析力學，會看到更深的結構。在拉格朗日表述中，平面中心力問題的拉格朗日量為 $L=\tfrac12\mu(\dot r^2+r^2\dot\theta^2)+k/r$，由於 $\theta$ 為循環座標（cyclic coordinate），其共軛動量 $p_\theta=\mu r^2\dot\theta$ 守恆，這正是角動量守恆的拉格朗日語言。轉到哈密頓表述，$H=\dfrac{p_r^2}{2\mu}+\dfrac{p_\theta^2}{2\mu r^2}-\dfrac{k}{r}$，其中第二項即為「離心位能」，使有效位能 $V_{\text{eff}}(r)$ 出現位能井，束縛軌道對應於井底附近的振盪。

克卜勒問題的真正特殊之處在於它是最大可積（maximally superintegrable）系統：除了能量與角動量，還存在第三個守恆向量——Laplace–Runge–Lenz 向量 $\vec{A}=\vec{p}\times\vec{L}-\mu k\,\hat r$。它指向近日點方向且大小正比於 $e$，其守恆正是「橢圓軌道不進動」的根源。一旦力場偏離嚴格的 $1/r^2$（例如廣義相對論修正、或其他天體的微擾），LRL 向量不再守恆，軌道便緩慢進動——水星近日點每世紀約 $43''$ 的反常進動，正是這個機制最著名的觀測證據，也是廣義相對論的早期勝利之一。

從對稱性看，這個額外守恆量對應一個隱藏的動力學對稱：束縛態的克卜勒問題具有 $SO(4)$ 對稱群（Fock 1935），其量子版本正是氫原子的能階簡併，能解釋為何主量子數 $n$ 相同而角動量 $\ell$ 不同的態能量相同。換言之，行星軌道的封閉性與氫原子光譜的簡併，源自同一個數學結構——這是中心力問題最迷人的跨領域連結。

在現代軌道力學中，二體解只是零階近似。真實太空任務要處理多體攝動、地球扁率（$J_2$ 項）、大氣阻力與輻射壓，常用方法是把軌道要素 $(a,e,i,\Omega,\omega,\nu)$ 視為慢變量，透過 Lagrange 或 Gauss 行星方程描述其演化。三體問題雖無一般解析解，但限制性三體問題的 Lagrange 點與 Jacobi 積分，已成為今日韋伯太空望遠鏡駐留 $L_2$、以及低能量轉移軌道（利用弱穩定邊界）設計的理論基礎。從克卜勒的橢圓到混沌動力學的不變流形，同一條運動方程不斷被重新詮釋。