超松弛迭代法matlab代码

以下是超松弛迭代法的 MATLAB 代码示例：

function [x, k] = sor(A, b, x, omega, tol, maxit)
% 超松弛迭代法求解线性方程组 Ax=b
% 输入参数：
% A: 系数矩阵
% b: 右端向量
% x: 初始向量
% omega: 松弛因子
% tol: 相对误差容限
% maxit: 最大迭代次数
% 输出参数：
% x: 方程组的解
% k: 实际迭代次数

n = length(b);
x = x;
k = ;
while k < maxit
x_old = x;
for i = 1:n
sigma = ;
for j = 1:n
if j ~= i
sigma = sigma + A(i, j) * x(j);
end
end
x(i) = (1 - omega) * x(i) + omega * (b(i) - sigma) / A(i, i);
end
k = k + 1;
if norm(x - x_old) / norm(x_old) < tol
break;
end
end

if k == maxit
warning('超过最大迭代次数，可能未收敛！');
end

end