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.

Sparse Signal Recovery from Phaseless Measurements Via Hard Thresholding Pursuit

A Fast and Provable Algorithm for Sparse Phase Retrieval

Phase Retrieval from Fourier Measurements with Masks

Phase Retrieval with PhaseLift Algorithm

PhaseMax: Stable Guarantees from Noisy Sub-Gaussian Measurements

Sparse Signals Recovery from Noisy Measurements by Orthogonal Matching Pursuit

Phase Retrieval from Quantized Measurements Via Approximate Message Passing

GraHTP: A Provable Newton-like Algorithm for Sparse Phase Retrieval

Sample-Efficient Sparse Phase Retrieval Via Stochastic Alternating Minimization

Sparse Phase Retrieval with Partial Convolutional Measurements

Sparse Signal Phase Retrieval for Phaseless Short-Time Fourier Transform Measurement Based on Local Search

A Stochastic Alternating Minimizing Method for Sparse Phase Retrieval

Projective Iterative Hard Thresholding Algorithm for Sparse Signal Recovery

One-Bit Phase Retrieval: Optimal Rates and Efficient Algorithms

A Pseudo-Inverse-Based Hard Thresholding Algorithm for Sparse Signal Recovery

An Alternating Optimization Approach For Phase Retrieval

Extended OMP Algorithm for Sparse Phase Retrieval.

Fast Compressive Phase Retrieval under Bounded Noise

Sparse representation utilizing tight frame for phase retrieval

A ReLU-based Hard-Thresholding Algorithm for Non-Negative Sparse Signal Recovery

Fundamental Limits of Weak Recovery with Applications to Phase Retrieval