Convex and nonconvex geometries of symmetric tensor factorization

Tensors provide natural representations for massive multi-mode datasets and tensor methods also form the backbone of many machine learning, signal processing, and statistical algorithms. This work develops theories and computational methods for guaranteed overcomplete, non-orthogonal symmetric tensor factorization using convex and nonconvex optimizations. In particular, we show when the symmetric tensor factors are uniformly sampled from the unit sphere, they are provably recoverable using convex atomic norm minimization. To design scalable polynomial-time algorithms, we apply low-rank parameterization to reformulate the atomic norm regularized tensor optimization as a nonconvex program. We analyze the optimization landscape of this nonconvex program to ensure (local) convergence of gradient descent algorithms.

[1]  Jieping Ye,et al.  Tensor Completion for Estimating Missing Values in Visual Data , 2013, IEEE Trans. Pattern Anal. Mach. Intell..

[2]  Annie Chen,et al.  Context-Aware Collaborative Filtering System: Predicting the User's Preference in the Ubiquitous Computing Environment , 2005, LoCA.

[3]  Michael Elad,et al.  Optimally sparse representation in general (nonorthogonal) dictionaries via ℓ1 minimization , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[4]  Zhihui Zhu,et al.  Global Optimality in Low-Rank Matrix Optimization , 2017, IEEE Transactions on Signal Processing.

[5]  M. Rudelson,et al.  Lp-moments of random vectors via majorizing measures , 2005, math/0507023.

[6]  Zhihui Zhu,et al.  Geometry of Factored Nuclear Norm Regularization , 2017, ArXiv.

[7]  Gongguo Tang,et al.  Approximate support recovery of atomic line spectral estimation: A tale of resolution and precision , 2016, 2016 IEEE Global Conference on Signal and Information Processing (GlobalSIP).

[8]  Pierre Comon,et al.  Tensor Decompositions, State of the Art and Applications , 2002 .

[9]  J. Kruskal Three-way arrays: rank and uniqueness of trilinear decompositions, with application to arithmetic complexity and statistics , 1977 .

[10]  Parikshit Shah,et al.  Guaranteed Tensor Decomposition: A Moment Approach , 2015, ICML.

[11]  Pierre Comon,et al.  Multiarray Signal Processing: Tensor decomposition meets compressed sensing , 2010, ArXiv.

[12]  Anima Anandkumar,et al.  Guaranteed Non-Orthogonal Tensor Decomposition via Alternating Rank-1 Updates , 2014, ArXiv.

[13]  Michael I. Jordan,et al.  Gradient Descent Converges to Minimizers , 2016, ArXiv.

[14]  Gongguo Tang,et al.  The nonconvex geometry of low-rank matrix optimizations with general objective functions , 2017, 2017 IEEE Global Conference on Signal and Information Processing (GlobalSIP).

[15]  Lixin Shen,et al.  Overcomplete tensor decomposition via convex optimization , 2015, 2015 IEEE 6th International Workshop on Computational Advances in Multi-Sensor Adaptive Processing (CAMSAP).

[16]  Christopher J. Hillar,et al.  Most Tensor Problems Are NP-Hard , 2009, JACM.

[17]  Gongguo Tang,et al.  The nonconvex geometry of low-rank matrix optimizations with general objective functions , 2016, 2017 IEEE Global Conference on Signal and Information Processing (GlobalSIP).

[18]  Renato D. C. Monteiro,et al.  A nonlinear programming algorithm for solving semidefinite programs via low-rank factorization , 2003, Math. Program..

[19]  Alexander Barvinok,et al.  A course in convexity , 2002, Graduate studies in mathematics.

[20]  Gongguo Tang,et al.  Resolution limits for atomic decompositions via Markov-Bernstein type inequalities , 2015, 2015 International Conference on Sampling Theory and Applications (SampTA).

[21]  Parikshit Shah,et al.  Compressed Sensing Off the Grid , 2012, IEEE Transactions on Information Theory.

[22]  Zhihui Zhu,et al.  The Global Optimization Geometry of Nonsymmetric Matrix Factorization and Sensing , 2017, ArXiv.