R相当于舒尔分解的Matlab重新排序？

Question

R 中是否有等效于 MATLAB 函数 ordschur（文档 here）的函数？

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

请注意，Matrix 包中有一个函数 Schur（参见 CRAN 文档 here）R，它执行方阵的 Schur 分解和特征值。更新：此函数还返回酉正交矩阵 U。

Answer 1

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