[論文レビュー] Estimation under uncertainties of acoustic and electromagnetic fields from noisy observations

Youri Shestopalov, Yuri Podlipenko+1 more|ArXiv.org|Oct 13, 2009

Numerical methods in inverse problems被引用数 8

ひとこと要約

本稿では、ノイズのある観測下で無限領域に定義されたヘルムホルツ方程式に従う音響場および電磁場のミニマックス推定手法を開発する。問題は有界領域における積分微分方程式または表面積分方程式に還元され、ノイズや不確実性が poorly 特徴付けられていても、解および右辺の関数型の推定がロバストに行える。明示的なミニマックス推定式と誤差上限が提供されている。

ABSTRACT

The creation and justification of the methods for minimax estimation of parameters of the external boundary value problems for the Helmholtz equation in unbounded domains are considered. When observations are distributed in subdomains, the determination of minimax estimates is reduced to the solution of integro-differential equations in bounded domains. When observations are distributed on a system of surfaces the problem is reduced to solving integral equations on an unclosed bounded surface which is a union of the boundary of the domain and this system of surfaces. Minimax estimation of the solutions to the boundary value problems from point observations is also studied.

研究の動機と目的

• 無限領域におけるヘルムホルツ方程式の外部境界値問題（BVP）の解およびパラメータの、不確実な観測下でのロバスト推定手法の開発を目的とする。
• 観測が未知のノイズによって汚されている状況で、解の線形関数型を推定する課題に取り組む。特に測定演算子が非可逆である場合を想定する。
• 最悪のノイズ実現の下でも最適な性能を保証するミニマックス推定の厳密な数学的枠組みを提供することを目的とする。
• 無限次元の推定問題を、有界領域における可解な積分微分方程式または表面積分方程式の系に還元することを目的とする。
• ミニマックス推定およびそれに付随する最大推定誤差の明示的表現を導出し、理論的・計算的妥当性を保証することを目的とする。

提案手法

• ヒルベルト空間における双対性および変分原理を用いて、ミニマックス推定問題を最適制御問題に還元する。
• 部分領域における分布観測からの関数型推定を、有界領域における積分微分方程式の解法に変換する。
• 表面に基づく観測（例：アンテナ）の場合は、領域境界および観測表面の和集合上での表面積分方程式の系に問題を還元する。
• 単層（S）、二重層（D）、特異的（H）ポテンシャルおよびその随伴作用素を用いて解を表現する。
• 基本解法およびグリーン関数を用いて、ミニマックス推定の明示的表現を構築する。
• 推定係数および観測点における関数値を決定するために、作用素 $I - K'_{-ar{k}} + i\eta T_{-ar{k}}$、$I - K_k - i\eta T_k$ およびその随伴作用素の連立系を解く。

実験結果

リサーチクエスチョン

• RQ1 観測が未知のノイズによって汚されている場合、ヘルムホルツ方程式の解の線形関数型のミニマックス推定をどのように構築できるか？
• RQ2 観測が部分領域に分布している場合と表面に分布している場合の、ミニマックス推定問題の数学的構造はそれぞれどのように異なるか？
• RQ3 ノイズおよびモデル不確実性の有界集合上での最悪のシナリオにおいて、推定誤差を最小化するにはどうすればよいか？
• RQ4 点観測または表面観測の場合、ミニマックス推定を計算するために解かなければならない積分方程式系は何か？
• RQ5 任意の波数 $k$ に対して、$ \operatorname{Im}k \geq 0$ を満たす場合、得られる方程式系が一意に解かれるための条件は何か？

主な発見

• 線形関数型 $l(\varphi)$ のミニマックス推定は、$\widehat{\widehat{l(\varphi)}} = \sum_{k=1}^{N} \bar{\hat{u}}_k y_k + \hat{c}$ で与えられ、係数 $\hat{u}_k$ および定数 $\hat{c}$ の明示的表現が得られている。
• 部分領域における分布観測の場合は、有界領域における積分微分方程式系の解法からミニマックス推定が導かれる。
• 表面観測の場合は、境界作用素 $S_k$、$K_k$、$T_k$ およびその随伴作用素を含む、表面積分方程式の連立系に還元される。
• 方程式系 (3.28)–(3.30) は、すべての $k$ に対して $ \operatorname{Re}k > 0$ かつ $ \operatorname{Im}k \geq 0$ を満たす場合に一意に解かれるため、推定の存在および一意性が保証される。
• 関数型の推定は、$\widehat{\widehat{l(\varphi)}} = \sum_{i=1}^{m} \bar{a}_i \hat{\varphi}(x_i)$ として表現され、ここで $\hat{\varphi}$ は境界積分ポテンシャルを用いて再構成される。
• 解は、観測点における値 $\hat{\varphi}(x'_l)$ を作用素 $\mathcal{W}_k$ および $\mathcal{V}_k$ を通じて決定する方程式系 (3.34)–(3.36) を解くことで得られる。