[論文レビュー] A Superexponentially Convergent Functional-Discrete Method for Solving the Cauchy Problem for Systems of Ordinary Differential Equations

Makarov Volodymyr, Denis Dragunov|arXiv (Cornell University)|Dec 30, 2010

Fractional Differential Equations Solutions参考文献 6被引用数 3

ひとこと要約

本稿では、常微分方程式（ODE）の1階系系（SODE）のコーシー問題を解くための、超指数的収束性を示す関数的離散（FD）法を提示する。この手法は、アドミアン多項式に基づく再帰的スキームと区分的定数離散化グリッドを用いて、FD法をSODEに一般化したものであり、系の非線形性および外力項に十分な条件が満たされれば、超指数的収束を達成する。主な貢献は、誤差の減衰率に関する定量的推定を含む、厳密な収束定理の確立である。

ABSTRACT

In the paper a new numerical-analytical method for solving the Cauchy problem for systems of ordinary differential equations of special form is presented. The method is based on the idea of the FD-method for solving the operator equations of general form, which was proposed by V.L. Makarov. The sufficient conditions for the method converges with a superexponential convergence rate were obtained. We have generalized the known statement about the local properties of Adomian polynomials for scalar functions on the operator case. Using the numerical examples we make the comparison between the proposed method and the Adomian Decomposition Method.

研究の動機と目的

• 1階系ODE系のコーシー問題を解くための、優れた収束特性を有する数値解析的手法の開発。
• 従来固有値問題や境界値問題に適用されてきた関数的離散（FD）法を、SODEの初期値問題に拡張すること。
• SODEのFD法の超指数的収束を保証する十分条件の確立。
• 特にADMが収束しない場合に、提案手法のFD法とAdomian分解法（ADM）を比較すること。

提案手法

• 解をベクトル関数の部分和の級数として近似し、各項を区分的グリッド上の線形コーシー問題により再帰的に解く。
• 基本問題（j=0）は、直前のグリッド点における非線形作用素の評価から得られる区分的定数係数線形ODEを解く。
• 高次項（j≥1）は、非線形作用素のアドミアン多項式を用いた再帰的公式により計算され、系を逐次線形化する。
• 離散化グリッドはステップサイズhの均等グリッドであり、各グリッド点でジャンプ条件を満たすことで、区間間での解の連続性を保証する。
• 非線形作用素N(t, u)の構造を活用し、アドミアン多項式によるテイラー型展開を用いて非線形性を効率的に扱う。
• 収束解析は、非線形作用素のヤコビアンおよび外力項のバインドに基づき、誤差の超指数的減衰のための十分条件を導出する。

実験結果

リサーチクエスチョン

• RQ1 関数的離散法は、非線形項を有する1階系ODE系に一般化可能であり、その際超指数的収束性を維持するか？
• RQ2 SODEのFD法の超指数的収束を保証する非線形作用素および外力項の十分条件は何か？
• RQ3 ADMと比較して、FD法の収束行動はどのように異なるか、特に発散する場合にどうか？
• RQ4 ADMが不安定性や発散のために失敗する問題に対しても、FD法を適用可能か？

主な発見

• SODEのFD法は、ヤコビアンの固有値のスペクトル的性質および外力項の有界性に依存する十分な条件下で、超指数的収束を示す。
• 数値的に、ADMが発散する非線形系（周期的解を有する）に対しても、FD法は超指数的収束を達成している。
• テスト問題において、ランクp=4のFD法は区間[0,6]で絶対誤差2×10⁻⁶未満を達成し、ランクを増加させることで誤差が指数的に減少した。
• 理論的誤差境界は、ヤコビアンの対称部の最小固有値および外力項のL∞-ノルムを用いて導出され、α≈0.47およびκ≈2.9が得られた。
• 数値結果から、FD法の収束速度がADMよりも顕著に速いことが確認され、同様の問題に対してADMは[0,2]で発散を示した。
• 正確な解u₁(t)=sin(t), u₂(t)=cos(t)が、FD法により高い精度で回復され、理論的収束主張の妥当性が検証された。