Matlab-1:jacobi迭代法工具箱

function [u,n]=Jacobi(A,b,u0,eps,varargin)
%Jacobi.m函数为用于雅可比迭代法求解线性方程组
%A为线性方程组的系数矩阵
%b为线性方程组的常数向量
%u0为迭代初始向量
%eps为解的精度控制
%varargin为迭代步数控制
%u为线性方程组的解
%n为求出所有精度的解实际的迭代步数
if nargin==3
    eps=1.0e-10;
    M=200;
elseif nargin<3
    error
    return
elseif nargin==5
    M=varargin{1};
end
D=diag(diag(A));%求A的对角矩阵
L=-tril(A,-1);%求A的下三角阵
U=-triu(A,1);%求A的上三角阵
B=D(L+U);
f=D;
u=B*u0+f;
n=1;  %迭代次数
while norm(u-u0)>=eps & n<=1000
    u0=u;
    u=B*u0+f;
    n=n+1;
%     if (n>=M)
%         disp('Warning:迭代次数太多，可能不收敛！');
%         return;
%     end
end

 算法推导：