三种迭代速度比较

1、5.三种迭代速度比较程序：function v,sN,vChain=gaussseidel(A,b,x0,errorBound,maxSp)%Gauss-Seidel迭代法求解线性方程组%A-系数矩阵b-右端向量x0-初始迭代点errorBound-近似精度maxSp-最大迭代次数%v-近似解sN-迭代次数vChain-迭代过程的所有值step=0;error=inf;s=size(A);D=zeros(s(1);vChain=zeros(15,3); %最多能记录15次迭代次数k=1;fx0=x0;for i=1:s(1)D(i,i)=A(i,i);end;L=-tril(A,-1);U=-

2、triu(A,1);while error=errorBound & stepmaxSpx0=inv(D)*(L+U)*x0+inv(D)*b;vChain(k,:)=x0;k=k+1;error=norm(x0-fx0);fx0=x0;step=step+1;endv=x0;sN=step;function x, k, index=jacobi(A, b, ep, it_max)%求解线性方程组的Jacobi迭代法，其中% A -方程组的系数矩阵% b -方程组的右端项% ep -精度要求。省缺为1e-5% it_max最大迭代次数，省缺为100% x -方程组的解% k -迭代次数% in

3、dexindex=1表示迭代收敛到指定要求；% index=0表示迭代失败if nargin 4 it_max=100; endif nargin 3 ep=1e-5; end n=length(A); k=0;x=zeros(n,1); y=zeros(n,1); index=1;while 1for i=1:ny(i)=b(i);for j=1:nif j=iy(i)=y(i)-A(i,j)*x(j);endendif abs(A(i,i)1e-10 | k=it_max index=0 return;endy(i)=y(i)/A(i,i);endif norm(y-x,inf)epbre

4、ak;endx=y; k=k+1;Endfunction x,k,index=SOR(A,b,ep,w,it_max)%求解线性方程组的SOR迭代法% A为系数矩阵% b为方程组右端项% ep为精度要求，缺省值1e-5% w为超松弛因子，缺省值为1；% it_max为最大迭代次数，缺省值100% x为方程组的解% k为迭代次数% index为指标变量 index=1表示迭代收敛到指定要求 index=0表示迭代失败。 if nargin5 it_max=100;end if nargin4 w=1;endif nargin3 ep=1e-5;endn=length(A);k=0;x=zeros

5、(n,1);y=zeros(n,1);index=1;while 1y=x;for i=1:nz=b(i);for j=1:nif j=iz=z-A(i,j)*x(j);endendif abs(A(i,i)1e-10|k=it_maxindex=return;endz=z/A(i,i);x(i)=(1-w)*x(i)+w*z;endif norm(y-x,inf) A=4 3 0;3 4 -1;-1 0 4;b=24;30;-24; x,k,index=jacobi(A,b,5e-8,100)x =3.09683.8710-5.2258k =80index =1 A=4 3 0;3 4 -1

6、;-1 0 4;b=24;30;-24;x0=1;1;1; v,sN,vChain=gaussseidel(A,b,x0,5e-008,11)v =3.31584.1326-5.3058sN =11vChain =5.25007.0000-5.75000.75002.1250-4.68754.40635.7656-5.81251.67582.7422-4.89843.94345.0186-5.58112.23613.1472-5.01423.63964.5694-5.44102.57303.4101-5.09013.44254.2978-5.35682.77673.5790-5.13943.31584.1326-5.3058000000000000 A=4 3 0;3 4 -1