数值分析5剖析.ppt

? Gauss消去法的舍入误差 ? 等价问题 舍入误差分析方法 1. 向前误差分析方法：按照所执行的运算次序而估计舍入误差积累的界限. 这种方法的好处是估计比较准确，但对复杂算法（如Gauss消去法）一般难以进行。 2. 向后误差分析方法：将实际计算过程的误差转换为关于原始数据的误差。 对Gauss消去法,其具体提法是： 对于系数矩阵扰动： 求解 时, 所得的计算解 是 如下扰动方程 的精确解。 Gauss消去法由两个独立算法所组成: 一是对A做LU分解; 二是求解三角方程组. 这两个独立运算均会产生舍入误差. 如果是选主元的Gauss消去法,由于只增加了 矩阵行、列的交换，并不产生新的舍入误差，因此并不影响误差分析。 由上可见，对Gauss消去法， 只要 不是很大，则Gauss消去法是 数值稳定的 结论：条件数越大，扰动对解的影响越大. 作业: P173-175 题17， 20, §4 解线性方程组的迭代法 将 等价 改写为 形式，建立迭代 　 。从初值 出发，得到序列 。 研究 内容： ? 如何建立迭代格式？　 ? 收敛速度？ ? 向量序列的收敛条件？ ? 误差估计？ 思 路： Jacobi 法和 Gauss - Seidel 法 ? Jacobi 迭代方法 写成矩阵形式： A = -L -U D B Jacobi 迭代阵 例 求 的Jacobi迭代格式。 解 算法: Jacobi迭代方法 求解 给定初始值 . Input: 方程和未知量的个数 n; 矩阵元素 a[ ][ ]; 元素 b[ ]; 初始近似值 X0[ ]; 误差容限 TOL; 最大迭代次数 Nmax. Output: 近似解 X[ ]或失败的信息. Step 1 Set k = 1; Step 2 While ( k ? Nmax) do steps 3-6 Step 3 For i = 1, …, n Set ; /* 计算 xk */ Step 4 If then Output (X[ ]); STOP; /* 成功 */ Step 5 For i = 1, …, n Set X 0[ ] = X [ ]; /* 更新 X0 */ Step 6 Set k ++; Step 7 Output (超过最大迭代次数); STOP. /* 失败 */ 如果 aii = 0，怎么办? 迭代过程中，A 的元素 不改变，故可以事先调整好 A 使得 aii ? 0，否则 A不可逆。 必须等X(k)完全计算 好了才能计算X(k+1)，因此 需要两组向量存储。 ? Gauss - Seidel 迭代方法 … … … … 只存一组向量即可。 写成矩阵形式： B Gauss-Seidel 迭代阵 例 求方程组 的Gauss—Seidel迭代格式。 ? 松弛法 换个角度看Gauss - Seidel 方法： 其中ri(k+1) = 残量 相当于在 的基础上加个余项生成 。 下面令 ，希望通过选取合适的 ? 来加速收敛，这就是松弛法 。 ii k i k i k i a r x x ) 1 ( ) ( ) 1 ( + + + = w 0 ? 1 低松弛法 ? = 1 Gauss - Seidel 法 ? 1 (渐次)超松弛法 写成矩阵形式： 松弛迭代阵 例 求方程组 的SOR迭代格式。 解 的收敛条件 充分条件: ||B|| 1 必要条件: →0 ? 迭代法的收敛性 对任意非零向量 成立 定理 设 存在唯一解，则从任意