PhaseCode: Fast and Efficient Compressive Phase Retrieval Based on Sparse-Graph Codes

• Published in: IEEE transactions on information theory Vol. 63; no. 6; pp. 3663 - 3691
• Main Authors: Pedarsani, Ramtin, Dong Yin, Kangwook Lee, Ramchandran, Kannan
• Format: Journal Article
• Language: English
• Published: New York IEEE 01.06.2017; The Institute of Electrical and Electronics Engineers, Inc. (IEEE)
• Subjects: Algorithms; coding theory; Complexity; Complexity theory; compressive sensing; Computation; Computer simulation; Diffraction patterns; Extraterrestrial measurements; Noise measurement; Optical diffraction; Optical variables measurement; Phase measurement; Phase retrieval; Sparse matrices; sparse-graph codes
• ISSN: 0018-9448; 1557-9654
• DOI: 10.1109/TIT.2017.2693287

We consider the problem of recovering a complex signal x ∈ ℂ n from m intensity measurements of the form |a i H x|, 1 ≤ i ≤ m, where a i H is the ith row of measurement matrix A ∈ ℂ m×n . Our main focus is on the case where the measurement vectors are unconstrained, and where x is exactly K-sparse, or the so-called general compressive phase retrieval problem. We introduce PhaseCode, a novel family of fast and efficient algorithms that are based on a sparsegraph coding framework. We show that in the noiseless case, the PhaseCode algorithm can recover an arbitrarily-close-toone fraction of the K nonzero signal components using only slightly more than 4K measurements when the support of the signal is uniformly random, with the order-optimal time and memory complexity of Θ(K). 1 It is known that the fundamental limit for the number of measurements in compressive phase retrieval problem is 4K - o(K) for the more difficult problem of recovering the signal exactly and with no assumptions on its support distribution. This shows that under mild relaxation of the conditions, our algorithm is the first constructive capacity-approaching compressive phase retrieval algorithm: in fact, our algorithm is also order-optimal in complexity and memory. Furthermore, we show that for any signal x, PhaseCode can recover a random (1 - p)-fraction of the nonzero components of x with high probability, where p can be made arbitrarily close to zero, with sample complexity m = c(p)K, where c(p) is a small constant depending on p that can be precisely calculated, with optimal time and memory complexity. As a result, assuming that the nonzero components of x are lower bounded by Θ(1) and upper bounded by Θ(K γ ) for some positive constant γ <; 1, we are able to provide a strong ℓ 1 guarantee for the estimated signal x̂ as follows: ||x̂ - x|| 1 ≤ p||x|| 1 (1+o(1)), where p can be made arbitrarily close to zero. As one instance, the PhaseCode algorithm can provably recover, with high probability, a random 1 - 10 -7 fraction of the significant signal components, using at most m = 14K measurements. Next, motivated by some important practical classes of optical systems, we consider a"Fourier-friendly" constrained measurement setting, and show that its performance matches that of the unconstrained setting, when the signal is sparse in the Fourier domain with uniform support. In the Fourier-friendly setting that we consider, the measurement matrix is constrained to be a cascade of Fourier matrices (corresponding to optical lenses) and diagonal matrices (corresponding to diffraction mask patterns). Finally, we tackle the compressive phase retrieval problem in the presence of noise, where measurements are in the form of y i = |a i H x| 2 + w i , and wi is the additive noise to the w i measurement. We assume that the signal is quantized, and each nonzero component can take L m possible magnitudes and L p possible phases. We consider the regime, where K = βn δ , δ ∈ (0, 1). We use the same architecture of PhaseCode for the noiseless case, and robustify it using two schemes: the almost-linear scheme and the sublinear scheme. We prove that with high probability, the almost-linear scheme recovers x with sample complexity O(K log(n)) and computational complexity Θ(L m L p n log(n)), and the sublinear scheme recovers x with sample complexity Θ(K log 3 (n)) and computational complexity Θ(L m L p K log 3 (n)). Throughout, we provide extensive simulation results that validate the practical power of our proposed algorithms for the sparse unconstrained and Fourier-friendly measurement settings, for noiseless and noisy scenarios.