ブロック行列の逆行列

ある本でブロック行列で表現された行列の逆行列を目にした。証明は簡単で、掛けて単位行列になることを確認すればよい。しかし、どのように逆行列を導いたかは書かれておらず、長い期間分からないままだった。このままでは必要になるたびに逆行列を解説しているテキストを探さねばならず、しかも取り扱っているテキストは少ないため非常に不便である。たまたま授業の準備で基本変形を調べていた時に導出法に思い至った。

ブロック行列 X=(ABCD) の行列式と逆行列を求める。もちろん X の逆行列を考えるからには X は正方行列である。

A が正則行列の場合

まず、A が正則行列の場合を考える。このとき D は正則であるかどうかは分からないが、正方行列ではある。B,C は正方行列とは限らない。通常の逆行列を求める時をイメージすると次のように基本変形される。

(ABCD)→(E∗CD)→(E∗O∗)→(E∗OE)→(EOOE)

上の手順を具体化する。

(A−1OOE)(ABCD)=(EA−1BCD)(EO−CE)(EA−1BCD)=(EA−1BOD−CA−1B)

ここで、(EO−CE)(A−1OOE)X=(EA−1BOD−CA−1B) の行列式を取って、|X|=|A||D−CA−1B| を得る。

基本変形を更に進めるためには F=D−CA−1B が正則であることを要請する。

(EOOF−1)(EA−1BOF)=(EA−1BOE)(EA−1BOE)(E−A−1BOE)=(EOOE)

上の基本変形をまとめて逆行列を求める。

X−1=(EA−1BOE)(EOOF−1)(EO−CE)(A−1OOE)=(A−1−A−1BF−1CA−1A−1BF−1−F−1CA−1F−1)

実際に計算してみると X−1X=E が確認できるだろう。結果を次にまとめる。

ブロック行列で表された正方行列 X=(ABCD) について、A を正則行列とし、F=D−CA−1B とすると次が成り立つ。
(1) |X|=|A||F|
(2) F が正則な場合、X−1=(A−1−A−1BF−1CA−1A−1BF−1−F−1CA−1F−1)

(1) より、F が正則でない場合は X は正則でないことが分かる。

B が正則行列の場合

各ブロックについて正則性を仮定したバージョンが存在する。ここでは B が正則行列の場合を考える。この時、C は正方行列であるが、A,D は正方行列とは限らない。後の都合で G=C−DB−1A と置いておく。

(OEEO)(ABCD)=(CDAB)(EOOB−1)(CDAB)=(CDB−1AE)(E−DOE)(CDB−1AE)=(GOB−1AE)

(E−DOE)(EOOB−1)(OEEO)X=(GOB−1AE) から |X|=|B||G| を得る。さらに基本変形を続ける。

(G−1OOE)(GOB−1AE)=(EOB−1AE)(EO−B−1AE)(EOB−1AE)=(EOOE)

基本変形をまとめると

X−1=(EO−B−1AE)(G−1OOE)(E−DOE)(EOOB−1)(OEEO)=(−G−1DB−1G−1B−1AG−1DB−1+B−1−B−1AG−1)

結果を次にまとめる。

ブロック行列で表された正方行列 X=(ABCD) について、B を正則行列とし、G=C−DB−1A とすると次が成り立つ。
(1) |X|=|B||G|
(2) G が正則な場合、X−1=(−G−1DB−1G−1B−1AG−1DB−1+B−1−B−1AG−1)