FRI Sampling With Arbitrary Kernels

This paper addresses the problem of sampling non-bandlimited signals within the Finite Rate of Innovation (FRI) setting. We had previously shown that, by using sampling kernels whose integer span contains specific exponentials (generalized Strang-Fix conditions), it is possible to devise non-iterative, fast reconstruction algorithms from very low-rate samples. Yet, the accuracy and sensitivity to noise of these algorithms is highly dependent on these exponential reproducing kernels - actually, on the exponentials that they reproduce. Hence, our first contribution here is to provide clear guidelines on how to choose the sampling kernels optimally, in such a way that the reconstruction quality is maximized in the presence of noise. The optimality of these kernels is validated by comparing with Cramér-Rao's lower bounds (CRB). Our second contribution is to relax the exact exponential reproduction requirement. Instead, we demonstrate that arbitrary sampling kernels can reproduce the “best” exponentials within quite a high accuracy in general, and that applying the exact FRI algorithms in this approximate context results in near-optimal reconstruction accuracy for practical noise levels. Essentially, we propose a universal extension of the FRI approach to arbitrary sampling kernels. Numerical results checked against the CRB validate the various contributions of the paper and in particular outline the ability of arbitrary sampling kernels to be used in FRI algorithms.

[1]  Thierry Blu,et al.  Sampling Curves With Finite Rate of Innovation , 2014, IEEE Transactions on Signal Processing.

[2]  Alex ChiChung Kot,et al.  2D Finite Rate of Innovation Reconstruction Method for Step Edge and Polygon Signals in the Presence of Noise , 2012, IEEE Transactions on Signal Processing.

[3]  Vivek K. Goyal,et al.  Estimating Signals With Finite Rate of Innovation From Noisy Samples: A Stochastic Algorithm , 2007, IEEE Transactions on Signal Processing.

[4]  Thierry Blu,et al.  Cardinal exponential splines: part I - theory and filtering algorithms , 2005, IEEE Transactions on Signal Processing.

[5]  Thierry Blu,et al.  Generalized L-spline wavelet bases , 2005, SPIE Optics + Photonics.

[6]  Pedro M. Crespo,et al.  A new stochastic algorithm inspired on genetic algorithms to estimate signals with finite rate of innovation from noisy samples , 2010, Signal Process..

[7]  Michael Unser,et al.  Polynomial spline signal approximations: Filter design and asymptotic equivalence with Shannon's sampling theorem , 1992, IEEE Trans. Inf. Theory.

[8]  James A. Cadzow,et al.  Signal enhancement-a composite property mapping algorithm , 1988, IEEE Trans. Acoust. Speech Signal Process..

[9]  Yonina C. Eldar,et al.  Innovation Rate Sampling of Pulse Streams With Application to Ultrasound Imaging , 2010, IEEE Transactions on Signal Processing.

[10]  Thierry Blu,et al.  MOMS: maximal-order interpolation of minimal support , 2001, IEEE Trans. Image Process..

[11]  Thierry Blu,et al.  ON THE EXPONENTIAL REPRODUCING KERNELS FOR SAMPLING SIGNALS WITH FINITE RATE OF INNOVATION , 2010 .

[12]  Martin Vetterli,et al.  Compressive Sampling of Multiple Sparse Signals Having Common Support Using Finite Rate of Innovation Principles , 2011, IEEE Signal Processing Letters.

[13]  B. Hofmann-Wellenhof,et al.  Introduction to spectral analysis , 1986 .

[14]  Pedro M. Crespo,et al.  Reconstruction of aperiodic FRI signals and estimation of the rate of innovation based on the state space method , 2011, Signal Process..

[15]  Thierry Blu,et al.  Extrapolation and Interpolation) , 2022 .

[16]  M. Vetterli,et al.  Sparse Sampling of Signal Innovations , 2008, IEEE Signal Processing Magazine.

[17]  Yonina C. Eldar,et al.  MMSE whitening and subspace whitening , 2003, IEEE Transactions on Information Theory.

[18]  Bhaskar D. Mo Model Based Processing of Signals: A State Space Approach , 1992 .

[19]  Michael Unser Cardinal exponential splines: part II - think analog, act digital , 2005, IEEE Transactions on Signal Processing.

[20]  Tapan K. Sarkar,et al.  Matrix pencil method for estimating parameters of exponentially damped/undamped sinusoids in noise , 1990, IEEE Trans. Acoust. Speech Signal Process..

[21]  Peyman Milanfar,et al.  A Stable Numerical Method for Inverting Shape from Moments , 1999, SIAM J. Sci. Comput..

[22]  Amos Ron,et al.  Factorization theorems for univariate splines on regular grids , 1990 .

[23]  M. Unser Sampling-50 years after Shannon , 2000, Proceedings of the IEEE.

[24]  Arie Feuer,et al.  On perfect conditioning of Vandermonde matrices on the unit circle , 2007 .

[25]  Thierry Blu,et al.  Sampling Piecewise Sinusoidal Signals With Finite Rate of Innovation Methods , 2010, IEEE Transactions on Signal Processing.

[26]  G. Strang,et al.  A Fourier Analysis of the Finite Element Variational Method , 2011 .

[27]  V. Koivunen,et al.  On the Cramér-Rao bound for the constrained and unconstrained complex parameters , 2008, 2008 5th IEEE Sensor Array and Multichannel Signal Processing Workshop.

[28]  Martin Vetterli,et al.  Exact sampling results for some classes of parametric nonbandlimited 2-D signals , 2004, IEEE Transactions on Signal Processing.

[29]  Martin Vetterli,et al.  Sampling and reconstruction of signals with finite rate of innovation in the presence of noise , 2005, IEEE Transactions on Signal Processing.

[30]  Thierry Blu,et al.  Sampling signals with finite rate of innovation , 2002, IEEE Trans. Signal Process..

[31]  Pier Luigi Dragotti,et al.  Sampling Schemes for Multidimensional Signals With Finite Rate of Innovation , 2007, IEEE Transactions on Signal Processing.