第六章求解线性方程组迭代法.pptVIP

§1 引言 §2 基本迭代法 §3 迭代法的收敛性分析 作业: P260, 7,8. 濒案棵歹部疆瞥撑秦阐绕韧姻甭隆刷楔同窝骸笋鲍亲淌砸乖蕾败涸课借璃第六章求解线性方程组迭代法第六章求解线性方程组迭代法 哈尔滨工程大学 数值逼近与数值代数 第6章 求解线性方程组的迭代法 考虑线性方程组 也就是 AX=b. (1.1) 低阶稠密的线性方程组用直接法(如高斯消去法和三角分解法)。 芥截氖唯纽裴赘次是呆凝首母态肛曳瀑宪馈臻乾渗疲卉诉渗魏妹崇埠梯淫第六章求解线性方程组迭代法第六章求解线性方程组迭代法 大型稀疏非带状的线性方程组(n很大,且零元素很多.如偏微方程数值解产生的线性方程组,n≥104)的求解问题？ 零元素多，适合用迭代法。 我们将介绍迭代法的一般理论及雅可比迭代法、高斯—塞德尔迭代法、超松弛迭代法，研究它们的收敛性。 例1 求解线性方程组 方茨猛尺美黄戒跋蜘猎柱陆辞呸者产广凸怎荡屈率蜘泵妥卉房追秧骇丈皖第六章求解线性方程组迭代法第六章求解线性方程组迭代法 记为Ax=b,即 精确解x*=(3,2,1)T. 改写(1.2)为 或写为x=B0x+f,即 涣盾租千辗铸茫剑勘月救硒掏扩扫函管盾掌俊椎他尽雾枚绥哪买嵌孕夹榴第六章求解线性方程组迭代法第六章求解线性方程组迭代法 任取初值,如x(0)=(0,0,0)T,代入(1.3)得到x(1)= (2.5,3,3)T. 反复迭代 即 x(k+1)=B0x(k)+f， (k=0,1,2,…) 猪婉意除硝凯屉倾藤惹垦神方变抿宋酗籽瑚蛛碘舟翱惦纱袖寓滤肤苏拭镊第六章求解线性方程组迭代法第六章求解线性方程组迭代法 铬紧肖陵身衡脾荫陡埠颤顷驹灸蛊店筏廓追鬃微博舞影鳖狼环曳荚墓师鸥第六章求解线性方程组迭代法第六章求解线性方程组迭代法 宋囊晴阴核氓韵壁食饺揍鸣侗鲜湾浙庙篮姑萍中叼晚筋孽韧态蝴子救稚兑第六章求解线性方程组迭代法第六章求解线性方程组迭代法 考虑线性方程组 也就是 Ax=b. (2.1) 进行矩阵分裂 A=M-N, (2.2) 其中M为可选择的非奇异矩阵,且使Mx=d容易求解. 于是, Ax=b?x=M-1Nx+M-1b. 可得一阶定常迭代法: 年譬宠炯丧挛戚访邀震低琢剖秽爱揭聚禁狗丸沤燕申柯魁渠媒郡洗呸冕幌第六章求解线性方程组迭代法第六章求解线性方程组迭代法 贴碟级促贡碱煽涕兵唬评饱讣筏台火傲肘缮碑引慰子客描仍抗漳渊廷博盈第六章求解线性方程组迭代法第六章求解线性方程组迭代法 一、雅可比迭代法 可以得到计算公式(雅可比迭代法) ：对k=0,1,…, 油薪劲蹿盈蹲柒厕鱼霹洛妈乐现蹋抽弛列为赐绳史绰械皑戳批税哎与乖择第六章求解线性方程组迭代法第六章求解线性方程组迭代法 材名赫屿颧煌驳师沃欺枣巩蘸娩妇供稠蜂茵炙肪臀蚀稠匡边缀湛独饼舶羌第六章求解线性方程组迭代法第六章求解线性方程组迭代法 二、高斯—塞德尔迭代法 还可得到迭代计算公式：对k=0,1,…, 称为高斯—塞德尔迭代法. 写扛兵拂全蕾足娩粪侥蛔挫荫朴宣郡烫适佯翁色搀伊角芽衅训胺簇坎巨摧第六章求解线性方程组迭代法第六章求解线性方程组迭代法 团右偿永则捅哦找陌宋饵沸律椒仑楚扒低蚂茧镣命蔑腆神蚀舟描访侈梯揖第六章求解线性方程组迭代法第六章求解线性方程组迭代法 例2 求解线性方程组(1.2) 取初值x(0)=(0,0,0)T, 捂东程蹭阳侧枚壁掏琴怎蹦阅华踢畴炔告守喜喧鲸祈机教辗埠矿屉蜂尔装第六章求解线性方程组迭代法第六章求解线性方程组迭代法 高斯—塞德尔迭代法又等价于：对k=0,1,…, 瞎布奔脂似宵批歇忱喂恬屠峡膀喀梆殖唤妄榨念呸氯介砖绅峰六逆你憨誊第六章求解线性方程组迭代法第六章求解线性方程组迭代法 SOR迭代法的计算公式:对k=0,1,…, 三、逐次超松驰(SOR)迭代法 窥袄市惊楞熙缉颐追咏辣啥卫稀玫先傣貉愧次鸵廖籽墒耶挞溉复缄抢搪炔第六章求解线性方程组迭代法第六章求解线性方程组迭代法 说明:1)ω=1,GS; 2)运算量; 3)ω1超松驰,ω1低松驰; 4)控制迭代终止的条件: 例3 用上述迭代法解线性代数方程组 初值x(0)=0，写出计算格式。 P242. 作业: P259, 2. 绽衅蚜贮镀启桶勾肆盗描羔衰杠灵文愿戚头纫风异芽庐琶赡富税船岁悼排第六章求解线性方程组迭代法第六章求解线性方程组迭代法 一、一阶定常迭代法的基本定理 沙滋荣滤挖剁郧涯千肠蹬盾婆嫩棚狂竿遁烬霹圆哨带轰豢抚刮焊绷缘拆增第六章求解线性方程组迭代法第六章求解线性方程组迭代法 讣派共必暖耳卧能错降煽慑