1.数值计算方法
Young, Jr. 和Stanley P. Frankel于1950年提出,主要为了实现在计算机上求解线性方程组。逐次超松弛迭代法从高斯-塞德尔迭代出发,加入了
w
w
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