Automatic Decomposition of Wigner Distribution and its Application to Heart Rate Variability

Summary

Objective: We introduce an algorithm for the automatic decomposition of Wigner Distribution (WD) and we applied it for the quantitative extraction of Heart Rate Variability (HRV) spectral parameters during non-stationary events. Early response to tilt was investigated.

Methods: Quantitative analysis of multi-components non-stationary signals is obtained through an automatic decomposition of WD based on least square (LS) fitting of the instantaneous autocorrelation function (ACF). Through this approach the different signal and interference terms which contributes to the ACF may be separated and their parameters (instantaneous frequency and amplitude) quantified. A beat-to-beat monitoring of HRV spectral components is obtained.

Results: Analysis of simulated signals demonstrated the capability of the proposed approach to track and separate the signal components. Analysis of HRV data evidenced different dynamics in the early Autonomic Nervous System (ANS) response to tilt.

Conclusions: The novel approach to the quantification of the beat-to-beat HRV spectral parameters obtained from decomposition of Wigner distribution was demonstrated to be effective in the analysis of HRV data. Relevant physiological information about the dynamics of the early sympathetic response to tilt were obtained. The method is a general approach which may be employed for a quantitative time-frequency analysis of non-stationary biological signals.

Keywords

Time-frequency distribution - autonomic nervous system - least square fitting

• References

• 1 Task Force of ESC/NASPE. Heart rate variability: standards of measurement, physiological interpretation and clinical use. Circulation 1996; 93 (05) 1043-65.
• 2 Novak P, Novak V. Time/frequency mapping of the heart rate, blood pressure and respiratory signals. Med Biol Eng Comput 1993; 31 (02) 103-10.
• 3 Bianchi AM, Mainardi L, Petrucci E, Signorini MG, Mainardi M, Cerutti S. Time-variant power spectrum analysis for the detection of transient episodes in HRV signal. IEEE Trans Biomed Eng 1993; 40 (02) 136-44.
• 4 Kumaresan R, Tufts DW. “Estimating the parameters of exponentially dumped sinusoids and pole zero modeling in noise”. IEEE Trans ASSP 1982; 30 (06) 833-40.
• 5 Hlawatsch F, Boudreaux-Bartels GF. Linear and Quadratic Time-frequency Signal Representations. IEEE SP Magazine (April) 1992; 9 (02) 21-67.
• 6 Montano, et al. Power spectrum analysis of heart rate variability to assess the changes in sympathovagal balance during graded orthostatic tilt. Circ 1994; 90 (04) 1826-31.