Counting faces of randomly projected polytopes when the projection radically lowers dimension
Top Cited Papers
Open Access
- 10 July 2008
- journal article
- Published by American Mathematical Society (AMS) in Journal of the American Mathematical Society
- Vol. 22 (1) , 1-53
- https://doi.org/10.1090/s0894-0347-08-00600-0
Abstract
Let Q = Q N Q = Q_N denote either the N N -dimensional cross-polytope C N C^N or the N − 1 N-1 -dimensional simplex T N − 1 T^{N-1} . Let A = A n , N A = A_{n,N} denote a random orthogonal projector A : R N ↦ b R n A: \mathbf {R}^{N} \mapsto bR^n . We compare the number of faces f k ( A Q ) f_k(AQ) of the projected polytope A Q AQ to the number of faces of f k ( Q ) f_k(Q) of the original polytope Q Q . We concentrate on the case where n n and N N are both large, but n n is much smaller than N N ; in this case the projection dramatically lowers dimension. We consider sequences of triples ( k , n , N ) (k,n,N) where N = N n N = N_n is not exponentially larger than n n . We identify thresholds of the form c o n s t ⋅ n log ( n / N ) const \cdot n \log (n/N) where the relationship of f k ( A Q ) f_k(AQ) and f k ( Q ) f_k(Q) changes abruptly. These properties of polytopes have significant implications for neighborliness of convex hulls of Gaussian point clouds, for efficient sparse solution of underdetermined linear systems, for efficient decoding of random error correcting codes and for determining the allowable rate of undersampling in the theory of compressed sensing. The thresholds are characterized precisely using tools from polytope theory, convex integral geometry, and large deviations. Asymptotics developed for these thresholds yield the following, for fixed ϵ > 0 \epsilon > 0 . With probability tending to 1 as n n , N N tend to infinity: (1a) for k > ( 1 − ϵ ) ⋅ n [ 2 e ln ( N / n ) ] − 1 k > (1-\epsilon ) \cdot n [2e\ln (N/n)]^{-1} we have f k ( A Q ) = f k ( Q ) f_k(AQ) = f_k(Q) , (1b) for k > ( 1 + ϵ ) ⋅ n [ 2 e ln ( N / n ) ] − 1 k > (1 +\epsilon ) \cdot n [2e\ln (N/n)]^{-1} we have f k ( A Q ) > f k ( Q ) f_k(AQ) > f_k(Q) , with E {\mathcal E} denoting expectation, (2a) for k > ( 1 − ϵ ) ⋅ n [ 2 ln ( N / n ) ] − 1 k > (1-\epsilon ) \cdot n [2\ln (N/n)]^{-1} we have E f k ( A Q ) > ( 1 − ϵ ) f k ( Q ) {\mathcal E} f_k(AQ) > (1-\epsilon ) f_k(Q) , (2b) for k > ( 1 + ϵ ) ⋅ n [ 2 ln ( N / n ) ] − 1 k > (1 +\epsilon ) \cdot n [2\ln (N/n)]^{-1} we have E f k ( A Q ) > ϵ f k ( Q ) {\mathcal E} f_k(AQ) > \epsilon f_k(Q) . These asymptotically sharp transitions in the behavior of face numbers as k k varies relative to n log ( N / n ) n \log (N/n) are proven, interpreted, and related to the above-mentioned applications.Keywords
All Related Versions
This publication has 25 references indexed in Scilit:
- For most large underdetermined systems of equations, the minimal 𝓁1‐norm near‐solution approximates the sparsest near‐solutionCommunications on Pure and Applied Mathematics, 2006
- Extensions of compressed sensingSignal Processing, 2006
- High-Dimensional Centrally Symmetric Polytopes with Neighborliness Proportional to DimensionDiscrete & Computational Geometry, 2005
- Decoding by Linear ProgrammingIEEE Transactions on Information Theory, 2005
- Geometric Representation of High Dimension, Low Sample Size DataJournal of the Royal Statistical Society Series B: Statistical Methodology, 2005
- On the inherent intractability of certain coding problems (Corresp.)IEEE Transactions on Information Theory, 1978
- Neighbourliness of centrally symmetric polytopes in high dimensionsMathematika, 1975
- Introduction to Approximation TheoryMathematics of Computation, 1969
- Diagrams for centrally symmetric polytopesMathematika, 1968
- Neighborly and cyclic polytopesProceedings of Symposia in Pure Mathematics, 1963