[論文レビュー] Exponential sum approximations and fast evaluation of fractional integrals

William McLean|arXiv (Cornell University)|Jun 1, 2016

Fractional Differential Equations Solutions被引用数 1

ひとこと要約

本稿では、$ t^{-\beta} $ をコン pact 区間上で指数和を用いて近似するための新規手法を提示する。Beylkin と Monzon の手法とは異なり、代替積分表現を用いることで、Prony の手法による項の削減前における精度を著しく向上させる。その後、Prony の手法を用いて項数を削減しても精度の損失は最小限に抑えられる。

ABSTRACT

Given $\beta>0$ and $\delta>0$, the function $t^{-\beta}$ may be approximated for $t$ in a compact interval $[\delta,T]$ by a sum of terms of the form $we^{-at}$, with parameters $w>0$ and $a>0$. One such an approximation, studied by Beylkin and Monzon, is obtained by applying the trapezoidal rule to an integral representation of $t^{-\beta}$, after which Prony's method is applied to reduce the number of terms in the sum with essentially no loss of accuracy. We review this method, and then describe a similar approach based on an alternative integral representation. The main difference is that the new approach achieves much better results before the application of Prony's method; after applying Prony's method the performance of both is much the same.

研究の動機と目的

• 区間 $[\delta, T]$ における $ t^{-\beta} $ の指数和近似を、より高精度な初期近似として開発すること。
• 台形則に基づく既存手法と比較して、Prony の手法による項の削減前の近似品質を向上させること。
• Prony の手法による項の削減後も高い精度を維持すること。
• 最適化された指数和表現を用いて、分数階積分の効率的かつ高精度な評価を可能とすること。

提案手法

• 本手法は、Beylkin と Monzon の定式化とは異なる、$ t^{-\beta} $ の代替積分表現を用いる。
• この新たな積分表現に対して台形則を適用し、初期の指数和近似を生成する。
• その後、Prony の手法を用いて和の項数を削減しつつ精度を保持する。
• 得られた指数和は $ w>0 $、$ a>0 $ を満たし、最小限の誤差で $ t^{-\beta} $ を近似する。
• Prony の手法による項の削減前の精度を、Beylkin と Monzon の手法と比較する。
• 本手法は、Prony の手法による項の削減前において、優れた初期近似品質を達成する。

実験結果

リサーチクエスチョン

• RQ1 代替積分表現が初期指数和近似の精度に与える影響は何か？
• RQ2 異なる積分表現は、Beylkin と Monzon の手法よりも優れたPronyの手法による項の削減前の近似品質をもたらすか？
• RQ3 Prony の手法による項の削減後に、指数和近似の性能がどの程度低下するか？
• RQ4 Prony の手法を適用する前後における、新規手法と既存手法の近似誤差はどのように比較されるか？
• RQ5 指数和に基づく分数階積分評価において、初期近似の精度と計算効率のトレードオフは何か？

主な発見

• Prony の手法を適用する前段階で、本手法は Beylkin と Monzon の手法よりも著しく高い近似精度を達成する。
• Prony の手法を適用後、両手法とも同等の精度を示すため、項の削減段階における効果は同等である。
• 代替積分表現により、初期指数和の精度が向上し、圧縮前の誤差が低減される。
• 項数を削減しても高い精度を維持できることから、分数階積分の効率的評価が可能になる。
• Prony の手法による前段階の性能向上は、積分表現の選択が初期近似精度に大きく影響することを示唆する。
• 本手法は、指数和近似を用いた分数階積分の高速かつ高精度な評価の実用的フレームワークを提供する。