Abstract
This paper develops a novel algorithm, termed SPARse Truncated Amplitude flow (SPARTA), to reconstruct a sparse signal from a small number of magnitude-only measurements. It deals with what is also known as sparse phase retrieval (PR), which is NP-hard in general and emerges in many science and engineering applications. Upon formulating sparse PR as an amplitude-based nonconvex optimization task, SPARTA works iteratively in two stages: In stage one, the support of the underlying sparse signal is recovered using an analytically well-justified rule, and subsequently a sparse orthogonality-promoting initialization is obtained via power iterations restricted on the support; and in the second stage, the initialization is successively refined by means of hard thresholding based gradient-type iterations. SPARTA is a simple yet effective, scalable, and fast sparse PR solver. On the theoretical side, for any n-dimensional k-sparse (k n) signal x with minimum (in modulus) nonzero entries on the order of (1/k)x2 , SPARTA recovers the signal exactly (up to a global unimodular constant) from about k2 log n random Gaussian measurements with high probability. Furthermore, SPARTA incurs computational complexity on the order of k2 n log n with total runtime proportional to the time required to read the data, which improves upon the state of the art by at least a factor of k. Finally, SPARTA is robust against additive noise of bounded support. Extensive numerical tests corroborate markedly improved recovery performance and speedups of SPARTA relative to existing alternatives.
Original language | English (US) |
---|---|
Pages (from-to) | 479-491 |
Number of pages | 13 |
Journal | IEEE Transactions on Signal Processing |
Volume | 66 |
Issue number | 2 |
DOIs | |
State | Published - Jan 15 2018 |
Bibliographical note
Funding Information:Manuscript received December 11, 2016; revised May 25, 2017 and August 25, 2017; accepted October 29, 2017. Date of publication November 8, 2017; date of current version December 22, 2017. The associate editor coordinating the review of this manuscript and approving it for publication was Prof. Tsung-Hui Chang. The work of G. Wang, L. Zhang, and G. B. Giannakis was supported by NSF Grants 1500713 and 1514056. The work of M. Akc¸akaya was supported by NSF Grant 1651825. The work of J. Chen was partially supported by the National Natural Science Foundation of China Grants U1509215, 61621063, and the Program for Changjiang Scholars and Innovative Research Team in University (IRT1208). This paper was presented in part at the 42nd IEEE International Conference on Acoustics, Speech, and Signal Processing, New Orleans, LA, USA, March 5–9, 2017. (Corresponding author: Georgios B. Giannakis.) G. Wang is with the Digital Technology Center and the Department of Electrical and Computer Engineering, University of Minnesota, Minneapolis, MN 55455 USA, and also with the State Key Laboratory of Intelligent Control and Decision of Complex Systems, Beijing Institute of Technology, Beijing 100081, China (e-mail: gangwang@umn.edu).
Funding Information:
The work of G. Wang, L. Zhang, and G. B. Giannakis was supported by NSF Grants 1500713 and 1514056. The work of M. Ak?akaya was supported by NSF Grant 1651825. The work of J. Chen was partially supported by the National Natural Science Foundation of China Grants U1509215, 61621063, and the Program for Changjiang Scholars and Innovative Research Team in University (IRT1208). This paper was presented in part at the 42nd IEEE International Conference on Acoustics, Speech, and Signal Processing, New Orleans, LA, USA, March 5?9, 2017. The authors would like to thank the anonymous reviewers for their thorough review and all constructive comments and suggestions, which helped to improve the quality of the manuscript. The authors also thank Prof. Xiaodong Li for sharing the codes of the thresholded Wirtinger flow algorithm.
Publisher Copyright:
© 2017 IEEE.
Keywords
- Compressive sampling
- Iterative hard thresholding
- Linear convergence to the global optimum.
- Nonconvex optimization
- Support recovery