逆問題のソースを表示

'''逆問題'''（ぎゃくもんだい、{{lang-en-short|inverse problem}}）とは、[[数学]]・[[物理学]]の一分野であり、入力（原因）から出力（結果、観測）を求める問題を'''順問題'''（じゅんもんだい、{{lang-en-short|direct problem}}）と呼び、その逆に出力から入力を推定する問題や入出力の関係性を推定する問題を逆問題と呼ぶ。

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

==概要==
順問題と逆問題は対になる概念であり、どちらが順でどちらが逆かというのは相対的な問題である。しかし対称的ではない。一般に、古くから問題として認識され研究が行われている方向のプロセスによるものを順問題とし、その逆方向のプロセスで解く方法は自明ではないのだが、それを解くことで何らかの[[工学]]的・その他の利用ができるような問題のことを逆問題と言う<ref group="注">物理現象の因果関係と、工学的応用では順方向と逆方向が逆の場合もある。わかりやすい例としては、3Dコンピュータグラフィックスにおいては、いわゆる「[[レイトレーシング]]」においてカメラ側から追跡したほうが容易なため、光源側からの追跡を「逆」と表現することがあり、ややこしい場合がある。</ref>。

単純な順問題・逆問題の例を示す。''f''(''x'') = ''x''<sup>2</sup> という関数について考える。''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'' <sup>-1</sup>(''x'') = √''x'' によって容易に得られる。しかし、一般には逆関数が容易にはわからない関数も多く、そういった場合を特に扱うのがこの分野である。

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

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

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

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

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

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

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

===ティホノフの正則化法===
非適切な問題の近似解を得る手法として最もよく使われるのが'''ティホノフの正則化法''' ({{lang-en-short|Tikhonov regularization method}}) である。

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

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

最適正則化パラメータの推定法として文献<ref>* 名取亮(編)：「数値計算法」、オーム社、ISBN 4-274-13153-X (1998年9月20日)の第8.6節「最適正則化パラメータ推定法」</ref>では、GCV法、L-カーブ法、Quasi-Optimal法を掲げている。

*関連項目： [[ヒルベルト空間]]、[[コンパクト (数学)|コンパクト]]

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

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

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

式で書けば、ノイズを最小にするには
:<math>J = \sum^N_{i=1} \sum^N_{j=1} n_i S^{-1}_{ij} n_j </math>
あるいは、行列表示して（上付き添字 <math>\intercal</math> は[[転置行列]]を表す）
:<math>
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})
</math>
なる ''J'' を最小にする '''''x''''' を決める問題になるが、行列 '''''K''''' は行より列が多く、'''''Kx'''''='''''y''''' の解が無数にあるという状況になる。そのため、正則化を行って解をひとつに定める。以下にいくつかの正則化の方法を紹介する。以下の議論で本質的に重要でないため、ノイズは分散 1 でそれぞれ無相関なものとする（つまり'''''S'''''は[[単位行列]]）。

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

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

ここに、'''''B''''' は
:<math>
\begin{bmatrix}
-1 & 1 & 0 & \ldots & 0 \\
0  & -1 & 1 & \ldots & 0 \\
\vdots & & \ddots & & \vdots \\
0 & 0 & \ldots & -1 & 1
\end{bmatrix}
</math>
なる成分を持つ。

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

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

:<math> \hat{\boldsymbol{x}} = \boldsymbol{L} \boldsymbol{y} </math>

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

:<math> \hat{\boldsymbol{x}} = \boldsymbol{L}(\boldsymbol{K}\boldsymbol{x} + \boldsymbol{n}) </math>

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

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

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

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

:<math> \boldsymbol{L} \mathrm{E}(\boldsymbol{n}\boldsymbol{n}^\intercal)\boldsymbol{L}^\intercal </math>

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

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

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

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

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

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

==逆問題を応用した分野==
* [[弾道学]]
* [[レーダー]]、[[電波探知機]]
*[[非破壊検査]](例 マイクロ波マンモグラフィ<ref>{{Cite web|和書
  |url=https://www.med-device.jp/repository/20200204-zenkoku-kaigi_27_igs.pdf
  |title= 散乱逆問題の解析解発見とマイクロ波マンモグラフィの実現
  |publisher= 神戸大学、(株)Integral Geometry Science
  |accessdate=2021-03-21
  |language=jp
  |format=pdf}}</ref><ref>{{Cite web|和書
  |url=https://www.amed.go.jp/content/000051771.pdf
  |title= 日本医療研究開発機構 医療分野研究成果展開事業（先端計測分析技術・機器開発プログラム）事後評価報告書 
  |publisher= 国立研究開発法人日本医療研究開発機構
  |accessdate=2021-03-21
  |language=jp
  |format=pdf}}</ref><ref>{{Cite web|和書
  |url=https://www.kobe-u.ac.jp/research_at_kobe/NEWS/people/researcher0012_ex01.html
  |title=【対談】木村建次郎教授 × ニュースキャスター 膳場貴子さん
  |publisher= 神戸大学
  |accessdate=2021-03-21
  |language=jp
  |format=pdf}}</ref>)、[[超音波]]探査、[[コンピュータ断層撮影|CT]]スキャン
*非鮮明な[[画像]]の復元・補完
*[[温度]]分布や[[熱伝導]]係数の推定
*[[地震学]]における[[震源]][[断層]]すべり分布履歴の推定
*[[彗星]]の[[ダストテイル]]の観測画像から、ダストの生成率、サイズ分布、放出速度を推定
*複数の反応率や検出器の波高分布から放射線のエネルギー分布の推定（アンフォールディング）
*[[翼型]]設計<ref>{{cite journal|和書
|title="温故知新（先人の教え）～低速の翼に関する話題～ 第５回 補遺１：翼型理論と逆問題設計法の考察"
|journal=ながれ : 日本流体力学会誌 = Nagare : journal of Japan Society of Fluid Mechanics
|publisher=日本流体力学会
|location=東京
|date=2023-08
| volume = 42
| issue = 4
|pages=259-268
|url=https://www.nagare.or.jp/download/noauth.html?d=42-4_259_rensai3.pdf&dir=190 
|author=吉田憲司
|language=jp
|format=pdf
|ISSN=0286-3154
|crid=1520860553742090752
|accessdate=2024-10-31
}}</ref>

== 脚注 ==
{{脚注ヘルプ}}

=== 注釈 ===
<references group="注"/>

=== 出典 ===
{{Reflist}}


==参考文献==
* 日本機械学会(編):「逆問題のコンピュータアナリシス」、コロナ社、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日).

==関連項目==
*[[正則化]]
*[[関数解析学]]
*[[数値解析]]

{{Normdaten}}
{{デフォルトソート:きやくもんたい}}
[[Category:応用数学]]
[[Category:数学の問題]]
[[Category:逆問題|*]]
[[Category:数学に関する記事]]