[論文レビュー] Recovery of a potential in a fractional diffusion equation

William Rundell, Masahiro Yamamoto|arXiv (Cornell University)|Nov 14, 2018

Fractional Differential Equations Solutions参考文献 11被引用数 8

ひとこと要約

本稿では、時間分数拡散方程式におけるポテンシャル $ q(x) $ を側面境界フラックスデータを用いて回復する問題について、一意性と安定性の結果を確立する。解析接続と線形化を用いて問題をスペクトル逆Sturm-Liouville問題に変換することで、逆問題が著しく悪条件であることが示され、$ \alpha $ が小さいほど悪条件性が増大する。一方、分数拡散では、Mittag-Leffler関数の遅い減衰のため、古典的拡散に比べて長時間における回復が良好であることが判明する。

ABSTRACT

We consider the determination of an unknown potential $q(x)$ form a fractional diffusion equation subject to overposed lateral boundary data. We show that this data allows recovery of two spectral sequences for the associated inverse Sturm-Liouville problem and these are sufficient to apply standard uniqueness results for this case. We also look at reconstruction methods and in particular examine the issue of stability of the solution with respect to the data. The outcome shows the inverse problem to be severely ill-conditioned and we consider the differences between the cases of fractional and of classical diffusion.

研究の動機と目的

• 時間分数拡散方程式におけるポテンシャル $ q(x) $ が、側面境界フラックスデータから一意に回復可能かどうかを確立すること。
• 逆問題の条件数を分析し、古典的放物型ケース（$ \alpha = 1 $）と比較すること。
• 逆係数問題において、分数拡散モデルが古典的拡散に比べて安定性や回復特性を向上させるかどうかを調査すること。
• スペクトル系列回復法と線形化写像逆問題に基づく再構成手法の開発と比較すること。
• サンプリングレートと測定タイミングが逆問題の条件数に与える影響を定量化すること。

提案手法

• 非同次問題を $ v(x,t) = u(x,t) - a(t) $ を用いて同次化し、源項 $ f(x,t) $ を含む修正された分数階PDEに変換する。
• $ t > T $ におけるフラックスデータ $ b(t) $ に対して解析接続を適用し、Sturm-Liouville作用素 $ \mathcal{L}u = -u_{xx} + qu $ に関連するスペクトル系列 $ \{ \lambda_j \} $ と $ \{ \gamma_j \} $ を回復する。
• 回復されたスペクトルデータを用いて、逆Sturm-Liouville問題における既知の一意性結果により $ q(x) $ を再構成する。
• $ q(x) $ からフラックスデータ $ b(t) $ への前向き写像を線形化し、特異値分解を用いて得られる線形逆問題を分析する。
• Mittag-Leffler関数 $ E_{\alpha,\beta}(z) $ を用いて解の時間領域挙動を特徴付け、悪条件性に影響する減衰特性を分析する。
• $ \alpha \in (0,1] $ の異なる値における線形化写像の特異値の減衰を比較し、$ \alpha = 1 $ では指数的減衰、$ \alpha < 1 $ では多項式的減衰が見られることを示す。

実験結果

リサーチクエスチョン

• RQ1 初期データがゼロで、$ u_x(0,t) = a(t) $ であると仮定したもとで、時間分数拡散方程式におけるポテンシャル $ q(x) $ が、境界フラックスデータ $ u_x(1,t) = b(t) $ から一意に回復可能かどうか。
• RQ2 逆問題の悪条件性は分数階数 $ \alpha $ にどのように依存するか。また、古典的放物型ケース（$ \alpha = 1 $）と比較するとどうなるか。
• RQ3 $ \alpha < 1 $ の場合、Mittag-Leffler関数の減衰が遅いため、長時間における解の安定性や再構成可能性が、$ \alpha = 1 $ の指数的減衰に比べて向上するか。
• RQ4 サンプリングレートや測定ウィンドウ（例：即時測定 vs. 遅延測定）が逆問題の条件数に果たす役割は何か。
• RQ5 スペクトル系列回復法や線形化写像逆問題を用いて、逆問題が効果的に再構成可能か。また、実用的な制限は何か。

主な発見

• 側面フラックスデータを用いた $ q(x) $ の逆問題は、関連するSturm-Liouville作用素の2つのスペクトル系列の回復が可能であるため、一意に解ける。
• すべての $ \alpha \in (0,1] $ において、逆問題は著しく悪条件である。$ \alpha $ が小さくなるほど、Mittag-Leffler関数の小引数における急速な減衰により、悪条件性が増大する。
• $ \alpha = 1 $ の場合、特異値は指数的に減衰するため、ノイズや測定誤差に対して非常に敏感であり、特にサンプリングが制限されている場合には回復が困難である。
• $ \alpha < 1 $ の場合、負の大きな引数におけるMittag-Leffler関数の減衰は多項式的であるため、長時間における特異値が相対的に大きく保たれ、古典的ケースよりも長時間測定が有効である。
• 測定が遅延している場合（例：$[1.5T, 3T]$）、$ \alpha = 1 $ の場合、再構成可能な固有モードの数は著しく減少するが、$ \alpha < 1 $ の場合、依然として有効である。$ \alpha = 3/4 $ の場合、最初の3つの特異値が $ 10^{-7} $ を上回る。
• $ \alpha < 1 $ では長時間における安定性が向上するが、特異値は依然として極めて小さい（例：$ \alpha = 3/4 $ の場合 $ \sim 10^{-7.5} $）ため、正則化なしでは実用的な再構成は不可能である。