Theorem 1.4: minimization recovers every -sparse vector when
ProvedCandesTao.Decoding.l1_recovers_sparse_vectorLet be a real matrix with columns , and suppose that is such that
Let be a real vector supported on a set obeying , and put . Then is the unique minimizer of
This is the paper's main result: under a deterministic "restricted orthonormality" hypothesis on , far weaker than orthonormality of its columns and compatible with much larger than , the convex program , which can be recast as a linear program, recovers every sufficiently sparse vector exactly and with no probability of failure. By Lemma 1.2, condition (1.10) implies (the hypothesis of Lemma 1.3) and is implied by .
Formalization Note "Unique minimizer" means that is feasible and every other feasible has strictly larger norm. The hypothesis is the domain on which Definition 1.1 defines .
import Definitions.Def_CandesTao_Decoding_RestrictedIsometry import Definitions.Def_CandesTao_Decoding_L1Minimization
namespace CandesTao.Decoding
theorem l1_recovers_sparse_vector {p m : ℕ} (F : Matrix (Fin p) (Fin m) ℝ) (S : ℕ)
(hS : 1 ≤ S) (hSm : 3 * S ≤ m)
(h : restrictedIsometryConst F S + restrictedOrthogonalityConst F S S +
restrictedOrthogonalityConst F S (2 * S) < 1)
(T : Finset (Fin m)) (c : Fin m → ℝ) (hT : T.card ≤ S) (hc : SupportedOn c T) :
IsUniqueL1Minimizer F (F.mulVec c) c := by sorry
end CandesTao.DecodingRead-back
What the Lean code literally says, in plain math · claude-fable-5-1
Read-back of l1_recovers_sparse_vector
Setting and notation. Throughout, and are natural numbers (a priori either may be ; see the edge cases below), and are the index sets, and is a real matrix. For , is the ordinary matrix–vector product , and for , is the ordinary dot product. Sets of indices are always finite, and is their cardinality. The following notions are the bundle's own definitions, unfolded here.
-
Support. For and , " is supported on " means: for every with , . So ; may vanish on part or all of , and forces . A vector supported on some with is therefore a vector with at most nonzero entries.
-
norm. .
-
norm. , using Mathlib's real square root (which returns on negative inputs; here the argument is a sum of squares, so this is the usual Euclidean norm). Hence . When the sum defining is empty, so for every .
-
Restricted isometry constant . For a natural number ,
Here ranges over all subsets of of size at most (including , for which only qualifies and both inequalities are trivial), and both inequalities are non-strict.
- Restricted orthogonality constant . For natural numbers ,
The bound on the right is times the product of the two (unsquared) Euclidean norms; and may be empty.
- Unique minimizer. For and , " is the unique minimizer for " means
The inequality is strict, so this says that minimizes over the affine set and that no other point of that set attains the same value.
About the two infima. In Mathlib, the infimum of a set is the genuine greatest lower bound when is nonempty and bounded below, and is defined to be when is empty or not bounded below. Both defining sets above consist of nonnegative reals, so they are bounded below by , and both are nonempty: by Cauchy–Schwarz, for every , so lies in the first set and lies in the second. Consequently and are honest infima and both are . Each set is moreover upward closed (if qualifies, so does every ; likewise for ) and cut out by non-strict inequalities, so each is a closed half-line and its infimum is attained: is the least for which the two-sided inequality holds for all vectors with at most nonzero entries, and is the least for which the dot-product bound holds for all disjointly supported pairs of the stated sparsities.
The statement. Let , let be a real matrix, and let . Assume:
- ;
- (ordinary multiplication and comparison of natural numbers);
where concerns disjoint index sets of size at most each, and concerns disjoint index sets of size at most and at most respectively; 4. is a set of indices with ; 5. is supported on , i.e. for every .
Then is the unique minimizer for . Unfolded, and noting that the first conjunct holds trivially, the conclusion is exactly
Here is exact equality of two vectors in , and ranges over all of with no sparsity restriction.
Edge cases and what the quantifiers silently include.
- Hypotheses 1 and 2 together force (so is excluded), , and , so disjoint index sets of the sizes named in and do exist inside . Nothing is assumed about relative to or .
- Since all three constants are , hypothesis 3 implies that each of , , is individually .
- If some nonzero vector with at most nonzero entries satisfies , then every in the defining set satisfies , hence , so and hypothesis 3 cannot hold: the theorem is vacuous for such . In particular this happens when ( has no rows, for every , and a nonzero -sparse vector exists because and ), so the hypotheses implicitly force .
- may be empty or have fewer than elements, and may vanish on some or all of ; in particular is allowed, in which case the conclusion reads "every with has ", which is automatic.
- The statement concerns one fixed, deterministic matrix : there is no probability, no normalization of the columns of , and no relation such as is imposed.
- The bundle also defines the -th column of and the linear span of all columns of ; neither is referenced by this statement.
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.