matlab雅可比迭代法函数

以下是使用Matlab实现雅可比迭代法的函数示例：

function [x, n] = jacobi(A, b, x0, eps)
%定义Jacobi迭代函数
%输入参数：
%A：系数矩阵
%b：常数向量
%x0：初始解向量
%eps：误差限
%输出参数：
%x：方程组的数值解
%n：迭代次数

%计算迭代矩阵
D = diag(diag(A));
L = -tril(A,-1);
U = -triu(A,1);
BJ = D\(L+U);
f = D\b;

%判断收敛性
a = max(abs(eig(BJ)));
if a >= 1
    disp('Jacobi迭代不收敛');
    return %不再向下执行
else
    n = 1;
    x = BJ*x0 + f;
    while norm(x-x0,inf)>=eps %无穷范数
        x0 = x;
        x = BJ*x0+f;
        n = n+1;
    end
end

该函数的输入参数包括系数矩阵A、常数向量b、初始解向量x0和误差限eps，输出参数包括方程组的数值解x和迭代次数n。函数中首先计算迭代矩阵BJ和常数向量f，然后判断迭代矩阵的收敛性，如果不收敛则输出提示信息并结束函数，否则进行迭代计算直到满足误差限为止。