Hadamard Wirtinger Flow for Sparse Phase Retrieval

by   Fan Wu, et al.

We consider the problem of reconstructing an n-dimensional k-sparse signal from a set of magnitude-only measurements. Formulating the problem as an unregularized empirical risk minimization task, we study the sample complexity performance of gradient descent with Hadamard parametrization, which we call Hadamard Wirtinger flow (HWF). Provided knowledge of the signal sparsity k, we prove that a single step of HWF is able to recover the support from O(k(x^*_max)^-2log n) samples, where x^*_max is the largest component of the signal in magnitude. This support recovery procedure can be used to initialize existing reconstruction methods and yields algorithms with total runtime proportional to the cost of reading the data and improved sample complexity, which is linear in k when the signal contains at least one large component. We numerically investigate the performance of HWF at convergence and show that, while not requiring any explicit form of regularization nor knowledge of k, HWF adapts to the signal sparsity and reconstructs sparse signals with fewer measurements than existing state-of-the-art methods.


Sample-Efficient Algorithms for Recovering Structured Signals from Magnitude-Only Measurements

We consider the problem of recovering a signal x^* ∈R^n, from magnitude-...

Sparse Phase Retrieval via Sparse PCA Despite Model Misspecification: A Simplified and Extended Analysis

We consider the problem of high-dimensional misspecified phase retrieval...

Signal Reconstruction from Modulo Observations

We consider the problem of reconstructing a signal from under-determined...

A Continuous-Time Mirror Descent Approach to Sparse Phase Retrieval

We analyze continuous-time mirror descent applied to sparse phase retrie...

Compressive Phase Retrieval: Optimal Sample Complexity with Deep Generative Priors

Advances in compressive sensing provided reconstruction algorithms of sp...

Implicit Regularization for Optimal Sparse Recovery

We investigate implicit regularization schemes for gradient descent meth...

A Universal Sampling Method for Reconstructing Signals with Simple Fourier Transforms

Reconstructing continuous signals from a small number of discrete sample...

Please sign up or login with your details

Forgot password? Click here to reset