Neighborly Polytopes And Sparse Solution Of Underdetermined Linear Equations
Neighborly Polytopes And Sparse Solution Of Underdetermined Linear Equations
复制标题
DOI:
--
复制
发表时间:
2005
期刊:
影响因子:
--
通讯作者:
D. Donoho
中科院分区:
文献类型:
--
作者:
D. Donoho
Consider a d × n matrix A, with d < n. The problem of solving for x in y = Ax is underdetermined, and has many possible solutions (if there are any). In several fields it is of interest to find the sparsest solution – the one with fewest nonzeros – but in general this involves combinatorial optimization. Let ai denote the i-th column of A, 1 ≤ i ≤ n. Associate to A the quotient polytope P formed by taking the convex hull of the 2n points (±ai) in R. P is centrosymmetric and is called (centrally) k-neighborly if every subset of k + 1 elements (±ilail) k+1 l=1 are the vertices of a face of P . We show that if P is k-neighborly, then if a system y = Ax has a solution with at most k nonzeros, that solution is also the unique solution of the convex optimization problem min ‖x‖1 subject to y = Ax; the converse holds as well. This complete equivalence between the study of sparse solution by ` minimization and neighborliness of convex polytopes immediately gives new results in each field. On the one hand, we get new families of neighborly centrosymmetric polytopes, by exploiting known results about sparsity of ` minimization; on the other, we get new limits on the ability of ` minimization to find sparse solutions, by exploiting known limits on neighborliness of centrally symmetric polytopes. Weaker notions of equivalence between ` and sparse optimization have also been studied recently. These are equivalent to other interesting properties of the quotient polytope. Thus, suppose the columns of A are in general position. Consider the vectors having k < d/2 nonzeros that are simultaneously the sparsest solution of y = Ax and the minimal ` solution. These make up a fraction 1 − of all vectors with k nonzeros if and only if the quotient polytope P has at least 1− as many k-dimensional faces as the regular cross polytope C. Combining this with recent work on face numbers of randomly-projected cross-polytopes, we learn that for large d, the overwhelming majority of systems of linear equations with d equations and 4d/3 unknowns have the following property: if there is a solution with fewer than .49d nonzeros, it is the unique minimum ` solution. A stylized application in digital communication is sketched; for large n, it is possible to transmit n/4 pieces of information using a codeword of length n with immunity to .49n gross errors in the received codeword, if the signs and sites of the gross errors are random, and with immunity to .11n gross errors chosen by a malicious opponent. The receiver uses ` minimization.