[論文レビュー] A Tunably-Accurate Laguerre Petrov-Galerkin Spectral Method for Multi-Term Fractional Differential Equations on the Half Line

Anna Lischke, Mohsen Zayernouri+1 more|arXiv (Cornell University)|Jul 28, 2016

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

ひとこと要約

本稿では、半直線上の導関数次数が1未満の多項式分数微分方程式を解くために、調整可能な精度を持つラゲル・ペトロフ・ガレルキンスペクトル法を提案する。一般化関連ラゲル関数（GALFs）を基底関数として用い、再帰関係を活用することで、スペクトル収束を達成し、良好に条件付けられたスパース行列を用いてO(N log N)の計算量で線形方程式系を解く。

ABSTRACT

We present a new tunably-accurate Laguerre Petrov-Galerkin spectral method for solving linear multi-term fractional initial value problems with derivative orders at most one and constant coefficients on the half line. Our method results in a matrix equation of special structure which can be solved in $\mathcal{O}(N \log N)$ operations. We also take advantage of recurrence relations for the generalized associated Laguerre functions (GALFs) in order to derive explicit expressions for the entries of the stiffness and mass matrices, which can be factored into the product of a diagonal matrix and a lower-triangular Toeplitz matrix. The resulting spectral method is efficient for solving multi-term fractional differential equations with arbitrarily many terms. We apply this method to a distributed order differential equation, which is approximated by linear multi-term equations through the Gauss-Legendre quadrature rule. We provide numerical examples demonstrating the spectral convergence and linear complexity of the method.

研究の動機と目的

• 分数微分の非局所的性質に起因する、多項式分数微分方程式を解く際の高い計算コストを低減すること。
• 半直線上の導関数次数≤1の線形多項式分数初期値問題（FIVPs）に対して、効率的なスペクトル法を開発すること。
• ガウス＝レジェンドル求積法を用いて分布型微分方程式を近似し、その結果得られる多項式FIVPsを解くことで、分布型微分方程式に対しても高い精度と効率を実現すること。
• 収束速度を制御するパラメータを導入し、精度を調整可能にする方法を提示すること。
• GALFsを基底関数として用いるペトロフ・ガレルキン定式化により、良好に条件付けられ、スパースな線形方程式系を達成すること。

提案手法

• 半直線上の分数階スツルム＝リウヴィル問題の固有関数として得られる一般化関連ラゲル関数（GALFs）を、試験関数および試行関数として用いる。
• 分数階積分による部分積分を用いて、試行関数から試験関数へ分数階微分を柔軟に移転できるペトロフ・ガレルキン変分法を定式化する。
• ラゲル多項式の再帰関係を活用し、剛性行列および質量行列の要素に対して明示的な表現を導出する。数値積分を回避することで、積分誤差を低減し、アセンブリコストを削減する。
• 行列要素を対角行列と下三角テープリッツ行列の積に因数分解し、高速な行列・ベクトル積を可能にする。
• 行列の構造的性質を活用して、高速変換を用いてO(N log N)の計算量で線形方程式系を解く。
• ガウス＝レジェンドル求積法を用いて、分布型微分方程式を多項式FIVPsに近似し、台形則よりも収束性を向上させる。

実験結果

リサーチクエスチョン

• RQ1 導関数次数≤1の半直線上の多項式分数初期値問題に対して、高い精度と計算効率を両立するペトロフ・ガレルキンスペクトル法を構築可能か？
• RQ2 基底関数をどのように選べば、解の特異性 t^α を自然に捉えられ、スペクトル収束を達成できるか？
• RQ3 収束速度を制御するパラメータを導入することで、この手法を調整可能な精度にできるか？
• RQ4 得られる線形方程式系のアセンブリおよび解法の計算量は何か？行列の構造を活用することで、O(N²)未満に削減可能か？
• RQ5 分布型微分方程式を近似する際、台形則と比較してガウス＝レジェンドル求積法の選択が全体の精度に与える影響は何か？

主な発見

• 重み付きL²ノルムにおける半直線上のスペクトル収束が達成され、収束速度は調整パラメータα₁に敏感である。
• 剛性行列および質量行列は、対角行列と下三角テープリッツ行列の積として構造化されており、O(N log N)の解法時間の実現が可能である。
• 数値積分を回避することで、明示的な行列要素が導出され、統合誤差が低減され、アセンブリコストが削減される。
• 数値実験の結果、ガウス＝レジェンドル求積法はK=10の求積点で機械精度に達するが、台形則はK=50であってもO(10⁻³)の誤差に飽和する。
• t^λ特異性を示す滑らかでない解に対しては、調整パラメータα₁を特異性指数に一致させる（例：λ = 2 + 1/3の場合はα₁ = 1/3）ことで、高い精度を維持する。
• 分布型微分方程式は、多項式FIVPsに近似することで効果的に解かれるが、Kが十分に大きくない場合、全体の誤差は求積誤差が支配的になる。