[论文解读] Linear Multistep Numerical Methods for Ordinary Differential Equations

Nikesh S. Dattani|ArXiv.org|Oct 28, 2008

Numerical methods for differential equations被引用 4

一句话总结

本文全面回顾了求解常微分方程的线性多步法，重点讨论了显式 Adams-Bashforth 方法、隐式 Adams-Moulton 方法以及后向微分公式（BDF）方法。通过右端项的多项式逼近，基于 Stone-Weierstrass 定理推导出这些方法，分析了其阶数、相容性及稳定性（包括 0-稳定性、A-稳定性与 L-稳定性），并证明显式方法无法实现 A-稳定性，且不存在阶数高于 2 的 A-稳定方法。

ABSTRACT

A review of the most popular Linear Multistep (LM) Methods for solving Ordinary Differential Equations numerically is presented. These methods are first derived from first principles, and are discussed in terms of their order, consistency, and various types of stability. Particular varieties of stability that may not be familiar, are briefly defined first. The methods that are included are the Adams-Bashforth Methods, Adams-Moulton Methods, and Backwards Differentiation Formulas. Advantages and disadvantages of these methods are also described. Not much prior knowledge of numerical methods or ordinary differential equations is required, although knowledge of basic topics from calculus is assumed.

研究动机与目标

• 系统推导并分析求解非刚性与刚性常微分方程中最广泛使用的线性多步法。
• 在多步格式背景下，建立阶数、相容性以及各种稳定性概念（0-稳定性、A-稳定性、L-稳定性）的理论基础。
• 从精度与稳定性特性角度，比较显式（Adams-Bashforth）、隐式（Adams-Moulton）与 BDF 方法的优缺点。
• 阐明线性多步法的局限性，如显式格式无法实现 A-稳定性，以及 A-稳定方法存在阶数障碍。

提出的方法

• 通过在前 k 个点上对函数 f(t,y) 进行次数 ≤k−1 的多项式插值，并在区间 [t_{n-1}, t_n] 上积分该多项式近似，推导 Adams-Bashforth 方法。
• 通过包含 f_n 的隐式多项式插值推导 Adams-Moulton 方法，从而实现更高阶精度（k 步方法的阶数为 k+1）。
• 应用测试方程 y' = λy 分析绝对稳定性，定义稳定性多项式 ζ(ξ) = ρ(ξ) − hλσ(ξ)，并将稳定性区域定义为满足 |ξ| ≤ 1 的 z = λh 的集合。
• 通过特征多项式 ρ(ξ) 的根条件定义 0-稳定性，要求 |ξ_i| ≤ 1，且单位圆上的根为单根。
• 将 A-稳定性定义为稳定性区域包含整个左半平面（Re(λh) < 0）的性质，L-稳定性则定义为当 Re(z) → −∞ 时 |R(z)| → 0。
• 利用误差常数 C_{p+1} 衡量局部截断误差，并证明阶数为 p 的方法其误差为 O(h^{p+1})。

实验结果

研究问题

• RQ1 Adams-Bashforth 与 Adams-Moulton 方法如何从 f(t,y) 在前序步骤上的多项式插值中导出？
• RQ2 系数 α_j 与 β_j 需满足何种条件，才能保证线性多步法的相容性与阶数 p？
• RQ3 为何显式线性多步法无法实现 A-稳定性？
• RQ4 L-稳定性与刚性衰减之间有何关系？为何尽管阶数较低，后向欧拉法仍优于梯形法则？
• RQ5 A-稳定线性多步法的最大阶数是多少？该限制存在的原因是什么？

主要发现

• Adams-Bashforth 方法为显式方法，k 步格式的阶数为 k，通过在前 k 个点插值 f 并积分插值多项式推导得出。
• Adams-Moulton 方法为隐式方法，k 步格式的阶数为 k+1，与相同步数的 Adams-Bashforth 方法相比具有更高精度。
• BDF 方法为 k 步格式，阶数为 k，当 k ≤ 2 时为 A-稳定，但当 k > 2 时不再 A-稳定。
• A-稳定的线性多步法阶数不能超过 2，且梯形法则（阶数为 2）是误差常数 C_3 = 1/12 的最精确 A-稳定方法。
• 当 Re(z) → −∞ 时 |R(z)| → 0 时实现 L-稳定性，这解释了为何后向欧拉法（L-稳定）在刚性问题中表现优于梯形法则（A-稳定但非 L-稳定）。
• 阶数为 p 的方法的局部截断误差为 O(h^{p+1})，误差常数 C_{p+1} 决定其主导系数。