把行變換記錄成一個矩陣
初等行變換最初是矩陣上的操作步驟:交換兩行、把一行乘以非零純量、把一行的倍數加到另一行。在增廣矩陣上,只有把操作施加到受影響的整行——包括增廣線右邊的常數項——才是在等價地改寫方程。此時操作可以逆轉,因而保留原方程組的解集。
但同一件事還有一個更有結構的讀法:每一個初等行變換,都等同於左乘某個特別的方陣。這個觀點很重要,因為它把一串初等行變換轉成普通的矩陣等式。之後討論行化簡、可逆性、秩、行列式和基底時,這正是程序性操作和代數結構之間的橋。
這個結構觀點亦回答一個實際問題:若同一個行變換要同時作用於許多欄,怎樣只記錄一次規則,而不逐項重做?單位矩陣恰好提供這種編碼。它的每一行會記錄新行由哪些舊行組成;隨後,矩陣乘法便把這套規則同時施加到相容矩陣的每一欄。
由單位矩陣建立初等矩陣
假設一個初等行變換 ρ 作用於有 p 行的矩陣。先從單位矩陣 Ip
開始,對 Ip 做同一個初等行變換,所得矩陣記作 Eρ。所有純量都屬於矩陣所採用的底域;特別地,每個非零倍數 c 在該域中都有倒數
1/c。
與 p 行矩陣上的初等行變換 ρ 對應的初等行變換矩陣,就是把
ρ 作用在 Ip 上所得的 p×p 矩陣。
這種矩陣亦簡稱初等矩陣;較長的名稱可以明確指出它與指定行變換的關係。
換句話說,
Eρ=ρ(Ip).
這個定義有用,是因為以下定理。
設 A 是任何有 p 行的矩陣,ρ 是一個初等行變換,並以
ρ(A) 表示把它作用於 A 後得到的矩陣。則
ρ(A)=EρA.所以,對 A 做這個初等行變換,等同於用同一個操作在 Ip 上得到的初等行變換矩陣左乘 A。
這裏必須是左乘。初等行變換改變的是行,會把新行寫成舊行的線性組合;左乘正是把
A 的新行寫成舊行的線性組合。右乘則會混合欄,不是同一件事。
維數亦說明了同一點。若 A 是 p×n 矩陣,則 Eρ 是
p×p 矩陣,所以 EρA 有定義,而且仍是 p×n。這裏完全不要求 A 是方陣。相反,對這 n 欄進行欄變換的矩陣大小為 n×n,並從右邊相乘。
為甚麼左乘會執行指定行變換
把任意 p×n 矩陣 A 的各行記作
a1,a2,…,ap,並以 ekT 表示 Ip 的第 k 行。基本等式
ekTA=ak表示左乘矩陣的一行會選取 A 的相應行。現在逐類檢查。
-
對 Ri↔Rj,Eρ 的第 i 行是
ejT,第 j 行是 eiT,其餘第 k 行仍是
ekT。因此 EρA 恰好交換 A 的第 i,j 行。
-
對 Ri←cRi(c=0),Eρ 的第 i 行是
ceiT,其餘行不變。因此乘積的第 i 行為
(ceiT)A=cai,
而每個 k=i 的第 k 行仍為 ak。
-
對 Rj←Rj+cRi(i=j),Eρ 的第 j 行是
ejT+ceiT,其餘行不變,所以新的第 j 行為
(ejT+ceiT)A=aj+cai,
正是指定的行加法。
三類情形窮盡了初等行變換,因此對每個相容矩陣 A 都有
EρA=ρ(A)。這個證明沒有聲稱任意改動行的方式都對應一個初等矩陣;
ρ 必須是剛才檢查的三類操作之一。
作用於增廣矩陣
假設 Ax=b 由增廣矩陣 [A∣b] 表示。初等行變換必須作用於整個增廣矩陣,不能只改其中一邊。相應的乘法等式是
Eρ[A∣b]=[EρA∣Eρb].
因此,係數和右端常數會接受完全相同的行組合。向量 x 滿足 Ax=b,若且唯若它滿足 EρAx=Eρb:正向可把原等式左乘 Eρ,反向則施加逆向初等行變換。這才是解集保持不變的嚴格理由。
若只改 A 而固定 b,或只改 b,一般會得到另一個方程組。增廣線只是視覺上的分隔,不是行變換可以中途停止的邊界。
分別記錄替換、縮放與交換
三種初等行變換,對應三類初等行變換矩陣。
對有三行的矩陣,考慮
ρ:R2←R2+3R1.把它作用在 I3 上:
I3=100010001⟼Eρ=130010001.因此,對任何有三行的矩陣 A,
EρA就是把 A 的第 2 行改成「第 2 行加上 3 倍第 1 行」後所得的矩陣。
對
ρ:R3←−2R3,對應的初等行變換矩陣是
Eρ=10001000−2.行倍乘要求純量非零,這一點在矩陣上也看得出來:如果倍數是 0,所得矩陣會出現零行,不能反向恢復。
對
ρ:R1↔R3,我們得到
Eρ=001010100.左乘這個矩陣,就會交換任何相容矩陣的第 1 行和第 3 行。
逐步理解初等矩陣的左乘
對單位矩陣施加行變換,便記錄了該變換的係數;再用所得矩陣左乘,就把同一個行變換施加到每一欄。
把行變換讀成左乘看清初等行變換如何變成左乘矩陣:先把操作作用在單位矩陣上,再用所得矩陣改變任何相容矩陣的行。
由 I 開始
對 ρ: R₂ ← R₂ + 3R₁,先把操作作用在 I₃ 上。
建立 Eρ
只有 I₃ 的第 2 行改變,得到第 2 行為 [3,1,0] 的 Eρ。
放在左邊相乘
同一個矩陣會對 A 執行該操作:ρ(A)=Eρ A。
讀新行
新的第 2 行等於舊第 2 行加三份舊第 1 行。
保留乘積次序
在 A₃=E₂E₁A₁ 中,E₁ 最靠近 A₁,因為 ρ₁ 最先作用。
反向得到逆矩陣
R₂ ← R₂ + 3R₁ 的反向操作是 R₂ ← R₂ − 3R₁,所以相應矩陣相乘得到 I。
行變換矩陣是把行變換先作用在單位矩陣上得到的。之後用它左乘任何相容矩陣,就會執行同一個行變換;反向行變換則給出逆矩陣。
比較行變換與矩陣乘積
用下面的步驟工具重溫一串行變換。每一步都要同時讀成:在左邊乘上對應的初等矩陣。
互動 · 逐步演示第 1 步,共 5 步
跟著走完一條行化簡路徑
每一步做一組基本行運算。留意主元上下的元素怎樣逐個變成 0。
先在第 1 列選主元。
1正在處理的主元0本步改動竪線右邊:常數
先由增廣矩陣開始。第一個主元的工作,是幫我們把它下面的元素清掉。
還未到行階梯形
第 2 行第 1 列仍是 1,沒有在第 1 行主元的右方。
要留意第 1 列的第一行已經有方便的主元 1,所以暫時不用換行。
一串初等行變換就是一個矩陣乘積
真正的用途不只是表示一個初等步驟,而是把整串初等行變換寫成一個矩陣乘積。
假設
A1ρ1A2ρ2A3ρ3⋯ρkAk+1.則
Ak+1=EρkEρk−1⋯Eρ2Eρ1A1.這裏每個 ρj 都必須是作用於同樣行數矩陣的初等行變換。這個乘積表示整串操作,但乘積本身未必仍是單一個初等行變換矩陣。
這個公式中的次序很重要。第一個行變換最靠近 A1,因為它最先作用:
A2=Eρ1A1,A3=Eρ2A2=Eρ2Eρ1A1.
令
A=101020110.依次做
ρ1: R2←R2+R1,ρ2: R1←R1+2R2.對應的行變換矩陣是
Eρ1=110010001,Eρ2=100210001.兩個操作完成後,結果就是
Eρ2Eρ1A=310210001101020110=311420520.若不想直接相乘 Eρ2Eρ1,也可以把同一串行變換按同一次序作用在
I3 上,所得矩陣就是合併後的左乘矩陣。最後的數值結果亦確認:第二步使用的是已經更新後的第 2 行。
讀懂較長的行變換乘積
把一長串初等行變換寫成矩陣乘積時,需要先分清三種對象:
- 正在被變換的矩陣 A1,A2,…,Ak+1;
- 每一步對應的初等行變換矩陣 H1,H2,…,Hk;
- 合併後的單一左乘矩陣 J=Hk⋯H2H1。
以下是一個四行矩陣的典型例子。假設
ρ1ρ2ρ3ρ4ρ5ρ6:R2←−21R2,:R1↔R2,:R3←R3−2R1,:R4←R4−2R3,:R1←R1−R3,:R1←R1−R2.
如果
A1ρ1A2ρ2A3ρ3A4ρ4A5ρ5A6ρ6A7,
那麼
A7=H6H5H4H3H2H1A1.
相應的行變換矩陣是
H1H3H5=10000−210000100001,=10−20010000100001,=10000100−10100001,H2H4H6=0100100000100001,=10000100001−20001,=1000−110000100001.
它們的乘積是
J=H6H5H4H3H2H1=−1100−2301−2−101−20001.
有效率的做法,是把六個行變換依次作用在 I4 上,而不是手動展開六個因子。矩陣 J 記錄了整串步驟對行的總作用:
A7=JA1.
這個等式也可用來檢查乘積次序。如果把第一步放在最左邊,所得乘積會描述另一串不同的操作。
快速檢查
(在上面六步行變換中,若 K 是由 A7 反向回到 A1 的行變換矩陣乘積,J 和 K 應滿足甚麼等式?)
解答 · 答案
矩陣 J 和 K 編碼的是每一個四行矩陣上的變換,不只是題中展示的
A1。對任意四行矩陣 X,先做正向步驟,再做反向步驟,會得到
K(JX)=X.取 X=I4,便有 KJ=I4。反過來,先做反向步驟,再做正向步驟,亦會恢復任意四行矩陣 Y,所以 J(KY)=Y;取 Y=I4,便有 JK=I4。因此
KJ=I4JK=I4.換句話說,K=J−1。
反向行變換與逆矩陣
每個初等行變換都有反向操作:
- Rj←Rj+cRi 的反向操作是
Rj←Rj−cRi;
- Ri←cRi(其中 c=0)的反向操作是
Ri←(1/c)Ri;
- 交換兩行的反向操作仍然是同一次交換。
這給出一個精確的矩陣敘述。
每個初等行變換矩陣都是可逆的。它的逆矩陣,就是反向初等行變換所對應的初等行變換矩陣。
設 Eρ 表示初等行變換 ρ,而 Eρ−1 表示它的反向初等操作。對每個有 p 行的矩陣 X,先做 ρ 再撤銷,都會恢復 X。由左乘定理,
Eρ−1EρX=X.這個等式對每個這樣的 X 都成立;特別取 X=Ip,得到
Eρ−1Eρ=Ip。先做反向操作再做 ρ,同理得到
EρEρ−1=Ip。因此
Eρ−1=Eρ−1.具體來說,交換兩行的逆仍是同一次交換;乘以 c=0 的逆是乘以 1/c;把第 i 行的 c 倍加到第 j 行,則由加上 −c 倍第 i 行來撤銷。
例如,若
E=130010001
表示 R2←R2+3R1,則
E−1=1−30010001
表示 R2←R2−3R1。
這就是初等行變換可逆的代數原因,也解釋了為甚麼行化簡和可逆矩陣有密切關係。
設 E1,…,Ek 按時間順序編碼同一組 p 行矩陣上的操作,合併左乘矩陣就是 J=Ek⋯E1。每個 Ej−1 表示對應初等行變換的反向操作。要恢復輸入,首先必須撤銷 Ek,因為它是最後一次對行所作的改變;然後依次撤銷 Ek−1,…,E1。因此,反向鏈的合併矩陣為
K=E1−1E2−1⋯Ek−1.這裏必須區分兩種次序:反向操作的時間順序從 Ek−1 開始;寫成乘積時,最先作用的因子卻放在最靠近待恢復矩陣的位置。以三步為例,由結合律可得
KJ=E1−1E2−1(E3−1E3)E2E1=E1−1(E2−1E2)E1=Ip.另一個乘積 JK 則從相鄰的 E1E1−1 開始消去,最後同樣得到 Ip。任意有限鏈都可這樣逐對消去。整個論證沒有交換因子,也不要求被變換的矩陣本身是方陣。
把這個讀法用於前面的兩步例子:合併矩陣是 Eρ2Eρ1,所以逆矩陣是 Eρ1−1Eρ2−1。從最終矩陣開始,先從第一行減去目前第二行的兩倍,第一行便恢復為 [1 0 1];再從第二行減去剛恢復的第一行,第二行便恢復為 [0 2 1]。撤銷操作是在逐個恢復中間狀態,乘積次序也就隨之確定,無須另背一個口訣。
為甚麼這個觀點之後會有用?
如果 B 與 A 行等價,則存在有限多個初等行變換把 A 變成 B。因此存在某個初等行變換矩陣的乘積 E,使得
B=EA.
由於每個初等行變換矩陣都可逆,乘積 E 也可逆。因此行等價既可以用初等操作步驟來描述,也可以用左邊乘上一個可逆矩陣的等式來描述。反向敘述則要小心:任意可逆左乘矩陣未必本身只代表一步初等操作。
這個觀點會在多個地方使用:
- 一個方陣若行等價於 In,它就是若干初等行變換矩陣的乘積;
- 初等行變換保留齊次方程組的解資訊,因為它等同於左乘可逆矩陣;
- 行列式關於初等行變換的規則,可以透過初等矩陣來表述;
- 秩和基底的論證可以使用行化簡,同時清楚知道原矩陣的欄本身通常已被改變。
常見錯誤
初等行變換由左乘表示。右乘會混合欄,而不是混合行。
如果 ρ1 先於 ρ2 作用,合併矩陣是
Eρ2Eρ1,不是 Eρ1Eρ2。
行倍乘要求純量非零。把一行乘以 0 不能反向恢復,也不會產生可逆的初等行變換矩陣。
從某個特定矩陣 A 滿足 KJA=A,一般不能推出 KJ=I,因為 A 的各欄未必能檢測所有可能輸入。應證明複合行變換會恢復每個相容矩陣,或直接把它們作用於 Ip。所以上面對 KJ=JK=Ip 的論證先使用任意 p 行矩陣,再代入單位矩陣。
按作用次序理解乘積
要編碼有 p 行矩陣上的一個初等行變換,只需把該操作施加到 Ip。所得初等矩陣 Eρ 會以左乘方式作用於每個相容矩陣;對交換、非零倍乘和行加法三類情形逐行檢查,就得到 EρA=ρ(A)。若對象是增廣矩陣,操作必須橫跨整行,包括右端常數,方程組才會保持完全相同的解集。
依次施加 ρ1,…,ρk 時,乘積按作用次序寫成
Eρk⋯Eρ1。每個因子都可逆,因為反向初等操作給出它的雙側逆矩陣。因此,有限的行化簡既可以看成逐步操作,亦可以看成一次可逆的左乘;各個因子則繼續記錄每一步的準確內容。
快速檢查
快速檢查
(對有三行的矩陣,E=105010001 表示哪個行變換?)
解答 · 答案
只有第 3 行被改變,而且它變成「第 3 行加上 5 倍第 1 行」。因此
R3←R3+5R1.
快速檢查
(R2←R2−4R1 的反向行變換是甚麼?)
解答 · 答案
反向操作是
R2←R2+4R1.
練習
快速檢查
(寫出有三行矩陣上 R1↔R2 的行變換矩陣。)
解答 · 引導解答
交換 I3 的前兩行,得到
010100001.這就是左乘時交換第 1 行和第 2 行的矩陣。
快速檢查
(假設 B 由 A 經以下兩步得到:先做 R2←R2+R1,再做 R3←2R3。把 B 寫成涉及 A 的矩陣乘積。)
解答 · 引導解答
令 E1 是 R2←R2+R1 的行變換矩陣,令 E2 是
R3←2R3 的行變換矩陣。因為 E1 先作用,而 E2 後作用,
B=E2E1A.
快速檢查
(令 β1,β2=0。假設五行矩陣 A 經以下行變換鏈得到 B:先做 R3←R3+α1R1,再做 R2←β1R2,再做 R1↔R4,再做 R3←R3+α2R2,最後做 R1←β2R1。寫出單一矩陣 G,使得 B=GA。)
把同一串操作作用在 I5 上,並記住後面的操作使用的是當時的行,不是原來的行。
解答 · 引導解答
把這些操作作用在 I5 上,得到
G=00α1100β1α2β10000100β2000000001.重點是 α2β1 這一項。操作
R3←R3+α2R2 發生在第 2 行已經被縮放成
β1R2 之後,所以來自原來第 2 行的貢獻是
α2β1R2。
快速檢查
(沿用上一題的矩陣 G,假設 D=GC。寫出矩陣 H,使得 C=HD。)
解答 · 引導解答
反向鏈是
GR1←β21R1R3←R3−α2R2R1↔R4R2←β11R2R3←R3−α1R1I5.因此
H=000β2100β11−α2000010010−α10000001.這就是 G−1,所以 D=GC 會推出 C=HD。
相關筆記
這一頁建立在
2.2 增廣矩陣與行變換
和
3.2 矩陣乘法、單位矩陣與線性方程組
之上。接著可讀
3.6 分塊矩陣,完成第三章。它也準備了
5.1 可逆矩陣
和
7.2 行變換、乘積與可逆性
中使用的行化簡代數觀點。
補充練習與完整解答
題 1
對(a)–(d),以下依次的行變換把A送到B。寫出從B返回A的操作鏈,並給出滿足B=HA的方陣H。大小分別為(a)、(b):3×4;(c):3×6;(d):4×5。
(a)
1. R2←R2+(−1)R12. R3←R3+(−2)R13. R3←R3+(−2)R24. R3←(−1)R35. R1←R1+(−2)R26. R2←R2+(−1)R3
(b)
1. R1↔R22. R1←(−1)R13. R3←R3+(−2)R14. R3←R3+(−3)R25. R1←R1+(−2)R2
(c)
1. R1↔R22. R3←R3+(2)R13. R3←R3+(−3)R24. R1←R1+(−2)R25. R1←R1+(2)R36. R2←R2+(−2)R3
(d)
1. R2←R2+(−1)R12. R3←R3+(−3)R13. R4←R4+(−1)R14. R2←(−1)R25. R3←R3+(4)R26. R4←R4+(3)R27. R3↔R48. R4←R4+(4)R39. R1←R1+(−2)R210. R1←R1+(−1)R3
完整解答
以下在實數上討論。行替換的來源行與目標行不同,縮放係數非零;操作作用於整行,包括所有增廣塊。有限鏈中第一步先作用,所以合併左乘矩陣中的初等矩陣按相反次序排列。
把正向操作逐一作用於Ir。以下Hj記錄每個中間乘數;對步數歸納,目前矩陣就是HjA,最後一個乘數即所求H。倒序撤銷各步:加行倍數改為減同一倍數,縮放改為乘倒數,交換仍作同一交換。每個逆操作對任意相容矩陣都撤銷原操作,因此所寫逆鏈確實返回A。
(a)
H0=100010001
H1=1−10010001
H2=1−1−2010001
H3=1−1001−2001
H4=1−1001200−1
H5=3−10−21200−1
H6=3−10−2−1201−1
逆鏈依以下次序:
1. R2←R2+(1)R32. R1←R1+(2)R23. R3←(−1)R34. R3←R3+(2)R25. R3←R3+(2)R16. R2←R2+(1)R1
(b)
H0=100010001
H1=010100001
H2=010−100001
H3=010−102001
H4=01−3−102001
H5=−21−3−102001
逆鏈依以下次序:
1. R1←R1+(2)R22. R3←R3+(3)R23. R3←R3+(2)R14. R1←(−1)R15. R1↔R2
(c)
H0=100010001
H1=010100001
H2=010102001
H3=01−3102001
H4=−21−3102001
H5=−81−3502201
H6=−87−35−422−21
逆鏈依以下次序:
1. R2←R2+(2)R32. R1←R1+(−2)R33. R1←R1+(2)R24. R3←R3+(3)R25. R3←R3+(−2)R16. R1↔R2
(d)
H0=1000010000100001
H1=1−100010000100001
H2=1−1−30010000100001
H3=1−1−3−1010000100001
H4=11−3−10−10000100001
H5=111−10−1−4000100001
H6=11120−1−4−300100001
H7=11210−1−3−400010010
H8=11290−1−3−1600010014
H9=−11292−1−3−1600010014
H10=−31295−1−3−160001−1014
逆鏈依以下次序:
1. R1←R1+(1)R32. R1←R1+(2)R23. R4←R4+(−4)R34. R3↔R45. R4←R4+(−3)R26. R3←R3+(−4)R27. R2←(−1)R28. R4←R4+(1)R19. R3←R3+(3)R110. R2←R2+(1)R1
題 3
設A,A′為m×p矩陣,B為p×q矩陣,且A′=ρ(A)。(a) 原題要求直接由定義證明AB=ρ(A′B),其中(i)Rk←Rk+αRi、i=k;(ii)Rk←βRk、β=0;(iii)Ri↔Rk。(b) 用歸納法證明:若A,A′行等價,則AB,A′B在同一有限操作序列下行等價。
完整解答
以下在實數上討論。行替換的來源行與目標行不同,縮放係數非零;操作作用於整行,包括所有增廣塊。有限鏈中第一步先作用,所以合併左乘矩陣中的初等矩陣按相反次序排列。
(a):(i)、(ii)原印方向不成立
對(i),取A=(1,0)T、B=[1],並作R2←R2+R1。則A′B=(1,1)T,再作同一操作得到(1,2)T=AB。對(ii),取A=B=[1]、β=2,則A′B=[2],而ρ(A′B)=[4]=[1]。所以原印的普遍命題不能被證明。正確方向是A′B=ρ(AB),等價地AB=ρ−1(A′B)。交換操作滿足ρ−1=ρ,所以(iii)原印方向成立。
直接逐項證明正向恆等式:在行替換情形,只有輸出第k行改變:
(A′B)kj=r=1∑p(akr+αair)brj=(AB)kj+α(AB)ij.
其餘輸出行逐項相同。縮放時(A′B)kj=β(AB)kj,其餘行不變。交換A的第i,k行,會恰好交換AB的第i,k行。由此直接證明三種正向傳遞恆等式;再作逆操作,便得到反向恆等式。
(b):對行變換序列長度歸納
長度為零時,兩矩陣相同,乘積也相同。假設結論對長度n成立。對長度n+1的序列,令An為前n步後的矩陣。歸納假設表明,同樣的前n步把AB送到AnB;再由(a)已逐項證明的正向恆等式,最後一步把AnB送到A′B。所以整個n+1步序列把AB送到A′B。歸納得任意有限長度的結論。所有乘積大小均為m×q,因此始終是同尺寸矩陣上的合法行變換。
題 4
3×5矩陣依次滿足A1R3+2R1A2R2↔R3A33R1A4。(a)–(c) 給出滿足Aj+1=JjAj的J1,J2,J3;(d) 給出A4=JA1中的J;(e) 給出A1=J′A4中的J′。
完整解答
以下在實數上討論。行替換的來源行與目標行不同,縮放係數非零;操作作用於整行,包括所有增廣塊。有限鏈中第一步先作用,所以合併左乘矩陣中的初等矩陣按相反次序排列。
每個Jj由對應的一步操作作用於I3得到。連續代入可得A4=J3J2J1A1。具體為:
J1=102010001
J2=100001010
J3=300010001
J=J3J2J1=320001010
J′=J−1=310−32001010
逆過程先把第1行乘1/3,再交換第2、3行,最後從第3行減去第1行的兩倍。把這些操作作用於所寫J得到I3;因此J′J=JJ′=I3,對任意A1都得到要求的反向等式。
題 5
下面六個4×6矩陣是連續階段A,A(1),…,A(5)=A′。(a) 指出五步行變換ρ1,…,ρ5;(b) 給出滿足A′=HA的明確4×4矩陣H。
A0stage=210−1131−2−101348403−4−308−4−4−1
A1stage=120−1311−20−1138440−43−30−48−4−1
A2stage=100−13−51−20−1138−1240−411−30−416−4−1
A3stage=10003−5110−1138−1248−411−3−4−416−4−5
A4stage=100031−5101−1384−128−4−311−4−4−416−5
A5stage=1000311101513845128−4−3−511−4−4−4−516−5
完整解答
以下在實數上討論。行替換的來源行與目標行不同,縮放係數非零;操作作用於整行,包括所有增廣塊。有限鏈中第一步先作用,所以合併左乘矩陣中的初等矩陣按相反次序排列。
(a) 比較整行,包括最後兩欄,依次得到:
1. R1↔R22. R2←R2+(−2)R13. R4←R4+(1)R14. R2↔R35. R3←(−51)R3
例如,交換後第2行減第1行的兩倍,得到(0,−5,−1,−12,11,16);第4行加第1行,得到(0,1,3,8,−4,−5)。第四步交換第2、3行;第五步把當時的第3行乘−1/5,最後兩項也同時縮放。(b) 把同樣操作作用於I4,逐步得到:
H0=1000010000100001
H1=0100100000100001
H2=01001−20000100001
H3=01001−20100100001
H4=001010−2101000001
H5=00−5101052101000001
取H=H5。每個記錄行恰為相應原行的線性組合,所以逐項證明了A′=HA。
題 6
考慮下面一串5×7矩陣,並將中間階段記為U0,U1,U2,U3。箭頭為Aρ1U0→⋯→U1ρ2U2→⋯→U3ρ3B1,其後是所示五步B1→⋯→B6。(a) 求ρ1,ρ2,ρ3;(b) 填寫B3(I)(II)(III)B1的三處空白;(c) 求滿足B5=KB4的K;(d) 是否存在從I5到某個滿足B6=HB1的H的行變換鏈?若有,須寫出全部操作和每個中間矩陣,包括H;若無,說明不存在。
A=311−14212−13113−12−20−1−2−2132−94−2−1−22−3−311−2−2
U0=131−14122−13113−120−2−1−2−2312−94−1−2−22−31−31−2−2
U1=100041−11031−22020−2−1−2−23−8−1−64−11−11−31−60−1−2
U2=100001−110−11−220−20−2−1−2−23−8−1−6−8−11−1111−60−1−6
U3=10000110001200002−3−2038−9−60−1−101016−6−10
B1=100001100012000021−20383−60−1−1010162−10
B2=1000011000120000210038300−1−101016230
B3=1000001000−12000−22100−583000−1010−56230
B4=1000001000−1200002100183000−1010−16230
B5=1000001000−1200000100123000−1010−12230
B6=1000001000−12000001001230000010−15230
從B1開始的五個箭頭依次為:
1. R4←R4+(2)R32. R1←R1+(−1)R23. R1←R1+(2)R34. R2←R2+(−2)R35. R2←R2+(1)R4
完整解答
以下在實數上討論。行替換的來源行與目標行不同,縮放係數非零;操作作用於整行,包括所有增廣塊。有限鏈中第一步先作用,所以合併左乘矩陣中的初等矩陣按相反次序排列。
(a) 比較三個標記箭頭兩側的整行:
ρ1:R1↔R2,ρ2:R5←R5−4R1,ρ3:R3←−31R3.
對ρ2,從(4,3,2,−2,4,−3,−2)減去(1,1,1,0,3,−1,1)的四倍,得到(0,−1,−2,−2,−8,1,−6)。對ρ3,縮放(0,0,0,−3,−9,0,−6)後得到(0,0,0,1,3,0,2)。(b) 把B1→B2→B3的兩步倒序撤銷:
(I):R1←R1+R2,(II)=B2,(III):R4←R4−2R3.
(c) 只有第2行改為第2行減第3行的兩倍,故:
K=10000010000−21000001000001
(d) 存在。對I5作以下同一組五步操作:
1. R4←R4+(2)R32. R1←R1+(−1)R23. R1←R1+(2)R34. R2←R2+(−2)R35. R2←R2+(1)R4
H0=1000001000001000001000001
H1=1000001000001200001000001
H2=10000−11000001200001000001
H3=10000−11000201200001000001
H4=10000−110002−21200001000001
H5=10000−11000201200101000001
取H=H5。上面已經寫出I5=H0→H1→⋯→H5的每個矩陣。每一步都以同一方式改變原輸入行的組合,因此把操作鏈作用於B1時,第j步有HjB1=Bj+1;特別地HB1=B6。這對矩形B1同樣成立,完全不需要把B1求逆。
題 7
3×4矩陣A1,…,A5按下面四步操作相連。(a) 求滿足Aj+1=JjAj的J1,…,J4;(b) 若增廣矩陣[A1∣B1],…,[A4∣B4],[A5∣I3]也按完全相同的四步相連,各Bj均為3×3,求全部四個Bj;(c) B1A1是否與A5行等價?若是,寫出操作鏈的每一步;若否,解釋原因。
1. R2←R2+(−1)R12. R3←R3+(−2)R13. R3←R3+(−1)R24. R1←R1+(−1)R2
完整解答
以下在實數上討論。行替換的來源行與目標行不同,縮放係數非零;操作作用於整行,包括所有增廣塊。有限鏈中第一步先作用,所以合併左乘矩陣中的初等矩陣按相反次序排列。
(a) 把各步指定操作作用於I3,得到:
J1=1−10010001
J2=10−2010001
J3=10001−1001
J4=100−110001
(b) 每個行變換也作用於右塊,因此J1B1=B2、J2B2=B3、J3B3=B4、J4B4=I3。從最後一個等式開始,依次施行逆操作:
B4=J4−1,B3=J3−1J4−1,B2=J2−1J3−1J4−1,B1=J1−1J2−1J3−1J4−1.
B1=112123001
B2=102113001
B3=100111001
B4=100110001
把這些明確矩陣代入上述四個右塊等式,也可逐一核對;並沒有對矩形Aj作非法消去。(c) 是。令P=J4J3J2J1。增廣操作鏈同時給出A5=PA1及PB1=I3。所以對B1A1施行一次四步操作,得到P(B1A1)=A1;再作一次四步,得到PA1=A5。因此完整操作鏈依次為:
1. R2←R2+(−1)R12. R3←R3+(−2)R13. R3←R3+(−1)R24. R1←R1+(−1)R25. R2←R2+(−1)R16. R3←R3+(−2)R17. R3←R3+(−1)R28. R1←R1+(−1)R2
題 8
設α1,α2,α3,β1,β2,β3為任意實數。四行矩陣A,B經下面六步操作相連。給出滿足B=GA及A=HB的明確4×4矩陣G,H。
1. R2←R2+(α1)R12. R3←R3+(α2)R13. R4←R4+(α3)R14. R1←R1+(β1)R25. R1←R1+(β2)R36. R1←R1+(β3)R4
完整解答
以下在實數上討論。行替換的來源行與目標行不同,縮放係數非零;操作作用於整行,包括所有增廣塊。有限鏈中第一步先作用,所以合併左乘矩陣中的初等矩陣按相反次序排列。
把原行記為a1,a2,a3,a4。前三步後各行為a1、a2+α1a1、a3+α2a1、a4+α3a1。最後三步使用這些已經更新的行,因此最終首行為(1+α1β1+α2β2+α3β3)a1+β1a2+β2a3+β3a4,其餘三行保持上述形式。這四條原行組合給出G;把六步倒序撤銷,則得到H:
G=α1β1+α2β2+α3β3+1α1α2α3β1100β2010β3001
H=1−α1−α2−α3−β1α1β1+1α2β1α3β1−β2α1β2α2β2+1α3β2−β3α1β3α2β3α3β3+1
具體地,依次從第1行減去β3R4、β2R3、β1R2,恢復a1;再從第4、3、2行分別減去α3R1、α2R1、α1R1,恢復全部原行。若最終各行為bj,這就是說a1=b1−β1b2−β2b3−β3b4,且ai+1=bi+1−αia1;展開便恰好得到所寫H。所以HG=GH=I4。這些操作全是行替換,沒有任何參數必須非零,即使G左上角等於零也不影響結論。
題 9
設實數α1,α2,β1,β2滿足β1β2=0。四行矩陣A,B經下面五步操作相連。(a) 給出一個逐項明確的4×4矩陣G,使B=GA,並寫出從I4到G的操作鏈;不必寫中間矩陣。(b) 寫出從G返回I4的操作鏈。
1. R3←R3+(α1)R22. R1←(β1)R13. R1↔R44. R1←R1+(α2)R35. R2←(β2)R2
完整解答
以下在實數上討論。行替換的來源行與目標行不同,縮放係數非零;操作作用於整行,包括所有增廣塊。有限鏈中第一步先作用,所以合併左乘矩陣中的初等矩陣按相反次序排列。
(a) 設原行為a1,a2,a3,a4。第一步使第3行成為a3+α1a2;縮放第1行再交換第1、4行,使第1行為a4、第4行為β1a1。把目前第3行的α2倍加到第1行,得到a4+α2a3+α1α2a2。最後第2行變為β2a2。因此:
G=000β1α1α2β2α10α20101000
從I4出發的所求操作鏈,就是依原次序施行題目五步;它們作用於標準行後,逐項得到所寫G。(b) 倒序撤銷:
1. R2←(β21)R22. R1←R1+(−α2)R33. R1↔R44. R1←(β11)R15. R3←R3+(−α1)R2
前兩步逆操作後,第2行為e2、第1行為e4。交換後第4行恢復為e4,第1行成為β1e1;乘倒數恢復e1,最後一步把第3行從e3+α1e2恢復為e3。終點因此為I4。β1,β2非零的條件,恰好用於這兩次倒數縮放。
題 10
設m×n矩陣A,A′行等價。證明對每個n×p矩陣B,AB與A′B行等價。
完整解答
以下在實數上討論。行替換的來源行與目標行不同,縮放係數非零;操作作用於整行,包括所有增廣塊。有限鏈中第一步先作用,所以合併左乘矩陣中的初等矩陣按相反次序排列。
證明中把內維數和輸出欄數分別記為p和q;這只是重新命名題目的維數符號。
(a):單步正向傳遞
直接逐項證明正向恆等式:在行替換情形,只有輸出第k行改變:
(A′B)kj=r=1∑p(akr+αair)brj=(AB)kj+α(AB)ij.
其餘輸出行逐項相同。縮放時(A′B)kj=β(AB)kj,其餘行不變。交換A的第i,k行,會恰好交換AB的第i,k行。由此直接證明三種正向傳遞恆等式;再作逆操作,便得到反向恆等式。
(b):對行變換序列長度歸納
長度為零時,兩矩陣相同,乘積也相同。假設結論對長度n成立。對長度n+1的序列,令An為前n步後的矩陣。歸納假設表明,同樣的前n步把AB送到AnB;再由(a)已逐項證明的正向恆等式,最後一步把AnB送到A′B。所以整個n+1步序列把AB送到A′B。歸納得任意有限長度的結論。所有乘積大小均為m×q,因此始終是同尺寸矩陣上的合法行變換。