Moment inversion problem for piecewise D-finite functions

We consider the problem of exact reconstruction of univariate functions with jump discontinuities at unknown positions from their moments. These functions are assumed to satisfy an a priori unknown linear homogeneous differential equation with polynomial coefficients on each continuity interval. Therefore, they may be specified by a finite amount of information. This reconstruction problem has practical importance in signal processing and other applications.It is somewhat of a 'folklore' that the sequence of the moments of such 'piecewise D-finite' functions satisfies a linear recurrence relation of bounded order and degree. We derive this recurrence relation explicitly. It turns out that the coefficients of the differential operator which annihilates every piece of the function, as well as the locations of the discontinuities, appear in this recurrence in a precisely controlled manner. This leads to the formulation of a generic algorithm for reconstructing a piecewise D-finite function from its moments. We investigate the conditions for solvability of resulting linear systems in the general case, as well as analyse a few particular examples. We provide results of numerical simulations for several types of signals, which test the sensitivity of the proposed algorithm to noise.

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

[2]  George Labahn,et al.  Computation of Numerical Padé-Hermite and Simultaneous Padé Systems I: Near Inversion of Generalized Sylvester Matrices , 1996, SIAM J. Matrix Anal. Appl..

[3]  Knut S. Eckhoff.  Accurate reconstructions of functions of finite regularity from truncated Fourier series expansions , 1995 .

[4]  Zhonggang Zeng.  Algorithm 835: MultRoot---a Matlab package for computing polynomial roots and multiplicities , 2004, TOMS.

[5]  Timothy A. Davis,et al.  MATLAB Primer , 1994 .

[6]  Knut S. Eckhoff.  Accurate and efficient reconstruction of discontinuous functions from truncated series expansions , 1993 .

[7]  Philippe Flajolet,et al.  Analytic Combinatorics , 2009 .

[8]  Tobin A. Driscoll,et al.  A Padé-based algorithm for overcoming the Gibbs phenomenon , 2004, Numerical Algorithms.

[9]  Richard P. Stanley,et al.  Differentiably Finite Power Series , 1980, Eur. J. Comb..

[10]  Riccardo March,et al.  Application of the Padé Method to Solving the Noisy Trigonometric Moment Problem: Some Initial Results , 1998, SIAM J. Appl. Math..

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

[12]  G. W. Stewart,et al.  The behavior of a multiplicity independent root-finding scheme in the presence of error , 1980 .

[13]  D. Zeilberger.  A holonomic systems approach to special functions identities , 1990 .

[14]  N. G. Parke,et al.  Ordinary Differential Equations. , 1958 .

[15]  G. Talenti.  Recovering a function from a finite number of moments , 1987 .

[16]  George Labahn,et al.  Inversion components of block Hankel-like matrices , 1992 .

[17]  P. Milanfar,et al.  Reconstructing planar domains from their moments , 2000 .

[18]  Pier Luigi Dragotti,et al.  Perfect Reconstruction Schemes for Sampling Piecewise Sinusoidal Signals , 2006, 2006 IEEE International Conference on Acoustics Speech and Signal Processing Proceedings.

[19]  Jonathan M. Borwein,et al.  A Proof of a Recurrence for Bessel Moments , 2007, Exp. Math..

[20]  George Kvernadze.  Approximating the jump discontinuities of a function by its Fourier-Jacobi coefficients , 2004, Math. Comput..

[21]  Riccardo March,et al.  Reconstruction of a Piecewise Constant Function from Noisy Fourier Coefficients by Padé Method , 2000, SIAM J. Appl. Math..

[22]  Ana C. Matos,et al.  Reduction of the Gibbs phenomenon for smooth functions with jumps by the ε-algorithm , 2008 .

[23]  O. Ore.  Theory of Non-Commutative Polynomials , 1933 .

[24]  Yosef Yomdin,et al.  Signal Acquisition from Measurements via Non-Linear Models , 2008 .

[25]  Yosef Yomdin,et al.  Linear versus Non-Linear Acquisition of Step-Functions , 2007, ArXiv.

[26]  S. Elaydi.  An introduction to difference equations , 1995 .

[27]  Jonathan M. Borwein,et al.  A Proof of a Recursion for Bessel Moments , 2007, ArXiv.

[28]  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.

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

[30]  R. E. Kalman,et al.  On minimal partial realizations of a linear input/output map , 1971 .

[31]  James F. Geer,et al.  Exponential Approximations Using Fourier Series Partial Sums , 1997 .

[32]  R. Gorenflo,et al.  Moment Theory and Some Inverse Problems in Potential Theory and Heat Conduction , 2002 .

[33]  C. Byrnes,et al.  On the partial stochastic realization problem , 1997, IEEE Trans. Autom. Control..

[34]  Zhonggang Zeng.  Computing multiple roots of inexact polynomials , 2005, Math. Comput..