gaosai.asv

来自「fit programme for a serials of data」· ASV 代码 · 共 37 行

function [X,n]=gaosai(A,b,X0,e,N0)
% input A is an N*N nonsingular matrix
%       b is an N*1 matrix
%       e is the tolerance for iterations
%       X0 is an N*1 matrix,the initial guess of AX=b
% output X is an N*1 matrix; the jacobi iterations approximation to 
%        the solution of AX=b
%         N0 is the maximum number of iterations
%rou=norm(eig(A));
%if rou>=1
%   disp(sprintf('the process of iteration is diverge'));
%  return;
%end
N=length(b);
X=X0;
for n=1:N0
    error=0;
    for i=1:N
        t=X(i);
        X(i)=b(i);
        for j=1:N
            if j~=i
               X(i)=X(i)-A(i,j)*X(j); 
            end
        end
        X(i)=X(i)/A(i,i);
        n
        X
        if abs(X(i)-t)>error
           error=abs(X(i)-t); 
        end
    end
    if error<e
       return; 
    end
end
disp(sprintf('false'))