songchidiedai.m

来自「线性方程组的迭代解法解法」· M 代码 · 共 27 行

function [x,n]=Songchidiedai(A,b,w,eps)
% w为松弛迭代因子，目前对于一般的矩阵A
% 还没有最好的w，只是对于一些特殊的矩阵才有
% w 为0～2
if nargin==3
    eps=1e-5;
end
D=diag(diag(A));
L=-tril(A,-1);
U=-triu(A,1);
Sw=inv(D-w*L)*((1-w)*D+w*U);
f=w*inv(D-w*L)*b;
if max(abs(eig(Sw)))<1
    x0=zeros(length(A),1);
    x1=Sw*x0+f;
    k=0;
    while norm(x1-x0,inf)>eps
        x0=x1;
        x1=Sw*x0+f;
        k=k+1;
    end
    x=x1;
    n=k;
else
    disp('不满足jacobi迭代条件！');
end
end