Basis pursuit denoising explained

In applied mathematics and statistics, basis pursuit denoising (BPDN) refers to a mathematical optimization problem of the form

minx\left(

1
2

\|y-

2
Ax\|
2

+λ\|x\|1\right),

where

λ

is a parameter that controls the trade-off between sparsity and reconstruction fidelity,

x

is an

N x 1

solution vector,

y

is an

M x 1

vector of observations,

A

is an

M x N

transform matrix and

M<N

. This is an instance of convex optimization.

Some authors refer to basis pursuit denoising as the following closely related problem:

minx\|x\|1subjectto\|y-

2
Ax\|
2

\le\delta,

which, for any given

λ

, is equivalent to the unconstrained formulation for some (usually unknown a priori) value of

\delta

. The two problems are quite similar. In practice, the unconstrained formulation, for which most specialized and efficient computational algorithms are developed, is usually preferred.

Either types of basis pursuit denoising solve a regularization problem with a trade-off between having a small residual (making

y

close to

Ax

in terms of the squared error) and making

x

simple in the

\ell1

-norm sense. It can be thought of as a mathematical statement of Occam's razor, finding the simplest possible explanation (i.e. one that yields

minx\|x\|1

) capable of accounting for the observations

y

.

Exact solutions to basis pursuit denoising are often the best computationally tractable approximation of an underdetermined system of equations. Basis pursuit denoising has potential applications in statistics (see the LASSO method of regularization), image compression and compressed sensing.

When

\delta=0

, this problem becomes basis pursuit.

Basis pursuit denoising was introduced by Chen and Donoho in 1994,[1] in the field of signal processing. In statistics, it is well known under the name LASSO, after being introduced by Tibshirani in 1996.

Solving basis pursuit denoising

The problem is a convex quadratic problem, so it can be solved by many general solvers, such as interior-point methods. For very large problems, many specialized methods that are faster than interior-point methods have been proposed.

Several popular methods for solving basis pursuit denoising include the in-crowd algorithm (a fast solver for large, sparse problems[2]), homotopy continuation, fixed-point continuation (a special case of the forward–backward algorithm[3]) and spectral projected gradient for L1 minimization (which actually solves LASSO, a related problem).

External links

Notes and References

  1. Book: Shaobing . Chen . D. . Donoho . Basis pursuit . Proceedings of 1994 28th Asilomar Conference on Signals, Systems and Computers . 1994 . 1 . 41–44 . 10.1109/ACSSC.1994.471413 . 0-8186-6405-3 . 96447294 .
  2. See Patrick R. . Gill . Albert . Wang . Alyosha . Molnar . The In-Crowd Algorithm for Fast Basis Pursuit Denoising . . 59 . 10 . 2011 . 4595–4605 . 10.1109/TSP.2011.2161292 . 15320645 . . demo MATLAB code available http://molnargroup.ece.cornell.edu/files/InCrowdBeta1.zip.
  3. Web site: Forward Backward Algorithm . February 16, 2014 . https://web.archive.org/web/20140216231347/http://www.ugcs.caltech.edu/~srbecker/wiki/Forward_Backward_Algorithm .