## Abstract

We discuss the relationship between the RR interval spectral analysis and the spectral analysis of the corresponding ECG signal from which the RR intervals were evaluated. The ECG signal spectrum is bounded below the frequency fB by using an electronic filter and sampled at rate larger than 2fB, thus excluding aliasing from spectral analysis. A similar procedure cannot be applied to the RR interval spectral analysis, and in this case aliasing is possible. One of our main efforts in this chapter is devoted to the problem of how to detect aliasing in the heart rate spectral analysis. In order to get an insight, we performed an experiment with an adult man, in which the ECG signal was detected in a case where the breathing rate was larger than half the heart rate. A constant breathing rate for time intervals exceeding 5 minutes was monitored with good accuracy using a special breathing procedure. The results show distinctively a very sharp peak in the spectral analysis of the ECG signal, and corresponding (diffused) aliasing peaks in the RR interval spectral analysis. A new method of dealing with unevenly sampled data was developed, which has interesting anti-aliasing properties. There are indications that the VLF peaks of the RR spectrum are originated by aliasing. Some of the LF peaks may have the same property. The chapter is fully based on the preprint arXiv:physics/9911017, submitted on 11 Nov 1999, by authors A. Gersten, O. Gersten, A. Ronen, and Y. Cassuto.

### Keywords

- heart rate
- ECG signal
- spectral analysis
- aliasing
- HRV

## 1. Introduction

The RR interval spectral analysis is usually based on heart rate data collected in two ways. In one method, the data are collected by analog to digital conversion of the ECG signal and computer evaluation of the RR intervals from the ECG signal. In the second method, devices are used whose output is the RR interval alone. The advantage of the first method is the control of accuracy and flexibility of the evaluations. The second method has the advantage of storing smaller amount of data, and it can be easily used online.

In the first method, usually the number of collected data (sampled ECG signal) is of two to three orders of magnitude larger than the RR interval data. Thus if only RR interval is analyzed, a large amount of data is unused. In this paper we are trying to take advantage of the ECG sampled signal and to derive new information in addition to the conventional RR interval analysis [1, 2, 3, 4, 5].

The ECG signal spectrum is bounded below the frequency f by using an electronic filter and sampled at rate larger than 2f

_{B}

*, thus excluding aliasing from spectral analysis [6]. A similar procedure cannot be applied to the RR interval spectral analysis, and in this case an aliasing is possible. One of our main efforts in this paper is devoted to the problem of how to detect aliasing in the RR interval spectral analysis.*

_{B}

In order to get an insight, we performed an experiment, in which the ECG signal of one of the authors (AG) was detected, while the breathing rate was larger than half the heart rate. A constant breathing rate for a time exceeding 5 minutes was monitored with good accuracy using a special breathing procedure with a metronome. The results show distinctively a very sharp peak in the spectral analysis of the ECG signal and corresponding (diffused) aliasing peaks in the RR interval spectral analysis.

The spectral analysis of the ECG signal was performed with the standard FFT procedures. The spectral analysis of the RR intervals was performed with several techniques in order to take into consideration that the data were unevenly sampled. This is presented in Section 2. In Section 3, we discuss the possibility of aliasing in the spectral analysis of the RR intervals. In Section 4, we compare power estimations of ECG’s and RR intervals of three experiments. In Section 5, we analyze the results. In Section 6, summary and conclusions are presented.

## 2. Spectral analysis of unevenly sampled data

The methods of spectral analysis are well developed for evenly sampled data [6, 7]. The RR interval data are unevenly sampled in time. In most cases an analysis is performed with respect to beat numbers which are evenly spaced. We will below justify this method using least square principles. But as was recently indicated by Laguna et al. [8], the resampling of data is causing the appearance of additional harmonics. They recommend to use a method developed by Lomb [9]. The errors of resampling the beats can, to large extent, be overcome by using a cubic spline interpolation. In this work we are suggesting a new method of treating unevenly sampled data, which, unexpectedly, gave good results beyond the Nyquist frequency.

### 2.1 Analysis according to beat numbers

Let us assume that the RR intervals are given at unevenly sampled times

and let us generate in the interval

We will use the discrete Fourier transform (DFT) for a basis formed from the evenly sampled times

The coefficient

with the result

Eqs. (3 and 5) can be handled easily with standard FFT programs. This is the usual procedure which is adopted in most of the papers dealing with RR interval analysis [4, 5].

### 2.2 Other methods

FFT can be applied more efficiently if the unevenly sampled data are interpolated at evenly spaced intervals of Eq. (2). The cubic spline interpolation is one of the good ways to do it.

The Lomb method [9] was extensively analyzed in Ref. [8]. We give here only the formulae in the form of the Lomb normalized periodogram:

where

### 2.3 Nonuniform discrete Fourier transform (NUDFT)

We present here a new method of treating unevenly spaced events which we call the “nonuniform discrete Fourier transform” (NUDFT).

Let us assume that

Our aim is to find a good approximation to this expression in terms of the unevenly sampled signal

We start with the Euler summation formula:

and make the following decomposition of the integral on the right hand side of Eq. (9)

and approximate each of the integrals on the right hand side with the trapezoidal rule

From Eqs. (9) and (11), we obtain:

When

and the final result, the approximation to Eq. (8), after rearranging the terms, becomes

where

with the inverse formula

which is an interpolation formula for

## 3. Aliasing

Aliasing is a result of undersampling and is a well-known phenomenon. In Ref. [10], aliasing was looked upon from the point of view of symmetry. It is an example of wrong symmetry and as such should be given more attention. It is the outcome of an incomplete basis. It was found in Ref. [10] that for evenly sampled data with a sampling rate

where

In order to avoid the aliasing symmetry of Eq. (17), the frequencies should be bounded by the Nyquist frequency (denoted here by

The ECG signal was sampled with sampling rate 250 Hz, and an electronic filter was applied, which have eliminated practically all frequencies above 32 Hz, thus aliasing cannot occur at frequencies below 125 Hz or even below 32 Hz. The RR intervals were calculated directly from the ECG signal. The sampling rate for RR intervals can be defined only for evenly sampled data and for the methods that interpolate the unevenly sampled data, or one can consider the average sampling rate from Eq. (1) in both cases:

where

Another possibility of detecting aliasing is by comparing the heart rate spectrum with the ECG signal spectrum. Marked differences below the Nyquist frequency for the power distribution of the RR intervals compared to the ECG signal power distribution in the same range may indicate aliasing. But we do not have yet a sound basis to treat this problem.

We have devised an experiment which definitely demonstrates the aliasing in the RR interval spectrum. To the best of our knowledge, this is the first experiment in which one can exactly know the correct frequency above the Nyquist frequency and can follow the development of the aliasing, which appears to be diffused to great extent because the symmetry of Eq. (17) is represented not by one sampling rate but by a distribution of sampling rates, as the RR interval is unevenly sampled.

Below we describe three experiments. One of them was devised to demonstrate aliasing and the other two for learning about the relations between the RR interval spectrum and the spectrum of the ECG signal.

## 4. Three experiments

We present below results of three experiments. In the first experiment, the ECG signal was collected in a normal resting state. The aim of this experiment was to compare the ECG spectrum with the RR interval spectrum. In the second experiment, very slow breathing was monitored at a rate of 0.04 Hz. Again the ECG and RR interval spectra were compared. In the third experiment, very fast breathing was accurately monitored at the rate of 74/min and 84/min. These respiratory rates were above half of the heart rates, thus allowing to observe in detail the development of aliasing.

### 4.1 The first experiment

In this experiment (linked with the names of Zahi and Ori, where the second is one of the authors: O.G), which was done in normal, resting conditions, we compare the power estimation of the RR interval and the ECG signal, from which the RR interval was obtained. The ECG signal was sampled at a rate of 250 Hz. Stable intervals of 7-minute duration were chosen for analysis.

In Figure 1a the power distribution of the ECG signal of Zahi is depicted. The attenuation of the power with increasing frequency above 12 Hz is due to the action of an electronic filter. Above 32 Hz the contribution is practically zero. The average heart rate was 0.97 Hz. The above results were zoomed to the interval [0–12] Hz in Figure 1b. One can see distinctively the peak around the average heart rate and the higher harmonics of this peak. The second harmonic is missing, but the third, fourth, fifth, and sixth are distinctively visible; higher harmonics became more and more smeared and indistinguishable above the sixth harmonic. One should also note the large difference in power in the heart rate range, below the Nyquist frequency of 0.49 Hz, which is much smaller than the peak around the average heart rate 0.97 Hz.

The power distribution of the RR intervals in the range {0–0.5} Hz was computed according to the methods discussed in Section 2 and is presented in Figure 2a (DFT, beat number analysis), Figure 2b (Spline interpolation), and Figure 2c (NUDFT). For comparison also the power distribution of the ECG signal in the above range is presented in Figure 2d.

The results of Figure 2a–c are quite similar, but the spline interpolation (Figure 2b) and the NUDFT (Figure 2c) are practically identical. The three graphs show the structure commonly found in the power estimation analysis of RR intervals, namely, the existence of the “high-frequency” (HF), “low-frequency” (LF), and the “very low-frequency” (VLF) peaks. The ECG spectrum shows qualitatively the same structure (but not a quantitative agreement), except that the ECG spectrum is highly suppressed below 0.04 Hz, in the VLF region, indicating a possibility of aliasing in this region in the RR analysis.

In Figures 3 and 4a–d, the results of Ori are presented. The conclusions are similar to those of Zahi, except that in the ECG spectrum, both VLF and LF peaks are missing, indicating the possibility of aliasing in these regions for the RR analysis. Also in the ECG spectrum of Ofek, VLF and LF, present in Figure 5a, are missing. VLF is missing in J.C.’s ECG spectrum (see Figure 6a–6b).

### 4.2 The second experiment

In this experiment (linked again with the name Ori), we have checked the ECG spectrum near the VLF region, as the VLF was absent in the ECG spectrum for the resting state in the first experiment. The question was whether such a result persists in all ECG spectra. Therefore we have probed the VLF region by monitoring very prolonged breathing with a rate of 0.04 Hz. For the spectrum of RR intervals, we found that the DFT, spline interpolation, and NUDFT give similar results, and again NUDFT was practically identical to the spline interpolation. Therefore we present only the results of NUDFT, which are presented in Figure 7a. For comparison the spectrum of the ECG signal is given in Figure 7b. In Figure 7a one can see a very clean pattern of a peak at 0.04 Hz and its higher harmonics. In Figure 7b one can see a similar but somewhat diffused pattern. Thus this experiment indicates that similar respiratory patterns exist in both the RR and in the ECG signals.

### 4.3 The third experiment

In this experiment (linked to the name Alex, who is one of the authors: AG), very fast breathing was accurately monitored at the rate of 74/min and 84/min, respectively. These rates were well above half of the average heart rate, thus allowing to observe in detail the development of aliasing. In Figure 8 the ECG spectrum is dominated by the very high and narrow peak at the frequency

In principle the NUDFT and the Lomb methods should not be used above the Nyquist frequency. Surprisingly enough we have found that both methods have a sharp peak at

Similar results for the breathing frequency 84/min are presented in Figures 11–12.

## 5. Further developments

Since our experiment, which demonstrated how aliasing is developing in human beings, nobody had performed experiments on human beings. The reason is that till now, nobody dared (except one of us, AG) to do extremely fast breathing of 74 breaths/min and 84/min, for more than 5 minutes. After reading our preprint, Campbell [17] and his colleagues found an aliasing in fish [17]. Other researchers were more concerned with preventing aliasing, observing the phenomenon in speeded heart rate, and in constructing aliasing filters [18, 19, 20, 21].

## 6. Summary and conclusions

The ECG signal spectrum is bounded below the Nyquist frequency f by using an electronic filter and sampled at rate larger than 2f

_{B}

*, thus excluding aliasing from spectral analysis. A similar procedure cannot be applied to the RR interval spectral analysis, and in this case an aliasing is possible. One of our main efforts in this paper was devoted to the problem of how to detect aliasing in the RR interval spectral analysis.*

_{B}

In order to get insight into this problem, three experiments have been analyzed. In the first experiment, the ECG signal was collected in a normal resting state. The aim of this experiment was to compare the ECG spectrum with the RR interval spectrum. In the second experiment very slow breathing was monitored at a rate of 0.04 Hz. Again the ECG and RR interval spectra were compared. In the third experiment, very fast breathing was accurately monitored at the rate of 74/min and 84/min, respectively. These respiratory rates were above half of the heart rates, thus allowing to observe in detail the development of aliasing.

The experiments which were described above led us to the following conclusions:

The spectral analysis of the ECG signal is more sensitive and accurate than the RR interval spectral analysis and is free from aliasing. Still in the present stage, it contains too much information to be of practical use. Efforts should be made to understand what will be the best way to extract information (not related to the heart condition alone as in the standard analysis of ECG) about the external influences on the heart signal.

We have conducted an experiment which gave a clear insight about the mechanism of aliasing in the RR interval spectrum. The very sharp peak in the spectrum of the ECG signal, which came as the result of enforced quick breathing, reappeared as a diffused signal in the RR spectrum. The extension of the diffuseness agrees with the extension of the sampling rates of unevenly sampled data.

The VLF peak observed in the RR interval spectrum is usually missing in the ECG spectrum. This leads us to suspect that the VLF observed in the RR spectrum has its origin in aliasing.

In some cases the LF peak does not show up in the ECG spectrum. This led us to suspect that part of the LF peak is of aliasing origin.

Unlike in electronic devices, it is very difficult to devise procedures to detect aliasing in humans. In electronic devices aliasing can be easily detected by changing the sampling rate. In humans the fluctuations of the heart rate are of the same order as the required changes in the sampling rates. It will be an important task to develop a proper procedure for detecting aliasing in humans.

We have developed a new technique for spectral analysis for unevenly sampled data called nonuniform discrete Fourier transform (NUDFT). When employed to the RR data, below the Nyquist frequency, it gave similar results as those obtained by interpolating the data with a cubic spline. Above the Nyquist frequency, the correct peak in the spectrum was detected with great accuracy. A similar result was obtained with the recently rediscovered Lomb method. We interpret this unexpected result by a partial destruction of aliasing symmetry in both methods. More efforts should be made in order to understand the anti-aliasing properties of the above methods.

We consider aliasing to be a wrong symmetry, resulting from the use of an incomplete basis, which has intrinsic symmetries inconsistent with the properties of the signal. Aliasing can be partially removed by reducing the symmetry of the basis.

## References

- 1.
Sayers BMA. Analysis of heart rate variability. Ergonomics . 1973;16 :85-97 - 2.
Akselrod SD et al. Power spectrum analysis of heart rate fluctuations: A quantitative probe of beat-to-beat cardiovascular control. Science . 1981;213 :213-220 - 3.
Kamath MV, Fallen EL. Power spectral analysis of HRV: A noninvasive signature of cardiac autonomic functions. Critical Reviews in Biomedical Engineering. 1993; 21 :245-311 - 4.
Malik M et al. Heart rate variability, standards of measurements, physiological interpretation, and clinical use. Circulation . 1998;93 :1043-1065 - 5.
Malik M, Camm AJ, editors. Heart Rate Variability . Armonk NY: Futura; 1995 - 6.
Oppenheim AV, Schafer RW. Discrete-Time Signal Processing . Englewood Cliffs NJ: Prentice-Hall; 1989 - 7.
Cohen L. Time-Frequency Analysis . Englewood Cliffs NJ: Prentice-Hall; 1995 - 8.
Laguna P, Moody GB, Mark RG. Power spectral density of unevenly sampled data by least-square analysis: Performance and application to heart rate signals. IEEE Transactions on Biomedical Engineering . 1998;45 :698-715 - 9.
Lomb NR. Least-squares frequency analysis of unequally spaced data. Astrophysical and Space Science . 1976;39 :447-462 - 10.
Gersten A. Dirac’s representation theory as a framework for signal theory. I. Discrete finite signals. Annals of Physics (N.Y.) . 1998;262 :47-72 - 11.
Witte H et al. Evidence of a previously undescribed form of respiratory sinus arrhythmia (RSA)–the physiological manifestation of cardiac aliasing. Pflügers Archiv . 1988;412 :442-444 - 12.
Rother M et al. Cardiac aliasing–a possible cause for the misinterpretation of cardiorespirographic data in neonates. Early Human Development . 1989;20 :1-12 - 13.
Nilsson J, Panizza M, Hallett M. Principles of digital sampling of a physiologic signal. Electroencephalography and Clinical Neurophysiology . 1993;89 :349-358 - 14.
Zwiener U et al. Heart rate fluctuations in rabbits during different behavioral states. Biomedica Biochimica Acta . 1990;49 :59-68 - 15.
Zwiener U et al. Forms of physiological aliasing within the heart rate fluctuations by higher frequent respiratory movements. Journal of Physiology and Pharmacology . 1994;45 :563-572 - 16.
Zwiener U et al. Heart rate fluctuations of lower frequencies than the respiratory rhythm but caused by it. Pflügers Archiv . 1995;429 :455-461 - 17.
Capbell HA, Kledpacki JZ, Eggington S. A new method in applying power spectral statistics to examine cardio-respiratory interactions in fish. Journal of Theoretical Biology . 2006;241 :410-419 - 18.
Morren G. Advanced signal processing applied to in-vivo spectroscopy and heart rate variability [thesis]. Katholieke Universiteit Leuven; 2004 - 19.
Keresnyei R. Selecting the optimal anti-aliasing filter for multichannel biosignal acquisition intended for inter-signal phase shift analysis. Physiological Measurement. 2015; 36 :N23-N24 - 20.
J. Pučik, O. Ondraček Heart rate variability spectrum: Physiologic aliasing and nonstationarity considerations. In: Trends in Biomedical Engineering; Bratislava. 16–18 Sept. 2009 - 21.
Toledo E, Pinhas I, Aravot D, Akselrod S. Very high frequency oscillations in the heart rate and blood pressure of heart transplant patients. Medical & Biological Engineering & Computing. 2003; 41 :432-438