[論文レビュー] Modified Tikhonov regularization for identifying several sources

Ole Løseth Elvetun, Bjørn Fredrik Nielsen|arXiv (Cornell University)|Nov 9, 2020

Numerical methods in inverse problems参考文献 15被引用数 4

ひとこと要約

本稿では、楕円型PDEにおけるディリクレ境界データから複数の局所的源を同定するための修正されたチホノフ正則化フレームワークを提案する。問題固有の正則化作用素Wを用いることで、源の区別可能性が向上する。この手法により、単一源の正確な回復が可能となり、三段階のアルゴリズムを用いることで複数源の局所化が改善される。理論的および数値的証拠から、同定可能性はドメインの幾何構造および境界からの源の位置に強く依存することが示されている。

ABSTRACT

We study whether a modified version of Tikhonov regularization can be used to identify several local sources from Dirichlet boundary data for a prototypical elliptic PDE. This paper extends the results presented in [5]. It turns out that the possibility of distinguishing between two, or more, sources depends on the smoothing properties of a second or fourth order PDE. Consequently, the geometry of the involved domain, as well as the position of the sources relative to the boundary of this domain, determines the identifiability. We also present a uniqueness result for the identification of a single local source. This result is derived in terms of an abstract operator framework and is therefore not only applicable to the model problem studied in this paper. Our schemes yield quadratic optimization problems and can thus be solved with standard software tools. In addition to a theoretical investigation, this paper also contains several numerical experiments.

研究の動機と目的

• 単一源同定に関する先行研究を、楕円型PDE逆問題における複数の局所的源への拡張すること。
• 複数の源が境界測定値から区別可能であるための理論的条件を確立すること、特にPDEの次数とドメインの幾何構造に基づくものである。
• チホノフ正則化と半径最適化による後処理を組み合わせることで、源局所化を向上させる数値アルゴリズムを開発すること。
• 特定のPDEモデルに限定されない抽象的作用素フレームワークに適用可能な、単一源同定の一意性結果を証明すること。
• 数値実験を通じて、源回復性能がドメイン形状および境界からの源の位置に著しく依存することを示すこと。

提案手法

• 修正されたチホノフ正則化アプローチを二次計画問題として定式化する：L²(∂Ω)におけるデータ不適合項と、源fに作用する線形作用素Wを含む正則化項の和を最小化する。
• 正則化作用素Wは、各基底関数ϕiに対してWϕi = ||Pϕi||X ϕiを満たすように定義される。ここでPは、前方作用素の随伴作用素の像への正射影である。
• 源を表すためにL²(Ω)の有限次元部分空間Fhを用い、局所的サポートを持つ基底関数を用いることで空間分解能を向上させる。
• 三段階のアルゴリズムを提案する：第一に、修正チホノフ問題を解いて初期推定値を得る。第二に、解の局所的最大値を特定する。第三に、これらの最大値を中心とする円形源の半径について最適化を行うことで、源の幾何形状を精緻化する。
• ε=0（ポアソン方程式）の場合、前方問題の適切性を保つために、源の平均を引く処理が施される。
• 前方作用素K_hは、源fを∂Ωにおける解uのディリクレトレースに写像する。解は境界上で平均がゼロになるように正規化され、定数解による一意性の欠如に対処する。

実験結果

リサーチクエスチョン

• RQ1 修正チホノフ正則化スキームは、ディリクレ境界データに基づいて複数の局所的源を区別できるか？
• RQ2 PDEの次数（第二次と第四次）は、複数源の同定可能性にどのように影響するか？
• RQ3 ドメインの幾何構造および源が境界に対して相対的にどこに位置するかが、源回復性能にどの程度影響を及えるか？
• RQ4 提案された三段階のアルゴリズムは、標準的チホノフ正則化と比較して源局所化を改善できるか？
• RQ5 最小ノルム最小二乗解が逆源問題における真の源と等価となる条件は何か？

主な発見

• 提案されたW作用素を用いた修正チホノフ正則化により、抽象的作用素フレームワークを用いた理論的証明により、単一局所源の正確な回復が可能である。
• 複数源を区別できる能力は、PDEの平滑化性に強く依存し、第四次問題の方が第二次問題よりも優れた源分解能を示す。
• 数値実験の結果、凸領域よりも非凸領域が、同程度のサイズの領域において、特にコーナーやエッジ付近に源がある場合、より良好な源回復をもたらすことが示された。
• 三段階のアルゴリズム（初期チホノフ解法、最大値の局所化、半径最適化）は、特に小さなまたは弱い源に対して、源再構成を顕著に改善する。
• ポアソン方程式（ε=0）の場合、源の平均を引いた後も、この手法は有効であり、最小ノルム解は自動的に∫Ω f dx = 0を満たす。
• Method Iと半径最適化法との間で一貫した局所化が得られることは、信頼性の重要な指標である。これより、初期源局所化にはMethod Iを最初に用いるべきである。