跳轉到

7.2 極限環與 Poincaré–Bendixson 定理#

本章導言#

線性系統的穩態振盪只有一種來源：無阻尼的中心（6.1 節），且它脆弱不堪——任何阻尼都毀掉它。然而真實世界充滿自我維持的振盪：心臟跳動、神經元放電、石英錶擺動、化學鐘。它們不靠外力、不怕擾動——擾動後自動回到同樣的振幅。這種軌道就是極限環（limit cycle）：一條被周圍軌道吸引（或排斥）的孤立閉軌。證明它存在的第一把理論鑰匙，是 Poincaré 提出、Ivar Bendixson（本迪克森，1861–1935） 於 1901 年嚴格化的Poincaré–Bendixson 定理——平面動力系統的鎮殿之寶。

極限環：孤立且穩定的閉軌#

定義：閉軌 \(\Gamma\) 是極限環，若附近存在非閉軌道螺旋趨向（穩定）或離開（不穩定） \(\Gamma\) ——即 \(\Gamma\) 在相空間中孤立（不像中心情形那樣被連續閉軌族包圍）。

類型 內側行為 外側行為 例子
穩定極限環 螺旋外擴 → 環 螺旋內縮 → 環 心跳、Van der Pol
不穩定極限環 內縮逃離 外擴逃離 分岔臨界的中間物
半穩定 內外同向 — 退化（結構不穩定）

與中心的區別是本節的概念核心：中心的閉軌不成孤（每條都緊鄰另一條），振幅由初值決定；極限環成孤，振幅由系統參數決定——自持振盪的「音量旋鈕」內建於方程，這就是為什麼心臟每次跳得一樣大。

Poincaré–Bendixson 定理#

定理（Poincaré 1881 思想，Bendixson 1901 證明）：平面自治系統 \(\dot{x} = f(x, y)\) 、 \(\dot{y} = g(x, y)\) （ \(f, g\) 光滑）的一條軌道若停留在緊集（有限區域）內且不趨於平衡點，則其 \(\omega\) -極限集必為閉軌（或連接平衡點的奇異圖）。

使用配方（Bendixson 袋）：找一個環形區域 \(R\) （內外兩條閉曲線圍成），使向量場在內邊界一律流出、外邊界一律流入（或反之）——軌道進入 \(R\) 就出不去；若 \(R\) 內沒有平衡點，則極限集只能是閉軌——極限環存在！

配角的理論根據：Jordan 曲線定理（平面上閉曲線分內外）——這是定理本質上是平面定理的原因。三維以上 Poincaré–Bendixson 失效（軌道可以打結、可以混沌，7.3 節）——平面的拓撲，是秩序的憲法。

否定面：Bendixson 判準#

判準：單連通區域 \(D\) 上若發散度 \(f_x + g_y\) 不變號，則 \(D\) 內無閉軌。

證明：若閉軌 \(\Gamma\) 存在，Green 定理（微積分！）給 \(\oint_\Gamma (f\,dy - g\,dx) = \iint_{\mathrm{int}} (f_x + g_y)\, dA\) 。左端 \(= \int (f \cdot \dot y - g \cdot \dot x)\, dt = \int (f g - g f)\, dt = 0\) （沿軌道 \(\dot x = f, \dot y = g\) ——向量場與切向量平行，法向通量零）；右端若被積函數不變號且非恆零則 \(\neq 0\) ——矛盾。∎

用法：想證極限環在這裡，用 Poincaré–Bendixson 袋；想證極限環不在那裡，用 Bendixson 判準——一對夾鉗，夾出極限環的位置。Dulac 判準（乘權函數 \(B(x, y)\) 後查 \(Bf_x + Bg_y + f B_x + g B_y\) ）是細化版。

範式：Van der Pol 振盪器#

Balthasar van der Pol（范德波爾，1889–1959） 1920 年代研究真空管電路（Philips 實驗室），方程化簡為

\[ \ddot{x} - \mu(1 - x^2)\,\dot{x} + x = 0. \]

等效負阻尼機制：寫成系統 \(\dot{x} = y\) 、 \(\dot{y} = -x + \mu(1 - x^2) y\) ，考察「能量」 \(E = \frac{x^2 + y^2}{2}\) ：

\[ \dot{E} = x\dot{x} + y\dot{y} = y(-x + y) + y\big(-x + \mu(1 - x^2) y\big)\cdot\frac{?}{} \quad\text{（用 } \dot y = -x + \mu(1 - x^2)y\text{）} = \mu(1 - x^2)\, y^2. \]

能量收支表： \(|x| < 1\) 時 \(\dot E > 0\) （注入）， \(|x| > 1\) 時 \(\dot E < 0\) （耗散）——小振幅被推、大振幅被拉，兩股力量匯合於一條平衡振幅的閉軌。這正是 Bendixson 袋的物理形狀：內圓外向、外環內向。

結構： \(\mu\) 小：近正弦振盪（諧波少，1927 年 van der Pol 與 van der Mark 觀察到的「鬆弛振盪」分頻是電子學史的著名實驗）； \(\mu\) 大：鬆弛振盪——慢積累（沿慢曲線爬）+ 快跳（沿快軌衝），波形像鋸齒/心電圖。心臟的竇房節放電正是鬆弛振盪（生物數學的標準對應）。

一個具體例題#

例題 7.2.1：證明系統

\[ \dot{x} = -y + x(1 - x^2 - y^2), \qquad \dot{y} = x + y(1 - x^2 - y^2) \]

存在穩定極限環，並求出它。

解：極座標是為此系統量身定做的： \(x = r\cos\theta\) 、 \(y = r\sin\theta\) 。標準變換（對 \(\dot x, \dot y\) 作 \(r\dot r = x\dot x + y \dot y\) 、 \(r^2 \dot\theta = x\dot y - y\dot x\) ）：

\[ r\dot{r} = x\dot{x} + y\dot{y} = x\big[-y + x(1 - r^2)\big] + y\big[x + y(1 - r^2)\big] = (x^2 + y^2)(1 - r^2) = r^2 (1 - r^2), \]
\[ r^2\dot{\theta} = x\dot{y} - y\dot{x} = x\big[x + y(1 - r^2)\big] - y\big[-y + x(1 - r^2)\big] = x^2 + y^2 = r^2. \]

故（ \(r > 0\) 時）：

\[ \dot{r} = r(1 - r^2), \qquad \dot{\theta} = 1. \]

讀解： \(\dot\theta = 1\) ——以恆定角速度旋轉； \(\dot r = r(1 - r^2)\) ——一維邏輯斯蒂型（1.3 節練習 5 的親戚）： \(r < 1\) 時 \(\dot r > 0\) （半徑增長）， \(r > 1\) 時 \(\dot r < 0\) （半徑收縮）， \(r = 1\) 是穩定平衡。故所有軌道螺旋趨向

\[ \boxed{\,r = 1\,} \qquad \text{即圓 } x^2 + y^2 = 1, \]

上以角速度 1 巡行——穩定極限環，振幅恆為 1，與初值無關（只要 \(r_0 > 0\) ）。 \(r = 0\) （原點）是不穩定焦點（線性化 \(\dot x = -y + x, \dot y = x + y\) ，特徵值 \(1 \pm i\) ，實部正）。

Bendixson 袋覆核：內邊界 \(r = 0.5\) 上 \(\dot r = 0.5 \times 0.75 > 0\) （向外）；外邊界 \(r = 2\) 上 \(\dot r = 2 \times (-3) < 0\) （向內）；環形域 \(0.5 < r < 2\) 內無平衡點——Poincaré–Bendixson 適用，閉軌必存在；本例更進一步顯式找到了它。∎

方法的教訓：極座標把二維系統「解耦」為一維徑向 + 一維角向——對稱性良好的系統可以用座標選擇直接馴服（1.3 節等角螺線的同族魔法）。一般系統沒有這份運氣，只能靠袋與判準的夾鉗。

歷史啟示#

  • 從「解方程」到「證存在」：Poincaré–Bendixson 不給公式，只給「某處有一條閉軌」的存在證明——定性理論的成年禮：數學家第一次理直氣壯地「知道而不算出」。
  • 生物的鐘錶：心跳（SA 節）、晝夜節律（Goodwin 1965、Kai 蛋白振盪 1990s）、化學鐘（Belousov–Zhabotinsky 反應，1951/1964——Belousov 的論文被拒五年，Zhabotinsky 復活它，最終 Zhabotinsky 獲 1980 列寧獎而 Belousov 已於 1970 鬱鬱而終）——生命就是一串極限環。
  • 平面是秩序的孤島：P–B 定理保證平面有限行為（閉軌或平衡）；三維一開，混沌進場（7.3）。維數 3 是秩序與混沌的國界——本定理是國界線上的海關。

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

  • 6.2 節：指標定理（閉軌內指標和 \(+1\) ）與 Bendixson 判準——排除閉軌的雙工具。
  • 6.3 節：LaSalle 原理與 Lyapunov 技巧——Bendixson 袋的構造常以 Lyapunov 函數等位線為牆。
  • 7.1 節：受迫共振的穩態振盪是「被動振盪」；極限環是「主動振盪」——兩者經 7.1 節能量收支視角統一。
  • 7.3 節：Lorenz 系統中極限環讓位於奇怪吸引子——本定理失效的現場。
  • 神經科學：FitzHugh–Nagumo（1961）是 Van der Pol 的生物版；動作電位 = 快慢系統的極限環；神經元「激發閾值」= 不穩定極限環的存在。
  • 工程：真空管振盪器、風奏音（von Kármán 渦列）、閥門水擊顫動——一切「自發嗡鳴」都是極限環；除振設計 = 破壞袋條件。
  • 分岔理論：Hopf 分岔（6.2 場景一）是極限環的出生方式（超臨界：穩定環；亞臨界：不穩定環——後者對應突發災難式失穩，如某些風振）。

練習#

  1. 對例題系統，寫出 \(r(t)\) 的顯式解（分離變數： \(\dot r = r(1 - r^2)\) ），並求 \(r_0 = 2\) 時 \(r(t) \to 1\) 的時間尺度。
  2. 用 Bendixson 判準證明： \(\dot{x} = y\) 、 \(\dot{y} = -x - y + x^2\) 在全平面（單連通）無閉軌。
  3. 對 Van der Pol（ \(\mu = 1\) ）：以 Liénard 變換（ \(u = x - x^3/3 - y/\mu\) 型）討論其「慢曲線」與鬆弛振盪的慢積累段（定性描述即可）。
  4. 設計袋：對 \(\dot{x} = -y + x(4 - x^2 - y^2)(x^2 + y^2 - 1)\) ，找出兩個極限環的位置並判斷各自的穩定性（提示：極座標下 \(\dot r = r(4 - r^2)(r^2 - 1)\) ）。
  5. 思考題：心臟的竇房節細胞若被「去極化封鎖」（參數改變使極限環消失、系統塌向穩定平衡），心臟停跳。請從 Hopf 分岔角度解釋：這是超臨界還是亞臨界分岔取決於什麼？臨床意義（可逆/突發）有何不同？

進階研讀：極限環的五個深化場景#

場景一：Liénard 方程與泛周期性。 \(\ddot{x} + f(x)\dot{x} + g(x) = 0\) （Liénard 1928）：條件（ \(F\) 奇函數形狀單峰、 \(g\) 奇且 \(x g > 0\) 等，Liénard–Levinson–Smith 定理 1943）保證唯一穩定極限環。Van der Pol 是 \(f = \mu(x^2 - 1)\) 的特例。「阻尼函數的形狀」決定極限環的存在與唯一——一套可驗證的工程判據。

場景二：化學振盪與其十九年的幽禁。Belousov 1951 年發現檸檬酸—溴酸鹽—鈰反應自行振盪（顏色週期變化），論文被拒（「違反熱力學第二定律」的審稿誤解——振盪不違反平衡態熱力學，只需遠離平衡的開放/准開放條件）；Zhabotinsky 1964 重啟研究。BZ 反應的 Oregonator 模型（Field–Nörös 1974）是標準極限環系統。科學共同體拒絕了它十九年——振盪理論的社會學教訓。

場景三：相鎖與混頻。受迫極限環（外部週期力 + 自持振盪）展現Arnold 舌（Arnold 1961–：鎖頻區 \(\frac{\omega}{\omega_{\rm ext}} \approx \frac{p}{q}\) 有理數附近）、準週期（兩頻不可通約）與通向混沌的路線。心臟起搏器的干擾、圓周映射 \(\theta \to \theta + \Omega - \epsilon \sin\theta\) 是其最小模型——極限環 + 強迫 = 從有序到混沌的旋鈕。

場景四：空間中的波——極限環的擴散版。反應—擴散方程（Turing 1952）中，極限環動力學 + 擴散產生行進波與螺旋波：BZ 反應的螺旋圖案、心臟纖顫（螺旋波失調——電擊除顫 = 重置相位，臨床理論 Winfree 1987）。一個點的極限環，鋪滿空間就是生命與疾病。

場景五：離散版的極限環。二維映射（如 Poincaré 截面）的不變閉曲線對應準週期運動；映射的不動點對應週期解。KAM 定理（Kolmogorov 1954、Arnold–Moser 1960s）證明：微擾哈密頓系統中大部分不變環面存活——「極限環」思想的哈密頓表親（雖無吸引性，因體積守恆）。太陽系長期穩定的數值研究（Laskar 1989）建立於 KAM 語言。

補充：Bendixson 其人#

Ivar Bendixson（1861–1935），瑞典斯德哥爾摩數學家（後任 Stockholm Högskola 校長、皇家科學院要職），1901 年論文嚴格化並推廣 Poincaré 的平面軌道理論（含極限集概念 \(\omega\) -limit set 的正式定義）。同代瑞典數學（Mittag-Leffler 的 Acta Mathematica 為歐洲三大期刊之一）使北歐成為定性理論的橋頭堡。一個「補完定理」的工作者，名字永遠與天才並列在定理名中——數學的署名制度對「整理者」相當慷慨。

補充：三維為何「失控」——亂流的第一次數學預言#

平面系統軌道不能交叉（3.2 節）+ Jordan 定理 ⇒ 緊軌道只有閉曲線的命運。三維中，軌道仍不交叉，但可以纏繞：不打結的線在一點的極限集可以是無窮複雜的集合（Smale 馬蹄 1967 的映射版、Lorenz 吸引子 1963 的流版）。「不交叉」在二維是緊箍咒，在三維是自由的開始——7.3 節就從這道門進入。

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

  • 極限環：孤立的閉軌；自持振盪的數學形式。
  • 穩定/不穩定/半穩定極限環：內外側軌道的趨近方向。
  • Poincaré–Bendixson 定理：平面緊軌道的極限集 = 閉軌或奇異圖。
  • Bendixson 袋：向量場進出不平衡的環域；極限環的存在構造。
  • Bendixson 判準：散度不變號 ⇒ 無閉軌。
  • Dulac 判準：加權散度版判準。
  • ** \(\omega\) -極限集**：軌道反覆訪問的點集。
  • Van der Pol 振盪器： \(\ddot x - \mu(1 - x^2)\dot x + x = 0\) ；負阻尼自限。
  • 鬆弛振盪：大 \(\mu\) 下的慢積累 + 快跳波形。
  • Liénard 方程： \(\ddot x + f(x)\dot x + g(x) = 0\) ；極限環判據框架。
  • Arnold 舌：受迫振盪的鎖頻區域。

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

  • 1881–1886：Poincaré 積分曲線理論（極限環概念）。
  • 1901：Bendixson 定理嚴格化。
  • 1920s：Van der Pol 真空管振盪（1926、1927 鬆弛振盪與分頻實驗）。
  • 1928：Liénard 定理。
  • 1951：Belousov 發現化學振盪（論文被拒）。
  • 1952：Turing 形態發生（反應—擴散）。
  • 1961：FitzHugh–Nagumo；Arnold 舌。
  • 1964：Zhabotinsky 復活化學振盪研究。
  • 1974：Oregonator 模型。
  • 1987：Winfree《心臟節律的幾何》（除顫理論）。

深化問答三則#

問：極限環與 7.1 節的「穩態受迫振盪」有什麼本質區別？ 答：受迫振盪的頻率與振幅由外力決定（方程右端），去掉外力即消失；極限環的頻率與振幅由系統自身決定（方程左端的結構），是方程的解。物理學家說：前者是響應，後者是節律——心臟不問外界頻率而自跳。數學上：前者是非自治系統的特解，後者是自治系統的閉軌。自治性（無外力）是生命節律的數學定義。

問：Poincaré–Bendixson 為何在本書同時是「天使」與「天花板」？ 答：天使：它給平面系統完整的漸近行為分類（平衡、閉軌、圖——沒有別的可能），是定性理論最強的完整性結果。天花板：它同時宣告平面不可能有混沌（極限集窮盡了），故尋找混沌必須上三維。一個定理既是答案也是邊界——數學理論的成熟往往如此。

問：如何在工程上「製造」或「消滅」極限環？ 答：製造：設計負阻尼區（真空管、雷射泵浦、Van der Pol 電路），並保證大振幅的正阻尼回壓（能量收支交點）。消滅：破壞 Bag 條件（加線性阻尼使外邊界也向外）、或使 Bendixson 判準成立（加發散度恆正的項）、或把 Hopf 分岔參數拉回穩定側。除顫器是消滅極限環的醫療儀器，石英振盪器是製造極限環的計時儀器——同一數學，兩個文明。

（本節完）