05.1 可積性與混沌#
淺顯說明#
到目前為止,我們解過的系統(自由質點、諧振子、單擺、自由落體)都有一個共同點: 它們都能被「完全解出」——運動方程可以積分,軌跡可以寫成公式。 這樣的系統叫可積系統(integrable system)。
但大自然不是處處可積。Poincaré 在 1890 年研究三體問題(太陽、地球、月亮互相吸引)時 發現:三體問題不可積!不存在足夠多的守恆量,軌跡寫不成公式。
更可怕的是,他發現不可積系統的相空間裡有「同宿纏結」—— 軌跡像亂成一團的毛線球,初始條件的微小差異,會隨時間放大成完全不同的命運。 這就是混沌(chaos):確定的方程,不可預測的未來。
而 1954 年的 KAM 定理帶來了一線希望:不可積系統並非全盤皆輸—— 大部分的可積結構會存活下來,只有「共振」的地方才會碎成混沌。 相空間是「可積的島嶼」與「混沌的海洋」交織的圖案。
正規理論#
可積性的定義#
\(n\) 自由度哈密頓系統若存在 \(n\) 個互相獨立且兩兩對合的守恆量
則稱可積(Liouville 可積)。此時:
- 運動方程可用正則變換「完全積出」(Liouville–Arnold 定理)。
- 相空間分解成不變環面(invariant tori): \(F_i\) = 常數定義的 \(n\) 維環面,軌跡在環面上纏繞。
- 可引入作用量–角變數 \((J_i, \theta_i)\) :
則 \(\dot J_i = 0\) 、 \(\dot\theta_i = \omega_i(J) = \dfrac{\partial H}{\partial J_i}\) —— 環面上的均勻直線運動。
一維系統都可積( \(H\) 本身就是那個守恆量); 二維以上不一定(三體問題是最著名的反例)。
微擾與共振#
近可積系統: \(H(J, \theta) = H_0(J) + \epsilon H_1(J, \theta)\) 。 用正則變換消去角依賴時,會遇到分母為零的項——共振:
共振處微擾被放大,環面被破壞——這是混沌的入口。
KAM 定理#
Kolmogorov–Arnold–Moser 定理(1954 提出,1963 完成證明):
若 (1) \(H_0(J)\) 非退化( \(\det \frac{\partial^2 H_0}{\partial J_i \partial J_j} \neq 0\) ), (2) 微擾夠小( \(|\epsilon|\) 足夠小),(3) 頻率向量 \(\omega\) 滿足丟番圖條件 (無理到「夠遠離共振」):
則對應的不變環面存活(僅微幅變形)。
而頻率有理數(共振)附近的環面碎裂——相空間變成:
- 活存的 KAM 環面(多數,測度隨 \(\epsilon \to 0\) 趨於全體)
- 共振區的混沌層與週期軌道鏈
混沌的特徵#
- 對初始條件敏感:相鄰軌跡指數分離 \(\|\delta x(t)\| \sim e^{\lambda t}\) , \(\lambda > 0\) 是** Lyapunov 指數**。
- 混合:相空間中的區域被拉長、摺疊(麵團拉伸–摺疊模型),長期後到處都是。
- 同宿纏結:Poincaré 發現,雙曲型週期軌道的穩定與不穩定流形無限次相交, 軌跡因此無限複雜——混沌的第一張照片。
簡易範例:二維標準映射#
離散化的轉子: \((\theta, p) \to (\theta + p',\ p + K\sin(\theta + p'))\) (模 \(2\pi\) )。 參數 \(K\) 控制混沌程度:
- \(K = 0\) :可積, \(p\) 守恆,相圖是水平直線。
- \(K = 0.5\) :微擾,多數 KAM 環面存活,相圖是「直線微彎」。
- \(K = 2\) :大範圍混沌海,只有零星島嶼(共振區的週期島)。
一張圖看盡 KAM 定理:可積 → 微擾存活 → 共振碎裂成混沌海。
歷史典故#
- Poincaré 發現三體問題不可積與同宿纏結,見 1890-Poincare三體問題.md——他為此贏得瑞典國王的大獎,卻在交稿後發現錯誤、重寫了整篇論文。
- KAM 定理:Kolmogorov 1954 年提出構想,見 1954-KAM定理.md;Arnold 與 Moser 1963 年完成證明,見 1963-HeonHeiles混沌.md。
- Hénon–Heiles 1963 年用電腦畫出典型混沌相圖(星系中恆星運動模型),見 1963-HeonHeiles混沌.md。
- FPUT 實驗(1955):電腦模擬非線性晶格,發現能量不均分與復發現象——最早的「數值實驗發現可積結構」,見 1955-FPUT實驗.md。
習題#
- 證明一維系統 \(H(q,p)\) 必可積(提示: \(H\) 本身就是守恆量, \(n=1\) )。
- 二維各向同性諧振子 \(H = \frac{1}{2}(p_x^2 + p_y^2) + \frac{1}{2}(x^2 + y^2)\) , 找出兩個互相對合的守恆量(提示: \(H\) 與角動量 \(L = xp_y - yp_x\) ,驗證 \(\{H, L\} = 0\) )。
- 在標準映射中取 \(K = 0.5\) ,說明為何大多數 KAM 環面存活(頻率滿足丟番圖條件)。