An Alternating Projection Algorithm for Approximate Simultaneous Diagonalization

In this paper, we present a novel formulation of Approximate Simultaneous Diagonalization (ASD) with a nonconvex feasibility problem to find a structured low rank matrix whose building blocks are the Kronecker sums of given multiple matrices. To tackle this feasibility problem, we propose an alternating projection algorithm that can generate a matrix sequence approaching monotonically to a solution. By this algorithm, simultaneously diagonalizable matrices are obtained in the neighborhood of the given matrices which are not necessarily diagonalizable simultaneously. By using further the Diagonalize-One-then-Diagonalize-the-Other (DODO) method, we can obtain finally a common similarity transformation which diagonalizes the simultaneously diagonalizable matrices. Numerical experiments show that, compared with a Jacobi-like method, the proposed algorithm achieves a better approximation to the desired common similarity transformation.

[1]  A. Bunse-Gerstner,et al.  Numerical Methods for Simultaneous Diagonalization , 1993, SIAM J. Matrix Anal. Appl..

[2]  Josef A. Nossek,et al.  Simultaneous Schur decomposition of several nonsymmetric matrices to achieve automatic pairing in multidimensional harmonic retrieval problems , 1998, IEEE Trans. Signal Process..

[3]  Joos Vandewalle,et al.  Computation of the Canonical Decomposition by Means of a Simultaneous Generalized Schur Decomposition , 2005, SIAM J. Matrix Anal. Appl..

[4]  Florian Roemer,et al.  A semi-algebraic framework for approximate CP decompositions via simultaneous matrix diagonalizations (SECSI) , 2013, Signal Process..

[5]  R. Plemmons,et al.  Structured low rank approximation , 2003 .

[6]  X. Luciani,et al.  Canonical Polyadic Decomposition based on joint eigenvalue decomposition , 2014 .

[7]  Ed F. Deprettere,et al.  Azimuth and elevation computation in high resolution DOA estimation , 1992, IEEE Trans. Signal Process..

[8]  Aude Rondepierre,et al.  On Local Convergence of the Method of Alternating Projections , 2013, Foundations of Computational Mathematics.

[9]  øöö Blockinø Phase retrieval, error reduction algorithm, and Fienup variants: A view from convex optimization , 2002 .

[10]  Laurent Condat,et al.  Cadzow Denoising Upgraded: A New Projection Method for the Recovery of Dirac Pulses from Noisy Linear Measurements , 2015 .

[11]  Isao Yamada,et al.  An Efficient Greedy Algorithm for finding the Nearest Simultaneous Diagonalizable Family , 2018, 2018 IEEE Statistical Signal Processing Workshop (SSP).

[12]  Adrian S. Lewis,et al.  Local Linear Convergence for Alternating and Averaged Nonconvex Projections , 2009, Found. Comput. Math..

[13]  Ed F. Deprettere,et al.  Analysis of joint angle-frequency estimation using ESPRIT , 2003, IEEE Trans. Signal Process..

[14]  P. Comon,et al.  ICAR: a tool for blind source separation using fourth-order statistics only , 2005, IEEE Transactions on Signal Processing.

[15]  Antoine Souloumiac,et al.  Jacobi Angles for Simultaneous Diagonalization , 1996, SIAM J. Matrix Anal. Appl..

[16]  Laurent Albera,et al.  Joint Eigenvalue Decomposition of Non-Defective Matrices Based on the LU Factorization With Application to ICA , 2015, IEEE Transactions on Signal Processing.