Under-Determined tensor diagonalization for decomposition of difficult tensors

Analysis of multidimensional arrays, usually called tensors, often becomes difficult in cases when the tensor rank (a minimum number of rank-one components) exceeds all the tensor dimensions. Traditional methods of canonical polyadic decomposition of such tensors, namely the alternating least squares, can be used, but a presence of a large number of false local minima can make the problem hard. Usually, multiple random initializations are advised in such cases, but the question is how many such random initializations are sufficient to get a good chance of finding the right solution. It appears that the number of the initializations can be very large. We propose a novel approach to the problem. The given tensor is augmented by some unknown parameters to the shape that admits ordinary tensor diagonalization, i.e., transforming the augmented tensor into an exact or nearly diagonal form through multiplying the tensor by non-orthogonal invertible matrices. Three possible constraints are proposed to make the optimization problem well defined. The method can be modified for an under-determined block-term decomposition.

[1]  Andrzej Cichocki,et al.  Non-orthogonal tensor diagonalization , 2014, Signal Process..

[2]  Eric Moreau,et al.  Joint Matrices Decompositions and Blind Source Separation , 2014 .

[3]  Julian D. Laderman,et al.  A noncommutative algorithm for multiplying $3 \times 3$ matrices using 23 multiplications , 1976 .

[4]  Andrzej Cichocki,et al.  Numerical CP decomposition of some difficult tensors , 2016, J. Comput. Appl. Math..

[5]  Zbynek Koldovský,et al.  Weight Adjusted Tensor Method for Blind Separation of Underdetermined Mixtures of Nonstationary Sources , 2011, IEEE Transactions on Signal Processing.

[6]  P. Comon Tensor Diagonalization, A useful Tool in Signal Processing , 1994 .

[7]  D. D. Morrison Methods for nonlinear least squares problems and convergence proofs , 1960 .

[8]  Markus Bläser,et al.  Fast Matrix Multiplication , 2013, Theory Comput..

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

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

[11]  P. Comon,et al.  TENSOR DIAGONALIZATION BY ORTHOGONAL TRANSFORMS , 2007 .

[12]  Martin Haardt,et al.  Extension of the semi-algebraic framework for approximate CP decompositions via non-symmetric simultaneous matrix diagonalization , 2016, 2016 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[13]  A. Smirnov,et al.  The bilinear complexity and practical algorithms for matrix multiplication , 2013 .

[14]  Lieven De Lathauwer,et al.  Blind Identification of Underdetermined Mixtures by Simultaneous Matrix Diagonalization , 2008, IEEE Transactions on Signal Processing.

[15]  Tamara G. Kolda,et al.  Tensor Decompositions and Applications , 2009, SIAM Rev..

[16]  Andrzej Cichocki,et al.  Two-sided diagonalization of order-three tensors , 2015, 2015 23rd European Signal Processing Conference (EUSIPCO).

[17]  L. Lathauwer,et al.  Canonical polyadic decomposition of third-order tensors: relaxed uniqueness conditions and algebraic algorithm , 2015, 1501.07251.