Schur分解找到酉矩陣？

Question

在庫'Matrix'中使用函數Schur執行Schur分解后，如何在R中找到關聯的酉矩陣？

我可以在MATLAB中用函數Schur（ 這里有文檔）來做這個，但是，R中的Schur包只提供了Schur分解中的三角形= $$ U * T * U'$$

Answer 1

看一下“矩陣”包的文檔，我注意到'Schur'類有一個Q的插槽，它是與分解相關的'Square orthogonal “Matrix”。 所以你想做：

Sch.A <- Schur(A);
U <- Sch.A@Q

這有點令人困惑，因為他們將分解引用為$ A = Q ^ {\\ top} TQ $，這也許就是你錯過它的原因。