解线性方程组的矩阵三角分解法31.pptVIP

* 第五章解线性方程组的直接方法 计算方法 —— 矩阵三角分解法 * 本讲内容 一般线性方程组 LU 分解与 PLU 分解 对称正定线性方程组 平方根法－－Cholesky 分解 对角占优三对角线性方程组 追赶法 * LU 分解 将一个矩阵分解成结构简单的三角形矩阵的乘积 矩阵的三角分解 矩阵的 LU(Doolittle) 分解 矩阵的 LDR 分解 克洛脱 (Crout) 分解 * 计算 LU 分解 利用矩阵乘法直接计算 LU 分解 L ? U = A ?比较等式两边的第一行得： u1j = a1j 比较等式两边的第一列得： ?比较等式两边的第二行得： 比较等式两边的第二列得： ( j = 1,…, n ) ( i = 2,…, n ) ( j = 2,…, n ) ( i = 3,…, n ) U 的第一行 L 的第一列 U 的第二行 L 的第二列 * 计算 LU 分解 第 k 步：此时 U 的前 k-1 行和 L 的前 k-1 列已经求出 直到第 n 步，便可求出矩阵 L 和 U 的所有元素。 比较等式两边的第 k 行得： ( j = k, …, n ) 比较等式两边的第 k 列得： ( i = k+1, …, n ) * LU 分解算法 算法 ：(LU 分解 ) for k = 1 to n end j =