sor.asv

为配套王能超版本数值分析简明教程编写,方便学生教师演示!

function SOR(A,b,e,w)
[i,j]=size(A);
[o,q]=size(b);
I=eye(i);
[l,u,p]=lu1(A);
if w<2 & w>0
    if i==j & i==o;
        if diag(p)~=zeros(i,1);
            B=inv(p-w*l)*((1-w)*p+w*u);
            a=abs(eig(B));
            a=a';
            %n=0;
            if a<ones(1,i);
                m=1;
                x1=zeros(i,1);
                x2=B*x1+w*inv(p-w*l)*b;
                while max(abs(x1-x2))>e;
                    x1=x2;
                    x2=B*x1+w*inv(p-w*l)*b;
                    fprintf('%d',m);
                    m=m+1;
                    for k=1:i;
                        fprintf('%13.8f',x1(k));
                    end
                    fprintf('\n');
                end
            else
                fprintf('此多项式方程组不能进行SOR迭代');
            end
        else
            fprintf('此多项式方程组的矩阵的主对角线上有零元,不能进行SOR迭代');
        end
    else
        fprintf('此多项式方程组不能进行SOR迭代');
    end
else
    disp('')