2011-11-11 14 views
4

R中是否有与MATLAB函数ordschur(文档here)等价的函数?R等效于Schur分解的Matlab重新排序?

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

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

回答

2

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