Least-Squares Joint Diagonalization of a matrix set by a congruence transformation

Marco Congedo (1) , Dinh-Tuan Pham (2)

1 GIPSA-SIGMAPHY - GIPSA - Signal Images Physique
2 SMS - Statistique et Modélisation Stochatisque

The approximate joint diagonalization (AJD) is an important analytic tool at the base of numerous independent component analysis (ICA) and other blind source separation (BSS) methods, thus finding more and more applications in medical imaging analysis. In this work we present a new AJD algorithm named SDIAG (Spheric Diagonalization). It imposes no constraint either on the input matrices or on the joint diagonalizer to be estimated, thus it is very general. Whereas it is well grounded on the classical leastsquares criterion, a new normalization reveals a very simple form of the solution matrix. Numerical simulations shown that the algorithm, named SDIAG (spheric diagonalization), behaves well as compared to state-of-the art AJD algorithms.

Mots clés: Approximate Joint Diagonalization

Domaines: Calcul [stat.CO]; Statistiques [math.ST]; Théorie [stat.TH]

Marco Congedo, Dinh-Tuan Pham. Least-Squares Joint Diagonalization of a matrix set by a congruence transformation. SinFra'09 - 2nd Singaporean-French IPAL Symposium, Feb 2009, Singapour, Singapore. pp.96-106, ⟨10.1142/9789814277563_0010⟩. ⟨hal-00371934⟩