Donoho David L, Tanner Jared
Department of Statistics, Stanford University, Stanford, CA 94305-4065, USA.
Proc Natl Acad Sci U S A. 2005 Jul 5;102(27):9452-7. doi: 10.1073/pnas.0502258102. Epub 2005 Jun 22.
Let A be a d x n matrix and T = T(n-1) be the standard simplex in Rn. Suppose that d and n are both large and comparable: d approximately deltan, delta in (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 Rn. We derive rhoN(delta) > 0 with the property that, for any rho < rhoN(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 < or = k < or = rhod. This implies that P is left floor rhod right floor-neighborly, and its skeleton Skel(left floor rhod right floor)(P) is combinatorially equivalent to Skel( left floor rhod right floor)(T). We also study a weaker notion of neighborliness where the numbers of k-dimensional faces f(k)(P) > or = f(k)(T)(1-epsilon). Vershik and Sporyshev previously showed existence of a threshold rhoVS(delta) > 0 at which phase transition occurs in k/d. We compute and display rhoVS and compare with rhoN. Corollaries are as follows. (1) The convex hull of n Gaussian samples in Rd, with n large and proportional to d, has the same k-skeleton as the (n-1) simplex, for k < rhoN (d/n)d(1 + oP(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 rhoVS(d/n)d(1 + o(1)) nonzeros, linear programming will find that solution.
设(A)为一个(d\times n)矩阵,(T = T(n - 1))为(\mathbb{R}^n)中的标准单纯形。假设(d)和(n)都很大且具有可比性:(d\approx\delta n),其中(\delta\in(0,1))。当从(\mathbb{R}^n)的(d)维正交投影算子的格拉斯曼流形中随机均匀地选择投影算子(A)时,我们对投影后的单纯形(AT)的面进行计数。我们推导出(\rho_N(\delta)>0),其性质为:对于任何(\rho<\rho_N(\delta)),当(d)很大时,以压倒性概率,对于(0\leq k\leq\rho d),(P = AT)的(k)维面的数量与(T)的完全相同。这意味着(P)是下取整(\rho d)邻接的,并且其骨架(Skel(\lfloor\rho d\rfloor)(P))与(Skel(\lfloor\rho d\rfloor)(T))组合等价。我们还研究了一种较弱的邻接概念,其中(k)维面的数量(f(k)(P)\geq f(k)(T)(1 - \epsilon))。韦尔希克和斯波里舍夫先前证明了存在一个阈值(\rho_{VS}(\delta)>0),在该阈值处(k/d)会发生相变。我们计算并展示(\rho_{VS})并与(\rho_N)进行比较。推论如下:(1) 在(\mathbb{R}^d)中(n)个高斯样本的凸包,当(n)很大且与(d)成比例时,对于(k < \rho_N(d/n)d(1 + o_P(1))),其(k)骨架与((n - 1))维单纯形相同。(2) 在求解欠定线性方程组的最稀疏非负解时,线性规划的能力存在“相变”。对于大多数具有少于(\rho_{VS}(d/n)d(1 + o(1)))个非零解的系统,线性规划将找到该解。