数值分析迭代法的收敛性分析.pptVIP

第1页，共18页，星期日，2025年，2月5日平面点列:(x1,y1),(x2,y2),···,(xk,yk),······Xk∈Rn:X1,X2,···,Xk,·······利用向量范数等价性,对任意范数||·||2/18第2页，共18页，星期日，2025年，2月5日原方程:AX=b记?(k)=X(k)–X*(k=0,1,2,3,······)则有?(k+1)=B?(k)?(k)=B?(k-1)(k=1,2,3,······)计算格式:X(k+1)=BX(k)+fX(k+1)–X*=B(X(k)–X*)设方程组的精确解为X*,则有X*=BX*+f?3/18第3页，共18页，星期日，2025年，2月5日(1)?(k)=B?(k-1)=B2?(k-2)=···=Bk?(0)(2)迭代格式X(k+1)=BX(k)+f收敛！4/18第4页，共18页，星期日，2025年，2月5日证:由?(k)=B?(k-1),得||?(k)||≤||B||||?(k-1)||(k=1,2,3,······)所以命题若||B||1,则迭代法X(k+1)=BX(k)+f收敛||?(k)||≤||B||k||?(0)||||B||15/18第5页，共18页，星期日，2025年，2月5日矩阵B的谱设n阶方阵B的n个特征值为:则称集合为B的谱.记为chB矩阵B的谱半径注1:当B是对称矩阵时,||B||2=?(B)注2:对Rn×n中的范数||·||,有?(B)≤||B||特征值取模最大6/18第6页，共18页，星期日，2025年，2月5日定理4.1迭代法X(k+1)=BX(k)+f收敛?谱半径ρ(B)1证:对任何n阶矩阵B都存在非奇矩阵P使B=P–1JP其中,J为B的Jordan标准型其中,Ji为Jordan块7/18第7页，共18页，星期日，2025年，2月5日其中,λi是矩阵B的特征值,由B=P–1JPBk=(P–1JP)(P–1JP)···(P–1JP)=P–1JkP迭代法x(k+1)=Bx(k)+f收敛=(i=1,2,···,r)(i=1,2,···,r)谱半径?(B)18/18第8页，共18页，星期日，2025年，2月5日注1:AX=b?X=BX+f?(I–B)X=f?X=(I–B)-1f注2:若则(I-B)-1=I+B+B2+······+Bk+······事实上(I-B)(I+B+B2+······+Bk)=I–Bk+1注3:X(k)=BX(k-1)+f=B(BX(k-2)+f)+f=····=BkX(0)+(I+B+····+Bk-1)f≈(I–B)-1f9/18第9页，共18页，星期日，2025年，2月5日Ans=1.2604e-005例线性方程组AX=b,分别取系数矩阵为试分析Jacobi迭代法和Seidel迭代法的敛散性D=diag(diag(A1));B1=D\(D-A1);max(abs(eig(B1)))(1)A1=[1,2,-2;1,1,1;2,2,1]10/18第10页，共18页，星期日，2025年，2月5日DL=tril(A1)B1=DL\(DL-A1)max(abs(eig(B1)))Ans=2(2)A2=[2,-1,1;1,1,1;1,1,-2]D=d