The problem of reconstructing a sparse signal vector from magnitude-only measurements (a.k.a. compressive phase retrieval) emerges in diverse applications, but it is NP-hard in general. Building on recent advances in nonconvex optimization, this paper puts forth a new algorithm that is termed compressive reweighted amplitude flow (CRAF) and abbreviated as CRAF, for compressive phase retrieval. Specifically, CRAF operates in two stages. The first stage seeks a sparse initial guess via a new spectral procedure. In the second stage, CRAF implements a few hard thresholding based iterations using reweighted gradients. When there are sufficient measurements, CRAF provably recovers the underlying signal vector exactly with high probability under suitable conditions. Moreover, its sample complexity coincides with that of state-of-the-art procedures. Finally, numerical tests showcase improved performance of the new spectral initialization, as well as improved exact recovery relative to competing alternatives. Empirically, the number of measurements required for exact recovery by CRAF is about 40% less than those of competing alternatives in our experiments using real images.
- Nonconvex optimization
- iteratively reweighting
- linear convergence to the global optimum
- model-based hard thresholding