Neighborliness of randomly projected simplices in high dimensions

被引:207
作者
Donoho, DL [1 ]
Tanner, J [1 ]
机构
[1] Stanford Univ, Dept Stat, Stanford, CA 94305 USA
关键词
neighborly polytopes; convex hull of Gaussian sample; underdetermined systems of linear equations; uniformly distributed random projections; phase transitions;
D O I
10.1073/pnas.0502258102
中图分类号
O [数理科学和化学]; P [天文学、地球科学]; Q [生物科学]; N [自然科学总论];
学科分类号
07 ; 0710 ; 09 ;
摘要
Let A be a d x n matrix and T=Tn-1 be the standard simplex in R-n. Suppose that d and n are both large and comparable: d approximate to delta n, delta epsilon (0, 1). We count the faces of the projected simplex AT when the projector A is chosen uniformly at random from the Grassmann manifold of d-dimensional orthoprojectors of R-n. We derive rho(N)(delta) > 0 with the property that, for any rho < rho(N)(delta), with overwhelming probability for large d, the number of k-dimensional faces of P = AT is exactly the same as for T, for 0 <= k <= rho d. This implies that P is vertical bar(-)rho d(-)vertical bar-neighborly, and its skeleton Skel vertical bar(-rho d-)vertical bar(P) is combinatorially equivalent to Skel vertical bar(-rho d-)vertical bar(T). We also study a weaker notion of neighborliness where the numbers of k-dimensional faces f(k)(P) >= f(k)(T)(1-epsilon). Vershik and sporyshev previously showed existence of a threshold rho vs(delta) > 0 at which phase transition occurs in k/d. We compute and display rho vs and compare with rho(N). Corollaries are as follows. (1) The convex hull of n Gaussian samples in R-d, with n large and proportional to d, has the same k-skeleton as the (n-1) simplex, for k < rho(N) (d/n)d(1 + o(p)(1)). (2) There is a "phase transition" in the ability of linear programming to find the sparsest nonnegative solution to systems of underdetermined linear equations. For most systems having a solution with fewer than rho vs(d/n)d(1 + o(1)) nonzeros, linear programming will find that solution.
引用
收藏
页码:9452 / 9457
页数:6
相关论文
共 14 条
[1]   RANDOM PROJECTIONS OF REGULAR SIMPLICES [J].
AFFENTRANGER, F ;
SCHNEIDER, R .
DISCRETE & COMPUTATIONAL GEOMETRY, 1992, 7 (03) :219-226
[2]  
[Anonymous], 2003, GRADUATE TEXTS MATH
[3]   REGULAR SIMPLICES AND GAUSSIAN SAMPLES [J].
BARYSHNIKOV, YM ;
VITALE, RA .
DISCRETE & COMPUTATIONAL GEOMETRY, 1994, 11 (02) :141-147
[4]  
Boroczky K, 1999, ARCH MATH, V73, P465
[5]  
DEMBO A, APPL MATH, V38
[6]   Sparse nonnegative solution of underdetermined linear equations by linear programming [J].
Donoho, DL ;
Tanner, J .
PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA, 2005, 102 (27) :9446-9451
[7]  
DONOHO DL, 2005, 200504 STANF U DEP S
[8]  
DONOHO DL, 2005, 200505 STANF U DEP S
[9]  
DONOHO DL, 200506 STANF U DEP S
[10]  
DONOHO DL, 2004, 200410 STANF U DEP S