高斯塞德尔迭代matlab

高斯塞德尔迭代是一种解线性方程组的迭代方法，可以用Matlab实现。下面是一个简单的高斯塞德尔迭代的Matlab代码实现：

function [x, iter] = gauss_seidel(A, b, x0, tol, max_iter)
% 高斯塞德尔迭代求解线性方程组 Ax = b
% A: 系数矩阵
% b: 右侧常数向量
% x0: 初始解向量
% tol: 相对误差容限
% max_iter: 最大迭代次数

n = length(b);
x = x0;
iter = 0;
err = tol + 1;

while err > tol && iter < max_iter
    x_old = x;
    for i = 1:n
        x(i) = (b(i) - A(i,1:i-1)*x(1:i-1) - A(i,i+1:n)*x(i+1:n))/A(i,i);
    end
    err = norm(x - x_old, inf)/norm(x_old, inf);
    iter = iter + 1;
end

其中，A为系数矩阵，b为右侧常数向量，x0为初始解向量，tol为相对误差容限，max_iter为最大迭代次数，x为求解出的解向量，iter为实际迭代次数。