Spectral Method for Phase Retrieval: an Expectation Propagation Perspective

Abstract

Phase retrieval refers to the problem of recovering a signal x∈Cn from its phaseless measurements yi=|aiHx|, where \ai\i=1m are the measurement vectors. Many popular phase retrieval algorithms are based on the following two-step procedure: (i) initialize the algorithm based on a spectral method, (ii) refine the initial estimate by a local search algorithm (e.g., gradient descent). The quality of the spectral initialization step can have a major impact on the performance of the overall algorithm. In this paper, we focus on the model where the measurement matrix A=[a1,…,am]H has orthonormal columns, and study the spectral initialization under the asymptotic setting m,n∞ with m/nδ∈(1,∞). We use the expectation propagation framework to characterize the performance of spectral initialization for Haar distributed matrices. Our numerical results confirm that the predictions of the EP method are accurate for not-only Haar distributed matrices, but also for realistic Fourier based models (e.g. the coded diffraction model). The main findings of this paper are the following: (1) There exists a threshold on δ (denoted as δweak) below which the spectral method cannot produce a meaningful estimate. We show that δweak=2 for the column-orthonormal model. In contrast, previous results by Mondelli and Montanari show that δweak=1 for the i.i.d. Gaussian model. (2) The optimal design for the spectral method coincides with that for the i.i.d. Gaussian model, where the latter was recently introduced by Luo, Alghamdi and Lu.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…