Matlab-1:jacobi迭代法工具箱

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

算法推导：