第三章 线性代数方程组的直接解法1.ppt

Cholesky分解公式 存储？？？ 因对称性无需存储 Step1 Step2 Step3 Stepn 的计算过程: 逐 列 计 算 元素 仍然存放在矩阵 的相应位置上 矩阵 分解的实际计算公式(算法3.3.1)： for for for for 例5：用Cholesky分 解法求解下列方程组 解： 系数矩阵为 Step1 Step2 Step3 求解方程组 求解方程组 ? Cholesky分解法求解方程组中需说明的几个问题 ?工作量：约为 分解的一半； ?不必选主元： 的正定性和算法的稳定性 ?稳定性：是数值稳定的； ?缺陷：存在开平方运算。 ?改进方法： 分解 改进的平方根法 ?改进的Cholesky法（改进的平方根法） 求解方程组 等价方程组 三对角线性方程组 其中 三. 三对角线性方程组的追赶法 给定矩阵 ，如果满足： 且 时， 则称 为上半带宽为 ，下半带宽为 的带状矩阵， 称为带状方程组； 如果 ，则称 为 的半带宽， 并称之为等带宽方程组； 为 的总带宽。 例如 上半带宽为2，下半带宽为1 总带宽为3 ?一些基本概念和理论： （保带状结构定理） 设 为上半带宽为 ，下半带宽为 的带状矩阵， 且其顺序主子式 ，则 有唯一的三角分解 ， 其中 是下半带宽 为 的单位下三角阵， 是上半带宽为 的上三角阵。 证明可根据前面讲过的三角分解公式 保带状结构定理说明：矩阵的三角分解中， 和 带外元素为零，因此不必计算，且不必参加求和运算 根据保带状结构定理，系数矩阵可作如下三角分解： ?三对角矩阵 分解的计算公式： ?方程组求解的计算公式： ?解方程组 ?解方程组 “追”的过程 “赶”的过程 例2：用追赶法求解三对角方程组 , 其中: 求解方程组 求解方程组 * 因为 * 因为 * 第三章 线性方程组的直接解法 ?线性方程组： ?线性方程组的求解： Cramer法则: 所需乘除法的运算量大约为(n+1)!+n n=20时，每秒1亿次运算速度的计算机要算30多万年！ 在没有舍入误差的情况下，经过有限次 运算可以得到方程组的精确解的方法。 直接法 选取初始值，逐步逼近的方法。迭代法 具有占存储单元少，程序设计简单， 原始系数矩阵在迭代过程中不变等优点， 但存在收敛性及收敛速度等问题 迭代法 基于分裂的经典迭代法； Krylov子空间迭代法。 三角形方程组 三角分解法 高斯消去法 向量和矩阵范数 线性方程组的扰动分析 本章主要内容： §3.1 三角形方程组 一、 三角形方程组的解法 ?考虑下三角形方程组 的计算公式为： 算法3.1.1 下三角形方程组的前代法： for end ?考虑上三角形方程组 的计算公式为： §3.2 高斯消去法 一、 高斯消去法 对方程组，作如下的变换，解不变 ①交换两个方程的次序 ②一个方程的两边同时乘以一个非0的数 ③一个方程的两边同时乘以一个非0数，加到另一个方程 ? 高斯消去法思想： 思路 首先将增广矩阵化为上三角阵，再回代求解 。 = 因此，对应的对增广矩阵(A，b)，作如下的变换，解不变 ①交换矩阵的两行 ②某一行乘以一个非0的数 ③某一个乘以一个非0数，加到另一行 高斯消去法就是对增广矩阵作上述行的变换，变为 上三角 方程组，而后求解 设 Ax=b. 记[A b]:=[A(1) b(1)] 1、第一次消元。设 ? 高斯消去法过程： 设第k-1次消元得A(k)x=b(k) 其中 最后 Gauss消元法的可行条件为： 1. 有些方程组，不能用Gauss消元求解 2. 如果某个 很小的话，会引入大的误差 ? 高斯消元法可行条件： 当且仅当A的所有顺序主子式不为0，即 二、 列选主元高斯消去法 在Gauss消元第k步之前，做如下的事情： 若 则交换k行和j行 注: 行的交换有效地克服了Gauss消去地缺陷 例:求解方程组 请用4位浮点数计算方程组的解。 记 为下三角矩阵, 为上三角矩阵. §3.3 三角分解法 ?基本思想: ?在矩阵计算中，矩阵还有很多其它的分解，比如QR分解，奇异值分解，满秩分解等。 ?分解算法推导: 比较第1行： 比较第1列： L为单位下三角，U为上三角 一. Doolittle分解法 ?分解形式: 存储？？？ 比较第2行： 比较第2列： 存储？？？ 比