第五章 线性方程组直接解法.pptVIP

第五章 线性方程组直接解法 矩阵三角分解 LU 分解 LU 分解的存在唯一性 类似分解 LU分解紧凑方式 LU分解紧凑算法（续） LU分解紧凑算法（续） Crout分解的紧凑算法 SPD 的 Cholesky 分解 Cholesky 分解实现算法 LU分解求解线性方程组 举例（一） 追赶法 追赶法（续） 作业 * * 第三节 矩阵三角分解 将一个矩阵分解成结构简单的三角形矩阵的乘积称为矩阵的三角分解。 高斯消去过程其实就是一个矩阵的三角分解过程。 则 A(k) 与 A(k+1) 之间的关系式可以表示为： 其中： ( i = k+1, …, n ) , 将 Gauss 消去过程中第 k-1 步消元后的系数矩阵记为： ( k = 1, …, n-1) 记： ，则 其中：L --- 单位下三角矩阵，U --- 上三角矩阵 LU 分解 于是有： 容易验证： ( k = 1, …, n-1) (杜利脱尔Doolittle分解) LU 分解存在 普通高斯消去法不被中断 ? 定理 普通高斯消去法求解方程组 Ax = b 时，主元 的充要条件是：A 的所有顺序主子式不为零。 定理 若 A的所有顺序主子式 ，则 A存在唯一的 LU分解。 ( L