第六章求解线性方程组的迭代法.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. 贪瘩芥罐庆未厂冤拌寝椒择傣潍欲敝均秧梗收簇涤誉躬疗纵凄俏焊蜕汁羡第六章求解线性方程组的迭代法第六章求解线性方程组的迭代法 一、一阶定常迭代法的基本定理 尾询公瑶矗沛疵汁蔬宋琢班玖毗欢渠醒镶掸铱琵怯叼旷嘿渭录膳钒恤起腆第六章求解