有限實矩陣乘法法則的證明
設 m , n , p , q m,n,p,q m , n , p , q 為正整數,A ∈ R m × n A\in\mathbb R^{m\times n} A ∈ R m × n 、B ∈ R n × p B\in\mathbb R^{n\times p} B ∈ R n × p 、C ∈ R p × q C\in\mathbb R^{p\times q} C ∈ R p × q 。以下只用行乘欄定義與實數的有限求和。對任意輸出位置 ( i , ℓ ) (i,\ell) ( i , ℓ ) ,
( ( A B ) C ) i ℓ = ∑ k = 1 p ( ∑ j = 1 n a i j b j k ) c k ℓ = ∑ j = 1 n a i j ( ∑ k = 1 p b j k c k ℓ ) = ( A ( B C ) ) i ℓ . \begin{aligned}
((AB)C)_{i\ell}
&=\sum_{k=1}^{p}\left(\sum_{j=1}^{n}a_{ij}b_{jk}\right)c_{k\ell}\\
&=\sum_{j=1}^{n}a_{ij}\left(\sum_{k=1}^{p}b_{jk}c_{k\ell}\right)
=(A(BC))_{i\ell}.
\end{aligned} (( A B ) C ) i ℓ = k = 1 ∑ p ( j = 1 ∑ n a ij b j k ) c k ℓ = j = 1 ∑ n a ij ( k = 1 ∑ p b j k c k ℓ ) = ( A ( B C ) ) i ℓ .
中間等號只交換兩個有限和的求和次序,並使用實數分配律。兩邊同為 m × q m\times q m × q 矩陣,每個位置相等,便證明結合律;沒有調換矩陣因子的次序。
對 D ∈ R n × p D\in\mathbb R^{n\times p} D ∈ R n × p ,同一定義給出
[ A ( B + D ) ] i k = ∑ j = 1 n a i j ( b j k + d j k ) = ( A B ) i k + ( A D ) i k . [A(B+D)]_{ik}=\sum_{j=1}^{n}a_{ij}(b_{jk}+d_{jk})
=(AB)_{ik}+(AD)_{ik}. [ A ( B + D ) ] ik = j = 1 ∑ n a ij ( b j k + d j k ) = ( A B ) ik + ( A D ) ik .
對 E ∈ R m × n E\in\mathbb R^{m\times n} E ∈ R m × n ,則有
[ ( A + E ) B ] i k = ∑ j = 1 n ( a i j + e i j ) b j k = ( A B ) i k + ( E B ) i k . [(A+E)B]_{ik}=\sum_{j=1}^{n}(a_{ij}+e_{ij})b_{jk}
=(AB)_{ik}+(EB)_{ik}. [( A + E ) B ] ik = j = 1 ∑ n ( a ij + e ij ) b j k = ( A B ) ik + ( E B ) ik .
所以在所寫的相容大小下,兩條分配律均成立。對任意實純量 c c c ,
[ ( c A ) B ] i k = ∑ j ( c a i j ) b j k = c ∑ j a i j b j k = ∑ j a i j ( c b j k ) = [ A ( c B ) ] i k , [(cA)B]_{ik}=\sum_j(ca_{ij})b_{jk}
=c\sum_j a_{ij}b_{jk}
=\sum_j a_{ij}(cb_{jk})=[A(cB)]_{ik}, [( c A ) B ] ik = j ∑ ( c a ij ) b j k = c j ∑ a ij b j k = j ∑ a ij ( c b j k ) = [ A ( c B ) ] ik ,
逐項相等便證明 ( c A ) B = c ( A B ) = A ( c B ) (cA)B=c(AB)=A(cB) ( c A ) B = c ( A B ) = A ( c B ) 。特別地,對 x , y ∈ R n x,y\in\mathbb R^n x , y ∈ R n 與實數 s , t s,t s , t ,
[ A ( s x + t y ) ] i = ∑ j a i j ( s x j + t y j ) = s ( A x ) i + t ( A y ) i . [A(sx+ty)]_i=\sum_j a_{ij}(sx_j+ty_j)
=s(Ax)_i+t(Ay)_i. [ A ( s x + t y ) ] i = j ∑ a ij ( s x j + t y j ) = s ( A x ) i + t ( A y ) i .
這證明了比較同一線性系統的解時所用的矩陣乘向量線性性質。所有恆等式都只針對大小相容的有限實矩陣,不聲稱任意矩陣因子可交換或可消去。
由單位矩陣建立初等矩陣
假設一個初等 行變換 ρ \rho ρ 作用於有 p p p 行的矩陣。先從單位矩陣 I p I_p I p
開始,對 I p I_p I p 做同一個初等行變換,所得矩陣記作 E ρ E_\rho E ρ 。
所有純量都屬於矩陣所採用的底域;特別地,每個非零倍數 c c c 在該域中都有倒數
1 / c 1/c 1/ c 。
定義
初等行變換矩陣 與 p p p 行矩陣上的初等行變換 ρ \rho ρ 對應的初等行變換矩陣 ,就是把
ρ \rho ρ 作用在 I p I_p I p 上所得的 p × p p \times p p × p 矩陣。
這種矩陣亦簡稱初等矩陣 ;較長的名稱可以明確指出它與指定行變換的關係。
換句話說,
E ρ = ρ ( I p ) . E_\rho = \rho(I_p). E ρ = ρ ( I p ) .
這個定義有用,是因為以下定理。
定理
初等行變換等同於左乘矩陣 設 A A A 是任何有 p p p 行的矩陣,ρ \rho ρ 是一個初等行變換,並以
ρ ( A ) \rho(A) ρ ( A ) 表示把它作用於 A A A 後得到的矩陣。則
ρ ( A ) = E ρ A . \rho(A) = E_\rho A. ρ ( A ) = E ρ A . 所以,對 A A A 做這個初等行變換,等同於用同一個操作在 I p I_p I p 上得到的
初等行變換矩陣左乘 A A A 。
這裏必須是左乘。初等行變換改變的是行,會把新行寫成舊行的線性組合;左乘正是把
A A A 的新行寫成舊行的線性組合。右乘則會混合欄,不是同一件事。
維數亦說明了同一點。若 A A A 是 p × n p \times n p × n 矩陣,則 E ρ E_\rho E ρ 是
p × p p \times p p × p 矩陣,所以 E ρ A E_\rho A E ρ A 有定義,而且仍是 p × n p \times n p × n 。
這裏完全不要求 A A A 是方陣。相反,對這 n n n 欄進行欄變換的矩陣大小為 n × n n \times n n × n ,並從右邊相乘。
為甚麼左乘會執行指定行變換
證明
對三類初等行變換逐行驗證 把任意 p × n p \times n p × n 矩陣 A A A 的各行記作
a 1 , a 2 , … , a p a_1,a_2,\ldots,a_p a 1 , a 2 , … , a p ,並以 e k T e_k^{\mathsf T} e k T 表示 I p I_p I p 的第 k k k 行。
基本等式
e k T A = a k e_k^{\mathsf T}A=a_k e k T A = a k 表示左乘矩陣的一行會選取 A A A 的相應行。現在逐類檢查。
對 R i ↔ R j R_i\leftrightarrow R_j R i ↔ R j ,E ρ E_\rho E ρ 的第 i i i 行是
e j T e_j^{\mathsf T} e j T ,第 j j j 行是 e i T e_i^{\mathsf T} e i T ,其餘第 k k k 行仍是
e k T e_k^{\mathsf T} e k T 。因此 E ρ A E_\rho A E ρ A 恰好交換 A A A 的第 i , j i,j i , j 行。
對 R i ← c R i R_i\leftarrow cR_i R i ← c R i (c ≠ 0 c\ne0 c = 0 ),E ρ E_\rho E ρ 的第 i i i 行是
c e i T c e_i^{\mathsf T} c e i T ,其餘行不變。因此乘積的第 i i i 行為
( c e i T ) A = c a i , (c e_i^{\mathsf T})A=c a_i, ( c e i T ) A = c a i ,
而每個 k ≠ i k\ne i k = i 的第 k k k 行仍為 a k a_k a k 。
對 R j ← R j + c R i R_j\leftarrow R_j+cR_i R j ← R j + c R i (i ≠ j i\ne j i = j ),E ρ E_\rho E ρ 的第 j j j 行是
e j T + c e i T e_j^{\mathsf T}+c e_i^{\mathsf T} e j T + c e i T ,其餘行不變,所以新的第 j j j 行為
( e j T + c e i T ) A = a j + c a i , (e_j^{\mathsf T}+c e_i^{\mathsf T})A=a_j+c a_i, ( e j T + c e i T ) A = a j + c a i ,
正是指定的行加法。
三類情形窮盡了初等行變換,因此對每個相容矩陣 A A A 都有
E ρ A = ρ ( A ) E_\rho A=\rho(A) E ρ A = ρ ( A ) 。這個證明沒有聲稱任意改動行的方式都對應一個初等矩陣;
ρ \rho ρ 必須是剛才檢查的三類操作之一。
作用於增廣矩陣
假設 A x = b Ax=b A x = b 由增廣矩陣 [ A ∣ b ] [A\mid b] [ A ∣ b ] 表示。初等行變換必須作用於整個增廣
矩陣,不能只改其中一邊。相應的乘法等式是
E ρ [ A ∣ b ] = [ E ρ A ∣ E ρ b ] . E_\rho[A\mid b]=[E_\rho A\mid E_\rho b]. E ρ [ A ∣ b ] = [ E ρ A ∣ E ρ b ] .
因此,係數和右端常數會接受完全相同的行組合。向量 x x x 滿足 A x = b Ax=b A x = b ,若且唯若
它滿足 E ρ A x = E ρ b E_\rho Ax=E_\rho b E ρ A x = E ρ b :正向可把原等式左乘 E ρ E_\rho E ρ ,反向則施加逆向
初等行變換。這才是解集保持不變的嚴格理由。
若只改 A A A 而固定 b b b ,或只改 b b b ,一般會得到另一個方程組。增廣線只是視覺
上的分隔,不是行變換可以中途停止的邊界。
一串初等行變換就是一個矩陣乘積
真正的用途不只是表示一個初等步驟,而是把整串初等行變換寫成一個矩陣乘積。
定理
一串初等行變換的乘積表示 假設
A 1 → ρ 1 A 2 → ρ 2 A 3 → ρ 3 ⋯ → ρ k A k + 1 . A_1 \xrightarrow{\rho_1} A_2 \xrightarrow{\rho_2} A_3
\xrightarrow{\rho_3} \cdots \xrightarrow{\rho_k} A_{k+1}. A 1 ρ 1 A 2 ρ 2 A 3 ρ 3 ⋯ ρ k A k + 1 . 則
A k + 1 = E ρ k E ρ k − 1 ⋯ E ρ 2 E ρ 1 A 1 . A_{k+1}
= E_{\rho_k}E_{\rho_{k-1}}\cdots E_{\rho_2}E_{\rho_1}A_1. A k + 1 = E ρ k E ρ k − 1 ⋯ E ρ 2 E ρ 1 A 1 . 這裏每個 ρ j \rho_j ρ j 都必須是作用於同樣行數矩陣的初等行變換。這個乘積表示整串
操作,但乘積本身未必仍是單一個初等行變換矩陣。
這個公式中的次序很重要。第一個行變換最靠近 A 1 A_1 A 1 ,因為它最先作用:
A 2 = E ρ 1 A 1 , A 3 = E ρ 2 A 2 = E ρ 2 E ρ 1 A 1 . A_2 = E_{\rho_1}A_1,
\qquad
A_3 = E_{\rho_2}A_2 = E_{\rho_2}E_{\rho_1}A_1. A 2 = E ρ 1 A 1 , A 3 = E ρ 2 A 2 = E ρ 2 E ρ 1 A 1 .
例題
合併兩個行變換 令
A = [ 1 0 1 0 2 1 1 0 0 ] . A =
\begin{bmatrix}
1 & 0 & 1\\
0 & 2 & 1\\
1 & 0 & 0
\end{bmatrix}. A = 1 0 1 0 2 0 1 1 0 . 依次做
ρ 1 : R 2 ← R 2 + R 1 , ρ 2 : R 1 ← R 1 + 2 R 2 . \rho_1:\ R_2 \leftarrow R_2 + R_1,
\qquad
\rho_2:\ R_1 \leftarrow R_1 + 2R_2. ρ 1 : R 2 ← R 2 + R 1 , ρ 2 : R 1 ← R 1 + 2 R 2 . 對應的行變換矩陣是
E ρ 1 = [ 1 0 0 1 1 0 0 0 1 ] , E ρ 2 = [ 1 2 0 0 1 0 0 0 1 ] . E_{\rho_1} =
\begin{bmatrix}
1 & 0 & 0\\
1 & 1 & 0\\
0 & 0 & 1
\end{bmatrix},
\qquad
E_{\rho_2} =
\begin{bmatrix}
1 & 2 & 0\\
0 & 1 & 0\\
0 & 0 & 1
\end{bmatrix}. E ρ 1 = 1 1 0 0 1 0 0 0 1 , E ρ 2 = 1 0 0 2 1 0 0 0 1 . 兩個操作完成後,結果就是
E ρ 2 E ρ 1 A = [ 3 2 0 1 1 0 0 0 1 ] [ 1 0 1 0 2 1 1 0 0 ] = [ 3 4 5 1 2 2 1 0 0 ] . E_{\rho_2}E_{\rho_1}A
=
\begin{bmatrix}
3&2&0\\
1&1&0\\
0&0&1
\end{bmatrix}
\begin{bmatrix}
1&0&1\\
0&2&1\\
1&0&0
\end{bmatrix}
=
\begin{bmatrix}
3&4&5\\
1&2&2\\
1&0&0
\end{bmatrix}. E ρ 2 E ρ 1 A = 3 1 0 2 1 0 0 0 1 1 0 1 0 2 0 1 1 0 = 3 1 1 4 2 0 5 2 0 . 若不想直接相乘 E ρ 2 E ρ 1 E_{\rho_2}E_{\rho_1} E ρ 2 E ρ 1 ,也可以把同一串行變換按同一次序作用在
I 3 I_3 I 3 上,所得矩陣就是合併後的左乘矩陣。最後的數值結果亦確認:第二步使用的
是已經更新後的第 2 2 2 行。
反向行變換與逆矩陣
每個初等行變換都有反向操作:
R j ← R j + c R i R_j \leftarrow R_j + cR_i R j ← R j + c R i 的反向操作是
R j ← R j − c R i R_j \leftarrow R_j - cR_i R j ← R j − c R i ;
R i ← c R i R_i \leftarrow cR_i R i ← c R i (其中 c ≠ 0 c \ne 0 c = 0 )的反向操作是
R i ← ( 1 / c ) R i R_i \leftarrow (1/c)R_i R i ← ( 1/ c ) R i ;
交換兩行的反向操作仍然是同一次交換。
這給出一個精確的矩陣敘述。
定理
初等行變換矩陣都是可逆矩陣 每個初等行變換矩陣都是可逆的。它的逆矩陣,就是反向初等行變換所對應的
初等行變換矩陣。
證明
為甚麼反向矩陣是雙側逆矩陣 設 E ρ E_\rho E ρ 表示初等行變換 ρ \rho ρ ,而 E ρ − 1 E_{\rho^{-1}} E ρ − 1 表示它的反向初等
操作。對每個有 p p p 行的矩陣 X X X ,先做 ρ \rho ρ 再撤銷,都會恢復 X X X 。由左乘
定理,
E ρ − 1 E ρ X = X . E_{\rho^{-1}}E_\rho X=X. E ρ − 1 E ρ X = X . 這個等式對每個這樣的 X X X 都成立;特別取 X = I p X=I_p X = I p ,得到
E ρ − 1 E ρ = I p E_{\rho^{-1}}E_\rho=I_p E ρ − 1 E ρ = I p 。先做反向操作再做 ρ \rho ρ ,同理得到
E ρ E ρ − 1 = I p E_\rho E_{\rho^{-1}}=I_p E ρ E ρ − 1 = I p 。因此
E ρ − 1 = E ρ − 1 . E_\rho^{-1}=E_{\rho^{-1}}. E ρ − 1 = E ρ − 1 . 具體來說,交換兩行的逆仍是同一次交換;乘以 c ≠ 0 c\ne0 c = 0 的逆是乘以 1 / c 1/c 1/ c ;
把第 i i i 行的 c c c 倍加到第 j j j 行,則由加上 − c -c − c 倍第 i i i 行來撤銷。
例如,若
E = [ 1 0 0 3 1 0 0 0 1 ] E =
\begin{bmatrix}
1 & 0 & 0\\
3 & 1 & 0\\
0 & 0 & 1
\end{bmatrix} E = 1 3 0 0 1 0 0 0 1
表示 R 2 ← R 2 + 3 R 1 R_2 \leftarrow R_2 + 3R_1 R 2 ← R 2 + 3 R 1 ,則
E − 1 = [ 1 0 0 − 3 1 0 0 0 1 ] E^{-1} =
\begin{bmatrix}
1 & 0 & 0\\
-3 & 1 & 0\\
0 & 0 & 1
\end{bmatrix} E − 1 = 1 − 3 0 0 1 0 0 0 1
表示 R 2 ← R 2 − 3 R 1 R_2 \leftarrow R_2 - 3R_1 R 2 ← R 2 − 3 R 1 。
這就是初等行變換可逆的代數原因,也解釋了為甚麼行化簡和可逆矩陣有密切關係。
證明
先撤銷最後一步操作 設 E 1 , … , E k E_1,\ldots,E_k E 1 , … , E k 按時間順序編碼同一組 p p p 行矩陣上的操作,合併左乘矩陣就是 J = E k ⋯ E 1 J=E_k\cdots E_1 J = E k ⋯ E 1 。每個 E j − 1 E_j^{-1} E j − 1 表示對應初等行變換的反向操作。要恢復輸入,首先必須撤銷 E k E_k E k ,因為它是最後一次對行所作的改變;然後依次撤銷 E k − 1 , … , E 1 E_{k-1},\ldots,E_1 E k − 1 , … , E 1 。因此,反向鏈的合併矩陣為
K = E 1 − 1 E 2 − 1 ⋯ E k − 1 . K=E_1^{-1}E_2^{-1}\cdots E_k^{-1}. K = E 1 − 1 E 2 − 1 ⋯ E k − 1 . 這裏必須區分兩種次序:反向操作的時間順序從 E k − 1 E_k^{-1} E k − 1 開始;寫成乘積時,最先作用的因子卻放在最靠近待恢復矩陣的位置。以三步為例,由結合律可得
K J = E 1 − 1 E 2 − 1 ( E 3 − 1 E 3 ) E 2 E 1 = E 1 − 1 ( E 2 − 1 E 2 ) E 1 = I p . KJ=E_1^{-1}E_2^{-1}(E_3^{-1}E_3)E_2E_1
=E_1^{-1}(E_2^{-1}E_2)E_1=I_p. K J = E 1 − 1 E 2 − 1 ( E 3 − 1 E 3 ) E 2 E 1 = E 1 − 1 ( E 2 − 1 E 2 ) E 1 = I p . 另一個乘積 J K JK J K 則從相鄰的 E 1 E 1 − 1 E_1E_1^{-1} E 1 E 1 − 1 開始消去,最後同樣得到 I p I_p I p 。任意有限鏈都可這樣逐對消去。整個論證沒有交換因子,也不要求被變換的矩陣本身是方陣。
把這個讀法用於前面的兩步例子:合併矩陣是 E ρ 2 E ρ 1 E_{\rho_2}E_{\rho_1} E ρ 2 E ρ 1 ,所以逆矩陣是 E ρ 1 − 1 E ρ 2 − 1 E_{\rho_1}^{-1}E_{\rho_2}^{-1} E ρ 1 − 1 E ρ 2 − 1 。從最終矩陣開始,先從第一行減去目前第二行的兩倍,第一行便恢復為 [ 1 0 1 ] [1\ 0\ 1] [ 1 0 1 ] ;再從第二行減去剛恢復的第一行,第二行便恢復為 [ 0 2 1 ] [0\ 2\ 1] [ 0 2 1 ] 。撤銷操作是在逐個恢復中間狀態,乘積次序也就隨之確定,無須另背一個口訣。
有限消元必停止,並保持解集
設 A ∈ R m × n A\in\mathbb R^{m\times n} A ∈ R m × n 、b ∈ R m b\in\mathbb R^m b ∈ R m ,其中 m , n m,n m , n 是有限正整數。允許的行變換為交換不同的兩行、以非零標量乘一行,以及把一行的標量倍加到另一行。逆變換分別為同一交換、乘以倒數,以及減去該倍數。因此對整個增廣陣施行變換恰好保持 A x = b Ax=b A x = b 的解集:新方程由舊方程推出,逆變換給出反向蘊含。每個行變換也等於左乘把同一變換施於單位矩陣所得的可逆初等矩陣。
從左到右掃描係數欄,並記錄下一未使用行。若當前欄在未使用行中有非零元素,就將其所在行交換到下一位置,除以非零主元,再從其餘各行減去適當倍數,清空該欄其他元素,然後推進未使用行指標。若該欄所有未使用位置都是零,就跳過該欄。先前主元欄始終是單位欄:當前未使用行在這些欄均為零,所以新主元行以及後續減法保留這些零與原主元。總共只有 n n n 次欄掃描及至多 min ( m , n ) \min(m,n) min ( m , n ) 個主元,每次只作有限個行變換,故算法必停止。最後非零係數行的主元欄依次遞增,各主元上下均為零,其餘係數行全零。這證明下面使用的簡化階梯形確實能得到,而非假設它最終出現。
解的分類與齊次自由度
令 r r r 為係數部分的主元數。若某個零係數行對應非零變換後右端,就出現不可能的方程 0 = c 0=c 0 = c ,因此無解。否則任意指定 n − r n-r n − r 個非主元變量的實數值。每個主元行都把其主元變量寫成變換後右端減去自由變量的已知線性組合,其中沒有其他主元變量。這些公式既構造解,也恢復每一個解。若 r = n r=n r = n ,解唯一;若 r < n r<n r < n ,改變自由變量就改變解向量,故有無窮多解。
對 A z = 0 Az=0 A z = 0 ,全部自由變量取零給出零解。若 r < n r<n r < n ,將一個自由變量設為一,其餘設為零,再按主元公式求解,就構造出非零齊次解;若 r = n r=n r = n ,公式迫使所有齊次變量為零。對相容系統固定一個解 x p x_p x p ,兩解相減表明每個解均為 x p + z x_p+z x p + z ,其中 A z = 0 Az=0 A z = 0 ;代入又證明反向包含。
當 A A A 為方陣且每欄都有主元時,係數階梯形為 I I I 。有限個初等矩陣之積 E E E 滿足 E A = I EA=I E A = I 。因為 E E E 可逆,A = E − 1 A=E^{-1} A = E − 1 ,故也有 A E = I AE=I A E = I ,所以 E E E 是雙側逆。若缺少主元,上面構造的非零齊次解排除可逆性。這也證明可逆時,對 [ A ∣ I ] [A\mid I] [ A ∣ I ] 作相同行變換會得到 [ I ∣ A − 1 ] [I\mid A^{-1}] [ I ∣ A − 1 ] 。
不預設計數結論的有限基存在與擴充
設實向量空間 V V V 已給定有限生成列表 t 1 , … , t N t_1,\ldots,t_N t 1 , … , t N 。先證明任何線性無關列表 u 1 , … , u k u_1,\ldots,u_k u 1 , … , u k 的長度不超過 N N N 。插入 j j j 個向量後,維持一個由 u 1 , … , u j u_1,\ldots,u_j u 1 , … , u j 與 N − j N-j N − j 個剩餘舊向量組成的生成列表。將 u j + 1 u_{j+1} u j + 1 按此列表展開,至少一個剩餘舊向量的係數非零;否則它在先前 u i u_i u i 的張成空間中,與無關性矛盾。解出這個舊向量,再以 u j + 1 u_{j+1} u j + 1 取代它。新列表仍張成同一空間,而且保留全部先前插入向量。若 k > N k>N k > N ,則 N N N 步後前 N N N 個向量已張成 V V V ,與 u N + 1 u_{N+1} u N + 1 的無關性矛盾。此即所需上界,也包括 N = 0 N=0 N = 0 。
每個子空間 W W W 都有有限基。從空的無關列表開始;若當前張成空間尚不是 W W W ,就在 W W W 中取一個不在該張成空間內的向量。若新向量在某關係中的係數非零,就能由舊向量表示它,矛盾,故新列表仍無關。已證上界允許至多添加 N N N 次,所以過程必停止於 W W W 的無關生成列表。若 W = { 0 } W=\{0\} W = { 0 } ,立即得到空基。改從 W W W 內任意給定無關列表開始,同樣證明該列表可擴充為 W W W 的基;只需有界有限次選擇。最後,若兩組基的長度為 p , q p,q p , q ,分別把它們用作上界證明中的有限生成列表,得到 p ≤ q p\le q p ≤ q 與 q ≤ p q\le p q ≤ p 。故基的長度有良定義,稱為維數。這些論證補足秩–零度定理所用的基存在、擴充與維數事實。
方陣可能把不同輸入送到同一個輸出,也可能無法產生某些輸出。逆矩陣要完全撤銷原矩陣的作用,就必須保證每個輸出都恰好對應一個輸入。因此,可逆性與方程組 A x = b Ax=b A x = b 的解集直接相關。除非另有說明,矩陣元素均取實數。
開始之前
矩陣乘法必須保持因子次序,每一步合法的初等行變換都可逆。方程組 A x = b Ax=b A x = b 有唯一解,當且僅當它一致,且化簡後每個變量欄都有主元。可逆性將保證每個維數適當的右端都滿足這個條件。
左逆與右逆
定義
左逆與右逆 設 A A A 是一個 p × q p \times q p × q 矩陣。
一個 q × p q \times p q × p 矩陣 H H H 如果滿足 H A = I q HA = I_q H A = I q ,就稱為 A A A 的左逆 。
一個 q × p q \times p q × p 矩陣 G G G 如果滿足 A G = I p AG = I_p A G = I p ,就稱為 A A A 的右逆 。
對長方形矩陣,左逆或右逆可以單獨存在;但若兩者都存在,它們必定相同。方陣還具有更強的性質:只要任何一側的逆存在,另一側的逆也存在。
定義
可逆矩陣 若 A A A 是一個 p × p p \times p p × p 方陣,而且存在另一個 p × p p \times p p × p 矩陣 B B B
使得
B A = A B = I p , BA = AB = I_p, B A = A B = I p , 我們就稱 A A A 為可逆 。矩陣 B B B 叫做 A A A 的逆矩陣 ,記作 A − 1 A^{-1} A − 1 。
定理
逆矩陣是唯一的 設 A A A 是 p × q p\times q p × q 矩陣,B , C B,C B , C 是 q × p q\times p q × p 矩陣。若 B B B 是 A A A 的左逆,C C C 是右逆,則 B = C B=C B = C 。特別地,可逆方陣恰好有一個逆矩陣。
證明
為何逆矩陣唯一 假設給出 B A = I q BA=I_q B A = I q 與 A C = I p AC=I_p A C = I p 。由結合律,
B = B I p = B ( A C ) = ( B A ) C = I q C = C . B=BI_p=B(AC)=(BA)C=I_qC=C. B = B I p = B ( A C ) = ( B A ) C = I q C = C . 這裏沒有交換任何因子。可逆方陣的任何兩個逆矩陣都滿足這些假設,因此必定相同。
甚麼叫可逆
逆矩陣的定義要求恢復每個輸入向量。對角矩陣尤其容易顯示這種作用。
例題
對角矩陣最容易看出逆矩陣 設
D = diag ( 2 , − 1 , 3 ) . D = \operatorname{diag}(2, -1, 3). D = diag ( 2 , − 1 , 3 ) . 那麼
D − 1 = diag ( 1 2 , − 1 , 1 3 ) . D^{-1} = \operatorname{diag}\left(\tfrac{1}{2}, -1, \tfrac{1}{3}\right). D − 1 = diag ( 2 1 , − 1 , 3 1 ) . 原因很直接:對角線上每個非零元素變成倒數,而其餘零位仍然是零。將 D D D 與 D − 1 D^{-1} D − 1 相乘,就會得到 I 3 I_3 I 3 。
等價表述
定理
識別可逆矩陣的等價條件 設 A A A 是實 p × p p\times p p × p 矩陣。以下條件等價:
A A A 可逆。
A A A 與 I p I_p I p 行等價。
A A A 是若干個 p p p 階初等行變換矩陣的乘積。
A A A 有左逆。
A A A 有右逆。
A A A 非奇異,即 A x = 0 Ax=0 A x = 0 只有零解。
對每個 b ∈ R p b\in\mathbb R^p b ∈ R p ,方程組 A x = b Ax=b A x = b 都一致。
對每個 b ∈ R p b\in\mathbb R^p b ∈ R p ,方程組 A x = b Ax=b A x = b 都恰好有一個解。
等價地,A A A 的欄向量線性獨立,或它們張成 R p \mathbb R^p R p 。這些條件成立時,唯一解為 A − 1 b A^{-1}b A − 1 b 。
證明
方陣假設在等價性證明中的作用 若 H H H 是左逆,則 A x = 0 Ax=0 A x = 0 蘊含 x = H A x = 0 x=HAx=0 x = H A x = 0 。反過來,設齊次系統只有零解。在 A A A 的任何最簡行階梯形中,自由變量都會產生一個非零齊次解,因此每欄都有主元。由於 A A A 是方陣,恰有 p p p 個主元行,所得最簡行階梯形必為 I p I_p I p 。
若行變換乘積 E E E 滿足 E A = I p EA=I_p E A = I p ,反向鏈說明 E E E 可逆,而且 A = E − 1 A=E^{-1} A = E − 1 。因此 A A A 可逆,並且本身是初等矩陣的乘積。反過來,這類乘積可逆,而逆矩陣尤其是左逆。這就連接了條件 1–4 和 6,並沒有假設 RREF 唯一。
現設 A G = I p AG=I_p A G = I p 。此時 A A A 是方陣 G G G 的左逆,剛證明的左逆判別使 G G G 可逆,從而 A = A ( G G − 1 ) = ( A G ) G − 1 = G − 1 A=A(GG^{-1})=(AG)G^{-1}=G^{-1} A = A ( G G − 1 ) = ( A G ) G − 1 = G − 1 可逆。這就證明了右逆判別,沒有擅自顛倒乘積。
若每個 A x = b Ax=b A x = b 都一致,對 I p I_p I p 的第 j j j 欄 e j e_j e j 選擇解 u j u_j u j ,使 A u j = e j Au_j=e_j A u j = e j 。令 U = [ u 1 ⋯ u p ] U=[u_1\ \cdots\ u_p] U = [ u 1 ⋯ u p ] ,便有 A U = I p AU=I_p A U = I p ,可用右逆判別。最後,若 A A A 可逆,則 A ( A − 1 b ) = b A(A^{-1}b)=b A ( A − 1 b ) = b 證明解存在;任何解 x x x 都滿足 x = A − 1 A x = A − 1 b x=A^{-1}Ax=A^{-1}b x = A − 1 A x = A − 1 b ,所以解唯一。每個右端都有唯一解,自然蘊含每個右端都有解。
欄向量的說法只是翻譯:線性獨立等價於 A x = 0 Ax=0 A x = 0 只有零系數向量;張成則表示每個 b b b 都有表示 A x = b Ax=b A x = b 。這裏沒有額外增加欄向量的假設。
行化簡與逆矩陣
每個初等行變換矩陣都可逆,逆矩陣編碼對應的反向操作,證明見 3.5 初等行變換矩陣 。若按時間順序施加的操作對應 E 1 , … , E k E_1,\ldots,E_k E 1 , … , E k ,合併左乘矩陣就是 E = E k ⋯ E 1 E=E_k\cdots E_1 E = E k ⋯ E 1 。
每一步都須作用於整個增廣矩陣,因此
E [ A ∣ I p ] = [ E A ∣ E ] . E[A\mid I_p]=[EA\mid E]. E [ A ∣ I p ] = [ E A ∣ E ] .
若左塊變成 I p I_p I p ,則 E A = I p EA=I_p E A = I p 。由於 E E E 可逆,有 A = E − 1 A=E^{-1} A = E − 1 ,右塊 E E E 正是 A − 1 A^{-1} A − 1 。若左塊化簡後缺少主元,方陣便不滿足上面的等價條件。這一方法同時判斷可逆性和計算逆矩陣;右塊再複雜,也不能補救左塊缺少主元的問題。
如何使用可逆性的等價字典
按已知資料選擇判別:主元、單側逆等式、非零零空間向量,或所有右端的可解性。
例題
非零零空間向量可以證明不可逆 設 A A A 是方陣,並且
v = [ 1 − 2 1 ] ≠ 0 , A v = 0. v =
\begin{bmatrix}
1 \\ -2 \\ 1
\end{bmatrix}
\ne 0,
\qquad
Av = 0. v = 1 − 2 1 = 0 , A v = 0. 那麼 A A A 不可能可逆。理由是:若 A − 1 A^{-1} A − 1 存在,把 A v = 0 Av = 0 A v = 0 左乘
A − 1 A^{-1} A − 1 ,就會得到
A − 1 A v = A − 1 0. A^{-1}Av = A^{-1}0. A − 1 A v = A − 1 0. 左邊等於 v v v ,右邊等於 0 0 0 ,於是會推出 v = 0 v = 0 v = 0 。這和已知的非零向量矛盾。
因此,只要找到一個非零的 A x = 0 Ax = 0 A x = 0 解,就已經證明 A A A 是奇異矩陣。
例題
用逆矩陣解方程組 設
A − 1 = [ 2 − 1 0 3 ] , b = [ 5 − 1 ] . A^{-1} =
\begin{bmatrix}
2 & -1 \\
0 & 3
\end{bmatrix},
\qquad
b =
\begin{bmatrix}
5 \\ -1
\end{bmatrix}. A − 1 = [ 2 0 − 1 3 ] , b = [ 5 − 1 ] . 若 A x = b Ax = b A x = b ,兩邊左乘 A − 1 A^{-1} A − 1 :
x = A − 1 b = [ 2 − 1 0 3 ] [ 5 − 1 ] = [ 11 − 3 ] . x = A^{-1}b
=
\begin{bmatrix}
2 & -1 \\
0 & 3
\end{bmatrix}
\begin{bmatrix}
5 \\ -1
\end{bmatrix}
=
\begin{bmatrix}
11 \\ -3
\end{bmatrix}. x = A − 1 b = [ 2 0 − 1 3 ] [ 5 − 1 ] = [ 11 − 3 ] . 逆矩陣不只是「把 A A A 反過來」的符號,它給出了唯一解的明確公式。
由可逆矩陣角度理解行等價
定理
行等價等同於左乘一個可逆矩陣 設 A A A 與 B B B 是有 p p p 行的矩陣。以下敘述等價:
A A A 與 B B B 行等價;
存在一個可逆的 p × p p \times p p × p 矩陣 G G G ,使得
B = G A . B = GA. B = G A . 而且一旦 B = G A B = GA B = G A ,就同時有
A = G − 1 B . A = G^{-1}B. A = G − 1 B .
行變換鏈給出這樣的乘積 G G G 。反過來,等價字典把任何可逆的 G G G 表成初等矩陣的乘積,所以 G A GA G A 可由相應行變換從 A A A 得到。左乘 G − 1 G^{-1} G − 1 就能恢復 A A A 。
例題
把行等價讀成一條矩陣等式 設
A = [ 1 0 1 1 2 3 0 1 1 ] , G = [ 1 0 0 − 1 1 0 0 0 1 ] . A = \begin{bmatrix}
1 & 0 & 1 \\
1 & 2 & 3 \\
0 & 1 & 1
\end{bmatrix},
\qquad
G = \begin{bmatrix}
1 & 0 & 0 \\
-1 & 1 & 0 \\
0 & 0 & 1
\end{bmatrix}. A = 1 1 0 0 2 1 1 3 1 , G = 1 − 1 0 0 1 0 0 0 1 . 矩陣 G G G 就是行變換
R 2 ← R 2 − R 1 R_2 \leftarrow R_2 - R_1 R 2 ← R 2 − R 1 對應的行變換矩陣。所以
G A = [ 1 0 1 0 2 2 0 1 1 ] . GA
=
\begin{bmatrix}
1 & 0 & 1 \\
0 & 2 & 2 \\
0 & 1 & 1
\end{bmatrix}. G A = 1 0 0 0 2 1 1 2 1 . 如果把這個新矩陣記作 B B B ,就有 B = G A B = GA B = G A 。一條等式已經完整記錄了這次
行變換。因為 G G G 可逆,所以 A A A 與 B B B 行等價。
定理
行變換保留對應欄向量之間的線性關係 設 A A A 與 B B B 是行等價的 p × q p \times q p × q 矩陣,並把它們的欄寫成
A = [ a 1 a 2 ⋯ a q ] , B = [ b 1 b 2 ⋯ b q ] . A = [a_1 \; a_2 \; \cdots \; a_q],
\qquad
B = [b_1 \; b_2 \; \cdots \; b_q]. A = [ a 1 a 2 ⋯ a q ] , B = [ b 1 b 2 ⋯ b q ] . 如果
a j = α 1 a k 1 + α 2 a k 2 + ⋯ + α n a k n , a_j = \alpha_1 a_{k_1} + \alpha_2 a_{k_2} + \cdots + \alpha_n a_{k_n}, a j = α 1 a k 1 + α 2 a k 2 + ⋯ + α n a k n , 就必然有
b j = α 1 b k 1 + α 2 b k 2 + ⋯ + α n b k n . b_j = \alpha_1 b_{k_1} + \alpha_2 b_{k_2} + \cdots + \alpha_n b_{k_n}. b j = α 1 b k 1 + α 2 b k 2 + ⋯ + α n b k n . 特別地,對應欄之間的線性相依與線性獨立,都會被行等價保留。
左乘 G G G 把關係向前傳遞,左乘 G − 1 G^{-1} G − 1 則恢復原關係。因此同一組係數關係在兩個方向都成立,也就證明了線性相依與線性獨立均被保留。欄向量本身通常會改變。完整的 RREF 唯一性證明及其與秩的關係,見 5.2 RREF 的唯一性與秩 。
乘積與逆矩陣的次序
定理
方陣乘積可逆當且僅當各因子可逆 設 A , B A,B A , B 是實 p × p p\times p p × p 矩陣。A B AB A B 可逆,當且僅當 A , B A,B A , B 都可逆。此時
( A B ) − 1 = B − 1 A − 1 . (AB)^{-1}=B^{-1}A^{-1}. ( A B ) − 1 = B − 1 A − 1 .
證明
先構造逆矩陣,再使用它的名稱 先假設兩個因子都可逆,並令 C = B − 1 A − 1 C=B^{-1}A^{-1} C = B − 1 A − 1 。用結合律核對兩側:
( A B ) C = A ( B B − 1 ) A − 1 = I p , C ( A B ) = B − 1 ( A − 1 A ) B = I p . (AB)C=A(BB^{-1})A^{-1}=I_p,\qquad
C(AB)=B^{-1}(A^{-1}A)B=I_p. ( A B ) C = A ( B B − 1 ) A − 1 = I p , C ( A B ) = B − 1 ( A − 1 A ) B = I p . 現在才能斷言 C = ( A B ) − 1 C=(AB)^{-1} C = ( A B ) − 1 。反向次序使互逆因子彼此相鄰;論證沒有使用交換律。
證明反向時,假設 A B AB A B 可逆,記 C = ( A B ) − 1 C=(AB)^{-1} C = ( A B ) − 1 。此時還不能使用 A − 1 A^{-1} A − 1 或 B − 1 B^{-1} B − 1 ,應先寫出
A ( B C ) = I p , ( C A ) B = I p . A(BC)=I_p,\qquad (CA)B=I_p. A ( B C ) = I p , ( C A ) B = I p . 於是 A A A 有右逆,B B B 有左逆。兩個因子都是同階方陣,所以前面已證明的單側逆判別保證它們可逆。這一步引用等價字典,正是反向證明的關鍵。
由歸納法,同階方陣的有限乘積可逆,當且僅當每個因子都可逆,而且
( A 1 ⋯ A k ) − 1 = A k − 1 ⋯ A 1 − 1 . (A_1\cdots A_k)^{-1}=A_k^{-1}\cdots A_1^{-1}. ( A 1 ⋯ A k ) − 1 = A k − 1 ⋯ A 1 − 1 .
轉置與冪次
定理
轉置、純量倍乘與整數冪 若 A A A 可逆,則 ( A − 1 ) − 1 = A (A^{-1})^{-1}=A ( A − 1 ) − 1 = A ,且 ( A T ) − 1 = ( A − 1 ) T (A^T)^{-1}=(A^{-1})^T ( A T ) − 1 = ( A − 1 ) T 。對非零實數 c c c ,有 ( c A ) − 1 = c − 1 A − 1 (cA)^{-1}=c^{-1}A^{-1} ( c A ) − 1 = c − 1 A − 1 。
每個整數冪 A n A^n A n 都可逆,且 ( A n ) − 1 = ( A − 1 ) n (A^n)^{-1}=(A^{-1})^n ( A n ) − 1 = ( A − 1 ) n 。這裏 A 0 = I p A^0=I_p A 0 = I p ;對正整數 m m m ,定義 A − m = ( A − 1 ) m A^{-m}=(A^{-1})^m A − m = ( A − 1 ) m 。
轉置公式來自對兩條互逆等式取轉置,此時因子次序反轉。純量倍乘公式可在左右兩個乘積中約去 c c c 與 c − 1 c^{-1} c − 1 ;正整數冪來自乘積定理,負整數冪則對 A − 1 A^{-1} A − 1 使用該定理,零次冪就是單位矩陣。對正整數 n n n ,A n A^n A n 可逆也能反推 A A A 可逆;但 n = 0 n=0 n = 0 時反推不成立,因為奇異方陣仍滿足 A 0 = I p A^0=I_p A 0 = I p 。
乘積等於單位矩陣時怎樣推理
定理
方陣情況下,一側單位等式已經足夠 設 P P P 和 Q Q Q 是同大小方陣。若
P Q = I , PQ=I, P Q = I , 則 P P P 和 Q Q Q 都可逆,而且
Q = P − 1 , Q P = I . Q=P^{-1},
\qquad
QP=I. Q = P − 1 , QP = I .
例題
由 A B C D = I ABCD=I A B C D = I 得到循環恆等式 假設 A , B , C , D A,B,C,D A , B , C , D 都是 5 × 5 5 \times 5 5 × 5 矩陣,且
A B C D = I 5 . ABCD=I_5. A B C D = I 5 . 先把乘積分組成
( A B C ) D = I 5 . (ABC)D=I_5. ( A B C ) D = I 5 . 由一側單位等式的定理,A B C ABC A B C 和 D D D 互為逆矩陣。因此
D A B C = I 5 . DABC=I_5. D A B C = I 5 . 再用已經證明的循環等式繼續分組,可以得到
C D A B = I 5 , B C D A = I 5 . CDAB=I_5,
\qquad
BCDA=I_5. C D A B = I 5 , B C D A = I 5 . 所以 B C D A BCDA B C D A 、C D A B CDAB C D A B 和 D A B C DABC D A B C 確實必然等於 I 5 I_5 I 5 。但像 D C B A DCBA D C B A 或
D B A C DBAC D B A C 這些任意重排則沒有被迫成立;矩陣乘法仍然不可交換。
用行化簡計算逆矩陣
例題
用行化簡求逆矩陣 設
A = [ 1 2 3 5 ] . A = \begin{bmatrix}
1 & 2 \\
3 & 5
\end{bmatrix}. A = [ 1 3 2 5 ] . 先寫出 [ A ∣ I 2 ] [A\mid I_2] [ A ∣ I 2 ] :
[ 1 2 1 0 3 5 0 1 ] . \left[\begin{array}{cc|cc}
1 & 2 & 1 & 0 \\
3 & 5 & 0 & 1
\end{array}\right]. [ 1 3 2 5 1 0 0 1 ] . 先清走左下角的 3 3 3 :
[ 1 2 1 0 0 − 1 − 3 1 ] . \left[\begin{array}{cc|cc}
1 & 2 & 1 & 0 \\
0 & -1 & -3 & 1
\end{array}\right]. [ 1 0 2 − 1 1 − 3 0 1 ] . 再把第二行除以 − 1 -1 − 1 ,然後清走第一行第二欄:
[ 1 0 − 5 2 0 1 3 − 1 ] . \left[\begin{array}{cc|cc}
1 & 0 & -5 & 2 \\
0 & 1 & 3 & -1
\end{array}\right]. [ 1 0 0 1 − 5 3 2 − 1 ] . 所以
A − 1 = [ − 5 2 3 − 1 ] . A^{-1} = \begin{bmatrix}
-5 & 2 \\
3 & -1
\end{bmatrix}. A − 1 = [ − 5 3 2 − 1 ] .
一族逆矩陣
例題
由行化簡得到參數化逆矩陣 對實數 α \alpha α ,令
A α = [ 1 2 1 2 1 3 1 2 0 0 1 α 1 2 1 3 ] . A_\alpha =
\begin{bmatrix}
1 & 2 & 1 & 2 \\
1 & 3 & 1 & 2 \\
0 & 0 & 1 & \alpha \\
1 & 2 & 1 & 3
\end{bmatrix}. A α = 1 1 0 1 2 3 0 2 1 1 1 1 2 2 α 3 . 把增廣矩陣 [ A α ∣ I 4 ] [A_\alpha \mid I_4] [ A α ∣ I 4 ] 行化簡,得到
[ 1 0 0 0 5 − α − 2 − 1 α − 2 0 1 0 0 − 1 1 0 0 0 0 1 0 α 0 1 − α 0 0 0 1 − 1 0 0 1 ] . \left[\begin{array}{cccc|cccc}
1 & 0 & 0 & 0 & 5-\alpha & -2 & -1 & \alpha-2 \\
0 & 1 & 0 & 0 & -1 & 1 & 0 & 0 \\
0 & 0 & 1 & 0 & \alpha & 0 & 1 & -\alpha \\
0 & 0 & 0 & 1 & -1 & 0 & 0 & 1
\end{array}\right]. 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 5 − α − 1 α − 1 − 2 1 0 0 − 1 0 1 0 α − 2 0 − α 1 . 由於左邊對每個實數 α \alpha α 都是 I 4 I_4 I 4 ,矩陣 A α A_\alpha A α 對每個實數
α \alpha α 都可逆,而且
A α − 1 = [ 5 − α − 2 − 1 α − 2 − 1 1 0 0 α 0 1 − α − 1 0 0 1 ] . A_\alpha^{-1} =
\begin{bmatrix}
5-\alpha & -2 & -1 & \alpha-2 \\
-1 & 1 & 0 & 0 \\
\alpha & 0 & 1 & -\alpha \\
-1 & 0 & 0 & 1
\end{bmatrix}. A α − 1 = 5 − α − 1 α − 1 − 2 1 0 0 − 1 0 1 0 α − 2 0 − α 1 . 請留意這個結論的意思。參數仍然留在逆矩陣中;它不是需要解出的未知數。
行化簡真正說明的是:沒有任何 α \alpha α 會令左邊失去主元。
常見錯誤
不要把每個參數都當成限制條件 有時參數出現,是因為答案是一族逆矩陣,而不是因為矩陣在某些值下失效。
限制來自缺少主元或非法除以零,不是來自符號本身。
讀取已給出的行化簡表格
例題
一對互逆矩陣與轉置系統 假設一個行化簡表格顯示 [ A ∣ I 4 ] [A \mid I_4] [ A ∣ I 4 ] 與 [ I 4 ∣ D ] [I_4 \mid D] [ I 4 ∣ D ] 行等價,其中
D = [ 3 − 2 3 − 2 0 1 − 2 1 − 3 1 − 2 2 1 0 1 − 1 ] . D=
\begin{bmatrix}
3 & -2 & 3 & -2 \\
0 & 1 & -2 & 1 \\
-3 & 1 & -2 & 2 \\
1 & 0 & 1 & -1
\end{bmatrix}. D = 3 0 − 3 1 − 2 1 1 0 3 − 2 − 2 1 − 2 1 2 − 1 . 這個最終形狀本身已經說明
A − 1 = D . A^{-1}=D. A − 1 = D . 對應的係數矩陣為
A = [ 1 1 1 1 1 1 2 3 2 1 3 3 3 2 4 3 ] . A=
\begin{bmatrix}
1 & 1 & 1 & 1 \\
1 & 1 & 2 & 3 \\
2 & 1 & 3 & 3 \\
3 & 2 & 4 & 3
\end{bmatrix}. A = 1 1 2 3 1 1 1 2 1 2 3 4 1 3 3 3 . 直接相乘可得 A D = I 4 AD=I_4 A D = I 4 ,驗證了這對矩陣互為逆矩陣。現在令
g = [ 1 0 3 0 ] . g=
\begin{bmatrix}
1 \\ 0 \\ 3 \\ 0
\end{bmatrix}. g = 1 0 3 0 . 要求解 A t x = g A^t x=g A t x = g ,不需要從頭再行化簡一次。因為 A A A 可逆,所以 A t A^t A t 也可逆,
並且
( A t ) − 1 = ( A − 1 ) t = D t . (A^t)^{-1}=(A^{-1})^t=D^t. ( A t ) − 1 = ( A − 1 ) t = D t . 因此
x = ( A t ) − 1 g = D t g = [ − 6 1 − 3 4 ] . x=(A^t)^{-1}g=D^t g
=
\begin{bmatrix}
-6 \\ 1 \\ -3 \\ 4
\end{bmatrix}. x = ( A t ) − 1 g = D t g = − 6 1 − 3 4 . 這個計算很短,是因為行化簡表格已經完成了主要工作。
例題
用主元欄讀四個逆矩陣計算 以下每個矩陣都用同一個方法處理:先行化簡 [ A ∣ I 4 ] [A \mid I_4] [ A ∣ I 4 ] ,先看左邊區塊的
主元模式;只有當左邊區塊變成 I 4 I_4 I 4 時,右邊區塊才可讀作逆矩陣。
情況 (a)。 對
A 1 = [ 1 0 3 3 2 1 2 3 − 1 3 0 3 4 2 2 4 ] , A_1=
\begin{bmatrix}
1 & 0 & 3 & 3 \\
2 & 1 & 2 & 3 \\
-1 & 3 & 0 & 3 \\
4 & 2 & 2 & 4
\end{bmatrix}, A 1 = 1 2 − 1 4 0 1 3 2 3 2 0 2 3 3 3 4 , 其中一個行階梯階段的左邊區塊是
A 1 ♯ = [ 1 0 3 3 0 1 − 4 − 3 0 0 1 1 0 0 0 0 ] . A_1^\sharp=
\begin{bmatrix}
1 & 0 & 3 & 3 \\
0 & 1 & -4 & -3 \\
0 & 0 & 1 & 1 \\
0 & 0 & 0 & 0
\end{bmatrix}. A 1 ♯ = 1 0 0 0 0 1 0 0 3 − 4 1 0 3 − 3 1 0 . 第四欄沒有主元,所以 A 1 A_1 A 1 不可逆。
情況 (b)。 取前一個互逆矩陣例子中的 A 2 = A A_2=A A 2 = A 。行化簡得到 [ I 4 ∣ D ] [I_4 \mid D] [ I 4 ∣ D ] ,四欄都有主元,所以 A 2 − 1 = D A_2^{-1}=D A 2 − 1 = D 。前面用來解轉置系統的同一個矩陣,在這裏給出了一個可逆的情形。
情況 (c)。 對
A 3 = [ 1 2 − 1 30 3 2 0 15 0 1 3 15 1 1 1 10 ] , A_3=
\begin{bmatrix}
1 & 2 & -1 & 30 \\
3 & 2 & 0 & 15 \\
0 & 1 & 3 & 15 \\
1 & 1 & 1 & 10
\end{bmatrix}, A 3 = 1 3 0 1 2 2 1 1 − 1 0 3 1 30 15 15 10 , 其中一個行階梯階段的左邊區塊是
A 3 ♯ = [ 1 2 − 1 30 0 1 3 15 0 0 5 − 5 0 0 0 0 ] . A_3^\sharp=
\begin{bmatrix}
1 & 2 & -1 & 30 \\
0 & 1 & 3 & 15 \\
0 & 0 & 5 & -5 \\
0 & 0 & 0 & 0
\end{bmatrix}. A 3 ♯ = 1 0 0 0 2 1 0 0 − 1 3 5 0 30 15 − 5 0 . 同樣有主元欄缺失,所以 A 3 A_3 A 3 不可逆。
情況 (d)。 對
A 4 = [ 1 1 2 1 2 3 4 1 3 3 3 1 1 2 3 1 ] , A_4=
\begin{bmatrix}
1 & 1 & 2 & 1 \\
2 & 3 & 4 & 1 \\
3 & 3 & 3 & 1 \\
1 & 2 & 3 & 1
\end{bmatrix}, A 4 = 1 2 3 1 1 3 3 2 2 4 3 3 1 1 1 1 , 行化簡可達到 [ I 4 ∣ A 4 − 1 ] [I_4 \mid A_4^{-1}] [ I 4 ∣ A 4 − 1 ] ,其中
A 4 − 1 = [ 1 1 0 − 2 − 2 − 2 1 3 1 2 − 1 − 2 0 − 3 1 3 ] . A_4^{-1}=
\begin{bmatrix}
1 & 1 & 0 & -2 \\
-2 & -2 & 1 & 3 \\
1 & 2 & -1 & -2 \\
0 & -3 & 1 & 3
\end{bmatrix}. A 4 − 1 = 1 − 2 1 0 1 − 2 2 − 3 0 1 − 1 1 − 2 3 − 2 3 . 在情況 (b) 與 (d) 中,左邊的單位矩陣確認了可逆性。在情況 (a) 與 (c) 中,缺少主元就已排除可逆性,右邊區塊的數值並不能改變這個結論。
例題
最後的單位區塊控制整條行操作鏈 假設一條行操作鏈把 [ A ∣ I 5 ] [A \mid I_5] [ A ∣ I 5 ] 經過中間矩陣 [ B ∣ C ] [B \mid C] [ B ∣ C ] ,最後化成
[ D ∣ E ] [D \mid E] [ D ∣ E ] 。在最後幾步之後,所得增廣矩陣是
[ D ∣ E ] = [ 1 0 0 0 0 1 0 − 2 0 0 0 1 0 0 0 9 1 − 18 − 3 − 2 0 0 1 0 0 − 6 2 11 2 0 0 0 0 1 0 − 3 0 6 1 0 0 0 0 0 1 − 6 0 12 2 1 ] . [D \mid E]
=
\left[
\begin{array}{ccccc|ccccc}
1 & 0 & 0 & 0 & 0 & 1 & 0 & -2 & 0 & 0 \\
0 & 1 & 0 & 0 & 0 & 9 & 1 & -18 & -3 & -2 \\
0 & 0 & 1 & 0 & 0 & -6 & 2 & 11 & 2 & 0 \\
0 & 0 & 0 & 1 & 0 & -3 & 0 & 6 & 1 & 0 \\
0 & 0 & 0 & 0 & 1 & -6 & 0 & 12 & 2 & 1
\end{array}
\right]. [ D ∣ E ] = 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 1 9 − 6 − 3 − 6 0 1 2 0 0 − 2 − 18 11 6 12 0 − 3 2 1 2 0 − 2 0 0 1 . 左邊區塊是 I 5 I_5 I 5 ,所以 D = I 5 D=I_5 D = I 5 ,而
A − 1 = E = [ 1 0 − 2 0 0 9 1 − 18 − 3 − 2 − 6 2 11 2 0 − 3 0 6 1 0 − 6 0 12 2 1 ] . A^{-1}=E=
\begin{bmatrix}
1 & 0 & -2 & 0 & 0 \\
9 & 1 & -18 & -3 & -2 \\
-6 & 2 & 11 & 2 & 0 \\
-3 & 0 & 6 & 1 & 0 \\
-6 & 0 & 12 & 2 & 1
\end{bmatrix}. A − 1 = E = 1 9 − 6 − 3 − 6 0 1 2 0 0 − 2 − 18 11 6 12 0 − 3 2 1 2 0 − 2 0 0 1 . 於是可以立即推出 A t A^t A t 可逆,且
( A t ) − 1 = E t = [ 1 9 − 6 − 3 − 6 0 1 2 0 0 − 2 − 18 11 6 12 0 − 3 2 1 2 0 − 2 0 0 1 ] . (A^t)^{-1}=E^t
=
\begin{bmatrix}
1 & 9 & -6 & -3 & -6 \\
0 & 1 & 2 & 0 & 0 \\
-2 & -18 & 11 & 6 & 12 \\
0 & -3 & 2 & 1 & 2 \\
0 & -2 & 0 & 0 & 1
\end{bmatrix}. ( A t ) − 1 = E t = 1 0 − 2 0 0 9 1 − 18 − 3 − 2 − 6 2 11 2 0 − 3 0 6 1 0 − 6 0 12 2 1 . 同一條行操作鏈也可以被讀成一次左乘。令 H H H 為整條鏈中所有行操作矩陣的
乘積。因為 H A = I 5 HA=I_5 H A = I 5 且 H I 5 = E HI_5=E H I 5 = E ,所以 H = A − 1 H=A^{-1} H = A − 1 。因此,如果同一批行操作
把 F + I 5 F+I_5 F + I 5 化成 A 3 + 3 A + I 5 A^3+3A+I_5 A 3 + 3 A + I 5 ,就有
A 3 + 3 A + I 5 = A − 1 ( F + I 5 ) . A^3+3A+I_5=A^{-1}(F+I_5). A 3 + 3 A + I 5 = A − 1 ( F + I 5 ) . 左乘 A A A 得到
F + I 5 = A ( A 3 + 3 A + I 5 ) = A 4 + 3 A 2 + A , F+I_5=A(A^3+3A+I_5)=A^4+3A^2+A, F + I 5 = A ( A 3 + 3 A + I 5 ) = A 4 + 3 A 2 + A , 所以
F = A 4 + 3 A 2 + A − I 5 . F=A^4+3A^2+A-I_5. F = A 4 + 3 A 2 + A − I 5 . 同一條行操作鏈對每個矩陣的作用,都是左乘同一個矩陣 H = A − 1 H=A^{-1} H = A − 1 。因此,這次數值化簡也確定了關於 F F F 的方程。
不用行列式證明或否定可逆性
例題
由向量恆等式證明可逆 假設 A A A 是 7 × 7 7 \times 7 7 × 7 矩陣,而且
A 2 x = A x + x A^2x = Ax + x A 2 x = A x + x 對每個 x ∈ R 7 x \in \mathbb R^7 x ∈ R 7 都成立。把含有 A A A 的項移到左邊:
( A 2 − A ) x = x 對每個 x . (A^2-A)x = x
\qquad\text{對每個 }x. ( A 2 − A ) x = x 對每個 x . 因為 A 2 − A = ( A − I 7 ) A A^2-A=(A-I_7)A A 2 − A = ( A − I 7 ) A ,這就是
( ( A − I 7 ) A ) x = I 7 x 對每個 x . ((A-I_7)A)x = I_7x
\qquad\text{對每個 }x. (( A − I 7 ) A ) x = I 7 x 對每個 x . 若兩個矩陣作用在每個向量上的結果都相同,它們就是同一個矩陣,所以
( A − I 7 ) A = I 7 . (A-I_7)A = I_7. ( A − I 7 ) A = I 7 . 因此 A − I 7 A-I_7 A − I 7 是 A A A 的左逆。由於 A A A 是方陣,一側單位等式已經足夠推出:
A A A 可逆,而且其逆矩陣是 A − I 7 A-I_7 A − I 7 。
例題
用反證法證明某個多項式矩陣不可逆 假設 A A A 和 B B B 是兩個不同的 6 × 6 6 \times 6 6 × 6 矩陣,且
A 3 = B 3 , A 2 B = B 2 A . A^3=B^3,
\qquad
A^2B=B^2A. A 3 = B 3 , A 2 B = B 2 A . 我們證明 A 2 + B 2 A^2+B^2 A 2 + B 2 不可逆。令
C = A 2 + B 2 . C=A^2+B^2. C = A 2 + B 2 . 先比較 C B CB C B 與 C A CA C A :
C B = ( A 2 + B 2 ) B = A 2 B + B 3 , CB=(A^2+B^2)B=A^2B+B^3, C B = ( A 2 + B 2 ) B = A 2 B + B 3 , 而
C A = ( A 2 + B 2 ) A = A 3 + B 2 A . CA=(A^2+B^2)A=A^3+B^2A. C A = ( A 2 + B 2 ) A = A 3 + B 2 A . 利用已知的 A 2 B = B 2 A A^2B=B^2A A 2 B = B 2 A 和 B 3 = A 3 B^3=A^3 B 3 = A 3 ,可得
C B = C A . CB=CA. C B = C A . 如果 C C C 可逆,左乘 C − 1 C^{-1} C − 1 就會得到 B = A B=A B = A ,這與 A A A 、B B B 不同矛盾。
所以 C = A 2 + B 2 C=A^2+B^2 C = A 2 + B 2 不可逆。
常見錯誤
常見錯誤
不要把長方形矩陣的一側逆當成另一側逆 對非方陣而言,有左逆並不保證有右逆,有右逆也不保證有左逆;若兩側逆都存在,則必定相同。對於方陣,前面的等價性定理保證:任何一側逆的存在,都能推出唯一的雙側逆存在。
常見錯誤
不要靠外表估計可逆性 一個矩陣表面上看來很簡單,也未必可逆。真正方法是行化簡,或者用上面列出的等價條件。
常見錯誤
不要在未證明時任意交換因子 由 A B C D = I ABCD=I A B C D = I 可以透過逆矩陣理論推出某些循環恆等式,例如 D A B C = I DABC=I D A B C = I 。但這
不代表可以把因子任意改排成 D C B A DCBA D C B A ;除非已知矩陣彼此可交換,否則不能
這樣做。
快速檢查
思考檢查
如果 A A A 可逆,A − 1 A A^{-1}A A − 1 A 是甚麼?
解答 · 答案 A − 1 A = I p A^{-1}A = I_p A − 1 A = I p 。
思考檢查
如果 A A A 可逆,齊次方程組 A x = 0 Ax = 0 A x = 0 會不會有非零解?
解答 · 答案 不會。因為可逆時,A x = b Ax = b A x = b 對每個 b b b 都有唯一解;取 b = 0 b = 0 b = 0 時,唯一解就是 x = 0 x = 0 x = 0 。
思考檢查
如果 A A A 可逆,A t A^t A t 是否可逆?
解答 · 答案 是。A t A^t A t 可逆,而且 ( A t ) − 1 = ( A − 1 ) t (A^t)^{-1} = (A^{-1})^t ( A t ) − 1 = ( A − 1 ) t 。
思考檢查
若 B = G A B = GA B = G A ,而 G G G 可逆,同時 A A A 的欄滿足 a 3 = 2 a 1 − a 2 a_3 = 2a_1 - a_2 a 3 = 2 a 1 − a 2 ,那麼 B B B 的對應欄必須滿足甚麼關係? 保留同一組係數,並用 b j = G a j b_j = Ga_j b j = G a j 。
解答 · 答案 B B B 的對應欄必須滿足
b 3 = 2 b 1 − b 2 . b_3 = 2b_1 - b_2. b 3 = 2 b 1 − b 2 . 行等價會改變欄向量本身,但不會改變線性關係所用的係數模式。
思考檢查
假設 A , B , C , D A,B,C,D A , B , C , D 都是方陣,且 A B C D = I ABCD=I A B C D = I 。哪一個循環等式必然成立:D A B C = I DABC=I D A B C = I 還是 D C B A = I DCBA=I D C B A = I ? 把原式分組成 ( A B C ) D = I (ABC)D=I ( A B C ) D = I 。
解答 · 答案 D A B C = I DABC=I D A B C = I 必然成立。因為 ( A B C ) D = I (ABC)D=I ( A B C ) D = I ,方陣的一側單位等式定理說明 D D D 是
A B C ABC A B C 的逆矩陣,所以 D ( A B C ) = I D(ABC)=I D ( A B C ) = I 。而 D C B A = I DCBA=I D C B A = I 是把因子倒序重排,通常不會
被原式保證。
思考檢查
在上面的參數化逆矩陣中,當 α = 3 \alpha=3 α = 3 時,A α − 1 A_\alpha^{-1} A α − 1 的 ( 1 , 4 ) (1,4) ( 1 , 4 ) 元素是多少? 把 α = 3 \alpha=3 α = 3 代入第一行第四欄。
解答 · 答案 ( 1 , 4 ) (1,4) ( 1 , 4 ) 元素是 α − 2 \alpha-2 α − 2 ,所以當 α = 3 \alpha=3 α = 3 時,它等於 1 1 1 。
思考檢查
若對每個向量 x x x 都有 A 2 x = A x + x A^2x=Ax+x A 2 x = A x + x ,哪個矩陣可作為 A A A 的逆矩陣? 把所有項移到同一邊,再尋找含有 A A A 的因式分解。
解答 · 答案 逆矩陣是 A − I A-I A − I 。向量恆等式給出 ( A − I ) A = I (A-I)A=I ( A − I ) A = I ;對方陣而言,這個一側單位
等式足以推出 A A A 可逆。
思考檢查
若行操作把 [ A ∣ I p ] [A \mid I_p] [ A ∣ I p ] 化成 [ I p ∣ D ] [I_p \mid D] [ I p ∣ D ] ,D D D 是甚麼?
解答 · 答案 D = A − 1 D=A^{-1} D = A − 1 。同一批把 A A A 化成 I p I_p I p 的行操作,會把 I p I_p I p 化成 A A A 的
逆矩陣。
思考檢查
行變換把 [ A ∣ I p ] [A\mid I_p] [ A ∣ I p ] 化成 [ I p ∣ E ] [I_p\mid E] [ I p ∣ E ] ,其中 A A A 是方陣。說明為甚麼除了 E A = I p EA=I_p E A = I p ,還必有 A E = I p AE=I_p A E = I p 。 區分行操作乘積直接給出的等式,與方陣單側逆定理補上的結論。
解答 · 答案 設初等矩陣的乘積為 H H H ,兩個分塊直接給出 H A = I p HA=I_p H A = I p 和 H I p = E HI_p=E H I p = E 。因此 H = E H=E H = E ,並且 E A = I p EA=I_p E A = I p 。由於 A A A 與 E E E 是同階方陣,左逆也必是右逆,所以 E = A − 1 E=A^{-1} E = A − 1 且 A E = I p AE=I_p A E = I p 。最後這個推論必須使用方陣假設。
練習
思考檢查
假設 A A A 可逆,而且 A B = I p AB = I_p A B = I p 。證明 B = A − 1 B = A^{-1} B = A − 1 。 用逆矩陣唯一性或者左乘 A − 1 A^{-1} A − 1 都可以。
解答 · 引導解答 因為 A A A 可逆,所以存在 A − 1 A^{-1} A − 1 。將等式 A B = I p AB = I_p A B = I p 左乘 A − 1 A^{-1} A − 1 ,得
A − 1 ( A B ) = A − 1 I p . A^{-1}(AB) = A^{-1}I_p. A − 1 ( A B ) = A − 1 I p . 由結合律,
( A − 1 A ) B = A − 1 . (A^{-1}A)B = A^{-1}. ( A − 1 A ) B = A − 1 . 而 A − 1 A = I p A^{-1}A = I_p A − 1 A = I p ,所以
B = A − 1 . B = A^{-1}. B = A − 1 . 因此任何可逆方陣的右逆都必定等於其逆矩陣。
思考檢查
假設 ( A − B ) A = I (A-B)A=I ( A − B ) A = I ,其中 A A A 和 B B B 是方陣。證明 A A A 和 B B B 可交換。 先使用方陣的一側單位等式定理,再展開 A B − B A AB-BA A B − B A 。
解答 · 引導解答 由 ( A − B ) A = I (A-B)A=I ( A − B ) A = I ,方陣的一側單位等式定理給出
A ( A − B ) = I . A(A-B)=I. A ( A − B ) = I . 現在計算交換子:
A B − B A = A B − A 2 + A 2 − B A . AB-BA
= AB-A^2+A^2-BA. A B − B A = A B − A 2 + A 2 − B A . 把兩部分重新分組:
A B − A 2 = A ( B − A ) = − A ( A − B ) , AB-A^2 = A(B-A) = -A(A-B), A B − A 2 = A ( B − A ) = − A ( A − B ) , 以及
A 2 − B A = ( A − B ) A . A^2-BA = (A-B)A. A 2 − B A = ( A − B ) A . 因此
A B − B A = − A ( A − B ) + ( A − B ) A = − I + I = 0. AB-BA = -A(A-B)+(A-B)A = -I+I=0. A B − B A = − A ( A − B ) + ( A − B ) A = − I + I = 0. 所以 A B = B A AB=BA A B = B A ,即 A A A 和 B B B 可交換。
相關筆記
這一頁特別依賴
2.3 高斯消元與 RREF 、
3.2 矩陣乘法、單位矩陣與線性方程組
和
3.3 轉置、對稱矩陣與反對稱矩陣 。
零空間只有零向量與解的唯一性之間的關係,見
4.1 齊次方程組與零空間 。