On Uniqueness and Computation of the Decomposition of a Tensor into Multilinear Rank-$(1,L_r,L_r)$ Terms

In this paper we focus on the decomposition of a tensor $\mathcal T$ into a sum of multilinear rank-$(1,L_r,L_r)$ terms, $r=1,\dots,R$. This particular decomposition type has already found applications in wireless communication, chemometrics and the blind signal separation of signals that can be modelled as exponential polynomials and rational functions. We find conditions on the terms which guarantee that the decomposition is unique and can be computed by means of the eigenvalue decomposition of a matrix even in the cases where none of the factor matrices has full column rank. We consider both the case where the decomposition is exact and the case where the decomposition holds only approximately. We show that in both cases the number of the terms $R$ and their `sizes' $L_1,\dots,L_R$ do not have to be known a priori and can be estimated as well. The conditions for uniqueness are easy to verify, especially for terms that can be considered `generic'. In particular, we obtain the following two generalizations of a well known result on generic uniqueness of the CPD (i.e., the case $L_1=\dots=L_R=1$): we show that the multilinear rank-$(1,L_r,L_r)$ decomposition of an $I\times J\times K$ tensor is generically unique if i) $L_1=\dots=L_R=:L$ and $R\leq \min((J-L)(K-L),I)$ or if ii) $\sum L_R\leq \min((I-1)(J-1),K)$ and $J\geq \max(L_i+L_j)$.