permut.m

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

%置换
function B=permut(A,b,p)
%A表示需要置换的矩阵，b表示置换的方式
%p为参数，p为1时对行进行置换，p为2时对列进行置换，默认值为1
if nargin==2
    p=1;
end
if p==2
    A=A';
end
len=length(b);
[m,n]=size(A);
if m~=len
    disp('A的行数必须与b的长度相等！');
else
    [a,c]=count(b);
    if max(c)>1|a(length(a))>m
        disp('向量b有错！');
    else
        B=A(b(1),:);
        for i=2:len
            B=[B;A(b(i),:)];
        end
        if p==2
            B=B';
        end
    end
end