Hx=b H为Hilbert矩阵.pptVIP

Hx=b H为Hilbert矩阵分别构造Jacobi迭代和Guass-Seidel迭代，求它们是否收敛 Jacobic迭代法 function [n,x]=jacobi(A,b,X,nm,w) %用雅克比迭代法求解方程组Ax=b %输入：A为方程组的系数矩阵，b为方程组右端的列向量，X为迭代初值构成的列向量，nm为最大迭代次数，w为误差精度 %输出：x为求得的方程组的解构成的列向量，n为迭代次数n=1;m=length(A);D=diag(diag(A)); %令A=D-L-U,计算矩阵DL=tril(-A)+D; %令A=D-L-U,计算矩阵LU=triu(-A)+D; %令A=D-L-U,计算矩阵UM=inv(D)*(L+U); %计算迭代矩阵g=inv(D)*b; %计算迭代格式中的常数项%下面是迭代过程while n=nm x=M*X+g; %用迭代格式进行迭代 if norm(x-X,2)w disp(迭代次数为);n disp(方程组的解为);x return; %上面：达到精度要求就结束程序，输出迭代次数和方程组的解 end X=x;n=n+1;end%下面：如果达到最大迭代次数仍不收敛，输出警告语句及迭代的最终结果（并不是方程组的解）disp(在最大迭代次数内不收敛!);disp(最大迭代次数后的结果为);x Gauss-Seidel迭代法 function [n,x]=gaussseidel(A,b,X,nm,w) %用高斯-赛德尔迭代法求解方程组Ax=b %输入：A为方程组的系数矩阵，b为方程组右端的列向量，X为迭代初值构成的列向量，nm为最大迭代次数，w为误差精度 %输出：x为求得的方程组的解构成的列向量，n为迭代次数n=1;m=length(A);I=eye(m); %生成m*m阶的单位矩阵D=diag(diag(A)); %令A=D-L-U,计算矩阵DL=tril(-A)+D; %令A=D-L-U,计算矩阵LU=triu(-A)+D; %令A=D-L-U,计算矩阵UM=inv(D-L)*U; %计算迭代矩阵g=inv(I-inv(D)*L)*(inv(D)*b); %计算迭代格式中的常数项%下面是迭代过程while n=nm x=M*X+g; %用迭代格式进行迭代 if norm(x-X,2)w disp(迭代次数为);n disp(方程组的解为);x return; %上面：达到精度要求就结束程序，输出迭代次数和方程组的解 end X=x;n=n+1;end %下面：如果达到最大迭代次数仍不收敛，输出警告语句及迭代的最终结果（并不是方程组的解）disp(在最大迭代次数内不收敛!);disp(最大迭代次数后的结果为);x 系数矩阵为H时Gauss-seidel迭代求解 H=hilb(3) H = 1.0000 0.5000 0.3333 0.5000 0.3333 0.2500 0.3333 0.2500 0.2000 X=[1,1,1]; b=[-12,20,3]; w=1.0e-4; nm=1000; gaussseidel(H,b,X,nm,w) 迭代次数为 n = 715 方程组的解为 x = 1.0e+003 * -0.7380 3.7320 -3.4200 ans = 715 求高斯迭代矩阵及普半径 D=diag(diag(H)); %令A=D-L-U,计算矩阵D L=tril(-H)+D; %令A=D-L-U,计算矩阵L U=triu(-H)+D; %令A=D-L-U,计算矩阵U M=inv(D-L)*U; %计算迭代矩阵 M M = 0 -0.5000 -0.3333 0 0.7500 -0.2500 0 -0.1042 0.8681 [x,C]=eig(M) x = 1.0000 0.6812 -0.1425 0