Matlab实现赛德尔(Seidel)迭代

 

A=[1,0.4,0.4;0.4,1,0.8;0.4,0.8,1];
b=[1;2;3];
n=length(b);
N=100;
ep=1e-6;
x0=zeros(n,1);
x=zeros(n,1);
k=0;
Result=[];
while k<N
    for i=1:n
        if i==1
            x(1)=(b(1)-A(1,2:n)*x0(2:n))/A(1,1);
        else if i==n
                x(n)=(b(n)-A(n,1:n-1)*x(1:n-1))/A(n,n);
            else
                x(i)=(b(i)-A(i,1:i-1)*x(1:i-1)-A(i,i+1:n)*x0(i+1:n))/A(i,i);
            end
        end
    end
    if norm(x-x0,inf)<ep
        break;
    end
%将每次的结果并入一个矩阵
    Result=[Result;x'];
    x0=x;
    k=k+1;
end
%将数据存表格中
xlswrite('b.xlsx',Result);
disp(['k=',num2str(k)])