【数值分析】超松弛迭代（SOR）

---

1.数值计算方法

Young, Jr. 和Stanley P. Frankel于1950年提出，主要为了实现在计算机上求解线性方程组。逐次超松弛迭代法从高斯-塞德尔迭代出发，加入了$w$作为松弛因子，加快了迭代的收敛速度。

2.算法实现

% 超松弛迭代（SOR）
function x = SOR(A, b, x0, maxIter, tol, w)
%input
%A 系数矩阵
%b 非齐次项
%x0 初始解
%maxIter 最大迭代次数
%tol 误差精度
%w 松弛因子  0<w<1，低松弛法；1<w<2，超松弛法
x = x0;
n = length(b);
for t = 1:maxIter
    for i = 1:n
        x(i) = (1 - w) * x(i) +  (b(i) - A(i, 1:i - 1) * x0(1:i - 1) - A(i, i+1:n) * x0(i+1:n)) * w / A(i,i);
    end
    
    if max(abs(x - x0)) < tol
        break;
    end
    x0 = x;
end