單擺運動
林琦焜
『Philosophy is written in this grand book — I mean the universe — which stands continually open to our gaze, but it cannot be understood unless one first learns to comprehend the language and interpret the characters in which it is written. It is written in the language of Mathematics, and its characters are triangles, circles, and other geometric figures, without which it is humanly impossible to understand a single word of it.』
— Galileo Galilei, Il Saggiatore (1623) —
§1 前言
近代科學之父伽利略 (Galileo Galilei 1564-1642), 能夠經由數學的運用, 將人類的推理 能力應用於自然界, 這是人類理智精華之結晶。 現在我們知道自然界的定律, 可以經過實驗, 隔 離某些不重要的因素, 再細心的觀察, 與數學的推導來加以了解, 這是科學家的世界觀, 也曾是 伽 利略表達得非常清楚的世界觀。 近代科學之所以成功, 主要在於科學活動中採取伽利略所揭 示的 定量方法來描述各種現象以取代神學哲學與玄學的解釋, 這也正是他超越古希臘眾博學多 才之士的地方 。
在所有數學領域中, 就屬微分方程 (differential equation) 與大自然的關係最密切, 法國 偉 大數學家 Henri Poincare(1854-1912) 在“科學的價值”一書中就曾說:
『The science of physics does not only give us (mathematicians) an opportu- nity to solve problems, but helps us to discover the means of solving them, and it does this in two ways: it leads us to anticipate the solution and suggests suitable lines of argument. 』
我 個人對於微分方程這門學問的 philosophy 正如英國物理學家 Dirac 所言真正理解物理問題 的意思是 不用解方程式就看出答案是甚麼。 “如果你相信微分方程是在描述大自然的現象, 那麼, 在你還沒有解方程式之前, 它就應該透露你一些秘密。” 學數學最好是從例子(example) 開始。
32
由於 個人的偏好, 所以選取單擺運動作為這篇文章所要闡述的對象, 聽說伽利略少年時因為教 堂的崇拜儀式十分乏味使人生厭, 轉而注視教堂吊燈之運動並發現:單擺完成一次擺動與其振幅 無關 。 這個方程式直接是牛頓第二運動定律之推論, 與虎克定律有相通處, 甚至到二十世紀的量 子力學也出現這個方程式, 因為它們都是在描述波(wave) 的現象。 所有的微分方程都是一種逼 近 (approximation), 因此要真正了解單擺運動則有賴於橢圓函數 (elliptic function)。 這門 十九世紀重要的學問一直到 1970年代的孤粒子理論 (soliton) 可積系統 (integrable system) 才發揚 光大, 我們在最後一節簡單提一下這個理論, 更深入的內容則留給讀者作更進一步的探 討。
§2 鐘擺之運動方程式
單擺是一理想化物體, 假設有一質點, 以不能伸長之細繩 (表示細繩之質量可以 忽 略) 懸之。 然後將擺錘拉往平衡位置之 一側後釋放, 單擺由於重力的影響而作左 右來回振盪之運動, 試問其運動方程式?
由 力的分解 ( 也就是平行四邊形法 則), 可以將擺錘之重力 mg 分解成法向量 mg cos θ與切向量mg sin θ 兩部分, 法向 量與細繩之張力互相平衡, 故唯一有作用 的 力是 mg sin θ , 因此回復力為
.
... ...
... .. ...
... ...
... .. ...
.. ...
... ...
... .. ...
.. ...
... .
.. . .. .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. ..
L
..
... .
...
S
.
...
•
........................... ...
.. .. .. .. .. . .. ..
. .. .. .. .. .. .. .. .. . .. . .. . .. .. .. .. .. .. .. .. . . .. . .. .. .. .. .. .. .. .. .. .. . .. .. .. .. .. .. .. .. .. .. . . .. .. . . . .. ... . ...
mg cos θ
. .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. .. . .. . .. .. . .. . . .. .. . .. .. . .. .. . .. .. . .. ...
mg
. ...
mg sin θ
...
θ
...
θ ◦
. ...
θ
圖一
F = −mg sin θ (2.1)
另 一方面擺錘之加速度可以這麼看, 因為位移就是弧長 x = Lθ , 故加速度 a 等於 a = d 2 x
dt 2 = d 2 (Lθ)
dt 2 = L d 2 θ
dt 2 (2.2)
由牛頓第二運動定律 (F = ma) 可以導出單擺之運動方程式 mL d 2 θ
dt 2 = −mg sin θ 或 d 2 θ dt 2 + g
L sin θ = 0 (2.3)
它 表示所考慮的量 θ 是隨時間 t 按規律 (2.3) 而變化, 這是一個二階非線性常微分方程 (sec-
ond order nonlinear ordinary differential equation)。 這個方程式看起來簡單, 但實際上卻
有很深 刻的數學內涵。 值得提醒的是方程式 (2.1) 告訴我們: 回復力 F 不與鐘擺移動的角度 θ
34 27 2 92 6
成正比, 而是與正弦 sin θ 成正比。 因此單擺運動並不是簡諧運動 (simple harmonic motion)。
但若擺動的角度 θ 很小, θ ≪ 1 , 單擺運動近乎直線, 則可視為簡諧運動, (sin θ ≈ θ)
F = −mgθ = −mg Lθ
L = − mg
L x (2.4)
故擺 動角度甚小時, 0 < θ ≪ 1 , 回復力與位移 x 成正比而方向相反, 這正是虎克定律 (也就 是簡諧運動), 方程式 (2.3) 則簡化為
d 2 θ dt 2 + g
L θ = 0 或 d 2 x dt 2 + g
L x = 0 (2.5)
因此方程式 (2.5) 的解是原鐘擺方程式 (2.3) 的近似解 (approximate solution)。 由 Taylor 展開式來看
sin θ = θ 1! − θ 3
3! + θ 5
5! − . . . (2.6)
所 以對充分小的 θ , 即單擺之擺動角度很小時方程式 (2.5) 是單擺運動方程 (2.3) 很好的近似。
§3 單擺方程的解
如何解二階微分方程 (2.5)? 這是一門很有趣的學問, 我們將從幾個方向來看這個問題。
3.1 觀察與猜測:
將 (2.5) 移項寫成
d 2 x
dt 2 = − g
L x (3.1)
這方程式告訴我們, 要找一個函數 x(t) 其二次微分, 除了常數 L g 外, 等於原函數之負值, 而且 單擺運動有最高與最低點, 這相當於函數 x(t) 有上界與下界。 由微積分的知識知道具有這樣性 質的就是正弦或餘弦函數。 因此可以假設 (2.5) 之解為
x(t) = A cos(ωt + δ) (3.2) 因為
cos(ωt + δ) = cos δ cos ωt − sin δ sin ωt = a cos ωt + b sin ωt 所 以常數 δ 之存在以允許 (3.2) 是正弦、 餘弦函數之組合。(當然 (3.2) 也可寫成 x(t) = A sin(ωt + δ) !) 微分兩次後代回 (3.1):
−ω 2 A cos(ωt + δ) = − g
L A cos(ωt + δ)
若取 ω 2 = L g , 則方程式 (3.1) 的解為
x(t) = A cos(ωt + δ), ω 2 = g
L (3.3)
A 與 δ 尚未決定, 為任意常數, 意即 (3.1) 有無窮多解, 至於為何會有兩個參數 A, δ , 那是完 全自然的, 因為 (3.1) 本來就是二階微分方程有兩組獨立解!
(3.2) 這個猜測可以換為
x(t) = e mt (3.2 ′ )
這理由是因為指數函數的任意次微分仍然是指數函數。 將 (3.2 ′ ) 代回 (2.5) m 2 + g
L = 0 ⇒ m = ±
r g L i 所 以 x(t) = e ± i √
g/Lt 由 Euler 公式取實部與虛部 (因為原方程式是實數值!) 因此 x(t) = sin q g/L t 或 cos q g/L t , 因為方程式是線性 (linear), 所以一般解是這兩個獨立解的線性 組合
x(t) = a cos
r g
L t + b sin
r g
L t (3.3 ′ )
3.2 降階法:
我們所關心的是鐘擺的擺動, 其速率是角速度, 令 v = ˙x = dx
dt (3.4)
則 (3.1) 成為
˙v = dv
dt = d 2 x
dt 2 = − g L x 但另一方面由連鎖律
dv dt = dv
dx dx
dt = v dv
dx (3.5)
因此我們將方程式 (3.1) 轉換成一階非線性微分方程 v dv
dx = − g
L x 或 vdv + g
L xdx = 0 (3.6)
這是一個全微分 (total differential) d
1
2 v 2 + 1 2 g L x 2
= 0 可以積分
1 2 v 2 + 1
2 g
L x 2 = C (動能 + 位能 = 常數) (3.7)
36 27 2 92 6
其中 C 是一常數, 這個等式之本質就是能量守恆律或者寫成 1
2 ( dx dt ) 2 + 1
2 g
L x 2 = C (3.8)
我 們可以透過微分的相反運算 (積分) 來解 (3.8), 為著方便, 我們令 C = 1 2 L g A 2 (事後孔明!) 開根號
dx dt =
r g L
√ A 2 − x 2 , −A ≤ x ≤ A (3.9)
這個方程的典型解法就是分離變數法 ( x 歸 x , t 歸 t , 凱撒的歸凱撒, 上帝的歸上帝)
Z dx
√ A 2 − x 2 =
r g L
Z
dt (3.10)
左邊這個積分正是反三角函數 sin −1
x A
=
r g
L t + C ′ 或 cos −1
x A
=
r g L t + C ′ 所以
x = A sin
r g L t + C ′
, x = A cos
r g L t + C ′
(3.11) 與 (3.3) 完全吻合。 而且由 (3.9) 更可得知 x 之最大值等於 A, 所以 A 就是振幅 (amplitude) 由初始值即最開始的高度所決定。
降階法的精神在於引進速度 v = dx dt 這個新變數, 在微分方程理論我們將 (x, v) 平面稱為 相空間 (phase space)。 因此二階微分方程 (2.5) 可轉換為一階聯立微分方程組
d 2 x dt 2 + g
L x = 0 ⇐⇒
( dx
dt = v
dv
dt = − L g x (3.12) 最後這兩式相除 (將變數 t 隱藏) 就是 (3.6)
dv
dx = − g L
x
v (3.13)
也因此可以在 (x, v) 平面上研究原微分方程 (2.5)。 我們稱 (x, v) 平面為相空間 (phase space)
由能量守恆來看, 可以更清楚明白為何相空間這個觀念對於研究微分方程會如此重要。
3.3 能量守恆:
既然方程式 (2.5) 是由牛頓定律而來, 我們就理所當然從力學的角度來思考。 首先介紹幾 個基本物理量
x : 位置函數 v = dx
dt : 速度 a = dv
dt = d 2 x
dt 2 : 加速度 T = 1
2 v 2 = 1 2 ( dx
dt ) 2 : 動能(質量視為1) 除此之外, 力與位能函數之關係 F = − dU dx , 因此 (2.5) 可以改寫為
d 2 x dt 2 = − g
L x = − d dx
1 2
g L x 2
= − d
dx U (3.14)
其中 U = 1 2 L g x 2 是 位能, 而 E = T + U = 1 2
dx dt
2
+ 1 2 L g x 2 就是總能量, 我們在 (3.7) 就已 推導得這個能量守恆律。
定理 :(能量守恆)
E(t) = 1 2
dx dt
2
+ 1 2
g
L x 2 = C(常數) (3.15) 證明 :直接微分
dE dt = dx
dt d 2 x
dt 2 + g L x dx
dt
= d dt
d 2 x dt 2 + g
L x
= 0
特 別注意的是, 整個證明過程並不依賴微分方程的精確解 (exact solution), 微分方程本身就 已經告訴我們這個秘密了。 從證明的過程反推回去, 如果直接乘 dx dt , 也可以得到能量守恆律, 並且分別得到動能 T 與位能 U
dx dt ( d 2 x
dt 2 + g
L x) = d dt
1 2
dx dt
2
+ 1 2
g L x 2
= d
dt (T + U)
在推導過程中, 也清楚看到動能 T 是由加速度而來, 位能 U 則是力的結果 d 2 x
dt 2 + g
L x = 0
↓ ↓
T U
38 27 2 92 6
方程式 (3.15) 可以視為古典力 學的 Hamilton-Jacobi 方程式, 它 除了說明能量守恆之外, 由於 C 是 任意常數, 它更代表了等位線 (level curve), 所以將 v = dx dt 視為另一個 變數, 則
E(t) = E(x, v) = 1 2 v 2 + 1
2 g L x 2
= C ≥ 0
其圖形在 (x, v) 平面是橢圓族, 若 g = L 則是圓。
.
...
x
. .. .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . .. .. . .. . .. .. .. . .. .. . .. .. . .. .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. . .. .. . ... .. . .. .. . .. .. .. .
...
U(x) = 1 2 L g x 2
0
.
...
x
.. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. . .. .. . ... .. . .. .. . .. .. .. .
v
.. .
... .
.. ..
(0, 0)
( q 2 LC g , 0) E(t) = C
. ...
. .. .. . .. . .. . . ..
•
..................... . .. .. .. .. . .. . ..•
.............................. ...
. . .. .. . .. . .. . ..
. .. .. .. . .. . .. . .. .. .. . .. .. . .. . .. .. . .. .. . .. . .. .. . .. .. . .. .. . .. . .. .. . .. .. . .. .. . .. . .. .. . .. .. . .. . .. .. . .. .. . .. . .. .. . .. .. . .. . .. .. .. . .. .. . . .. .. .. . .. .. .. . .. . .. . .. .. .. . .. . .. .. . .. ..
...
. .. .. . .. .. . .. . .. .. . .. .. . .. . .. .. . .. .. . .. .. . .. . .. .. . .. .. . .. .. . .. . .. .. . .. .. . .. .. . .. . .. .. . .. .. . .. .. . .. . .. .. . .. .. . .. .. . .. . .. .. . .. .. .. . . .. .. . .. .. .. . . .. .. .. . .. .. . . ..
.. .. .. . .. . . .. .. .. .. . .. .. . .. . .. .. . .. .. . .. .. . .. . .. .. . .. .. . .. .. . .. . .. .. . .. .. . .. .. . .. . .. .. . .. .. . .. .. . .. . .. .. . .. .. . .. .. . .. . .. .. . .. .. .
圖二
3.4 幾何觀點 — 相空間:
從幾 何的角度來探討微分方程這個觀念可以追溯至法國數學家H. Poincare(1854-1912) 與蘇俄數學家 A. Liapunov。 我們考慮比 (3.12) 更一般的方程式
dx
dt = ˙x = p(x, v) , dv
dt = ˙v = q(x, v) (3.16) 在物理上 (x, v) 代表一系統的狀態即位置 (position) 與速度 (velocity), 這個平面也稱為相 平面 (phase plane), 因為是描述狀態之變化例如詩人所言 「月有陰晴圓缺」就是在描述月相 (moon phase) 之變化。
... . . . .. . . . .. . . . .
...............
(x, v)
. .. . . .. . . .. . . .. . .. . . .. . . .. . . .. . . .. . .. . . . .. . .. . . .. . . .. . . .. . .. . . .. . . .. . . .. . . .. . .. . . . .. . .. . . .. . . .. . . .. . . .. . .. . . .. . . .. . . .. . .. . . . .. . .. . . .. . . .. . . .. . . .. . .. . . .. . . .. . . .. . .. . . .. . . . .. . .. . . .. . . .. . . .. . .. . . .. . . .. . . .. . . .. . .. . . . .. . .. . . .. . . .. . . .. . .. . . .. . . .. . . .. . . . . ... . ...
.. . ...
. .. .. .. . .. .. .. .
(x + ∆x, v + ∆v)
(∆x, ∆v)
. .. . .. .. . .. ...
.. .. . .. .. . .. ... .
...
... . . . .. . . . .. . . . . .
(x, v)
.. . .. . .. . .. .. .. .. .. .. .. .. . . .. . .. .. .. .. .. .. .. .. . .. .. .. .. .. .. .. .. .. .. . .. .. .. .. .. .. .. .. .. .. . .. .. .. .. .. .. .. .. . . .. .. . .. .. .. .. .. .. .. .. . . .. . .. .. .. .. .. .. .. ...
...
(dx, dv)
...
圖三
因為 (dx, dv) 代表 (x, v) -平面之切向量 (tangent vector), 如果將微分方程解之曲線
視為水流 (flow), 則 (3.16) 告訴我們在每一點的速度 V = (p, q) 其水平分量是 p(x, v) 垂直
分量是 q(x, v), 透過向量場 V = (p, q) 我們可以看出微分方程之積分曲線, 將 (3.16) 兩式相 除 (參考 (3.13))
dv
dx = q(x, v)
p(x, v) ≡ f(x, v) (3.17) 這個方程式告訴我們曲線 v 在該點之切線斜率。 這正是微分方程之義意 — 未知函數與其微分 滿足某關係式 (等式), 求此函數? 從歷史而言常微分方程是伴隨微積分發展起來的分析之一分 支, 如果常微分方程可以被積分, 則由微積分基本定理就達到目的, 也因此這方法稱為面積解法 (quadrature), 但是可以積分的微分方程實在是少之又少。 直到 1881-1886 年間 Poincare 將 常微 分方程的研究方法從分析方法轉化為更直觀的幾何方法, 從此開創了常微分方程定性理論 這一學科, 並為後來的動態系統 (dynamical system) 預測了美好的遠景。
回 到 (3.12) 或 (3.13) (dx, dv) = (v, −gx/L), dv
dx = − gx Lv 如果你覺得會困惑, 則不妨先將 g/L 視為 1, (dx, dv) = (v, −x) 與向量 (x, v) 比 較正好 x, v 互換差一個負號:
(dx, dv) ⊥ (x, v) 平面上具有這性質的圖 形 就是圓, 同理可推平面上任意向量 (x, v) 其切向量為 (dx, dv) = (v,−gx/L) 的 圖形必定是橢圓 (圓的變形), 或者根據 方程式可以描繪各點之向量場, 則隱約可
.
...
.. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . .. . .. .. . .. .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. .. . .. . . .. .. . ... . .. .. . .. .. .. . ..
.. .. .. .
... .
.. .. . .. .. .. .. .
... .
.. .. . ..
. .. . .. .. . .. .. .. . .. . .. .. . .. .. .. . .. . .. .. .. . .. .. . .. . .. .. . .. .. .. . .. . .. .. .. . .. .. . .. . .. . .. . .. .. .. . .. .. . .. . ...
(0, − L g )
.. . .. . .. .. . .. .. .. . .. . .. .. .. . .. .. . .. . .. .. . .. .. .. . .. . .. .. .. . .. .. . .. . .. .. .. . .. .. . . . .. . .. ... . .. . .. .. .. . .. ..
(0, L g )
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ... . ...
(−1, 0)
. ...