Power spectral density of unevenly sampled data by least-square analysis: performance and application to heart rate signals
- 1 June 1998
- journal article
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Transactions on Biomedical Engineering
- Vol. 45 (6) , 698-715
- https://doi.org/10.1109/10.678605
Abstract
This work studies the frequency behavior of a least-square method to estimate the power spectral density of unevenly sampled signals. When the uneven sampling can be modeled as uniform sampling plus a stationary random deviation, this spectrum results in a periodic repetition of the original continuous time spectrum at the mean Nyquist frequency, with a low-pass effect affecting upper frequency bands that depends on the sampling dispersion. If the dispersion is small compared with the mean sampling period, the estimation at the base band is unbiased with practically no dispersion. When uneven sampling is modeled by a deterministic sinusoidal variation respect to the uniform sampling the obtained results are in agreement with those obtained for small random deviation. This approximation is usually well satisfied in signals like heart rate (HR) series. The theoretically predicted performance has been tested and corroborated with simulated and real HR signals. The Lomb method has been compared with the classical power spectral density (PSD) estimators that include resampling to get uniform sampling. The authors have found that the Lomb method avoids the major problem of classical methods: the low-pass effect of the resampling. Also only frequencies up to the mean Nyquist frequency should be considered (lower than 0.5 Hz if the HR is lower than 60 bpm). It is concluded that for PSD estimation of unevenly sampled signals the Lomb method is more suitable than fast Fourier transform or autoregressive estimate with linear or cubic interpolation. In extreme situations (low-HR or high-frequency components) the Lomb estimate still introduces high-frequency contamination that suggest further studies of superior performance interpolators. In the case of HR signals the authors have also marked the convenience of selecting a stationary heart rate period to carry out a heart rate variability analysis.Keywords
This publication has 28 references indexed in Scilit:
- ECG-based indices of physical activityPublished by Institute of Electrical and Electronics Engineers (IEEE) ,2003
- The European ST-T database: development, distribution and usePublished by Institute of Electrical and Electronics Engineers (IEEE) ,2002
- Pole-tracking algorithms for the extraction of time-variant heart rate variability spectral parametersIEEE Transactions on Biomedical Engineering, 1995
- An Efficient Algorithm for Spectral Analysis of Heart Rate VariabilityIEEE Transactions on Biomedical Engineering, 1986
- Coherent averaging technique: A tutorial review Part 1: Noise reduction and the equivalent filterJournal of Biomedical Engineering, 1986
- Spectrum of a series of point events, generated by the integral pulse frequency modulation modelMedical & Biological Engineering & Computing, 1985
- The Measurement of Heart Rate Variability Spectra with the Help of a Personal ComputerIEEE Transactions on Biomedical Engineering, 1982
- The Shannon sampling theorem—Its various extensions and applications: A tutorial reviewProceedings of the IEEE, 1977
- Spectral Analysis of Pulse Frequency Modulation in the Nervous SystemsIEEE Transactions on Biomedical Engineering, 1968
- On Nonuniform Sampling of Bandwidth-Limited SignalsIRE Transactions on Circuit Theory, 1956