Я борюсь с разложением QZ двух особых матриц.
В частности, я выполнил следующую команду в Matlab
[T,S,Q,Z] = qz(V,W)
Это возвращает дватреугольные матрицы T и S, а также две унитарные матрицы Q и Z, такие что QVZ = T и QWZ = S.
Теперь я заинтересован в разбиении Q и Z на квадратный блок так, чтобы левый верхнийблок в Q, Z, S и T связан с наименьшими обобщенными собственными значениями, заданными диагональными элементами T, умноженными на обратное значение S.
Помимо того факта, что S не обратимо в моем приложенииЯ боюсь, что я также не вижу, как изменить порядок собственных значений, оставив нетронутой треугольную форму S и T.
Математическое изложение моей проблемы см. В: https://math.stackexchange.com/questions/3380997/generalized-schur-decompostion-reordering-unitary-matrices
Было бы очень полезно, если бы кто-то мог дать мне представление об общей логике фрагмента кода, который мог бы выполнить вышеупомянутое переупорядочение - учитывая матрицуes T, S, Q, Z, полученные из QZ-разложения V и W.
Большое спасибо.