線形多段階法

線形多段階法は、いくつかの以前の解の値と導関数の線形結合から新しい解の値を計算し、過去の作業を再利用することで、ステップあたりのコストを低く抑えつつ高次の精度を達成します。

Definition

線形多段階法は、いくつかの以前の解の値と右辺の評価値の間の固定された線形関係を通じて次の解の値を決定する常微分方程式の解法です。

Scope

このトピックでは、Adams-Bashforth（陽的）およびAdams-Moulton（陰的）ファミリー、スティッフ問題に対する後退差分公式、予測子-修正子実装、ゼロ安定性を定義する特性多項式と根の条件、およびこれらの方法が達成できる限界を定めるDahlquistの次数障壁について扱います。

Core questions

• 多段階法は、どのようにして過去の値を再利用し、ステップごとに1回の新しい関数評価で高次精度を達成するのでしょうか？
• ゼロ安定性とは何ですか？また、特性多項式の根の条件はそれをどのように表現するのでしょうか？
• 予測子-修正子ペアは、陽的公式と陰的公式を実際にどのように組み合わせるのでしょうか？
• Dahlquistの次数障壁は、多段階法の精度と安定性の限界について何を述べているのでしょうか？

Key theories

ゼロ安定性と根の条件

多段階法は、その第一特性多項式の根が閉単位円盤内にあり、境界上の根が単純根である場合に限り、ゼロ安定であり、したがって一貫性がある場合に収束します。この根の条件は、安定性の多段階法における類似物です。

Dahlquistの障壁

Dahlquistの第一障壁は、ゼロ安定なk段階法の次数に上限を設け、第二障壁は、A安定な線形多段階法が2次を超える次数を持つことはできないことを示しています。このため、高次のスティッフソルバーは、絶対安定性ではなく相対安定性というBDFの妥協に依存しています。

Mechanisms

Adams法は、過去の導関数値を通る補間多項式を積分します。Adams-Bashforthは既知の値のみを使用し（陽的）、Adams-Moultonは未知の新しい値を含み（陰的）、より高い精度と安定性を実現します。実際には、これら2つは予測子-修正子としてペアで使用されます。陽的公式が予測を行い、陰的公式が修正を行います（通常は1〜2回の反復）。一方、後退差分公式は、過去の解の値を差分して新しい点での導関数を近似し、スティッフODEコードの中心となるスティッフ安定な方法を提供します。多段階法はいくつかの初期値を必要とするため、一段階法によってブートストラップされます。

Clinical relevance

線形多段階法、特に後退差分公式は、化学反応速度論、電子回路シミュレーション、および大規模な微分代数系で使用される実用的なスティッフODEソルバーの基盤となっています。これらの分野では、右辺の評価が高価であり、多段階公式を通じて過去の評価を再利用することで、大幅な効率向上が得られます。

History

AdamsとBashforthは19世紀に多段階公式を導入し、Moultonが陰的変種を追加しました。Dahlquistの1950年代から60年代の分析は、この分野を支配する安定性理論と次数障壁を確立し、C. William Gearの1970年代の研究により、後退差分公式コードがスティッフ問題の標準となりました。

Key figures

• John Couch Adams
• Francis Bashforth
• Forest Ray Moulton
• Germund Dahlquist
• C. William Gear

Seminal works

• hairer1993
• iserles2008

Frequently asked questions

多段階法はRunge-Kutta法とどう異なりますか？

Runge-Kutta法は、各ステップ内でいくつかの新しい導関数評価を行いますが、その後それらを破棄します。一方、多段階法は、以前のステップからの導関数値を再利用します。したがって、多段階法はステップあたりのコストは低いですが、追加の初期値とステップサイズ変更の特別な処理が必要です。

根の条件とは何ですか？

これは、方法の第一特性多項式の根が単位円の内側または上にあり、境界上の根が単純であるという要件です。これにより、ステップが累積しても小さな誤差が増幅されないことが保証され、方法がゼロ安定であり、したがって収束することが保証されます。