第六章SOR法及加.pptVIP

华长生制作 * 第六章 线性方程组迭代法 迭代法的加速 迭代法的加速 无论是解线性方程组的Jacobi迭代法和G—S迭代法 还是解非线性方程Newton系列迭代法 都涉及到收敛速度问题 如何加快迭代法的速度呢？ 也涉及到初值的选取问题 如何改善迭代法的适用范围呢？ 由G-S迭代法的矩阵形式 加速 加速法主要思想 ------(5) 两边同时乘D ----(6) 上式为逐次超松弛法(SOR迭代法)的矩阵形式 令 ----(7) SOR法化为 G-S迭代法 G-S法为SOR法的特例, SOR法为G-S法的加速 例1. 用G-S法和SOR法求下列方程组的解, 要求精度1e-6 解: (1)G-S迭代法 [x,k]=gauss_seidel(a,b,[1,1,1],1e-6) 1 1 1 0.7500000 0.3750000 1.5000000 0.5625000 0.5312500 1.5416667 0.6510417 0.5963542 1.6145833 0.7018229 0.6582031 1.6727431 ………………………………………. 0.9999933 0.9999923 1.9999926 0.999994