擬似逆行列

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

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

定義[編集]

m × n 行列 A に対し、A の随伴行列（複素共軛かつ転置行列）を A^* とするとき、以下の4条件を満足する n × m 行列 A⁺ はただ一つ定まる：

{\begin{aligned}AA^{+}A&=A,\\A^{+}AA^{+}&=A^{+}.\end{aligned}}

{\begin{aligned}(AA^{+})^{*}&=AA^{+},\\(A^{+}A)^{*}&=A^{+}A.\end{aligned}}

この行列 A⁺ を A の擬似逆行列と呼ぶ。A が正則でなくとも A⁺ は定まるが、A が正則ならば逆行列 A⁻¹ はこの条件を満たす。ゆえに擬似逆行列の概念は逆行列の概念の一般化を与えていることがわかる。

擬似逆行列は以下のような性質を持つ。

$(A^{+})^{+}=A$
$(A^{\intercal })^{+}=(A^{+})^{\intercal },(AA^{\intercal })^{+}=(A^{+})^{\intercal }A^{+}$
$\operatorname {rank} A=\operatorname {rank} B=n\implies (AB)^{+}=B^{+}A^{+}$
$A^{+}=(A^{\intercal }A)^{+}A^{\intercal }=A^{\intercal }(AA^{\intercal })^{+}$
$m\times n$ 行列 A に対して
${\begin{aligned}\operatorname {rank} A=m&\implies A^{+}=A^{\intercal }(AA^{\intercal })^{-1}\\\operatorname {rank} A=n&\implies A^{+}=(A^{\intercal }A)^{-1}A^{\intercal }\end{aligned}}$
A の特異値分解を $A=U\Sigma V^{\intercal }$ とすると、
$A=U\Sigma V^{\intercal }\implies A^{+}=V\Sigma ^{+}U^{\intercal }$

が成立する。（ $\Sigma ^{+}$ の成分は $\sigma _{ii}^{+}$ 、 $\Sigma$ の成分は $\sigma _{ii}$ とすると、 $\sigma _{ii}^{+}=1/\sigma _{ii}$ である。）

$m\times n$ $m\times n$ 行列 $A$ $A$ に対して
- n 次正方行列 $A^{+}A$ は、 $A$ の零空間の直交補空間 $(\ker A)^{\bot }$ への直交射影である。
- n 次正方行列 $I_{n}-A^{+}A$ は、 $\ker A$ への直交射影である。
$A$ $A$ を $m\times n$ $m\times n$ 行列とする。連立一次方程式 $Ax=b$ $Ax=b$ に対して
- 方程式が解を持つとき
  $k$ を任意の $n$ 次元列ベクトルとして、すべての解は $x=A^{+}b+(I_{n}-A^{+}A)k$ と表せる。ノルム $\|x\|$ が最小の解は $A^{+}b$ で与えられる。 $A$ が正則なら $A^{+}=A^{-1}$ で、ただ一つの解を持つ。
- 方程式が解を持たないとき
  前述の $x$ は $\|Ax-b\|^{2}$ を最小にするベクトル（最小２乗解）である。