The problem of recovering a signal from a quadratic system with full-rank matrices frequently arises in applications such as unassigned distance geometry and sub-wavelength imaging. With i.i.d. standard Gaussian matrices , this paper addresses the high-dimensional case where by incorporating prior knowledge of . First, we consider a -sparse and introduce the thresholded Wirtinger flow (TWF) algorithm that does not require the sparsity level . TWF comprises two steps: the spectral initialization that identifies a point sufficiently close to (up to a sign flip) when , and the thresholded gradient descent (with a good initialization) that produces a sequence linearly converging to with measurements. Second, we explore the generative prior, assuming that lies in the range of an -Lipschitz continuous generative model with -dimensional inputs in an -ball of radius . We develop the projected gradient descent (PGD) algorithm that also comprises two steps: the projected power method that provides an initial vector with -error given measurements, and the projected gradient descent that refines the -error to at a geometric rate when . Experimental results corroborate our theoretical findings and show that: (i) our approach for the sparse case notably outperforms the existing provable algorithm sparse power factorization; (ii) leveraging the generative prior allows for precise image recovery in the MNIST dataset from a small number of quadratic measurements.
View on arXiv