Polar decomposition for quaternionic matrices? A non-zero complex number can be uniquely written in polar form as $re^{i\theta}$.  There is an analogous result for complex matrices: any invertible complex matrix can be uniquely written as $UP$, where $U$ is a unitary matrix and $P$ is a positive definite Hermitian matrix (see e.g. the description on Wikipedia)
Suppose we consider $n\times n$ invertible matrices over the quaternions.  Is there an analogous polar decomposition?
Or to ask the question a little more colorfully, can we complete the following list?


*

*roots of unity: positive reals

*unitary matrices: positive definite matrices

*compact symplectic matrices: ???


PS: Note that there exists a polar decomposition for quaternions; I'm interested in the result for matrices of quaternions.
 A: For all matrices, even non-invertible matrices, one gets a polar decomposition  $X=UP$  with  $U$  unitary and $P$  positive semidefinite.  Of course $U$ is not unique unless $X$ is invertible.  I found this result in "Quaternions and matrices of quaternions" by F. Zhang, Linear Algebra Appl., 251 (1997), pp. 21–57.
There also the Jordan canonical form, Schur factorization and the spectral theorem.  I have a survey on this that needs some polishing:  http://arxiv.org/abs/1107.0500 "Factorization of Matrices of Quaternions."
A: A key wordphrase here is "Cartan decomposition". Since $G=SL_n(\mathbb H)$ is a semisimple group, there is a diffeomorphism 
$$K\times \mathfrak p\rightarrow G$$ 
taking $(k,X)$ to $k\cdot {\rm exp}(X)$, where $K$ is a maximal compact subgroup of $G$ (i.e. the compact symplectic group) and $\mathfrak p$ is the vector subspace of ${\mathfrak sl}_n(\mathbb H)$ fixed by the Cartan involution $X\rightarrow \bar X^t$ (the quaternionic version of Hermitian). 
A: The good news is that this is addressed in:
Zhuang, Wa Jin(PRC-ZNU)
Polar decomposition and GL partial ordering for quaternion rectangular matrices. (Chinese. English, Chinese summary) 
Adv. Math. (China) 34 (2005), no. 2, 187–193.
The bad news is that I can't seem to find the paper (but the math review is reasonably informative).
A: I guess just as for normal SVD, you can recover the polar decomposition for quaternion matrices from their corresponding SVD, for which two immediate papers that I found were:


*

*Singular value decomposition of quaternion matrices

*Quaternion SVD computation
