Abstract:In phase retrieval, we want to recover an unknown signal x∈Cd\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$${{\varvec{x}}}\in {{\mathbb {C}}}^d$$\end{document} from n quadratic measurements of the form yi=|⟨ai,x⟩|2+wi\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$y_i = |\langle {{\varvec{a}}}_i,{{\varvec{x}}}\rangle |^2+w_i$$\end{document}, where ai∈Cd\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$${{\varvec{a}}}_i\in {{\mathbb {C}}}^d$$\end{document} are known sensing vectors and wi\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$w_i$$\end{document} is measurement noise. We ask the following weak recovery question: What is the minimum number of measurements n needed to produce an estimator x^(y)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$${\hat{{{\varvec{x}}}}}({{\varvec{y}}})$$\end{document} that is positively correlated with the signal x\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$${{\varvec{x}}}$$\end{document}? We consider the case of Gaussian vectors ai\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$${{\varvec{a}}}_i$$\end{document}. We prove that—in the high-dimensional limit—a sharp phase transition takes place, and we locate the threshold in the regime of vanishingly small noise. For n≤d-o(d)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$n\le d-o(d)$$\end{document}, no estimator can do significantly better than random and achieve a strictly positive correlation. For n≥d+o(d)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$n\ge d+o(d)$$\end{document}, a simple spectral estimator achieves a positive correlation. Surprisingly, numerical simulations with the same spectral estimator demonstrate promising performance with realistic sensing matrices. Spectral methods are used to initialize non-convex optimization algorithms in phase retrieval, and our approach can boost the performance in this setting as well. Our impossibility result is based on classical information-theoretic arguments. The spectral algorithm computes the leading eigenvector of a weighted empirical covariance matrix. We obtain a sharp characterization of the spectral properties of this random matrix using tools from free probability and generalizing a recent result by Lu and Li. Both the upper bound and lower bound generalize beyond phase retrieval to measurements yi\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$y_i$$\end{document} produced according to a generalized linear model. As a by-product of our analysis, we compare the threshold of the proposed spectral method with that of a message passing algorithm.

SUPER: Sparse Unknown Phases Easily Recovered

Phase Retrieval from Fourier Measurements with Masks

PhaseMax: Stable Guarantees from Noisy Sub-Gaussian Measurements

A Fast and Provable Algorithm for Sparse Phase Retrieval

Phase Retrieval from Quantized Measurements Via Approximate Message Passing

Super Resolution Phase Retrieval for Sparse Signals

Sparse Phase Retrieval with Partial Convolutional Measurements

Subspace Phase Retrieval

Phase Transitions for Greedy Sparse Approximation Algorithms

Sample-Efficient Sparse Phase Retrieval Via Stochastic Alternating Minimization

Fast Compressive Phase Retrieval under Bounded Noise

Block Sparse Recovery With Redundant Measurement Matrices and Its Application in Frequency Agile Radar

Efficient Compressive Phase Retrieval with Constrained Sensing Vectors

Phase Retrieval By Projections

One-Bit Phase Retrieval: Optimal Rates and Efficient Algorithms

Sparse Affine Sampling: Ambiguity-Free and Efficient Sparse Phase Retrieval

Phase Retrieval From the Magnitudes of Affine Linear Measurements

Sparse Support Recovery with Phase-Only Measurements

Phase retrieval in high dimensions: Statistical and computational phase transitions

Uniform Exact Reconstruction of Sparse Signals and Low-rank Matrices from Phase-Only Measurements

Fundamental Limits of Weak Recovery with Applications to Phase Retrieval