Jacobi迭代的matlab程序

Jacobi迭代是一种用于求解方程组的迭代方法，其基本思想是将矩阵分解为对角矩阵和非对角矩阵的和，然后逐次求解对角矩阵对应的线性方程组，得到近似解，不断迭代直到满足精度要求。以下是一个使用Matlab实现的Jacobi迭代程序的示例：

% Jacobi iterative method to solve linear equations
function x = jacobi(A, b, N, tol)
% Initialize solution
x = zeros(size(b));
% Iterate for N times or until desired error tolerance
for k = 1:N
x_new = zeros(size(x));
for i = 1:length(b)
% Use the formula for Jacobi iteration
x_new(i) = (b(i) - A(i,[1:i-1,i+1:end]) * x([1:i-1,i+1:end])) / A(i,i);
end
% Check for convergence
if norm(x - x_new) < tol
x = x_new;
return;
end
x = x_new;
end
% Return the best result obtained so far
return;

注意：以上内容仅供参考，请根据实际情况进行修改调整。