[論文レビュー] IFOHAM-an iterative algorithm based on the first-order equation of HAM: exploratory preliminary results

Miguel Moreira|arXiv (Cornell University)|Jan 1, 2017

Fractional Differential Equations Solutions被引用数 1

ひとこと要約

本稿では、ホモトピー摂動法（HAM）の一次方程式に基づく反復アルゴリズムであるIFOHAMを提案する。この手法はピカード＝リンドェルフの反復を一般化し、制御パラメータ$c_0$を用いて収束を保証する。非線形初期値問題の解法において、収束速度が優れており、CPU時間も低く抑えられ、特に$c_0 = -1.2$のとき、標準HAMおよびピカード＝リンドェルフ法を上回る性能を示す。

ABSTRACT

In this work we present and study an iterative algorithm used to asymptotically solve nonlinear differential equations. This algorithm (Iterative First Order HAM or IFOHAM) is based on the first order equation of the Homotopy Analysis Method, HAM. We show that IFOHAM generalizes Picard-Lindeloff's iteration algorithm. Moreover, IFOHAM shares with HAM the possibility of ensuring convergence by adequately choosing c0, a convergence control parameter. Preliminary results show that IFOHAM exhibits a very good performance both in aspects related to the speed of convergence and in aspects related to the CPU calculation time. It should also be noted that the IFOHAM is a very low complexity algorithm easily programmable in a symbolic computing environment.

研究の動機と目的

• 非線形常微分方程式を解くための、一次HAM式に基づく新しい反復アルゴリズムIFOHAMの開発。
• 収束制御パラメータ$c_0$を組み込むことで、ピカード＝リンドェルフの反復を一般化し、より広い収束特性を実現すること。
• 標準HAMおよびピカード＝リンドェルフ法と比較して、IFOHAMの収束速度、計算効率、実装の容易さを評価すること。
• $c_0$が収束行動に与える影響を調査し、特定の非線形問題に対して最適なパラメータ選定を同定すること。
• 記号計算環境におけるIFOHAMの実用的妥当性を評価し、解の解析的近似に向けた可能性を検討すること。

提案手法

• 線形作用素$L$、収束制御パラメータ$c_0$、埋め込みパラメータ$q \in [0,1]$を用いて、零次近似の変形方程式を定式化する。
• 一次HAM方程式を導出：$m \geq 0$に対して$L[u_{m+1}(t)] = c_0 \left[ N\left( \sum_{k=0}^m u_k(t) \right) \right]$ であり、これは反復更新則を定義する。
• 初期推定値$u_0(t)$を、問題の初期条件を満たすように選び、順次高次項$u_m(t)$を再帰的に計算する。
• 多項式項の導出および評価を自動化するために、記号計算ツール（例：Mathematica, Maple）を適用する。
• パラメータ$c_0$による収束制御を行い、理論的および数値的証拠から、最適な値が$] -1, 0[$内、あるいはすなわち$c_0 < -1$である可能性が示唆される。
• 標準HAMおよびピカード＝リンドェルフ法との性能比較として、$x' = 1 + x^2$, $x(0) = 0$というベンチマーク非線形IVPを用いる。この問題の正確な解は$\tan t$である。

実験結果

リサーチクエスチョン

• RQ1 IFOHAMは、HAMフレームワーク内でピカード＝リンドェルフの反復アルゴリズムをどのように一般化するか？
• RQ2 収束制御パラメータ$c_0$は、IFOHAMにおける収束の保証および収束速度の向上にどのような役割を果たすか？
• RQ3 ベンチマーク非線形IVPに対して、IFOHAMは標準HAMおよびピカード＝リンドェルフ法と比較して、収束速度および計算コストの点でどのように異なるか？
• RQ4 パラメータ$c_0$を$-1$未満にチューニングすることで、IFOHAMはHAMおよびピカード＝リンドェルフ法を上回る速い収束を達成できるか？
• RQ5 特に高次項における式の複雑さに起因する、IFOHAMの実用的限界は何か？

主な発見

• IFOHAMは、$c_0 = -1$のとき、ピカード＝リンドェルフの反復に正確に一致し、この古典的手法の一般化であることが示された。
• ベンチマークIVP $x' = 1 + x^2$, $x(0) = 0$に対して、$c_0 = -1.2$のIFOHAMは4次近似で残差誤差$5.45 \times 10^{-6}$を達成し、標準HAM（$c_0 = -1$）およびピカード＝リンドェルフ法を上回った。
• IFOHAMを用いた4次近似解の計算に要したCPU時間は2.078秒であり、複雑な式であっても計算オーバーヘッドが低く抑えられていた。
• 収束は$c_0 \in ]-1, 0[$に対して保証され、非線形作用素$f$の構造に応じて$c_0 < -1$のときより速い収束が観察された。
• アルゴリズムの実装が非常に簡単であり、低計算複雑性であるため、記号計算環境に非常に適している。
• 高次反復における式の長さが増加しても、IFOHAMは収束速度および計算効率を強く維持しており、実用的妥当性が裏付けられた。