偏微分方程的數值方法:有限差分

拉克斯等價定理(Lax equivalence theorem)

/ lahks /

用數值法解 PDE 時,你真正在乎的是「收斂」:當網格細化時,數值答案是否趨近真解?收斂難以直接檢查——你得知道真解才行。拉克斯等價定理是讓你繞過這點的美麗結果:它把困難的問題(會收斂嗎?)換成兩個你能分別檢查的容易問題(相容嗎?穩定嗎?)。

定理陳述:對一個適定的線性初值問題與一個相容的有限差分格式,該格式「收斂」,若且唯若它「穩定」。化成口號就是:相容加穩定等於收斂。相容性用泰勒級數驗證(格式在 h、k -> 0 時趨向 PDE)。穩定性用馮諾伊曼分析驗證(誤差不被無界放大)。兩者單獨都不夠:相容但不穩定的格式發散,穩定但不相容的格式收斂到錯誤答案。唯有兩者並存才保證數值解收斂到真解,而且它們互相升級——一旦你知道相容性,穩定性就「等價於」收斂。

這是整個學科的組織原則。它正是為何實務者把心力花在證明穩定性(真正微妙的部分),而非嘗試直接的收斂證明:依此定理,對相容的格式,穩定性即收斂。誠實的細則:它陳述於「線性」、適定的問題,並假設相容性在合適的範數下成立。對非線性 PDE(震波、紊流)沒有這麼乾淨的等價,穩定性概念必須加強(全變差界、熵條件)。但在其線性範圍內,拉克斯定理是把有限差分中其他一切繫在一起的橋。

熱方程的 FTCS 相容(階 k + h^2)。馮諾伊曼分析說它恰當 r = alpha*k/h^2 <= 1/2 時穩定。依拉克斯定理,FTCS 因此「收斂」,若且唯若你保持 r <= 1/2——不需另作收斂證明。違反此界則格式不穩定,故依定理也不收斂。

相容 + 穩定 = 收斂:證明兩個容易的部分,免費得到困難的那個。

此等價是對「線性」、適定的問題。對非線性守恆律(震波、可壓縮流),相容加線性穩定「不」保證收斂到正確(滿足熵)的解;你需要更強的工具,如 Lax-Wendroff 定理與 TVD/熵條件。別過度套用那乾淨的線性口號。

又称
Lax-Richtmyer theoremconsistency + stability = convergence拉克斯-里希特邁耶定理