Nested Sampling and its Applications in Stable Compressive Covariance Estimation and Phase Retrieval with Near-Minimal Measurements

Title of thesis: Nested Sampling and its Applications in Stable Compressive Covariance Estimation and Phase Retrieval with Near-Minimal Measurements Heng Qiao, Master of Science, 2016 Thesis directed by: Professor Piya Pal Department of Electrical and Computer Engineering Compressed covariance sensing using quadratic samplers is gaining increasing interest in recent literature. Covariance matrix often plays the role of a sufficient statistic in many signal and information processing tasks. However, owing to the large dimension of the data, it may become necessary to obtain a compressed sketch of the high dimensional covariance matrix to reduce the associated storage and communication costs. Nested sampling has been proposed in the past as an efficient sub-Nyquist sampling strategy that enables perfect reconstruction of the autocorrelation sequence of Wide-Sense Stationary (WSS) signals, as though it was sampled at the Nyquist rate. The key idea behind nested sampling is to exploit properties of the difference set that naturally arises in quadratic measurement model associated with covariance compression. In this thesis, we will focus on developing novel versions of nested sampling for low rank Toeplitz covariance estimation, and phase retrieval, where the latter problem finds many applications in high resolution optical imaging, X-ray crystallography and molecular imaging. The problem of low rank compressive Toeplitz covariance estimation is first shown to be fundamentally related to that of line spectrum recovery. In absence if noise, this connection can be exploited to develop a particular kind of sampler called the Generalized Nested Sampler (GNS), that can achieve optimal compression rates. In presence of bounded noise, we develop a regularization-free algorithm that provably leads to stable recovery of the high dimensional Toeplitz matrix from its order-wise minimal sketch acquired using a GNS. Contrary to existing TV-norm and nuclear norm based reconstruction algorithms, our technique does not use any tuning parameters, which can be of great practical value. The idea of nested sampling idea also finds a surprising use in the problem of phase retrieval, which has been of great interest in recent times for its convex formulation via PhaseLift, By using another modified version of nested sampling, namely the Partial Nested Fourier Sampler (PNFS), we show that with probability one, it is possible to achieve a certain conjectured lower bound on the necessary measurement size. Moreover, for sparse data, an l1 minimization based algorithm is proposed that can lead to stable phase retrieval using order-wise minimal number of measurements. Nested Sampling and its Applications in Stable Compressive Covariance Estimation and Phase Retrieval with Near-Minimal Measurements by Heng Qiao Thesis submitted to the Faculty of the Graduate School of the University of Maryland, College Park in partial fulfillment of the requirements for the degree of Master of Science 2016 Advisory Committee: Professor Piya Pal, Chair/Advisor Professor Rama Chellappa Professor Min Wu

[1]  D. Donoho 50 Years of Data Science , 2017 .

[2]  Emmanuel J. Candès,et al.  Super-Resolution of Positive Sources: The Discrete Setup , 2015, SIAM J. Imaging Sci..

[3]  Piya Pal,et al.  Generalized nested sampling for compression and exact recovery of symmetric Toeplitz matrices , 2014, 2014 IEEE Global Conference on Signal and Information Processing (GlobalSIP).

[4]  Richard G Baraniuk,et al.  More Is Less: Signal Processing and the Data Deluge , 2011, Science.

[5]  Jeffrey A. Fessler,et al.  Undersampled Phase Retrieval With Outliers , 2014, IEEE Transactions on Computational Imaging.

[6]  Hamid Krim,et al.  Two Decades of Array Signal Processing , 1997 .

[7]  Jian Li,et al.  Computationally efficient maximum-likelihood estimation of structured covariance matrices , 1998, Proceedings of the 1998 IEEE International Conference on Acoustics, Speech and Signal Processing, ICASSP '98 (Cat. No.98CH36181).

[8]  U. Grenander,et al.  Toeplitz Forms And Their Applications , 1958 .

[9]  Prateek Jain,et al.  Phase Retrieval Using Alternating Minimization , 2013, IEEE Transactions on Signal Processing.

[10]  Yonina C. Eldar,et al.  Sub-Nyquist Sampling for Power Spectrum Sensing in Cognitive Radios: A Unified Approach , 2013, IEEE Transactions on Signal Processing.

[11]  Yuejie Chi,et al.  Guaranteed Blind Sparse Spikes Deconvolution via Lifting and Convex Optimization , 2015, IEEE Journal of Selected Topics in Signal Processing.

[12]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

[13]  Babak Hassibi,et al.  Recovery of sparse 1-D signals from the magnitudes of their Fourier transform , 2012, 2012 IEEE International Symposium on Information Theory Proceedings.

[14]  Benjamin Recht,et al.  Superresolution without separation , 2015, 2015 IEEE 6th International Workshop on Computational Advances in Multi-Sensor Adaptive Processing (CAMSAP).

[15]  Dustin G. Mixon,et al.  Saving phase: Injectivity and stability for phase retrieval , 2013, 1302.4618.

[16]  Babak Hassibi,et al.  On robust phase retrieval for sparse signals , 2012, 2012 50th Annual Allerton Conference on Communication, Control, and Computing (Allerton).

[17]  Emmanuel J. Candès,et al.  PhaseLift: Exact and Stable Signal Recovery from Magnitude Measurements via Convex Programming , 2011, ArXiv.

[18]  Chandra Sekhar Seelamantula,et al.  Fienup Algorithm With Sparsity Constraints: Application to Frequency-Domain Optical-Coherence Tomography , 2014, IEEE Transactions on Signal Processing.

[19]  Yonina C. Eldar,et al.  On conditions for uniqueness in sparse phase retrieval , 2013, 2014 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[20]  J. Makhoul,et al.  Linear prediction: A tutorial review , 1975, Proceedings of the IEEE.

[21]  Piya Pal,et al.  Sparse phase retrieval using partial nested fourier samplers , 2015, 2015 IEEE Global Conference on Signal and Information Processing (GlobalSIP).

[22]  Holger Boche,et al.  Fast compressive phase retrieval from Fourier measurements , 2015, 2015 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[23]  Brian M. Sadler,et al.  Cyclic Feature Detection With Sub-Nyquist Sampling for Wideband Spectrum Sensing , 2012, IEEE Journal of Selected Topics in Signal Processing.

[24]  Xiaodong Li,et al.  Solving Quadratic Equations via PhaseLift When There Are About as Many Equations as Unknowns , 2012, Found. Comput. Math..

[25]  P. Bickel,et al.  Regularized estimation of large covariance matrices , 2008, 0803.1909.

[26]  Emmanuel J. Candès,et al.  Exact Matrix Completion via Convex Optimization , 2008, Found. Comput. Math..

[27]  Petre Stoica,et al.  Maximum likelihood methods for direction-of-arrival estimation , 1990, IEEE Trans. Acoust. Speech Signal Process..

[28]  Yaniv Plan,et al.  Robust 1-bit Compressed Sensing and Sparse Logistic Regression: A Convex Programming Approach , 2012, IEEE Transactions on Information Theory.

[29]  Yonina C. Eldar,et al.  Phase Retrieval via Matrix Completion , 2011, SIAM Rev..

[30]  Piya Pal,et al.  Gridless Line Spectrum Estimation and Low-Rank Toeplitz Matrix Compression Using Structured Samplers: A Regularization-Free Approach , 2017, IEEE Transactions on Signal Processing.

[31]  R. Balan,et al.  On signal reconstruction without phase , 2006 .

[32]  Andrea J. Goldsmith,et al.  Exact and Stable Covariance Estimation From Quadratic Sampling via Convex Programming , 2013, IEEE Transactions on Information Theory.

[33]  Gongguo Tang,et al.  Near minimax line spectral estimation , 2013, 2013 47th Annual Conference on Information Sciences and Systems (CISS).

[34]  Gongguo Tang,et al.  Atomic Norm Denoising With Applications to Line Spectral Estimation , 2012, IEEE Transactions on Signal Processing.

[35]  Thomas Strohmer,et al.  Self-calibration and biconvex compressive sensing , 2015, ArXiv.

[36]  David Gross,et al.  Recovering Low-Rank Matrices From Few Coefficients in Any Basis , 2009, IEEE Transactions on Information Theory.

[37]  Emmanuel J. Candès,et al.  Near-Optimal Signal Recovery From Random Projections: Universal Encoding Strategies? , 2004, IEEE Transactions on Information Theory.

[38]  Piya Pal,et al.  Generalized Nested Sampling for Compressing Low Rank Toeplitz Matrices , 2015, IEEE Signal Processing Letters.

[39]  Zhi Ding,et al.  Blind MIMO channel identification from second order statistics using rank deficient channel convolution matrix , 2003, IEEE Trans. Signal Process..

[40]  B. Moore Principal component analysis in linear systems: Controllability, observability, and model reduction , 1981 .

[41]  Sundeep Rangan,et al.  Compressive Phase Retrieval via Generalized Approximate Message Passing , 2014, IEEE Transactions on Signal Processing.

[42]  Nanda Kambhatla,et al.  Dimension Reduction by Local Principal Component Analysis , 1997, Neural Computation.

[43]  Martin J. Wainwright,et al.  A unified framework for high-dimensional analysis of $M$-estimators with decomposable regularizers , 2009, NIPS.

[44]  Yonina C. Eldar,et al.  GESPAR: Efficient Phase Retrieval of Sparse Signals , 2013, IEEE Transactions on Signal Processing.

[45]  Piya Pal,et al.  Sparse phase retrieval with near minimal measurements: A structured sampling based approach , 2016, 2016 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[46]  Tomas Bäckström,et al.  Vandermonde Factorization of Toeplitz Matrices and Applications in Filtering and Warping , 2013, IEEE Transactions on Signal Processing.

[47]  Geert Leus,et al.  Compressive Wideband Power Spectrum Estimation , 2012, IEEE Transactions on Signal Processing.

[48]  Piya Pal,et al.  Stable compressive low rank Toeplitz covariance estimation without regularization , 2016, 2016 50th Asilomar Conference on Signals, Systems and Computers.

[49]  Yonina C. Eldar,et al.  Phase Retrieval with Application to Optical Imaging: A contemporary overview , 2015, IEEE Signal Processing Magazine.

[50]  Yonina C. Eldar,et al.  Phase Retrieval from 1D Fourier Measurements: Convexity, Uniqueness, and Algorithms , 2016, IEEE Transactions on Signal Processing.

[51]  Justin K. Romberg,et al.  Efficient Compressive Phase Retrieval with Constrained Sensing Vectors , 2015, NIPS.

[52]  Xiaodong Li,et al.  Phase Retrieval via Wirtinger Flow: Theory and Algorithms , 2014, IEEE Transactions on Information Theory.

[53]  Emmanuel J. Candès,et al.  The Power of Convex Relaxation: Near-Optimal Matrix Completion , 2009, IEEE Transactions on Information Theory.

[54]  Babak Hassibi,et al.  Sparse Phase Retrieval: Uniqueness Guarantees and Recovery Algorithms , 2013, IEEE Transactions on Signal Processing.

[55]  Phillipp Meister,et al.  Statistical Signal Processing Detection Estimation And Time Series Analysis , 2016 .

[56]  Petre Stoica,et al.  Spectral Analysis of Signals , 2009 .

[57]  Emmanuel J. Candès,et al.  Decoding by linear programming , 2005, IEEE Transactions on Information Theory.

[58]  Xiaodong Li,et al.  Phase Retrieval from Coded Diffraction Patterns , 2013, 1310.3240.

[59]  Geert Leus,et al.  Recovering second-order statistics from compressive measurements , 2011, 2011 4th IEEE International Workshop on Computational Advances in Multi-Sensor Adaptive Processing (CAMSAP).

[60]  Parikshit Shah,et al.  Compressed Sensing Off the Grid , 2012, IEEE Transactions on Information Theory.

[61]  Peter Jung,et al.  Stable recovery from the magnitude of symmetrized fourier measurements , 2014, 2014 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[62]  James R Fienup,et al.  Phase retrieval algorithms: a personal tour [Invited]. , 2013, Applied optics.

[63]  Björn E. Ottersten,et al.  Structured covariance matrix estimation: a parametric approach , 2000, 2000 IEEE International Conference on Acoustics, Speech, and Signal Processing. Proceedings (Cat. No.00CH37100).

[64]  Emmanuel J. Candès,et al.  Super-Resolution from Noisy Data , 2012, Journal of Fourier Analysis and Applications.

[65]  Geert Leus,et al.  Compressive covariance sampling , 2013, 2013 Information Theory and Applications Workshop (ITA).

[66]  Zhiqiang Xu,et al.  Phase Retrieval for Sparse Signals , 2013, ArXiv.

[67]  Richard G. Baraniuk,et al.  Kronecker Compressive Sensing , 2012, IEEE Transactions on Image Processing.

[68]  Huaiyu Zhu On Information and Sufficiency , 1997 .

[69]  Stephen P. Boyd,et al.  A rank minimization heuristic with application to minimum order system approximation , 2001, Proceedings of the 2001 American Control Conference. (Cat. No.01CH37148).

[70]  Emmanuel J. Candès,et al.  Towards a Mathematical Theory of Super‐resolution , 2012, ArXiv.

[71]  Robert Bregovic,et al.  Multirate Systems and Filter Banks , 2002 .

[72]  Allen Y. Yang,et al.  CPRL -- An Extension of Compressive Sensing to the Phase Retrieval Problem , 2012, NIPS.

[73]  Yonina C. Eldar,et al.  Phase Retrieval: Stability and Recovery Guarantees , 2012, ArXiv.

[74]  Geert Leus,et al.  Compression Limits for Random Vectors with Linearly Parameterized Second-Order Statistics , 2013, IEEE Transactions on Information Theory.

[75]  R. Gerchberg A practical algorithm for the determination of phase from image and diffraction plane pictures , 1972 .

[76]  Yonina C. Eldar,et al.  Sparse Phase Retrieval from Short-Time Fourier Measurements , 2014, IEEE Signal Processing Letters.

[77]  Yonina C. Eldar,et al.  Simultaneously Structured Models With Application to Sparse and Low-Rank Matrices , 2012, IEEE Transactions on Information Theory.

[78]  P. P. Vaidyanathan,et al.  Nested Arrays: A Novel Approach to Array Processing With Enhanced Degrees of Freedom , 2010, IEEE Transactions on Signal Processing.

[79]  S. Frick,et al.  Compressed Sensing , 2014, Computer Vision, A Reference Guide.

[80]  Yuejie Chi,et al.  Off-the-Grid Line Spectrum Denoising and Estimation With Multiple Measurement Vectors , 2014, IEEE Transactions on Signal Processing.

[81]  Parikshit Shah,et al.  Sketching Sparse Matrices, Covariances, and Graphs via Tensor Products , 2015, IEEE Transactions on Information Theory.

[82]  F. Gamboa,et al.  Spike detection from inaccurate samplings , 2013, 1301.5873.

[83]  Holger Rauhut,et al.  Low rank matrix recovery from rank one measurements , 2014, ArXiv.

[84]  Xiaodong Li,et al.  Sparse Signal Recovery from Quadratic Measurements via Convex Programming , 2012, SIAM J. Math. Anal..

[85]  Holger Rauhut,et al.  A Mathematical Introduction to Compressive Sensing , 2013, Applied and Numerical Harmonic Analysis.

[86]  Martin Vetterli,et al.  Phase Retrieval for Sparse Signals: Uniqueness Conditions , 2013, ArXiv.

[87]  Yohann de Castro,et al.  Exact Reconstruction using Beurling Minimal Extrapolation , 2011, 1103.4951.

[88]  Piya Pal,et al.  Finite sample analysis of covariance compression using structured samplers , 2016, 2016 IEEE Sensor Array and Multichannel Signal Processing Workshop (SAM).

[89]  Holger Rauhut,et al.  Stable low-rank matrix recovery via null space properties , 2015, ArXiv.

[90]  I. S. Gál On the Representation of 1, 2, . . . , N by Differences , 2004 .