スポンサーリンク

「余因子行列の性質〜逆行列が存在することの必要十分条件〜」【線型代数学の基礎シリーズ】行列式編 その6

線型代数学

本記事の内容

本記事は、余因子行列の性質である逆行列が存在することの必要十分条件について解説する記事です。

本記事を読むにあたり、余因子、余因子行列について知っている必要があるため、以下の記事も合わせてご覧ください。

本記事で言いたいこと

本記事で言いたい定理はシンプルです。
何を言いたいかというと、

定理0.

正方行列AAが正則であるための必要十分条件は、det(A)0det(A)0である。このとき、AAの逆行列A1A1A1=1det(A)˜AA1=1det(A)~A であたえられる。ただし、˜A~AAAの余因子行列である。

です。

要するに、det(A)det(A)を計算して00でなければ必ず逆行列があって、余因子行列でもって逆行列が求まる、という話です。

逆行列のチャラい復習

逆行列をチャラく復習します。
逆行列とは以下でした。

逆行列

nNnNAAnn次正方行列とする。このとき、 AB=BA=InAB=BA=In を満たすようなBBが存在するとき、AA正則である、または正則行列であるという。 また、このときBBAA逆行列といい B=A1B=A1 と書く。

詳しくは、【線型代数学の基礎シリーズ】行列編 その3を御覧ください。

余因子、余因子行列のチャラい復習

余因子および余因子行列をチャラく復習します。
余因子および余因子行列とは以下でした。

余因子

nNnNnn次正方行列A=(aij)A=(aij)から第ii行と第jj列を取り除いて得られたたn1n1次正方行列をAijAijと書く。すなわち、AijAijjA=(a11a1ja1nai1aijainan1anjann)i としたとき、第i行と第j列を取り除いて得られた行列    j   Aij=(a11a1 j1a1 j+1a1nai1 1ai1 j1ai1 j+1ai1 nai+1 1ai+1 j1ai+1 j+1ai+1 nan1an j1an j+1ann)i である。
 このとき、 ~aij=(1)i+jdet(Aij) とおき、これを行列Aにおけるaij余因子という。

余因子行列

nNA=(aij)n次正方行列とする。このとき、aijの余因子˜aijを成分とする行列 (˜a11˜a12˜a1n˜a21˜a22˜a2n˜an1˜an2˜ann) の転置行列 ˜A=(˜a11˜a21˜an1˜a12˜a22˜an2˜a1n˜a2n˜ann)A余因子行列という。

詳しくは【線型代数学の基礎シリーズ】行列式編 その5を御覧ください。

定理0.の証明のために必要な事実

ここでは本記事で言いたい定理0.を示すために必要な事実を証明します。

一言でいえば、

行列の積の行列式は、それぞれの行列の行列式の積と等しい。

ということです。
主張はシンプルですが、実は結構重要だったりします。

では、定理を明示しましょう。

定理1.

nNAおよびBが共にn次正方行列とする。 このとき、次が成り立つ。 det(AB)=det(A)det(B)

(定理7.へジャンプ)

定理1.の証明は割と真正直にやれば証明できます。

定理1.の証明

nNA=(aij)およびB=(bij)n次正方行列とします。
このとき、
A=(a11a1nam1amn),B=(b1bn)
と書いたとします。
ただし、b1,,bnBの行ベクトルです。
すなわち、
bj=(bj1 bj2  bjn)
です。

このとき、AB(n,n)型なので、ABが定まります。
従って、
AB=(a11a1nam1amn)(b1bn)=(a11b1++a1nbnan1b1++annbn)

ここで、以下の2つの定理を繰り返し使います。

定理2.

nNA=(aij)n次正方行列とする。Aの第i行が、2つの行ベクトルの和ならば、行列式は他の行は同じで、第i行は各々のベクトルを取った行列の行列式の和になる。 すなわち、 (iN;1in) s.t. (ai1,,ain)=(bi1+ci1,,bin+cin) としたとき、 |a11a1nbi1+ci1bin+cinan1ann|=|a11a1nbi1binan1ann|+|a11a1nci1cinan1ann| が成り立つ。

定理2.の証明は【線型代数学の基礎シリーズ】行列式編 その3を御覧ください。

定理3.

n次正方行列A=(aij)の1つの行をcC倍すると、行列式はc倍となる。すなわち、 |a11a1ncai1cainan1ann|=c|a11a1nai1ainan1ann| が成り立つ。

定理3.の証明は【線型代数学の基礎シリーズ】行列式編 その3を御覧ください。

まず、
(a11b1++a1nbnan1b1++annbn)
第1行に目をつけてみます。
すると、定理2.および定理3.から
|a11b1++a1nbnan1b1++annbn|=a11|b1an1b1++annbn|++a1n|bnan1b1++annbn|=nj1=1a1j1|bj1an1b1++annbn|
が成り立ちます。

この操作をすべての行に対して繰り返し行います。
すると、
|a11b1++a1nbnan1b1++annbn|=njn=1an jnnjn1=1an1 jn1nj2=1a2 j2nj1=1a1j1|bj1bj2bjn1bjn|=njn=1njn1=1nj2=1nj1=1a1j1a2j2an1 jn1an jn|bj1bj2bjn1bjn|

となります。
ここで、和はj1,j2,,jnがそれぞれ1からnまでを動くわけですので、nn個の項が出現します。
ここで、detに出現する和はσSnでの和ですので、項数としてはn!個です。
nNのときnn>n!ですので、今回計算したものの方が項数が多いことに注意です。
何が言いたいか、というと、「余分なもの(というとちょっと語弊があるけれど)まで足しているZE」ということです。

実は、次の事実が成り立っているのでした。

命題4.

2つの行が等しい行列の行列式は0である。

この命題の証明は【線型代数学の基礎シリーズ】行列式編 その3を御覧ください。

命題4.により、bj1,bj2,,bjn1,bjn,のうち、同じものがあれば、
|bj1bj2bjn1bjn|=0
ということになります。
すなわち、j1,j2,,jn1,jnがすべて相異なる場合の和を考えれば良いことになります。
先程述べた「余分なもの」というのは「j1,j2,,jn1,jnのうち少なくとも1ペアが等しい場合」を指しています。

さて、j1,j2,,jn1,jnがすべて相異なる場合というのは、j1,j2,,jn1,jnがちょうど1,2,,nの順列となる場合ということになります。
というのも、j1,j2,,jn1,jnがすべて相異なって、それぞれが1からnを動くからです。
しかも、和を取れば、ちょうどすべての順列での和ということになります。
従って、
σ=(12nj1j2jn)
としたとき、

det(AB)=njn=1njn1=1nj2=1nj1=1a1j1a2j2an1 jn1an jn|bj1bj2bjn1bjn|=σSna1j1anjn|bj1bjn|=σSna1σ(1)anσ(n)|bσ(1)bσ(n)|

ここで、
|bσ(1)bσ(n)|
に対して、以下の定理を適用させます。

定理5.

nNA=(aij)n次正方行列とする。Aの各行の順序を置換τによって変更すると、行列式はsgn(τ)になる。すなわち、 τ=(12nk1k2kn) としたとき、 |ak11ak12ak1nak21ak22ak2nakn1akn2aknn|=sgn(τ)|a11a12a1na21a22a2nan1an2ann|=sgn(τ)det(A) が成り立つ。

この定理5.の証明は【線型代数学の基礎シリーズ】行列式編 その3を御覧ください。

定理5.により、
|bσ(1)bσ(n)|=sgn(σ)|b1bn|
です。
以上のことから、
det(AB)=σSna1σ(1)anσ(n)|bσ(1)bσ(n)|=σSna1σ(1)anσ(n)sgn(σ)|b1bn|=(σSnsgn(σ)a1σ(1)anσ(n))|b1bn|
が成り立ち、
σSnsgn(σ)a1σ(1)anσ(n)=det(A),|b1bn|=det(B)
なわけですので、
det(AB)=det(A)det(B)
です。

定理1.の証明終わり

某五番隊隊長から「ひゃあ、こら長いわ」とツッコミが来そうですね。

逆行列が存在することの必要十分条件

準備が整いましたので、いよいよ本題に入ります。
とはいえ、今まで証明したことをまとめて使えばなんてこと無いんですね、実は。

定理0.

正方行列Aが正則であるための必要十分条件は、det(A)0である。このとき、Aの逆行列A1A1=1det(A)˜A であたえられる。ただし、˜AAの余因子行列である。

定理0.の証明

nNAn次正方行列だとします。

(「Aが正則det(A)0」の証明)

Aが正則行列だとします。
このとき、Aが正則であることから、ある行列n次正方行列Bが存在して、
AB=In
(ただし、Inn次単位行列とします)を満たします。

先程示した定理1.により、行列の積の行列式はそれぞれの行列式の積と等しいので、
det(A)det(B)=det(AB)=det(In)=1
となるので、det(A)0です。
実際、仮にdet(A)=0ならば、det(A)det(B)=0となって矛盾です。

(「det(A)0Aは正則」の証明)

det(A)0だとします。
このとき、1det(A)˜Aという行列を考えることができます。
そこで、
B=1det(A)˜A
とします。
このBAの逆行列であれば、証明完了です。
すなわち、
AB=BA=In
であれば証明完了です。

AB=A(1det(A)˜A)=1det(A)A˜A
となります。
ここで、次の事実を使います。

定理6.

nNA=(aij)n次正方行列、˜AAの余因子行列とする。 このとき、次が成り立つ。 A˜A=A˜A=(det(A)OOdet(A))=det(A)(1OO1)=det(A)In ただし、Inn次の単位行列である。

定理6.の証明は【線型代数学の基礎シリーズ】行列式編 その5を御覧ください。

定理6.から、A˜A=˜AA=det(A)Inなのですから、
AB=1det(A)A˜A=1det(A)det(A)In=In
です。

同様にして、定理6.を使うことで
BA=(1det(A)˜A)A=1det(A)˜AA=1det(A)det(A)In=In
です。

従って、Aは正則行列で、かつA1
A1=B=1det(A)˜A
です。

定理0.の証明終わり

次に紹介するのが、逆行列の行列式に対する事実です。

定理7.

正則行列Aに対して、以下が成り立つ。 det(A1)=det(A)1

定理1.を使えば一瞬です。

定理7.の証明

定理1.から、
1=det(In)=det(AA1)=det(A)det(A1)
であるので、
det(A)det(A1)=1
だから、
det(A1)=1det(A)=det(A)1
となるから成り立ちます。

定理7.の証明終わり

いっちょ計算してみっか

いっちょ計算してみましょう。

例8.
A=(100201154)
が正則行列かを判定して、正則だったら逆行列を求めてみましょう。

det(A)=|100201154|=|0154|=0×4(5×(1))=50
なので、Aは正則行列です。
では、逆行列を求めてみましょう。

  • a11の余因子˜a11(第1行と第1列を抜いたとき)
    A11=(0154),˜a11=(1)2det(A11)=5
  • a12の余因子˜a12(第1行と第2列を抜いたとき)
    A12=(2114),˜a12=(1)3det(A12)=7
  • a13の余因子˜a13(第1行と第3列を抜いたとき)
    A13=(2015),˜a13=(1)4det(A13)=10
  • a21の余因子˜a21(第2行と第1列を抜いたとき)
    A21=(0054),˜a21=(1)3det(A21)=0
  • a22の余因子˜a22(第2行と第2列を抜いたとき)
    A22=(1014),˜a22=(1)4det(A22)=4
  • a23の余因子˜a23(第2行と第3列を抜いたとき)
    A23=(1015),˜a23=(1)5det(A23)=5
  • a31の余因子˜a31(第3行と第1列を抜いたとき)
    A31=(0001),˜a31=(1)4det(A31)=0
  • a32の余因子˜a32(第3行と第2列を抜いたとき)
    A32=(1021),˜a32=(1)5det(A32)=1
  • a33の余因子˜a33(第3行と第3列を抜いたとき)
    A33=(1020),˜a33=(1)6det(A33)=0

故に、Aの余因子行列˜A
˜A=(˜a11˜a21˜a31˜a12˜a22˜a32˜a13˜a23˜a33)=(5007411050)
となるので、Aの逆行列A1
A1=1det(A)˜A=15(5007411050)=(100754515210)
です。
本当にこの行列がAの逆行列かを確かめてみます。
AA1=(100201154)(100754515210)=(1+0+00+0+00+0+02+020+0+10+0+017+80+440+1+0)=(100010001)=I3

A1A=(100754515210)(100201154)=(1+0+00+0+00+0+075+85150+0+1045+4522+00+0+00+1+0)=(100010001)=I3
確かに、逆行列です。

今回は、余因子行列の性質である「逆行列の存在の必要十分条件」について解説しました。
「余因子行列なんてそんな面倒なモノがなんで必要なんだ?」と思ったかもしれませんが、逆行列を求める際に誠に有用だからです。
そして、行列式は逆行列が存在するかどうかを判定する際に非常に有用である、ということが今回わかりました。

次回から線型空間(ベクトル空間)の話をします。

乞うご期待!質問、コメントなどお待ちしております!

コメントをする

タイトルとURLをコピーしました