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

1第五章解线性方程组的直接方法计算方法—— 矩阵三角分解法第一页，共十六页。 2本讲内容 一般线性方程组 LU 分解与 PLU 分解 对称正定线性方程组 平方根法－－Cholesky 分解 对角占优三对角线性方程组 追赶法第二页，共十六页。 3LU 分解将一个矩阵分解成结构简单的三角形矩阵的乘积矩阵的三角分解矩阵的 LU(Doolittle) 分解矩阵的 LDR 分解克洛脱 (Crout) 分解第三页，共十六页。 4计算 LU 分解利用矩阵乘法直接计算 LU 分解L ? U = A?比较等式两边的第一行得：u1j = a1j比较等式两边的第一列得：?比较等式两边的第二行得：比较等式两边的第二列得：( j = 1,…, n )( i = 2,…, n )( j = 2,…, n )( i = 3,…, n )U 的第一行 L 的第一列 U 的第二行 L 的第二列 第四页，共十六页。 5计算 LU 分解第 k 步：此时 U 的前 k-1 行和 L 的前 k-1 列已经求出直到第 n 步，便可求出矩阵 L 和 U 的所有元素。比较等式两边的第 k 行得：( j = k, …, n )比较等式两边的第 k 列得：( i = k+1, …, n )第五页，共十六页。 6LU 分解算法算法 ：(LU 分解 )for k = 1 to nendj = k, …, ni = k+1, …, nMatlab程序参见：ex51.m乘除法运算量：(n3 - n)/3为了节省存储空间，通常用 A 的绝对下三角部分来存放 L (对角线元素无需存储)，用 A 的上三角部分来存放 U 第六页，共十六页。 7PLU 分解矩阵的 PLU 分解for k = 1 to nendi = k, k+1, …, nj = 1, 2, …, ni = k+1, …, nj = k+1, …, nMatlab程序：上机练习 第七页，共十六页。 8Cholesky 分解 对称正定矩阵的三角分解－－Cholesky 分解定理：设 A 是对称矩阵，若 A 的所有顺序主子式都不为 0，则 A 可唯一分解为其中 L 为单位下三角阵，D 为对角矩阵A = LDLT定理：（Cholesky分解）若 A 对称正定，则 A 可唯一分解为其中 L 为下三角实矩阵，且对角元素都大于 0A = LLT第八页，共十六页。 9计算 Cholesky 分解 Cholesky 分解的计算直接比较等式两边的元素 计算公式第九页，共十六页。 10Cholesky 分解算法for j = 1 to nendi = j +1, …, n算法 ：(Cholesky 分解 )第十页，共十六页。 11平方根法A 对称正定算法 ：(解对称正定线性方程组的平方根法 )计算 A 的 Cholesky 分解解方程：Ly = b 和 LTx = yi = 2, 3, …, ni = n-1, …, 2, 1 第十一页，共十六页。 12改进的 Cholesky 分解 计算公式 改进的 Cholesky 分解第十二页，共十六页。 13改进的 Cholesky 分解for j = 1 to nendi = j +1, …, n算法 ：(改进的 Cholesky 分解 ) 优点：避免开方运算第十三页，共十六页。 14改进的平方根法A 对称正定算法 ：(解对称正定线性方程组的改进的平方根法 )计算 改进的 Cholesky 分解解方程：Ly = b 和 DLTx = yi = 2, 3, …, ni = n-1, …, 2, 1 第十四页，共十六页。 15追赶法 对角占优的三对角矩阵的 LU 分解 计算公式i = 2, 3, …, n-1第十五页，共十六页。 16追赶法A 三对角矩阵（对角占优）算法 ：(追赶法 )i = 2, 3, …, ni = n-1, …, 2, 1 i = 2, 3, …, n-1 运算量：5n-42n – 3 次2n 次n – 1 次第十六页，共十六页。