已知方阵A=[1 2 0 1,0 2 1 0,2 0 3 1,2 0 3 K] 的行最简形矩阵为H （K为学生自己的序号，注意：不是学号），求可逆矩阵P 使得PA=H 。写一段matlab代码来解决这个问题

在MATLAB中，我们可以使用`rref`函数来计算矩阵A的行最简形矩阵H，然后通过Householder反射或初等行变换找到可逆矩阵P。这里假设矩阵A的秩等于其阶数（即4），因为有非零元素构成了一组线性无关的行向量，这通常意味着A是满秩的。

% 假设矩阵A存储在一个变量中（例如A = ...）
[n, m] = size(A); % 获取矩阵的维度

% 使用rref函数计算行最简形矩阵H
[H, P] = rref(A); % H为行最简形矩阵，P为对应的行操作矩阵

% 初等行变换可能不会直接得到单位矩阵形式，所以我们需要检查P是否已经是可逆的
if rank(P) == n && isequal(diag(diag(P)), ones(n, 1)) % 如果P是对角阵且对角线都是1，它就是可逆的
    disp('矩阵P已经是可逆的');
else
    error('矩阵P不可逆，无法找到满足条件的P');
end

disp(['可逆矩阵P为:\n', num2str(P)]); % 输出P

运行上述代码后，如果能得到可逆矩阵P，那么它将满足PA=H。若代码返回错误，说明可能存在一些特殊情况，比如K导致的秩不足，此时可能需要处理或者调整输入矩阵A。