Abstract:This paper deals with finding an n-dimensional solution x to a system of quadratic equations of the form y i = |〈a i , x〉| 2 for 1 ≤ i ≤ m, which is also known as the generalized phase retrieval problem. For this NP-hard problem, a novel approach is developed for minimizing the amplitude-based leastsquares empirical loss, which starts with a weighted maximal correlation initialization obtainable through a few power or Lanczos iterations, followed by successive refinements based on a sequence of iteratively reweighted gradient iterations. The two stages (initialization and gradient flow) distinguish themselves from prior contributions by the inclusion of a fresh (re)weighting regularization procedure. For certain random measurement models, the novel scheme is shown to be able to recover the true solution x in time proportional to reading the data {(a i ; y i )} 1 ≤i≤m. This holds with high probability and without extra assumption on the signal vector x to be recovered, provided that the number m of equations is some constant c > 0 times the number n of unknowns in the signal vector, namely m > cn. Empirically, the upshots of this contribution are: first, (almost) 100% perfect signal recovery in the high-dimensional (say n ≥ 2000) regime given only an information-theoretic limit number of noiseless equations, namely m = 2n - 1, in the real Gaussian case; and second, (nearly) optimal statistical accuracy in the presence of additive noise of bounded support. Finally, substantial numerical tests using both synthetic data and real images corroborate markedly improved recovery performance and computational efficiency of the novel scheme relative to the state-of-the-art approaches.

Phase retrieval from incomplete data via weighted nuclear norm minimization

PhaseMax: Stable Guarantees from Noisy Sub-Gaussian Measurements

Single-frame Noisy Interferogram Phase Retrieval Using an End-to-end Deep Learning Network with Physical Information Constraints

Phase Retrieval from Quantized Measurements Via Approximate Message Passing

Geophysical Data Sparse Reconstruction Based on L0-Norm Minimization

Phase Retrieval Via Nonlocal Complex-Domain Sparsity.

Phase Retrieval: A Data-Driven Wavelet Frame Based Approach

Undersampled Phase Retrieval via Majorization-Minimization

Phase Retrieval Via Reweighted Amplitude Flow

Low Complexity Regularized Phase Retrieval

Phase Retrieval from Multiple FRFT Measurements Based on Nonconvex Low-Rank Minimization

Seismic Data Reconstruction Via Weighted Nuclear-Norm Minimization

Robust Compressive Phase Retrieval via L1 Minimization With Application to Image Reconstruction

Weighted Phase Retrieval of Fourier Measurement With Outliers: Measurement Structure and Reconstruction Algorithm

Undersampled Phase Retrieval with Outliers

Robust Image Retrieval From Phaseless Data

Complex Noise-Based Phase Retrieval Using Total Variation and Wavelet Transform Regularization

Phase retrieval in high dimensions: Statistical and computational phase transitions

Weighted Null Vector Initialization and Its Application to Phase Retrieval

Image Recovery With Data Missing In The Presence Of Salt-And-Pepper Noise

Recovery Conditions in Weighted Sparse Phase Retrieval Via Weighted $$\ell _q\, (0<q\le 1)$$ Minimization