第六章 线性与非线性方程组的迭代解法2.pptVIP

* §6.3 超松弛(SOR)迭代法/*Overrelaxation Iteration*/ 一、超松弛(SOR)迭代法 思想 类似于G-S迭代法的改进方法，利用第k次 迭代值和第k +1次的G-S迭代值作加权平均 G-S迭代法的计算公式 ： 作加权平均 ?超松弛(SOR)迭代法的分量形式： 其中 称为松弛因子； 时即为G-S迭代 ?SOR迭代法的迭代矩阵： 记 超松弛(SOR)迭代法的迭代公式 迭代矩阵 其中 例1：写出SOR迭代法求解下列方程组的迭代格式 解： SOR迭代法的迭代公式 取初始向量 选取不同的 值进行计算，结果见下表 (2.9999993 4.0000001 -4.9999996） 12 0.00001 (2.9999853 4.0000031 -4.9999935） 10 0.0001 (2.9997451 4.0000653 -4.9998924） 8 0.001 (3.0000210 3.9999820 -5.0000046) 23 0.00001 (3.0001673 3.9998567 -5.0000368) 18 0.0001 (3.0020191 3.9982705 -5.0004444) 12 0.001 (3.0000186 3.9999845 -5.0000039) 21 0.00001 (3.0001952 3.9998374 -5.0000407) 16 0.0001 (3.0012790 3.9989342 -5.0002665) 12 0.001 (2.9995106 4.0017780 -5.0027919） 151 0.001 (3.0006104 4.0001741 -5.0007434) 13 0.001 迭代 次数 要求 精度 方 程 组 的 近 似 解 的值 ?向后的SOR迭代法： 记 SOR迭代法的向量形式： 迭代矩阵： 向后SOR迭代法的分量形式： ?对称超松弛迭代法（SSOR） ： 例2：写出下列方程组向后的SOR和SSOR的迭代格式 解： 向后的SOR迭代法的迭代格式 SSOR迭代法的迭代格式 (3.0000012 3.9999988 -5.0000027） 11 0.00001 (2.9999991 4.0000051 -4.9999831） 10 0.0001 (2.9998426 4.0003635 -4.9995660） 8 0.001 bSOR SSOR 所用 方法 (3.0000099 3.9999843 -5.0000035） 28 0.00001 (3.0000939 3.9998514 -5.0000334) 23 0.0001 (3.0008900 3.9985916 -5.0003161) 18 0.001 迭代 次数 要求 精度 方 程 组 的 近 似 解 取初始向量 选取 进行计算，结果见下表 二、SOR迭代法的收敛性： 其中 对 ， 设其对角元皆非零， 则对所有实数 ，有 证明： 设迭代矩阵 的特征值为 推论6.3.1 如果求解方程组 的SOR法收敛，则 设 为对称正定矩阵，且 则求解方程组 的SOR法收敛。 证明： 迭代矩阵 设 的任意特征值为 ，相应的特征向量为 ，则 上式两边与 作内积 因为矩阵 正定，所以 也正定： 记 则 分子减去分母： 因为矩阵 正定， 设 为对称矩阵，且所有对角元 ，若求解方程组 的SOR法 收敛，则 正定，且 。 证明： 记 由 得 考察二次型 注意到 说明二次型 是单调递减趋于零 如果 不是正定的，则存在非零向量 ，满足 且有 矛盾