Adaptive Filter-bank Approach to Restoration and Spectral Analysis of Gapped Data*

The main topic of this paper is the nonparametric estimation of complex (both amplitude and phase) spectra from gapped data, as well as the restoration of such data. The focus is on the extension of the APES (amplitude and phase estimation) approach to data sequences with gaps. APES, which is one of the most successful existing nonparametric approaches to the spectral analysis of full data sequences, uses a bank of narrowband adaptive (both frequency and data dependent) filters to estimate the spectrum. A recent interpretation of this approach showed that the filterbank used by APES and the resulting spectrum minimize a least-squares (LS) fitting criterion between the filtered sequence and its spectral decomposition. The extended approach, which is called GAPES for somewhat obvious reasons, capitalizes on the aforementioned interpretation: it minimizes the APES-LS fitting criterion with respect to the missing data as well. This should be a sensible thing to do whenever the full data sequence is stationary, and hence the missing data have the same spectral content as the available data. We use both simulated and real data examples to show that GAPES estimated spectra and interpolated data sequences have excellent accuracy. We also show the performance gain achieved by GAPES over two of the most commonly used approaches for gapped-data spectral analysis, viz., the periodogram and the parametric CLEAN method.

[1]  Mauricio D. Sacchi,et al.  Interpolation and extrapolation using a high-resolution discrete Fourier transform , 1998, IEEE Trans. Signal Process..

[2]  T. J. Ulrych,et al.  A new method for estimating the power spectrum of gapped data , 1982 .

[3]  Petre Stoica,et al.  Introduction to spectral analysis , 1997 .

[4]  P. R. Swan Discrete Fourier transforms of nonuniformly spaced data , 1982 .

[5]  Jeff R. Kuhn Recovering spectral information from unevenly sampled data - Two machine-efficient solutions , 1982 .

[6]  N. Lomb Least-squares frequency analysis of unequally spaced data , 1976 .

[7]  B. Anderson,et al.  Greatest common divisor via generalized Sylvester and Bezout matrices , 1978 .

[8]  D. H. Roberts,et al.  Time Series Analysis with Clean - Part One - Derivation of a Spectrum , 1987 .

[9]  Jian Li,et al.  A new derivation of the APES filter , 1999, IEEE Signal Processing Letters.

[10]  W. Press,et al.  Interpolation, realization, and reconstruction of noisy, irregularly sampled data , 1992 .

[11]  J. Scargle Studies in astronomical time series analysis. II - Statistical aspects of spectral analysis of unevenly spaced data , 1982 .

[12]  J. Scargle Studies in astronomical time series analysis. III - Fourier transforms, autocorrelation functions, and cross-correlation functions of unevenly spaced data , 1989 .

[13]  Jian Li,et al.  An adaptive filtering approach to spectral estimation and SAR imaging , 1996, IEEE Trans. Signal Process..

[14]  W. Press,et al.  Fast algorithm for spectral analysis of unevenly sampled data , 1989 .

[15]  Thomas P. Bronez,et al.  Spectral estimation of irregularly sampled multidimensional processes by generalized prolate spheroidal sequences , 1988, IEEE Trans. Acoust. Speech Signal Process..

[16]  Andreas Jakobsson,et al.  Matched-filter bank interpretation of some spectral estimators , 1998, Signal Process..

[17]  Fionn Murtagh,et al.  In Astronomical Data Analysis Software and Systems IV , 1995 .

[18]  Timothy M. Brown,et al.  A technique for estimating complicated power spectra from time series with gaps , 1990 .

[19]  Steven Kay,et al.  Modern Spectral Estimation: Theory and Application , 1988 .

[20]  Alf Isaksson,et al.  Identification of ARX-models subject to missing data , 1993, IEEE Trans. Autom. Control..

[21]  D. D. Meisel Fourier transforms of data sampled in unequally spaced segments. , 1979 .

[22]  Yingbo Hua,et al.  A subspace method for the computation of the GCD of polynomials , 1997, Autom..