On uniqueness and computation of the decomposition of a tensor into multilinear rank-$(1, L_{r},L_{r})$ terms

TL;DR

This paper investigates the uniqueness and computation of tensor decompositions into multilinear rank-$(1,L_r,L_r)$ terms, extending classical results on CPD and providing practical conditions for applications in signal processing and chemometrics.

Contribution

It introduces new conditions for the uniqueness and computability of multilinear rank-$(1,L_r,L_r)$ tensor decompositions, generalizing known CPD results.

Findings

Decomposition is unique under specific rank and size conditions.

Eigenvalue decomposition can compute the decomposition even without full column rank.

Number of terms and their sizes can be estimated without prior knowledge.

Abstract

Canonical Polyadic Decomposition (CPD) represents a third-order tensor as the minimal sum of rank-1 terms. Because of its uniqueness properties the CPD has found many concrete applications in telecommunication, array processing, machine learning, etc. On the other hand, in several applications the rank-1 constraint on the terms is too restrictive. A multilinear rank-$(M,N,L)$ constraint (where a rank-1 term is the special case for which $M=N=L=1$) could be more realistic, while it still yields a decomposition with attractive uniqueness properties. 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,...,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,...,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=...=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=...=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)$.