线性方程数值解中jacobi迭代matlab实现

%jacobi迭代
%author:野渡无人
%create date:2007-11-23
function res=jacobi(a,b,x0,eps)
%a为系数矩阵,b为方程组常数向量，x0为初值
if nargin==3
    eps=1.0e-6;
elseif nargin<3
    error
    return;
end
D=diag(diag(a));%求矩阵对角阵
D=inv(D);
L=tril(a,-1);%求严格下三角矩阵
U=triu(a,1);%求严格上三角矩阵
B=-D*(L+U);
if norm(B,inf)>=1
    fprintf('B范数>=1.');
    return;
end
f=D*b;
res=B*x0+f;
i=0;
while norm(res-x0)>=eps
    x0=res;
    res=B*x0+f;
    i=i+1;
end