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

      Low rank approximation and decomposition of large matrices using error correcting codes

      Preprint
      , ,

      Read this article at

          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

          Low rank approximation is an important tool used in many applications of signal processing and machine learning. Recently, randomized sketching algorithms were proposed to effectively construct low rank approximations and obtain approximate singular value decompositions of large matrices. Similar ideas were used to solve least squares regression problems. In this paper, we show how matrices from error correcting codes can be used to find such low rank approximations and matrix decompositions, and extend the framework to linear least squares regression problems. The benefits of using these code matrices are the following: (i) They are easy to generate and they reduce randomness significantly. (ii) Code matrices with mild properties satisfy the subspace embedding property, and have a better chance of preserving the geometry of an entire subspace of vectors. (iii) For parallel and distributed applications, code matrices have significant advantages over structured random matrices and Gaussian random matrices. (iv) Unlike Fourier or Hadamard transform matrices, which require sampling \(O(k\log k)\) columns for a rank-\(k\) approximation, the log factor is not necessary for certain types of code matrices. That is, \((1+\epsilon)\) optimal Frobenius norm error can be achieved for a rank-\(k\) approximation with \(O(k/\epsilon)\) samples. (v) Fast multiplication is possible with structured code matrices, so fast approximations can be achieved for general dense input matrices. (vi) For least squares regression problem \(\min\|Ax-b\|_2\) where \(A\in \mathbb{R}^{n\times d}\), the \((1+\epsilon)\) relative error approximation can be achieved with \(O(d/\epsilon)\) samples, with high probability, when certain code matrices are used.

          Related collections

          Author and article information

          Journal
          2015-12-30
          2016-08-19
          Article
          1512.09156
          276f6b99-6442-4e70-88fe-4141ec6a95c4

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

          History
          Custom metadata
          cs.IT cs.LG cs.NA math.IT

          Numerical & Computational mathematics,Numerical methods,Information systems & theory,Artificial intelligence

          Comments

          Comment on this article