R 相当于 Matlab 对 Schur 分解的重新排序？

Question

在 R 中是否有与 MATLAB 函数 ordschur(文档 here)等价的函数？

该函数对 schur 函数生成的 Schur 分解 X = U*T*U' 重新排序，并返回重新排序的 Schur 矩阵 TS 和累积正交变换 US，使得 X = US*TS*US'。我对“lhp”方法特别感兴趣 - 在 MATLAB 文档链接中也有描述。

请注意，矩阵包中有一个函数 Schur(请参阅 CRAN 文档 here)R，它执行方矩阵的 Schur 分解和特征值。更新:该函数还返回酉正交矩阵 U。

Answer 1

据我所知，MATLAB 使用 LAPACK 中的 ?TRSEN 函数来执行重新排序。你可以看看一些有限的实现here .为了将此功能引入 R，您可以自己实现此例程。