Canonical Polyadic decomposition of complex-valued multi-way arrays based on Simultaneous Schur Decomposition

In this paper, we propose a new semi-algebraic algorithm to compute the Canonical Polyadic (CP) decomposition of complex-valued multi-way arrays. The proposed algorithm is based on the Simultaneous Schur Decomposition (SSD) of particular matrices derived from the array to process. This CP algorithm solves some convergence problems of classical iterative techniques and its identifiability assumptions are less restrictive than those of other semi-algebraic methods. We also propose a new Jacobi-like algorithm to calculate the SSD of several complex-valued matrices. Finally the usefulness of the proposed method is illustrated in the context of fluorescence spectroscopy and epileptic source localization.

[1]  Rasmus Bro,et al.  MULTI-WAY ANALYSIS IN THE FOOD INDUSTRY Models, Algorithms & Applications , 1998 .

[2]  Laurent Albera,et al.  Multiway space-time-wave-vector analysis for source localization and extraction , 2010, 2010 18th European Signal Processing Conference.

[3]  J. Kruskal.  Three-way arrays: rank and uniqueness of trilinear decompositions, with application to arithmetic complexity and statistics , 1977 .

[4]  Patrick Dupont,et al.  Canonical decomposition of ictal scalp EEG reliably detects the seizure onset zone , 2007, NeuroImage.

[5]  Nikos D. Sidiropoulos,et al.  Kruskal's permutation lemma and the identification of CANDECOMP/PARAFAC and bilinear models with constant modulus constraints , 2004, IEEE Transactions on Signal Processing.

[6]  Lieven De Lathauwer,et al.  New Simultaneous Generalized Schur Decomposition methods for the computation of the Canonical Polyadic decomposition , 2010, 2010 Conference Record of the Forty Fourth Asilomar Conference on Signals, Systems and Computers.

[7]  Rasmus Bro,et al.  A comparison of algorithms for fitting the PARAFAC model , 2006, Comput. Stat. Data Anal..

[8]  Laurent Albera,et al.  Multi-way space-time-wave-vector analysis for EEG source separation , 2012, Signal Process..

[9]  Daniel M. Dunlavy,et al.  A scalable optimization approach for fitting canonical tensor decompositions , 2011 .

[10]  Laurent Albera,et al.  Semi-algebraic canonical decomposition of multi-way arrays and Joint Eigenvalue Decomposition , 2011, 2011 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[11]  Joos Vandewalle,et al.  A Multilinear Singular Value Decomposition , 2000, SIAM J. Matrix Anal. Appl..

[12]  J. Chang,et al.  Analysis of individual differences in multidimensional scaling via an n-way generalization of “Eckart-Young” decomposition , 1970 .

[13]  A. Villringer,et al.  Non-invasive optical spectroscopy and imaging of human brain function , 1997, Trends in Neurosciences.

[14]  Lieven De Lathauwer,et al.  A Link between the Canonical Decomposition in Multilinear Algebra and Simultaneous Matrix Diagonalization , 2006, SIAM J. Matrix Anal. Appl..

[15]  R. Harshman,et al.  PARAFAC: parallel factor analysis , 1994 .

[16]  R. Bro.  PARAFAC. Tutorial and applications , 1997 .

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

[18]  P. Comon,et al.  Tensor decompositions, alternating least squares and other tales , 2009 .

[19]  Ben H. Jansen,et al.  Electroencephalogram and visual evoked potential generation in a mathematical model of coupled cortical columns , 1995, Biological Cybernetics.

[20]  M. Haardt,et al.  A closed-form solution for multilinear PARAFAC decompositions , 2008, 2008 5th IEEE Sensor Array and Multichannel Signal Processing Workshop.

[21]  Arogyaswami Paulraj,et al.  An analytical constant modulus algorithm , 1996, IEEE Trans. Signal Process..

[22]  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..

[23]  Laurent Albera,et al.  Iterative methods for the canonical decomposition of multi-way arrays: Application to blind underdetermined mixture identification , 2011, Signal Process..