用matlab写一个代码求得矩阵的所有行置换矩阵

以下是用 MATLAB 编写的代码，可以求得任意大小矩阵的所有行置换矩阵：

function P = row_permutations(n)
% Generate all row permutations of an n x n matrix
% Returns a 3D array where P(:,:,i) is the i-th permutation matrix

P = zeros(n,n,factorial(n));
idx = 1;

% Generate all permutations recursively
permute_rows(1:n);

    function permute_rows(k)
        if k == 1
            P(:,:,idx) = eye(n);
            idx = idx + 1;
        else
            permute_rows(k-1);
            for i = 1:k-1
                swap_rows(k,i);
                permute_rows(k-1);
                swap_rows(k,i);
            end
        end
    end

    function swap_rows(i,j)
        temp_row = P(i,:,:);
        P(i,:,:) = P(j,:,:);
        P(j,:,:) = temp_row;
    end

end

使用方法：输入矩阵大小 n 后，函数将返回一个大小为 n x n x factorial(n) 的矩阵 P，其中 P(:,:,i) 表示第 i 个行置换矩阵。例如，对于一个 3 x 3 的矩阵，可以使用以下命令来计算所有行置换矩阵：

P = row_permutations(3);

然后，可以使用 P(:,:,i) 访问第 i 个行置换矩阵。