一條定律,兩個舞台:單擺與行星
在見過 Lotka-Volterra 與 SIR 模型之後,你已經明白:一個由兩個耦合速率組成的系統,能夠承載一整個故事。單擺與繞行的行星是下一批重要角色,而它們共用同一個劇本:牛頓第二定律——力等於質量乘以加速度。由於加速度是二階導數,它永遠是一條二階方程。一個長度為 L、擺過角度 theta 的單擺,遵守 theta'' = -(g/L) sin(theta);一顆位置為 r、被太陽吸引的行星,遵守 r'' = -(G M / |r|^3) r。字母不同,文法相同——而且用的正是你一路打造起來的那套工具箱。
請留意單擺方程中關鍵的誠實之處:右邊是 sin(theta),而不是 theta。你先前熟悉的「簡諧」單擺——theta'' = -(g/L) theta,配著俐落的正弦與餘弦——只是小角度近似,僅在 theta 小到 sin(theta) 接近 theta 時才成立。真正的單擺是非線性的,而那一個 sin 使它無法以初等函數求解。我們正好身處「定性方法」那一階梯為我們準備好的地帶:多數有趣的 ODE 沒有閉式解,因此我們改從幾何上把行為讀出來。
把二階定律化為相圖
第一步是你已經做過許多次的動作:把二階方程化為一階系統,做法是替速度命名。令 v = theta',於是單擺變成 theta' = v 與 v' = -sin(theta)。現在我們有了一個二維自治系統——正是你學過如何描繪相圖的那一類:水平軸放角度,鉛直軸放角速度。單擺的每一個狀態,都是這個平面上的一個點,而它的整個未來,就是從那一點穿引而出的曲線。
平衡點在哪裡?要兩個速率同時為零,需要 v = 0 且 sin(theta) = 0,於是 theta = 0、pi、2pi 等等。其中兩個在物理上截然不同:theta = 0 是單擺筆直下垂,theta = pi 是它完美地倒立平衡。在底部附近做線性化,得到 theta'' = -theta,是純粹的振盪,因此那一點是個中心——四周環繞著閉軌,也就是來回的擺動。在頂端附近線性化,則得到 theta'' = +theta,有一個增長方向與一個衰減方向,正是鞍點的標誌。
能量:組織起一切的守恆量
這裡最深刻的念頭是:單擺擁有一個守恆量——它在擺動時總能量永不改變。對 theta'' = -sin(theta) 而言,能量為 E = (1/2) v^2 + (1 - cos(theta)):來自速度的動能項,加上來自高度的位能項。你可以直接驗證:dE/dt = v v' + sin(theta) theta' = v(-sin(theta)) + sin(theta) v = 0。由於 E 保持固定,每一條軌跡都必須落在單一條等位曲線 E = 常數上,而那些等位曲線*就是*相圖。守恆量不必求解方程,便把每一條軌道免費交到我們手裡。
至此,整幅圖像驟然對焦。低能量意味著繞著底部的小閉環——輕柔的擺動。隨著 E 增大,環逐漸鼓脹,直到恰好在 E = 2 時,曲線不再閉合,反而撞向倒立的鞍點:這條特殊的軌跡就是分界線(separatrix),是一條分隔界線,在此單擺*剛好*有足夠的能量緩緩爬向倒立、卻又永遠到不了。在 E = 2 之上,曲線完全不再回頭——單擺能量太大,便一次又一次翻過頂端,一圈接著一圈。三種性質迥異的命運,全都直接從能量等高線上讀出。
energy level what the pendulum does phase-plane shape ---------------- ----------------------------- ------------------- E small (E < 2) small back-and-forth swings closed loops (center) E = 2 exactly creeps toward upright forever separatrix (to saddle) E large (E > 2) whirls over the top, round wavy open curves
哈密頓觀點,及它為何能推廣
這種能量守恆的結構,並非單擺碰巧的好運;它是一整類系統。當方程能寫成 theta' = dH/dv 與 v' = -dH/dtheta(對某個單一函數 H 而言——這裡 H = E,即能量),我們就稱它為哈密頓系統。那兩個偏導數巧妙的反對稱配對,正是逼出 dH/dt = 0 的緣由,因此哈密頓系統*永遠*守恆 H。這也是為什麼它的平衡點是中心與鞍點,而絕不會是螺旋點或結點:沒有摩擦項讓能量洩漏,軌跡便永遠不會盤旋向內收進一點。
若加上真實世界的摩擦——空氣阻力、僵硬的樞軸——你就附上一個阻尼項:theta'' + b theta' + sin(theta) = 0。如今 dE/dt = -b v^2,只要單擺在運動就為負,於是能量穩定地流失。守恆量的論證失效,底部的中心變成一個吸引的螺旋點,每一次擺動都衰減向靜止下垂。這正是理想化、無摩擦模型,與你真正會打造的有阻尼模型之間誠實的分野——也正是為何一座真正的落地擺鐘,需要一根發條不斷把能量餵送回去。
克卜勒軌道,源自同一套工具箱
現在把同一套機器抬升到天上。克卜勒二體問題——一顆行星繞著重得多的太陽——由 r'' = -(G M / |r|^3) r 支配,又是牛頓定律化成的二階 ODE。它看起來比單擺糟糕得多,但同樣的兩個念頭就能馴服它。首先,這裡有兩個守恆量,而非一個:總能量 E,以及對太陽的角動量。每一條守恆律都是運動必須永遠遵守的關係,兩者合起來,幾乎把軌道完全釘死。
角動量守恆正是克卜勒第二定律——行星在相等時間內掃過相等面積,靠近太陽時疾馳、遠離時則拖沓。利用它消去角度,可把運動化簡為半徑在一個*有效位能*中的單一方程,這位能由真正的重力吸引,加上守恆律變出的一個「離心」項所組成。那條徑向方程又一次由它的能量所組織:束縛軌道安坐於位能井中、成為閉軌——橢圓——而恰好在逃逸能量、剛好未被束縛的軌道則描出拋物線,正是單擺分界線在天上的表親。
當單擺彼此耦合:簡正模態
最後一座橋,通回你先前精通的線性世界。掛起兩個單擺、以一根弱彈簧相連,並讓擺幅保持微小,使每個 sin(theta) 都退化為 theta,你便得到一對耦合振子——一個由兩條二階方程組成的線性系統。把它對角化,正如你在特徵值階梯裡對角化系統那樣,會揭示一些特殊的運動樣式,稱為簡正模態:在這些組合中,兩個單擺以單一共享的頻率振盪。任何運動,無論看起來多麼糾纏,都是這些模態的疊加。