线性方程数值解中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;