Sign Up to our social questions and Answers Engine to ask questions, answer people’s questions, and connect with other people.
Login to our social questions & Answers Engine to ask questions answer people’s questions & connect with other people.
Lost your password? Please enter your email address. You will receive a link and will create a new password via email.
Please briefly explain why you feel this question should be reported.
Please briefly explain why you feel this answer should be reported.
Please briefly explain why you feel this user should be reported.
The polar decomposition of a square complex matrix A is a matrix decomposition of the form
A = UP,
where U is a unitary matrix and P is a positive-semidefinite Hermitian matrix. This differs from the QR Decomposition implemented as scipy.linalg.qr.
How can I find P with numpy/scipy?
Well, the wikipedia page on Polar decomposition you linked to contains formulas to compute it from the SVD.
Secondly, a few paragraphs up from the one you linked to from the scipy linalg tutorial, it is explained how to compute the SVD.
So, by combining these two, you should have the polar decomposition, no?