arXiv:1207.3094v2 [cs.IT] 15 Jul 2013 eehleTrust. Leverhulme ( UK | oyehiu ´deaed asne(PL,Luan,Sw Lausanne, (EPFL), Lausanne de ( F´ed´erale Polytechnique ie nDfiiin11adterilsrto nFg 1. Fig. in i illustration graph their expander and lossless 1.1 a Definition of in given definition adjacen precise sparse a very with matrices, Expand graphs connected graphs. highly expander are lossless graphs of matrices adjacency are here). considered (as and column Gaussian per dense nonzeros fixe small of to very well number with compared even how matrices do Transform Cosine implementation matrices Discrete of GPU [ type a Tanner these and in Blanchard demonstrated fact, In recently matrices counterparts. sensing dense matric good as their sparse be as can shows graphs that expander from results coming numerical guarante exhibiting [4 performance with also theoretical both sensing, giving made [7] compressed [6], been in [5], has matrices, progress matrices dense significant sparse late, to incorporate Of compared [3]. [2], as [1], implementation see comple fast storage and low their ity of because mathematics putational n,sas arcs xadrgraphs. expander matrices, sparse ing, esn apigterm hnuigsas o mean-zero non sparse using spars compress when the matrices. quantitative measurement theorems hence also and sampling theor nove and graphs consider quantitative sensing a we expander for is matrices lossless allow analysis random of that this existence analysis constants to on order detailed Key better sets. more of of a splitting unions through dyadic in derived collisions th are of for expe bounds bounds the similar These than are less bounds the be these expansion from will of Deducible cardinality cardinality value. this expected that the probability for neighbors gra expander of formulae degree left of present corresponde fixed of number We one-to-one matrices a adjacency fixed the have wit a with random matrices at has uniformly These drawn replacement. column are indices each row whose nonzeros where matrices random xadri ti iatt rp with graph bipartite a is it if expander [email protected] V a ihnl preMtie n xadrGraphs, Expander and Matrices Sparse Vanishingly oyih c 02IEEE 2012 (c) Copyright ahmtc nttt n xtrClee nvriyo Ox of University College, Exeter and Institute Mathematics aoaoyfrIfrainadIfrneSses(LIONS), Systems Inference and Information for Laboratory ento 1.1: Definition that matrices sparse random consider we manuscript this In com- and applied in useful particularly are matrices Sparse ne Terms Index Abstract | = [email protected] n ih etcsadhsarglrlf degree left regular a has and vertices right ftegahwihi fitrs nmn applications. many in interest of is which graph the of W eii h rbblsi osrcino sparse of construction probabilistic the revisit —We o hs rps n rsn albud nthe on bounds tail present and graphs, these for ihApiaint opesdSensing Compressed to Application With Agrtm,cmrse esn,sga process- signal sensing, compressed —Algorithms, ehiu.Teaayi e otederivation the to led analysis The technique. G .I I. ( = ). NTRODUCTION ,V E V, U, .JT cnwegsspotfo the from support acknowledges J.T. ). ) salossless a is | U | uaarBh n ae Tanner Jared and Bah, Bubacarr = N etvertices, left od Oxford, ford, ( ,d ǫ d, k, d itzerland such , sand es Ecole ´ cted phs. ems nce set cy ed 8] es x- to er ) ], d h e s e - l neighbors. any that i.1 lutaino lossless a of Illustration 1. Fig. d ntefr fabudo h alpoaiiyo h ieof size the of probability tail the on bound obje a these of of form construction the probabilistic in the on guarantees [12]. achiev al. et. only Guruswami constructions see parameters, Deterministic sub-optimal here. consider optim we with constructions probabil Probabilistic high parameters graph. with is, expander graph an a bipartite [11] that Bassylago showed left-regular and and random expanders Pinsker lossless survey. of detailed existence [10] the more proved see Storag a Complexity, for Proof Distributed [9] in or a Tautologies Hard and and Time of Method, Fault-tolerance Linear Complexity Subsets, Bitprobe Networks, Storing applications Codes, in many Error-Correcting have Routing Decodable and Distributed mathematics including: pure and science osesexpanders lossless conductors 2 = 1 )Te r called are They 2) )The 3) eak1.2: Remark u ancnrbto,i h rsnaino quantitative of presentation the is contribution, main Our computer theoretical in studied well been have graphs Such osesexpanders Lossless . (1 − X ihprmtr htaebs oaihso h parameters the of logarithms 2 base are that parameters with expansion ǫ 1 ,N n, ) ⊂ d . e 9,[]adterfrne therein. references the and [5] [9], see U xs u r o utbefrteapplications the for suitable not are but exist )Tegah are graphs The 1) with ihparameters with falossless a of naacdexpanders unbalanced | X ≤ | ( ,d ǫ d, k, k ,k ,N n, k, d, ) has epne rpswith graphs -expander ( ,d ǫ d, k, lossless | Γ( r qiaetto equivalent are X ) epne rp is graph -expander ) | when because (1 = n − k ≪ ǫ 4 = ǫ ) lossless ≪ d N | ity, X and cts ; ny 1 al of e e ; 1 | 2 the set of neighbors, Γ(X) for a given X U, of a randomly theory and linear algebra we define the set of neighbors in generated left-degree bipartite graph. Moreover,⊂ we provide both notation. deducible bounds on the tail probability of the expansion of Definition 1.4: Consider a bipartite graph G = (U,V,E) the graph, Γ(X) / X . We derive quantitative guarantees for where E is the set of edges and e = (x ,y ) is the edge that | | | | ij i j randomly generated non-mean zero sparse binary matrices to connects vertex x to vertex y . For a set of left vertices S U i j ⊂ be adjacency matrices of expander graphs. In addition, we its set of neighbors is Γ(S) = yj xi S and eij E . In derive the first phase transitions showing regions in parameter terms of the adjacency matrix, {A,| of G∈ = (U,V,E∈) the} set space that depict when a left-regular bipartite graph with a of neighbors of AS for S = s, denoted by As, is the set of given set of parameters is guaranteed to be a lossless expander rows with at least one nonzero.| | with high probability. The most significant contribution, which Definition 1.5: Using Definition 1.4 the expansion of the is the key innovation, of this paper is the use of a novel graph is given by the ratio Γ(S) / S , or equivalently, As /s. technique of dyadic splitting of sets. We derived our bounds By the definition of a lossless| | expander,| | Definition 1.1,| | we using this technique and apply them to derive ℓ1 restricted need Γ(S) to be large for every small S U. In terms of | | ⊂ isometry constants (RIC1). the class of matrices defined by Definition 1.3, for every AS Numerous compressed sensing algorithms have been de- we want to have As as close to n as possible, where n is signed for sparse matrices [4], [5], [6], [7]. Another contribu- the number of rows.| Henceforth,| we will only use the linear tion of our work is the derivation of sampling theorems, pre- algebra notation As which is equivalent to Γ(S). Note that sented as phase transitions, comparing performance guarantees A is a random variable depending on the draw of the set of | s| for some of these algorithms as well as the more traditional ℓ1 columns, S, for each fixed A. Therefore, we can ask what is minimization compressed sensing formulation. We also show the probability that A is not greater than a , in particular | s| s how favorably ℓ1 minimization performance guarantees for where a is smaller than the expected value of A . This is s | s| such sparse matrices compared to what ℓ2 restricted isome- the question that Theorem 1.6 to answers. We then use this try constants (RIC2) analysis yields for the dense Gaussian theorem with RIC1 to deduce the corollaries that follow which matrices. For this comparison, we used sampling theorems are about the probabilistic construction of expander graphs, the and phase transitions from related work by Blanchard et. al. matrices we consider, and sampling theorems of some selected [13] that provided such theorems for dense Gaussian matrices compressed sensing algorithms. based on RIC2 analysis. Theorem 1.6: For fixed s,n,N and d, let an n N matrix, The outline of the rest of this introduction section goes as A be drawn from either of the classes of matrices× defined in follows. In Section I-A we present our main results in Theorem Definition 1.3, then 1.6 and Corollary 1.7. In Section I-B we discuss RIC1 and its implication for compressed sensing, leading to two sampling Prob ( A a )
where pmax(s, d) is given by A. Main results 2 Our main results is about a class of sparse matrices coming pmax(s, d)= , and (2) 25√2πs3d3 from lossless expander graphs, a class which include non-mean zero matrices. We start by defining the class of matrices we consider and a key concept of a set of neighbors used in the s/2 1 ⌈ ⌉ s a2i ai derivation of the main results of the manuscript. Firstly, we Ψ (as,...,a1)= (n ai) H − n 2i − · n a denote H(p) := p log(p) (1 p) log(1 p) as the Shannon " i=1 − i − − − − X entropy function of base e logarithm. a a a +a H 2i i n H i +3s log (5d) (3) Definition 1.3: Let be an matrix with nonzeros i − A n N d · ai − · n # in each column. We refer to A as× a random where . 1) sparse expander (SE) if every nonzero has value 1 a1 := d 2) sparse signed expander (SSE) if every nonzero has value If no restriction is imposed on as then the ai for i> 1 take from 1, 1 on their expected value ai given by {− } and the support set of the d nonzeros per column are drawn ai a2i = ai 2 b for i =1, 2, 4,..., s/2 . (4) uniformly at random, with each column drawn independently. − n ⌈ ⌉ SE matrices are adjacency matrices of lossless (k,d,ǫ)- b If a is restricted to be less than a , then the a for i > 1 expander graphs while SSE matrices have random sign pat- bs b s i are the unique solutions to the following polynomial system terns in the nonzeros of an adjacency matrix of a lossless (k,d,ǫ)-expander graph. If A is either an SE or SSE it will 3 2 2 2 b a2i 2aia2i +2ai a2i ai a4i =0 for i =1, 2,..., s/4 (5) have only d nonzeros per column and since we fix d n, A is − − ⌈ ⌉ therefore “vanishingly sparse.” We denote A as a≪ submatrix with a a for each i. S 2i ≥ i of A composed of columns of A indexed by the set S with Corollary 1.7: For fixed s,n,N,d and 0 <ǫ< 1/2, let an S = s. To aid translation between the terminology of graph n N matrix, A be drawn from the class of matrices defined | | × 3 in Definition 1.3, then Theorem 1.9: If an n N matrix A is either SE or SSE defined in Definition 1.3,× then A/d satisfies RIP-1 with Prob ( ASx 1 (1 2ǫ)d x 1)