第二章解线性方程组的迭代法.pptVIP

理学院UniversityofShanghaiforScienceandTechnologyCollegeofScience上海理工大学第1页，共39页，星期日，2025年，2月5日对方程组做等价变换从某一初值x(0)出发，我们可以构造序列若同时：所以，序列收敛与初值的选取无关如令A=D-L-U，于是x=D-1(L+U)x+D-1b,第2页，共39页，星期日，2025年，2月5日定义5.1：设G为n阶方阵，若Gk?0，则称G为收敛矩阵定理：即矩阵G为收敛矩阵，当且仅当G的谱半径1由知，若有某种范数则，迭代收敛第3页，共39页，星期日，2025年，2月5日迭代法的收敛性定理：迭代法X(m+1)=GX(m)+g收敛的充分必要条件是迭代矩阵G为收敛矩阵，即G的谱半径?(G)1。定理：迭代法X(m+1)=GX(m)+g的迭代矩阵G的某种范数||G||＝q1,那么：1）对任意初值X(0)及g右端向量，迭代格式收敛于X*;2)||X(m)-X*||?qm||X(1)–X(0)||/(1-q);3)||X(m)-X*||?q||X(m)–X(m-1)||/(1-q).第4页，共39页，星期日，2025年，2月5日Jacobi迭代第5页，共39页，星期日，2025年，2月5日格式很简单：第6页，共39页，星期日，2025年，2月5日1、输入系数矩阵A和向量b，和误差控制eps2、x1={0,0,…..,0},x2={1,1,…..,1}//赋初值3、while(||A*x2-b||eps){x1=x2;for(i=0;i=n;i++){x2[i]=0;for(j=0;ji;j++){x2[i]+=A[i][j]*x1[j]}for(j=i+1;jn;j++){x2[i]+=A[i][j]*x1[j]}x2[i]=-(x2[i]-b[i])/A[i][i]}}4、输出解x2Jacobi迭代算法第7页，共39页，星期日，2025年，2月5日迭代矩阵记Jacobi迭代法的收敛性第8页，共39页，星期日，2025年，2月5日易知，Jacobi迭代有第9页，共39页，星期日，2025年，2月5日练习讨论用雅可比(Jacobi)迭代法求解下列线性方程组的收敛性。若收敛，求其解；若发散，作适当变换使其收敛并求解。第10页，共39页，星期日，2025年，2月5日G的谱半径?(G)=4.01971.Jacobi迭代不收敛。迭代矩阵为G的特征值为：?1=4.02408,?2=-2.01204?3.10115i,??1?=4.02408；??2,3?=3.69668第11页，共39页，星期日，2025年，2月5日将方程组变形，化为：第12页，共39页，星期日，2025年，2月5日G的谱半径?(G)=0.3085071.Jacobi迭代收敛。此时迭代矩阵为G的特征值分别为：0.308507,-0.154254+0.18304i,-0.154254-0.18304i第13页，共39页，星期日，2025年，2月5日收敛条件迭代格式收敛的充要条件是G的谱半径1。对于Jacobi迭代，我们有一些保证收敛的充分条件定理：若线性方程组AX=b的系数矩阵A满足下列条件之一，则Jacobi迭代收敛。①A为行对角占优阵②A为列对角占优阵③A满足④若A对称正定阵，且2D-A也为对称正定阵，则Jacobi迭代收敛。第14页，共39页，星期日，2025年，2月5日证明：②A为列对角占优阵，则AT为行对角占优阵，有＃证毕第15页，共39页，星期日，2025年，2月5日在Jacobi迭代中，使用最新计算出的分量值Gauss－Seidel迭代第16页，共39页，星期日，2025年，2月5日1、输入系数矩阵A和向量b，和误差控制eps2、x2={1,1,…..,1}//赋初值3、while(||A*x2-b||eps){