Matlab SOR方法实现

Question

使用零向量的初始近似值并且不考虑容差，我将代码缩短为只需要 4 个参数。使得 x1 总是等于 c，以此类推 x(k+1)=x(k)T+c。

但是，代码似乎没有产生您期望的正确近似值。有没有人注意到我在哪里搞砸了？假设 DLU_decomposition(A) 返回正确的矩阵。

function x = sor2(A,b,omega,kmax)
[D,L,U] = DLU_decomposition(A);
T=inv(D-omega*L)*(((1-omega)*D)+(omega*U));
c= (omega*inv(D-omega*L))*b;

for k=1:kmax,

    if(k==1),
        x=c;
    end
    x=T*x+c;

end
norm(A*x-b)
end

Answer 1

好吧，我猜所有的困惑都可能来自乘法。您需要按元素计算矩阵 --> 使用 .* 而不是普通的 *。这会提供正确的近似值吗？