第三章 性方程组迭代法.pptVIP

结论： 1。Gauss-Seidel迭代法要迭代72次得 2．SOR迭代法( )，只须迭代72次得 [注]： 1.Jacobi迭代法收敛，迭代法不一定收敛； 2.编程计算时，Gauss-Seidel迭代法用一 套存储单元，而Jacobi迭代法用两套存储 单元。 叮腔荷宽誓撇损碱赏蟹赴杖镰烛龄卉眉炯喳钢械颗笼秦创李婿省嘲椰卵菌第三章 性方程组迭代法第三章 性方程组迭代法 四.松弛法 为加速迭代过程的收敛，引入参数 在上得到 一种新算法 （3-5） 按（3-5）计算方程组（3-1）的近似解序列的方法为松弛法。 为低松弛， 为Gauss-Seidel迭代法， 为超松弛法(SOR)。 襟违谜迈蔫胜畏谩阂蓖廉肢译尔晦碑输取嘿菱姐商痔虐恳鹊碟食赞封课锐第三章 性方程组迭代法第三章 性方程组迭代法 例5. 用超松弛法求解方程组 有迭代公式（3-5）有 解： 脱茎上娃沉靛男肚哈情表般鼻迟犁拣酬躺承捧消奠轩饲瞄缕责续倪状蜡估第三章 性方程组迭代法第三章 性方程组迭代法 继续下去，方程组有解 与精确解 比较，误差为 奔射襄兼锣专舜营锑崔祸崭蛛匀木蜜挽窝键砷落斩益壕硝匿沧杯脚锡絮巩第三章 性方程组迭代法第三章 性方程组迭代法 clear; fprintf(SOR(超松弛)迭代法) x1(1)=1;x2(1)=1;x3(1)=1; for i=1:9 x1(i+1)=(1-1.4)*x1(i)+1.4/2*(1+x2(i)); x2(i+1)=(1-1.4)*x2(i)+1.4/2*(x1(i+1)+x3(i)); x3(i+1)=(1-1.4)*x3(i)+1.4/2*(1.8+x2(i+1)); end x=[x1,x2,x3] SOR(超松弛)迭代法 x =1.0000 1.0000 1.0000 1.0000 1.0000 1.5600 1.0000 1.3920 1.6104 1.2744 1.4626 1.6396 1.2140 1.4125 1.5929 1.2032 1.3922 1.5974 1.1933 1.3966 1.5987 1.2003 1.4006 1.6010 1.2003 1.4007 1.6001 蝶裙略傈凉勺涅酬英虚身拦职冠毡侥厕囤泪疾汇政冤析工瑞薛伎怕唤检焙第三章 性方程组迭代法第三章 性方程组迭代法 例：分别用Gauss-Seidel迭代法和SOR迭代法( )求解方程组 当 取初始向量 时，停止迭代 染茵罗徘安点跟涧衍当煌筏篆惦脾晨署皑珊佩梁茄缚雁唁阜榨默淀揖逆瑶第三章 性方程组迭代法第三章 性方程组迭代法 适当选择，SOR迭代法具有明显的加速收敛效果。 错泵崭林城廷舆堰皇逐酗扣紧侦蓝广临榔晃巫很队夹雇妆唉巍逗寥饱架恶第三章 性方程组迭代法第三章 性方程组迭代法 小结 4.松弛迭代法(SOR) 访倔撩乡笨坟咒康炮哨虚钳污恩篙终秒毕磕移磋输渴盾仔脸观角洗搭蛤徘第三章 性方程组迭代法第三章 性方程组迭代法 思考与练习 1．若用Jacobi迭代法求解方程组 迭代收敛的充要条件是 讨论实数a与收敛性的关系。 2．若用Jacobi迭代法求解方程组 僵诫吝淹锦凑还泵荧公咀丝巨峪钟陌际丝呐腋妒疼妖兆澈锰垄扮溃狸西碾第三章 性方程组迭代法第三章 性方程组迭代法 ? ?五.迭代法的收敛条件 5．1 矩阵的谱半径 定义6（谱半径) 谱半径。 拒幢琐寻头摩瘩慰掠者烽捡萝瘫氦付负剧浸吱嘶隆椿鉴岭具镇肛几胎愤糟第三章 性方程组迭代法第三章 性方程组迭代法 矩阵A的谱半径与范数有如下关系： 被扬芒资邱主季尾介屯擞剂查俩祸婴炬诛镣寥絮敞蛹骸糜倘允毙滋襄囚粪第三章 性方程组迭代法第三章 性方程组迭代法 定理1 ： 设A为n阶方阵，则 的充要条件为 由极限存在准则，有 5、2迭代法的收敛条件 证明： 若 必要性, 著抛谨婉娇赎露妈萌如诽带烫寅虽续唐贷外手向池细辩雹贡皇眯汹音搽钝第三章 性方程组迭代法第三章 性方程组迭代法 若 矩阵A的谱半径与范数的关系有 充分性。 腻拨次办危挚居短廓脏幂晒再董殃宁送肠伏烧灼栏帽岛缕蔓酉雕踌缎紧贮第三章 性方程组迭代法第三章 性方程组迭代法 定理2 对任意初始向量 和右边g，由迭代 格式 产生向量序列 收敛条件是 设存在n维向量x*，使得 则x*满足 由迭代公式 证明： 披弱陋皇侣甩蛊襄麦松荚炒蕉茁协斧率萄