10. 消元法与三角分解

数值分析



周维祺

引入:消元法

例子

  • 孙子算经:今有雉兔同笼,上有三十五头,下有九十四足,问雉兔各几何?



消元法的策略

利用初等(行)变换将矩阵化为上三角形式,再自下而上依次解出未知量

对角阵和单位阵

  • 对角阵:主对角线以外的元素全为
    为对角阵,若时,,例:

  • 单位阵:主对角线元素全为的对角阵,记作

三角阵

  • 上三角阵:主对角线以下的元素全为
    为上三角阵,若时,,例:

  • 下三角阵:主对角线以上的元素全为
    为上三角阵,若时,,例:

练习

证明下三角阵的乘积仍是下三角阵

初等变换1

  • 交换矩阵中的两行/两列

  • 为交换单位阵中第行和第行所得的矩阵

  • 交换中第行和第行所得的矩阵即

  • 交换中第列和第列所得的矩阵即 A

初等变换1: 交换矩阵中的两行/两列

初等变换2

  • 将第行(列)乘以常数

  • 为将中第行乘以常数所得的矩阵

  • 中第行乘以常数所得的矩阵即

  • 中第列乘以常数所得的矩阵即

初等变换2: 将第行/列乘以常数

初等变换3

  • 将第行(列)中每个元素都加上第行(列)中对应元素的

  • 为将中第行中每个元素都加上第行中对应元素的倍所得的矩阵

  • 加上倍所得矩阵即

  • 加上倍所得的矩阵即

初等变换3: 将第行/列中每个元素都加上第行/列中对应元素的

练习

写出上述三种初等变换的逆矩阵

消元法的分解

  • 的逆也均为三角阵,因此其乘积也仍是三角阵
    特别地,是下三角阵

  • 不断运用初等行变换使之成为上三角阵

三角分解(分解):

其中可逆,是下三角阵,是上三角阵

练习

作以下矩阵的三角分解:

运用初等行变换(1)的场合

例:

练习

  • 分别是对角阵、下三角阵、上三角阵,证明

  • 都是下三角阵

  • 都是上三角阵

  • 若D可逆,则也是对角阵

三角分解的唯一性

  • 可逆,的一个三角分解

  • 是任意可逆对角阵,也是的一个三角分解

  • 因此,一般要求三角分解中的主对角线是,此时该分解唯一

练习

  • 是主对角线元素都为的下三角阵,证明也是主对角线元素都是的下三角阵

  • 是上三角阵,证明也是上三角阵

唯一性的证明

  • 可逆,的两个三角分解,且的主对角线元素都为



  • 是对角阵,直接计算得其主对角线元素都为

直接作三角分解

  • 例:

  • 第一行:
    代入第一列:
    代入第二行,再代入第二列....

小结

  • 消元法的概念

  • 三种初等行变换所对应的矩阵及其逆矩阵

  • 三角分解的两种计算方法

  • 三角分解具有唯一性时的形式