Explicit constructions of RIP matrices and related problems
Jean Bourgain, S. J. Dilworth, Kevin Ford, Sergei Konyagin, Denka Kutzarova
Introduction
While most authors work with real signals and matrices, in this paper we work with complex matrices for convenience. Given a complex matrix satisfying (1.1), the real matrix , formed by replacing each element of by the matrix , also satisfies (1.1) with the same parameters .
Given , we wish to find RIP matrices of order with constant , and with as large as possible. If the entries of are independent Bernoulli random variables with values , then with high probability, will have the required properties forFor convenience, we utilize the Vinogradov notation , which means , and the Hardy notation , which means .
See ; also for a proof based on the Johnson-Lindenstrauss lemma . The first result of similar type for these matrices is due to Kashin . See also for RIP matrices with rows randomly selected from the rows of a discrete Fourier transform matrix and for other random constructions of RIP matrices. The parameter cannot be taken larger; in fact
It is an open problem to find good explicit constructions of RIP matrices; see T. Tao’s Weblog for a discussion of the problem. We mention here that all known explicit examples of RIP matrices are based on constructions of systems of unit vectors (the columns of the matrix) with small coherence.
Matrices whose columns are unit vectors with small coherence are connected to a number of well-known problems, a few of which we describe below. Systems of vectors with small coherence are also known as spherical codes. Some other applications of matrices with small coherence may be found in .
Suppose that are the columns of a matrix and have coherence . Then satisfies RIP of order with constant .
For any -sparse vector ,
All explicit constructions of matrices with small coherence are based on number theory. There are many constructions producing matrices with
In particular, such examples have been constructed by Kashin , Alon, Goldreich, Håstad and Peralta , DeVore , and Nelson and Temlyakov . By Proposition 1, these matrices satisfy RIP with constant and order
It follows from random constructions of Erdős and Rényi for Turán’s problem (see Proposition 2 and (1.15) below) that for any there are vectors with coherence
By contrast, there is a universal lower bound
valid for and all , due to Levenshtein (see also and ). Therefore, by estimating RIP parameters in terms of the coherence parameter we cannot construct RIP matrices of order larger than and constant .
Using methods of additive combinatorics, we construct RIP matrices of order with .
There is an effective constant and an explicit number such that for any positive integers and , there is an explicit RIP matrix of order with constant .
For application to sparse signal recovery, it is sufficient to take fixed , and one needs an upper bound on in terms of . By Theorem 1, for some , large and , we construct explicit RIP matrices with .
The proof of Theorem 1 uses a result on additive energy of sets (Corollary 2, Theorem 4), estimates for sizes of sumsets in product sets (Theorem 5), and bounds for exponential sums over products of sets possessing special additive structure (Lemma 10).
We now return to the problem of constructing matrices with small coherence. By (1.6), the bound (1.4) cannot be improved if , but there is a gap between bounds (1.6) and (1.4) when . For example, (1.4) is nontrivial only for . Of particular interest in coding theory is the range for fixed , where there have been some improvements made to (1.4). A construction obtained by concatenating algebraic-geometric codes with Hadamard codes (see e.g. [23, Corollary 3] and Section 3 of ) produces matrices with coherence
which is nontrivial for , and is better than (1.4) when . In the range , Ben-Aroya and Ta-Shma improved both (1.4) and (1.7) by constructing binary codes (vectors with entries ) with coherence
In this paper, we introduce very elementary constructions of matrices with coherence which matches (up to a factor) the bound (1.7). Our constructions, which are based on a method of Ajtai, Iwaniec, Komlós, Pintz and Szemerédi , have the added utility of applying to Turán’s power-sum problem and to the problem of finding thin sets with small Fourier coefficients. For the last two problems, our construction gives better estimates than existing explicit constructions in certain ranges of the parameters.
Roughly speaking, a set with small Fourier coefficients can be used to construct a set of numbers for Turán’s problem, and a set of numbers in Turán’s problem can be used to produce a matrix with small coherence. This is made precise below.
We next describe the problem of explicitly constructing thin sets with small Fourier coefficients. If is a positive integer and is a set (or multiset) of residues modulo , we let
Given , we wish to find a small set with also small.
Turán’s problem concerns the estimation of the function
where are positive integers. There is a vast literature related to Turán’s problem; see, e.g., , , (chapter 5), , .
If is a multiset of integers modulo and for , we see that
We also have the following easy connection between Turán’s problem and coherence.
Given any vector with for all , the coherence of the matrix with the columns
satisfies .
Combining (1.9) and Proposition 2, for any multiset of residues modulo , the vectors (1.10) satisfy
An application of Dirichlet’s approximation theorem shows that a set with must have . In , sets which are not much larger are explicitly constructed so that is small. Specifically, by [1, (1),(2)], for each primeA corresponding result when is composite is given in . there is a set with and
where is the integer so that the -th iterate of the logarithm of lies in . The proof uses an iterative procedure. By modifying this procedure, and truncating after two steps, we prove the following. To state our results, for brevity write
For sufficiently large prime and such that
a set of residues modulo can be explicitly constructed so that
The method from , if applied without modification (with two iterations of the basic lemma), produces a conclusion in Theorem 2 with
The bound on in Theorem 2 is better than (1.12) for .
Together, the construction for Theorem 2 and (1.9) give explicit sets for Turán’s problem. By further modifying the construction, we can do better.
For sufficiently large positive integer and such that
a multiset such that , can be explicitly constructed so that
To put Theorem 3 in context, we briefly review what is known about . P. Erdős and A. Rényi used probabilistic methods to prove an upper estimate
Using the character sum bound of Katz , J. Andersson gave explicit examples of sets which give
One can see that (1.16) supersedes (1.15) for . Also, combining (1.16) with Proposition 2 provides yet another construction of matrices with coherence satisfying (1.4). On the other hand, by (1.6) and Proposition 2, we have the lower estimate
By comparison, the constructions in Theorem 3 are better than (1.16) in the range , that is, throughout the range (1.14) (our constructions require to be prime, however).
The constructions in Theorem 3 also produce, by Proposition 2, explicit examples of matrices with coherence
which is close to the bound (1.7). By Proposition 1, these matrices satisfy RIP with constant and order
We prove Theorem 1 in Sections 2–6, Theorem 2 in Section 7 and Theorem 3 in Section 8.
Construction of the matrix in Theorem 1
We fix a large even number . A value of can be specified; it depends on the constant in an estimate from additive combinatorics (Proposition 3, Section 4). Also, the value can be reduced if one proves a better version of the Balog–Szemerédi–Gowers lemma (Lemma 6 below).
and the sets will be defined below. Notice that the matrix can be extended to a matrix by adding zero rows. Clearly, the matrices and have the same RIP parameters.
We notice that all elements of are at most , and
For , take to be the matrix formed by the first columns of , padded with rows of zeros.
In the next four sections, we show that has the required properties for Theorem 1. First, in Section 3, we show that in (1.1) we need only consider vectors whose components are 0 or 1 (emphflat vectors). We prove the following.
Let and be a positive integer. Assume that the coherence parameter of the matrix is . Also, assume that for some and any disjoint with we have
Then satisfies the RIP of order with constant .
and, for and ,
with some , where is the number of solutions of with each .
holds where .
The proof of Lemma 2 is quite involved, and will be handled in three subsequent sections. We next demonstate how Theorem 1 may be deduced from it.
We first prove (2.4) for the specific set defined in (2.1), provided that (and thus ). We have to show that for any distinct and any nonzero integers such that and the sum
All summands in the right-hand side of (2.6) but the first one are divisible by . For the first summand we have
This shows that . Therefore, . By assumption, , and
Condition (2.5) is satisfied due to Corollary 4 of Section 5 with . If then Lemma 2 gives a nontrivial estimate with . Thus, satisfies the conditions of Corollary 1 with and (using for large , which follows from the prime number theorem). Let . Let , and let be the matrix formed by taking the first columns of , then adding rows of zeros. Clearly, satisfies the conditions of Corollary 1 with the same parameters as . By Lemma 1 with , Theorem 1 follows.
In Section 4 we introduce some notation and recall standard estimates in additive combinatorics, which will be applied to subsets of . Section 5 is devoted to the sumset theory of , from which we deduce (2.5). The completion of the proof of Lemma 2 is in Section 6. We give some preliminaries here.
where the summands with are excluded from the summation. We next break into balanced sets. For and , let
whenever are powers of two and, for and for any ,
Indeed, there are choices for . To prove the cancellation in (2.8), we basically split into two cases: (i) some has additive structure (that is, is large), where the cancellation comes from the sum over (with fixed), and (ii) when does not have additive structure, in which case one gets dispersion of the phases from the dilation weights (taking a large moment and using (2.4)). Incidentally, oscillations of the factor play no role in the argument.
The Flat-RIP property
Let be the columns of an matrix . Suppose that for every , . We say that satisfies the flat RIP of order with constant if for any disjoint with we have
For technical reasons, it is more convenient to work with the flat-RIP than with the RIP. However, flat-RIP implies RIP with an increase in . The flat-RIP property is closely related to the property that (1.1) holds for any with entries which are zero or one and at most ones (see the calculation at the end of this section).
Let and be a positive integer. Suppose that satisfies flat-RIP of order with constant . Then satisfies RIP of order with constant .
First, by a convexity-type argument and our assumption,
provided that , for all . Next, suppose , and for all . Without loss of generality assume that , where denotes the norm. For a positive integer let
Applying (3.2) to sets , we get
Let . By the Cauchy–Schwarz inequality we infer that
For the next step, suppose take arbitrary complex values, and . We partition and into subsets of cardinality at most each: , Next, for any we have
where are non-negative. By (3.4) and the Cauchy–Schwarz inequality,
For any disjoint with we have
Using the assumptions of the Lemma 1 directly rather than reducing it to Lemma 3, one can get a better constant for RIP; However, we do not need a stronger version of the corollary for our purposes.
Some definitions and results from additive combinatorics
For an (additive) abelian group we define the sum and the difference of subsets :
We will use the following lemma which is a particular case of Plünecke – Ruzsa estimates (, Exercise 6.5.15).
For any nonempty set we have .
If , we define the (additive) energy of the sets and as the number of solutions of the equation
Next, let . The -restricted sum of and is defined as
Trivially If is close to then must have a special additive structure.
(, Lemma 2.30) If then there exists such that and .
The following lemma is a version of the Balog–Szemerédi–Gowers lemma which plays a very important role in additive combinatorics.
If , and . Then there exists a set such that and .
Combining Lemma 5 and Lemma 6 gives the following.
If then there exists a set such that and .
By we denote the indicator function of the set . With this notation, we have
An explicit version of Proposition 3, with , is given in .
Note that if , we may decompose as a disjoint union of at most sets with and apply (4.2) for each . Hence
Applying the Cauchy–Schwarz inequality we get
It would be interesting to find best possible value for in Proposition 3. The example shows that .
Put , and let be a permutation of such that . By (4.3), for we have , where
Denote . Notice that since . Separately considering and and using the Cauchy–Schwarz inequality, we get
Using a parameter which will be specified later we define the sets
Decompose where
The contribution to the sum in the theorem from and is negligible. First,
Using Young’s inequality (cf , Theorem 4.8), we find that
So, it suffices to estimate the contribution of . We have
Hence, . Now we can use Corollary 2:
Combining the last inequality with (4.6) – (4.8) we get
Taking completes the proof of the theorem. ∎
A sumset estimate in product sets
The main result of this section is the following.
Then for any subsets we have
Observe that for we have where
By Theorem 5, . On the other hand, . If then
So, the asymptotic behavior of as is sharp. Likely, inequality (5.1) holds with . This was proved in the case by Woodall .
For positive integers we define an path as a sequence of pairs of integers such that for any either , or .
Let , , , . Then there exists an path such that
We proceed by induction on . For or the assertion is obvious. We prove it for with , supposing that it holds for replaced by and . Without loss of generality we assume that
By the induction supposition, there exists an path such that and
Similarly, Thus, where
The function has negative third derivative on $f(0)=f(1/M)=f(1)=0ff(u)>0uf(x)\geq 01/M\leq x\leq 1f(w)\geq 0$ as desired. ∎
We will need Lemma 7 only for (although for the proof it was convenient to have varying ).
Let , be non-negative numbers, and . Then
Lemma 8 has some similarity with inequality (2.1) from .
We order and in the descending order and , respectively, where for some permutations and of the set we have . We consider an arbitrary path with . Since and ,
Consequently, there is a permutation of so that
Thus, for some and we have
But for some . Recalling that and we obtain . Similarly, . Therefore,
Now we are ready to prove Theorem 5. We proceed by induction on . For the set is a singleton, and there is nothing to prove. Now suppose that the assertion holds for replaced by . We consider arbitrary subsets . For we denote
Let . For we denote
By the induction supposition, . Hence,
The set is a translate of some set , and is Freiman isomorphic to . Hence, for any we have . If then . By (5.2) and a short calculation using , . ∎
Let . By Corollary 1, there is a set such that and . If and is so large that then we get contradiction with Corollary 3. ∎
The proof of Lemma 2
We may assume , otherwise there is nothing to prove. Adopt the notation () from Section 2. If , then by (2.9), and (2.8) holds (recall that , hence ). Thus, we can assume that , which implies, by (2.3), that
Let denote the double sum over . By the Cauchy–Schwarz inequality,
Another application of the Cauchy–Schwarz inequality gives
A third application of the Cauchy–Schwarz inequality, followed by Parseval’s identity yields a well-known inequality (cf. , Problem 14(a) for Chapter 6)
By (6.1), , and by Lemma 9 and (2.5),
Thus, if and , then and (2.8) follows. Otherwise, without loss of generality we may assume that
The following lemma gives the necessary estimates to complete the proof of Lemma 2. For , set
The proof of Lemma 10 applies to more general sums, e.g. in one may replace the Legendre symbol with arbitrary complex numbers with modulus , and one may replace with different quantities having the dissociative property (the analog of (2.4) holds).
Postponing the proof of Lemma 10, we show first how to deduce Lemma 2.
We take a maximal subset so that (6.5) holds for . Denote . By Lemma 9, (2.9), and (2.3) we have
Now assume that (6.6) does not hold. By (2.9), we get
Applying now Corollary 1 and (2.9) we obtain the existence of a set such that
and . Using Lemma 10 we get inequality (6.5) for . Therefore, (6.5) is also satisfied for , contradicting the choice of .
Thus, we have shown that (6.6) must hold. Using (6.5) for and (6.7) we get
Summing on and using (2.3) and (2.9), we obtain
Hence, for some complex numbers of modulus ,
Then . By Hölder’s inequality,
As , we have by the triangle inequality,
Define the probability measure by
with . By (2.4), this has only trivial solutions and thus
On the other hand, it follows from (6.3) and (6.4) that
Subsequent application of (6.9), (6.10) and (6.11) gives
By (6.3), . Recalling , (2.9), (6.2) and (6.4), we conclude that
Plugging the last estimate into (6.8), we get
Thin sets with small Fourier coefficients
Denote by the inverse of modulo . It is easy to see for relatively prime integers that
Let , , and be a positive integer. Suppose that for every prime , is a set of integers in . Suppose is a prime satisfying . Then the numbers , where , are distinct modulo .
Multiplying both sides by gives
The right side is divisible by and the absolute value of the right side is , hence both sides are zero, , and . ∎
For brevity, we write for is what follows.
Let , , and be a positive integer. Suppose that for every prime , is a multiset of integers in , and . Suppose is a prime satisfying . Then the multiset
where is the number of primes in .
Since , we may assume without loss of generality that . We have
If , we use the trivial bound and conclude
Now assume . If , then . When , by (7.1),
Since there are primes with , we have
Combining our estimates for and , we arrive at
For a specific choice of , the inequality (7.2) can be strengthened.
Let and be a positive integer. For every prime denote by the set of all integers in . Suppose is a prime satisfying . Then the multiset
where is the number of primes in and .
Again, we may assume without loss of generality that . We use notation from the proof of Lemma 12. If , we use the trivial estimate . Now there are primes with . When , by (7.1),
where it is assumed that . For we denote
Taking into account that for we get
If then is divisible by . But . Therefore, the number of prime divisors of any number is at most and for any we get
Combining our estimates for and ((7.3) and (7.5)), we arrive at
Applying Lemma 12 for all primes in a dyadic interval, we can then feed these multisets back into the lemma and iterate.
Using explicit estimates for counts of prime numbers , we have
For , there are more than primes in . For any , there are at most primes in .
Using Proposition 4 we obtain a more convenient version of Lemma 13.
Let . For every prime denote by the set of all nonzero integers in . Suppose is a prime satisfying and suppose is a positive integer. Then the multiset
We use the notation of Lemma 13. By Proposition 4 we have
On the other hand, using Proposition 4 again we get
Now the inequality (7.6) follows from (7.7) and (7.4). ∎
Using just one iteration one can get the following effective result on thin sets with small Fourier coefficients, of nearly the same strength as (1.12).
For sufficiently large prime and such that there is a set of residues modulo so that
Clearly, . Let be the multiset constructed in Lemma 14. We have . By Lemma 11, is a set. Moreover,
We choose real parameters , and positive integers , so that
For , let be the set of integers in . By Lemmas 11, 14 and (7.8), for each prime , there is a set of residues modulo such that
By an application of Lemmas 11 and 12 with , , , and , together with (7.9), there is a set of residues modulo so that
so that (7.9) follows immediately. The condition (1.13) implies (7.8) for large enough . ∎
Theorem 2 supersedes Corollary 5 for .
An explicit construction for Turán’s problem
We follow the proof of Theorem 2 and Lemma 12. We choose real parameters , and a positive integer , so that
For , let be the set of integers in . By Lemma 14 and (8.1), for each prime , there is a multiset of residues modulo such that
We have for all , where . Now define a multiset as a union of multisets . We have, for ,
If , then . When , by (8.3), . Therefore,
The sum over is estimated at the same way as in Lemma 12:
Combining (8.4), (8.5) and using Proposition 4 we arrive at
as required. Moreover, by Proposition 4 we have
Now we take the same as in the proof of Theorem 2 so that (8.2) follows immediately. The condition (1.14) implies (8.1) for large enough . ∎
As in , one can construct thin sets modulo with and small, by iterating Lemma 12. Roughly speaking, applying Lemma 14 followed by iterations of Lemma 12 produces sets , with small , as small as , where is the -th iterate of the logarithm of . We omit the details.
Acknowledgments. The authors thank Ronald DeVore, Zeev Dvir, Venkatesan Guruswami, Piotr Indyk, Sina Jafarpour, Boris Kashin, Howard Karloff, Imre Leader, Igor Shparlinski and Avi Wigderson for helpful conversations.