求解线性方程组的方法.ppt

设方程组AX=b, 其中，A=(aij) 为非奇异阵， x=(x1, x2, …, xn)T, b=(b1, b2, …, bn)T. 假设已算出 x(k) ， §8.3 超松驰迭代法SOR方法 2. SOR算法的构造 ω称为松弛因子 利用高斯-塞德尔迭代法得: §8.3 超松驰迭代法SOR方法 2. SOR算法的构造 (基于G-S迭代) 解方程组AX=b的逐次超松弛迭代公式： 显然，当取ω=1时，上式就是高斯-塞德尔迭代公式. §8.3 超松驰迭代法SOR方法 2. SOR算法的构造(基于Jacobi迭代) 得到解方程组 AX=b 的逐次超松弛迭代公式： 显然，上式就是 基于Jacobi 迭代的 SOR 方法. 下面令 ， 希望通过选取合适的 ? 来加速收敛，这就是松弛法 。 3. SOR算法的进一步解释 SOR方法 其中ri(k+1) = 相当于在 的基础上加个余项生成 。 0 ? 1 低松弛法 ? = 1 Gauss - Seidel 法 2? 1 (渐次)超松弛法 利用SOR方法解方程组 SOR例题分析: 其准确解为x*={1