6.1迭代法和敛散性及其MATLAB程序6.1.1迭代法敛散性的判别及其MATLAB程序用谱半径判别迭代法产生的迭代序列的敛散性的MATLAB主程序functionH=ddpbj(B)H=eig(B);mH=norm(H,inf);ifmH>=1disp('请注意：因为谱半径不小于1，所以迭代序列发散,谱半径mH和B的所有的特征值H如下：')elsedisp('请注意：因为谱半径小于1，所以迭代序列收敛,谱半径mH和B的所有的特征值H如下：')endmH6.2雅可比（Jacobi）迭代及其MATLAB程序6.2.1雅可比迭代的收敛性及其MATLAB程序判别雅可比迭代收敛性的MATLAB主程序functiona=jspb(A)[nm]=size(A);forj=1:ma(j)=sum(abs(A(:,j)))-2*(abs(A(j,j)));endfori=1:nifa(i)>=0disp('请注意：系数矩阵A不是严格对角占优的，此雅可比迭代不一定收敛')returnendendifa(i)<0disp('请注意：系数矩阵A是严格对角占优的，此方程组有唯一解，且雅可比迭代收敛')end例6.2.1用判别雅可比迭代收敛性的MATLAB主程序，判别由下列方程组的雅可比迭代产生的序列是否收敛？（1）（2）解（1）首先保存名为jspb.m的M文件，然后在MATLAB工作窗口输入程序>>A=[10-1-2;-110-2;-1-15];a=jspb(A)37.运行后输出结果请注意：系数矩阵A是严格对角占优的，此方程组有唯一解，且雅可比迭代收敛a=-8-8-1（2）在MATLAB工作窗口输入程序>>A=[10-1-2;-110-2;-1-10.5];a=jspb(A)运行后输出结果请注意：系数矩阵A不是严格对角占优的，此雅可比迭代不一定收敛a=-8.0000e+000-8.0000e+0003.5000e+0006.2.2雅可比迭代的两种MATLAB程序（一）雅可比迭代公式的MATLAB程序用雅可比迭代解线性方程组的MATLAB主程序functionX=jacdd(A,b,X0,P,wucha,max1)[nm]=size(A);forj=1:ma(j)=sum(abs(A(:,j)))-2*(abs(A(j,j)));endfori=1:nifa(i)>=0disp('请注意：系数矩阵A不是严格对角占优的，此雅可比迭代不一定收敛')returnendendifa(i)<0disp('请注意：系数矩阵A是严格对角占优的，此方程组有唯一解，且雅可比迭代收敛')endfork=1:max1kforj=1:mX(j)=(b(j)-A(j,[1:j-1,j+1:m])*X0([1:j-1,j+1:m]))/A(j,j);endX,djwcX=norm(X'-X0,P);xdwcX=djwcX/(norm(X',P)+eps);X0=X';X1=A\b;if(djwcXwucha)&(xdwcX>wucha)disp('请注意：雅可比迭代次数已经超过最大迭代次数max1')enda,X=X;jX=X1',例6.2.2用范数和判别雅可比迭代的MATLAB主程序解例6.2.1中的方程组，解的精度为0.0...
