relax.m

SOR方法解线性方程组的matlab程序例子（李庆扬数值分析例子）

Nmax=1000;
eps=5e-6;
w=0.5;
x0=ones(3,1);

a=[4 -1 0;-1 4 -1;0 -1 4];
b=[1 4 -3]';
D=diag(diag(a));
L=-tril(a,-1);
U=-triu(a,1);
C=inv(D-w*L);
B=C*((1-w)*D+w*U);
f=w*C*b;
x1=B*x0+f;
k=1;
while norm(x1-x0)>=eps
x0=x1;
x1=B*x0+f;
k=k+1;
  if k>Nmax
   ps='迭代失败'
   break;
  end
end

ss=x1