r - ユニタリ行列を見つけるためのシューア因数分解?

Question

ライブラリ 'Matrix' の関数 Schur を使用して Schur 分解を実行した後、R で関連するユニタリ行列を見つけるにはどうすればよいですか?

MATLAB で関数 Schur (ドキュメントはこちら) を使用してこれを行うことができますが、R の Schur パッケージは Schur 分解の三角形のみを提供します = $$ U * T * U' $$

Answer 1

「Matrix」パッケージのドキュメントを見て、「Schur」クラスに、分解に関連付けられたQ「正方直交「行列」」である のスロットがあることに気付きました。だからあなたがしたい：

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

これは、分解を $A = Q^{\top} TQ$ として引用しているため、やや混乱しています。