この章の目標
行列の和・スカラー倍・積を正確に計算でき、行列の積が「写像の合成」を表すことを理解する
行基本変形と基本行列の関係を理解し、任意の行列を簡約階段行列に変形できる(簡約階段行列の一意性の証明を含む)
階数を用いて連立一次方程式の解の存在を判定し、解の全体を「特殊解+同次方程式の解」として記述できる
正則行列の同値条件を証明し、逆行列を行基本変形で計算できる
ブロック行列の積を使いこなせる
前提 :数学の言葉 第2章 集合 、第3章 写像
1.1 連立一次方程式と行列
線形代数学の出発点は、次のような連立一次方程式 (system of linear equations) である。
{ x 1 + 2 x 2 + x 3 = 4 2 x 1 + 4 x 2 + x 3 + x 4 = 7 3 x 1 + 6 x 2 − x 3 + 4 x 4 = 8 \left\lbrace
\begin{aligned}
x_1 + 2x_2 + x_3 &= 4 \\
2x_1 + 4x_2 + x_3 + x_4 &= 7 \\
3x_1 + 6x_2 - x_3 + 4x_4 &= 8
\end{aligned}
\right. ⎩ ⎨ ⎧ x 1 + 2 x 2 + x 3 2 x 1 + 4 x 2 + x 3 + x 4 3 x 1 + 6 x 2 − x 3 + 4 x 4 = 4 = 7 = 8
未知数が 4 個、方程式が 3 個ある。高校までは 2〜3 元の連立方程式を代入法や加減法で解いたが、未知数や方程式の個数が増えると、手順を系統立てなければ見通しが立たない。本章の目的は、次の二つの問いに完全に答えることである。
解はいつ存在するか。
解が存在するとき、解の全体はどのような形をしているか。
鍵になるのは、係数だけを長方形に並べた行列 に注目することである。未知数の文字 x 1 , x 2 , … x_1, x_2, \dots x 1 , x 2 , … は計算の途中で何の役割も果たさず、係数の並び方だけが本質だからである。
本章では、数の範囲 K K K は R \mathbb{R} R または C \mathbb{C} C とする。実際には四則演算が自由にできる数の体系(体 、第2章 で定義する)であれば、本章の議論はすべてそのまま成り立つ。K K K の元をスカラー (scalar) という。
定義 1.1 (行列, matrix)m , n ∈ N m, n \in \mathbb{N} m , n ∈ N とする。K K K の元 a i j a_{ij} a ij ( 1 ≤ i ≤ m , 1 ≤ j ≤ n ) (1 \leq i \leq m,\ 1 \leq j \leq n) ( 1 ≤ i ≤ m , 1 ≤ j ≤ n ) を縦 m m m 行・横 n n n 列の長方形に並べたもの
A = ( a 11 a 12 ⋯ a 1 n a 21 a 22 ⋯ a 2 n ⋮ ⋮ ⋱ ⋮ a m 1 a m 2 ⋯ a m n ) A = \begin{pmatrix} a_{11} & a_{12} & \cdots & a_{1n} \\ a_{21} & a_{22} & \cdots & a_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ a_{m1} & a_{m2} & \cdots & a_{mn} \end{pmatrix} A = a 11 a 21 ⋮ a m 1 a 12 a 22 ⋮ a m 2 ⋯ ⋯ ⋱ ⋯ a 1 n a 2 n ⋮ a mn
を m × n m \times n m × n 行列 (m m m 行 n n n 列の行列)という。a i j a_{ij} a ij を A A A の ( i , j ) (i, j) ( i , j ) 成分 (entry) といい、A = ( a i j ) A = (a_{ij}) A = ( a ij ) と略記する。m × n m \times n m × n 行列全体を M m , n ( K ) \operatorname{M}_{m,n}(K) M m , n ( K ) 、n × n n \times n n × n 行列(n n n 次正方行列 , square matrix)全体を M n ( K ) \operatorname{M}_n(K) M n ( K ) と書く。
横の並び ( a i 1 , … , a i n ) (a_{i1}, \dots, a_{in}) ( a i 1 , … , a in ) を第 i i i 行 (row)、縦の並びを第 j j j 列 (column) という。m × 1 m \times 1 m × 1 行列を列ベクトル 、1 × n 1 \times n 1 × n 行列を行ベクトル という。n n n 成分の列ベクトル全体を K n K^n K n と書き、数ベクトル空間 と呼ぶ。紙面の節約のため、列ベクトルを t ( x 1 , … , x n ) {}^t(x_1, \dots, x_n) t ( x 1 , … , x n ) と横に書くことがある(t {}^t t は後で定義する転置)。
特別な行列に名前をつけておく。すべての成分が 0 の行列を零行列 O O O (サイズを明示するときは O m , n O_{m,n} O m , n )という。n n n 次正方行列で対角成分が 1、それ以外が 0 のものを単位行列 (identity matrix) I n I_n I n という。クロネッカーのデルタ δ i j \delta_{ij} δ ij (i = j i = j i = j なら 1、そうでなければ 0)を使えば I n = ( δ i j ) I_n = (\delta_{ij}) I n = ( δ ij ) である。I n I_n I n の第 j j j 列を e j e_j e j と書き、基本ベクトル (標準単位ベクトル)という。また ( i , j ) (i, j) ( i , j ) 成分だけが 1 で他が 0 の行列を行列単位 E i j E_{ij} E ij という。
例 1.2 冒頭の連立方程式の係数を並べた行列 A A A 、右辺を並べた列ベクトル b b b 、それらを横に並べた行列 ( A ∣ b ) (A \mid b) ( A ∣ b ) はそれぞれ
A = ( 1 2 1 0 2 4 1 1 3 6 − 1 4 ) , b = ( 4 7 8 ) , ( A ∣ b ) = ( 1 2 1 0 4 2 4 1 1 7 3 6 − 1 4 8 ) A = \begin{pmatrix} 1 & 2 & 1 & 0 \\ 2 & 4 & 1 & 1 \\ 3 & 6 & -1 & 4 \end{pmatrix}, \quad
b = \begin{pmatrix} 4 \\ 7 \\ 8 \end{pmatrix}, \quad
(A \mid b) = \left(\begin{array}{cccc|c} 1 & 2 & 1 & 0 & 4 \\ 2 & 4 & 1 & 1 & 7 \\ 3 & 6 & -1 & 4 & 8 \end{array}\right) A = 1 2 3 2 4 6 1 1 − 1 0 1 4 , b = 4 7 8 , ( A ∣ b ) = 1 2 3 2 4 6 1 1 − 1 0 1 4 4 7 8
である。A A A を係数行列 (coefficient matrix)、( A ∣ b ) (A \mid b) ( A ∣ b ) を拡大係数行列 (augmented matrix) という。次節で定義する行列の積を用いると、連立方程式は A x = b Ax = b A x = b (x = t ( x 1 , x 2 , x 3 , x 4 ) x = {}^t(x_1, x_2, x_3, x_4) x = t ( x 1 , x 2 , x 3 , x 4 ) )と一行で書ける。
ここで、本教材全体を貫く視点を先に述べておく。m × n m \times n m × n 行列 A A A は、列ベクトル x ∈ K n x \in K^n x ∈ K n に列ベクトル A x ∈ K m Ax \in K^m A x ∈ K m を対応させる写像を定める。連立方程式 A x = b Ax = b A x = b を解くとは、この写像による b b b の逆像を求めることにほかならない。行列は線形写像を座標で表示したものである ——この見方は第3章 で完全な形になり、以後のすべての章の土台になる。
1.2 行列の演算
和・スカラー倍・積
定義 1.3 (行列の演算)
A = ( a i j ) , B = ( b i j ) ∈ M m , n ( K ) A = (a_{ij}), B = (b_{ij}) \in \operatorname{M}_{m,n}(K) A = ( a ij ) , B = ( b ij ) ∈ M m , n ( K ) と c ∈ K c \in K c ∈ K に対し、和 A + B A + B A + B とスカラー倍 c A cA c A を成分ごとに ( A + B ) i j = a i j + b i j (A + B)_{ij} = a_{ij} + b_{ij} ( A + B ) ij = a ij + b ij 、( c A ) i j = c a i j (cA)_{ij} = c a_{ij} ( c A ) ij = c a ij で定める。
A = ( a i j ) ∈ M m , n ( K ) A = (a_{ij}) \in \operatorname{M}_{m,n}(K) A = ( a ij ) ∈ M m , n ( K ) 、B = ( b j k ) ∈ M n , p ( K ) B = (b_{jk}) \in \operatorname{M}_{n,p}(K) B = ( b j k ) ∈ M n , p ( K ) に対し、積 A B ∈ M m , p ( K ) AB \in \operatorname{M}_{m,p}(K) A B ∈ M m , p ( K ) を
( A B ) i k = ∑ j = 1 n a i j b j k = a i 1 b 1 k + a i 2 b 2 k + ⋯ + a i n b n k (AB)_{ik} = \sum_{j=1}^{n} a_{ij} b_{jk} = a_{i1} b_{1k} + a_{i2} b_{2k} + \cdots + a_{in} b_{nk} ( A B ) ik = j = 1 ∑ n a ij b j k = a i 1 b 1 k + a i 2 b 2 k + ⋯ + a in b nk
で定める。積は「A A A の列の数 = = = B B B の行の数」のときにだけ定義される。
積の ( i , k ) (i, k) ( i , k ) 成分は「A A A の第 i i i 行と B B B の第 k k k 列の対応する成分を掛けて足したもの」である。この一見人工的な定義の理由は、後で見るように「写像の合成」を表すからである。
命題 1.4 (演算法則)以下、積が定義されるサイズの行列について、次が成り立つ。
( A + B ) + C = A + ( B + C ) (A + B) + C = A + (B + C) ( A + B ) + C = A + ( B + C ) 、A + B = B + A A + B = B + A A + B = B + A 、A + O = A A + O = A A + O = A 、A + ( − 1 ) A = O A + (-1)A = O A + ( − 1 ) A = O
(結合法則)( A B ) C = A ( B C ) (AB)C = A(BC) ( A B ) C = A ( B C )
(分配法則)A ( B + C ) = A B + A C A(B + C) = AB + AC A ( B + C ) = A B + A C 、( A + B ) C = A C + B C (A + B)C = AC + BC ( A + B ) C = A C + B C
c ( A B ) = ( c A ) B = A ( c B ) c(AB) = (cA)B = A(cB) c ( A B ) = ( c A ) B = A ( c B )
A ∈ M m , n ( K ) A \in \operatorname{M}_{m,n}(K) A ∈ M m , n ( K ) に対し I m A = A = A I n I_m A = A = A I_n I m A = A = A I n
証明. 1, 3, 4, 5 は成分を書けば直ちにわかる。例えば 3 の前半は ∑ j a i j ( b j k + c j k ) = ∑ j a i j b j k + ∑ j a i j c j k \sum_j a_{ij}(b_{jk} + c_{jk}) = \sum_j a_{ij} b_{jk} + \sum_j a_{ij} c_{jk} ∑ j a ij ( b j k + c j k ) = ∑ j a ij b j k + ∑ j a ij c j k から、5 は ∑ j δ i j a j k = a i k \sum_j \delta_{ij} a_{jk} = a_{ik} ∑ j δ ij a j k = a ik から従う。2 を示す。A ∈ M m , n ( K ) A \in \operatorname{M}_{m,n}(K) A ∈ M m , n ( K ) 、B ∈ M n , p ( K ) B \in \operatorname{M}_{n,p}(K) B ∈ M n , p ( K ) 、C ∈ M p , q ( K ) C \in \operatorname{M}_{p,q}(K) C ∈ M p , q ( K ) とすると、有限和の順序を交換して
( ( A B ) C ) i l = ∑ k = 1 p ( ∑ j = 1 n a i j b j k ) c k l = ∑ j = 1 n a i j ( ∑ k = 1 p b j k c k l ) = ( A ( B C ) ) i l ((AB)C)_{il} = \sum_{k=1}^{p} \left( \sum_{j=1}^{n} a_{ij} b_{jk} \right) c_{kl} = \sum_{j=1}^{n} a_{ij} \left( \sum_{k=1}^{p} b_{jk} c_{kl} \right) = (A(BC))_{il} (( A B ) C ) i l = k = 1 ∑ p ( j = 1 ∑ n a ij b j k ) c k l = j = 1 ∑ n a ij ( k = 1 ∑ p b j k c k l ) = ( A ( B C ) ) i l
となる。□ \square □
結合法則により、A B C ABC A B C のように括弧を省いて書いてよい。一方、積は可換ではない 。
例 1.5 次の 2 次正方行列を考える。
A = ( 0 1 0 0 ) , B = ( 0 0 1 0 ) , A B = ( 1 0 0 0 ) , B A = ( 0 0 0 1 ) A = \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix}, \quad B = \begin{pmatrix} 0 & 0 \\ 1 & 0 \end{pmatrix}, \quad
AB = \begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix}, \quad BA = \begin{pmatrix} 0 & 0 \\ 0 & 1 \end{pmatrix} A = ( 0 0 1 0 ) , B = ( 0 1 0 0 ) , A B = ( 1 0 0 0 ) , B A = ( 0 0 0 1 )
であるから A B ≠ B A AB \neq BA A B = B A である。さらに A 2 = O A^2 = O A 2 = O である。つまり A ≠ O A \neq O A = O なのに A 2 = O A^2 = O A 2 = O となる(このような行列を冪零行列 , nilpotent matrix という)。また A E 11 = O A E_{11} = O A E 11 = O も直接計算で確かめられる。このように、行列の世界には「X Y = O XY = O X Y = O だが X ≠ O X \neq O X = O かつ Y ≠ O Y \neq O Y = O 」という零因子 (zero divisor) がある。
注意
行列の計算では、数の計算の習慣が通用しない場面がある。A B = A C AB = AC A B = A C から B = C B = C B = C は一般に導けない(A A A が正則なら導ける)。( A + B ) 2 = A 2 + A B + B A + B 2 (A + B)^2 = A^2 + AB + BA + B^2 ( A + B ) 2 = A 2 + A B + B A + B 2 であり、A B = B A AB = BA A B = B A でない限り A 2 + 2 A B + B 2 A^2 + 2AB + B^2 A 2 + 2 A B + B 2 にはならない。
積の列ごとの見方
積の定義を列ごとに読み直すと、非常に役に立つ見方が得られる。
命題 1.6 A ∈ M m , n ( K ) A \in \operatorname{M}_{m,n}(K) A ∈ M m , n ( K ) の列を a 1 , … , a n ∈ K m a_1, \dots, a_n \in K^m a 1 , … , a n ∈ K m とし、A = ( a 1 a 2 ⋯ a n ) A = (a_1 \ a_2 \ \cdots \ a_n) A = ( a 1 a 2 ⋯ a n ) と書く。
x = t ( x 1 , … , x n ) ∈ K n x = {}^t(x_1, \dots, x_n) \in K^n x = t ( x 1 , … , x n ) ∈ K n に対し、A x = x 1 a 1 + x 2 a 2 + ⋯ + x n a n Ax = x_1 a_1 + x_2 a_2 + \cdots + x_n a_n A x = x 1 a 1 + x 2 a 2 + ⋯ + x n a n 。特に A e j = a j A e_j = a_j A e j = a j 。
B ∈ M n , p ( K ) B \in \operatorname{M}_{n,p}(K) B ∈ M n , p ( K ) の列を b 1 , … , b p b_1, \dots, b_p b 1 , … , b p とすると、A B = ( A b 1 A b 2 ⋯ A b p ) AB = (Ab_1 \ Ab_2 \ \cdots \ Ab_p) A B = ( A b 1 A b 2 ⋯ A b p ) 。
A B AB A B の第 i i i 行は、(A A A の第 i i i 行)⋅ B \cdot B ⋅ B に等しい。
証明. 1:A x Ax A x の第 i i i 成分は ∑ j a i j x j \sum_j a_{ij} x_j ∑ j a ij x j であり、これは ∑ j x j a j \sum_j x_j a_j ∑ j x j a j の第 i i i 成分に等しい。2:A B AB A B の第 k k k 列の第 i i i 成分は ∑ j a i j b j k \sum_j a_{ij} b_{jk} ∑ j a ij b j k で、これは A b k A b_k A b k の第 i i i 成分である。3 も同様である。□ \square □
命題 1.6 の 1 は、「A x = b Ax = b A x = b が解をもつ」ことが「b b b が A A A の列の一次結合 (定数倍の和)で書ける」ことと同じだと言っている。この言い換えは第2章 以降で繰り返し使う。
行列の積と写像の合成
定義 1.7 A ∈ M m , n ( K ) A \in \operatorname{M}_{m,n}(K) A ∈ M m , n ( K ) に対し、写像 L A : K n → K m L_A\colon K^n \to K^m L A : K n → K m を L A ( x ) = A x L_A(x) = Ax L A ( x ) = A x で定める。
分配法則と命題 1.4 の 4 から、L A L_A L A は
L A ( x + y ) = L A ( x ) + L A ( y ) , L A ( c x ) = c L A ( x ) ( x , y ∈ K n , c ∈ K ) L_A(x + y) = L_A(x) + L_A(y), \qquad L_A(cx) = c L_A(x) \qquad (x, y \in K^n,\ c \in K) L A ( x + y ) = L A ( x ) + L A ( y ) , L A ( c x ) = c L A ( x ) ( x , y ∈ K n , c ∈ K )
を満たす。この二つの性質をもつ写像を線形写像 と呼ぶ(第3章 )。
命題 1.8
A ∈ M m , n ( K ) A \in \operatorname{M}_{m,n}(K) A ∈ M m , n ( K ) 、B ∈ M n , p ( K ) B \in \operatorname{M}_{n,p}(K) B ∈ M n , p ( K ) に対し L A ∘ L B = L A B L_A \circ L_B = L_{AB} L A ∘ L B = L A B 。
A , A ′ ∈ M m , n ( K ) A, A' \in \operatorname{M}_{m,n}(K) A , A ′ ∈ M m , n ( K ) について、L A = L A ′ L_A = L_{A'} L A = L A ′ ならば A = A ′ A = A' A = A ′ 。
証明. 1:x ∈ K p x \in K^p x ∈ K p に対し ( L A ∘ L B ) ( x ) = A ( B x ) = ( A B ) x = L A B ( x ) (L_A \circ L_B)(x) = A(Bx) = (AB)x = L_{AB}(x) ( L A ∘ L B ) ( x ) = A ( B x ) = ( A B ) x = L A B ( x ) (結合法則)。2:命題 1.6 より A A A の第 j j j 列は L A ( e j ) L_A(e_j) L A ( e j ) であるから、L A = L A ′ L_A = L_{A'} L A = L A ′ なら A A A と A ′ A' A ′ の列はすべて一致する。□ \square □
つまり行列 A A A と写像 L A L_A L A は一対一に対応し、行列の積は写像の合成に対応する。行列の積の定義は「合成がこうなるように」決められているのである。実は結合法則もこの観点から説明できる:写像の合成は常に結合的であり、L ( A B ) C = ( L A ∘ L B ) ∘ L C = L A ∘ ( L B ∘ L C ) = L A ( B C ) L_{(AB)C} = (L_A \circ L_B) \circ L_C = L_A \circ (L_B \circ L_C) = L_{A(BC)} L ( A B ) C = ( L A ∘ L B ) ∘ L C = L A ∘ ( L B ∘ L C ) = L A ( B C ) から命題 1.8 の 2 により ( A B ) C = A ( B C ) (AB)C = A(BC) ( A B ) C = A ( B C ) が従う(ただしこの議論は 1 の証明で結合法則 A ( B x ) = ( A B ) x A(Bx) = (AB)x A ( B x ) = ( A B ) x を使っているので、C C C が列ベクトルの場合の結合法則から一般の場合を導いたことになる)。
例 1.9 (回転行列)角 θ \theta θ の回転を表す行列
R ( θ ) = ( cos θ − sin θ sin θ cos θ ) R(\theta) = \begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix} R ( θ ) = ( cos θ sin θ − sin θ cos θ )
を考える。R ( θ ) e 1 = t ( cos θ , sin θ ) R(\theta) e_1 = {}^t(\cos\theta, \sin\theta) R ( θ ) e 1 = t ( cos θ , sin θ ) 、R ( θ ) e 2 = t ( − sin θ , cos θ ) R(\theta) e_2 = {}^t(-\sin\theta, \cos\theta) R ( θ ) e 2 = t ( − sin θ , cos θ ) であり、L R ( θ ) L_{R(\theta)} L R ( θ ) は平面の原点のまわりの角 θ \theta θ の回転である。「角 β \beta β 回してから角 α \alpha α 回す」ことは「角 α + β \alpha + \beta α + β 回す」ことだから、命題 1.8 より R ( α ) R ( β ) = R ( α + β ) R(\alpha) R(\beta) = R(\alpha + \beta) R ( α ) R ( β ) = R ( α + β ) が期待される。実際に計算すると
R ( α ) R ( β ) = ( cos α cos β − sin α sin β − ( cos α sin β + sin α cos β ) sin α cos β + cos α sin β cos α cos β − sin α sin β ) R(\alpha) R(\beta) = \begin{pmatrix} \cos\alpha\cos\beta - \sin\alpha\sin\beta & -(\cos\alpha\sin\beta + \sin\alpha\cos\beta) \\ \sin\alpha\cos\beta + \cos\alpha\sin\beta & \cos\alpha\cos\beta - \sin\alpha\sin\beta \end{pmatrix} R ( α ) R ( β ) = ( cos α cos β − sin α sin β sin α cos β + cos α sin β − ( cos α sin β + sin α cos β ) cos α cos β − sin α sin β )
であり、これが R ( α + β ) R(\alpha + \beta) R ( α + β ) に等しいことは三角関数の加法定理そのものである。逆に言えば、加法定理は「回転の合成は回転」という幾何学的事実の座標表示である。
転置
定義 1.11 (転置, transpose)A = ( a i j ) ∈ M m , n ( K ) A = (a_{ij}) \in \operatorname{M}_{m,n}(K) A = ( a ij ) ∈ M m , n ( K ) に対し、( j , i ) (j, i) ( j , i ) 成分が a i j a_{ij} a ij である n × m n \times m n × m 行列を A A A の転置行列 といい、t A {}^t A t A と書く(A T A^{\mathsf{T}} A T と書く本も多い)。t A = A {}^t A = A t A = A を満たす正方行列を対称行列 (symmetric matrix)、t A = − A {}^t A = -A t A = − A を満たすものを交代行列 (歪対称行列, skew-symmetric matrix)という。
命題 1.12 t ( t A ) = A {}^t({}^t A) = A t ( t A ) = A 、t ( A + B ) = t A + t B {}^t(A + B) = {}^t A + {}^t B t ( A + B ) = t A + t B 、t ( c A ) = c t A {}^t(cA) = c\ {}^t A t ( c A ) = c t A 、t ( A B ) = t B t A {}^t(AB) = {}^t B\ {}^t A t ( A B ) = t B t A 。
証明. 最後の式のみ示す。( t ( A B ) ) k i = ( A B ) i k = ∑ j a i j b j k = ∑ j ( t B ) k j ( t A ) j i = ( t B t A ) k i ({}^t(AB))_{ki} = (AB)_{ik} = \sum_j a_{ij} b_{jk} = \sum_j ({}^t B)_{kj} ({}^t A)_{ji} = ({}^t B\ {}^t A)_{ki} ( t ( A B ) ) k i = ( A B ) ik = ∑ j a ij b j k = ∑ j ( t B ) k j ( t A ) j i = ( t B t A ) k i 。□ \square □
順序が入れ替わることに注意する。これは「靴下を履いてから靴を履く」の逆操作が「靴を脱いでから靴下を脱ぐ」であるのと同じ構造で、逆行列についても同様の現象が起こる。
正則行列
定義 1.13 (正則行列, 逆行列)A ∈ M n ( K ) A \in \operatorname{M}_n(K) A ∈ M n ( K ) に対し、A X = X A = I n AX = XA = I_n A X = X A = I n を満たす X ∈ M n ( K ) X \in \operatorname{M}_n(K) X ∈ M n ( K ) が存在するとき、A A A は正則 (invertible, regular, nonsingular) であるといい、X X X を A A A の逆行列 (inverse matrix) という。n n n 次正則行列全体を GL n ( K ) \operatorname{GL}_n(K) GL n ( K ) と書く。
命題 1.14
逆行列は存在すればただ一つである。これを A − 1 A^{-1} A − 1 と書く。
A , B ∈ GL n ( K ) A, B \in \operatorname{GL}_n(K) A , B ∈ GL n ( K ) ならば A B ∈ GL n ( K ) AB \in \operatorname{GL}_n(K) A B ∈ GL n ( K ) で、( A B ) − 1 = B − 1 A − 1 (AB)^{-1} = B^{-1} A^{-1} ( A B ) − 1 = B − 1 A − 1 。
A ∈ GL n ( K ) A \in \operatorname{GL}_n(K) A ∈ GL n ( K ) ならば A − 1 , t A ∈ GL n ( K ) A^{-1}, {}^t A \in \operatorname{GL}_n(K) A − 1 , t A ∈ GL n ( K ) で、( A − 1 ) − 1 = A (A^{-1})^{-1} = A ( A − 1 ) − 1 = A 、( t A ) − 1 = t ( A − 1 ) ({}^t A)^{-1} = {}^t(A^{-1}) ( t A ) − 1 = t ( A − 1 ) 。
証明. 1:X , Y X, Y X , Y がともに逆行列なら X = X I n = X ( A Y ) = ( X A ) Y = I n Y = Y X = X I_n = X(AY) = (XA)Y = I_n Y = Y X = X I n = X ( A Y ) = ( X A ) Y = I n Y = Y 。2:( A B ) ( B − 1 A − 1 ) = A ( B B − 1 ) A − 1 = A A − 1 = I n (AB)(B^{-1}A^{-1}) = A(BB^{-1})A^{-1} = AA^{-1} = I_n ( A B ) ( B − 1 A − 1 ) = A ( B B − 1 ) A − 1 = A A − 1 = I n 、逆順の積も同様。3:前半は定義の対称性から。後半は命題 1.12 より t A t ( A − 1 ) = t ( A − 1 A ) = t I n = I n {}^t A\ {}^t(A^{-1}) = {}^t(A^{-1} A) = {}^t I_n = I_n t A t ( A − 1 ) = t ( A − 1 A ) = t I n = I n 、逆順も同様。□ \square □
命題 1.14 は GL n ( K ) \operatorname{GL}_n(K) GL n ( K ) が行列の積について群 をなすことを示している(代数学 第2章 )。GL n ( K ) \operatorname{GL}_n(K) GL n ( K ) を一般線形群 (general linear group) という。
1.3 行基本変形と基本行列
連立方程式を解くとき、私たちは「二つの式を入れ替える」「式を 0 でない数で割る」「ある式の定数倍を別の式に加える」という操作を行う。これらの操作は解の集合を変えない(逆操作で元に戻せるからである)。係数だけに注目すると、これは拡大係数行列の行に対する操作になる。
定義 1.15 (行基本変形, elementary row operations)行列に対する次の 3 種類の操作を行基本変形 という。
(I) 第 i i i 行と第 j j j 行を入れ替える ( i ≠ j ) (i \neq j) ( i = j ) 。
(II) 第 i i i 行を c c c 倍する ( c ≠ 0 ) (c \neq 0) ( c = 0 ) 。
(III) 第 i i i 行に第 j j j 行の c c c 倍を加える ( i ≠ j , c ∈ K ) (i \neq j,\ c \in K) ( i = j , c ∈ K ) 。
行を列に置き換えたものを列基本変形 という。
行基本変形は、次の行列を左から掛けることで実現できる。
定義 1.16 (基本行列, elementary matrix)n n n 次正方行列として次を定める。
P ( i , j ) P(i, j) P ( i , j ) :I n I_n I n の第 i i i 行と第 j j j 行を入れ替えた行列 ( i ≠ j ) (i \neq j) ( i = j ) 。
Q ( i ; c ) = I n + ( c − 1 ) E i i Q(i; c) = I_n + (c - 1) E_{ii} Q ( i ; c ) = I n + ( c − 1 ) E ii :I n I_n I n の ( i , i ) (i, i) ( i , i ) 成分を c c c に置き換えた行列 ( c ≠ 0 ) (c \neq 0) ( c = 0 ) 。
R ( i , j ; c ) = I n + c E i j R(i, j; c) = I_n + c E_{ij} R ( i , j ; c ) = I n + c E ij ( i ≠ j ) (i \neq j) ( i = j ) 。
これらを基本行列 という。
例えば n = 3 n = 3 n = 3 のとき
P ( 1 , 2 ) = ( 0 1 0 1 0 0 0 0 1 ) , Q ( 2 ; c ) = ( 1 0 0 0 c 0 0 0 1 ) , R ( 1 , 3 ; c ) = ( 1 0 c 0 1 0 0 0 1 ) P(1, 2) = \begin{pmatrix} 0 & 1 & 0 \\ 1 & 0 & 0 \\ 0 & 0 & 1 \end{pmatrix}, \quad
Q(2; c) = \begin{pmatrix} 1 & 0 & 0 \\ 0 & c & 0 \\ 0 & 0 & 1 \end{pmatrix}, \quad
R(1, 3; c) = \begin{pmatrix} 1 & 0 & c \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix} P ( 1 , 2 ) = 0 1 0 1 0 0 0 0 1 , Q ( 2 ; c ) = 1 0 0 0 c 0 0 0 1 , R ( 1 , 3 ; c ) = 1 0 0 0 1 0 c 0 1
である。
命題 1.17 A ∈ M m , n ( K ) A \in \operatorname{M}_{m,n}(K) A ∈ M m , n ( K ) とし、基本行列は m m m 次とする。
P ( i , j ) A P(i, j) A P ( i , j ) A 、Q ( i ; c ) A Q(i; c) A Q ( i ; c ) A 、R ( i , j ; c ) A R(i, j; c) A R ( i , j ; c ) A は、A A A にそれぞれ行基本変形 (I), (II), (III) を施した行列である。
同様に、n n n 次の基本行列を右から掛けると列基本変形になる:A P ( i , j ) A P(i, j) A P ( i , j ) は第 i i i 列と第 j j j 列の入れ替え、A Q ( i ; c ) A Q(i; c) A Q ( i ; c ) は第 i i i 列の c c c 倍、A R ( i , j ; c ) A R(i, j; c) A R ( i , j ; c ) は第 j j j 列に第 i i i 列の c c c 倍を加えたものである。
証明. 命題 1.6 の 3 より、X A XA X A の第 k k k 行は(X X X の第 k k k 行)⋅ A \cdot A ⋅ A である。E i j A E_{ij} A E ij A は第 i i i 行が A A A の第 j j j 行で、他の行が 0 の行列だから、R ( i , j ; c ) A = A + c E i j A R(i, j; c) A = A + c E_{ij} A R ( i , j ; c ) A = A + c E ij A は A A A の第 i i i 行に第 j j j 行の c c c 倍を加えたものである。(I), (II) も同様に確かめられる。列については、A E i j A E_{ij} A E ij が第 j j j 列が A A A の第 i i i 列で他の列が 0 の行列であることから同様である。□ \square □
命題 1.18 基本行列は正則であり、P ( i , j ) − 1 = P ( i , j ) P(i, j)^{-1} = P(i, j) P ( i , j ) − 1 = P ( i , j ) 、Q ( i ; c ) − 1 = Q ( i ; c − 1 ) Q(i; c)^{-1} = Q(i; c^{-1}) Q ( i ; c ) − 1 = Q ( i ; c − 1 ) 、R ( i , j ; c ) − 1 = R ( i , j ; − c ) R(i, j; c)^{-1} = R(i, j; -c) R ( i , j ; c ) − 1 = R ( i , j ; − c ) 。
証明. 命題 1.17 により、左辺と右辺の積は「行基本変形とその逆の変形を続けて行うこと」を I n I_n I n に施したものであり、I n I_n I n に戻る。例えば R ( i , j ; − c ) R ( i , j ; c ) = ( I + c E i j ) − c E i j ( I + c E i j ) = I − c 2 E i j E i j = I R(i, j; -c) R(i, j; c) = (I + c E_{ij}) - c E_{ij}(I + c E_{ij}) = I - c^2 E_{ij} E_{ij} = I R ( i , j ; − c ) R ( i , j ; c ) = ( I + c E ij ) − c E ij ( I + c E ij ) = I − c 2 E ij E ij = I である(i ≠ j i \neq j i = j なので E i j E i j = O E_{ij} E_{ij} = O E ij E ij = O )。□ \square □
定義 1.19 (行同値)行列 B B B が A A A に有限回の行基本変形を施して得られるとき、A A A と B B B は行同値 (row equivalent) であるという。これは、基本行列の積 P P P によって B = P A B = PA B = P A と書けることと同じである。
命題 1.18 により行基本変形は逆に戻せるので、行同値は同値関係である(数学の言葉 第4章 )。
命題 1.20 拡大係数行列 ( A ∣ b ) (A \mid b) ( A ∣ b ) と ( A ′ ∣ b ′ ) (A' \mid b') ( A ′ ∣ b ′ ) が行同値ならば、連立方程式 A x = b Ax = b A x = b と A ′ x = b ′ A'x = b' A ′ x = b ′ の解の集合は一致する。
証明. ( A ′ ∣ b ′ ) = P ( A ∣ b ) (A' \mid b') = P(A \mid b) ( A ′ ∣ b ′ ) = P ( A ∣ b ) 、すなわち A ′ = P A A' = PA A ′ = P A 、b ′ = P b b' = Pb b ′ = P b となる正則行列 P P P (基本行列の積)がある。A x = b Ax = b A x = b なら A ′ x = P A x = P b = b ′ A'x = PAx = Pb = b' A ′ x = P A x = P b = b ′ 。逆に A ′ x = b ′ A'x = b' A ′ x = b ′ なら A x = P − 1 A ′ x = P − 1 b ′ = b Ax = P^{-1}A'x = P^{-1}b' = b A x = P − 1 A ′ x = P − 1 b ′ = b 。□ \square □
1.4 簡約階段行列
行基本変形を繰り返して、方程式を「一目で解ける形」まで変形したい。その到達点が簡約階段行列である。
定義 1.21 (簡約階段行列, reduced row echelon form)行列 R R R が次の条件を満たすとき、簡約階段行列 であるという。
零でない行(0 でない成分を含む行)は、零行(すべての成分が 0 の行)より上にある。
零でない各行の最も左にある 0 でない成分(主成分 , pivot)は 1 である。
零でない行が r r r 本あり、第 i i i 行の主成分が第 j i j_i j i 列にあるとすると、j 1 < j 2 < ⋯ < j r j_1 < j_2 < \cdots < j_r j 1 < j 2 < ⋯ < j r 。
主成分を含む列(主列 , pivot column)では、主成分以外の成分はすべて 0 である。
条件 1〜3 のみを満たすものを階段行列 (row echelon form) という。
例 1.22 次の行列はいずれも簡約階段行列である。
( 1 2 0 1 0 0 1 − 1 0 0 0 0 ) , ( 1 0 0 0 1 0 0 0 1 ) , ( 0 1 3 0 0 0 0 1 ) \begin{pmatrix} 1 & 2 & 0 & 1 \\ 0 & 0 & 1 & -1 \\ 0 & 0 & 0 & 0 \end{pmatrix}, \quad
\begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix}, \quad
\begin{pmatrix} 0 & 1 & 3 & 0 \\ 0 & 0 & 0 & 1 \end{pmatrix} 1 0 0 2 0 0 0 1 0 1 − 1 0 , 1 0 0 0 1 0 0 0 1 , ( 0 0 1 0 3 0 0 1 )
一方、次の行列は簡約階段行列ではない(括弧内は破れている条件)。
( 1 2 3 0 0 1 ) ( 条件 4 ) , ( 0 1 1 0 ) ( 条件 3 ) , ( 2 0 0 1 ) ( 条件 2 ) , ( 0 0 1 0 ) ( 条件 1 ) \begin{pmatrix} 1 & 2 & 3 \\ 0 & 0 & 1 \end{pmatrix}\ (\text{条件 4}), \quad
\begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}\ (\text{条件 3}), \quad
\begin{pmatrix} 2 & 0 \\ 0 & 1 \end{pmatrix}\ (\text{条件 2}), \quad
\begin{pmatrix} 0 & 0 \\ 1 & 0 \end{pmatrix}\ (\text{条件 1}) ( 1 0 2 0 3 1 ) ( 条件 4 ) , ( 0 1 1 0 ) ( 条件 3 ) , ( 2 0 0 1 ) ( 条件 2 ) , ( 0 1 0 0 ) ( 条件 1 )
定理 1.23 (簡約化の存在)任意の行列は、有限回の行基本変形によって簡約階段行列に変形できる。
証明. A ∈ M m , n ( K ) A \in \operatorname{M}_{m,n}(K) A ∈ M m , n ( K ) として、行の数 m m m に関する帰納法で示す。A = O A = O A = O ならすでに簡約階段行列である。A ≠ O A \neq O A = O とし、0 でない成分を含む最も左の列を第 j 1 j_1 j 1 列とする。
第 j 1 j_1 j 1 列の 0 でない成分がある行を (I) で第 1 行に移し、(II) でその成分を 1 にする。次に (III) で他の行から第 1 行の定数倍を引いて、第 j 1 j_1 j 1 列の第 1 行以外の成分をすべて 0 にする。得られた行列を A 1 A_1 A 1 とすると、A 1 A_1 A 1 の第 1 列〜第 j 1 − 1 j_1 - 1 j 1 − 1 列は 0、第 j 1 j_1 j 1 列は e 1 e_1 e 1 である。
m = 1 m = 1 m = 1 ならこれで終わりである。m ≥ 2 m \geq 2 m ≥ 2 のとき、A 1 A_1 A 1 の第 2 行〜第 m m m 行からなる ( m − 1 ) × n (m-1) \times n ( m − 1 ) × n 行列 B B B に帰納法の仮定を適用し、第 2 行〜第 m m m 行だけに対する行基本変形で B B B を簡約階段行列 B ′ B' B ′ にする。これらの変形は第 1 行を変えない。B B B の第 1 列〜第 j 1 j_1 j 1 列は 0 なので、B ′ B' B ′ でもそうであり、B ′ B' B ′ の主成分はすべて第 j 1 j_1 j 1 列より右にある。最後に、B ′ B' B ′ の各主成分(行列全体では第 i i i 行・第 j i j_i j i 列、i ≥ 2 i \geq 2 i ≥ 2 )について、第 1 行から第 i i i 行の a 1 j i a_{1 j_i} a 1 j i 倍(a 1 j i a_{1j_i} a 1 j i はその時点での第 1 行の第 j i j_i j i 列成分)を引く。第 i i i 行は第 j i j_i j i 列より左の成分と、他の主列の成分が 0 なので、この操作は第 1 行の第 j i j_i j i 列の成分だけを 0 にし、第 1 行の主成分や他の主列の成分を変えない。こうして得られた行列は簡約階段行列である。□ \square □
この証明の手順は掃き出し法 (ガウス–ジョルダンの消去法, Gauss–Jordan elimination)と呼ばれる。次に、簡約階段行列の形の連立方程式がどう解けるかを確認する。
補題 1.24 R ∈ M m , n ( K ) R \in \operatorname{M}_{m,n}(K) R ∈ M m , n ( K ) を簡約階段行列とし、零でない行の数を r r r 、主列を j 1 < ⋯ < j r j_1 < \cdots < j_r j 1 < ⋯ < j r 、主列でない列の番号の集合を F F F とする。c = t ( c 1 , … , c m ) ∈ K m c = {}^t(c_1, \dots, c_m) \in K^m c = t ( c 1 , … , c m ) ∈ K m が c r + 1 = ⋯ = c m = 0 c_{r+1} = \cdots = c_m = 0 c r + 1 = ⋯ = c m = 0 を満たすとする。このとき、任意に与えた t j ∈ K t_j \in K t j ∈ K ( j ∈ F ) (j \in F) ( j ∈ F ) に対し、R x = c Rx = c R x = c の解 x x x で x j = t j x_j = t_j x j = t j ( j ∈ F ) (j \in F) ( j ∈ F ) を満たすものがただ一つ存在し、
x j i = c i − ∑ j ∈ F r i j t j ( i = 1 , … , r ) x_{j_i} = c_i - \sum_{j \in F} r_{ij} t_j \qquad (i = 1, \dots, r) x j i = c i − j ∈ F ∑ r ij t j ( i = 1 , … , r )
で与えられる。
証明. 第 i i i 行 ( i ≤ r ) (i \leq r) ( i ≤ r ) の成分は、主成分 r i j i = 1 r_{i j_i} = 1 r i j i = 1 を除けば、他の主列では 0(条件 4)である。よって R x = c Rx = c R x = c の第 i i i 式は x j i + ∑ j ∈ F r i j x j = c i x_{j_i} + \sum_{j \in F} r_{ij} x_j = c_i x j i + ∑ j ∈ F r ij x j = c i となる。第 r + 1 r+1 r + 1 式以降は 0 = 0 0 = 0 0 = 0 である。したがって x j x_j x j ( j ∈ F ) (j \in F) ( j ∈ F ) を自由に決めると、x j i x_{j_i} x j i は上の式でただ一通りに決まり、そのとき全式が成り立つ。□ \square □
x j x_j x j ( j ∈ F ) (j \in F) ( j ∈ F ) を自由変数 という。補題の特別な場合として、次が重要である:主列でない列 j j j があれば、R x = 0 Rx = 0 R x = 0 の解で x j = 1 x_j = 1 x j = 1 となるものが存在する (t j = 1 t_j = 1 t j = 1 、他の自由変数を 0 とすればよい)。
定理 1.25 (簡約階段行列の一意性)行列 A A A を行基本変形で変形して得られる簡約階段行列は、変形の手順によらずただ一つに定まる。
証明. R , S R, S R , S がともに A A A と行同値な簡約階段行列であるとする。行同値は同値関係なので R R R と S S S は行同値であり、命題 1.20(b = 0 b = 0 b = 0 の場合)より
R x = 0 ⟺ S x = 0 ( x ∈ K n ) (1) Rx = 0 \iff Sx = 0 \qquad (x \in K^n) \tag{1} R x = 0 ⟺ S x = 0 ( x ∈ K n ) ( 1 )
が成り立つ。列の数 n n n に関する帰納法で R = S R = S R = S を示す。
n = 1 n = 1 n = 1 のとき、m × 1 m \times 1 m × 1 の簡約階段行列は O O O か e 1 e_1 e 1 のどちらかである。R x = 0 Rx = 0 R x = 0 が解 x = 1 x = 1 x = 1 をもつのは R = O R = O R = O のときに限るから、(1) より R = O ⟺ S = O R = O \iff S = O R = O ⟺ S = O であり、R = S R = S R = S となる。
n ≥ 2 n \geq 2 n ≥ 2 とし、n − 1 n - 1 n − 1 列の場合は正しいとする。R , S R, S R , S から最後の列を除いた行列を R ′ , S ′ R', S' R ′ , S ′ とする。簡約階段行列から最後の列を除いたものはまた簡約階段行列である(主成分が最後の列にある行は零行になるが、その行は最後の主成分をもつ行なので零でない行の最も下にあり、条件 1 は保たれる。他の条件は明らかに保たれる)。また S = P R S = PR S = P R なら S ′ = P R ′ S' = PR' S ′ = P R ′ なので R ′ R' R ′ と S ′ S' S ′ は行同値である。よって帰納法の仮定から R ′ = S ′ R' = S' R ′ = S ′ である。
R ≠ S R \neq S R = S と仮定して矛盾を導く。このとき R R R と S S S の最後の列 r n , s n r_n, s_n r n , s n が異なる。R x = 0 Rx = 0 R x = 0 を満たす任意の x x x について (1) より S x = 0 Sx = 0 S x = 0 でもあるから ( R − S ) x = 0 (R - S)x = 0 ( R − S ) x = 0 である。R − S R - S R − S は最後の列 r n − s n r_n - s_n r n − s n 以外が 0 なので ( R − S ) x = x n ( r n − s n ) (R - S)x = x_n (r_n - s_n) ( R − S ) x = x n ( r n − s n ) であり、r n − s n ≠ 0 r_n - s_n \neq 0 r n − s n = 0 から x n = 0 x_n = 0 x n = 0 が従う。つまり R x = 0 Rx = 0 R x = 0 のすべての解は x n = 0 x_n = 0 x n = 0 を満たす。補題 1.24 の後の注意により、第 n n n 列は R R R の主列でなければならない。同じ議論で第 n n n 列は S S S の主列でもある。
R ′ R' R ′ の零でない行の数を k k k とする。R R R の第 n n n 列の主成分がある行は、第 n n n 列以外の成分がすべて 0 の行、すなわち R ′ R' R ′ で零行になる行であり、条件 1, 3 からそれは第 k + 1 k+1 k + 1 行である。条件 2, 4 より r n = e k + 1 r_n = e_{k+1} r n = e k + 1 となる。S ′ = R ′ S' = R' S ′ = R ′ なので同様に s n = e k + 1 s_n = e_{k+1} s n = e k + 1 であり、r n = s n r_n = s_n r n = s n となって矛盾する。□ \square □
注意
一意なのは簡約 階段行列であって、階段行列は一意ではない。例えば第 1 行が ( 1 , 2 ) (1, 2) ( 1 , 2 ) 、第 2 行が ( 0 , 1 ) (0, 1) ( 0 , 1 ) の行列は階段行列だが、さらに変形して I 2 I_2 I 2 にもなる。どちらも同じ行列と行同値な階段行列である。
例 1.26 例 1.2 の拡大係数行列を簡約化する。矢印の上に行った変形を書く(R i R_i R i は第 i i i 行を表す)。
( 1 2 1 0 4 2 4 1 1 7 3 6 − 1 4 8 ) → R 3 − 3 R 1 R 2 − 2 R 1 ( 1 2 1 0 4 0 0 − 1 1 − 1 0 0 − 4 4 − 4 ) → ( − 1 ) × R 2 R 3 − 4 R 2 ( 1 2 1 0 4 0 0 1 − 1 1 0 0 0 0 0 ) → R 1 − R 2 ( 1 2 0 1 3 0 0 1 − 1 1 0 0 0 0 0 ) \left(\begin{array}{cccc|c} 1 & 2 & 1 & 0 & 4 \\ 2 & 4 & 1 & 1 & 7 \\ 3 & 6 & -1 & 4 & 8 \end{array}\right)
\xrightarrow[R_3 - 3R_1]{R_2 - 2R_1}
\left(\begin{array}{cccc|c} 1 & 2 & 1 & 0 & 4 \\ 0 & 0 & -1 & 1 & -1 \\ 0 & 0 & -4 & 4 & -4 \end{array}\right)
\xrightarrow[(-1) \times R_2]{R_3 - 4R_2}
\left(\begin{array}{cccc|c} 1 & 2 & 1 & 0 & 4 \\ 0 & 0 & 1 & -1 & 1 \\ 0 & 0 & 0 & 0 & 0 \end{array}\right)
\xrightarrow{R_1 - R_2}
\left(\begin{array}{cccc|c} 1 & 2 & 0 & 1 & 3 \\ 0 & 0 & 1 & -1 & 1 \\ 0 & 0 & 0 & 0 & 0 \end{array}\right) 1 2 3 2 4 6 1 1 − 1 0 1 4 4 7 8 R 2 − 2 R 1 R 3 − 3 R 1 1 0 0 2 0 0 1 − 1 − 4 0 1 4 4 − 1 − 4 R 3 − 4 R 2 ( − 1 ) × R 2 1 0 0 2 0 0 1 1 0 0 − 1 0 4 1 0 R 1 − R 2 1 0 0 2 0 0 0 1 0 1 − 1 0 3 1 0
最後の行列が簡約階段行列である。主列は第 1, 3 列、自由変数は x 2 , x 4 x_2, x_4 x 2 , x 4 である。これについては 1.6 節で解を書き下す。
1.5 階数
定理 1.25 によって、次の定義は変形の手順によらない。
定義 1.27 (階数, rank)行列 A A A と行同値な簡約階段行列の零でない行の数(=主成分の数)を、A A A の階数 といい rank A \operatorname{rank} A rank A と書く。
命題 1.28 A ∈ M m , n ( K ) A \in \operatorname{M}_{m,n}(K) A ∈ M m , n ( K ) とする。
rank A ≤ min ( m , n ) \operatorname{rank} A \leq \min(m, n) rank A ≤ min ( m , n ) 。
A A A と行同値な行列 B B B について rank B = rank A \operatorname{rank} B = \operatorname{rank} A rank B = rank A 。
証明. 1:主成分は各行に高々 1 個、各列に高々 1 個である。2:行同値は同値関係なので、A A A と B B B は同じ簡約階段行列と行同値であり、定理 1.25 から主張が従う。□ \square □
例 1.26 の係数行列 A A A の階数は 2、拡大係数行列 ( A ∣ b ) (A \mid b) ( A ∣ b ) の階数も 2 である。
階数は「行列が本当に含んでいる独立な情報の量」を測る量である。例 1.26 では、第 3 式が ( − 5 ) × (-5) \times ( − 5 ) × (第 1 式)+ 4 × + 4 \times + 4 × (第 2 式)として他の式から導かれる(係数と右辺を比べて確かめよ)ため、実質的な方程式は 2 本しかない。この直観は第3章 で「階数 = = = 像の次元」として正確になり、さらに行の階数と列の階数が等しい という非自明な定理(第3章)につながる。
1.6 連立一次方程式の解の存在と構造
定理 1.29 (解の存在条件)A ∈ M m , n ( K ) A \in \operatorname{M}_{m,n}(K) A ∈ M m , n ( K ) 、b ∈ K m b \in K^m b ∈ K m とする。連立方程式 A x = b Ax = b A x = b が解をもつための必要十分条件は
rank A = rank ( A ∣ b ) \operatorname{rank} A = \operatorname{rank}(A \mid b) rank A = rank ( A ∣ b )
である。
証明. ( A ∣ b ) (A \mid b) ( A ∣ b ) を簡約階段行列 ( R ∣ c ) (R \mid c) ( R ∣ c ) に変形する。定理 1.25 の証明中で注意したように R R R は簡約階段行列であり、R R R は A A A と行同値なので r : = rank A r := \operatorname{rank} A r := rank A は R R R の零でない行の数である。命題 1.20 より、A x = b Ax = b A x = b と R x = c Rx = c R x = c の解は一致する。
( R ∣ c ) (R \mid c) ( R ∣ c ) の零でない行の数は r r r か r + 1 r + 1 r + 1 である。r + 1 r + 1 r + 1 の場合、第 r + 1 r+1 r + 1 行は ( 0 ⋯ 0 ∣ 1 ) (0 \ \cdots \ 0 \mid 1) ( 0 ⋯ 0 ∣ 1 ) であり、方程式 0 = 1 0 = 1 0 = 1 を表すので解はない。r r r の場合、c r + 1 = ⋯ = c m = 0 c_{r+1} = \cdots = c_m = 0 c r + 1 = ⋯ = c m = 0 であり、補題 1.24 により解が存在する。rank ( A ∣ b ) \operatorname{rank}(A \mid b) rank ( A ∣ b ) は ( R ∣ c ) (R \mid c) ( R ∣ c ) の零でない行の数だから、主張が従う。□ \square □
定理 1.30 (解の構造)A x = b Ax = b A x = b の解 x 0 x_0 x 0 が一つ与えられたとする(これを特殊解 , particular solution という)。このとき A x = b Ax = b A x = b の解全体は
{ x 0 + y ∣ y ∈ K n , A y = 0 } \lbrace x_0 + y \mid y \in K^n,\ Ay = 0 \rbrace { x 0 + y ∣ y ∈ K n , A y = 0 }
である。すなわち「一般解 = = = 特殊解 + + + 同次方程式 A y = 0 Ay = 0 A y = 0 の一般解」である。
証明. A y = 0 Ay = 0 A y = 0 なら A ( x 0 + y ) = A x 0 + A y = b A(x_0 + y) = Ax_0 + Ay = b A ( x 0 + y ) = A x 0 + A y = b 。逆に A x = b Ax = b A x = b なら y : = x − x 0 y := x - x_0 y := x − x 0 は A y = b − b = 0 Ay = b - b = 0 A y = b − b = 0 を満たし、x = x 0 + y x = x_0 + y x = x 0 + y である。□ \square □
b = 0 b = 0 b = 0 の方程式 A x = 0 Ax = 0 A x = 0 を同次 (斉次)方程式 (homogeneous equation) という。同次方程式は常に自明な解 x = 0 x = 0 x = 0 をもつ。
系 1.31 A ∈ M m , n ( K ) A \in \operatorname{M}_{m,n}(K) A ∈ M m , n ( K ) 、r = rank A r = \operatorname{rank} A r = rank A とする。
A x = 0 Ax = 0 A x = 0 の解は、n − r n - r n − r 個の自由変数を任意に選ぶことで一意に定まる。特に A x = 0 Ax = 0 A x = 0 が自明でない解をもつための必要十分条件は r < n r < n r < n である。
m < n m < n m < n (方程式の数 < < < 未知数の数)ならば、A x = 0 Ax = 0 A x = 0 は自明でない解をもつ。
証明. 1 は補題 1.24(c = 0 c = 0 c = 0 )で、自由変数が n − r n - r n − r 個あることから従う。r < n r < n r < n なら自由変数の一つを 1 とする解は 0 でない。r = n r = n r = n なら自由変数がなく、解は x = 0 x = 0 x = 0 だけである。2 は r ≤ m < n r \leq m < n r ≤ m < n と 1 から従う。□ \square □
系 1.31 の 2 は単純だが、第2章 で次元の理論を展開するときの決定的な道具になる。
例 1.32 例 1.26 の結果から、例 1.2 の連立方程式は
x 1 + 2 x 2 + x 4 = 3 , x 3 − x 4 = 1 x_1 + 2x_2 + x_4 = 3, \qquad x_3 - x_4 = 1 x 1 + 2 x 2 + x 4 = 3 , x 3 − x 4 = 1
と同値である。自由変数を x 2 = s x_2 = s x 2 = s 、x 4 = t x_4 = t x 4 = t とおくと x 1 = 3 − 2 s − t x_1 = 3 - 2s - t x 1 = 3 − 2 s − t 、x 3 = 1 + t x_3 = 1 + t x 3 = 1 + t であり、解の全体は
x = ( 3 0 1 0 ) + s ( − 2 1 0 0 ) + t ( − 1 0 1 1 ) ( s , t ∈ K ) x = \begin{pmatrix} 3 \\ 0 \\ 1 \\ 0 \end{pmatrix} + s \begin{pmatrix} -2 \\ 1 \\ 0 \\ 0 \end{pmatrix} + t \begin{pmatrix} -1 \\ 0 \\ 1 \\ 1 \end{pmatrix} \qquad (s, t \in K) x = 3 0 1 0 + s − 2 1 0 0 + t − 1 0 1 1 ( s , t ∈ K )
である。第 1 項が特殊解、残りが同次方程式の一般解である。検算 :特殊解を元の式に代入すると 3 + 0 + 1 = 4 3 + 0 + 1 = 4 3 + 0 + 1 = 4 、6 + 0 + 1 + 0 = 7 6 + 0 + 1 + 0 = 7 6 + 0 + 1 + 0 = 7 、9 + 0 − 1 + 0 = 8 9 + 0 - 1 + 0 = 8 9 + 0 − 1 + 0 = 8 。t ( − 2 , 1 , 0 , 0 ) {}^t(-2, 1, 0, 0) t ( − 2 , 1 , 0 , 0 ) を係数行列に掛けると ( − 2 + 2 , − 4 + 4 , − 6 + 6 ) = ( 0 , 0 , 0 ) (-2 + 2,\ -4 + 4,\ -6 + 6) = (0, 0, 0) ( − 2 + 2 , − 4 + 4 , − 6 + 6 ) = ( 0 , 0 , 0 ) 、t ( − 1 , 0 , 1 , 1 ) {}^t(-1, 0, 1, 1) t ( − 1 , 0 , 1 , 1 ) では ( − 1 + 1 , − 2 + 1 + 1 , − 3 − 1 + 4 ) = ( 0 , 0 , 0 ) (-1 + 1,\ -2 + 1 + 1,\ -3 - 1 + 4) = (0, 0, 0) ( − 1 + 1 , − 2 + 1 + 1 , − 3 − 1 + 4 ) = ( 0 , 0 , 0 ) となり、正しい。
一方、第 3 式の右辺を 8 から 9 に変えると、同じ変形で最後の行が ( 0 0 0 0 ∣ 1 ) (0 \ 0 \ 0 \ 0 \mid 1) ( 0 0 0 0 ∣ 1 ) になる(9 − 3 ⋅ 4 = − 3 9 - 3 \cdot 4 = -3 9 − 3 ⋅ 4 = − 3 、− 3 − 4 ⋅ ( − 1 ) = 1 -3 - 4 \cdot (-1) = 1 − 3 − 4 ⋅ ( − 1 ) = 1 )。このとき rank A = 2 < 3 = rank ( A ∣ b ) \operatorname{rank} A = 2 < 3 = \operatorname{rank}(A \mid b) rank A = 2 < 3 = rank ( A ∣ b ) であり、解は存在しない。
1.7 逆行列
正方行列 A A A が正則であるための条件を、これまでの道具で完全に特徴づける。
定理 1.34 (正則性の同値条件)A ∈ M n ( K ) A \in \operatorname{M}_n(K) A ∈ M n ( K ) について、次は同値である。
A A A は正則である。
rank A = n \operatorname{rank} A = n rank A = n 。
A A A の簡約階段行列は I n I_n I n である。
A A A は基本行列の積として書ける。
A x = 0 Ax = 0 A x = 0 の解は x = 0 x = 0 x = 0 のみである。
B A = I n BA = I_n B A = I n を満たす B ∈ M n ( K ) B \in \operatorname{M}_n(K) B ∈ M n ( K ) が存在する(左逆行列の存在)。
A C = I n AC = I_n A C = I n を満たす C ∈ M n ( K ) C \in \operatorname{M}_n(K) C ∈ M n ( K ) が存在する(右逆行列の存在)。
さらに 6, 7 の B , C B, C B , C は A − 1 A^{-1} A − 1 に等しい。
証明. (2 ⇒ 3):n × n n \times n n × n の簡約階段行列で主成分が n n n 個あれば、条件 3 より j i = i j_i = i j i = i ( i = 1 , … , n ) (i = 1, \dots, n) ( i = 1 , … , n ) であり、条件 2, 4 からそれは I n I_n I n である。
(3 ⇒ 4):基本行列 E 1 , … , E k E_1, \dots, E_k E 1 , … , E k で E k ⋯ E 1 A = I n E_k \cdots E_1 A = I_n E k ⋯ E 1 A = I n となる。命題 1.18 より A = E 1 − 1 ⋯ E k − 1 A = E_1^{-1} \cdots E_k^{-1} A = E 1 − 1 ⋯ E k − 1 であり、各 E i − 1 E_i^{-1} E i − 1 も基本行列である。
(4 ⇒ 1):基本行列は正則であり(命題 1.18)、正則行列の積は正則である(命題 1.14)。
(1 ⇒ 6), (1 ⇒ 7):B = C = A − 1 B = C = A^{-1} B = C = A − 1 とすればよい。
(6 ⇒ 5):A x = 0 Ax = 0 A x = 0 なら x = I n x = B A x = B 0 = 0 x = I_n x = BAx = B0 = 0 x = I n x = B A x = B 0 = 0 。
(5 ⇒ 2):系 1.31 の 1 より、自明な解しかもたないなら rank A = n \operatorname{rank} A = n rank A = n 。
ここまでで 1〜6 の同値性が示された。
(7 ⇒ 1):A C = I n AC = I_n A C = I n は、C C C が左逆行列 A A A をもつことを意味する。すでに示した (6 ⇒ 1) を C C C に適用すると C C C は正則であり、A = A C C − 1 = C − 1 A = A C C^{-1} = C^{-1} A = A C C − 1 = C − 1 も正則である。
最後に、B A = I n BA = I_n B A = I n かつ A A A が正則なら B = B A A − 1 = A − 1 B = B A A^{-1} = A^{-1} B = B A A − 1 = A − 1 。C C C についても同様である。□ \square □
定理 1.34 の 6, 7 は、正方行列については片側だけ確かめれば逆行列である ことを保証する。これは有限次元性に依存する事実であり、無限次元では成り立たない(第3章 の例 3.8、関数解析 第1章 のシフト作用素の例を参照)。
逆行列の計算
定理 1.34 の証明は、そのまま計算法を与える。A A A に行基本変形を施して I n I_n I n に到達したとすると、E k ⋯ E 1 A = I n E_k \cdots E_1 A = I_n E k ⋯ E 1 A = I n より E k ⋯ E 1 = A − 1 E_k \cdots E_1 = A^{-1} E k ⋯ E 1 = A − 1 である。同じ変形を I n I_n I n に施すと E k ⋯ E 1 I n = A − 1 E_k \cdots E_1 I_n = A^{-1} E k ⋯ E 1 I n = A − 1 が得られる。そこで n × 2 n n \times 2n n × 2 n 行列 ( A ∣ I n ) (A \mid I_n) ( A ∣ I n ) を簡約化すれば
( A ∣ I n ) ⟶ ( I n ∣ A − 1 ) (A \mid I_n) \longrightarrow (I_n \mid A^{-1}) ( A ∣ I n ) ⟶ ( I n ∣ A − 1 )
となる。途中で左半分の簡約階段行列が I n I_n I n にならない(零行が現れる)ことがわかれば、A A A は正則でない。
例 1.35 次の行列の逆行列を求める。
A = ( 1 2 1 2 5 3 1 3 3 ) A = \begin{pmatrix} 1 & 2 & 1 \\ 2 & 5 & 3 \\ 1 & 3 & 3 \end{pmatrix} A = 1 2 1 2 5 3 1 3 3
( 1 2 1 1 0 0 2 5 3 0 1 0 1 3 3 0 0 1 ) → R 3 − R 1 R 2 − 2 R 1 ( 1 2 1 1 0 0 0 1 1 − 2 1 0 0 1 2 − 1 0 1 ) → R 3 − R 2 ( 1 2 1 1 0 0 0 1 1 − 2 1 0 0 0 1 1 − 1 1 ) \left(\begin{array}{ccc|ccc} 1 & 2 & 1 & 1 & 0 & 0 \\ 2 & 5 & 3 & 0 & 1 & 0 \\ 1 & 3 & 3 & 0 & 0 & 1 \end{array}\right)
\xrightarrow[R_3 - R_1]{R_2 - 2R_1}
\left(\begin{array}{ccc|ccc} 1 & 2 & 1 & 1 & 0 & 0 \\ 0 & 1 & 1 & -2 & 1 & 0 \\ 0 & 1 & 2 & -1 & 0 & 1 \end{array}\right)
\xrightarrow{R_3 - R_2}
\left(\begin{array}{ccc|ccc} 1 & 2 & 1 & 1 & 0 & 0 \\ 0 & 1 & 1 & -2 & 1 & 0 \\ 0 & 0 & 1 & 1 & -1 & 1 \end{array}\right) 1 2 1 2 5 3 1 3 3 1 0 0 0 1 0 0 0 1 R 2 − 2 R 1 R 3 − R 1 1 0 0 2 1 1 1 1 2 1 − 2 − 1 0 1 0 0 0 1 R 3 − R 2 1 0 0 2 1 0 1 1 1 1 − 2 1 0 1 − 1 0 0 1
続けて第 3 列、第 2 列の順に上を掃き出す。
→ R 1 − R 3 R 2 − R 3 ( 1 2 0 0 1 − 1 0 1 0 − 3 2 − 1 0 0 1 1 − 1 1 ) → R 1 − 2 R 2 ( 1 0 0 6 − 3 1 0 1 0 − 3 2 − 1 0 0 1 1 − 1 1 ) \xrightarrow[R_1 - R_3]{R_2 - R_3}
\left(\begin{array}{ccc|ccc} 1 & 2 & 0 & 0 & 1 & -1 \\ 0 & 1 & 0 & -3 & 2 & -1 \\ 0 & 0 & 1 & 1 & -1 & 1 \end{array}\right)
\xrightarrow{R_1 - 2R_2}
\left(\begin{array}{ccc|ccc} 1 & 0 & 0 & 6 & -3 & 1 \\ 0 & 1 & 0 & -3 & 2 & -1 \\ 0 & 0 & 1 & 1 & -1 & 1 \end{array}\right) R 2 − R 3 R 1 − R 3 1 0 0 2 1 0 0 0 1 0 − 3 1 1 2 − 1 − 1 − 1 1 R 1 − 2 R 2 1 0 0 0 1 0 0 0 1 6 − 3 1 − 3 2 − 1 1 − 1 1
よって
A − 1 = ( 6 − 3 1 − 3 2 − 1 1 − 1 1 ) A^{-1} = \begin{pmatrix} 6 & -3 & 1 \\ -3 & 2 & -1 \\ 1 & -1 & 1 \end{pmatrix} A − 1 = 6 − 3 1 − 3 2 − 1 1 − 1 1
である。検算 :A A − 1 A A^{-1} A A − 1 の第 1 行は ( 6 − 6 + 1 , − 3 + 4 − 1 , 1 − 2 + 1 ) = ( 1 , 0 , 0 ) (6 - 6 + 1,\ -3 + 4 - 1,\ 1 - 2 + 1) = (1, 0, 0) ( 6 − 6 + 1 , − 3 + 4 − 1 , 1 − 2 + 1 ) = ( 1 , 0 , 0 ) 、第 2 行は ( 12 − 15 + 3 , − 6 + 10 − 3 , 2 − 5 + 3 ) = ( 0 , 1 , 0 ) (12 - 15 + 3,\ -6 + 10 - 3,\ 2 - 5 + 3) = (0, 1, 0) ( 12 − 15 + 3 , − 6 + 10 − 3 , 2 − 5 + 3 ) = ( 0 , 1 , 0 ) 、第 3 行は ( 6 − 9 + 3 , − 3 + 6 − 3 , 1 − 3 + 3 ) = ( 0 , 0 , 1 ) (6 - 9 + 3,\ -3 + 6 - 3,\ 1 - 3 + 3) = (0, 0, 1) ( 6 − 9 + 3 , − 3 + 6 − 3 , 1 − 3 + 3 ) = ( 0 , 0 , 1 ) となり正しい。定理 1.34 により、片側の積を確かめれば十分である。A A A が対称行列なので A − 1 A^{-1} A − 1 も対称になっている(命題 1.14 の 3)ことにも注意する。
例 1.36 (2 次の逆行列)2 次正方行列 A A A と、それから作る行列 A ~ \tilde{A} A ~ を
A = ( a b c d ) , A ~ = ( d − b − c a ) A = \begin{pmatrix} a & b \\ c & d \end{pmatrix}, \qquad \tilde{A} = \begin{pmatrix} d & -b \\ -c & a \end{pmatrix} A = ( a c b d ) , A ~ = ( d − c − b a )
とし、Δ = a d − b c \Delta = ad - bc Δ = a d − b c とおく。直接計算で A A ~ = Δ I 2 A\tilde{A} = \Delta I_2 A A ~ = Δ I 2 となる。よって Δ ≠ 0 \Delta \neq 0 Δ = 0 なら A A A は正則で A − 1 = Δ − 1 A ~ A^{-1} = \Delta^{-1} \tilde{A} A − 1 = Δ − 1 A ~ である(定理 1.34 の 7)。Δ = 0 \Delta = 0 Δ = 0 なら A A A は正則でない:もし正則なら A A ~ = O A \tilde{A} = O A A ~ = O に左から A − 1 A^{-1} A − 1 を掛けて A ~ = O \tilde{A} = O A ~ = O となり、a = b = c = d = 0 a = b = c = d = 0 a = b = c = d = 0 すなわち A = O A = O A = O となるが、O O O は正則でないので矛盾する。量 Δ \Delta Δ は 2 2 2 次の行列式 であり、一般の n n n 次への拡張が第4章 の主題である。
ヒント
逆行列を求めたら、必ず元の行列と掛けて単位行列になるか確かめる習慣をつけよう。手計算での符号の誤りは非常によく起こるが、検算はすぐにできる。
1.8 ブロック行列
大きな行列を縦横の線で区切り、小さな行列(ブロック , block)を成分とする行列とみなすと、計算や議論の見通しがよくなることが多い。
命題 1.37 (ブロック積)A ∈ M m , n ( K ) A \in \operatorname{M}_{m,n}(K) A ∈ M m , n ( K ) 、B ∈ M n , p ( K ) B \in \operatorname{M}_{n,p}(K) B ∈ M n , p ( K ) とし、m = m 1 + ⋯ + m s m = m_1 + \cdots + m_s m = m 1 + ⋯ + m s 、n = n 1 + ⋯ + n t n = n_1 + \cdots + n_t n = n 1 + ⋯ + n t 、p = p 1 + ⋯ + p u p = p_1 + \cdots + p_u p = p 1 + ⋯ + p u と分割して
A = ( A 11 ⋯ A 1 t ⋮ ⋮ A s 1 ⋯ A s t ) , B = ( B 11 ⋯ B 1 u ⋮ ⋮ B t 1 ⋯ B t u ) A = \begin{pmatrix} A_{11} & \cdots & A_{1t} \\ \vdots & & \vdots \\ A_{s1} & \cdots & A_{st} \end{pmatrix}, \quad
B = \begin{pmatrix} B_{11} & \cdots & B_{1u} \\ \vdots & & \vdots \\ B_{t1} & \cdots & B_{tu} \end{pmatrix} A = A 11 ⋮ A s 1 ⋯ ⋯ A 1 t ⋮ A s t , B = B 11 ⋮ B t 1 ⋯ ⋯ B 1 u ⋮ B t u
と書く(A I J A_{IJ} A I J は m I × n J m_I \times n_J m I × n J 行列、B J K B_{JK} B J K は n J × p K n_J \times p_K n J × p K 行列)。このとき A B AB A B を同様に分割したブロックは
( A B ) I K = ∑ J = 1 t A I J B J K (AB)_{IK} = \sum_{J=1}^{t} A_{IJ} B_{JK} ( A B ) I K = J = 1 ∑ t A I J B J K
で与えられる。つまり、ブロックを成分とみなして通常の積の規則で計算してよい。
証明. A B AB A B の ( i , k ) (i, k) ( i , k ) 成分 ∑ j = 1 n a i j b j k \sum_{j=1}^{n} a_{ij} b_{jk} ∑ j = 1 n a ij b j k で、i i i が第 I I I ブロック、k k k が第 K K K ブロックに属するとする。和を j j j の属するブロック J = 1 , … , t J = 1, \dots, t J = 1 , … , t ごとに分けると、第 J J J ブロックについての部分和は A I J B J K A_{IJ} B_{JK} A I J B J K の対応する成分である。□ \square □
ただし、ブロックどうしの積は非可換なので、A I J B J K A_{IJ} B_{JK} A I J B J K の順序を入れ替えてはいけない。
例 1.38 (ブロック三角行列の逆行列)A ∈ GL k ( K ) A \in \operatorname{GL}_k(K) A ∈ GL k ( K ) 、D ∈ GL l ( K ) D \in \operatorname{GL}_l(K) D ∈ GL l ( K ) 、B ∈ M k , l ( K ) B \in \operatorname{M}_{k,l}(K) B ∈ M k , l ( K ) とする。ブロック積で確かめられるように
( A B O D ) ( A − 1 − A − 1 B D − 1 O D − 1 ) = ( I k − B D − 1 + B D − 1 O I l ) = I k + l \begin{pmatrix} A & B \\ O & D \end{pmatrix} \begin{pmatrix} A^{-1} & -A^{-1} B D^{-1} \\ O & D^{-1} \end{pmatrix}
= \begin{pmatrix} I_k & -BD^{-1} + BD^{-1} \\ O & I_l \end{pmatrix} = I_{k+l} ( A O B D ) ( A − 1 O − A − 1 B D − 1 D − 1 ) = ( I k O − B D − 1 + B D − 1 I l ) = I k + l
であるから、定理 1.34 より左の行列は正則で、右の行列がその逆行列である。特に B = O B = O B = O の場合、ブロック対角行列 について
( A O O D ) − 1 = ( A − 1 O O D − 1 ) \begin{pmatrix} A & O \\ O & D \end{pmatrix}^{-1} = \begin{pmatrix} A^{-1} & O \\ O & D^{-1} \end{pmatrix} ( A O O D ) − 1 = ( A − 1 O O D − 1 )
であり、ブロック対角行列の積・べきはブロックごとに計算できる。この事実は、固有値やジョルダン標準形の理論(第5章 、第6章 )で、空間を不変部分空間の直和に分けて考えるときの基本的な道具になる。
命題 1.39 (シューア補行列による分解)A ∈ GL k ( K ) A \in \operatorname{GL}_k(K) A ∈ GL k ( K ) 、B ∈ M k , l ( K ) B \in \operatorname{M}_{k,l}(K) B ∈ M k , l ( K ) 、C ∈ M l , k ( K ) C \in \operatorname{M}_{l,k}(K) C ∈ M l , k ( K ) 、D ∈ M l ( K ) D \in \operatorname{M}_l(K) D ∈ M l ( K ) とし、S = D − C A − 1 B S = D - CA^{-1}B S = D − C A − 1 B (シューア補行列 , Schur complement)とおく。このとき
( A B C D ) = ( I k O C A − 1 I l ) ( A O O S ) ( I k A − 1 B O I l ) \begin{pmatrix} A & B \\ C & D \end{pmatrix} = \begin{pmatrix} I_k & O \\ CA^{-1} & I_l \end{pmatrix} \begin{pmatrix} A & O \\ O & S \end{pmatrix} \begin{pmatrix} I_k & A^{-1}B \\ O & I_l \end{pmatrix} ( A C B D ) = ( I k C A − 1 O I l ) ( A O O S ) ( I k O A − 1 B I l )
が成り立つ。特に、左辺が正則であることと S S S が正則であることは同値である。
証明. 右辺の 3 つの行列を左から L , Δ , U L, \Delta, U L , Δ , U とおく。ブロック積により
Δ U = ( A B O S ) , L ( Δ U ) = ( A B C A − 1 A C A − 1 B + S ) = ( A B C D ) \Delta U = \begin{pmatrix} A & B \\ O & S \end{pmatrix}, \qquad L(\Delta U) = \begin{pmatrix} A & B \\ CA^{-1}A & CA^{-1}B + S \end{pmatrix} = \begin{pmatrix} A & B \\ C & D \end{pmatrix} Δ U = ( A O B S ) , L ( Δ U ) = ( A C A − 1 A B C A − 1 B + S ) = ( A C B D )
となる。後半:L , U L, U L , U は例 1.38 により正則(逆行列は非対角ブロックの符号を変えたもの)なので、左辺が正則であることは Δ \Delta Δ が正則であることと同値である。S S S が正則なら Δ \Delta Δ は例 1.38 より正則である。逆に Δ \Delta Δ が正則なら、その逆行列 X X X を同じサイズのブロック X 11 , X 12 , X 21 , X 22 X_{11}, X_{12}, X_{21}, X_{22} X 11 , X 12 , X 21 , X 22 に区切ると、Δ X = I \Delta X = I Δ X = I の右下ブロックを比較して S X 22 = I l S X_{22} = I_l S X 22 = I l となり(X 12 X_{12} X 12 に掛かるのは Δ \Delta Δ の左下ブロック O O O である)、定理 1.34 より S S S は正則である。□ \square □
この分解は、ガウスの消去法をブロック単位で行ったものにほかならない。第4章 ではこれを使ってブロック行列の行列式を計算し、第8章 では正定値性の判定に使う。
まとめ
行列の積 ( A B ) i k = ∑ j a i j b j k (AB)_{ik} = \sum_j a_{ij} b_{jk} ( A B ) ik = ∑ j a ij b j k は、写像 L A : x ↦ A x L_A\colon x \mapsto Ax L A : x ↦ A x の合成 L A ∘ L B = L A B L_A \circ L_B = L_{AB} L A ∘ L B = L A B に対応するように定められている。積は結合的だが非可換で、零因子や冪零行列がある。
A x Ax A x は A A A の列の一次結合 ∑ j x j a j \sum_j x_j a_j ∑ j x j a j である。
行基本変形は基本行列を左から掛けることに等しく、基本行列は正則である。行基本変形は連立方程式の解を変えない。
任意の行列は行基本変形で簡約階段行列に変形でき、その結果は一意である。階数は簡約階段行列の零でない行の数として定義される。
A x = b Ax = b A x = b が解をもつ ⟺ \iff ⟺ rank A = rank ( A ∣ b ) \operatorname{rank} A = \operatorname{rank}(A \mid b) rank A = rank ( A ∣ b ) 。解の全体は「特殊解+同次方程式の解」であり、自由変数は n − rank A n - \operatorname{rank} A n − rank A 個ある。
正方行列 A A A について、正則 ⟺ \iff ⟺ rank A = n \operatorname{rank} A = n rank A = n ⟺ \iff ⟺ 基本行列の積 ⟺ \iff ⟺ A x = 0 Ax = 0 A x = 0 が自明な解のみ ⟺ \iff ⟺ 片側逆行列が存在する。
逆行列は ( A ∣ I ) → ( I ∣ A − 1 ) (A \mid I) \to (I \mid A^{-1}) ( A ∣ I ) → ( I ∣ A − 1 ) の簡約化で計算できる。計算後は必ず検算する。
ブロック行列はブロックを成分とみなして掛けてよい(順序には注意)。シューア補行列による分解はブロック単位の消去法である。
演習問題
問題 1.1 ★ 次の A , B A, B A , B について A B AB A B と B A BA B A を計算し、それぞれ A A A の行・列にどのような操作をした結果になっているか説明せよ。
A = ( 1 2 3 4 ) , B = ( 0 1 1 0 ) A = \begin{pmatrix} 1 & 2 \\ 3 & 4 \end{pmatrix}, \qquad B = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} A = ( 1 3 2 4 ) , B = ( 0 1 1 0 )
解答
A B = ( 2 1 4 3 ) , B A = ( 3 4 1 2 ) AB = \begin{pmatrix} 2 & 1 \\ 4 & 3 \end{pmatrix}, \qquad BA = \begin{pmatrix} 3 & 4 \\ 1 & 2 \end{pmatrix} A B = ( 2 4 1 3 ) , B A = ( 3 1 4 2 )
B = P ( 1 , 2 ) B = P(1, 2) B = P ( 1 , 2 ) は基本行列であり、命題 1.17 のとおり、B A BA B A は A A A の第 1 行と第 2 行を入れ替えたもの、A B AB A B は第 1 列と第 2 列を入れ替えたものである。特に A B ≠ B A AB \neq BA A B = B A である。
問題 1.2 ★ 次の連立一次方程式を行基本変形で解け。
{ x + y + 2 z = 1 2 x + 3 y + 5 z = 3 x + 2 y + 4 z = 1 \left\lbrace
\begin{aligned}
x + y + 2z &= 1 \\
2x + 3y + 5z &= 3 \\
x + 2y + 4z &= 1
\end{aligned}
\right. ⎩ ⎨ ⎧ x + y + 2 z 2 x + 3 y + 5 z x + 2 y + 4 z = 1 = 3 = 1
解答
( 1 1 2 1 2 3 5 3 1 2 4 1 ) → R 3 − R 1 R 2 − 2 R 1 ( 1 1 2 1 0 1 1 1 0 1 2 0 ) → R 3 − R 2 ( 1 1 2 1 0 1 1 1 0 0 1 − 1 ) \left(\begin{array}{ccc|c} 1 & 1 & 2 & 1 \\ 2 & 3 & 5 & 3 \\ 1 & 2 & 4 & 1 \end{array}\right)
\xrightarrow[R_3 - R_1]{R_2 - 2R_1}
\left(\begin{array}{ccc|c} 1 & 1 & 2 & 1 \\ 0 & 1 & 1 & 1 \\ 0 & 1 & 2 & 0 \end{array}\right)
\xrightarrow{R_3 - R_2}
\left(\begin{array}{ccc|c} 1 & 1 & 2 & 1 \\ 0 & 1 & 1 & 1 \\ 0 & 0 & 1 & -1 \end{array}\right) 1 2 1 1 3 2 2 5 4 1 3 1 R 2 − 2 R 1 R 3 − R 1 1 0 0 1 1 1 2 1 2 1 1 0 R 3 − R 2 1 0 0 1 1 0 2 1 1 1 1 − 1
さらに R 2 − R 3 R_2 - R_3 R 2 − R 3 、R 1 − 2 R 3 R_1 - 2R_3 R 1 − 2 R 3 、R 1 − R 2 R_1 - R_2 R 1 − R 2 の順に変形すると
( 1 0 0 1 0 1 0 2 0 0 1 − 1 ) \left(\begin{array}{ccc|c} 1 & 0 & 0 & 1 \\ 0 & 1 & 0 & 2 \\ 0 & 0 & 1 & -1 \end{array}\right) 1 0 0 0 1 0 0 0 1 1 2 − 1
となる。よって ( x , y , z ) = ( 1 , 2 , − 1 ) (x, y, z) = (1, 2, -1) ( x , y , z ) = ( 1 , 2 , − 1 ) がただ一つの解である。検算:1 + 2 − 2 = 1 1 + 2 - 2 = 1 1 + 2 − 2 = 1 、2 + 6 − 5 = 3 2 + 6 - 5 = 3 2 + 6 − 5 = 3 、1 + 4 − 4 = 1 1 + 4 - 4 = 1 1 + 4 − 4 = 1 。
問題 1.3 ★ 問題 1.2 の係数行列 A A A の逆行列を ( A ∣ I 3 ) (A \mid I_3) ( A ∣ I 3 ) の簡約化で求め、A − 1 b A^{-1} b A − 1 b が問題 1.2 の解と一致することを確かめよ。
解答
( 1 1 2 1 0 0 2 3 5 0 1 0 1 2 4 0 0 1 ) → R 3 − R 1 R 2 − 2 R 1 ( 1 1 2 1 0 0 0 1 1 − 2 1 0 0 1 2 − 1 0 1 ) → R 3 − R 2 ( 1 1 2 1 0 0 0 1 1 − 2 1 0 0 0 1 1 − 1 1 ) \left(\begin{array}{ccc|ccc} 1 & 1 & 2 & 1 & 0 & 0 \\ 2 & 3 & 5 & 0 & 1 & 0 \\ 1 & 2 & 4 & 0 & 0 & 1 \end{array}\right)
\xrightarrow[R_3 - R_1]{R_2 - 2R_1}
\left(\begin{array}{ccc|ccc} 1 & 1 & 2 & 1 & 0 & 0 \\ 0 & 1 & 1 & -2 & 1 & 0 \\ 0 & 1 & 2 & -1 & 0 & 1 \end{array}\right)
\xrightarrow{R_3 - R_2}
\left(\begin{array}{ccc|ccc} 1 & 1 & 2 & 1 & 0 & 0 \\ 0 & 1 & 1 & -2 & 1 & 0 \\ 0 & 0 & 1 & 1 & -1 & 1 \end{array}\right) 1 2 1 1 3 2 2 5 4 1 0 0 0 1 0 0 0 1 R 2 − 2 R 1 R 3 − R 1 1 0 0 1 1 1 2 1 2 1 − 2 − 1 0 1 0 0 0 1 R 3 − R 2 1 0 0 1 1 0 2 1 1 1 − 2 1 0 1 − 1 0 0 1
R 2 − R 3 R_2 - R_3 R 2 − R 3 と R 1 − 2 R 3 R_1 - 2R_3 R 1 − 2 R 3 で第 3 列を掃き出すと第 1 行は ( 1 , 1 , 0 ∣ − 1 , 2 , − 2 ) (1, 1, 0 \mid -1, 2, -2) ( 1 , 1 , 0 ∣ − 1 , 2 , − 2 ) 、第 2 行は ( 0 , 1 , 0 ∣ − 3 , 2 , − 1 ) (0, 1, 0 \mid -3, 2, -1) ( 0 , 1 , 0 ∣ − 3 , 2 , − 1 ) となり、最後に R 1 − R 2 R_1 - R_2 R 1 − R 2 で
A − 1 = ( 2 0 − 1 − 3 2 − 1 1 − 1 1 ) A^{-1} = \begin{pmatrix} 2 & 0 & -1 \\ -3 & 2 & -1 \\ 1 & -1 & 1 \end{pmatrix} A − 1 = 2 − 3 1 0 2 − 1 − 1 − 1 1
を得る。検算:A A − 1 AA^{-1} A A − 1 の第 1 行は ( 2 − 3 + 2 , 0 + 2 − 2 , − 1 − 1 + 2 ) = ( 1 , 0 , 0 ) (2 - 3 + 2,\ 0 + 2 - 2,\ -1 - 1 + 2) = (1, 0, 0) ( 2 − 3 + 2 , 0 + 2 − 2 , − 1 − 1 + 2 ) = ( 1 , 0 , 0 ) 、第 2 行は ( 4 − 9 + 5 , 6 − 5 , − 2 − 3 + 5 ) = ( 0 , 1 , 0 ) (4 - 9 + 5,\ 6 - 5,\ -2 - 3 + 5) = (0, 1, 0) ( 4 − 9 + 5 , 6 − 5 , − 2 − 3 + 5 ) = ( 0 , 1 , 0 ) 、第 3 行は ( 2 − 6 + 4 , 4 − 4 , − 1 − 2 + 4 ) = ( 0 , 0 , 1 ) (2 - 6 + 4,\ 4 - 4,\ -1 - 2 + 4) = (0, 0, 1) ( 2 − 6 + 4 , 4 − 4 , − 1 − 2 + 4 ) = ( 0 , 0 , 1 ) 。また A − 1 t ( 1 , 3 , 1 ) = t ( 2 − 1 , − 3 + 6 − 1 , 1 − 3 + 1 ) = t ( 1 , 2 , − 1 ) A^{-1}\ {}^t(1, 3, 1) = {}^t(2 - 1,\ -3 + 6 - 1,\ 1 - 3 + 1) = {}^t(1, 2, -1) A − 1 t ( 1 , 3 , 1 ) = t ( 2 − 1 , − 3 + 6 − 1 , 1 − 3 + 1 ) = t ( 1 , 2 , − 1 ) で、問題 1.2 の解に一致する。
問題 1.4 ★★ a ∈ R a \in \mathbb{R} a ∈ R を定数とする。次の連立方程式が (i) ただ一つの解をもつ、(ii) 無数の解をもつ、(iii) 解をもたない、ための a a a の条件を求め、解があるときは解をすべて求めよ。
{ a x + y + z = 1 x + a y + z = 1 x + y + a z = 1 \left\lbrace
\begin{aligned}
ax + y + z &= 1 \\
x + ay + z &= 1 \\
x + y + az &= 1
\end{aligned}
\right. ⎩ ⎨ ⎧ a x + y + z x + a y + z x + y + a z = 1 = 1 = 1
解答
第 1 式と第 3 式を入れ替えてから掃き出す。
( 1 1 a 1 1 a 1 1 a 1 1 1 ) → R 3 − a R 1 R 2 − R 1 ( 1 1 a 1 0 a − 1 1 − a 0 0 1 − a 1 − a 2 1 − a ) → R 3 + R 2 ( 1 1 a 1 0 a − 1 1 − a 0 0 0 − ( a − 1 ) ( a + 2 ) − ( a − 1 ) ) \left(\begin{array}{ccc|c} 1 & 1 & a & 1 \\ 1 & a & 1 & 1 \\ a & 1 & 1 & 1 \end{array}\right)
\xrightarrow[R_3 - aR_1]{R_2 - R_1}
\left(\begin{array}{ccc|c} 1 & 1 & a & 1 \\ 0 & a - 1 & 1 - a & 0 \\ 0 & 1 - a & 1 - a^2 & 1 - a \end{array}\right)
\xrightarrow{R_3 + R_2}
\left(\begin{array}{ccc|c} 1 & 1 & a & 1 \\ 0 & a - 1 & 1 - a & 0 \\ 0 & 0 & -(a - 1)(a + 2) & -(a - 1) \end{array}\right) 1 1 a 1 a 1 a 1 1 1 1 1 R 2 − R 1 R 3 − a R 1 1 0 0 1 a − 1 1 − a a 1 − a 1 − a 2 1 0 1 − a R 3 + R 2 1 0 0 1 a − 1 0 a 1 − a − ( a − 1 ) ( a + 2 ) 1 0 − ( a − 1 )
(( 1 − a ) + ( 1 − a 2 ) = 2 − a − a 2 = − ( a − 1 ) ( a + 2 ) (1 - a) + (1 - a^2) = 2 - a - a^2 = -(a-1)(a+2) ( 1 − a ) + ( 1 − a 2 ) = 2 − a − a 2 = − ( a − 1 ) ( a + 2 ) を用いた。)
a = 1 a = 1 a = 1 のとき:第 2, 3 行は 0 になり、rank A = rank ( A ∣ b ) = 1 \operatorname{rank} A = \operatorname{rank}(A \mid b) = 1 rank A = rank ( A ∣ b ) = 1 。解は x = 1 − s − t x = 1 - s - t x = 1 − s − t 、y = s y = s y = s 、z = t z = t z = t ( s , t ∈ R ) (s, t \in \mathbb{R}) ( s , t ∈ R ) で無数にある。
a = − 2 a = -2 a = − 2 のとき:第 3 行は ( 0 0 0 ∣ 3 ) (0 \ 0 \ 0 \mid 3) ( 0 0 0 ∣ 3 ) となり、rank A = 2 < 3 = rank ( A ∣ b ) \operatorname{rank} A = 2 < 3 = \operatorname{rank}(A \mid b) rank A = 2 < 3 = rank ( A ∣ b ) で解はない(実際、3 式を足すと 0 = 3 0 = 3 0 = 3 になる)。
a ≠ 1 , − 2 a \neq 1, -2 a = 1 , − 2 のとき:rank A = 3 \operatorname{rank} A = 3 rank A = 3 で解はただ一つ。第 3 行から z = 1 / ( a + 2 ) z = 1/(a+2) z = 1/ ( a + 2 ) 、第 2 行を a − 1 a - 1 a − 1 で割って y = z y = z y = z 、第 1 行から x = 1 − y − a z = 1 − ( 1 + a ) / ( a + 2 ) = 1 / ( a + 2 ) x = 1 - y - az = 1 - (1 + a)/(a + 2) = 1/(a + 2) x = 1 − y − a z = 1 − ( 1 + a ) / ( a + 2 ) = 1/ ( a + 2 ) 。よって x = y = z = 1 / ( a + 2 ) x = y = z = 1/(a+2) x = y = z = 1/ ( a + 2 ) 。
問題 1.5 ★★ A ∈ M n ( K ) A \in \operatorname{M}_n(K) A ∈ M n ( K ) がすべての B ∈ M n ( K ) B \in \operatorname{M}_n(K) B ∈ M n ( K ) と可換(A B = B A AB = BA A B = B A )ならば、A = c I n A = cI_n A = c I n ( c ∈ K ) (c \in K) ( c ∈ K ) であることを示せ。
解答
行列単位 E i j E_{ij} E ij と可換であることを使う。成分を計算すると
( A E i j ) k l = ∑ p a k p δ p i δ j l = a k i δ j l , ( E i j A ) k l = ∑ p δ k i δ j p a p l = δ k i a j l (AE_{ij})_{kl} = \sum_{p} a_{kp} \delta_{pi} \delta_{jl} = a_{ki} \delta_{jl}, \qquad (E_{ij}A)_{kl} = \sum_{p} \delta_{ki} \delta_{jp} a_{pl} = \delta_{ki} a_{jl} ( A E ij ) k l = p ∑ a k p δ p i δ j l = a k i δ j l , ( E ij A ) k l = p ∑ δ k i δ j p a pl = δ k i a j l
である。l = j l = j l = j とおくと a k i = δ k i a j j a_{ki} = \delta_{ki} a_{jj} a k i = δ k i a j j がすべての i , j , k i, j, k i , j , k で成り立つ。k ≠ i k \neq i k = i なら a k i = 0 a_{ki} = 0 a k i = 0 、k = i k = i k = i なら a i i = a j j a_{ii} = a_{jj} a ii = a j j 。よって A A A は対角成分がすべて等しい対角行列、すなわち A = a 11 I n A = a_{11} I_n A = a 11 I n である。逆にスカラー行列がすべての行列と可換なのは明らかである。
問題 1.6 ★★ 正方行列 A = ( a i j ) A = (a_{ij}) A = ( a ij ) の対角成分の和をトレース (trace) といい tr A \operatorname{tr} A tr A と書く。A ∈ M m , n ( K ) A \in \operatorname{M}_{m,n}(K) A ∈ M m , n ( K ) 、B ∈ M n , m ( K ) B \in \operatorname{M}_{n,m}(K) B ∈ M n , m ( K ) に対し tr ( A B ) = tr ( B A ) \operatorname{tr}(AB) = \operatorname{tr}(BA) tr ( A B ) = tr ( B A ) を示せ。また、A B − B A = I n AB - BA = I_n A B − B A = I n を満たす A , B ∈ M n ( R ) A, B \in \operatorname{M}_n(\mathbb{R}) A , B ∈ M n ( R ) は存在しないことを示せ。
解答
tr ( A B ) = ∑ i = 1 m ∑ j = 1 n a i j b j i = ∑ j = 1 n ∑ i = 1 m b j i a i j = tr ( B A ) \operatorname{tr}(AB) = \sum_{i=1}^{m} \sum_{j=1}^{n} a_{ij} b_{ji} = \sum_{j=1}^{n} \sum_{i=1}^{m} b_{ji} a_{ij} = \operatorname{tr}(BA) tr ( A B ) = ∑ i = 1 m ∑ j = 1 n a ij b j i = ∑ j = 1 n ∑ i = 1 m b j i a ij = tr ( B A ) 。後半:tr \operatorname{tr} tr は線形なので tr ( A B − B A ) = tr ( A B ) − tr ( B A ) = 0 \operatorname{tr}(AB - BA) = \operatorname{tr}(AB) - \operatorname{tr}(BA) = 0 tr ( A B − B A ) = tr ( A B ) − tr ( B A ) = 0 だが、tr I n = n ≠ 0 \operatorname{tr} I_n = n \neq 0 tr I n = n = 0 である。
(注:無限次元では A B − B A = I AB - BA = I A B − B A = I となる線形写像が存在する。例えば多項式の空間で A A A を微分、B B B を x x x 倍とすると ( x p ) ′ − x p ′ = p (xp)' - xp' = p ( x p ) ′ − x p ′ = p 。これは量子力学の正準交換関係の数学的な背景である。)
問題 1.7 ★★ N ∈ M n ( K ) N \in \operatorname{M}_n(K) N ∈ M n ( K ) が N k = O N^k = O N k = O を満たすとき、I n − N I_n - N I n − N は正則で ( I n − N ) − 1 = I n + N + N 2 + ⋯ + N k − 1 (I_n - N)^{-1} = I_n + N + N^2 + \cdots + N^{k-1} ( I n − N ) − 1 = I n + N + N 2 + ⋯ + N k − 1 であることを示せ。これを用いて次の行列の逆行列を求めよ。
U = ( 1 a b 0 1 c 0 0 1 ) U = \begin{pmatrix} 1 & a & b \\ 0 & 1 & c \\ 0 & 0 & 1 \end{pmatrix} U = 1 0 0 a 1 0 b c 1
解答
( I − N ) ( I + N + ⋯ + N k − 1 ) = I − N k = I (I - N)(I + N + \cdots + N^{k-1}) = I - N^k = I ( I − N ) ( I + N + ⋯ + N k − 1 ) = I − N k = I で、逆順の積も同様。U = I + N U = I + N U = I + N とおくと
N = ( 0 a b 0 0 c 0 0 0 ) , N 2 = ( 0 0 a c 0 0 0 0 0 0 ) , N 3 = O N = \begin{pmatrix} 0 & a & b \\ 0 & 0 & c \\ 0 & 0 & 0 \end{pmatrix}, \qquad N^2 = \begin{pmatrix} 0 & 0 & ac \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}, \qquad N^3 = O N = 0 0 0 a 0 0 b c 0 , N 2 = 0 0 0 0 0 0 a c 0 0 , N 3 = O
である。U = I − ( − N ) U = I - (-N) U = I − ( − N ) に前半を適用して
U − 1 = I − N + N 2 = ( 1 − a a c − b 0 1 − c 0 0 1 ) U^{-1} = I - N + N^2 = \begin{pmatrix} 1 & -a & ac - b \\ 0 & 1 & -c \\ 0 & 0 & 1 \end{pmatrix} U − 1 = I − N + N 2 = 1 0 0 − a 1 0 a c − b − c 1
検算:U U − 1 UU^{-1} U U − 1 の第 1 行は ( 1 , − a + a , a c − b − a c + b ) = ( 1 , 0 , 0 ) (1,\ -a + a,\ ac - b - ac + b) = (1, 0, 0) ( 1 , − a + a , a c − b − a c + b ) = ( 1 , 0 , 0 ) 、第 2 行は ( 0 , 1 , − c + c ) = ( 0 , 1 , 0 ) (0, 1, -c + c) = (0, 1, 0) ( 0 , 1 , − c + c ) = ( 0 , 1 , 0 ) 。
問題 1.8 ★★ A ∈ M m , n ( K ) A \in \operatorname{M}_{m,n}(K) A ∈ M m , n ( K ) 、B ∈ M n , m ( K ) B \in \operatorname{M}_{n,m}(K) B ∈ M n , m ( K ) で m > n m > n m > n ならば、A B ∈ M m ( K ) AB \in \operatorname{M}_m(K) A B ∈ M m ( K ) は正則でないことを示せ。
解答
B B B は n × m n \times m n × m 行列で n < m n < m n < m だから、系 1.31 の 2 により B x = 0 Bx = 0 B x = 0 を満たす x ∈ K m x \in K^m x ∈ K m , x ≠ 0 x \neq 0 x = 0 が存在する。このとき ( A B ) x = A ( B x ) = 0 (AB)x = A(Bx) = 0 ( A B ) x = A ( B x ) = 0 であり、定理 1.34 の 5 より A B AB A B は正則でない。(第3章 の言葉では、L A B L_{AB} L A B は n n n 次元の空間 K n K^n K n を経由するので、像の次元が n < m n < m n < m 以下になる。)
問題 1.9 ★★ 命題 1.39 の記号で、A A A と S = D − C A − 1 B S = D - CA^{-1}B S = D − C A − 1 B がともに正則であるとき、
( A B C D ) − 1 = ( A − 1 + A − 1 B S − 1 C A − 1 − A − 1 B S − 1 − S − 1 C A − 1 S − 1 ) \begin{pmatrix} A & B \\ C & D \end{pmatrix}^{-1} = \begin{pmatrix} A^{-1} + A^{-1} B S^{-1} C A^{-1} & -A^{-1} B S^{-1} \\ -S^{-1} C A^{-1} & S^{-1} \end{pmatrix} ( A C B D ) − 1 = ( A − 1 + A − 1 B S − 1 C A − 1 − S − 1 C A − 1 − A − 1 B S − 1 S − 1 )
であることを示せ。
解答
命題 1.39 の証明の記号で、左辺の行列を M = L Δ U M = L \Delta U M = L Δ U と分解すると M − 1 = U − 1 Δ − 1 L − 1 M^{-1} = U^{-1} \Delta^{-1} L^{-1} M − 1 = U − 1 Δ − 1 L − 1 である。例 1.38 より
U − 1 = ( I − A − 1 B O I ) , Δ − 1 = ( A − 1 O O S − 1 ) , L − 1 = ( I O − C A − 1 I ) U^{-1} = \begin{pmatrix} I & -A^{-1}B \\ O & I \end{pmatrix}, \quad \Delta^{-1} = \begin{pmatrix} A^{-1} & O \\ O & S^{-1} \end{pmatrix}, \quad L^{-1} = \begin{pmatrix} I & O \\ -CA^{-1} & I \end{pmatrix} U − 1 = ( I O − A − 1 B I ) , Δ − 1 = ( A − 1 O O S − 1 ) , L − 1 = ( I − C A − 1 O I )
である。まず
Δ − 1 L − 1 = ( A − 1 O − S − 1 C A − 1 S − 1 ) \Delta^{-1} L^{-1} = \begin{pmatrix} A^{-1} & O \\ -S^{-1}CA^{-1} & S^{-1} \end{pmatrix} Δ − 1 L − 1 = ( A − 1 − S − 1 C A − 1 O S − 1 )
であり、これに左から U − 1 U^{-1} U − 1 を掛けると、左上ブロックは A − 1 + A − 1 B S − 1 C A − 1 A^{-1} + A^{-1}BS^{-1}CA^{-1} A − 1 + A − 1 B S − 1 C A − 1 、右上ブロックは − A − 1 B S − 1 -A^{-1}BS^{-1} − A − 1 B S − 1 、下の行のブロックは変わらず、主張の式を得る。
問題 1.10 ★★★ (シャーマン–モリソンの公式)A ∈ GL n ( K ) A \in \operatorname{GL}_n(K) A ∈ GL n ( K ) 、u , v ∈ K n u, v \in K^n u , v ∈ K n とし、α = 1 + t v A − 1 u \alpha = 1 + {}^t v A^{-1} u α = 1 + t v A − 1 u とおく。
α ≠ 0 \alpha \neq 0 α = 0 ならば A + u t v A + u\ {}^t v A + u t v は正則で、( A + u t v ) − 1 = A − 1 − α − 1 A − 1 u t v A − 1 (A + u\ {}^t v)^{-1} = A^{-1} - \alpha^{-1} A^{-1} u\ {}^t v A^{-1} ( A + u t v ) − 1 = A − 1 − α − 1 A − 1 u t v A − 1 であることを示せ。
α = 0 \alpha = 0 α = 0 ならば A + u t v A + u\ {}^t v A + u t v は正則でないことを示せ。
解答
t v A − 1 u {}^t v A^{-1} u t v A − 1 u はスカラー(1 × 1 1 \times 1 1 × 1 行列)であることに注意する。β = t v A − 1 u \beta = {}^t v A^{-1} u β = t v A − 1 u とおくと α = 1 + β \alpha = 1 + \beta α = 1 + β で
( A + u t v ) ( A − 1 − α − 1 A − 1 u t v A − 1 ) = I + u t v A − 1 − α − 1 u t v A − 1 − α − 1 u ( t v A − 1 u ) t v A − 1 = I + ( 1 − α − 1 − α − 1 β ) u t v A − 1 = I + ( 1 − α − 1 ( 1 + β ) ) u t v A − 1 = I \begin{aligned}
(A + u\ {}^t v)(A^{-1} - \alpha^{-1} A^{-1} u\ {}^t v A^{-1})
&= I + u\ {}^t v A^{-1} - \alpha^{-1} u\ {}^t v A^{-1} - \alpha^{-1} u ({}^t v A^{-1} u) {}^t v A^{-1} \\
&= I + \left(1 - \alpha^{-1} - \alpha^{-1}\beta\right) u\ {}^t v A^{-1} = I + \left(1 - \alpha^{-1}(1 + \beta)\right) u\ {}^t v A^{-1} = I
\end{aligned} ( A + u t v ) ( A − 1 − α − 1 A − 1 u t v A − 1 ) = I + u t v A − 1 − α − 1 u t v A − 1 − α − 1 u ( t v A − 1 u ) t v A − 1 = I + ( 1 − α − 1 − α − 1 β ) u t v A − 1 = I + ( 1 − α − 1 ( 1 + β ) ) u t v A − 1 = I
定理 1.34 より右側の行列が逆行列である。
α = 0 \alpha = 0 α = 0 なら β = − 1 \beta = -1 β = − 1 なので u ≠ 0 u \neq 0 u = 0 であり、x = A − 1 u ≠ 0 x = A^{-1}u \neq 0 x = A − 1 u = 0 とおくと ( A + u t v ) x = u + u ( t v A − 1 u ) = u ( 1 + β ) = 0 (A + u\ {}^t v)x = u + u({}^t v A^{-1} u) = u(1 + \beta) = 0 ( A + u t v ) x = u + u ( t v A − 1 u ) = u ( 1 + β ) = 0 。定理 1.34 の 5 より正則でない。
(この公式は、行列の一部だけが変化したときに逆行列を計算し直す手間を省くために、数値計算や統計学で使われる。)