Abstract

Alternating least squares (ALS) and its variations are the most commonly used algorithms for the PARAFAC decomposition of a tensor. However, it is still troubled for one how to accelerate the ALS algorithm with the reduced computational complexity. In this paper, a new acceleration method for the ALS with a matrix polynomial predictive model (MPPM) is proposed. In the MPPM, a matrix-valued function is first approximated by a matrix polynomial. It is shown that the future value of the function can be predicted by an FIR filter with the coefficients determined offline. By viewing each factor matrix of a tensor as a matrix-valued function, a new ALS algorithm, the ALS-MPPM algorithm, is then given. Analyses show that our ALS-MPPM algorithm is of low computational complexity and a close relation with the existing ALS algorithms. Moreover, to further accelerate the convergence of the proposed algorithm, a new technique called the multi-model (MM) prediction is also introduced. While the analytical results are verified by the numerical simulations, it is also shown that our ALS-MPPM outperforms the existing ALS-based algorithms in terms of the rate of convergence.

Highlights

  • It is well known that tensors, or multidimensional arrays, are commonly used in many areas such as the Direct SequenceCode Division Multiple Access (DS-CDMA) [1], Multiple Input Multiple Output (MIMO) radar [2], [3], Space Time Adaptive Processing (STAP) [4], and Multidimensional Harmonic Retrieval (MHR) [5]

  • The Matrix Polynomial Predictive Model (MPPM) approximates a matrix-valued function by a matrix polynomial

  • It is shown that the future value of a matrix polynomial can be predicted using an FIR filter, whose coefficients can be determined offline

Read more

Summary

INTRODUCTION

It is well known that tensors, or multidimensional arrays, are commonly used in many areas such as the Direct SequenceCode Division Multiple Access (DS-CDMA) [1], Multiple Input Multiple Output (MIMO) radar [2], [3], Space Time Adaptive Processing (STAP) [4], and Multidimensional Harmonic Retrieval (MHR) [5]. To make use of the structural advantages brought by the multidimensional data, tensor based methods for modeling and/or processing them have been reported for such applications in the recent years [6] [7], [8] Among those methods, the PARAllel FACtor (PARAFAC) decomposition [9], known as the Canonical Polyadic Decomposition (CPD) [10], or the CANonical DECOMPosition (CANDECOMP) [11], is one of the most widely used tools. The most commonly used algorithm for the PARAFAC decomposition is the Alternating Least Squares (ALS) [9], [11]. A new acceleration technique for the ALS algorithm, called as the Matrix Polynomial Predictive Model (MPPM), is proposed. For other forms of the PARAFAC decomposition, please refer to [12], [34] for details

ALTERNATING LEAST SQUARES
MATRIX-VALUED POLYNOMIAL MODEL
COEFFICIENTS DERIVATION
MULTI-MODEL PREDICTION
COMPUTATIONAL COMPLEXITY
NUMERICAL SIMULATIONS
DIFFERENCE CONFIGURATIONS
CONCLUSION
COMPUTATIONAL COMPLEXITY OF FAST ALS UPDATE
COMPUTATIONAL COMPLEXITY OF MPPM
Findings
COMPUTATIONAL COMPLEXITY OF ELS
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call