matlab.zip_GS法迭代MATLAB_iterative_迭代

function [x,k,e ] = sor( A,b,y,w ) %y是精确解 %sor迭代 % Detailed explanation goes here N=100000000; tol = 1e-4; q=1-1/(100*100); D = diag(diag(A)); L = triu(A)-A; U = tril(A)-A; % 生成矩阵B, 向量g B = (D-w*L)\((1-w)*D+w*U); g = (D-w*L)\b*w; %迭代首项 n = length(b); x0=eye(n,1); x1=eye(n,1); for i=1:n x0(i)=1; x1(i)=0; end r=1; k=0; while r > tol * (1-q)/q x1=B*x0+g; % 迭代公式 r = norm(x1-x0, inf); % 使用无穷范数 x0 = x1; k = k+1; if k>=N % 如果达到最大迭代次数 disp( '达到最大迭代次数，迭代失败'); break ; end end x = x0; e = norm(y-x, inf); end