4. 稳定性定义与Lax等价定理

在有限差分方法的理论体系中，稳定性不是一个可有可无的附加条件，而是连接“格式能算”与“算得对”的桥梁。前面我们讨论了截断误差与相容性，说明了一个格式在步长趋于零时是否逼近原微分方程；但即便格式相容，如果误差在推进过程中被无限放大，数值解依然会偏离真解。稳定性正是刻画“误差是否被控制”的核心概念。本节严格区分零稳定性与绝对稳定性，陈述Lax等价定理，并厘清稳定性、相容性与收敛性三者的逻辑关系。

一、零稳定性：格式对初始扰动的敏感程度

零稳定性（zero-stability）考察的是当步长趋于零时，差分格式本身对初始数据扰动的放大行为。考虑线性多步法的一般形式：

Σ_{j=0}^{k} α_j u^{n+j} = Δt · Σ_{j=0}^{k} β_j f^{n+j}

其中 α_j、β_j 为格式系数，α_k ≠ 0。零稳定性只与左端系数 α_j 有关，与右端 β_j 及右端项 f 无关。其严格定义是：存在常数 C 和 Δt₀，使得对任意初始扰动序列，当 0 < Δt ≤ Δt₀ 时，数值解扰动的范数不超过 C 乘以初始扰动范数。等价地，第一特征多项式 ρ(ζ) = Σ α_j ζ^j 的所有根满足根条件：|ζ| ≤ 1，且 |ζ| = 1 的根必须是单根。

根条件给出了非常直观的判据。以显式Euler格式 u^{n+1} = u^n + Δt f^n 为例，ρ(ζ) = ζ − 1，根为 ζ = 1，是单根且模为1，故零稳定。而格式 u^{n+1} = −u^n + 2Δt f^n 的 ρ(ζ) = ζ + 1，根为 ζ = −1，模为1且为单根，同样零稳定，但它的解会以振荡方式传播扰动。若出现 |ζ| > 1 的根，则扰动随 n 指数增长，格式零不稳定。需要强调：零稳定性是格式的“结构性”属性，与具体 Δt 取值无关，只要求 Δt 足够小。

二、绝对稳定性：固定步长下的误差衰减

零稳定性回答的是“Δt → 0 时格式是否可靠”，而实际计算中 Δt 是有限值。绝对稳定性（absolute stability）针对固定 Δt，考察格式对误差的衰减能力。标准做法是代入测试方程 u' = λu（λ 为复数），令 u^n = ζ^n，得到放大因子或放大矩阵，要求 |ζ| ≤ 1。使 |ζ| ≤ 1 成立的 λΔt 在复平面上的集合称为绝对稳定区域。

以显式Euler格式为例，放大因子 ζ = 1 + λΔt。当 λ 为负实数时，要求 |1 + λΔt| ≤ 1，即 −2 ≤ λΔt ≤ 0，稳定区间为 [−2, 0]。这意味着对于刚性问题（|λ| 很大），显式Euler必须取极小步长才能稳定，这正是其绝对稳定性受限的体现。隐式Euler格式的放大因子为 ζ = 1/(1 − λΔt)，对任意 Re(λ) < 0 都有 |ζ| < 1，属于A稳定格式。绝对稳定性是“格式 + 步长”的联合属性，同一格式在不同 Δt 下可能稳定也可能不稳定。

三、Lax等价定理及其意义

Lax等价定理（Lax Equivalence Theorem）是有限差分理论中最著名的结论之一，其表述为：对于一个适定的线性初值问题，若差分格式是相容的，则稳定性与收敛性等价。用逻辑式表达即：

相容 + 稳定  ⇔  收敛

该定理的证明思路是“误差方程”技巧：设 e^n = u^n − U^n（U^n 为真解在网格点上的值），代入格式后得到误差满足与格式同构的齐次方程，右端为截断误差 τ^n。于是 e^n 可表示为截断误差的累积传播：e^n = Σ G^{n−m} τ^m，其中 G 为放大算子。相容性保证 τ → 0，稳定性保证 G 的幂一致有界，两者结合即得 e^n → 0。反过来，若格式收敛，则可构造性证明其稳定。

Lax等价定理的意义在于把“证明收敛”这一困难任务拆解为两个相对可操作的子任务：相容性通过Taylor展开即可验证，稳定性通过根条件或放大因子分析即可判断。它同时给出了负面警示：相容但不稳定的格式一定不收敛，此时加密网格不仅不能改善精度，反而可能使解发散。因此在实际编程中，先验证相容性，再检查稳定性，是保证计算结果可信的基本流程。

四、三者的逻辑关系与工程判读

综合来看，相容性、稳定性与收敛性构成一条清晰的逻辑链。相容性是“格式逼近方程”的必要条件，回答“方程对不对”；稳定性是“误差不放大”的保证，回答“计算稳不稳”；收敛性则是最终目标，回答“结果准不准”。Lax等价定理指出，在适定线性问题的前提下，相容性与稳定性联合起来恰好等价于收敛性，三者并非彼此独立，而是通过误差传播方程紧密耦合。

在工程实践中，这一逻辑关系可以直接转化为检查流程。拿到一个格式后，第一步做相容性分析：将真解代入格式，做Taylor展开，确认截断误差阶数，若 τ = O(Δt^p) 则相容。第二步做稳定性分析：写出第一特征多项式检验根条件，判断零稳定性；再代入测试方程求出放大因子，画出绝对稳定区域，结合实际问题中 λΔt 的取值范围判断是否落在稳定域内。第三步才谈收敛阶估计：在稳定域内，全局误差阶通常与截断误差阶一致。若跳过第二步直接加密网格，遇到零不稳定格式时误差会随步数指数增长，表现为数值解出现非物理的高频振荡或溢出，此时应回退检查格式系数与根条件，而非盲目减小 Δt。

以中心差分格式 u^{n+1} = u^{n−1} + 2Δt f^n 为例，其 ρ(ζ) = ζ² − 1，根为 ±1，均为单根，零稳定成立；但代入测试方程后放大因子满足 ζ² − 1 − 2λΔt ζ = 0，对纯虚 λ 其模恒为1，对 Re(λ) < 0 则存在模大于1的根，绝对不稳定。这说明零稳定并不蕴含绝对稳定，实际使用时必须结合问题类型与步长共同判断。掌握这一层次关系，才能在格式选择、步长设计与结果校验中做出有依据的决策。