橢圓的曲率公式和萬有引力的 平方反比規律
張海潮 · 莊正良
一、 導言
牛頓 (1642∼1727) 在 1687 年出版 《自然哲學之數學原理》 (或稱 《原理》)。 牛頓先在書 中第一卷第二章的命題 1 和 2 證明了克卜勒的面積律等價於行星繞日是因為受到太陽的向心 (吸引) 力。 接著在同卷第三章的命題 11 證明從克卜勒的橢圓律可以推得行星受太陽的吸引力 (或向心加速度) 與行星到太陽距離的平方成反比, 比例常數是 4π2a3
p2, 式中 a 是橢圓軌道的半 長軸, p 是行星繞日的週期。 由於克卜勒的週期律說 a3
p2 對太陽系中任一行星均為定值, 因此而 有萬有引力之萬有一詞 (註一) 。
本文最主要的目的在從橢圓的曲率公式重新說明橢圓軌道與距離平方反比規律的關聯, 有 別於牛頓在 《原理》 中的方法, 見本文註二及三、 四節。
二、 平面曲線的曲率、 曲率半徑及物理意義
本節定義平面曲線的曲率, 並透過計算解釋曲率半徑的物理意義, 如圖一:
圖一
從 P 到 Q, 曲線的切線方向代表質點瞬間的速度方向, 切線之間的夾角對弧長的變率即 為曲率, 以半徑為 R 的圓為例, 如圖二:
圓從 P 點循逆時針回到 P 點, 角度總變化量是 2π, 經歷的弧長是 2πR, 兩者相除而得
24
圖二
曲率等於 1
R, 此為對稱的情形。 若是一般曲線, 曲率會逐點變化; 今以 κ 表曲率, 仿圓的情形, 另定義曲率半徑 ρ 為 1
κ。
現以 (x(t), y(t)) 表平面參數曲線, 速度 (切) 向量為 V = dx dt,dy
dt
, 加速度向量為 A = dV
dt =d2x dt2,d2y
dt2
。 弧長微量單位
ds2 = dx2+ dy2 或 ds dt
2
=dx dt
2
+dy dt
2
= V · V 或 ds
dt = |V | = 向量 V 的長度 (1) 今考慮單位切向量 T = V
|V |, 並以 (cos θ, sin θ) 表 T , θ 即 T 的方向 (角度) , 曲率 κ 即 dθ ds。 T = (cos θ, sin θ)
dT ds =
− sin θdθ
ds, cos θdθ ds
亦即
dT ds =
dθ
ds = κ (取正) (2)
將 T = V
|V | 代入,
dT
ds = d V
|V |
.ds = |V |dVds − V d|V |ds
|V |2 (3)
但 dV ds = dV
dt · dt
ds = Adt
ds, 由 (1) dV
ds = A/|V | 代入 (3) 得 dT
ds = A − V d|V |ds
|V |2 (4)
注意到 dT
ds 和 T 垂直, 因此 dT
ds 也和 V 垂直, 所以由 (2)、(4)
dT ds × V
= |V |κ =
|A × V |
|V |2 或 κ = |A × V |
|V |3 (5)
接著, 從 κ = |A × V |
|V |3 來看曲率半徑的物理意義 如圖三, V 為速度 (切) 向量, A 為加速度
圖三 向量, A, V 夾角為 ε, 虛線垂直 V , 為法線方向, κ = 1
ρ, ρ 為曲率半徑, 由 (5) ρ = |V |3
|A × V | = |V |3
|A||V | sin ε 或
|A| sin ε = |V |2
ρ (6)
亦即速度的平方除以曲率半徑會等於法線方向的加速度, 等速圓周運動時, ε = 90◦, (6) 相當 於 |A| = |V |2
ρ , ρ 是圓運動的半徑。 (註三)
三、 橢圓曲率與平方反比規律
如圖四, 行星在橢圓上運動, 服從面積律:
圖四
行星受到的向心加速度為 A, 假定 A 服從平方反比規律:
A = C · 1
SP2, C 為比例常數 兩邊同乘 sin ε, 得
A sin ε = C · sin ε SP2
= C · sin3εv2 SP2v2sin2ε 分母 SP2v2sin2ε 由於面積律, 會是一個常數 D, 因此
A sin ε = C
D · sin3ε · v2 (7) 但是因為由上一節公式 (6), A sin ε = v2
ρ , ρ 為曲率半徑 或
A sin ε = v2κ, κ 為曲率 與 (7) 比較發現, 如果 A 滿足距離平方反比, 則
κ = C
D · sin3ε (8)
反之, 如果知道橢圓的曲率 κ 正比於 sin3ε, 則同樣由上一節公式 (6), 曲率半徑的物理意義, A sin ε = v2
ρ = v2κ ∼ v2sin3ε 因此 A ∼ v2sin2ε =v2sin2εSP2
SP2 同樣由面積律得知, v2sin2εSP2 是常數, 故
A ∼ 1 SP2
此即平方反比規律, 由此可見, 平方反比規律和橢圓律的關係其實是在於對橢圓曲率的掌握, 亦 即曲率與 sin3ε 成正比。 (註四)
四、 橢圓曲率公式的計算
雖然在不同的場合, 我們均看到橢圓曲率公式的證明 (同見註四) κ ∼ sin3ε
但是這些證明多半比較迂迴。 例如在 《古代天文學中的幾何方法 pp.190-192》, 作者利用 (a cos t, b sin t) 這個參數式來作計算, 過程冗長。 在本節中, 我們要從曲率的定義, 光學性質以及焦切距 乘積定理來計算 κ, 我們認為這個辦法比較直接, 並且較能彰顯橢圓曲率的本質。 橢圓的光學性 質, 大家都很熟悉, 不必多言。 以下先介紹焦切距乘積定理。
圖五 焦切距乘積定理:
如圖五, 過橢圓上一點 P 的切線 l, 自焦點 F1 和 F2 分向 l 作垂線段, 其長度分別為 d1, d2, 則 d1d2 = b2, 式中 b 是橢圓的半短軸。 (證明見註五)
我們現在從定義來看橢圓的曲率, 如圖六, F1, F2 為焦點,
圖六
從 P 點走到 Q 點, 切線 l1 和 l2 有一角度差, 此一角差除以弧長P Q , 再取極限 Q → P ,⌢ 即為 P 點之曲率。 但 l1, l2 之角差可以用在 Q, P 點之法線 (與切線垂直) 角差來計算。 由光 學性質可知, 此法線實乃相關角度之分角線, 如圖七。
圖七
兩條虛線分別為 P, Q 角之分角線, 則
π = x + γ + ε2+ e + δ π = α + 2γ + ε2+ e π = β + 2δ + e + ε2
由上三式可見 2x = α + β 或 x = 12(α + β), x 是兩分角線之角差, 亦為過 P, Q 兩切線之角 差。 我們要計算 limQ→P x
P Q⌢ = κ, P Q 表 P , Q 間的弧長。⌢
圖七中的角 e 和角 ε2 在 Q → P 時, 均趨近圖六中的角 ε。 如圖七, 在 △F1P Q 中, 由 正弦定律, 有
sin α
P Q = sin e P F1
(9) 而在 △F2P Q 中, 亦由正弦定律, 有
sin β
P Q = sin(2δ + e) P F2
(10) (9)+(10)
sin α + sin β
P Q = sin e P F1
+sin(2δ + e) P F2
(11) 由於
sin α + sin β
α + β = 2 sin α+β2 cosα−β2
α + β 和 P Q P Q⌢
在 Q → P (x → 0, 或 α, β → 0) 時均趨近於 1, 而同時圖七的角 e 趨近於圖六的角 ε, 角 2δ + e 趨近於圖六的角 π − ε, 因此在 Q → P 時
(11) 的極限是
Plim→Q
α + β
⌢
P Q
= sin ε P F1
+sin ε P F2
或
2κ = sin ε 1 P F1
+ 1 P F2
= sin ε · 2a P F1· P F2
或
κ = a sin ε P F1· P F2
, a 為半長軸 (12)
如圖八
圖八
由 P F1 = d1
sin ε 及 P F2 = d2
sin ε, 代入 (12) 得 κ = a sin3ε
d1d2
再由焦切距乘積定理 (註五), d1d2 = b2, 因此得到曲率公式 κ = a
b2 · sin3ε (13)
由第三節, 此一公式可推導出太陽對行星的吸引力服從平方反比的規律。
五、 結語
牛頓在 《原理》 中從橢圓律推出引力的平方反比規律, 並未使用橢圓的曲率公式, 但卻用 了許多橢圓的性質, 例如: 共軛直徑乘積定理。
對橢圓而言, 任一過中心的弦均稱為直徑, 如圖九
圖九
P OQ 為一直徑, L1 和 L2 為過 P, Q 兩點的切線, 則可證得 L1//L2, 並且定過 O 與 L1(L2) 平行之直徑 RS 為直徑 P Q 之共軛直徑。 反之, 亦可證得 P Q 為 RS 之共軛直徑。
所謂共軛直徑乘積定理指的是, 如圖十
圖十
凡與 P Q 之共軛直徑 RS 平行之弦, 如 UV W , 則必有 UV = V W , 且 UV2
P V · V Q 是一定 值, 與 V 之位置無關。 現在利用伸縮線性變換來看此一定理並不困難, 但在牛頓的時代, 想必 難倒許多學者。 牛頓在 《原理》 第一卷第三章命題 11, 從橢圓律推得平方反比的論證, 用了包 括共軛直徑和橢圓其他的性質, 非常難懂。 所以後人多有注釋或另起爐灶企圖簡化牛頓的證明。
其中值得一提的是馬克斯威爾 (Clerk Maxwell, 1831-1879) 在著作 Matter and Motion 一 書的 108, 109 頁利用本文提到焦切距乘積定理 (見註五) 證明了橢圓律可以推得平方反比規 律, 他結論如下 (Matter and Motion, p.109)
“Hence the acceleration of the planet is in the direction of the sun, and is inversely as the square of the distance from the sun.”
另一位終其生寫 《Newton’s Principia for the Common Reader》 的學者 S. Chandrasekhar, (1910-1995), 在書中的 110 頁認為牛頓知道橢圓曲率扮演的角色, 他說 (ρ 為曲率半徑) :
ρ ∼ csc3ε
That Newton must have known this relation requires no argument!
在本文中, 我們並沒有考證牛頓是否知道這個曲率 (半徑) 公式, 而是提出了一個看起來更 直接的證明, 希望讀者欣賞。
註一: 克卜勒 (1571∼1630) 發現的行星三大定律是 1. 橢圓律: 行星繞日的軌道是一橢圓, 太陽位居一焦點。
2. 面積律: 行星繞日時, 在單位時間, 行星與太陽連線段所掃過的面積是一常數。
3. 週期律: 在太陽系中, 任一行星繞日軌道半長軸的立方和繞日週期的平方之比是一常數, 與 個別的行星無關。
本文引用 《原理》 的中譯本 (譯者王克迪, 台北大塊文化出版社)。 第一卷第二章之命題 1, 2 見
《原理》, pp.57-59, 第一卷第三章的命題 11 見 《原理》, pp.71-72。 由於力 F 與加速度 A 成 正比, 所以本文均以 (向心) 加速度 A 代表 (向心) 力, 牛頓所證乃
A = 4π2a3 p2
1 r2
式中 r 為行星與太陽之距離, a 為橢圓軌道之半長軸, p 為繞日之週期。 以今日習用的公式而言 F = mA = 4π2a3
p2m 1
r2 = GMm r2
其中 4π2ap32 換成 GM, G 為萬有引力常數, m, M 分別為行星和太陽之質量。
註二: 關於面積律等價於向心力, 現在的解釋是角動量守恆, 如圖:
S, P 分別是太陽和行星。 SP = ~−→ X, P Q = d ~−→ X 代表微量位移, | ~X × d ~X| 是 ~X 和 d ~X 決定 的外積向量長度, 代表微量掃過面積。 克卜勒的面積律即:
Z H 0
|X × dX| = H · C, H 為經歷的時間, C 為常數 或
Z H 0
X ×~ d ~X dt
dt = H · C 兩邊對 H 微分得
X ×~ d ~X dt
= C 或
| ~X × ~V | = C
式中 ~V 是 P 的速度向量。 由於是平面運動, 所以上式等價於 ~X × ~V 是一常數向量, 此向量 若與質量相乘即為角動量。
若角動量守恆, 將 ~X × ~V 對時間微分, 得 d ~X
dt × ~V + ~X ×d~V dt = ~0 但 d ~X
dt = ~V , d~dtV = ~A (加速度向量) , 因此得 ~X × ~A = 0, 易見受力 ~A 為 ~X 的反方向, 即受 到太陽之向心吸引力。
在牛頓的時代, 普遍相信行星繞日是受到太陽的吸引力, 也同意此一向心吸引力等價於面 積律。 因此下一個任務就是理解此一向心力的大小。 當時確有不少人, 如: 胡克, 猜測此力與距 離的平方成反比, 但是只有牛頓用嚴謹的數學從橢圓的幾何性質推得了平方反比, 此即原理的 第一卷第三章命題 11。
註三: 曲率概念由牛頓提出, 見 《原理》 引理 11, p.53。 在牛頓的時代, 大家先理解了等速圓周 運動的向心加速度公式, A = v2
ρ, 式中 A 為向心加速度, v 為 (等) 速度, ρ 為圓半徑。 大家也 理解到一般運動的加速度有兩個作用, 一是沿切線方向加速, 另一是轉彎, 負責轉彎的加速度就 是加速度在垂直切線方向的投影, 而轉彎的半徑就是曲率半徑, 此即本節所得公式 (6) :
|A| sin ε = |V |2 ρ
註四: 1983 年諾貝爾物理獎得主錢卓斯卡 (S. Chandrasekhar, 1910∼1995) 在 1990 年 以 80 高齡發憤註釋牛頓的 《原理》。 1994 年完稿交由牛津大學於次年出版, 書名 《Newton’s Principia for the Common Reader》 同時錢氏又在 1994 年 《當代科學 Current Sci., 67(7)(1994), 495-496》 發表 《On reading Newton’s Principia at age past eighty》, 在這 些著作中錢氏均提及橢圓的曲率半徑公式, 及此公式與平方反比規律的關聯, 但是錢氏對曲率 公式的證明比較不像本文那麼直接, 詳本文第四節。
註五: 如下圖
從焦點 F1, F2 到切線 L 的垂線段長分別為 d1, d2, 則由橢圓的光學性質看出有下面這個等腰 梯形:
圖中 F1′, F2′ 分別是 F1, F2 對切線 l 的鏡射點, F2, P, F1′ 和 F1, P, F2′ 均滿足三點共線, 圖中 F1′F2 = F1F2′ = 2a, a 為半長軸, F1F2 = F1′F2′ = 2c, c 為焦點到橢圓中心的長度。 令 h 為 此等腰梯形的高, 則有
4a2= F2H2+ h2 = (d1+ d2)2+ h2 4c2= F2′H2+ h2 = (d2− d1)2+ h2 4b2= 4a2− 4c2 = 4d1d2
因此, d1d2 = b2, 式中 b 為橢圓之半短軸, 此即焦切距乘積定理。
參考資料
1. 牛頓。 自然哲學之數學原理。 台北: 大塊文化, 2005。
2. 張海潮, 沈貽婷。 古代天文學中的幾何方法。 台北: 三民書局, 2015。
3. 項武義, 張海潮, 姚珩。 千古之謎。 台北: 商務印書館, 2010。
4. Chandrasekhar, S., Newton’s Principia for the Common Reader, Oxford Univ. Press, Oxford, 1997.
5. Maxwell, C., Matter and Motion, Dover Pub.,New York, [1877] 1991.
—本文作者張海潮及莊正良為台大數學系退休教授—
