[論文レビュー] Superconvergence Points For The Spectral Interpolation Of Riesz Fractional Derivatives

Beichuan Deng, Zhimin Zhang+1 more|arXiv (Cornell University)|Sep 29, 2017

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

ひとこと要約

本稿では、Lobatto型多項式と一般化ヤコビ関数（GJF）を用いたリーマン・フーリエ微分のスペクトル補間において、超収束点を同定し、$ \alpha \in (0,1) $ の場合に $ O(N^{-2}) $、$ \alpha > 1 $ の場合に $ O(N^{-(\alpha+3)/2}) $ の収束速度の向上を示している。超収束点はGJFを用いる場合、補間ノードと一致し、分数階微分方程式（FDE）の解法において高次精度を実現できる。

ABSTRACT

In this paper, superconvergence points are located for the approximation of the Riesz derivative of order $α$ using classical Lobatto-type polynomials when $α\in (0,1)$ and generalized Jacobi functions (GJF) for arbitrary $α> 0$, respectively. For the former, superconvergence points are zeros of the Riesz fractional derivative of the leading term in the truncated Legendre-Lobatto expansion. It is observed that the convergence rate for different $α$ at the superconvergence points is at least $O(N^{-2})$ better than the optimal global convergence rate. Furthermore, the interpolation is generalized to the Riesz derivative of order $α> 1$ with the help of GJF, which deal well with the singularities. The well-posedness, convergence and superconvergence properties are theoretically analyzed. The gain of the convergence rate at the superconvergence points is analyzed to be $O(N^{-(α+3)/2})$ for $α\in (0,1)$ and $O(N^{-2})$ for $α> 1$. Finally, we apply our findings in solving model FDEs and observe that the convergence rates are indeed much better at the predicted superconvergence points.

研究の動機と目的

• 古典的Lobatto多項式を用いて、$ \alpha \in (0,1) $ のリーマン・フーリエ微分のスペクトル補間のための超収束点を特定すること。
• $ \alpha > 1 $ への分析の拡張。強い特異性を扱うために、一般化ヤコビ関数（GJF）を用いて安定かつ高精度な補間を可能にする。
• 任意の $ \alpha > 0 $ に対して、GJFに基づく補間の適切な定式化、収束性、および超収束特性を理論的に分析すること。
• 超収束点の実用的有用性を、高精度な解を得るためのモデル分数階微分方程式（FDE）の解法において示すこと。
• 予測された点での超収束を達成するための補間法、Petrov-Galerkin法、およびスペクトルコロケーション法の性能を比較すること。

提案手法

• $ \alpha \in (0,1) $ の場合、超収束点は、ルジャンドル-Lobatto展開における主項のリーマン・フーリエ微分の零点として特定される。
• $ \alpha > 1 $ の場合、$ x = \pm 1 $ における特異性を扱うために一般化ヤコビ関数（GJF）を基底関数として用い、安定なスペクトル補間を可能にする。
• 補間は $ u_N(x) = \sum_{j=0}^N \hat{\ell}_j(x) v_j $ として構築され、ここで $ \hat{\ell}_j(x) = (1 - x^2)^{\alpha/2} \ell_j(x) $ であり、$ \ell_j $ はGJFノード上のラグランジュ基底関数である。
• 超収束点は、$ P^{\alpha/2, \alpha/2}_{N+1}(x) $ の根として定義され、これが補間ノードとして機能する。
• FDEをこれらの超収束点で強制することでスペクトルコロケーション法を定式化し、線形方程式系 $ (D + \Lambda)V = F $ を得る。ここで $ D $ は分数階微分行列である。
• GJF基底関数のリーマン・フーリエ微分の既知の公式により、分数階微分行列 $ D $ の解析的計算が可能になる。

実験結果

リサーチクエスチョン

• RQ1 $ \alpha \in (0,1) $ のスペクトル補間におけるリーマン・フーリエ微分の超収束点はどこに位置し、その点での収束速度の向上はどの程度か？
• RQ2 $ \alpha > 1 $ の場合、標準的なLobatto多項式が強い特異性のため失敗する中で、どのようにして超収束を達成できるか？
• RQ3 一般化ヤコビ関数（GJF）を用いた場合、$ \alpha > 1 $ における超収束点での理論的収束速度はどの程度か？
• RQ4 補間法、Petrov-Galerkin法、スペクトルコロケーション法の間で、超収束性能はどのように比較できるか？
• RQ5 予測された超収束点は、モデルFDEの数値解法の精度向上に実際に有効に利用できるか？

主な発見

• $ \alpha \in (0,1) $ の場合、超収束点はルジャンドル-Lobatto多項式のリーマン・フーリエ微分の零点であり、グローバル最適レート比で $ O(N^{-2}) $ の収束速度向上が達成される。
• $ \alpha > 1 $ の場合、一般化ヤコビ関数（GJF）の使用により安定な補間が可能となり、超収束点は補間ノードと一致する。具体的には、$ P^{\alpha/2, \alpha/2}_{N+1}(x) $ の根である。
• 超収束点における収束速度の向上は、$ \alpha \in (0,1) $ では $ O(N^{-(\alpha+3)/2}) $、$ \alpha > 1 $ では $ O(N^{-2}) $ であり、両者ともグローバルレートより顕著に優れている。
• スペクトルコロケーション法では、超収束レートはグローバル最適レート比で $ O(N^{-3}) $ から $ O(N^{-4}) $ の高速化を示し、補間法およびPetrov-Galerkin法を上回る性能を発揮する。
• 数値実験により、超収束点における誤差はグローバル最大誤差よりも著しく小さく、特にコロケーション法で顕著な改善が確認された。
• 理論は超収束点の位置を正確に予測できており、その利用により、GJFに基づく補間を用いた場合、特にモデルFDEの解法において顕著な精度向上が達成される。