用直接三角分解法解线性方程组.ppt

§ 5 用直接三角分解法解线性方程组 5.1 矩阵的三角分解 列主元高斯消去法实质上是对方程组进行等价变形，即是对 定理 10 1 、矩阵的杜里特尔（ Doolittle ）分解 系数矩阵施行 行 初等变换，这些 初等变换 又可以用矩阵表示。因此 矩阵的 三角分解 是列主元高斯消去法的另一种表示方法，或着说 高斯消去法的变形，即是高斯消去法 紧凑格式 的矩阵表示（矩阵的 n n R A ? ? 设 ， ) 1 , , 2 , 1 ( 0 ? ? ? n k D k ? ，如果 A 顺序主子式 其中 L 为 单位下三角 阵， U 为 上三角 矩阵。 说明： 则 A 可 唯一 分解为两个三角矩阵相乘，即 。 LU A ? L 为单位下三角阵 , 则 A 可 唯一 分解为两个三角矩阵相乘 , 否则， 分解不是唯一的。 三角分解），它在解方程组的直接法中起着重要的作用。 ： 存在性 因为 A 的顺序主子式非零，则由定理 8 ，使 阵 1 2 1 , , , ? n L L L ? U A L L L n ? ? 1 2 1 ? , , , 2 , 1 , 0 ) ( n k a k kk ? ? ? 得 ， 其中 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1 1 1 1 , 1 nk k k k m m L ? ? ? 证明： ， 对于 b x A ? ? ? 再由高斯消去法定理 6 知，存在初等下三角 ) ( 上三角阵 k k ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1 1 1 1 , , 1 1 k n k k k m m L ? ? ? U L L L A n ) ( 1 1 1 2 1 1 ? ? ? ? ? ? ? k k ， ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1 1 1 1 2 1 32 31 21 ? ? ? ? ? ? n n m m m m m ， LU ? ? ? 1 1 1 2 1 1 ? ? ? ? ? n L L L L ? 其中 n n ij a A ? ? ) ( ) 1 ( ， ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ) ( ) 2 ( 2 ) 2 ( 22 ) 1 ( 1 ) 1 ( 12 ) 1 ( 11 n nn n n a a a a a a U ? ? ? ? ， ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ) ( ) 2 ( 2 ) 2 ( 22 ) 1 ( 1 ) 1 ( 12 ) 1 ( 11 n nn n n a a a a a a U ? ? ? ? ， ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1 1 1 1 2 1 32 31 21 ? ? ? ? ? ? n n m m m m m ? ? 1 1 1 2 1 1 , , , ? ? ? ? ? n L L L L ? 其中 U L L L A n ) ( 1 1 1 2 1 1 ? ? ? ? ? ? ? ， LU ? 也就是 A 分解为一个单位下三角阵 L 与一个上三角阵 U 的乘积，即 。 LU A ? 为 其中 L L , 1 唯一性： . 1 1 LU U L A A ? ? 有两个三角分解，即 设 1 1 1 1 ? ? ? UU L L 由于上三角阵的逆仍为上三角阵，单位下三角阵的逆仍为单位下三 ) , ( 1 1 1 1 1 ? ? ? U L LU U L 右乘 两边同时左乘 角阵；上三角阵的乘积是上三角阵，单位下三角阵的乘积为单位下 三角阵。等式左边是单位下三角阵，右边是上三角阵，要使等式 成立，只能等于单位矩阵 I 。于是 , 1 1 1 1 I UU L L ? ? ? ? ，即此三角分解唯一。 则 1 1 , U U L L ? ? 上三角阵，于是有 非奇异 为 单位下三角阵， ) ( , 1 U U ) ( 为非奇异阵 仅 A ? ? ? ? ? ? ? ? ? ? ? ? ? 3 1 1 1 7 3 1 2 1 A 例 7 设 ，试将 A 进行三角分解。 解： 由高斯消去法得到 ? ? ? ? ? ? ? ? ? ? ? ? ? 1 1 1 31 21 1 m m L ? ? ? ? ? ? ? ? ? ? ? ? 1 1 1 32 2 m L ， ? ? ? ? ? ? ? ? ? ? ? ? ? 1 1 1 3 1 ? ? ? ? ? ? ? ? ? ? ? 1 1 1 1 1 1 1 3 1 3 31 21 ? ? ? ? m m ， 1 1 1 32 ? ? ? ? m ? ? ? ? ?