jacobi_function.m

来自「该函数文件是通过雅克比迭代格式求解线性方程的近似解」· M 代码 · 共 19 行

function J=Jacobi(A,b,x0,eps)
%该函数文件是通过雅克比迭代格式求解线性方程的近似解
%其中，A为方程组的系数矩阵，b为方程组的右端向量，x0为迭代初始值（向量）
if nargin==3
    eps=1.0e-6;
else error
    return
end
D=diag(diag(A));
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