第六章sor.txt

来自「包括牛顿法」· 文本 代码 · 共 26 行

function y=SOR(a,b,w,x0)
D=diag(diag(a));
U=-triu(a,1);
L=-tril(a,-1);
B=(D-w*L)\((1-w)*D+w*U);
g=(D-w*L)\b*w;
y=B*x0+g;
n=1;
while norm(y-x0)>=1.0e-6
    x0=y;
    y=B*x0+g;
    n=n+1;
end
y
n

a=[2 -1 0;-1 3 -1;0 -1 2];
b=[1;8;-5];
SOR(a,b,1.1,[0;0;0])

y =
    2.0000
    3.0000
   -1.0000
n =
    10