5.3超松弛迭代法.pptVIP

第五章线性方程组迭代解法 5.3 超松弛迭代法 5.3.2 超迭代法的收敛性 5.3.1 超迭代法的构造 5.3.1 超松弛迭代法的构造 经整理得 （5. 3. 1) 称此式为逐次超松弛迭代法，简记为SOR(Successive Over – Relaxation)法，其中 称为超松弛因子。当 时，（5. 3.1)就是GS法。 在很多情况小，J法和GS法收敛较慢，所以考虑GS法的改进。设计算第 个近似解 时，分量 已经算好。按GS法给出辅助量 再用参数 将 与 做加权平均，即 记A=D－L－U，(5. 3. 1)可写成矩阵形式 再整理得 (5.3.2) 其中，迭代矩阵为 (5.3.3) 例5.4 方程组 得准确解为 ，如果用 SOR跌代法（即GS法），计算公式是 如果用 的SOR迭代法，计算公式是 设 的特征值为 ，则 证 取 ，迭代7次，则 时得 时得 若继续算下去，要达到7位数字的精度， 时，要迭代34次，而 时，只需要迭代14次，显然选 收敛要快些。 按一般的迭代法收敛的理论，SOR迭代法收敛的充分必要条件是 而 与松弛因子 有关。下面讨论松弛因子 在什么范围内取值，SOR迭代法可能收敛。 定理5.7 如果解方程组 的SOR法收敛，则有 。 5.3.2 超松弛迭代法的收敛性 由于SOR法收敛，所以有 定理得证。 该定理说明，只有当松弛因子 在区间 内取值时，SOR法才可能收敛。下面给出SOR法收敛的充分条件。 定理5. 8 如果A为对称正定矩阵，且 ，则解 的SOR法收敛。 证 设 是 的一个特征值，对应特征向量 。由（5. 3. 3）可得 这里， 是实对称矩阵，所以有 。上式两边与 作内积得 （5. 3. 4） 因为A正定，D亦正定，记 ，有 。又记 ，则有 由（5. 3. 4）有 因A正定， ， ，所以 即 ，从而 ，SOR方法收敛。定理得证。 当 时，SOR法就是GS法，所以上面的定理说明，当系数矩阵是对称正定矩阵时，GS法收敛。 对于例5. 4所给出的方程组，其系数矩阵是对称正定的，因此对 和 的SOR迭代法都收敛。 例5. 5 设矩阵A非奇异，求证用GS法求解方程组 时是收敛的。 对 ，由A非奇异知 ，从而 即 是对称正定的，因此，用GS法求解方程组 是收敛。 证 引入超松弛迭代法的想法是希望能找到最优的松