On compressed blind de-convolution of filtered sparse processes

Suppose the signal x ∈ 葷n is realized by driving a k-sparse signal z ∈ 葷n through an arbitrary unknown stable discrete-linear time invariant system H, namely, x(t) = (h * z)(t), where h(·) is the impulse response of the operator H. Is x(·) compressible in the conventional sense of compressed sensing? Namely, can x(t) be reconstructed from small set of measurements obtained through suitable random projections? For the case when the unknown system H is auto-regressive (i.e. all pole) of a known order it turns out that x can indeed be reconstructed from O(k log(n)) measurements. We develop a novel LP optimization algorithm and show that both the unknown filter H and the sparse input z can be reliably estimated.

[1]  I. Rubin,et al.  Random point processes , 1977, Proceedings of the IEEE.

[2]  Holger Rauhut,et al.  Circulant and Toeplitz matrices in compressed sensing , 2009, ArXiv.

[3]  Enders A. Robinson,et al.  Seismic time‐invariant convolutional model , 1985 .

[4]  Erick Baziw,et al.  Principle phase decomposition: a new concept in blind seismic deconvolution , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[5]  Jerry D. Gibson,et al.  Handbook of Image and Video Processing , 2000 .

[6]  Joel A. Tropp,et al.  Signal Recovery From Random Measurements Via Orthogonal Matching Pursuit , 2007, IEEE Transactions on Information Theory.

[7]  D. Donoho For most large underdetermined systems of linear equations the minimal 𝓁1‐norm solution is also the sparsest solution , 2006 .

[8]  Venkatesh Saligrama,et al.  Deterministic Designs with Deterministic Guarantees: Toeplitz Compressed Sensing Matrices, Sequence Designs and System Identification , 2008, ArXiv.

[9]  Justin K. Romberg,et al.  Compressive Sensing by Random Convolution , 2009, SIAM J. Imaging Sci..

[10]  Martin Vetterli,et al.  Annihilating filter-based decoding in the compressed sensing framework , 2007, SPIE Optical Engineering + Applications.

[11]  R. DeVore,et al.  A Simple Proof of the Restricted Isometry Property for Random Matrices , 2008 .

[12]  Kjetil F. Kaaresen,et al.  Deconvolution of sparse spike trains by iterated window maximization , 1997, IEEE Trans. Signal Process..

[13]  E. Candès,et al.  Near-ideal model selection by ℓ1 minimization , 2008, 0801.0345.

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

[15]  Christophe Andrieu,et al.  Bayesian deconvolution of noisy filtered point processes , 2001, IEEE Trans. Signal Process..

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

[17]  J WainwrightMartin Sharp thresholds for high-dimensional and noisy sparsity recovery using l1-constrained quadratic programming (Lasso) , 2009 .

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

[19]  E. Candès,et al.  Stable signal recovery from incomplete and inaccurate measurements , 2005, math/0503066.

[20]  Wulfram Gerstner,et al.  SPIKING NEURON MODELS Single Neurons , Populations , Plasticity , 2002 .