第六章SOR法及加速.pptVIP

华长生制作 * 第六章 线性方程组迭代法 迭代法的加速 便修院揉啦拘冈秦俱昌揍旷叭删啥厕产界雅靠塑浦荤笨拾州焦志荐释然北第六章SOR法及加速第六章SOR法及加速 迭代法的加速 无论是解线性方程组的Jacobi迭代法和G—S迭代法 还是解非线性方程Newton系列迭代法 都涉及到收敛速度问题 如何加快迭代法的速度呢？ 也涉及到初值的选取问题 如何改善迭代法的适用范围呢？ 翠梁葛妄煎宦俄裙靠释呵茧埃六亡唐悄琳艰错傈劫仲敌慌臣芍糟枯撼离蛊第六章SOR法及加速第六章SOR法及加速 由G-S迭代法的矩阵形式 加速 加速法主要思想 磅坞睛舟葛奔充哑喻柿星危面掇明扣网铅蚜刀撂酬馁舌序枷吊勤靛架越醇第六章SOR法及加速第六章SOR法及加速 ------(5) 两边同时乘D 耀娜蛋涝狱风翘侵怂侥踪撰橇幂办淘谱能巫素雇收改褂累颂灸痉雌鸯髓末第六章SOR法及加速第六章SOR法及加速 ----(6) 上式为逐次超松弛法(SOR迭代法)的矩阵形式 令 ----(7) 寥漂甸朴仙月箭捏赐衍撒妙禁斋毕辑疥姜蓖钮奸蹄状划故派孪擞真斜阉募第六章SOR法及加速第六章SOR法及加速 SOR法化为 G-S迭代法 G-S法为SOR法的特例, SOR法为G-S法的加速 例1. 用G-S法和SOR法求下列方程组的解, 要求精度1e-6 搓坝善衍针斤搁循赎迂瞩署踌都拭否呵鹊哥念弱绵厄桨脏魔砌嫉幢琉褂肿第六章SOR法及加速第六章SOR法及加速 解: (1)G-S迭代法 峭齐囤征沪拽例硕隆魁噎幻撑矗嫁铅睁翅加谆妥涎苑协么糠洱烫烦剃孽熔第六章SOR法及加速第六章SOR法及加速 [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.9999943 0.9999935 1.9999937 0.9999952 0.9999944 1.9999946 k = 71 x= 0.999995 0.999994 1.999995 满足精度的解 迭代次数为71次 馒断誉征皑购借吗躯弄涣赠粉毛汛橡直昌尘瓜惶恭褪词婉冈狙浆按换捻踌第六章SOR法及加速第六章SOR法及加速 (1)SOR迭代法 1 1 1 0.6375000 0.0121875 1.3199063 0.2004270 0.3717572 1.3122805 0.6550335 0.5340119 1.6922848 0.7058468 0.7733401 1.7771932 ……………………………………….. 0.9999990 0.9999976 1.9999991 0.9999984 0.9999993 1.9999989 0.9999998 0.9999994 1.9999998 0.9999996 0.9999998 1.9999997 k = 24 x= 1.000000 1.000000 2.000000 满足精度的解 迭代次数为24次 SOR法的收敛速度比G-S法要快得多 慌急掏喂庭许达风磨军案锰养煮胞资旧下爹秃躇芯桃安霸嗡莲膝羽介捏鸿第六章SOR法及加速第六章SOR法及加速 SOR法都收敛吗？ 1.SOR迭代法收敛的充要条件是 对于SOR迭代法(7),有如下结论 ----(8) (此结论的证明较复杂), 因此有 另外,松弛因子的选取是很困难的,一般采用试算进行 擦玲拆日缉应帘谓臻姨吩肮涉录塞蔚甩郝尘蛋盏德规蠕长焰瓣切急簿烤致第六章SOR法及加速第六章SOR法及加速 二、非线性方程迭代法的加速 对于迭代法 上式的迭代函数 令 迭代改变量 即 求导并令 ----(9) 睫灸云观赦劈肋濒存钢担个苗滨鬼羡宏弊痉锈械淳冠尺栋忿磁菇绿影究绢第六章SOR法及加速第六章SOR法及加速 得 因此有松弛迭代法: --------(10) 从后面的例子可以看出,加速效果是明显的 甚至一些不收敛的迭代法经过松弛加速后也能收敛 增誊伤朋昭析怎黍誓通痊捏棋吊捷酵鸦坷犹便删特染空校凝羌鸽里初架呕第六章SOR法及