A Regularized Nonnegative Third Order Tensor decomposition Using a Primal-Dual Projected Gradient Algorithm: Application to 3D Fluorescence Spectroscopy

This paper investigates the use of Primal-Dual optimization algorithms on multidimensional signal processing problems. The data blocks interpreted in a tensor way can be modeled by means of multi-linear decomposition. Here we will focus on the Canonical Polyadic Decomposition (CPD), and we will present an application to fluorescence spectroscopy using this decomposition. In order to estimate the factors or latent variables involved in these decompositions, it is usual to use criteria optimization algorithms. A classical cost function consists of a measure of the modeling error (fidelity term) to which a regularization term can be added if necessary. Here, we consider one of the most efficient optimization methods, Primal-Dual Projected Gradient.

[1]  Jorge N. Tendeiro,et al.  First and second-order derivatives for CP and INDSCAL , 2011 .

[2]  J. Chang,et al.  Analysis of individual differences in multidimensional scaling via an n-way generalization of “Eckart-Young” decomposition , 1970 .

[3]  R. Harshman,et al.  PARAFAC: parallel factor analysis , 1994 .

[4]  F. L. Hitchcock The Expression of a Tensor or a Polyadic as a Sum of Products , 1927 .

[5]  Pierre Comon,et al.  Computing the polyadic decomposition of nonnegative third order tensors , 2011, Signal Process..

[6]  J Möcks,et al.  Topographic components model for event-related potentials and some biophysical considerations. , 1988, IEEE transactions on bio-medical engineering.

[7]  I. Ekeland,et al.  Convex analysis and variational problems , 1976 .

[8]  Richard A. Harshman,et al.  Foundations of the PARAFAC procedure: Models and conditions for an "explanatory" multi-model factor analysis , 1970 .

[9]  J. Mocks,et al.  Topographic components model for event-related potentials and some biophysical considerations , 1988, IEEE Transactions on Biomedical Engineering.

[10]  Antonin Chambolle,et al.  A First-Order Primal-Dual Algorithm for Convex Problems with Applications to Imaging , 2011, Journal of Mathematical Imaging and Vision.

[11]  Andrzej Cichocki,et al.  Nonnegative Matrix and Tensor Factorization T , 2007 .

[12]  Mingqiang Zhu,et al.  An Efficient Primal-Dual Hybrid Gradient Algorithm For Total Variation Image Restoration , 2008 .

[13]  C. G. Bollini,et al.  On the Reduction Formula of Feinberg and Pais , 1965 .

[14]  Stefan Kindermann,et al.  News Algorithms for tensor decomposition based on a reduced functional , 2014, Numer. Linear Algebra Appl..

[15]  Tamara G. Kolda,et al.  Scalable Tensor Factorizations with Missing Data , 2010, SDM.

[16]  David E. Booth,et al.  Multi-Way Analysis: Applications in the Chemical Sciences , 2005, Technometrics.

[17]  P. Coble Characterization of marine and terrestrial DOM in seawater using excitation-emission matrix spectroscopy , 1996 .