SOR.rar_高阶矩

function [x,k,flag]=SOR(A,b,ep,w,max1) %SOR迭代法 %A为方程组的系数矩阵 %b为方程组的右端项 %delta为精度要求，缺省值为1e-5 %max1为最大迭代次数，缺省值为100 %x为方程组的解 %k为迭代次数 %w为超松弛因子，缺省值为1 %flag为指标变量 flag='OK'表示迭代收敛到指标要求 % flag='fail'表示迭代失败 if nargin<5 max1=100; end if nargin<4 w=1; end if nargin<3 ep=1e-5; end n=length(A); k=0; x=zeros(n,1); y=zeros(n,1); flag='OK'; while 1 y=x; for i=1:n z=b(i); for j=1:n if j~=i z=z-A(i,j)*x(j); end end if abs(A(i,i))<1e-10|k==max1 flag='fail'; return; end z=z/A(i,i); x(i)=(1-w)*x(i)+w*z; end if norm(y-x,inf)<ep break; end k=k+1; end