18
views
0
recommends
+1 Recommend
0 collections
    0
    shares
      • Record: found
      • Abstract: found
      • Article: found
      Is Open Access

      Canonical polyadic decomposition of third-order tensors: relaxed uniqueness conditions and algebraic algorithm

      Preprint
      ,

      Read this article at

      Bookmark
          There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.

          Abstract

          Canonical Polyadic Decomposition (CPD) of a third-order tensor is a minimal decomposition into a sum of rank-1 tensors. We find new mild deterministic conditions for the uniqueness of individual rank-1 tensors in CPD and present an algorithm to recover them. We call the algorithm "algebraic" because it relies only on standard linear algebra. It does not involve more advanced procedures than the computation of the null space of a matrix and eigen/singular value decomposition. Simulations indicate that the new conditions for uniqueness and the working assumptions for the algorithm hold for a randomly generated I×J×K tensor of rank RKJI2 if R is bounded as R(I+J+K2)/2+(K(IJ)2+4K)/2 at least for the dimensions that we have tested. This improves upon the famous Kruskal bound for uniqueness R(I+J+K2)/2 as soon as I3. In the particular case R=K, the new bound above is equivalent to the bound R(I1)(J1) which is known to be necessary and sufficient for the generic uniqueness of the CPD. An existing algebraic algorithm (based on simultaneous diagonalization of a set of matrices) computes the CPD under the more restrictive constraint R(R1)I(I1)J(J1)/2 (implying that R<(J12)(I12)/2+1). On the other hand, optimization-based algorithms fail to compute the CPD in a reasonable amount of time even in the low-dimensional case I=3, J=7, K=R=12. By comparison, in our approach the computation takes less than 1 sec. We demonstrate that, at least for R24, our algorithm can recover the rank-1 tensors in the CPD up to R(I1)(J1).

          Related collections

          Author and article information

          Journal
          2015-01-28
          2016-07-19
          Article
          1501.07251
          6e51798d-60da-4d5a-8817-9fe94e315680

          http://arxiv.org/licenses/nonexclusive-distrib/1.0/

          History
          Custom metadata
          15A69, 15A23
          math.SP

          Functional analysis
          Functional analysis

          Comments

          Comment on this article