数值分析课件-num_3.4超松弛迭代法.pdfVIP

第三章 线性方程组迭代解法 § 3.4 超松弛迭代法(SOR) § 3.4超松弛迭代法(SOR) § 3.4 SOR) SOR （Successive Over-Relaxation） 法，即超松弛迭代法，是目前解大型线性 方程组的一种最常用的方法，是Gauss- Seidel迭代法的一种加速方法。  一、SOR 法迭代公式  例3.6 用SOR法求解线性方程组  二、SOR法的收敛性  SOR法收敛与收敛速度有关定理  SOR法分类与现状 返回引用 一 SOR法迭代公式 一 SOR法迭代公式 设线性方程组 AX=b 其中 A 非奇异，且a  0 （i=1,2,,n ） 。 ii 如果已经得到第k次迭代量x (k) 及第k+1次迭代量x (k+1) 的前i-1个 分量 （x 1 (k+1),x2 (k+1) ,,xi-1 (k+1) )， 在计算xi (k+1) 时，先用Gauss-Seidel迭代法得到 (1) 返回引用 选择参数 ω，取 (2) 把 式（1）代入式（2 ）可以综合写成： 即得超松弛法或逐次超松弛迭代法（Successive Over- Relaxation Method），简称SOR法。 返回引用 或可表示成增量的形式： 其中，  参数 ω叫做松弛因子；  若 ω=1，它就是Gauss-Seidel迭代法。 令A=D-L-U, SOR法(2)式可写成： 再整理成： (k1) 1 (k) 1 X  （D-L) [(1)DU]X (DL) b 于是可导出SOR法的矩阵形式： (k1) (k) X BX  f  其中，迭代矩阵和f为: B  （D-L)1[(1)DU]   1 f (D L) b 返回引用 例3.6 用SOR法求解线性方程组 例3.6 用SOR法求解线性方程组 解 方程组的精确解为 x=(3,4,-5) T，为了进行比较， (0) T 利用同一初值 x =(1,1,1) ，分别取 ω=1