第3章 解线性方程组及直接方法2.ppt

* * 2．2 三角分解法 　　2.2.1 杜里特尔分解法 　　求解线性代数议程组的三角分解法，起源于高斯消去法的矩阵形式。 　　高斯消去法消去过程中，将变换后增广矩阵的第k行-c倍加于第i行，相当于左乘初等矩陈 它们都是单位下三角矩阵，即对角元全为1、对角线上方元素全为零的矩阵。因此不选主元的高斯消去法消去过程，实质是增广矩陈 被左乘一系列倍加矩阵，变成上三角形矩阵 ，即 此式称为高斯消去法的矩阵形式。由此显然 注意 是将单位矩阵 的第 行倍数加于第 行， ，将第一行的倍数加于第 行、 、 第二行，可见 是单位下三角矩阵。故 这说明，高斯消去法的消去过程，实质上是把系数矩阵 分解为 单位下三角矩阵 与上三角矩阵 的乘积，并且求解议程组 的过程。回代过程就是求解上三角形方程组 矩阵 和 也可直接算出。事实上，比较等式 两 边等 行、第 列元素可知 注意 是单位下三角矩阵， 便知 从而 同样，因 为上三角阵， 知 可见 公式（2-2）和（2-3）就是计算 和 各元素的计算公式。 实际计算时 的对角元 不必存放， 和 中 肯定为零的元素也不必存放，因此 的 可共同存放在增广 矩阵 的位置： 此时公式（2-2）、（2-3）表明， 或 都是原始矩阵 对 应元素，减去同行左边 的元素与同列上边 的元素乘积；只是对 的元素，然后需除以 的对角元。计算顺序，通常先算 的第 行，再算 的第 列；也可先算 的第 列，再算 的第 行， 如图2—1所示： 图2—1 计算顺序 例2—1 分解 ，并解方程组 ，其中 解 按计算公式（2-2）和（2-3） 详细计算过程如下（下文不再写出）： 从而 回代（解方程组 ），得 分解 且 为单位下三角阵、 为上三角阵，称为杜里 特尔Dolittlse)分解。利用杜里特尔分解求解方程组 或 ， 相当于解两个三角形方程组 解下三角方程组 可以在分解 时同时完成（如例2—1）， 也可独立完成。这是因为，把 写成分量形式，就是 由此可见， 用杜里特尔分解求解方程组（2-1），所需乘除次数与高斯消 去法完全一样。其中分解 需 次，解 需 次，解 需 次，共计 次。 它们都是单位下三角矩阵，即对角全为1、对角线上方元素全为零的矩阵。因此不选主元的高斯消去过程，实质是增广矩阵 被左乘一系列倍加矩阵，变成上三角形矩阵 ,即 此式称为高斯消去法的矩阵形式。由此显然 注意 是将单位矩阵 三角分解法常用于求解系数矩阵都是 的若干方程式组 这是因为，一旦完成分解 ，只需再解 个三角形方程组 解这种三角形方程组每组只需 次乘除法，远比重复使用高斯消 去法节省工作量。 为保证三角分解顺序、稳定进行，与高斯消去法一样，也可选 主元。常用列主元法。 ２.２.２ 克洛特分解法 当矩阵 可作杜里特尔分解 时,令 为 对角元构 成的对角阵 则 再算第 行；或者先算第 行，再算第 列， 如图 2—2所示。克洛