JOVANA
Explore Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

不動點、穩定性與相圖

非線性系統的定性工具箱:找出流動靜止之處、判斷它是否穩定,並在不精確求解任何一條方程式的情況下,讀懂整張相圖。

不必求解也能讀懂流動

大多數非線性微分方程沒有封閉解——你寫不出一條整潔的 x(t) 公式。龐加萊的革命,是不再苦求公式,改問幾何問題:流動在哪裡靜止?在附近往哪個方向走?長期會如何?這些答案構成一張 相圖,而它們無需精確解。

任何相圖的骨架,是它的不動點連同其穩定性,再加上連接它們的特殊軌跡。掌握這些,你便定性地理解了系統的行為——而這往往正是物理學家所需的全部。

不動點與一維流動

不動點 \mathbf{x}^* 是流動消失之處,\mathbf{f}(\mathbf{x}^*)=0:把系統恰好放在那裡,它便永不移動。這些就是平衡態。先從最簡單的情形入手——直線上的單一變數。

\dot{x} = f(x), \qquad f(x^*) = 0

在一維中,不動點就是 f 的零點。兩零點之間,f 的正負決定狀態滑向何方——f > 0 向右,f < 0 向左。

穩定性由一次擾動決定。把狀態擾到 x = x^* + \eta 並線性化:因為 f(x^*)=0,微小位移滿足 \dot\eta = f'(x^*)\,\eta。若斜率 f'(x^*) 為負,擾動衰減,該點穩定;為正則成長,該點不穩定。這就是線性穩定性分析最純粹的形式。

\dot{\eta} = f'(x^*)\,\eta \quad\Longrightarrow\quad \eta(t) = \eta_0\, e^{\,f'(x^*)\,t}

擾動以速率 f'(x*) 指數成長或衰減。斜率為負 ⇒ 穩定平衡;斜率為正 ⇒ 不穩定。

高維的線性穩定性分析

同樣的想法適用於任意維度——只要把導數換成在不動點取值的雅可比矩陣(偏導數構成的矩陣)。它的特徵值 \lambda_k 決定每個擾動方向如何成長或收縮。

\dot{\boldsymbol{\eta}} = J\,\boldsymbol{\eta}, \qquad J_{ij} = \left.\frac{\partial f_i}{\partial x_j}\right|_{\mathbf{x}^*}, \qquad \boldsymbol{\eta}(t) = \sum_k c_k\, \mathbf{v}_k\, e^{\lambda_k t}

在不動點附近,流動由雅可比矩陣 J 主宰。當且僅當每個特徵值的實部皆為負時,不動點才穩定。

二維不動點的分類

在二維中,你很少需要親手算特徵值:兩個數字就夠了—— \tau = \operatorname{tr}J行列式 \Delta = \det J,因為 \tau = \lambda_+ + \lambda_-\Delta = \lambda_+\lambda_-

\lambda_\pm = \tfrac{1}{2}\left(\tau \pm \sqrt{\tau^2 - 4\Delta}\,\right), \qquad \tau = \operatorname{tr}J,\ \ \Delta = \det J

跡–行列式平面能一眼分類每個二維不動點。

規則如下:若 \Delta < 0,兩特徵值為異號實數——鞍點(永遠不穩定)。若 \Delta > 0,該點為節點或焦點:\tau < 0穩定\tau > 0不穩定\tau^2 - 4\Delta 的正負決定是節點(實根,>0)還是焦點(複根,<0)。邊界 \tau = 0,\ \Delta > 0中心——封閉環,如無阻尼振子。

範例:阻尼單擺

讓我們把整套方法用在一個真實系統上:一個阻尼單擺,由 \ddot\theta + \gamma\dot\theta + \omega_0^2\sin\theta = 0 支配。保留完整的 \sin\theta(不用小角度取巧)使它成為非線性,因此我們用幾何方法分析。

  1. 寫成一階系統。v = \dot\theta。二階方程化為一對式子 \dot\theta = v\dot v = -\gamma v - \omega_0^2\sin\theta——即 (\theta, v) 相空間中的流動。
  2. 找不動點。\dot\theta = \dot v = 0:需要 v = 0\sin\theta = 0,得 \theta = 0(垂直下垂)與 \theta = \pi(垂直倒立平衡)。
  3. 建立雅可比矩陣。 求導得 J = \begin{pmatrix} 0 & 1 \\ -\omega_0^2\cos\theta & -\gamma \end{pmatrix},於是跡 \tau = -\gamma、行列式 \Delta = \omega_0^2\cos\theta
  4. 分類底部(\theta = 0)。 此處 \Delta = \omega_0^2 > 0\tau = -\gamma < 0穩定點——若 \gamma < 2\omega_0 為焦點(阻尼振盪),過阻尼則為節點。單擺最終停在這裡,與預期完全相符。
  5. 分類頂部(\theta = \pi)。 此處 \cos\pi = -1,故 \Delta = -\omega_0^2 < 0:永遠是鞍點,因而不穩定。最輕微的一推便使它離開——正符合「要讓單擺筆直倒立平衡是無望的」這一日常真理。

極限環,以及為何平面禁止混沌

有些系統最終停在的不是一個點,而是一個孤立的封閉環——極限環,一種擁有自身偏好振幅的自持振盪(心跳、雷射、范德波耳電路)。不同於中心點周圍那一圈中性軌道,極限環是孤立的,並吸引鄰近軌跡,因此它的振幅對擾動具有韌性。

接著是一個深刻的結構性事實,龐加萊–本迪克森定理:在二維連續自治流動中,有界軌跡只能趨近一個不動點或一個極限環——不會有更奇異的東西。軌跡不能相交,而平面實在太過侷促,無法讓一條有界路徑永遠遊走而不重複。因此平面上不可能有混沌