richardson迭代法求解线性方程组的matlab代码

以下是使用Richardson迭代法求解线性方程组Ax=b的Matlab代码：

function [x, err] = richardson(A, b, x0, w, tol, maxit)

% A: 系数矩阵
% b: 常数向量
% x0: 初始解
% w: 松弛因子
% tol: 相对误差容限
% maxit: 最大迭代次数

n = length(b);
x = x0;
err = zeros(maxit, 1);

for k = 1:maxit
    x_old = x;
    x = x + w * (b - A * x);
    err(k) = norm(x - x_old) / norm(x);
    if err(k) < tol
        break;
    end
end

err = err(1:k);

其中，x表示方程的解，err表示每次迭代后的相对误差。在实际应用中，需要根据具体问题设定系数矩阵A、常数向量b、初始解x0、松弛因子w、相对误差容限tol和最大迭代次数maxit等参数。