擬似逆行列のソースを表示

ムーア-ペンローズの'''擬似逆行列'''（ぎじぎゃくぎょうれつ、''pseudo-inverse matrix''）は[[線型代数学]]における[[逆行列]]の概念の一般化である。'''擬逆行列'''、'''一般化逆行列'''、'''一般逆行列'''（{{lang-en-short|generalized inverse}}）ともいう。また擬は疑とも書かれる。

[[連立一次方程式]]の解を簡潔に表現するものとして[[逆行列]]の概念は重要であり、逆行列を持つ行列は、[[可逆]]あるいは[[正則行列|正則]]であると言われる。正則でない行列の場合にも逆行列のような都合のよい行列として擬逆の概念を導入する。[[ロボット工学]]に関していうならば、動特性の同定や冗長ロボットの制御などで良く用いられている。

== 定義 ==
{{see also|ムーア・ペンローズ逆行列}}

''m'' &times; ''n'' 行列 ''A'' に対し、''A'' の[[随伴行列]]（複素共軛かつ[[転置行列]]）を ''A''<sup>*</sup> とするとき、以下の4条件を満足する ''n'' &times; ''m'' 行列 ''A''<sup>+</sup> はただ一つ定まる：
* ''A'' と ''A''<sup>+</sup> は互いに広義可逆元である：
:<math>\begin{align}
A A^+ A &= A, \\
A^+ A A^+ &= A^+.
\end{align}</math>
*  ''A'' ''A''<sup>+</sup> および ''A''<sup>+</sup>''A'' は[[エルミート行列]]である：
:<math>\begin{align}
(A A^+)^* &= A A^+, \\
(A^+ A)^* &= A^+ A.
\end{align}</math>
この行列 ''A''<sup>+</sup> を ''A'' の'''擬似逆行列'''と呼ぶ。''A'' が正則でなくとも ''A''<sup>+</sup> は定まるが、''A'' が正則ならば逆行列 ''A''<sup>&minus;1</sup> はこの条件を満たす。ゆえに擬似逆行列の概念は逆行列の概念の一般化を与えていることがわかる。

== 性質 ==
擬似逆行列は以下のような性質を持つ。
* <math>(A^+)^+ = A</math>
* <math>(A^{\intercal})^+ = (A^+)^{\intercal}, (A A^{\intercal})^+ = (A^+)^{\intercal} A^+</math>
* <math>\operatorname{rank}A=\operatorname{rank}B=n \implies (AB)^+ = B^+ A^+</math>
* <math>A^+ = (A^{\intercal} A)^+ A^{\intercal} = A^{\intercal} (AA^{\intercal})^+</math>
* <math>m \times n</math> 行列 ''A'' に対して
*: <math>\begin{align}
\operatorname{rank}A = m &\implies A^+ = A^{\intercal} (AA^{\intercal})^{-1} \\
\operatorname{rank}A = n &\implies A^+ = (A^{\intercal}A)^{-1} A^{\intercal}
\end{align}</math>
* ''A'' の[[特異値分解]]を <math>A = U \Sigma V^{\intercal}</math>とすると、
*: <math>A = U \Sigma V^{\intercal} \implies A^+ = V \Sigma^+ U^{\intercal}</math>
が成立する。
（<math>\Sigma^+</math> の対角成分を <math>\sigma_{ii}^+</math>、<math>\Sigma</math> の対角成分を <math>\sigma_{ii}</math> とすると、
<math>\sigma_{ii}^+ = 1/\sigma_{ii}</math> である。）
* <math>m \times n</math> 行列 <math>A</math> に対して
** ''n'' 次正方行列  <math>A^+A</math> は、<math>A</math> の[[零空間]]の[[直交補空間]] <math>(\ker A)^\bot</math>への直交射影である。
** ''n'' 次正方行列 <math>I_n-A^+A</math> は、<math>\ker A</math> への直交射影である。
* <math>A</math>を<math>m \times n</math> 行列とする。連立一次方程式 <math>Ax=b</math> に対して
** 方程式が解を持つとき<br /><math>k</math>を任意の<math>n</math>次元列ベクトルとして、すべての解は<math> x=A^+b + (I_n - A^+ A)k</math>と表せる。ノルム <math>\|x\|</math> が最小の解は<math>A^+ b</math> で与えられる。<math>A</math> が正則なら<math>A^+=A^{-1}</math>で、ただ一つの解を持つ。
** 方程式が解を持たないとき<br />前述の <math>x</math> は<math>\|Ax-b\|^2</math>を最小にするベクトル（[[最小二乗法|最小２乗解]]）である。

== 例 ==
===スカラー===
スカラーの場合にも擬似逆行列を定義できる。スカラーを行列として扱うことになる。{{math|''&lambda;''}} が0ならば、その擬似逆行列は0であり、{{math|''&lambda;''}} がそれ以外の数ならば、 その擬似逆行列は {{math|''&lambda;''}} の逆数になる:
:<math> \lambda^+ = 
\begin{cases}
0 & (\lambda=0),\\
\lambda^{-1} & (\lambda\ne 0). 
\end{cases} </math>
===ベクトル===
零ベクトルの擬似逆行列は転置された零ベクトルである。零ベクトルでないベクトルの擬似逆行列はそのベクトルの大きさの2乗で割られた、随伴ベクトルである:
:<math>x^+ = 
\begin{cases} 
0^* & (x = 0),\\
 (x^* x)^{-1} x^* & (x\ne 0). 
\end{cases} </math>

===列が線形独立である場合===
<math>A</math> の各列が[[線形独立]]（このとき <math>m \ge n</math> である）ならば、<math>A^*A</math> は可逆である。この場合、擬似逆行列は次のようになる:
:<math>A^+ = (A^*A)^{-1}A^*</math>.
これから <math>A^+</math> が <math>A</math> の左逆元であることがわかる: つまり <math> A^+ A = I_n</math> 

===行が線形独立である場合===
<math>A</math> の各行が[[線形独立]]（このとき <math>m \le n</math> である）ならば、<math>A A^*</math> は可逆である。この場合、擬似逆行列は次のようになる:
:<math>A^+ = A^*(AA^*)^{-1}</math> 
これから <math>A^+</math> が <math>A</math> の右逆元であることがわかる: つまり <math> AA^+  = I_m</math> 

===2次正方行列===
2次正方行列
:<math>A=
\begin{bmatrix}
a&b\\
c&d
\end{bmatrix}
</math>
の擬似逆行列は <math>ad-bc\ne 0</math> のとき、
:<math>A^+=A^{-1}=
\frac{1}{ad-bc}
\begin{bmatrix}
d&-b\\
-c&a
\end{bmatrix}
</math>
である。 <math>ad-bc=0</math> のとき、 <math>A\ne O</math> のときは
:<math>A^+=
\frac{1}{|a|^2+|b|^2+|c|^2+|d|^2}
\begin{bmatrix}
\bar{a}&\bar{c}\\
\bar{b}&\bar{d}
\end{bmatrix}
</math>
となる。 <math>A=O</math> のときは
:<math>A^+=O=
\begin{bmatrix}
0&0\\
0&0
\end{bmatrix}
</math> 
である。

== 参考文献 ==
{{参照方法|date=2016年3月}}
* 「ロボット制御基礎論」(著者：吉川恒夫)
* {{cite book
|last1      = Harville
|first1     = David A
|year       = 1997
|title      = Matrix algebra from a statistician's perspective
|url        = {{google books|kZGBQijgGV8C|plainurl=yes|page=497}}
|publisher  = Springer-Verlag
|isbn       = 0-387-94978-X 
|mr         = 1467237
|zbl        = 0881.15001 
|ref        = harv
}}
* {{cite book
|和書
|author      = 岩井斉良
|year       = 1995
|title      = 基礎課程線形代数
|publisher  = 学術図書出版社
|isbn       = 978-4-87361-194-5
}}

== 関連項目 ==
* [[線型代数学]]、[[逆行列]]
* [[逆元]]（半群における擬逆性）

{{DEFAULTSORT:きしきやくきようれつ}}
[[Category:行列]]
[[Category:数学に関する記事]]