partialpivoting.m

来自「本文件提供了几个数值分析中的Matlab小程序。水平有限」· M 代码 · 共 19 行

function X=partialpivoting(A,b)
[n,n]=size(A);
X=zeros(n,1);
Aug=[A,b];
for k=1:n-1
    for i=k+1:n
        if abs(Aug(i,k))>abs(Aug(k,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))