Table of Contents
Fetching ...

Fast Orthogonal Matching Pursuit through Successive Regression

Huiyuan Yu, Jia He, Maggie Cheng

TL;DR

The paper tackles the computational bottlenecks of Orthogonal Matching Pursuit in sparse recovery by introducing OMP-SR, a fast successive-regression variant that preserves the same atom selection as OMP but avoids costly LS inversions via a backtracking coefficient update. It extends the idea to blocks with Blocked Successive Regression (BSR), reducing iterations and enabling faster recovery comparable to or superior to gOMP while retaining recovery guarantees. The authors provide both strong and weak exact-recovery conditions, connecting these guarantees to dictionary coherence and cumulative coherence measures, and show that under conditions like $\mu(2k-1) < 1$ or $\mu_1(l)+\mu_1(n) < 1$, BSR achieves recovery in $\lceil k/c \rceil$ iterations. Experiments on real and synthetic data demonstrate significant speedups with comparable accuracy under noise, highlighting the practical utility of the proposed fast OMP-based algorithms for large-scale sparse representation problems.

Abstract

Orthogonal Matching Pursuit (OMP) has been a powerful method in sparse signal recovery and approximation. However, OMP suffers computational issues when the signal has a large number of non-zeros. This paper advances OMP and its extension called generalized OMP (gOMP) by offering fast algorithms for the orthogonal projection of the input signal at each iteration. The proposed modifications directly reduce the computational complexity of OMP and gOMP. Experiment results verified the improvement in computation time. This paper also provides sufficient conditions for exact signal recovery. For general signals with additive noise, the approximation error is at the same order as OMP (gOMP), but is obtained within much less time.

Fast Orthogonal Matching Pursuit through Successive Regression

TL;DR

The paper tackles the computational bottlenecks of Orthogonal Matching Pursuit in sparse recovery by introducing OMP-SR, a fast successive-regression variant that preserves the same atom selection as OMP but avoids costly LS inversions via a backtracking coefficient update. It extends the idea to blocks with Blocked Successive Regression (BSR), reducing iterations and enabling faster recovery comparable to or superior to gOMP while retaining recovery guarantees. The authors provide both strong and weak exact-recovery conditions, connecting these guarantees to dictionary coherence and cumulative coherence measures, and show that under conditions like or , BSR achieves recovery in iterations. Experiments on real and synthetic data demonstrate significant speedups with comparable accuracy under noise, highlighting the practical utility of the proposed fast OMP-based algorithms for large-scale sparse representation problems.

Abstract

Orthogonal Matching Pursuit (OMP) has been a powerful method in sparse signal recovery and approximation. However, OMP suffers computational issues when the signal has a large number of non-zeros. This paper advances OMP and its extension called generalized OMP (gOMP) by offering fast algorithms for the orthogonal projection of the input signal at each iteration. The proposed modifications directly reduce the computational complexity of OMP and gOMP. Experiment results verified the improvement in computation time. This paper also provides sufficient conditions for exact signal recovery. For general signals with additive noise, the approximation error is at the same order as OMP (gOMP), but is obtained within much less time.
Paper Structure (22 sections, 11 theorems, 24 equations, 2 figures, 5 tables, 2 algorithms)

This paper contains 22 sections, 11 theorems, 24 equations, 2 figures, 5 tables, 2 algorithms.

Key Result

theorem 1

A sufficient condition for BSR to recover a $k$-sparse signal within $\lceil k/c\rceil$ iterations is that holds for all iterations.

Figures (2)

  • Figure 1: Images used for experiments, (a) phantom, (b) transaxial CT, (c) trees, (d) letters.
  • Figure 2: Running time (top) and iterations (bottom) used by the algorithms to recover $k$ non-zeros in the signal. OMP and OMP-SR use the same number of iterations, and gOMP and BSR use the same number of iterations. The datasets used: (a),(c) phantom; (b),(d) MNIST dataset handwritten digit '7'.

Theorems & Definitions (11)

  • theorem 1: The strong exact recovery condition for BSR
  • theorem 2: The Weak Exact Recovery Condition for BSR
  • lemma 1
  • lemma 2
  • theorem 3: The strong exact recovery condition for BSR
  • theorem 1: The strong exact recovery condition for BSR
  • theorem 2: The Weak Exact Recovery Condition for BSR
  • lemma 1
  • lemma 2
  • theorem 3: The strong exact recovery condition for BSR
  • ...and 1 more