線性化止步之處
前兩篇交給你一把銳利的工具,外加一個誠實的警告。建好雅可比矩陣、讀出它的特徵值,在雙曲平衡點處,哈特曼—葛羅曼定理保證線性判決確實就是局部真相。但只要有一個特徵值落到虛軸上——也就是臨界的、非雙曲的情形——這份保證便煙消雲散;此時被丟棄的彎曲項主宰了一切,而線性化只能聳聳肩。
無阻尼的單擺正是這類典型的頭痛問題。它底部靜止點的雅可比特徵值是 +/- i——一個純粹的中心,實部為零——所以線性化預測出無止盡的圓圈,卻說不出真正的非線性軌道究竟是緩慢向內盤旋、緩慢向外盤旋,還是恰好閉合。我們需要一種完全不倚賴線性化的方法。這方法就是李雅普諾夫直接法,而它的核心想法物理得近乎令人臉紅:盯著能量看。
你早已信賴的那份能量
想像一顆彈珠落在碗裡。無論你把它放在何處,總機械能——高度加上運動——只能維持不變,或因摩擦而一點一滴流失;它絕不可能自發地往上爬。因此彈珠無法一路向上越出碗外,而只要有任何摩擦,它終將停在碗底。請注意,我們得出這一切,從未求解過彈珠的路徑。我們純粹只是推論一個純量——能量——如何隨時間變化。
李雅普諾夫的躍進,在於看出彈珠的能量並非神聖不可替代——真正要緊的只是它的形狀與它的趨勢。任何在平衡點周圍呈碗狀、且當系統流動時絕不增大的純量函數 V(x, y),都會以完全相同的方式把軌跡困住。這樣的 V 稱為李雅普諾夫函數,它根本不必是真正的物理能量;只要具備那兩個正確的性質,它可以是我們自己發明的任何巧妙表達式。
兩個條件,與連鎖律的妙招
把「碗狀」與「絕不增大」說精確。其一,正定性:V 應在平衡點處恰為零,而在附近處處嚴格為正——這就是那只碗,其唯一的最低點恰好落在靜止狀態上。其二,是 V 沿著軌跡的變化率所滿足的趨勢條件,這個量寫作 V' 或 dV/dt。我們希望它非正,使得隨著時間推移,V 只能停在原處或往下滑。
接下來這一步,正是讓整套方法無須求解任何東西就能運作的關鍵。依多變數連鎖律,V 沿軌跡的變化率為 V' = (dV/dx) x' + (dV/dy) y'。但系統本身已告訴我們 x' = f(x, y)、y' = g(x, y)——於是我們直接把向量場代進去。未知的 x(t) 從頭到尾都不出現。我們純粹從 V 的公式與場的公式算出 V',判斷它的正負號,便能讀出系統的命運。
conditions on V(x,y) near an equilibrium (x*, y*):
V(x*, y*) = 0 and V > 0 elsewhere nearby (positive definite: a bowl)
V' = (dV/dx) x' + (dV/dy) y' (chain rule)
= (dV/dx) f(x,y) + (dV/dy) g(x,y) (substitute the field, no x(t) needed)
V' <= 0 ==> stable (energy never grows)
V' < 0 ==> asymptotically stable (energy strictly drains away)
讀出判決
現在來看那兩種情形。若在整個鄰域內 V' <= 0,軌跡便永遠爬不上更高的碗位,於是被困在平衡點附近:這就是李雅普諾夫穩定——起初靠近,便始終靠近。若趨勢更強一些,即除了平衡點本身之外處處 V' < 0,那麼這個值會持續嚴格流失,軌跡必定一路盤旋而下,抵達靜止點:這就是漸近穩定——起初靠近,並且真的抵達。
證明相反結論也有一個鏡像版本。若你能找到一個正定、但其 V' 嚴格為正的函數——也就是被迫增長的能量——那麼軌跡會被驅離,平衡點便不穩定。這有時稱為切塔耶夫式(Chetaev)論證,它讓同一套記帳法既能乾淨俐落地證明穩定,也能同樣俐落地證明不穩定。一個純量函數,便足以把問題往任一方向了結。
當 V' 只是觸碰零:拉薩爾
阻尼單擺裡藏著一個微妙之處。有摩擦時,只要擺錘在動,V' 就為負,但每當速度 y 經過零的瞬間——也就是每一次擺動的最高點——它會暫時等於零。因此 V' 並非嚴格為負;它只是在途中觸碰了零。單純的漸近穩定判準要求除平衡點外處處 V' < 0,而這個 V 在字面上沒過關,儘管單擺顯然確實會慢慢滑向靜止。
拉薩爾不變性原理(LaSalle)正好補上這道縫隙。它說:設 V' <= 0,那就去看 V' = 0 的那個集合,並追問有哪些完整的軌跡能永遠停留在這個集合裡。系統最終必定落在這種不變集裡最大的那一塊上。對阻尼單擺而言,V' = 0 要求始終 y = 0;但若擺錘永遠保持 y = 0,它便不可能在擺動,這就把它逼向底部靜止點。被困在零集裡的唯一軌跡,就是平衡點本身——於是拉薩爾把僅僅的穩定,重新提升為真正的漸近穩定。
誠實的難處:如何找到 V
要坦白這份能力的代價。這方法並未給出尋找李雅普諾夫函數的食譜——這正是它唯一真正的弱點。你必須猜一個候選的 V,再檢驗那兩個條件;若 V' 算出來符號不對,你的猜測什麼也沒告訴你,你就只能再換一個。對力學系統與電路系統而言,物理能量是自然的第一猜測;而在一般平衡點附近,像 V = a x^2 + b y^2 這樣的二次型則是主力,藉由調整常數 a 與 b 使 V' 算出來為非正。
還有兩個誠實的限度。這方法是單向的:找不到一個有效的 V,並不能證明不穩定——也許只是你還不夠聰明。而且它的結論通常是局部的,只在你的條件確實成立的那個鄰域上有效;要主張整體穩定,就得讓那只碗與向下的趨勢觸及處處,而這是個強得多、得仔細驗證的條件。然而在這些界限之內,回報非凡:對穩定性的完整判決——連那些雅可比矩陣啞口無言的臨界情形也涵蓋在內——僅憑一個精選的純量函數與一次連鎖律的運用便提取了出來,而那個解不出來的系統始終未被求解。