麻省理工線性代數筆記(二)

第二講主要內容:

  • 高斯消元法(Elimination)來解線性方程組(success or failure)
  • Back substitution(後向回代)
  • Elimination Matrices(消元矩陣)

1. 高斯消元法求解下面方程組

第一步:假設x,y,z滿足第一個方程,然後用該方程乘以某個數去消除第二個及第三個方程中的x;

本例中是第一個方程乘以-3加到第二個方程中,因為第三個方程本來就沒有x項,不變,因此最終方程組變成:

第二步:第二個方程不動,第二個方程乘以-2加到第三個方程中,得到以下方程組:

第三步,Back substitution(後向回代)

肯定是先算第三個方程,得到z;然後把z的值回代到第二個方程,得到y值;最後把y,z值回代到第一個方程,從而得到x值。至此我們已用高斯消元法得到了方程組的解。

2. 用矩陣來表示高斯消元法過程

把方程組右邊的向量放到矩陣豎線後面,稱為增廣矩陣。

反向回代簡單,計算機很快能做完,我們目前主要研究消元之後的矩陣。

其實大家可以看到,高斯消元法一直在做初等行變換!

首先,主元(pivot)不能為0,有一種情況就是0佔據了主元的位置,此時需要做行交換,找到主元不為0的行交換過來(還是在做初等行變換)。

其次,若已經找不到不為0的主元了,也就是不能通過行變換來得到不為0的主元了,比如高斯消元後得到下面的矩陣:

僅有一個主元1,其餘主元的位置均為0,此時高斯消元法失效,方程組可能有解,可能無解(此例無解),需要考察方程組右邊的向量情況(後面會涉及)。

3. Elimination Matrices(消元矩陣)

回憶上面消元的過程,第一個方程乘以-3加到第二個方程中,寫成相當於矩陣左乘一個消元矩陣:

提示:左乘是進行行變換,右乘是列變換!

我們把此消元矩陣記為E21,因為消掉的是2行1列上面的元素。以此類推,需要在原矩陣左乘一系列的消元矩陣,得到最終的上三角形式,這個例子是再左乘一個E32。

將這些消元矩陣乘在一起:

消元矩陣是下三角,乘積也是下三角,即得到右側完美等式

但是還是有一點小瑕疵,剛才提到如果主元位置上有0,可能需要行交換,行交換矩陣稱為置換矩陣(Permutation)。

若想交換矩陣的第一行和第二行,左乘置換矩陣(通過單位矩陣交換相應行即可):

也就是說上面提到的完美等式,少數情況需要左乘一些置換矩陣,但對於計算機來講並不是什麼大問題。

4. 逆矩陣

把消元矩陣

拿過來,它做的事情是把第二行加上第一行乘以-3,反向操作是第二行加上第一行乘以3,相當於沒加沒減。

這裡引出了逆矩陣的定義,先感性認識一下。

本節課內容結束。


推薦閱讀:

數學分析筆記整理(一)
定積分
經典求極限例題
傅里葉級數(三角)
物理化學公式閑談(1)

TAG:線性代數 | 微積分 | 高等數學 |