jacobi.m

来自「这是机械工业张德丰《MATLAB与外部程序接口编程》随书源码」· M 代码 · 共 16 行

%jacobi.m
function s=jacobi(A,b,X0,epx)
%A为系数矩阵，b即为Ax=b中的b,X0为迭代初值，eps为控制精度
D=diag(diag(A));  %求出对角阵
D=inv(D);  %求出对较阵的逆矩阵
L=tril(A,-1);  %求出严格下三角矩阵
U=triu(A,1);  %示出严格上三角矩阵
B=-D*(L+U);
f=D*b;
s=B*X0+f;
while norm(s-X0)>=eps
    X0=s;
    s=B*X0+f;  %迭代的矩阵表达形式
end
return