逆問題

**逆問題**とは、ある系（物理現象や数学モデル）において、観測や結果（出力）からその原因や内部構造（入力・パラメータ）を推定・復元する問題のことを指す。

## 一般的な定義
***順問題**（じゅんもんだい、）（正問題）：入力や原因が与えられたときに、その結果や応答を計算・予測する問題。
例）地震波の伝播の計算から地表の振動を求める。
***逆問題**（ぎゃくもんだい、）：観測された結果や応答が与えられたときに、その原因や内部パラメータを推定する問題。
例）観測された地震波から地下の地質構造を推定する。

## 特徴
*不適定性（ill-posedness）：逆問題は、多くの場合「解が存在しない」「解が一意でない」「解がデータのわずかな誤差に非常に敏感」という性質を持つ（アダマールの条件に反する）。
*正則化が必要：不安定な解の振る舞いを抑えるために、追加の条件（正則化）を用いて解を安定化させる手法が使われる。
*応用範囲が広い：地球物理学、医療画像（CT、MRI）、材料評価、音響・電磁気学、金融工学など様々な分野に応用される。

逆問題とは、「観測された結果から原因を推定する問題」であり、しばしば不安定で不確定な解決が求められる。

## 歴史
逆関数の問題であると解釈すると、紀元前から扱われている問題である。しかし歴史的には物理において順問題と逆問題は今の使われ方とは異なっていた。例えばニュートンの時代では物体の動きからその作用する力を導くことが順問題だとされ、作用する力から物体の軌道を導くことが逆問題だとされていた。順問題と逆問題の定義は実際曖昧で、時代や学問分野によって異なることが多い。一般的には1820年代にニールス・アーベルがヤコビの逆問題を研究したのが、逆問題の最初の研究とされる。アーベルは方程式の解の公式の研究でも有名だが、方程式の解の公式自体も逆問題である。1929年にヴィクトル・アンバルツミャンも逆問題に関する論文を発表している。第二次世界大戦中に、弾道計算やレーダー探査など軍事上の目的により急速に発展した。現在では、非破壊検査や医療を目的とした利用も盛んに研究されている。

## 概要
順問題と逆問題は対になる概念であり、どちらが順でどちらが逆かというのは相対的な問題である。しかし対称的ではない。一般に、古くから問題として認識され研究が行われている方向のプロセスによるものを順問題とし、その逆方向のプロセスで解く方法は自明ではないのだが、それを解くことで何らかの工学的・その他の利用ができるような問題のことを逆問題と言う。

単純な順問題・逆問題の例を示す。*f*(*x*) = *x*2 という関数について考える。*f*(2) や f(3) を計算して 4 や 9 と求めるのが順問題である。逆問題は2通りある。1つ目は、*f*(*x*) = 25 という問題で、*x* = 5 と解く問題である。2つ目は、関数が未知で、*f*(1) = 1, *f*(2) = 4, *f*(3) = 9 という情報から、*f*(*x*) がいかなるものかを推測する問題である。

この例において、特にひとつめは逆関数 *f* -1(*x*) = √*x* によって容易に得られる。しかし、一般には逆関数が容易にはわからない関数も多く、そういった場合を特に扱うのがこの分野である。

逆問題は入力を求める、と一口に言っても、ここでの「入力」とは単に入力信号のようなものだけを指すのではない。例えば、物理学・工学で材料に関する問題においては、扱う材料に作用している外力を求める逆問題だけでなく、
*材料の境界・領域形状を求める
*材料を支配している方程式を求める
*材料についての境界値あるいは初期値を求める
*材料の物性値を求める
といった、複数の逆問題が存在する。様々な問題設定があるように、様々な有益な用途があり、理論・実用の両面から研究が行われている。

### 問題の種類
逆問題としては、以下の2つのパターンがある。
# 既知：モデル（関数）と出力未知：入力
# 既知：入力と出力未知：モデル（関数）

順問題は、入力とモデル（関数）が既知で、出力が未知である。

## 適切性と非適切な問題
逆問題を解く際によく問題になるのが**適切性** (良設定問題、) である。次の3つの条件が満たされるとき、**[[ジャック・アダマール|アダマール]]の意味で適切である**という。

#解の存在性： 解が存在すること
#解の一意性： 解がただ一つであること
#解の安定性： 入力に微小な変動を与えたときに、出力の変動も微小であること

上に挙げた *f*(1) = 1, *f*(2) = 4, *f*(3) = 9 から *f*(*x*) を推測する例で、逆問題の答えとしては f(x) = x^2 のほか、例えば f(x) = x^3 - 5 x^2 + 11x - 6 も解となり、解の一意性が満たされない。よって、非適切 (ill-posed) な問題といえる。

その他、微分方程式、積分方程式などに関する逆問題では解の安定性が得られず非適切な問題となることが多い。

### ティホノフの正則化法
非適切な問題の近似解を得る手法として最もよく使われるのが**ティホノフの正則化法** () である。

線形有界作用素 *K* : *X*→*Y* についての方程式 *Kx* = *y* の近似解を得るために、**ティホノフ[[汎関数]]**：
:J_\alpha(x) = \|Kx-y\|^2+\alpha\|x\|^2　　for *x*∈*X*　　*α*：正則化パラメータ
を導入し、これを最小にする *x*α∈*X* を求める。

近似解を真の解に近づけるためには、正則化パラメータ *α* を誤差 *η* = (*δ*, *h*) に応じて次のように設定すればよいといわれている：
:\mu_\eta(K_h, y_\delta) = \inf_{x \in D} \|K_h x-y_\delta\|_Y
という汎関数を設定し、
:\rho_\eta^\kappa(\alpha) = \|K_h x_\eta^\alpha-y_\delta\|_Y^2 - (\delta+h\|x_\eta^\alpha\|_X)^2 - \mu_\eta^\kappa(K_h, y_\delta)^2
について、*&rho;*&eta;&kappa; (*&alpha;**) = 0 となるような *α** (*η*) を選ぶ。

最適正則化パラメータの推定法として文献では、GCV法、L-カーブ法、Quasi-Optimal法を掲げている。

*関連項目： ヒルベルト空間、コンパクト

## 入力が未知な線形モデルでの正則化
ここでは、モデルと出力が既知で、入力が未知の問題を扱う。モデルは線形モデルである。*N* 個の誤差のある観測値（出力） *y**1*, *y**2*, \dotsc, *y**N* から、*M* 個のパラメタ（入力） *x**1*, *x**2*, \dotsc, *x**M* を推定するという問題を扱う。

観測不可能な真の値 *x**i* と、観測値 *y**&mu;* は、線形の関係がある（線形モデル）と仮定する。
: y_i = \sum_\mu K_{i\mu} x_\mu + n_i 
ここで、*K**i&mu;* は分かっているものとする。ノイズ *n**i* は観測不可能だが、その統計的性質として平均 0 と、共分散
: S_{ij} = \mathrm{E}(n_i n_j) 
は分かっているものとする。ここで、E() は統計平均を取る操作。

もし、観測が全て独立でその数 *N* が、パラメタの数 *M* より多ければ、最小自乗法で *x* の推定値を求めることができる。しかし、観測が独立でなかったりその数がパラメタの数より少ないとき、*x* を求める問題は劣決定となり、上記*適切性*のうち*解の一意性*が満たされない非適切な問題となる。よって、その問題に即した適当な正則化を行って、解を求める必要がある。

式で書けば、ノイズを最小にするには
:J = \sum^N_{i=1} \sum^N_{j=1} n_i S^{-1}_{ij} n_j 
あるいは、行列表示して（上付き添字 \intercal は転置行列を表す）
:
J = \boldsymbol{n}^\intercal \boldsymbol{S}^{-1}\boldsymbol{n}
= (\boldsymbol{K}\boldsymbol{x} - \boldsymbol{y})^\intercal \boldsymbol{S}^{-1}(\boldsymbol{K}\boldsymbol{x} - \boldsymbol{y})

なる *J* を最小にする ***x*** を決める問題になるが、行列 ***K*** は行より列が多く、***Kx***=***y*** の解が無数にあるという状況になる。そのため、正則化を行って解をひとつに定める。以下にいくつかの正則化の方法を紹介する。以下の議論で本質的に重要でないため、ノイズは分散 1 でそれぞれ無相関なものとする（つまり***S***は単位行列）。

### 零次の正則化
正則化パラメタ *&alpha;* を用いて、
: J = (\boldsymbol{K}\boldsymbol{x} - \boldsymbol{y})^\intercal (\boldsymbol{K}\boldsymbol{x}-\boldsymbol{y}) + \alpha \boldsymbol{x}^\intercal \boldsymbol{x}
と取る。つまり、無数の解のうち ***x*** の大きさを小さくにするものを推定値として採用する。*&alpha;* が小さいとき、これは ***Kx***=***y*** を特異値分解で解いた解と一致する。パラメタ *&alpha;* の取り方は問題設定によって異なる。一例としては、観測誤差が正規分布に近いと期待される場合、第一項は自由度 *N* のカイ二乗分布となることが期待され、その平均値は *N* となる。よって、第一項が *N* に近くなるように *&alpha;* を調整する。

### 線形の正則化
**x** の大きさより、滑らかさが重要なときは、x_{i+1}-x_i を第二項にした
: J = (\boldsymbol{K}\boldsymbol{x} - \boldsymbol{y})^\intercal (\boldsymbol{K}\boldsymbol{x}-\boldsymbol{y}) + \alpha (\boldsymbol{Bx})^\intercal \boldsymbol{Bx}
を最小にするような ***x*** を定める。

ここに、***B*** は
:
\begin{bmatrix}
-1 & 1 & 0 & \ldots & 0 \\
0  & -1 & 1 & \ldots & 0 \\
\vdots & & \ddots & & \vdots \\
0 & 0 & \ldots & -1 & 1
\end{bmatrix}

なる成分を持つ。

同様に ***x*** が線形に増加すると期待されるとき、***x*** が二次関数的に増加すると期待されるとき、なども適当な ***B*** を設定することで解くことができる。

### バッカス=ギルバート法
上の二つの正則化もそうであったが、***x*** の推定値 \hat{\boldsymbol{x}} は観測値の線型結合で表されている。

: \hat{\boldsymbol{x}} = \boldsymbol{L} \boldsymbol{y} 

観測値 **y** は、真の値をノイズ付きで観測したものだから、***y*** の定義式を代入して

: \hat{\boldsymbol{x}} = \boldsymbol{L}(\boldsymbol{K}\boldsymbol{x} + \boldsymbol{n}) 

***K*** は *N* 行 *M* 列（ N ）だから、逆行列は存在しないが、ノイズがなければ、よい観測は \hat{\boldsymbol{x}}=\boldsymbol{x} となるはずである。そこで、***LK***=***I*** と ***L*** を定めればよさそうである。すなわち、行列 ***LK*** の成分 {*LK*}*ij* がクロネッカーのデルタ *&delta;**ij* になれば理想的である。しかし、実際にはノイズがあるからこのようにはならない。そこで、クロネッカーのデルタにできるだけ形の近いものになるようにする。バッカスとギルバートは ***LK*** の行ベクトルのクロネッカーのデルタからのずれ、

: \sum_{j=1}^N (i-j)^2(\{LK\}_{ij}-\delta_{ij})^2 

を最小にすれば良いと考えた。これが最小化関数 *J* の第一項となる。

正規化のための第二項は、多数の観測で得られたパラメタ推定値のばらつきが少ないという条件を用いる。観測値のばらつき（ノイズ）は **n** だから、

: \boldsymbol{L} \mathrm{E}(\boldsymbol{n}\boldsymbol{n}^\intercal)\boldsymbol{L}^\intercal 

が最小化関数 *J* の第二項となる。

### 正則化項の意味
劣決定な逆問題では、与えられたデータ ***y*** だけでは拘束条件 ***Kx***=***y*** を満たす推定パラメタ ***x*** が一意に決まらないため、正則化項 ***R*** を第二項に加えた

: J = (\boldsymbol{K}\boldsymbol{x} - \boldsymbol{y})^\intercal (\boldsymbol{K}\boldsymbol{x}-\boldsymbol{y}) + \boldsymbol{R} 

を最小にする ***x*** を求めた。第二項には第一項に含まれていない情報が付加されている。上記の例では、「推定パラメタはほとんど零である」や「推定パラメタはばらつきが少ない」などである。その情報は、問題が与えられる前にすでに期待されていることだから、先験情報（あるいは事前情報）と呼ばれる。これに対応させて、第一項を事後情報ととらえ、逆問題をベイズ統計学の観点から考えることもできる。

一般に、第一項はデータのノイズに敏感で、推定パラメタに大きな変動や空間スケールの小さな構造をもたらす。一方、第二項は推定パラメタを滑らか・安定にする働きをもつ。問題に応じて、第一項と第二項の比は（ときに主観的に）決められる。

## 逆問題でよく扱われる方程式・解析手法
*第一種フレドホルム積分方程式
*有限要素法
*境界要素法
*フーリエ変換
*ラプラス変換
*随伴作用素

## 逆問題を応用した分野
* 弾道学
* レーダー、電波探知機
*非破壊検査(例 マイクロ波マンモグラフィ)、超音波探査、CTスキャン
*非鮮明な画像の復元・補完
*温度分布や熱伝導係数の推定
*地震学における震源断層すべり分布履歴の推定
*彗星のダストテイルの観測画像から、ダストの生成率、サイズ分布、放出速度を推定
*複数の反応率や検出器の波高分布から放射線のエネルギー分布の推定（アンフォールディング）
*翼型設計

## 脚注

### 注釈


### 出典

## 参考文献
* 日本機械学会(編):「逆問題のコンピュータアナリシス」、コロナ社、ISBN 978-4-339-04140-8 (1991年11月20日).
* 久保司郎：「逆問題」、培風館（計算力学とCAEシリーズ 10）、ISBN 4-563-03385-5 (1992年).
* 田中博、岡部政之、鈴木貴：「逆問題」、岩波書店（岩波講座 応用数学 3）、ISBN 978-4-00010513-2 (1993年6月8日).
* 登坂宣好、大西和榮、山本昌宏：「逆問題の数理と解法：偏微分方程式の逆解析」、東京大学出版会、ISBN 4-13-062906-9 (1999年).
* W.メンケ(著)、柳谷俊・塚田和彦(訳)：「離散インバース理論」、古今書院、ISBN 4-7722-1558-1 (1997年).
* Andreas Kirsch：*An Introduction to the Mathematical Theory of Inverse Problems*, Springer(Applied Mathematical Sciences 120)、ISBN 0-387-94530-X (1996年).
* W. Press *et al*：*Numerical Recipes in C* (2nd Ed.), Camridge University Press, ISBN 0-521-43108-5 (1992年) ※ 特に 18 章。
* 小島史男、上坂充(編著)：「電磁現象と逆問題」、養賢堂、ISBN 4-8425-9907-3 (1999年).
* 堤正義：「逆問題の数学」、共立出版、ISBN 4-320-01656-4 (2000年3月25日).
* 村上章、登坂宣好、堀宗朗、鈴木誠：「有限要素法・境界要素法による 逆問題解析：カルマンフィルタと等価介在物法の応用」、コロナ社、ISBN 978-4-339-05212-1 (2002年4月26日).
* 小國健二：「応用例で学ぶ逆問題と計測」、オーム社、ISBN 978-4-274-06829-4 (2011年2月19日).
* 堤正義：「逆問題：理論および数理科学への応用」、朝倉書店、ISBN 978-4-254-11824-7 (2012年10月25日).

## 関連項目
*正則化
*関数解析学
*数値解析