Факторизация Шура для нахождения унитарной матрицы?

#r #matlab #matrix

#r #matlab #матрица

Вопрос:

После выполнения факторизации по Шуру с использованием функции Schur в библиотеке ‘Matrix’, как мне найти соответствующую унитарную матрицу в R?

Я могу сделать это в MATLAB с помощью функции Schur (документация здесь), однако пакет Schur в R предоставляет только треугольную в факторизации Шура = $$ U * T * U ‘ $$

Комментарии:

1. В классе Schur, похоже, есть слот для Q «Квадратной ортогональной «матрицы»» (из документов). Таким образом, я думаю, вы бы просто сделали Sch.A <- Schur(A); U <- Sch.A@Q .

2. Спасибо! Рад ответить на это, если вы хотите скопировать / вставить свой комментарий

Ответ №1:

Просматривая документы для пакета «Matrix», я заметил, что в классе ‘Schur’ есть слот для Q , который является ‘квадратной ортогональной «матрицей»‘, связанной с разложением. Итак, вы хотите сделать:

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

Это слегка сбивает с толку, потому что они цитируют разложение как $ A = Q ^ { top} T Q $ , возможно, поэтому вы это пропустили.