<< Chapter < Page Chapter >> Page >

We say that an n × N n times N matrix Φ Φ has the restricted isometry property (RIP) for k k if for each T { 1 , , N } T subseteq { lbrace {1 , dotslow , N} rbrace} such that # T k italic "#" T<= k , Φ T Φ_T (the matrix formed by choosing the columns of Φ Φ whose indices are in T T ) has the property

( 1 δ k ) x T 2 Φ T ( x ) 2 ( 1 + δ k ) x T 2 (RIP)

where 0 < δ k < 1 0<δ_k<1 . This useful definition is by Candes and Tao. The idea is that the embedding of a k k -dimensional space in M M -dimensional space almost preserves norm – like an isometry. Another way of looking at it is to consider the matrix Φ T t Φ T Φ_T^t Φ_T , of size k × k k times k . This matrix is symmetric, positive definite, and it’s eigen-values are between 1 δ k 1 - δ_k and 1 + δ k 1 +δ_k .

I prefer the following modified condition (dubbed the MIRP), which is more convenient for mathematical analysis:

( c 1 ) 1 x T 2 Φ T ( x ) 2 c 1 x T 2 (MRIP)

We can now state the following theorem.

If Φ Φ satisfies MRIP for 2 k 2 k then Δ exists Δ s.t. ( Φ , Δ ) \( {Φ , Δ} \) is instance optimal for 1 N ℓ_1^N for K K .

This shows that whenever we have a matrix Φ Φ satisfying the MRIP for 2 k 2 k then it will perform well on encoding vectors (at least in the sense of 1 N ℓ_1^N accuracy). The question is how can we construct measurement matrices with this property? We can construct Φ Φ using Gaussian entries and then normalizing the columns.

exists constant c > 0 c>0 s.t. if k c n log ( N n ) k<= c n over {l { \( {N ∕ n} \)}} then with high probability Φ Φ satisfies RIP and MRIP for k k .

Given N N and n n , the range of k k in the above results reflects how accurately we can recover data. There is another constant c c^′ that serves as a converse bound for Theorem 3. This converse can be derived using Gluskin widths.

The following generic problem is of great interest: Consider the class of matrices = { Φ M × N , Φ has some prescribed property(eg. Toeplitz, circulant, etc.) } . What is the largest k k for which such a matrix can have the MRIP.

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Compressive sensing. OpenStax CNX. Sep 21, 2007 Download for free at http://cnx.org/content/col10458/1.1
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Compressive sensing' conversation and receive update notifications?

Ask