正則行列の判定と逆行列の導出
正方行列\(A\in M_{n,n}\left( \mathbb{R} \right) \)に対して、以下の条件\begin{equation*}AB=BA=I_{n}
\end{equation*}を満たす正方行列\(B\in M_{n,n}\left( \mathbb{R} \right) \)が存在する場合には、\(A\)を正則行列と呼びます。ただし、\(I_{n}\in M_{n,n}\left( \mathbb{R} \right) \)は単位行列です。
正方行列\(A\in M_{n,n}\left( \mathbb{R} \right) \)が正則である場合には、定義より、以下の条件\begin{equation*}AB=BA=I_{n}
\end{equation*}を満たす正方行列\(B\in M_{n,n}\left( \mathbb{R} \right) \)が存在します。これを\(A\)の逆行列と呼び、\begin{equation*}A^{-1}
\end{equation*}で表記します。つまり、正方行列\(A\)の逆行列\(A^{-1}\)とは、以下の条件\begin{equation*}AA^{-1}=A^{-1}A=I_{n}
\end{equation*}を満たす正方行列として定義されます。
正方行列が正則であることは様々な形で表現可能です。これまで明らかになったものをまとめます。
- \(A\)は正則行列である。
- \(A\)の行標準形が単位行列\(I_{n}\in M_{n,n}\left( \mathbb{R} \right) \)である。
- \(n=\mathrm{rank}\left( A\right) \)が成り立つ。
- \(A\)を係数行列とする変数\(x_{1},x_{2},\cdots ,x_{n}\)に関する同次連立1次方程式\begin{equation*}\left\{ \begin{array}{c}
a_{11}x_{1}+a_{12}x_{2}+\cdots +a_{1n}x_{n}=0 \\
a_{21}x_{1}+a_{22}x_{2}+\cdots +a_{2n}x_{n}=0 \\
\vdots \\
a_{n1}x_{1}+a_{n2}x_{2}+\cdots +a_{nn}x_{n}=0\end{array}\right.
\end{equation*}が与えられているものとする。ただし、\(a_{ij}\in \mathbb{R} \)は係数であり\(x_{i}\in \mathbb{R} \)は変数である。この同次連立1次方程式の解は自明解\begin{equation*}\left(
\begin{array}{c}
x_{1} \\
x_{2} \\
\vdots \\
x_{n}\end{array}\right) =\left(
\begin{array}{c}
0 \\
0 \\
\vdots \\
0\end{array}\right)
\end{equation*}だけである。
正則行列\(A\)が正則であるものとします。この場合、\(A\)の行標準形は単位行列\(I_{n}\)であるため、\(A\)に対してガウス・ジョルダンの消去法を適用すれば単位行列\(I_{n}\)が得られます。消去法を構成する個々の行基本操作に対応する行基本行列を\(E_{1},E_{2},\cdots ,E_{p}\)で表記します。つまり、\begin{equation*}I_{n}=E_{p}\cdots E_{2}E_{1}A
\end{equation*}が成り立つということです。このとき、\(A\)の逆行列が、\begin{equation*}A^{-1}=E_{p}\cdots E_{2}E_{1}I_{n}
\end{equation*}として定まることが保証されるため、\(A\)から\(I_{n}\)を得るために適用した行基本操作をそのまま\(I_{n}\)に対して適用すれば\(A^{-1}\)が得られます。
A=\begin{pmatrix}
2 & 5 \\
1 & 3\end{pmatrix}\end{equation*}にガウス・ジョルダンの消去法を適用すると、\begin{eqnarray*}
A &=&\begin{pmatrix}
2 & 5 \\
1 & 3\end{pmatrix}
\\
&\rightarrow &\begin{pmatrix}
-5+6 & -15+15 \\
1 & 3\end{pmatrix}\quad \because R_{1}\rightarrow -5R_{2}+3R_{1} \\
&=&\begin{pmatrix}
1 & 0 \\
1 & 3\end{pmatrix}
\\
&\rightarrow &\begin{pmatrix}
1 & 0 \\
1-1 & 3-0\end{pmatrix}\quad \because R_{2}\rightarrow R_{2}-R_{1} \\
&=&\begin{pmatrix}
1 & 0 \\
0 & 3\end{pmatrix}
\\
&\rightarrow &\begin{pmatrix}
1 & 0 \\
0 & 1\end{pmatrix}\quad \because R_{2}\rightarrow \frac{1}{3}R_{2} \\
&=&I_{2}
\end{eqnarray*}を得ます。つまり、\(A\)の行標準形は単位行列\(I_{2}\)であるため\(A\)は正則です。そこで、単位行列\(I_{2}\)に対しても同様の行基本変形\begin{eqnarray*}R_{1} &\rightarrow &-5R_{2}+3R_{1} \\
R_{2} &\rightarrow &R_{2}-R_{1} \\
R_{2} &\rightarrow &\frac{1}{3}R_{2}
\end{eqnarray*}を適用すると、\begin{eqnarray*}
I_{2} &=&\begin{pmatrix}
1 & 0 \\
0 & 1\end{pmatrix}
\\
&\rightarrow &\begin{pmatrix}
3 & -5 \\
0 & 1\end{pmatrix}\quad \because R_{1}\rightarrow -5R_{2}+3R_{1} \\
&\rightarrow &\begin{pmatrix}
3 & -5 \\
-3 & 6\end{pmatrix}\quad \because R_{2}\rightarrow R_{2}-R_{1} \\
&\rightarrow &\begin{pmatrix}
3 & -5 \\
-1 & 2\end{pmatrix}\quad \because R_{2}\rightarrow \frac{1}{3}R_{2}
\end{eqnarray*}となるため、\begin{equation*}
A^{-1}=\begin{pmatrix}
3 & -5 \\
-1 & 2\end{pmatrix}\end{equation*}を得ます。実際、\begin{eqnarray*}
AA^{-1} &=&\begin{pmatrix}
2 & 5 \\
1 & 3\end{pmatrix}\begin{pmatrix}
3 & -5 \\
-1 & 2\end{pmatrix}=\begin{pmatrix}
1 & 0 \\
0 & 1\end{pmatrix}
\\
A^{-1}A &=&\begin{pmatrix}
3 & -5 \\
-1 & 2\end{pmatrix}\begin{pmatrix}
2 & 5 \\
1 & 3\end{pmatrix}=\begin{pmatrix}
1 & 0 \\
0 & 1\end{pmatrix}\end{eqnarray*}であるため、\begin{equation*}
AA^{-1}=A^{-1}A=I_{2}
\end{equation*}が成立しています。
行列式を用いた正則行列の判定
正方行列が正則であることは行列式を用いて判定することもできます。順番に解説します。
行列の行基本操作とは以下の3種類の操作\begin{eqnarray*}
&&\left( R_{1}\right) \ R_{i}\leftrightarrow R_{j} \\
&&\left( R_{2}\right) \ R_{i}\rightarrow kR_{i}\quad \left( k\not=0\right)
\\
&&\left( R_{3}\right) \ R_{i}\rightarrow R_{i}+kR_{j}
\end{eqnarray*}の総称です。行列\(A\)に対して行基本操作を有限回適用することで行列\(B\)が得られる場合、\(A\)は\(B\)と行同値であると言います。
単位行列\(I_{n}\in M_{n,n}\left( \mathbb{R} \right) \)に対して行基本操作\(\left( R_{1}\right) ,\left( R_{2}\right) ,\left( R_{3}\right) \)の中の1つを適用することにより行列\(E\in M_{n,n}\left( \mathbb{R} \right) \)が得られる場合、このような行列\(E\)を行基本行列と呼びます。
正方行列\(A\in M_{n,n}\left( \mathbb{R} \right) \)と行基本行列\(E\in M_{n,n}\left( \mathbb{R} \right) \)をそれぞれ任意に選んだとき、以下の関係\begin{equation*}\left\vert EA\right\vert =\left\vert E\right\vert \left\vert A\right\vert
\end{equation*}が成り立ちます。つまり、正方行列\(A\)に対して行基本行列\(E\)に対応する行基本操作を行うことにより得られる行列の行列式の値\(\left\vert EA\right\vert \)は、個々の行列の行列式の値どうしの積\(\left\vert E\right\vert \left\vert A\right\vert \)と一致します。
\end{equation*}であるとともに、\begin{equation*}
\left\vert EA\right\vert =\left\vert E\right\vert \left\vert A\right\vert
\end{equation*}が成り立つ。
以上の命題を踏まえた上で以下を示します。
\end{equation*}と表されるものとする。このとき、\begin{eqnarray*}
&&\left( a\right) \ \left\vert B\right\vert =\left\vert E_{p}\right\vert
\cdots \left\vert E_{2}\right\vert \left\vert E_{1}\right\vert \left\vert
A\right\vert \\
&&\left( b\right) \ \left\vert B\right\vert \not=0\Leftrightarrow \left\vert
A\right\vert \not=0
\end{eqnarray*}がともに成り立つ。
以上の命題を踏まえた上で以下を示します。
\end{equation*}が成り立つことと、\(A\)が正則行列であることは必要十分である。
A=\begin{pmatrix}
2 & 5 \\
1 & 3\end{pmatrix}\end{equation*}が正則行列であることは先に確認した通りです。同じことを行列式を用いて示します。具体的には、\begin{eqnarray*}
\left\vert A\right\vert &=&\begin{vmatrix}
2 & 5 \\
1 & 3\end{vmatrix}
\\
&=&2\cdot 3-5\cdot 1 \\
&=&1 \\
&\not=&0
\end{eqnarray*}であるため、先の命題より\(A\)は正則行列です。
余因子行列
行列式を用いて正方行列が正則であることを判定する方法が明らかになりました。では、正方行列が正則であることが判明した場合、行列式を用いて逆行列を具体的に求めることはできるのでしょうか。まずは必要な概念を整備します。
次数\(n\)の正方行列\(A=\left(a_{ij}\right) \in M_{n,n}\)が与えられたとき、それぞれの成分\(a_{ij}\ \left( i,j=1,\cdots ,n\right) \)の余因子は、\begin{equation*}A_{ij}=\left( -1\right) ^{i+j}\left\vert M_{ij}\right\vert
\end{equation*}と定義されます。ただし、\(M_{ij}\)は\(A\)から第\(i\)行と第\(j\)列を削除することにより得られる次数\(n-1\)の正方行列です。\(A\)のそれぞれの成分\(a_{ij}\)に対して余因子\(A_{ij}\)を計算することにより、\(A_{ij}\)を\(ij\)成分とする次数\(n\)の正方行列\begin{equation*}\left( A_{ij}\right) =\begin{pmatrix}
A_{11} & A_{11} & \cdots & A_{1n} \\
A_{21} & A_{22} & \cdots & A_{2n} \\
\cdots & \cdots & \cdots & \cdots \\
A_{n1} & A_{n2} & \cdots & A_{nn}\end{pmatrix}\end{equation*}が得られます。さらにこの正方行列の転置行列を、\begin{equation*}
\mathrm{adj}\left( A\right) =\left( A_{ij}\right)
\end{equation*}で表記し、これを\(A\)の余因子行列(cofactor matrix)と呼びます。つまり、\begin{equation*}\mathrm{adj}\left( A\right) =\begin{pmatrix}
A_{11} & A_{21} & \cdots & A_{n1} \\
A_{12} & A_{22} & \cdots & A_{n2} \\
\cdots & \cdots & \cdots & \cdots \\
A_{1n} & A_{2n} & \cdots & A_{nn}\end{pmatrix}\end{equation*}です。
a_{11} & a_{12} \\
a_{21} & a_{22}\end{pmatrix}\end{equation*}のそれぞれの成分の余因子は、\begin{eqnarray*}
A_{11} &=&\left( -1\right) ^{1+1}\left\vert M_{11}\right\vert =1\cdot
a_{22}=a_{22} \\
A_{12} &=&\left( -1\right) ^{1+2}\left\vert M_{12}\right\vert =\left(
-1\right) a_{21}=-a_{21} \\
A_{21} &=&\left( -1\right) ^{2+1}\left\vert M_{21}\right\vert =\left(
-1\right) a_{12}=-a_{12} \\
A_{22} &=&\left( -1\right) ^{2+2}\left\vert M_{22}\right\vert =1\cdot
a_{11}=a_{11}
\end{eqnarray*}であるため、\(A\)の余因子行列は、\begin{eqnarray*}\mathrm{adj}A &=&\begin{pmatrix}
A_{11} & A_{21} \\
A_{12} & A_{22}\end{pmatrix}
\\
&=&\begin{pmatrix}
a_{22} & -a_{12} \\
-a_{21} & a_{11}\end{pmatrix}\end{eqnarray*}です。
a_{11} & a_{12} & a_{13} \\
a_{21} & a_{22} & a_{23} \\
a_{31} & a_{32} & a_{33}\end{pmatrix}\end{equation*}のそれぞれの成分の余因子は、\begin{eqnarray*}
A_{11} &=&\left( -1\right) ^{1+1}\left\vert M_{11}\right\vert =\begin{vmatrix}
a_{22} & a_{23} \\
a_{32} & a_{33}\end{vmatrix}=a_{22}a_{33}-a_{23}a_{32} \\
A_{12} &=&\left( -1\right) ^{1+2}\left\vert M_{12}\right\vert =-\begin{vmatrix}
a_{21} & a_{23} \\
a_{31} & a_{33}\end{vmatrix}=a_{23}a_{31}-a_{21}a_{33} \\
A_{13} &=&\left( -1\right) ^{1+3}\left\vert M_{13}\right\vert =\begin{vmatrix}
a_{21} & a_{22} \\
a_{31} & a_{32}\end{vmatrix}=a_{21}a_{32}-a_{22}a_{31} \\
A_{21} &=&\left( -1\right) ^{2+1}\left\vert M_{21}\right\vert =-\begin{vmatrix}
a_{12} & a_{13} \\
a_{32} & a_{33}\end{vmatrix}=a_{13}a_{32}-a_{12}a_{33} \\
A_{22} &=&\left( -1\right) ^{2+2}\left\vert M_{22}\right\vert =\begin{vmatrix}
a_{11} & a_{13} \\
a_{31} & a_{33}\end{vmatrix}=a_{11}a_{33}-a_{13}a_{31} \\
A_{23} &=&\left( -1\right) ^{2+3}\left\vert M_{23}\right\vert =-\begin{vmatrix}
a_{11} & a_{12} \\
a_{31} & a_{32}\end{vmatrix}=a_{12}a_{31}-a_{11}a_{32} \\
A_{31} &=&\left( -1\right) ^{3+1}\left\vert M_{31}\right\vert =\begin{vmatrix}
a_{12} & a_{13} \\
a_{22} & a_{23}\end{vmatrix}=a_{12}a_{23}-a_{13}a_{22} \\
A_{32} &=&\left( -1\right) ^{3+2}\left\vert M_{32}\right\vert =-\begin{vmatrix}
a_{11} & a_{13} \\
a_{21} & a_{23}\end{vmatrix}=a_{13}a_{21}-a_{11}a_{23} \\
A_{33} &=&\left( -1\right) ^{3+3}\left\vert M_{33}\right\vert =\begin{vmatrix}
a_{11} & a_{12} \\
a_{21} & a_{22}\end{vmatrix}=a_{11}a_{22}-a_{12}a_{21}
\end{eqnarray*}であるため、\(A\)の余因子行列は、\begin{eqnarray*}\mathrm{adj}A &=&\begin{pmatrix}
A_{11} & A_{21} & A_{31} \\
A_{12} & A_{22} & A_{32} \\
A_{13} & A_{23} & A_{33}\end{pmatrix}
\\
&=&\begin{pmatrix}
a_{22}a_{33}-a_{23}a_{32} & a_{23}a_{31}-a_{21}a_{33} &
a_{12}a_{23}-a_{13}a_{22} \\
a_{13}a_{32}-a_{12}a_{33} & a_{11}a_{33}-a_{13}a_{31} &
a_{12}a_{31}-a_{11}a_{32} \\
a_{21}a_{32}-a_{22}a_{31} & a_{13}a_{21}-a_{11}a_{23} &
a_{11}a_{22}-a_{12}a_{21}\end{pmatrix}\end{eqnarray*}です。
余因子行列を用いた逆行列の導出
正方行列が正則である場合、その逆行列は余因子行列を用いて以下のように導出可能です。
\end{equation*}が成り立つことと、\(A\)が正則行列であることは必要十分である。さらにこのとき、\begin{equation*}A^{-1}=\frac{\mathrm{adj}\left( A\right) }{\left\vert A\right\vert }
\end{equation*}が成り立つ。ただし、\(\left\vert A\right\vert \)は\(A\)の行列式であり、\(\mathrm{adj}\left( A\right) \)は\(A\)の余因子行列である。
a_{11} & a_{12} \\
a_{21} & a_{22}\end{pmatrix}\end{equation*}の行列式は、\begin{equation*}
\left\vert A\right\vert =a_{11}a_{22}-a_{12}a_{21}
\end{equation*}であるため、\begin{equation*}
\left\vert A\right\vert \not=0
\end{equation*}すなわち、\begin{equation*}
a_{11}a_{22}-a_{12}a_{21}\not=0
\end{equation*}である場合には\(A\)は正則です。さらに、先に明らかにしたように\(A\)の余因子行列は、\begin{eqnarray*}\mathrm{adj}A &=&\begin{pmatrix}
A_{11} & A_{21} \\
A_{12} & A_{22}\end{pmatrix}
\\
&=&\begin{pmatrix}
a_{22} & -a_{12} \\
-a_{21} & a_{11}\end{pmatrix}\end{eqnarray*}であるため、先の命題より、\(A\)の逆行列は、\begin{eqnarray*}A^{-1} &=&\frac{\mathrm{adj}\left( A\right) }{\left\vert A\right\vert } \\
&=&\frac{1}{a_{11}a_{22}-a_{12}a_{21}}\begin{pmatrix}
a_{22} & -a_{12} \\
-a_{21} & a_{11}\end{pmatrix}
\\
&=&\begin{pmatrix}
\dfrac{a_{22}}{a_{11}a_{22}-a_{12}a_{21}} & -\dfrac{a_{12}}{a_{11}a_{22}-a_{12}a_{21}} \\
-\dfrac{a_{21}}{a_{11}a_{22}-a_{12}a_{21}} & \dfrac{a_{11}}{a_{11}a_{22}-a_{12}a_{21}}\end{pmatrix}\end{eqnarray*}となります。
a_{11} & a_{12} & a_{13} \\
a_{21} & a_{22} & a_{23} \\
a_{31} & a_{32} & a_{33}\end{pmatrix}\end{equation*}の行列式は、\begin{equation*}
\left\vert A\right\vert
=a_{11}a_{22}a_{33}+a_{12}a_{21}a_{31}+a_{13}a_{21}a_{32}-a_{11}a_{23}a_{32}-a_{12}a_{21}a_{33}-a_{13}a_{22}a_{31}
\end{equation*}であるため、\begin{equation*}
\left\vert A\right\vert \not=0
\end{equation*}すなわち、\begin{equation*}
a_{11}a_{22}a_{33}+a_{12}a_{21}a_{31}+a_{13}a_{21}a_{32}-a_{11}a_{23}a_{32}-a_{12}a_{21}a_{33}-a_{13}a_{22}a_{31}\not=0
\end{equation*}である場合には\(A\)は正則です。さらに、先に明らかにしたように\(A\)の余因子行列は、\(A\)の余因子行列は、\begin{eqnarray*}\mathrm{adj}A &=&\begin{pmatrix}
A_{11} & A_{21} & A_{31} \\
A_{12} & A_{22} & A_{32} \\
A_{13} & A_{23} & A_{33}\end{pmatrix}
\\
&=&\begin{pmatrix}
a_{22}a_{33}-a_{23}a_{32} & a_{23}a_{31}-a_{21}a_{33} &
a_{12}a_{23}-a_{13}a_{22} \\
a_{13}a_{32}-a_{12}a_{33} & a_{11}a_{33}-a_{13}a_{31} &
a_{12}a_{31}-a_{11}a_{32} \\
a_{21}a_{32}-a_{22}a_{31} & a_{13}a_{21}-a_{11}a_{23} &
a_{11}a_{22}-a_{12}a_{21}\end{pmatrix}\end{eqnarray*}であるため、先の命題より、\(A\)の逆行列は、\begin{eqnarray*}A^{-1} &=&\frac{\mathrm{adj}\left( A\right) }{\left\vert A\right\vert } \\
&=&\frac{1}{\left\vert A\right\vert }\begin{pmatrix}
a_{22}a_{33}-a_{23}a_{32} & a_{23}a_{31}-a_{21}a_{33} &
a_{12}a_{23}-a_{13}a_{22} \\
a_{13}a_{32}-a_{12}a_{33} & a_{11}a_{33}-a_{13}a_{31} &
a_{12}a_{31}-a_{11}a_{32} \\
a_{21}a_{32}-a_{22}a_{31} & a_{13}a_{21}-a_{11}a_{23} &
a_{11}a_{22}-a_{12}a_{21}\end{pmatrix}
\\
&=&\begin{pmatrix}
\dfrac{a_{22}a_{33}-a_{23}a_{32}}{\left\vert A\right\vert } & \dfrac{a_{23}a_{31}-a_{21}a_{33}}{\left\vert A\right\vert } & \dfrac{a_{12}a_{23}-a_{13}a_{22}}{\left\vert A\right\vert } \\
\dfrac{a_{13}a_{32}-a_{12}a_{33}}{\left\vert A\right\vert } & \dfrac{a_{11}a_{33}-a_{13}a_{31}}{\left\vert A\right\vert } & \dfrac{a_{12}a_{31}-a_{11}a_{32}}{\left\vert A\right\vert } \\
\dfrac{a_{21}a_{32}-a_{22}a_{31}}{\left\vert A\right\vert } & \dfrac{a_{13}a_{21}-a_{11}a_{23}}{\left\vert A\right\vert } & \dfrac{a_{11}a_{22}-a_{12}a_{21}}{\left\vert A\right\vert }\end{pmatrix}\end{eqnarray*}となります。ただし、\begin{equation*}
\left\vert A\right\vert
=a_{11}a_{22}a_{33}+a_{12}a_{21}a_{31}+a_{13}a_{21}a_{32}-a_{11}a_{23}a_{32}-a_{12}a_{21}a_{33}-a_{13}a_{22}a_{31}
\end{equation*}です。
演習問題
A=\begin{pmatrix}
1 & -2 \\
3 & 5\end{pmatrix}\end{equation*}が正則であることを示した上で、その逆行列\(A^{-1}\)を求めてください。
A=\begin{pmatrix}
3 & -1 & 2 \\
2 & 3 & -3 \\
3 & 1 & 3\end{pmatrix}\end{equation*}が正則であることを示した上で、その逆行列\(A^{-1}\)を求めてください。
プレミアム会員専用コンテンツです
【ログイン】【会員登録】