chapter2解线性方程组的直接法.ppt

第二章 解线性方程组的直接法 深圳大学计算机系 线性方程组的应用 递归方程求解 — 通解代入初始数据，解线性方程组 实际应用中 — 经济（线性规划） 、科学计算等 对线性方程组的研究，中国比欧洲至少早1500年，记载在公元初《九章算术》方程章中。 下三角形方程组的解—回代法 三角形方程组求解例 Guass消元法 Guass消元法 Guass消元法 Guass消元法 Guass消元法例 Guass消元法效率及不足 程序实现：二维数组存储增广矩阵 一维数组x 复杂度分析：O(n3)? O(n3/3) Guass消元法的不足: 2.3 Gauss列主元消元法 Guass列主元消元法例 Guass列主元消元法例 Guass列主元消元法的程序实现 消元法的计算中断 关于主元消元法 列主元算法在实际中应用较广。 “全面选主元”算法： 不带行交换的消元过程 不带行交换的消元过程 2.5 直接三角分解法 直接三角分解法 连续消元等于LA=U LU分解 Doolittle分解：L — 单位下三角矩阵， U — 上三角矩阵 Crout分解： L — 下三角矩阵， U — 单位上三角矩阵 定理2.1 任给矩阵A?Rn?n,若A的各阶顺序主子式不 为零，则A的LU分解存在且唯一（充要条件， 见P37）。 Doolittle分解中LU的计算 Doolittle分解 Doolittle分解 Doolittle分解程序实现 程序实现：二维数组存储增广矩阵(A，b)； L、U存储在A的对应位置； Ly = b, y存储在b的相应位置； Ux = y, x存储在b的相应位置； 复杂度分析：O(n3) 紧凑格式的Doolittle: Doolittle分解程序实现 说明:若计算没有舍入误差,用Gauss消去法和Doolittle法求得的结果应相同。但在计算机上计算,能产生不同的结果。 程序中可以用双精度运算.当完成内积计算后再舍入成单精度. 部分选主元的Doolittle分解 A = LU的计算公式： 若uii(i=1,2,..,n)(主元素）为小主元，则应加 入选主元过程，即列主元的Doolittle分解。 紧凑格式的列主元Doolittle法 紧凑格式的列主元Doolittle法步骤 紧凑格式的列主元Doolittle法例 紧凑格式的列主元Doolittle法例 Crout分解中LU的计算 Crout分解 2.6 追赶法 特殊方程组—带状矩阵：n阶矩阵A=(aij)的元素aij 若满足对小于n的正整数p和q ，有ji+q及 ij+p时，aij=0，则称A为具有上带宽为q,下带宽 为p的带状矩阵。 带状矩阵的LU分解 定理2.3 如果上带宽为q、下带宽为p的n阶带状矩 阵有Doolittle分解A=LU ，则L是下带宽为p的单位 下三角矩阵，U是上带宽为q的上三角矩阵。 三对角矩阵的LU分解—追赶法 AX=b，假设A为三对角矩阵，即p=q=1的带状矩阵。 根据定理2.3知A的Crout分解格式为： 三对角矩阵的LU分解—追赶法 根据矩阵乘积有： 三对角矩阵的LU分解—追赶法 解方程Ly=b有： 三对角矩阵的LU分解—追赶法 三对角矩阵的LU分解的唯一性 如果A为三对角线矩阵，且满足： 追赶法例 例8：用追赶法求解线性方程组。并说明分解是否唯一，为什么？ 2.7 平方根法 特殊方程组Ax=b,其中A为对称正定阵（A的各阶顺序 主子式全部大于0）。 定理2.4 设A是对称正定阵，则存在对角线元素全 为正数的下三角形矩阵L，使A=LLT存在且唯一。称这 种分解为Cholesky分解。 证明略：根据定理2.1及A的对称性质证。 参见P51~52。 Cholesky分解中L的计算 Cholesky分解中y,x的计算 Cholesky分解例 例9：用平方根法解对称正定方程组： 平方根法实现 平方根法中由于L只有下三角元素，且计算只与当前 aij有关，因此同一般的LU分解，也可采用紧凑格 式。 若A为对称阵，方程组Ax=b的求解思路同对称正定 阵，分解格式为A=LDLT。详细介绍请参见其它参考 书。 本章总结 本章主要学习求解线性方程组Ax=b的各种方法。 要求掌握： 解下三角或上三角