跳轉到

6.2 相平面與平衡點分類#

本章導言#

線性系統的平衡點已被特徵值分類（6.1 節）；非線性系統的平衡點能否照辦？答案是幾乎可以——在平衡點處線性化，只要線性化後的特徵值「不踩在虛軸上」（雙曲情形），非線性平衡點與其線性像拓撲等價（Hartman–Grobman 定理，1959–1960）。本節建立相平面的完整分析工具：平衡點圖譜、線性化判據的邊界（臨界情況）、以及 Poincaré 指標——並以捕食者—獵物系統為例展示「不求解、只讀圖」的定性藝術。

平衡點與線性化#

系統 \(\dot{x} = f(x, y)\) 、 \(\dot{y} = g(x, y)\) 的平衡點（奇點/不動點） \((x^*, y^*)\) 滿足 \(f = g = 0\) ——向量場在此歸零，軌道停駐。

線性化：在平衡點展開 \(f \approx f_x (x - x^*) + f_y (y - y^*)\) ，得 Jacobian 矩陣

\[ J = \begin{pmatrix} f_x & f_y \\ g_x & g_y \end{pmatrix}_{(x^*, y^*)}, \]

其特徵值 \(\lambda_{1,2}\) 決定局部行為。判據表（記 \(\tau = \operatorname{tr} J\) 、 \(\Delta = \det J\) ）：

條件 類型 穩定性
\(\Delta > 0,\ \tau^2 - 4\Delta < 0\) 焦點（螺旋） \(\tau < 0\) 穩定， \(\tau > 0\) 不穩
\(\Delta > 0,\ \tau^2 - 4\Delta > 0\) 結點 同上
\(\Delta < 0\) 鞍點 不穩定（總是）
\(\Delta = 0\) 或 \(\tau^2 = 4\Delta\) 踩線 非雙曲 需要更多分析！

** \(\tau\) – \(\Delta\) 平面（trace–determinant diagram，Strogatz 的著名圖）把上述判據畫成一張圖：拋物線 \(\tau^2 = 4\Delta\) 是結點與焦點的分界，下半平面全是鞍點——兩個純量決定一張地圖**。

Hartman–Grobman：線性化何時可信#

定理（Hartman 1960、Grobman 1959）：若平衡點雙曲（ \(J\) 的特徵值實部皆非零），則在平衡點鄰域內，非線性流與線性化流拓撲共軛（存在保定向的連續變換把一者的軌道映成另一者的軌道）。

含義：雙曲平衡點的定性行為（結點/鞍/焦、穩定與否）由線性化完全決定——高階項無話語權。臨界情況（虛軸上有特徵值）例外：純虛根時線性化說「中心」，非線性可能真給中心（Lotka–Volterra），也可能給穩定或不穩定螺旋——必須看二階項。線性化是雙曲世界的通行證，出了雙曲世界就要自帶證據。

Poincaré 指標：平衡點的拓撲指紋#

定義：沿相平面上包圍平衡點的閉曲線走一圈，向量場方向角總變化除以 \(2\pi\) ，即指標（index）。

  • 結點、焦點、中心：指標 \(+1\) ；
  • 鞍點：指標 \(-1\) 。

Poincaré–Hopf 指標定理：閉軌（週期解）的指標為 \(+1\) ；含若干平衡點的區域邊界上，向量場總指標 = 內部平衡點指標之和。立即推論：閉軌內部平衡點指標之和為 \(+1\) ——包著兩個鞍點的閉軌不存在（ \(-1 - 1 \neq 1\) ）；包一個鞍點與一個結點可以（ \(-1 + 1 = 0 \neq 1\) ，也不行！必須和為 \(+1\) ，如一個鞍點加兩個結點，或單一結點/焦點/中心）。指標是閉軌的「出生許可證」——不解方程就排除大量可能性，這是 Poincaré 1881 年定性理論的招牌絕技。

案例研究：Lotka–Volterra 捕食者—獵物系統#

\[ \dot{x} = x(a - b y), \qquad \dot{y} = y(-c + d\, x) \qquad (a, b, c, d > 0). \]

（獵物 \(x\) ：無捕食者時指數增長 \(a\) ；捕食者 \(y\) ：無獵物時指數衰減 \(-c\) ；相遇率 \(bxy\) 損獵物、 \(dxy\) 餵捕食者。Volterra 1925–26 為解釋亞得里亞海漁獲量波動提出；美國生化學家 Alfred Lotka 1920 年獨立提出於化學振盪研究。）

平衡點： \((0, 0)\) 與 \(\big(\frac{c}{d}, \frac{a}{b}\big)\) 。

線性化：Jacobian \(J = \begin{pmatrix} a - by & -bx \\ dy & -c + dx \end{pmatrix}\) 。

  • 在原點： \(\begin{pmatrix} a & 0 \\ 0 & -c \end{pmatrix}\) ，特徵值 \(a > 0, -c < 0\) ——鞍點（雙曲，結論可靠）：滅絕與爆炸的分水嶺。
  • 在 \(\big(\frac{c}{d}, \frac{a}{b}\big)\) ： \(J = \begin{pmatrix} 0 & -\frac{bc}{d} \\ \frac{ad}{b} & 0 \end{pmatrix}\) ，特徵值 \(\pm i\sqrt{ac}\) ——純虛根，臨界情況！線性化只說「中心候選」，需要二階證據。

二階證據（不變量）：考慮

\[ V(x, y) = d\, x - c \ln x + b\, y - a \ln y. \]

沿軌道求導： \(\dot{V} = \big(d - \frac{c}{x}\big)\dot{x} + \big(b - \frac{a}{y}\big)\dot{y} = \frac{d x - c}{x}\, x(a - by) + \frac{b y - a}{y}\, y(-c + dx) = (dx - c)(a - by) - (a - by)(dx - c) = 0.\)

** \(V\) 恆定——軌道是 \(V\) 的等位線（閉合的卵形線，因 \(V\) 在 \(\big(\frac{c}{d}, \frac{a}{b}\big)\) 嚴格凸且趨 \(\infty\) 於邊界）。故內平衡點是真中心，軌道全是閉曲線：捕食者與獵物以恆定振幅永續振盪**。

歷史後記：Volterra 的「閉軌預言」在漁業管理中名聲起伏——真實生態系統有阻尼（年齡結構、空間異質），振盪常衰減或放大；May（1970s）證明更複雜的食物網傾向不穩定。教科書模型的閉軌是結構不穩定的（參數微擾即毀）——這個「結構穩定性」概念（Andronov–Pontryagin 1937）本身成為動力系統的核心思想。

一個具體例題#

例題 6.2.1：分析系統 \(\dot{x} = 1 - xy\) 、 \(\dot{y} = x - y^3\) 的平衡點與類型。

解：

求平衡點： \(xy = 1\) 且 \(x = y^3\) 。代入： \(y^4 = 1 \Rightarrow y = 1\) （實數情形， \(y = -1\) 給 \(xy = 1\) ？ \((-1)(-1) = 1\) ✓， \(x = -1\) ✓——兩個平衡點： \((1, 1)\) 與 \((-1, -1)\) ）。

Jacobian： \(J = \begin{pmatrix} -y & -x \\ 1 & -3y^2 \end{pmatrix}\) 。

在 \((1, 1)\) ： \(J = \begin{pmatrix} -1 & -1 \\ 1 & -3 \end{pmatrix}\) ， \(\tau = -4\) ， \(\Delta = 3 + 1 = 4\) ，判別式 \(\tau^2 - 4\Delta = 16 - 16 = 0\) ——臨界（重特徵根 \(\lambda = -2\) 兩次）。跡為負、行列式為正，必是某種吸引子；重根情形結點（退化結點/結狀結點）與焦點可能——直接檢驗：矩陣非純量矩陣且只有一個特徵方向，軌道沿該方向切向流入——退化穩定結點（attracting）。

在 \((-1, -1)\) ： \(J = \begin{pmatrix} 1 & 1 \\ 1 & -3 \end{pmatrix}\) ， \(\tau = -2\) ， \(\Delta = -3 - 1 = -4 < 0\) ——鞍點（ \(\Delta < 0\) 無條件）。

指標檢查：若有閉軌，其內部指標和須為 \(+1\) ：包住 \((1,1)\) （指標 \(+1\) ）可行；包住 \((-1,-1)\) （指標 \(-1\) ）必須再補兩個正指標平衡點——本系統沒有，故包圍鞍點的閉軌不可能。這個「不解方程的排除法」正是指標理論的實戰價值。

歷史啟示#

  • 定性先於定量：Lotka–Volterra 的結論（振盪、分水嶺、中心）無一需要顯式解——Poincaré 範式的勝利：「畫出場，讀出行為」。
  • 臨界情況是科學的縫隙：純虛根處線性化啞火，恰好是生物振盪、工程失穩（Hopf 分岔，1942）的棲息地——分岔理論（誕生於這些「例外」）在二十世紀成為顯學。
  • 結構穩定性觀念：Volterra 閉軌一擾即毀；結點與鞍點則堅固。「什麼現象在擾動下存活」比「什麼現象存在」更接近自然——Andronov 學派的座右銘。

與其他章節及領域的關聯#

  • 6.1 節：線性分類字典（本節的理論基礎）。
  • 6.3 節：Lyapunov 函數——本節 Lotka–Volterra 的 \(V\) 是其非嚴格版（ \(\dot V = 0\) 而非 \(< 0\) ）；嚴格版給漸近穩定。
  • 7.2 節：閉軌與極限環——指標理論的下一步；本節「閉軌許可證」在那裡升級為「存在性定理」。
  • 7.3 節：三維、鞍點成對、混沌——指標理論的高維擴充（Smale 的指標配對）。
  • 10.2 節：Lotka–Volterra 與其後繼（邏輯斯蒂捕食、化學振盪 Brusselator、電路 Van der Pol）。
  • 神經科學：Hodgkin–Huxley（1952）神經元模型的相平面分析（二維簡化 FitzHugh–Nagumo 1961）——動作電位 = 極限環的誕生（Hopf 分岔）；一次神經衝動是一張相平面圖。
  • 經濟學：IS–LM 動態、Goodwin 成長循環（1967）直接借用 Lotka–Volterra 數學。

練習#

  1. 對系統 \(\dot{x} = 2x(1 - x) - xy\) 、 \(\dot{y} = -y + 2xy\) （邏輯斯蒂捕食模型）：求全部平衡點並分類。
  2. 用 \(\tau\) – \(\Delta\) 平面圖：畫出 \(\tau = -3\) 、 \(\Delta = 2\) 的平衡點類型；再討論沿直線 \(\Delta = \frac{\tau^2}{4}\) 移動時類型如何過渡。
  3. 對 \(\dot{x} = y\) 、 \(\dot{y} = \mu y + x - x^2\) （含參數 \(\mu\) ）：求平衡點，討論 \(\mu\) 變化時類型如何變（特徵值遷移），指出 \(\mu\) 的臨界值。
  4. 驗證 Lotka–Volterra 的 \(V\) 等位線在 \(\big(\frac{c}{d}, \frac{a}{b}\big)\) 處取得嚴格最小值（Hessian 正定），從而證明等位線閉合。
  5. 思考題：指標定理禁止「閉軌只包一個鞍點」。請設計一個向量場示意圖（畫箭頭），直觀解釋為什麼軌道繞鞍點一圈後方向角淨變化是 \(-2\pi\) ，而繞結點一圈是 \(+2\pi\) 。

進階研讀：定性方法的五個深化場景#

場景一：Hopf 分岔。參數 \(\mu\) 掃過臨界值時，平衡點由穩定焦點變不穩定焦點、誕生一個極限環（Hopf 1942；Poincaré 1892 已知）。例： \(\dot{r} = \mu r - r^3\) 、 \(\dot{\theta} = \omega\) —— \(\mu < 0\) 吸引到原點， \(\mu > 0\) 收縮到圓 \(r = \sqrt\mu\) 。Van der Pol 振盪器（7.2 節）、神經元放電、飛機翼振（flutter）都是 Hopf 分岔。臨界情況不只是麻煩，是現象的產房。

場景二：分離軌跡（separatrix）與分水嶺。鞍點的穩定/不穛定流形（invariant manifolds）是相平面的分水嶺：兩側軌道奔向不同命運。擺的相圖中，能量恰等於「最高點能量」的軌道是分離軌跡——它上方是旋轉、下方是振盪。Poincaré 1899 年發現分離軌跡橫截相交（同宿纏結）時行為極端複雜——混沌的第一次現形（7.3 節的史前史）。

場景三：Bendixson 判準（否定閉軌）。若單連通區域上 \(\nabla \cdot \mathbf{f} = f_x + g_y\) 不變號，則區域內無閉軌（Green 定理：閉軌內流量 = 環積分 \(\neq 0\) 矛盾）。與指標定理配合，構成「排除閉軌」的雙工具。Dulac 判準（以函數 \(B\) 加權的推廣）是十九–二十世紀之交俄法學派的細化。

場景四：中心 vs 細焦點的判定。純虛根時，區分「真中心」與「微弱螺旋」（細焦點，focus of multiplicity \(k\) ）需要計算Lyapunov 量（Poincaré 首項法）：焦點級數展開中第一個非零的徑向項 \(\dot{r} = c_k r^{2k+1}\) 。Lotka–Volterra 的所有 Lyapunov 量為零（可積系統），故真中心；一般系統不零。「中心是奇蹟，螺旋是常態」——可積性是測度零的美好。

場景五：從相平面到高維。三維以上，平衡點分類照樣線性化（特徵值實部符號），但全局結構質變：Poincaré–Bendixson 失效、馬蹄映射（Smale 1967）出現、混沌合法化（7.3 節）。維數的分水嶺：2 維有序到無聊、3 維有序到瘋狂——這一「維數門檻」是二十世紀數學最深的驚喜之一。

補充：Vito Volterra 與漁獲數據#

Vito Volterra（1860–1940）是義大利參議員、羅馬大學校長、實變函數論大家（Volterra 積分方程）。1925 年，女婿 Umberto D'Ancona（海洋生物學家）給他看亞得里亞海漁獲統計：一戰期間（捕撈減少）捕食性魚類比例反而大增。Volterra 以他的方程解釋：捕撈壓力減小對捕食者（繁殖慢、受捕食制約）相對有利——「Volterra 原理」今日稱為競爭釋放，是保育生物學的基準推理。一組漁獲數字，催生了數理生態學。1931 年 Volterra 因拒絕向法西斯政權宣誓效忠被逐出大學——數理生態學之父以科學良知落幕。

補充：τ–Δ 平面的記憶法#

把 \((\tau, \Delta)\) 想成「體溫與體重」： \(\Delta\) 管「平衡與否」（ \(\Delta < 0\) 必鞍——「分裂人格」）， \(\tau\) 管「吸與呼」（ \(\tau < 0\) 吸、 \(> 0\) 呼）。拋物線 \(\tau^2 = 4\Delta\) 是結點—焦點的國界：越線時兩特徵根由相異實數（結點）在重合處「結婚」變複數（焦點）。兩個純量講完二維線性世界的故事——這種「降維到不變量」的品味貫穿數學（判別式、Euler 示性數、譜）。

名詞速查表（本節出現的核心術語）#

  • 平衡點： \(\mathbf{f} = 0\) ；流的停駐點。
  • 線性化 / Jacobian：平衡點處的一階近似矩陣。
  • 雙曲平衡點：特徵值實部皆非零；Hartman–Grobman 可用。
  • 結點 / 焦點 / 鞍點 / 中心：平衡點四類型。
  • ** \(\tau\) – \(\Delta\) 圖**：跡—行列式分類平面。
  • Poincaré 指標：向量場繞平衡點的轉角圈數；結點 \(+1\) 、鞍 \(-1\) 。
  • 分離軌跡：鞍點流形構成的命運分界。
  • Bendixson 判準：散度不變號 ⇒ 無閉軌。
  • Hopf 分岔：平衡失穩誕生極限環（1942）。
  • 結構穩定性：小擾動下定性行為不變（Andronov–Pontryagin 1937）。

年表速查（本節涉及的關鍵年份）#

  • 1881–1886：Poincaré 相平面理論與指標。
  • 1920：Lotka 化學振盪模型。
  • 1925–26：Volterra 漁獲分析。
  • 1937：Andronov–Pontryagin 結構穩定性。
  • 1942：Hopf 分岔定理。
  • 1952：Hodgkin–Huxley 神經方程。
  • 1959–1960：Grobman–Hartman 定理。
  • 1961：FitzHugh–Nagumo 相平面神經模型。
  • 1967：Goodwin 經濟週期模型。

深化問答三則#

問：為什麼「臨界情況」這麼難，卻這麼重要？ 答：難，因為高階項接管話語權，沒有通用定理；重要，因為系統演化與相變恰恰發生在臨界處——參數掃過臨界面時，結點變焦點（Hopf）、穩定變失穩、中心誕生或毀滅。工程師的失敗模式（顫振、共振、失穩）全都住在臨界情況裡。定理管常態，直覺與計算管臨界——這是定性方法的分工。

問：指標理論「不解方程」的威力，邊界在哪？ 答：指標只給必要性（閉軌內指標和必為 \(+1\) ），不給充分性——指標合格不代表閉軌存在。存在性需要 7.2 節的 Poincaré–Bendixson 定理或構造（Lyapunov 型函數）。指標是安檢站，不是產房——它能攔住不可能的，不能保證可能的。

問：Lotka–Volterra 中心「結構不穩定」意味着模型沒用嗎？ 答：意味着逐條結論要謹慎：閉軌振幅恆定這一結論會被任何小擾動（密度制約、時滯）破壞；但定性現象（共存平衡存在、擾動後振盪回歸或發散、捕撈的相對效應）在廣泛的模型族中存活。模型的價值在「魯棒結論」，不在「脆弱定量」——這是數學建模的第一課（10.2 節展開）。

（本節完）