srref.m

Matlab中虽然有很多解方程的函数

%用来求简化行阶梯矩阵和他的列的变换向量
function [S,t]=srref(a)
%a是一个矩阵m*n
%S是简化行阶梯型的且前面是Ｅ的形式
[m,n]=size(a);
A=rref(a);
r=rank(A);
A=A(1:r,:);
t=1:n-1;
for k=2:r
    for p=k:n-1
        if A(k,p)==1
            temp=A(:,k);A(:,k)=A(:,p);A(:,p)=temp;
            tt=t(k);t(k)=t(p);t(p)=tt;
        end
    end
end
S=[A;zeros(m-r,n)];