4.8 行满秩方程
4.8 行滿秩方程
高斯約當消元法還可以應用于行滿秩矩陣,行滿秩矩陣進行高斯約當消元法,最終矩陣變為單位矩陣和自由矩陣。例如方程
2x+4y+6z=124x+9y+13z=362x + 4y + 6z = 12 \\ 4x + 9y + 13z= 36 2x+4y+6z=124x+9y+13z=36
系數矩陣為
A=[2464913]A= \left[ \begin{matrix} 2 & 4 & 6\\ 4 & 9 & 13 \end{matrix} \right] A=[24?49?613?]
是行滿秩矩陣。
方程有 333 個未知數,但是只有 222 個方程,有無窮多解。
增廣矩陣進行高斯消元法
[24612491336]?[2461201112]\left[ \begin{matrix} 2 & 4 & 6 & 12\\ 4 & 9 & 13 & 36 \end{matrix} \right] \Rightarrow \left[ \begin{matrix} 2 & 4 & 6 & 12\\ 0 & 1 & 1 & 12 \end{matrix} \right] [24?49?613?1236?]?[20?41?61?1212?]
方程變為
2x+4y+6z=12y+z=122x + 4y + 6z = 12 \\ y + z = 12 2x+4y+6z=12y+z=12
未知數 zzz 是自由變量,可取任意實數,方程解可以表示為
x=?18?zy=12?zz=0+zx = -18 - z \\ y = 12 - z \\ z = 0 + z x=?18?zy=12?zz=0+z
寫成向量形式
[xyz]=[?18120]+z[?1?11]\left[ \begin{matrix} x \\ y \\ z \end{matrix} \right]= \left[ \begin{matrix} -18 \\ 12 \\ 0 \end{matrix} \right]+ z\left[ \begin{matrix} -1 \\ -1 \\ 1 \end{matrix} \right] ???xyz????=????18120????+z????1?11????
zzz 取任意實數。
這是只有一個自由變量的情況,如果有多個呢?舉例說明。
2x+4y+6z+w=124x+9y+13z+4w=362x + 4y + 6z + w = 12 \\ 4x + 9y + 13z + 4w = 36 2x+4y+6z+w=124x+9y+13z+4w=36
系數矩陣為
A=[246149134]A= \left[ \begin{matrix} 2 & 4 & 6 & 1\\ 4 & 9 & 13 & 4 \end{matrix} \right] A=[24?49?613?14?]
是行滿秩矩陣。
方程有 444 個未知數,但是只有 222 個方程,有無窮多解。
增廣矩陣進行高斯消元法
[2461124913436]?[246112011212]\left[ \begin{matrix} 2 & 4 & 6 & 1 & 12\\ 4 & 9 & 13 & 4 & 36 \end{matrix} \right] \Rightarrow \left[ \begin{matrix} 2 & 4 & 6 & 1 & 12\\ 0 & 1 & 1 & 2 & 12 \end{matrix} \right] [24?49?613?14?1236?]?[20?41?61?12?1212?]
方程變為
2x+4y+6z+w=12y+z+2w=122x + 4y + 6z + w = 12 \\ y + z + 2w= 12 2x+4y+6z+w=12y+z+2w=12
未知數 z,wz,wz,w 是自由變量,可取任意實數,方程解可以表示為
x=?18?z?wy=12?z?2wz=0+z+0ww=0+0z+wx = -18 - z - w \\ y = 12 - z - 2w \\ z = 0 + z + 0w \\ w = 0 + 0z + w x=?18?z?wy=12?z?2wz=0+z+0ww=0+0z+w
寫成向量形式
[xyzw]=[?181200]+z[?1?110]+w[?1?201]\left[ \begin{matrix} x \\ y \\ z \\ w \end{matrix} \right]= \left[ \begin{matrix} -18 \\ 12 \\ 0 \\ 0 \end{matrix} \right]+ z\left[ \begin{matrix} -1 \\ -1 \\ 1 \\ 0 \end{matrix} \right]+ w\left[ \begin{matrix} -1 \\ -2 \\ 0 \\ 1 \end{matrix} \right] ?????xyzw??????=??????181200??????+z??????1?110??????+w??????1?201??????
z,wz,wz,w 取任意實數。
有幾點需要說明,第 111 點,自由變量不一定非取 z,wz,wz,w ,有多種取法,比如可取 y,wy,wy,w 。
方程為
2x+4y+6z+w=12y+z+2w=122x + 4y + 6z + w = 12 \\ y + z + 2w= 12 2x+4y+6z+w=12y+z+2w=12
未知數 y,wy,wy,w 是自由變量,可取任意實數,方程解可以表示為
x=?60+2y+11wz=12?y?2wy=0+y+0ww=0+0y+wx = -60 + 2y + 11w \\ z = 12 - y - 2w \\ y = 0 + y + 0w \\ w = 0 + 0y + w x=?60+2y+11wz=12?y?2wy=0+y+0ww=0+0y+w
寫成向量形式
[xzyw]=[?601200]+y[2?110]+w[11?201]\left[ \begin{matrix} x \\ z \\ y \\ w \end{matrix} \right]= \left[ \begin{matrix} -60 \\ 12 \\ 0 \\ 0 \end{matrix} \right]+ y\left[ \begin{matrix} 2 \\ -1 \\ 1 \\ 0 \end{matrix} \right]+ w\left[ \begin{matrix} 11 \\ -2 \\ 0 \\ 1 \end{matrix} \right] ?????xzyw??????=??????601200??????+y?????2?110??????+w?????11?201??????
y,wy,wy,w 取任意實數。
這兩組解的形式雖然差別很大,但其實都是方程的解。
第 222 點,自由變量雖然有多種取法,但自由變量的數量是固定的,為 n?mn-mn?m ,是矩陣 AAA 零空間的維度。
第 333 點,向量 [2?110]\left[ \begin{matrix} 2 \\ -1 \\ 1 \\ 0 \end{matrix} \right]?????2?110?????? 和向量 [11?201]\left[ \begin{matrix} 11 \\ -2 \\ 0 \\ 1 \end{matrix} \right]?????11?201?????? 是方程 Ax=0A\mathbf{x} = \mathbf{0}Ax=0 解,稱為零解,它們線性組合構成矩陣 AAA 零空間。一般來說,每個零解中自由變量組均取單位向量 ei\mathbf{e}_iei? ,即只有一個自由變量取 111 ,其它均為 000 ,這樣計算比較方便。
第 444 點,向量 [?601200]\left[ \begin{matrix} -60 \\ 12 \\ 0 \\ 0 \end{matrix} \right]??????601200?????? 是方程 Ax=bA\mathbf{x} = \mathbfAx=b 解,稱為特解。一般來說,特解中自由變量均取 000 ,這樣計算比較方便,當然也可取任意值。
所以行滿秩方程解的結構為:特解 +++ 零解的線性組合。
第 555 點,非自由變量對應的列向量,是矩陣列空間的基。本例中,[2449]\left[ \begin{matrix} 2 & 4 \\ 4 & 9 \end{matrix} \right][24?49?] 這兩個列向量是列空間的基。由于自由變量有多種取法,故基也有多種。
前面都是非自由變量對應的列向量在矩陣的最前列,即最前面的上三角陣是可逆的,但有時不可逆。例如方程
2x+y+3z+2w=54x+2y+7z+5w=122x + y + 3z + 2w = 5 \\ 4x + 2y + 7z + 5w = 12 2x+y+3z+2w=54x+2y+7z+5w=12
增廣矩陣進行高斯消元法
[21325427512]?[2132500112]\left[ \begin{matrix} 2 & 1 & 3 & 2 & 5\\ 4 & 2 & 7 & 5 & 12 \end{matrix} \right] \Rightarrow \left[ \begin{matrix} 2 & 1 & 3 & 2 & 5\\ 0 & 0 & 1 & 1 & 2 \end{matrix} \right] [24?12?37?25?512?]?[20?10?31?21?52?]
前面兩列是矩陣
[2100]\left[ \begin{matrix} 2 & 1 \\ 0 & 0 \end{matrix} \right] [20?10?]
對角元素有 000 ,上三角陣不可逆。對應方程為
2x+y+3z+2w=50x+0y+z+w=22x + y + 3z + 2w = 5 \\ 0x + 0y + z + w = 2 2x+y+3z+2w=50x+0y+z+w=2
那怎么選非自由變量呢?我們把非自由變量稱為主元。每行選一個,從左到右選首個非零元素對應的變量為主元。第一行首個非零元素是 222 ,對于主元是 xxx ,第二行首個非零元素是 111 ,對于主元是 zzz ;則自由變量就是 y,wy,wy,w ,可取任意實數,方程解可以表示為
x=?1?y+wz=2?0y?wy=0+y+0ww=0+0y+wx = -1 - y + w \\ z = 2 - 0y - w \\ y = 0 + y + 0w \\ w = 0 + 0y + w x=?1?y+wz=2?0y?wy=0+y+0ww=0+0y+w
寫成向量形式
[xzyw]=[?1200]+y[?1010]+w[1?101]\left[ \begin{matrix} x \\ z \\ y \\ w \end{matrix} \right]= \left[ \begin{matrix} -1 \\ 2 \\ 0 \\ 0 \end{matrix} \right]+ y\left[ \begin{matrix} -1 \\ 0 \\ 1 \\ 0 \end{matrix} \right]+ w\left[ \begin{matrix} 1 \\ -1 \\ 0 \\ 1 \end{matrix} \right] ?????xzyw??????=??????1200??????+y??????1010??????+w?????1?101??????
y,wy,wy,w 取任意實數。
矩陣 AAA 經過高斯消元變換為 [21320011]\left[ \begin{matrix} 2 & 1 & 3 & 2 \\ 0 & 0 & 1 & 1 \end{matrix} \right][20?10?31?21?] ,但我們希望前兩列向量是上三角陣,這樣更有利于理論分析和對比。這可通過列對調操作實現,即對調 2,32,32,3 列向量,則前兩列向量是上三角陣。列對調操作類似行對調操作,可以通過置換矩陣實現,差別在于:行對調是左乘,列對調是右乘,即 PijAP_{ij}APij?A 對調矩陣 i,ji,ji,j 兩行, APijAP_{ij}APij? 對調矩陣 i,ji,ji,j 兩列。
總結如下,行滿秩矩陣 AmnA_{mn}Amn?,高斯消元法變換為 [Umm,Fm,n?m]\left[ \begin{matrix} U_{mm} , F_{m,n-m} \end{matrix} \right][Umm?,Fm,n?m??] ,UmmU_{mm}Umm? 是 mmm 階上三角陣,其對角元素是矩陣 AAA 的主元,Fm,n?mF_{m,n-m}Fm,n?m? 是自由矩陣。用矩陣乘法表示為 LmmA=[Umm,Fm,n?m]L_{mm}A = \left[ \begin{matrix} U_{mm} , F_{m,n-m} \end{matrix} \right]Lmm?A=[Umm?,Fm,n?m??] ,LmmL_{mm}Lmm? 是 mmm 階單位下三角陣。上三角陣對應的變量是主元,自由矩陣對應的變量是自由變量。方程解結構為:特解加零解。如有必要,需進行列對調。矩陣乘法表示,即對任意行滿秩矩陣 AAA ,存在可逆矩陣 P,QP,QP,Q ,使 PAQ=[Umm,Fm,n?m]PAQ=\left[ \begin{matrix} U_{mm} , F_{m,n-m} \end{matrix} \right]PAQ=[Umm?,Fm,n?m??] 成立,進一步對 UmmU_{mm}Umm? 進行高斯約當消元,則可表示為,存在可逆矩陣 P,QP,QP,Q ,使 PAQ=[Emm,Fm,n?m]PAQ=\left[ \begin{matrix} E_{mm} , F_{m,n-m} \end{matrix} \right]PAQ=[Emm?,Fm,n?m??] 成立。
前面章節介紹了,行滿秩矩陣的列向量組是相關組,其極大無關組是 mmm 維空間的基。利用高斯消元法可以找到行滿秩矩陣的極大無關組,主元對應的列向量即是極大無關組。
總結
- 上一篇: 4.7 列满秩方程
- 下一篇: 4.9 行列均不满秩方程