---
title: Restricted Isometry Constant (RIC) in Compressed Sensing
url: https://www.emergentmind.com/topics/restricted-isometry-constant-ric
type: topic
---

# Restricted Isometry Constant (RIC) in Compressed Sensing

The restricted isometry constant (RIC) is a central quantitative metric in compressed sensing, sparse signal recovery, and high-dimensional data analysis. It measures to what extent a matrix acts nearly isometrically when restricted to sparse vectors, providing necessary and sufficient conditions for the success of various recovery algorithms, both convex and nonconvex. RIC thresholds govern phase transitions, inform algorithm selection, and underlie sample complexity analyses in modern compressed sensing theory. The RIC framework has shaped both theoretical understanding and algorithmic design, but is also marked by significant challenges in its computation, estimation, and sharpness for diverse matrix classes and recovery methods.

## 1. Formal Definition and Basic Properties

Given a matrix $A\in\mathbb{R}^{m\times n}$ and integer $s\geq 1$, the $s$-th order restricted isometry constant (RIC) $\delta_s(A)$ is defined as the smallest $\delta\geq0$ such that, for all $s$-sparse vectors $x$,
\[
(1-\delta)\|x\|_2^2 \leq \|Ax\|_2^2 \leq (1+\delta)\|x\|_2^2
\]
Equivalently,
\[
1-\delta \leq \frac{\|Ax\|_2^2}{\|x\|_2^2} \leq 1+\delta \qquad \forall x: \|x\|_0 \leq s
\]
RICs are monotone in order: if $s_1\leq s_2$ then $\delta_{s_1}(A)\leq \delta_{s_2}(A)$. They can also be expressed in eigenvalue language: for each subset $T\subset[n]$, $|T|=s$, let $A_T$ denote the associated submatrix. Then,
\[
\delta_s(A) = \max \left\{ 1 - \min_{|T|=s} \lambda_{\min}(A_T^T A_T),\ \max_{|T|=s} \lambda_{\max}(A_T^T A_T) - 1 \right\}
\]
The minimal $\delta$ for which $A$ is said to have the restricted isometry property (RIP) of order $s$ and level $\delta$.

## 2. RIC in Sparse Recovery: Thresholds and the Null-Space Property

RICs constitute a critical metric guaranteeing exact and stable recovery of $k$-sparse vectors. For recovery via (potentially nonconvex) $\ell_q$-minimization ($0 < q \leq 1$):
\[
x^* = \arg\min_{x\in \mathbb{R}^n} \|x\|_q^q \quad \text{subject to} \quad A x = y
\]
Necessary and sufficient conditions are formulated via the null-space property (NSP):
\[
\|h_T\|_q^q < \|h_{T^c}\|_q^q \quad \forall\ h\in\ker(A)\setminus\{0\},\ |T| \leq k
\]
RIC lower bounds are obtained by relating the NSP to the RIC via sharp inequalities between $\ell_q$ and $\ell_1$ norms. For $q=1$ or $q=1/2$ (including both convex and certain archetypal nonconvex relaxations), the exact-recovery threshold is:
\[
\delta_{tk} < \sqrt{\frac{t-1}{t}} \quad \forall\, t > 1
\]
This threshold is sharp for $\ell_1$-minimization and, as recently established, for $\ell_{1/2}$-minimization as well, showing an equivalence in recovery capability between these relaxations under the same RIC bound [1308.0455].

For more general $q \in (0,1)$, nontrivial improvements are obtained on $\delta_{2k}$. Explicit examples:
- For $0<q<1/2$, if $\delta_{2k} < 0.5547$ then recovery holds for all $k \geq 2$.
- For $1/2<q<1$, if $\delta_{2k}<0.6782$ then recovery holds for all $k\geq1$ [1308.0455].
Furthermore, in the minimal-sparsity regime, for $q=1/2, 1$, RIC bounds specialize to $\delta_k<1/3$ (even $k$) and $\delta_k<0.3203$ (odd $k$), pinpointing the sharpest-known thresholds.

## 3. RIC Bounds for Recovery Algorithms

The RIC framework enables the comparison of various recovery algorithms:

| Algorithm                   | Sufficient RIC Threshold ($\delta$)                      | Reference                  |
|-----------------------------|----------------------------------------------------------|----------------------------|
| $\ell_1$-minimization       | $\delta_{tk}<\sqrt{(t-1)/t},\ t>1$                       | [1308.0455]                |
| OMP                         | $\delta_{k+1}<1/(\sqrt{k}+1)$                            | [1201.2733], [1501.01708]  |
| gOMP                        | $\delta_{N K+1}<1/\sqrt{K/N + 1}$                        | [1604.03306]               |
| SP                          | $\delta_{3k}<0.4859$                                     | [1309.6073]                |
| CoSaMP                      | $\delta_{4k}<0.5$                                        | [1309.6073]                |
| $\ell_q$-min ($0<q<1$)      | see explicit bounds in [1308.0455], [1312.3379], [1310.2410] |                            |

These thresholds often match, or are provably sharp for, the breakdown points of the corresponding algorithms. For instance, OMP’s bound is sharp: a matrix can have $\delta_{k+1}=1/\sqrt{k+1}$ and yet OMP will fail on some $k$-sparse vector. The same tightness applies to gOMP and certain block-sparse reductions [1501.01708, 1604.03306].

## 4. Computational Complexity and Practical Estimation

Determining whether a given matrix $A$ satisfies the RIP for specified $k$ and $\delta$, or computing $\delta_k(A)$ exactly, is strongly intractable. Both the exact computation and decision versions of the RIC are NP-hard, and even coNP-complete for some formulations [1205.2081]. No polynomial-time (or even pseudo-polynomial time) approximation within arbitrarily small factors exists unless P=NP. This computational barrier arises via direct reduction from the spark problem, which is also NP-complete: determining whether $A$ has a dependent subset of $\leq k$ columns.

As a result, practical and theoretical work on RIC estimation focuses on:
- Efficient upper and lower bounds, often via semidefinite relaxations, convex programming, or probabilistic methods.
- Probabilistic guarantees for random matrix ensembles, such as Gaussian or subgaussian matrices, where high-probability uniform RIP is attainable with $\mathcal O(k\log n)$ rows.

Due to NP-hardness, almost all practical recovery and analysis methods rely on proxies or replace worst-case RICs with high-probability bounds for structured or random $A$.

## 5. RIC Bounds for Random and Structured Matrices

The typical mechanism for bounding RICs in random or structured ensembles is via concentration inequalities for the extreme singular values (or equivalently, eigenvalues) of submatrices.

Notable classes and scaling laws:
- **Gaussian/Random Matrices**: With probability at least $1-\epsilon$, $A\sim \mathcal{N}(0,1)$ satisfies $\delta_k(A) < \delta$ when $m \gtrsim C(\delta) k\log(n/k)$. Sharper large-deviation and replica-based bounds further optimize $C(\delta)$ and asymptotic tightness [1207.4883, 1306.3779, 1501.06281, 1410.1956].
- **Partial Random Circulant/Gabor**: For such structured matrices, $m \gtrsim c \delta^{-2} s\log^2 s\log^2 n$ rows suffice to ensure $\delta_s \leq \delta$ with overwhelming probability [1207.0235].
- **Khatri-Rao/Kronecker Products**: The RIC of a Khatri-Rao product $A\odot A$ is at most $\delta_k(A)^2$, making such product matrices strictly stronger isometries. For Kronecker products, $\delta_s(H_1 \otimes H_2)\leq (1+\delta_s(H_1))(1+\delta_s(H_2))-1$ [1709.05789, 2409.08699].

Extreme-value theory and random matrix concentration provide not only the scaling laws but the precise finite-$n$ performance, including sharp constants and the distributions of left/right RICs. Notably, in the Gaussian ensemble, left and right RICs asymptotically follow Weibull (min) and Gumbel (max) limit laws, respectively [1410.1956].

## 6. Role in Compressed Sensing Phase Transitions and Algorithmic Design

The RIC determines the algorithmic boundary (“phase transition”) for which classes of sparse recovery succeed uniformly. For example, in compressed sensing:
- Sharp RICs imply that for $\delta_{2k}<\sqrt{2}-1$, $\ell_1$-minimization exactly recovers all $k$-sparse vectors.
- For OMP, exact (and only exact) recovery for every $k$-sparse signal in $k$ steps is achievable when $\delta_{k+1}<1/(\sqrt{k}+1)$ [1201.2733, 1101.4458, 1501.01708].
- Extensions to block-sparsity, joint sparsity, and group-structured dictionaries similarly depend on structured RICs.

In practical design, knowledge of attainable RIC levels in candidate matrix ensembles guides the selection of $m$ (row count) required for reliable recovery, and dictates when higher computational-complexity algorithms (e.g., convex programming) are justified versus fast greedy algorithms (e.g., OMP, gOMP, SP, CoSaMP).

## 7. Open Problems and Future Directions

Despite rigorous advances, several open questions remain:
- **Sharpness for General Classes**: While exact sharp bounds are available for some algorithms and ensembles, closing gaps between sufficiency and necessity for $\delta_{sk}$ (with $s>2$ or for nonconvex $q$) remains a prominent goal.
- **RICs in Highly Structured or Deterministic Matrices**: Achieving RIP with near-optimal $m$ for deterministic constructions or those with significant algebraic structure is significantly less understood than for random matrices.
- **Beyond RIC: Alternative Metrics**: Given NP-hardness and extremality, complementary notions (mutual coherence, statistical RIP, etc.) are actively explored for practical verifiability.

Recent techniques leveraging advances in statistical mechanics (replica symmetric and replica symmetric breaking analysis), concentration of measure, chaos processes, and probabilistic order-statistics provide ongoing progress and increasingly tight upper and lower bounds for RICs in both random and structured settings [1501.06281, 1306.3779, 1410.1956].

---

References: [1308.0455], [1205.2081], [1201.2733], [1501.01708], [1604.03306], [1101.4458], [1312.3379], [1310.2410], [1309.6073], [1207.4883], [1410.1956], [1709.05789], [2409.08699], [1306.3779], [1501.06281], [1604.01171], [1207.0235], [1003.3299], [1802.03195].

Source: https://www.emergentmind.com/topics/restricted-isometry-constant-ric