跳轉到

1948 — Lax 等價定理：差分格式的第一定律#

案件檔案表#

欄位 內容
案件名稱 收斂性之謎：格點越細解越準，靠的是什麼保證？
主嫌（主角） Peter D. Lax（與 Richtmyer 於 1955-56 定稿發表）
案發年份 1948 年提出，1955-56 與 Richtmyer 發表證明
案發地點 洛斯阿拉莫斯（Los Alamos）國家實驗室與 NYU
被害人 天真信念：格式「看起來對」（一致性）就該收斂
兇器 不穩定格式：誤差隨步數指數放大
關鍵證物 一致性 \(E(\Delta t) \to 0\) 、穩定性 \(\|A^n\| \le C\) 、收斂性 \(u_n \to u\)
結案判決 對適定線性初值問題：一致性＋穩定性 \(\iff\) 收斂性
遺留線索 von Neumann 穩定性分析、CFL 條件、數值 PDE 的法典化

案發現場#

二戰之後，數值計算迎來爆發期：ENIAC 上線（1945-46），Los Alamos 的流體與激波計算、氣象預報（Charney–von Neumann 計畫）都要求把 PDE 離散化。差分格式是主角：把 \(\partial_t u = \partial_x u\) 換成

\[ \frac{u_j^{n+1} - u_j^n}{\Delta t} = \frac{u_{j+1}^n - u_{j-1}^n}{2 \Delta x}, \]

然後在機器上一步步推進。但現場立刻浮出三具疑屍：

  1. 格式 A（FTCS 對流方程）：截斷誤差 \(O(\Delta t + \Delta x^2)\) ，看起來「一致」，卻對任何 \(\Delta t\) 都不穩定——誤差像瘟疫一樣指數傳播。收斂？從未發生。
  2. 格式 B（跳蛙 leapfrog）：條件穩定，必須 \(\lambda = c \Delta t / \Delta x \le 1\) ；違反 CFL 條件立即爆炸。
  3. 理論真空：1948 年之前，「穩定性」一詞人人會說、無人能證；收斂與否全靠數值實驗的試錯。

Lax 在 Los Alamos 與 von Neumann 的討論中把案子理清：一致性只是必要條件，穩定性才是生死線；兩者合起來，恰好等價於收斂。這就是數值分析的第一定律。

偵查過程#

第一步：設定三項罪證的精確定義#

設 \(X\) 為 Banach 空間（例如 \(L^2\) 的格點副本），真解 \(u(t) = A(t) u_0\) 適定；離散推進算子 \(A_{\Delta t}\) 把 \(u^n\) 推到 \(u^{n+1}\) ， \(t_n = n \Delta t\) 。

證物 定義 檢驗方式
一致性 截斷誤差 \(E(\Delta t) = \frac{A_{\Delta t} u(t) - u(t+\Delta t)}{\Delta t} \to 0\) Taylor 展開
穩定性 \(\|A_{\Delta t}^n\| \le C(T)\) 對 \(n \Delta t \le T\) 一致有界 von Neumann 分析或能量法
收斂性 \(n \to \infty\) ， \(\Delta t \to 0\) ： \(\|u^n - u(t_n)\| \to 0\) 定理保證

第二步：收斂 ⇒ 穩定（半邊推理）#

若格式收斂，則對每個初值 \(v\) ， \(A_{\Delta t}^n v\) 有界（收斂到 \(u(t_n; v)\) ）。由 Banach–Steinhaus（一致有界原理，1929），算子族 \(\{A_{\Delta t}^n\}\) 一致有界——穩定性。這半邊用的是泛函分析的共振方法：逐點有界推出一致有界。1948 年的 Lax 深諳此道。

第三步：一致＋穩定 ⇒ 收斂（主推理）#

把誤差按時間切片。 \(u^{n+1} - u(t_{n+1})\) 可拆成兩段：

\[ u^{n+1} - u(t_{n+1}) = A_{\Delta t}\big( u^n - u(t_n) \big) + E(\Delta t). \]

反覆套用並用穩定性控制放大：

\[ \|u^n - u(t_n)\| \le \sum_{k=0}^{n-1} \|A_{\Delta t}^{\,n-1-k}\|\ \|E(\Delta t)\| \le T\, C(T)\, \frac{\|E(\Delta t)\|}{\Delta t} \to 0. \]

推理骨架一目瞭然：穩定性提供「放大不超過常數倍」，一致性提供「每步注入的毒量趨零」；兩者缺一，病人必死。

第四步：von Neumann 穩定性分析（驗屍工具）#

對常係數格式，取 Fourier 模式 \(u_j^n = \xi^n e^{i j \theta}\) ，代回格式得增長因子 \(g(\theta)\) 。穩定性判據為 von Neumann 條件

\[ |g(\theta)| \le 1 + C \Delta t \quad \forall \theta. \]
格式 增長因子 \(g(\theta)\) 判決
FTCS（對流） \(1 - i \lambda \sin\theta\) $
向後歐拉 \((1 + i\lambda \sin\theta)^{-1}\) 的形式 \(\|g\| \le 1\) ：無條件穩定
跳蛙 \(g^2 + 2 i \lambda g \sin\theta - 1 = 0\) \(\lambda \le 1\) 時 $
Lax–Friedrichs \(\cos\theta - i\lambda\sin\theta\) \(\lambda \le 1\) 時 $

注意 FTCS 的教訓： \(|g| = 1 + O(\lambda^2 \sin^2\theta)\ \Delta t^2\) 型的「每步微膨脹」累積 \(n \sim T/\Delta t\) 步後變成 \(e^{C T / \Delta t}\) ——指數爆炸。穩定性必須是 \(\le 1 + C\Delta t\) 這種剛剛好的界。

第五步：1955-56 定稿與適定性前提#

Lax 與 Richtmyer 在《Comm. Pure Appl. Math.》發表完整證明，處理變係數與抽象設定，並明確前提：問題必須適定（well-posed）且格式為線性。對非線性問題，等價性退化为「穩定性 ⇒ 收斂」的單向猜測（後由 Lax–Wendroff、Strang 等案續審）。

結案報告#

判決：Lax 等價定理成立——對適定線性初值問題的相容格式：一致性＋穩定性 \(\iff\) 收斂性。檢驗收斂不必真收斂，只需檢驗兩個更容易抓的嫌犯。

遺產：

  1. 數值分析法典化：從此任何差分格式的論文必附「一致性驗證 + 穩定性證明」，證明義務被制度化。這是本書因果鏈中「Lax等價 (1948) → 數值分析（Lax等價 1948）」的法律基礎。
  2. von Neumann 穩定性分析成為工程標準配備：寫下 \(g(\theta)\) 、畫出 CFL 三角、檢查 \(|g| \le 1\) 。
  3. 與 PDE 理論接軌：收斂性結論需適定性（Hadamard），而橢圓正則性與能量估計（Gårding 1956）提供穩定性的能量證法，兩案互為表裡（見 1950-LaxMilgram定理）。
  4. 現代延伸：非線性守恆律的 Lax–Wendroff 定理、算子分裂的 Strang 收斂理論、機器學習中的神經 ODE 穩定性分析，都是本案的續集。

證據與工具#

工具 用途
一致性（截斷誤差） 每步毒量測量，Taylor 展開檢驗
穩定性 \(\|A_{\Delta t}^n\| \le C\) 放大不失控，格式生死線
Banach–Steinhaus（1929） 收斂 ⇒ 穩定的半邊證明
von Neumann 分析 傅立葉增長因子，常係數格式驗屍刀
適定性（Hadamard） 整案前提：解存在、唯一、連續依賴

範例：對流方程 \(u_t = u_x\) ， \(\lambda = \Delta t / \Delta x\) 。取 \(\lambda = 0.9\) ：Lax–Friedrichs 穩定且收斂；跳蛙穩定且收斂；FTCS 一致但不穩定，數值解在幾百步後被捨入誤差的放大產物淹沒。把 \(\lambda\) 加到 \(1.1\) ：跳蛙爆炸，違反 CFL。同一個方程，三種結局，全由穩定性一項證物決定——這就是第一定律的教科書演示。

一句話結案：一致性說「每步走對方向」，穩定性說「別被誤差絆倒」；兩者同時成立，收斂便是必然的判決。