[論文レビュー] Probability distributions generated by fractional diffusion equations

Francesco Mainardi, Paolo Paradisi+1 more|ArXiv.org|Apr 3, 2007

Fractional Differential Equations Solutions参考文献 46被引用数 149

ひとこと要約

本稿は、時間的または空間的分数階微分を用いて標準的拡散方程式を一般化する分数階拡散方程式が、安定分布族に属する確率密度関数（pdf）を生成することを示している。特殊関数としての Wright 関数やMittag-Leffler 関数を用いた分数階微分の理論により、これらの方程式の基本解が、パワー則的または指数的尾部を示す対称的かつ安定なpdfを生成することを示しており、古典的ブラウン運動モデルを、経済物理学における重尾的・非ガウス的挙動を含むように拡張している。

ABSTRACT

Fractional calculus allows one to generalize the linear, one-dimensional, diffusion equation by replacing either the first time derivative or the second space derivative by a derivative of fractional order. The fundamental solutions of these equations provide probability density functions, evolving on time or variable in space, which are related to the class of stable distributions. This property is a noteworthy generalization of what happens for the standard diffusion equation and can be relevant in treating financial and economical problems where the stable probability distributions play a key role.

研究の動機と目的

• 分数階微分を用いた分数階微分法を用いて、非ガウス的確率過程をモデル化する。
• 時間的および空間的分数階拡散方程式の基本解が、安定確率分布を生成することを確立する。
• Wright 関数やMittag-Leffler 関数などの特殊関数を通じて、分数階拡散過程と安定法則を結びつける。
• 重尾的金融・経済データを分数階拡散を用いてモデル化する理論的基盤を提供する。
• ブラウン運動およびLévy飛行の古典的結果を、分数階動力学に拡張する。

提案手法

• Riemann-LiouvilleおよびCaputoによる分数階微分の定義を用いて、時間的分数階拡散方程式を定式化する。
• Fourier変換を適用して初期値問題を解き、基本解をWright関数を用いて表現する。
• ラプラス変換を用いて信号伝達問題を解き、時間における片側安定pdfを導出する。
• 安定分布の理論と、Fox H関数およびMellin-Barnes積分を用いた表現に依存する。
• 特殊関数の級数および積分表現を用いて、解の漸近的挙動を導出する。
• 積分恒等式およびラプラス変換を通じて、Wright関数と安定pdfとの関係を確立する。

実験結果

リサーチクエスチョン

• RQ1 拡散方程式における分数階微分が、得られる確率密度関数に与える影響は何か？
• RQ2 時間的分数階および空間的分数階拡散方程式の基本解を、安定確率分布として解釈できるか？
• RQ3 Wright関数およびMittag-Leffler関数は、これらの一般化された拡散過程を特徴付ける上で果たす役割は何か？
• RQ4 解の漸近的挙動（例：パワー則的または指数的減衰）が、基礎となる分布の安定性指数αとどのように関係するか？
• RQ5 これらの分数階拡散モデルは、確率過程における古典的ガウス分布およびLévy分布をどのように一般化するか？

主な発見

• 時間的分数階拡散方程式の初期値問題の基本解は、空間において対称的かつ安定な確率密度関数を生成し、すべてのモーメントが有限である。時間経過に伴い分散が t^β に比例して変化する。ここで β は分数階微分の次数である。
• 時間的分数階拡散方程式の信号伝達問題の解は、時間における片側安定pdfを生成し、α = 1/2 のLévy法に一致する。t^{-3/2} の代数的減衰を示す。
• 基本解は、Wright関数 M(z; ν) を用いて表現され、α = 2 の場合にガウス分布、α = 1/2 の場合にLévy分布を一般化する。
• α = 1/3 の場合、解はAiry関数および修正ベッセル関数 K_{1/3} に関連する安定pdfを生成し、安定分布理論における既知の解析的形と一致する。
• 片側安定pdf Φ₁(y) のラプラス変換は exp(−s^α) であり、Φ₂(y) のラプラス変換は (1/α)E_{1/α}(−s) である。これにより、解がMittag-Leffler関数と関連づけられる。
• 解の漸近的挙動は、Φ₁(y) ではゼロ近傍で指数的減衰を示し、Φ₂(y) では無限大で重尾的減衰を示す。これは安定分布の性質と整合的である。