Application of the Padé Method to Solving the Noisy Trigonometric Moment Problem: Some Initial Results

The noisy trigonometric moment problem for a finite linear combination of box functions is considered, and a research program, possibly leading to a superresolving method, is outlined and some initial steps are performed. The method is based on the remark that the poles of the Pade approximant to the Z-transform of the noiseless moments show, asymptotically, a regular pattern in the complex plane. The pattern can be described by a set of arcs, connecting points on the unit circle, and a pole density function defined on the arcs. When a moderate noise affects the moments, more arcs are needed to describe the pole pattern, but the noiseless pattern is slightly deformed, still allowing its identification. When this identification is possible, a very effective noise filter and moment extrapolator should be easily constructed. In this paper only some preliminary steps of the above research program are performed. Specifically, the case of one box function is considered. A method for computing the pole patterns,...

[1]  Fabrice Gamboa,et al.  The maximum entropy method on the mean: Applications to linear programming and superresolution , 1994, Math. Program..

[2]  Enrico Massaro,et al.  The segmented Prony method for the analysis of non-stationary time series. , 1989 .

[3]  J. Nuttall,et al.  Asymptotic Form of Hermite-Padé Polynomials , 1981 .

[4]  P. Barone,et al.  Hormonally induced modulation in the phosphate metabolites of breast cancer: Analysis of in vivo 31P MRS signals with a modified prony method , 1997, Magnetic resonance in medicine.

[5]  P. Barone,et al.  Modified Prony Method to Resolve and Quantify in Vivo31P NMR Spectra of Tumors , 1994 .

[6]  Piero Barone,et al.  A new method of magnetic resonance image reconstruction with short acquisition time and truncation artifact reduction , 1992, IEEE Trans. Medical Imaging.

[7]  J Nuttall,et al.  Asymptotics of diagonal Hermite-Padé polynomials , 1984 .

[8]  B. Logan,et al.  Signal recovery and the large sieve , 1992 .

[9]  D. Donoho Superresolution via sparsity constraints , 1992 .

[10]  C. Brezinski Padé-type approximation and general orthogonal polynomials , 1980 .

[11]  J. Nuttall,et al.  Convergence of Padé approximants to quasianalytic functions beyond natural boundaries , 1973 .

[12]  Gene H. Golub,et al.  Matrix computations , 1983 .

[13]  P. Jansson Deconvolution : with applications in spectroscopy , 1984 .

[14]  J. Skilling,et al.  Maximum entropy signal processing in practical NMR spectroscopy , 1984, Nature.

[15]  J. Nuttall,et al.  THE CONVERGENCE OF PADÉ APPROXIMANTS TO FUNCTIONS WITH BRANCH POINTS , 1977 .

[16]  D. Slepian,et al.  Prolate spheroidal wave functions, fourier analysis and uncertainty — II , 1961 .

[17]  J Nuttall,et al.  Orthogonal polynomials and Padé approximants associated with a system of arcs , 1977 .

[18]  W. Hackbusch Singular Integral Equations , 1995 .

[19]  I. Johnstone,et al.  Maximum Entropy and the Nearly Black Object , 1992 .

[20]  G. A. Baker Essentials of Padé approximants , 1975 .

[21]  E. Kosarev Shannon's superresolution limit for signal recovery , 1990 .

[22]  Michel Loève,et al.  Probability Theory I , 1977 .

[23]  Edward J. Wegman,et al.  Statistical Signal Processing , 1985 .