[論文レビュー] Stability result for elliptic inverse periodic coefficient problem by partial Dirichlet-to-Neumann map

Mourad Choulli, Yavar Kian+1 more|arXiv (Cornell University)|Jan 20, 2016

Numerical methods in inverse problems参考文献 26被引用数 7

ひとこと要約

この論文は、無限にのびる円柱型波導における1周期的ポテンシャルを、境界のわずかに半分以上に相当する部分で測定された部分的ディリクレ・トゥ・ノイマン写像から回復する逆問題に対して、対数対数安定性推定を確立する。フーリエ解析と複素幾何学的オプティクス解を用いて、測定誤差に二重対数的依存性を示す $ H^{-1} $-ノルムにおけるポテンシャルの安定回復を証明している。これは、ノイマンデータが境界の部分集合に制限されている場合でも成り立つ。

ABSTRACT

We study the inverse problem of identifying a periodic potential perturbation of the Dirichlet Laplacian acting in an infinite cylindrical domain, whose cross section is assumed to be bounded. We prove log-log stable determination of the potential with respect to the partial Dirichlet-to-Neumann map, where the Neumann data is taken on slightly more than half of the boundary of the domain.

研究の動機と目的

• 境界の一部からの測定値から、断面が有界な波導における1周期的ポテンシャルを特定する逆問題に取り組む。
• 無限領域、特に無限円柱型波導における逆係数問題の安定性推定を確立する。
• 有界領域から無限領域への部分的データ安定性結果の拡張、特に周期的係数の場合を対象とする。
• 境界の一部（全境界のわずかに半分以上）で測定されたノイマンデータから、ポテンシャルが安定に回復可能であることを証明する。
• 測定誤差に依存する再構成誤差の定量的評価を行う対数対数型安定性推定を導出する。

提案手法

• 周期的方向 $ x_1 $ におけるフーリエ分解を用い、3次元問題を $ k \in \mathbb{Z} $ をパrameterとする2次元問題の族に還元する。
• 波導の構造と周期的構造に適合した、シュレーディンガー方程式 $ (-\Delta + V)v = 0 $ の複素幾何学的オプティクス（CGO）解を構成する。
• 境界の部分集合 $ G \subset \partial\Omega $ で定義される部分的ディリクレ・トゥ・ノイマン写像 $ \Lambda_V $ を用いる。ここで $ G $ は境界の $ \xi $-影法線面の近傍である。
• カルレマン推定とエネルギー推定を用いて、二つのポテンシャル $ V_1 $ と $ V_2 $ に対応する解の差に対する $ L^2 $-推定を導出する。
• 双対性の議論を適用し、ポテンシャル差の $ H^{-1} $-ノルムを測定領域における解の差の $ L^2 $-ノルムに関連付ける。
• CGOアンザッツのパrameter（周波数 $ \rho $、スケーリング $ r $）を最適化することで、指数関数的増大と減衰の項のバランスを取る。これにより、二重対数的安定性推定が得られる。

実験結果

リサーチクエスチョン

• RQ1 無限波導において、境界のわずかに半分以上に相当する部分で測定された部分的ディリクレ・トゥ・ノイマンデータから、$ V \in L^\infty(\Omega) $ の周期的ポテンシャルを安定に回復できるか？
• RQ2 ノイマンデータが境界の部分集合に制限されている場合に、周期的ポテンシャルを特定する逆問題の最良の安定性推定は何か？
• RQ3 境界データが不完全で領域が無限である場合、安定性は測定誤差にどのように依存するか？
• RQ4 有界領域からの結果を拡張して、無限円柱領域における周期的係数に対して対数対数安定性推定を確立できるか？
• RQ5 照明領域と影領域の幾何的構造は、逆問題の安定性にどのような役割を果たすか？

主な発見

• 論文は、対数対数型安定性推定を確立した：$ \|V\|_{H^{-1}(\check{\Omega})} \leq C (\ln |\ln \gamma|)^{-1} $、ここで $ \gamma = \|\Lambda_{V_2,\theta} - \Lambda_{V_1,\theta}\| $ であり、測定誤差が小さい場合に成り立つ。
• 安定性推定は、ノイマンデータが $ \xi $-影法線面の近傍である集合 $ G $ で測定されているという仮定のもとで成立する。この $ G $ は全境界のわずかに半分以上を占める。
• 波導構造と周期性のおかげで、無限領域における解の振る舞いが有効に制御できるため、領域の非有界性に対しても結果はロバストである。
• 主な技術的革新は、測定領域における $ L^2 $-ノルムが制御されたCGO解の構成と、指数的重みにおける最適パラメータ選択にある。
• 最終的な安定性推定は、カルレマン推定における指数的項のバランスを最適化し、周波数パラメータ $ \rho $ について最適化することで得られ、測定誤差に二重対数的依存性が現れる。
• 安定性推定における定数 $ C $ は、断面 $ \omega $ の幾何、ポテンシャルの $ L^\infty $-有界値 $ M_+ $、測定領域 $ G' $ のみに依存し、ポテンシャル自体には依存しない。