Guass---Seidel迭代算法的matlab代码

以下是使用matlab实现Guass-Seidel迭代算法的示例代码：

function [x, iter] = gauss_seidel(A, b, x0, tol, max_iter)
% A: 系数矩阵
% b: 常数向量
% x0: 初始解向量
% tol: 容差
% max_iter: 最大迭代次数
% x: 迭代结果
% iter: 实际迭代次数

[n, ~] = size(A);

% 初始化
x = x0;
iter = 0;
err = tol + 1;

while err > tol && iter < max_iter
    for i = 1 : n
        % 更新x(i)
        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(A * x - b);
    iter = iter + 1;
end

end

使用示例：

A = [4 -1 0; -1 4 -1; 0 -1 4];
b = [10; 10; 10];
x0 = [0; 0; 0];
tol = 1e-6;
max_iter = 1000;

[x, iter] = gauss_seidel(A, b, x0, tol, max_iter);
disp(['迭代结果：', num2str(x')]);
disp(['实际迭代次数：', num2str(iter)]);

注意：这只是一个示例代码，Guass-Seidel迭代算法可能不适合所有的线性方程组，需要根据具体情况选择算法。