jacobi.m

来自「利用高斯消元法、雅克比消元法解决求解线性方程组的解的程序」· M 代码 · 共 22 行

%%本程序使用Jacobi迭代法解线性代数方程组
%%jacobi.m
function x=jacobi(A,b)
err = 0.0001;

[n m] = size(A);
D = diag(diag(A));
L = (tril(A)-D)*(-1);
U = (triu(A)-D)*(-1);
x(:,1) = zeros(n,1);    %%产生初始向量
invD = diag(1.0./diag(A));
k = 1;
while 1
   x(:,k+1)=invD*L*x(:,k)+invD*U*x(:,k)+invD*b;
   if max(abs(x(:,k+1)-x(:,k)))>err
      k = k+1;
   else
      break
   end
end