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

石家庄铁道学院数理系第1页，共28页，星期日，2025年，2月5日对方程组做等价变换如：令，则则，我们可以构造序列若同时：所以，序列收敛与初值的选取无关第2页，共28页，星期日，2025年，2月5日定义：（收敛矩阵）定理：即：矩阵B为收敛矩阵当且仅当B的谱半径1由知，若有某种范数则迭代收敛.第3页，共28页，星期日，2025年，2月5日1.Jacobi迭代法第4页，共28页，星期日，2025年，2月5日格式很简单：第5页，共28页，星期日，2025年，2月5日2.Gauss－Seidel迭代法在Jacobi迭代中，使用最新计算出的分量值第6页，共28页，星期日，2025年，2月5日迭代矩阵记A=-L-UD第7页，共28页，星期日，2025年，2月5日易知，Jacobi迭代有第8页，共28页，星期日，2025年，2月5日迭代矩阵第9页，共28页，星期日，2025年，2月5日JacobiiterationGauss-Seideliteration计算x(k+1)时需要x(k)的所有分量，因此需开两组存储单元分别存放x(k)和x(k+1)计算xi(k+1)时只需要x(k)的i+1~n个分量，因此x(k+1)的前i个分量可存贮在x(k)的前i个分量所占的存储单元，无需开两组存储单元.第10页，共28页，星期日，2025年，2月5日迭代公式:例用Gauss-seidel迭代法解方程组Ax=b计算结果:第11页，共28页，星期日，2025年，2月5日3逐次超松弛迭代法(SOR)记则可以看作在前一步上加一个修正量。若在修正量前乘以一个因子，有对Gauss－Seidel迭代格式整理得引入松弛因子第12页，共28页，星期日，2025年，2月5日写成分量形式，有第13页，共28页，星期日，2025年，2月5日迭代矩阵第14页，共28页，星期日，2025年，2月5日SOR方法收敛的快慢与松弛因子?的选择有密切关系.但是如何选取最佳松弛因子,即选取?=?*,使?(B?)达到最小,是一个尚未很好解决的问题.实际上可采用试算的方法来确定较好的松弛因子.经验上可取1.4?1.6.第15页，共28页，星期日，2025年，2月5日第16页，共28页，星期日，2025年，2月5日第17页，共28页，星期日，2025年，2月5日第18页，共28页，星期日，2025年，2月5日4迭代法的收敛性定义设有矩阵序列及，如果则称收敛于，记为第19页，共28页，星期日，2025年，2月5日一些关于收敛的定义及定理定理定理设，则其中为的谱半径。定理(迭代法基本定理）设有方程组对于任意初始向量及任意，解此方程组的迭代法收敛的充要条件是第20页，共28页，星期日，2025年，2月5日石家庄铁道学院数理系