Restricted isometry property

From testwiki
Jump to navigation Jump to search

In linear algebra, the restricted isometry property (RIP) characterizes matrices which are nearly orthonormal, at least when operating on sparse vectors. The concept was introduced by Emmanuel Candès and Terence Tao[1] and is used to prove many theorems in the field of compressed sensing.[2] There are no known large matrices with bounded restricted isometry constants (computing these constants is strongly NP-hard,[3] and is hard to approximate as well[4]), but many random matrices have been shown to remain bounded. In particular, it has been shown that with exponentially high probability, random Gaussian, Bernoulli, and partial Fourier matrices satisfy the RIP with number of measurements nearly linear in the sparsity level.[5] The current smallest upper bounds for any large rectangular matrices are for those of Gaussian matrices.[6] Web forms to evaluate bounds for the Gaussian ensemble are available at the Edinburgh Compressed Sensing RIC page.[7]

Definition

Let A be an m × p matrix and let 1 ≤ s ≤ p be an integer. Suppose that there exists a constant δs(0,1) such that, for every m × s submatrix As of A and for every s-dimensional vector y,

(1δs)y22Asy22(1+δs)y22.

Then, the matrix A is said to satisfy the s-restricted isometry property with restricted isometry constant δs.

This condition is equivalent to the statement that for every m × s submatrix As of A we have

As*AsIs×s22δs,

where Is×s is the s×s identity matrix and X22 is the operator norm. See for example [8] for a proof.

Finally this is equivalent to stating that all eigenvalues of As*As are in the interval [1δs,1+δs].

Restricted Isometric Constant (RIC)

The RIC Constant is defined as the infimum of all possible δ for a given An×m.

δK=inf[δ:(1δ)y22Asy22(1+δ)y22], |s|K,yR|s|

It is denoted as δK.

Eigenvalues

For any matrix that satisfies the RIP property with a RIC of δK, the following condition holds:[1]

1δKλmin(Aτ*Aτ)λmax(Aτ*Aτ)1+δK.

The tightest upper bound on the RIC can be computed for Gaussian matrices. This can be achieved by computing the exact probability that all the eigenvalues of Wishart matrices lie within an interval.

See also

References

Template:Reflist

  1. 1.0 1.1 E. J. Candes and T. Tao, "Decoding by Linear Programming," IEEE Trans. Inf. Th., 51(12): 4203–4215 (2005).
  2. E. J. Candes, J. K. Romberg, and T. Tao, "Stable Signal Recovery from Incomplete and Inaccurate Measurements," Communications on Pure and Applied Mathematics, Vol. LIX, 1207–1223 (2006).
  3. A. M. Tillmann and M. E. Pfetsch, "The Computational Complexity of the Restricted Isometry Property, the Nullspace Property, and Related Concepts in Compressed Sensing," IEEE Trans. Inf. Th., 60(2): 1248–1259 (2014)
  4. Abhiram Natarajan and Yi Wu, "Computational Complexity of Certifying Restricted Isometry Property," Approximation, Randomization, and Combinatorial Optimization. Algorithms and Techniques (APPROX/RANDOM 2014) (2014)
  5. F. Yang, S. Wang, and C. Deng, "Compressive sensing of image reconstruction using multi-wavelet transform", IEEE 2010
  6. B. Bah and J. Tanner "Improved Bounds on Restricted Isometry Constants for Gaussian Matrices"
  7. Template:Cite web
  8. Template:Cite web
  9. Template:Cite web
  10. Template:Cite journal