Рассчитать полярную декомпозицию матовой матрицы?

Question

Полярное разложение квадратной комплексной матрицы A является матричным разложением формы

A = UP,

где U - унитарная матрица, а P - положительно-полуопределенная эрмитова матрица. Это отличается от QR-разложения , реализованного как scipy.linalg.qr .

Как я могу найти P с numpy / scipy?

Answer 1

Сейчас есть "scipy.linalg.polar" :) http://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.linalg.polar.html

Answer 2

Ну, страница в Википедии о полярной декомпозиции , на которую вы ссылались, содержит формулы для ее вычисления из SVD.

Во-вторых, в нескольких абзацах от того, на который вы ссылались из учебника Scipy Linalg, объясняется, как вычислить SVD .

Итак, объединяя эти два, вы должны получить полярное разложение, нет?