jacobi.m

有关matlab r2007的实用教程

function J=Jacobi(A,b,x0,eps)
%该函数文件是通过Jacobi迭代格式求解线性方程组的近似解
%其中，A为方程组的系数矩阵，b为方程组的右端向量，x0为选代初值
if nargin==3
    eps=1.0e-6;
elseif nargin<3
    error
    return
end
D=diag(diag(A));    %求解对角阵D
D=inv(D);         %求解对角矩阵的逆阵
L=tril(A,-1);       %求解出严格下三角阵
U=triu(A,1);       %求解出严格上三角阵
B=-D*(L+U);
f=D*b;
J=B*x0+f;
while norm(J-x0)>=eps
    x0=J;
    J=B*x0+f;
end