scaledpivoting.m

本文件提供了几个数值分析中的Matlab小程序。水平有限

function X=scaledpivoting(A,b)
[n,n]=size(A);
X=zeros(n,1);
Aug=[A,b];
for i=1:n
    S(i)=max(abs(Aug(i,1:n)));
end
for k=1:n-1
    for i=k+1:n
        if abs(Aug(i,k)/S(i))>abs(Aug(k,k)/S(k))
            t=Aug(i,k:n+1);
            Aug(i,k:n+1)=Aug(k,k:n+1);
            Aug(k,k:n+1)=t;
        end
    end
    for i=k+1:n
        m=Aug(i,k)/Aug(k,k);
        Aug(i,k:n+1)=Aug(i,k:n+1)-m*Aug(k,k:n+1);
    end
end
X=backsub(Aug(1:n,1:n),Aug(1:n,n+1));