天體力學中的橢圓積分：克卜勒問題為什麼是例外？

\(\require{physics}\)

上一篇從橢圓弧長推出了第二類橢圓積分 $E(k)$，順便把「橢圓積分」這個名字的來源講清楚了：它源自橢圓的周長 $C=4aE(e)$。

但如果你讀過克卜勒系列第五篇，可能會覺得奇怪——那篇文章的結論是：行星沿橢圓軌道運行的週期，完全不需要橢圓積分，只靠一個三角代換就化成初等函數，得到克卜勒方程 $M=\psi-e\sin\psi$。

那麼，「橢圓積分」跟「橢圓軌道」到底有沒有關係？天體力學裡真的用得到橢圓積分嗎？

答案是：克卜勒問題本身是一個特例。只要稍微偏離這個特例——換一個量（弧長而不是時間）、換一個力場（不是純 $1/r$）、或換一個系統（不是兩體）——橢圓積分（甚至更複雜的特殊函數）就會冒出來。這篇文章把這些地方整理出來。

---

一、回顧：克卜勒問題為什麼是初等的

克卜勒系列第五篇把週期寫成積分

\[T=2\int_{r_{\min}}^{r_{\max}}\frac{\dd{r}}{\sqrt{2(\varepsilon-V_{\text{eff}}(r))}},\qquad V_{\text{eff}}(r)=\frac{h^2}{2r^2}-\frac{GM}{r}\]

展開後配方，得到

\[\dot r^2=\frac{GM}{ar^2}\left[a^2e^2-(r-a)^2\right]\]

關鍵的一步是離心近點角代換 $r=a(1-e\cos\psi)$，它讓 $a^2e^2-(r-a)^2=a^2e^2\sin^2\psi$，根號裡的東西變成完全平方，根號直接開掉，積分變成初等函數 $\psi-e\sin\psi$。

對比單擺的橢圓積分：能量守恆給出 $\dot\theta^2\propto\cos\theta-\cos\theta_0$，半角代換後根號裡剩下 $1-k^2\sin^2\phi$——這個東西不是完全平方，無法再開根，這就是 $K(k)$ 誕生的地方。

兩者外觀相似（都是「能量減位能」開根號再積分），但克卜勒問題的有效位能恰好是 $1/r$ 與 $1/r^2$ 的組合，這種特殊的代數結構讓配方之後自動出現完全平方。這不是運氣，下一節會說明背後的原因。

---

二、橢圓本身：軌道弧長與 $E(k)$ 的回歸

克卜勒方程處理的是「時間」這個量。但如果問的是另一個量——行星沿軌道走過的實際距離（弧長）——答案就完全不同了。

行星軌道在空間中就是一個離心率為 $e$ 的橢圓，半長軸 $a$。上一篇已經證明任意橢圓的周長是

\[C=4a\,E(e)\]

所以：行星繞太陽一圈，時間是初等公式 $T=2\pi\sqrt{a^3/GM}$，但走過的路程是 $C=4aE(e)$，需要第二類橢圓積分。

這是一個很好的對照組：同一個橢圓軌道，問「花多久」是初等問題，問「走多遠」卻是橢圓積分問題。兩個問題的差別在於：時間的積分變數是離心近點角 $\psi$（克卜勒方程的代換對象），而弧長的積分變數是真正的幾何弧長元素 $\dd{s}=\sqrt{\dot x^2+\dot y^2}\,\dd{t}$——這正是 $E(k)$ 定義時的那個積分，跟克卜勒問題的代數結構無關，純粹是「橢圓是橢圓」這個幾何事實。

---

三、其他會出現橢圓積分的天體力學問題

除了弧長，天體力學裡還有幾類經典問題會用到橢圓積分或更高階的特殊函數。

3.1 扁球重力場

地球不是完美的球體，而是扁球體（oblate spheroid）——赤道略凸、兩極略扁，這就是地球的 $J_2$ 項。計算扁球體（或更一般的三軸橢球體）外部的重力位能時，積分中會出現橢圓積分。這也是衛星軌道攝動分析中最基本的修正項之一。

3.2 受擾動的軌道

真實太陽系不是乾淨的二體問題：其他行星的引力、行星自身的扁率、潮汐效應都會對軌道造成微小擾動。這些擾動項展開後，往往不再具有克卜勒問題那種「恰好完全平方」的結構，計算擾動積分時自然會遇到橢圓積分或更高階的超橢圓積分。

3.3 限制性三體問題

三體問題本身不可解析求解，但在某些特殊情形（例如限制性三體問題中的特定能量曲線、週期軌道族）會出現橢圓積分或橢圓函數，用來描述軌道的週期與形狀。

3.4 廣義相對論：近日點進動

在廣義相對論中，行星軌道由 Schwarzschild 時空中的測地線方程描述。徑向方程的形式是

\[\left(\frac{\dd{r}}{\dd{\tau}}\right)^2 = (\text{關於 } r \text{ 的三次或四次多項式})\times\frac{1}{r^4}\]

牛頓力學的徑向方程是二次多項式（也就是前面那個「完全平方」的來源），廣義相對論多了一個 $1/r^3$ 的修正項，把多項式的次數從二次推到三次。二次多項式開根號可以用三角代換消掉；三次或四次多項式開根號，一般就是橢圓積分的定義式本身。 這正是水星近日點進動這類問題在精確處理時會碰到橢圓積分的原因（雖然實務上常用攝動展開取近似值）。

---

四、共通結構：對稱性決定一切

把以上幾個例子放在一起看，會發現一個共通的模式：

問題	有效位能 / 徑向方程	根號內	結果
克卜勒週期（時間）	$1/r,\ 1/r^2$ 組合 → 二次多項式	完全平方 $\sin^2\psi$	初等函數
橢圓弧長	（與位能無關，純幾何）	$1-e^2\sin^2\theta$	$E(k)$
扁球重力場	多了高階項	三次以上	橢圓積分
GR 近日點進動	多了 $1/r^3$ 項 → 三次多項式	三次多項式開根	橢圓積分

克卜勒問題的特殊性，本質上來自牛頓重力 $V\propto 1/r$ 這個特定的冪次。 根據 Bertrand 定理，在所有中心力場中，只有 $V\propto 1/r$（克卜勒）和 $V\propto r^2$（簡諧振子）這兩種情況，所有有界軌道都是封閉曲線——這背後藏著一個額外的守恆量，Laplace–Runge–Lenz 向量（克卜勒系列第五篇結尾已經預告過）。這個額外的對稱性，正是讓週期積分的根號自動變成完全平方的深層原因。

只要偏離 $1/r$（多一個 $1/r^2$ 的扁率項、$1/r^3$ 的相對論修正、或來自第三個天體的擾動），這份「多出來的對稱性」就被打破，根號下的多項式次數就會升高，橢圓積分（或更高階的特殊函數）就會出現。

---

五、銜接下一篇：轉動也有自己的「克卜勒問題」

到目前為止，我們看到的橢圓積分都跟軌道有關——位置 $r$ 隨時間或角度的變化。但天體力學還有另一類運動：自轉。

一個不受外力矩的剛體（例如一塊不規則形狀、自由翻滾的小行星），它的自轉軸方向會隨時間演化，由歐拉方程式（Euler’s equations）描述。這個系統同樣有兩個保守量（能量與角動量），結構上跟克卜勒問題的「兩個保守量決定一條軌道」非常相似——但它的解不是初等函數，而是上一篇 $K(k)$、這一篇 $E(k)$ 的「親戚」：雅可比橢圓函數 $\mathrm{sn},\mathrm{cn},\mathrm{dn}$。

真實世界中，像是翻滾的小行星（例如 4179 Toutatis）、彗星核（67P/Churyumov–Gerasimenko）的非主軸自轉，都可以用這套語言描述。下一篇將正式介紹 $\mathrm{sn},\mathrm{cn},\mathrm{dn}$——它們是 $K(k)$ 積分的反函數，也是整個「單擺與橢圓積分」系列從第一篇 $K(k)$ 出發，最終要抵達的目的地。

---

小結

量	是否需要橢圓積分/函數	原因
克卜勒週期（時間）	❌ 初等	$1/r,1/r^2$ 結構 → 完全平方（Bertrand 定理的特殊性）
橢圓軌道弧長	✅ $E(e)$	橢圓弧長的幾何本質，與動力學無關
扁球重力場	✅	重力位能展開出現高階項
受擾動軌道	✅（常見）	擾動破壞完全平方結構
GR 近日點進動	✅	徑向方程多出 $1/r^3$，多項式次數升高
剛體自由轉動	✅ Jacobi $\mathrm{sn,cn,dn}$	下一篇

一句話：「橢圓軌道」和「橢圓積分」這兩個「橢圓」其實是兩件不同的事——前者是牛頓力學在 $1/r$ 力場下的幾何結果，後者是橢圓弧長積分的特殊函數。克卜勒問題剛好讓兩者「脫鉤」，但只要系統稍微複雜一點，它們就會重新糾纏在一起。

還有問題嗎？
→ 歡迎預約家教課，直接針對你的問題討論