[論文レビュー] Hermite functions with discontinuous coefficients for the solution of fractal diffusion retrospective problems

Oleg Yaremko|arXiv (Cornell University)|Sep 18, 2013

Mathematical and Theoretical Analysis参考文献 3被引用数 3

ひとこと要約

本稿では、空間的に不連続な拡散係数を考慮するため、分岐的均質媒質におけるフラクタル拡散の逆問題を解くための新規な手法を提案する。古典的ヘルミート多項式を空間的不連続な拡散係数を反映させるように一般化することで、分数階時間微分と作用素変換に基づく級数解が導出され、逆問題の解の存在および一意性が証明される。

ABSTRACT

In this article we study the retrospective inverse problem. The retrospective inverse problem consists of in the reconstruction of a priori unknown initial condition of the dynamic system from its known final condition. Existence and uniqueness of the solution is proved.

研究の動機と目的

• フラクタル拡散における逆問題の不適定性に対処し、最終状態の観測から初期状態を再構築することを目的とする。
• 古典的ヘルミート関数法を、分岐的均質媒質をモデル化する不連続拡散係数を有する系へと拡張することを目的とする。
• 作用素変換と一般化されたヘルミート多項式を用いて、逆問題の形式的解を確立することを目的とする。
• 基礎となる境界値問題が無限に可解である条件の下で、解の存在および一意性を証明することを目的とする。
• 分数階時間微分（α）と双曲的極限（α=2）を含む解の枠組みを一般化することを目的とする。

提案手法

• 逆問題を、分岐的均質領域におけるスツルム＝リウヴィル問題の固有関数を用いてスペクトル領域に写像する作用素変換法を用いる。
• スペクトル領域におけるミットャグ＝レフラー関数のテイラー展開から導かれる不連続係数を有するヘルミート関数を導入する。
• 逆作用素変換を適用して、t=0における最終状態データの微分を含む級数として初期状態を再構築する。
• 一般化されたヘルミート多項式の形を導出する：$ H_j(x) = \Sigma_{k=0}^{[j/2]} \frac{(-1)^k j!}{\Gamma(k\alpha+1)(j-2k)!} x^{j-2k} $、ここでαはフラクタル挙動を支配する。
• 解の表現を確立する：$ f(x) = \sum_{j=0}^{\infty} \frac{u_j(\tau) \tau^{\beta j}}{j!} H_j\left(\frac{x}{\tau^\beta}\right) $、ここで$ \beta = \alpha/2 $であり、フラクタルおよび双曲的ケースに有効である。
• 双曲的ケース（α=2）との類似性を通じて手法の妥当性を検証し、既知の解$ f(x) = \frac{u(\tau, \tau+x) + u(\tau, \tau-x)}{2} $が回復されることを示す。

実験結果

リサーチクエスチョン

• RQ1 フラクタル拡散問題において、不連続な拡散係数を扱えるようにヘルミート関数をどのように一般化できるか？
• RQ2 分岐的均質媒質における逆問題の解の存在および一意性を保証する条件は何か？
• RQ3 分数階時間微分（αを介して）の導入が、解の構造および収束性にどのように影響を与えるか？
• RQ4 作用素変換法を不連続係数を有する逆問題に適応し、形式的解を導出できるか？
• RQ5 一般化されたヘルミート関数と古典的ヘルミート多項式との間の関係は、α→2の極限でどのように表れるか？

主な発見

• 逆問題の解は、不連続係数を有する一般化されたヘルミート関数を含む無限級数として表現される。
• 基礎となる境界値問題が無限に可解である条件の下で、解の存在および一意性が厳密に証明される。
• 双曲的ケース（α=2）では、解は$ f(x) = \frac{u(\tau, \tau+x) + u(\tau, \tau-x)}{2} $に簡略化され、既知の結果が回復される。
• 一般化されたヘルミート多項式は$ H_j(x) = \sum_{k=0}^{[j/2]} \frac{(-1)^k j!}{\Gamma(k\alpha+1)(j-2k)!} x^{j-2k} $として定義され、α=2の極限で古典的ヘルミート多項式に還元される。
• 標準的ヘルミート関数に代えて不連続係数に適合した関数を用いることで、古典的逆問題解法が成功裏に一般化される。
• 半平面における逆ディリクレ問題は、類似の枠組みを用いて解かれる。その結果、$ f(y) = \sum_{j=0}^{\infty} \frac{u^{(j)}(l,0)}{j!} \cdot \frac{(y+li)^j + (y-li)^j}{2} $が得られ、解析接続の場合は$ f(y) = \text{Re}\, u(l, y+li) $となる。