[論文レビュー] A Second Order Approximation for the Caputo Fractional Derivative

Yuri Dimitrov|arXiv (Cornell University)|Feb 2, 2015

Fractional Differential Equations Solutions参考文献 34被引用数 15

ひとこと要約

本稿では、リーマン・ゼータ関数を用いて、古典的な1次精度の公式の最初の3つの係数を修正することにより、キャプート・分数階微分の2次精度近似を提案する。新しいスキームは、すべての $\alpha \in (0,1)$ に対して $O(h^2)$ 収束を達成し、特に $\alpha \approx 1$ の近傍で顕著に精度が向上し、分数階緩和方程式および分数階拡散方程式の数値的解法において優れた性能を示す。

ABSTRACT

When $0

研究の動機と目的

• 当該 $\alpha \approx 1$ のとき、収束が $O(h)$ に低下する古典的キャプート微分近似における精度損失を是正すること。
• すべての $\alpha \in (0,1)$ の範囲で高精度を維持する2次精度の数値スキームを、キャプート分数階微分に対して開発すること。
• 時間方向離散化の次数を向上させることで、特に分数階拡散方程式および緩和方程式の数値的解法を改善すること。
• 漸近展開と畳み込み・クワドラチャ・重みへのゼータ関数補正を用いて、改善された収束速度の理論的基盤を提供すること。

提案手法

• 標準的な畳み込み・クワドラチャにおける係数 $\sigma_k^{(\alpha)}$ の最初の3つを、リーマン・ゼータ関数 $\zeta(\alpha-1)$ を用いて修正した $\delta_k^{(\alpha)}$ に置き換えることで、キャプート微分の2次近似を導出する。
• 和 $\sum_{k=0}^n k^{\alpha}$ の漸近展開を用いて、一次誤差項を消去する補正を導出する。
• 修正された重みを $\delta_0^{(\alpha)} = \sigma_0^{(\alpha)} - \zeta(\alpha-1)$, $\delta_1^{(\alpha)} = \sigma_1^{(\alpha)} + 2\zeta(\alpha-1)$, $\delta_2^{(\alpha)} = \sigma_2^{(\alpha)} - \zeta(\alpha-1)$ と定義し、$k \geq 3$ では $\delta_k^{(\alpha)} = \sigma_k^{(\alpha)}$ とする。
• 新しい近似を用いて、均等グリッド上の有限差分スキームを用いて、分数階緩和方程式および分数階拡散方程式を解く。
• 古典的 $O(h^{2-\alpha})$ スキームと新しい $O(h^2)$ スキームを比較し、誤差表および収束順序推定を用いて収束率を検証する。

実験結果

リサーチクエスチョン

• RQ1 すべての $\alpha \in (0,1)$ において、特に $\alpha \approx 1$ の近傍で、古典的一次精度キャプート微分近似の精度を2次精度に向上させることは可能か？
• RQ2 標準的な畳み込み・クワドラチャ近似におけるキャプート微分の一次誤差項を消去するために必要な補正項は何か？
• RQ3 修正された近似は、分数階拡散方程式および緩和方程式の数値的解法における収束速度および安定性にどのように影響を与えるか？
• RQ4 リーマン・ゼータ関数を用いて畳み込み重みを体系的に補正し、分数階数範囲全体にわたって一様な2次収束を達成することは可能か？

主な発見

• 提案された2次精度近似 $y^{(\alpha)}(x) = \frac{1}{\Gamma(2-\alpha)h^{\alpha}} \sum_{k=0}^n \delta_k^{(\alpha)} y(x-kh) + O(h^2)$ は、すべての $\alpha \in (0,1)$ に対して一様に2次収束を達成し、$\alpha = 1$ の近傍での $O(h)$ 劣化を解消する。
• 分数階緩和方程式（$\alpha = 0.6$）に対する数値実験では、新しいスキームが $\tau = h/2$ 時に収束順序 $\approx 1.98$ を達成し、古典的スキームの $\approx 1.4$ より顕著に向上している。
• 時間分数階拡散方程式（$\alpha = 0.4$）において、新しいスキームは $\tau = h/2$ 時に収束順序 $1.997$ を達成したが、古典的スキームは $1.56$ にとどまっていた。
• $\alpha = 0.6$ 時に、新しいスキームは拡散方程式に対して収束順序 $2.035$ を達成したのに対し、古典的スキームは $1.396$ にとどまり、一貫した2次収束の挙動を示した。
• 最初の3係数に対するリーマン・ゼータ関数の使用により、修正された重みは和がゼロとなる性質を維持し、漸近誤差展開の改善が図られた。
• 特に長時間シミュレーションや高階時間微分を含む方程式では、古典的スキームの誤差が急速に蓄積されるため、その精度向上は顕著である。