4.7 列满秩方程
4.7 列滿秩方程
對于可逆矩陣,高斯約當(dāng)消元法可以求得方程的解和逆矩陣,同樣也可以應(yīng)用于列滿秩矩陣。列滿秩矩陣進(jìn)行高斯約當(dāng)消元法,最終矩陣變?yōu)閱挝痪仃嚭土憔仃嚒@绶匠?/p>
2x+4y=24x+9y=86x+13y=102x + 4y = 2 \\ 4x + 9y = 8 \\ 6x + 13y = 10 2x+4y=24x+9y=86x+13y=10
系數(shù)矩陣為
A=[2449613]A= \left[ \begin{matrix} 2 & 4 \\ 4 & 9 \\ 6 & 13 \end{matrix} \right] A=???246?4913????
是列滿秩矩陣。
方程有 222 個未知數(shù),但是有 333 個方程,一般來說是矛盾方程,無解。但當(dāng) b\mathbf{b}b 能被矩陣 AAA 向量組表示時,有唯一解。
增廣矩陣進(jìn)行高斯消元法
[24249861310]?[242014014]?[242014000]\left[ \begin{matrix} 2 & 4 & 2\\ 4 & 9 & 8\\ 6 & 13 & 10 \end{matrix} \right] \Rightarrow \left[ \begin{matrix} 2 & 4 & 2\\ 0 & 1 & 4\\ 0 & 1 & 4 \end{matrix} \right]\Rightarrow \left[ \begin{matrix} 2 & 4 & 2\\ 0 & 1 & 4\\ 0 & 0 & 0 \end{matrix} \right] ???246?4913?2810????????200?411?244????????200?410?240????
注意此時最后一個方程變?yōu)?0x+0y=00x+0y=00x+0y=0 ,是永遠(yuǎn)成立的平凡方程!真正有效的方程數(shù)量是 222 個。方程解為 x=?7,y=4x=-7,y=4x=?7,y=4 。
如果第3個方程系數(shù)改變
2x+4y=24x+9y=86x+13y=112x + 4y = 2 \\ 4x + 9y = 8 \\ 6x + 13y = 11 2x+4y=24x+9y=86x+13y=11
最后變換為
[242014001]\left[ \begin{matrix} 2 & 4 & 2\\ 0 & 1 & 4\\ 0 & 0 & 1 \end{matrix} \right] ???200?410?241????
注意此時最后一個方程變?yōu)?0x+0y=10x+0y=10x+0y=1 ,無解!
總結(jié)如下,列滿秩矩陣 AmnA_{mn}Amn?,高斯消元法變換為 [UnnOm?n,n]\left[ \begin{matrix} U_{nn} \\ \mathbf{O}_{m-n,n} \end{matrix} \right][Unn?Om?n,n??] ,UnnU_{nn}Unn? 是 nnn 階上三角陣,其對角元素是矩陣 AAA 的主元且均不為零,Om?n,n\mathbf{O}_{m-n,n}Om?n,n? 是零矩陣。矩陣乘法表示為 LmmA=[UnnOm?n,n]L_{mm}A = \left[ \begin{matrix} U_{nn} \\ \mathbf{O}_{m-n,n} \end{matrix} \right]Lmm?A=[Unn?Om?n,n??] ,LmmL_{mm}Lmm? 是 mmm 階單位下三角陣。
對向量 b\mathbf{b}b ,如果 Lmmb=[b′0]L_{mm}\mathbf{b}=\left[ \begin{matrix} \mathbf{b'} \\ \mathbf{0} \end{matrix} \right]Lmm?b=[b′0?] ,即后 m?nm-nm?n 個分量都為 000 ,則方程 Amnx=bA_{mn}\mathbf{x}=\mathbf{b}Amn?x=b 有唯一解;只要后 m?nm-nm?n 個分量有一個不為 000 ,則方程無解。
UnnU_{nn}Unn? 繼續(xù)進(jìn)行高斯約當(dāng)消元,則變?yōu)閱挝痪仃?EnnE_{nn}Enn? ,用矩陣乘法表示為 PmmA=[EnnOm?n,n]P_{mm}A = \left[ \begin{matrix} E_{nn} \\ \mathbf{O}_{m-n,n} \end{matrix} \right]Pmm?A=[Enn?Om?n,n??] ,PmmP_{mm}Pmm? 是 mmm 階可逆矩陣。
上面的方法是沒有進(jìn)行行對調(diào),如果需要則可以先進(jìn)行行對調(diào),再進(jìn)行消元。
矩陣乘法表示,即對任意列滿秩矩陣 AAA ,存在可逆矩陣 P,QP,QP,Q ,使 PAQ=[UnnOm?n,n]PAQ=\left[ \begin{matrix} U_{nn} \\ \mathbf{O}_{m-n,n} \end{matrix} \right]PAQ=[Unn?Om?n,n??] 成立,進(jìn)一步對 UnnU_{nn}Unn? 進(jìn)行高斯約當(dāng)消元,則可表示為,存在可逆矩陣 P,QP,QP,Q ,使 PAQ=[EnnOm?n,n]PAQ=\left[ \begin{matrix} E_{nn} \\ \mathbf{O}_{m-n,n} \end{matrix} \right]PAQ=[Enn?Om?n,n??] 成立。
前面章節(jié)介紹了,列滿秩矩陣的行向量組是相關(guān)組,其極大無關(guān)組是 nnn 維空間的基。利用高斯消元法可以找到列滿秩矩陣的極大無關(guān)組,變換后的矩陣 UnnU_{nn}Unn? 對應(yīng)到矩陣 AAA 的行向量即是極大無關(guān)組。本例中
[2449]\left[ \begin{matrix} 2 & 4 \\ 4 & 9 \\ \end{matrix} \right][24?49?] 這兩個行向量就是矩陣 AAA 行向量的極大無關(guān)組。
總結(jié)
- 上一篇: 4.6 高斯约当消元法
- 下一篇: 4.8 行满秩方程