[論文レビュー] Linear Multistep Numerical Methods for Ordinary Differential Equations

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

Numerical methods for differential equations被引用数 4

ひとこと要約

この論文は、常微分方程式を解くための線形多段法について包括的なレビューを提供しており、アダムス・バッシュフォースト（陽的）法、アダムス・ムールトン（陰的）法、および後退差分公式（BDF）法に焦点を当てる。これらの方法は、右辺の多項式近似によってストーン＝ヴァイエルシュトラスの定理から導出され、次数、一貫性、安定性（0-安定性、A-安定性、L-安定性）を分析し、陽的法ではA-安定性が達成できないこと、およびA-安定性を満たす次数が2より大きい方法は存在しないことを確立する。

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.

研究の動機と目的

• 非剛性および剛性ODEを解くために広く用いられる線形多段法を体系的に導出し、分析すること。
• 線形多段スキームの文脈において、次数、一貫性、およびさまざまな安定性概念（0-安定性、A-安定性、L-安定性）の理論的基盤を確立すること。
• 精度と安定性特性の観点から、陽的（アダムス・バッシュフォースト）、陰的（アダムス・ムールトン）、およびBDF法の利点と欠点を比較すること。
• 線形多段法の制限、たとえば陽的スキームではA-安定性が達成できないこと、およびA-安定性を満たす方法の次数の上限があることについて明確にすること。

提案手法

• k個の過去の点で f(t,y) を次数 ≤k−1 の多項式で補間し、その近似を区間 [t_{n-1}, t_n] で積分することで、アダムス・バッシュフォースト法を導出する。
• f_n を含む陰的多項式補間を用いることで、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 アダムス・バッシュフォースト法とアダムス・ムールトン法は、過去のステップにおける f(t,y) の多項式補間によってどのように導出されるか？
• RQ2 線形多段法における係数 α_j と β_j が一貫性および次数 p を満たすために満たすべき条件は何か？
• RQ3 なぜ陽的線形多段法では A-安定性を達成できないのか？
• RQ4 L-安定性と剛性減衰の関係は何か？ なぜ、次数が低いにもかかわらず、後退オイラー法がトラペゾイド則よりも剛性問題に対して優れているのか？
• RQ5 A-安定な線形多段法の最大次数は何か？ そして、この上限が存在する理由は何か？

主な発見

• アダムス・バッシュフォースト法は陽的で、k段法では次数 k である。k 個の過去の点における f の補間を行い、その補間多項式を [t_{n-1}, t_n] で積分することで導出される。
• アダムス・ムールトン法は陰的で、k段法では次数 k+1 である。同じステップ数でも、アダムス・バッシュフォースト法よりも高い精度を達成する。
• BDF法はk段法で次数kであり、k ≤ 2 の場合にA-安定であるが、k > 2 の場合にはそうではない。
• A-安定な線形多段法は次数が2より大きくはならず、次数2のトレイプゾイド則（誤差定数 C_3 = 1/12）が、A-安定性を満たす中で最も精度の高い方法である。
• L-安定性は、Re(z) → −∞ のとき |R(z)| → 0 となる条件を満たす。このため、後退オイラー法（L-安定）は、次数が低いにもかかわらず、剛性問題に対してトレイプゾイド則（A-安定だがL-安定でない）よりも優れた性能を示す。
• 次数 p の方法の局所切り捨て誤差は O(h^{p+1}) であり、誤差定数 C_{p+1} が主係数を決定する。