Robust Low-rank Tensor Train Recovery

TL;DR

This paper introduces robust methods for recovering low-rank tensor train tensors from measurements contaminated with arbitrary outliers, using $ ext{l}_1$ loss and theoretical guarantees for linear measurement complexity.

Contribution

It establishes the $ ext{l}_1/ ext{l}_2$-RIP for Gaussian measurements, proves sharpness of the $ ext{l}_1$ loss, and develops two efficient algorithms with convergence guarantees.

Findings

The $ ext{l}_1/ ext{l}_2$-RIP holds with linear measurements.

Both proposed algorithms converge linearly with proper initialization.

Factorized method reduces memory cost significantly.

Abstract

Tensor train (TT) decomposition represents an $N$-order tensor using $O(N)$ matrices (i.e., factors) of small dimensions, achieved through products among these factors. Due to its compact representation, TT decomposition has found wide applications, including various tensor recovery problems in signal processing and quantum information. In this paper, we study the problem of reconstructing a TT format tensor from measurements that are contaminated by outliers with arbitrary values. Given the vulnerability of smooth formulations to corruptions, we use an $\ell_1$ loss function to enhance robustness against outliers. We first establish the $\ell_1/\ell_2$-restricted isometry property (RIP) for Gaussian measurement operators, demonstrating that the information in the TT format tensor can be preserved using a number of measurements that grows linearly with $N$. We also prove the sharpness property for the $\ell_1$ loss function optimized over TT format tensors. Building on the $\ell_1/\ell_2$-RIP and sharpness property, we then propose two complementary methods to recover the TT format tensor from the corrupted measurements: the projected subgradient method (PSubGM), which optimizes over the entire tensor, and the factorized Riemannian subgradient method (FRSubGM), which optimizes directly over the factors. Compared to PSubGM, the factorized approach FRSubGM significantly reduces the memory cost at the expense of a slightly slower convergence rate. Nevertheless, we show that both methods, with diminishing step sizes, converge linearly to the ground-truth tensor given an appropriate initialization, which can be obtained by a truncated spectral method.