matlab jacobi迭代法代码

好的，以下是 MATLAB 中 Jacobi 迭代法的代码：

function [x, iter] = jacobi(A, b, x0, tol, maxiter)
% Jacobi 迭代法求解线性方程组 Ax=b
% 输入参数：
% A：系数矩阵
% b：常数向量
% x0：初始解向量
% tol：容差
% maxiter：最大迭代次数
% 输出参数：
% x：方程组的解向量
% iter：实际迭代次数

n = length(b);
x = x0;
iter = 0;
while iter < maxiter
    x_old = x;
    for i = 1:n
        sigma = 0;
        for j = 1:n
            if j ~= i
                sigma = sigma + A(i,j) * x_old(j);
            end
        end
        x(i) = (b(i) - sigma) / A(i,i);
    end
    if norm(x - x_old) < tol
        break;
    end
    iter = iter + 1;
end

if iter == maxiter
    warning('Jacobi 迭代法未收敛！');
end

end

其中，`A` 是系数矩阵，`b` 是常数向量，`x0` 是初始解向量，`tol` 是容差，`maxiter` 是最大迭代次数。函数的输出参数为 `x` 和 `iter`，分别表示方程组的解向量和实际迭代次数。