Abstract
We present a novel method to iteratively calculate discrete Fourier transforms for discrete time signals with sample time intervals that may be widely nonuniform. The proposed recursive Fourier transform (RFT) does not require interpolation of the samples to uniform time intervals, and each iterative transform update of N frequencies has computational order N. Because of the inherent non-uniformity in the time between successive heart beats, an application particularly well suited for this transform is power spectral density (PSD) estimation for heart rate variability. We compare RFT based spectrum estimation with Lomb–Scargle Transform (LST) based estimation. PSD estimation based on the LST also does not require uniform time samples, but the LST has a computational order greater than Nlog(N). We conducted an assessment study involving the analysis of quasi-stationary signals with various levels of randomly missing heart beats. Our results indicate that the RFT leads to comparable estimation performance to the LST with significantly less computational overhead and complexity for applications requiring iterative spectrum estimations.
Similar content being viewed by others
References
Alexander S (1986) Fast adaptive filters: a geometrical approach 3(4):18–28
Cadzow J (1990) Signal processing via least squares error modeling 7(4):12–31. doi:10.1109/53.62941
Clifford G, Tarassenko L (2005) Quantifying errors in spectral estimates of hrv due to beat replacement and resampling. Biomed Eng IEEE Trans 52(4):630–638
Cooke WH, Salinas J, Convertino VA, Ludwig DA, Hinds D, Duke JH, Moore FA, Holcomb JB (2006) Heart rate variability and its association with mortality in prehospital trauma patients. J Trauma 60:363–370
Ellenby MS, McNames J, Lai S, McDonald BA, Krieger D, Sclabassi RJ, Goldstein B (2001) Uncoupling and recoupling of autonomic regulation of the heart beat in pediatric septic shock. Shock 16(4):274–377
Glentis GO, Berberidis K, Theodoridis S (1999) Efficient least squares adaptive algorithms for fir transversal filtering 16(4):13–41. doi:10.1109/79.774932
Hilton MF, Bates RA, Godfrey KR, Chappell MJ, Cayton RM, Gil E, Mantaras C, Aiolfi S, Cerutti S (1999) Evaluation of frequency and time–frequency spectral analysis of heart rate variability as a diagnostic marker of the sleep apnoea syndrome. Med Biol Eng Comput 37(1):760–769
HRV (1996) Task Force of the European Society of Cardiology and the North American Society of Pacing & Electrophysiology: Heart rate variability standards of measurement, physiological interpretation, and clinical use. Eur Heart J 17:354–381
Laguna P, Moody G, Mark R (1998) Power spectral density of unevenly sampled data by least-square analysis: performance and application to heart rate signals. IEEE Trans Biomed Eng 45(6):698–715
Lomb N (1976) Least-squares frequency analysis of unequally spaced data. Astrophys Space Sci 39:447–462
Mendez MO, Bianchi AM, Montano N, Patruno V, Gil E, Mantaras C, Aiolfi S, Cerutti S (2008) On arousal from sleep: time–frequency analysis. Med Biol Eng Comput 46(4):341–351
Plemmons R (1993) Fft-based rls in signal processing. In: Proc. IEEE International Conference on Acoustics, Speech, and Signal Processing ICASSP-93, vol 3, pp 571–574. doi:10.1109/ICASSP.1993.31956210.1109/ICASSP.1993.319562
Press WH, Teukolsky SA, Vetterling WT, Flannery BP (2002) Numerical Recipes in C++: the art of scientific computing. Cambridge University Press, Cambridge
Rajendra Acharya U, Paul Joseph K, Kannathal N, Lim CM, Suri JS (2006) Heart rate variability: a review. Med Biol Eng Comput 44(12):1031–1051. http://dx.doi.org/10.1007/s11517-006-0119-0
Sayed AH, Kailath T (1994) A state-space approach to adaptive RLS filtering. Signal Process Mag IEEE 11:18–60, 3 July. doi:10.1109/79.295229. ISSN: 1053-5888
Scargle JD (1982) Studies in astronomical time series analysis ii, statistical aspects of spectral analysis of unevenly spaced data. Astrophy J 263:835–853
Schwab K, Eiselt M, Putsche P, Helbig M, Witte H (2006) Time-variant parametric estimation of transient quadratic phase couplings between heart rate components in healthy neonates. Med Biol Eng Comput 44(12):1077–1083
Strang G (1986) Introduction to applied mathematics. Wellesley-Cambridge Press
Thong T (2006) Real-time evaluation of spectral heart rate variability. In: Proceedings of the IEEE EMBS annual international conference (EMBS’06), New York City, USA
Thong T, Raitt MH (2005) Ventricular tachyarrhythmia prediction. In: Proceedings of the 2005 IEEE EMBS 27th annual conference
Thong T, McNames J, Aboy M (2004) Lomb–Welch periodogram for non-uniform sampling. In: Proceedings of the 26th annual international conference of the IEEE EMBS, pp 271–274
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Holland, A., Aboy, M. A novel recursive Fourier transform for nonuniform sampled signals: application to heart rate variability spectrum estimation. Med Biol Eng Comput 47, 697–707 (2009). https://doi.org/10.1007/s11517-009-0461-0
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11517-009-0461-0