还原行梯形形式的转换矩阵

Question

i可以使用命令R = rref(C)计算MATLAB中矩阵C的简化行echelon表格R。

但是，我还想跟踪执行的步骤，即获得给我TC = R的转换矩阵T。据我所知，此矩阵应该在使用高斯时被隐式计算 - 犹太人消除。

有没有办法获得T？也许是解决方法？在MATLAB文档中，我找不到任何信息。其他编程语言中是否有RREF函数会返回t？

Answer 1

您可以使用一个事实，即基本行操作等效于繁殖左侧的基本矩阵。令c为大小的矩阵(MXN(;

z= rref([c eye(m)]); % [c I] is multiplied by some matrix T
% the result is [rref(c) T]
r= z(:,1:n); % the reduced row echelon form of c
t= z(:,n+1:end); % now we have T