The RR Interval Spectrum, the ECG Signal, and Aliasing

Alexander Gersten, Ori Gersten, Adi Ronen and Yair Cassuto

## 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 f<sup>B</sup> 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–5].

The ECG signal spectrum is bounded below the frequency f<sup>B</sup> by using an electronic filter and sampled at rate larger than 2fB, 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.

σ ¼ ∑ N n¼1

analysis [4, 5].

PXð Þ¼ ω<sup>k</sup>

defined as

∑ N n¼1

Eq. (9)

41

fð Þ¼ τ<sup>n</sup>

2.2 Other methods

the good ways to do it.

1 2σ<sup>2</sup> 8 ><

>:

∑<sup>N</sup>

Eq. (2). The corresponding DFT is

unevenly sampled signal s tð Þ<sup>n</sup> .

1 Δτ ð<sup>τ</sup><sup>N</sup> τ1

∑<sup>N</sup>

with the result

s tð Þ� <sup>n</sup>

1 <sup>N</sup> <sup>∑</sup> N k¼1

DOI: http://dx.doi.org/10.5772/intechopen.85327

The RR Interval Spectrum, the ECG Signal, and Aliasing

Sk exp ð Þ iωkτ<sup>n</sup> � � s tð Þ� <sup>n</sup>

> Sk ¼ ∑ N n¼1

the formulae in the form of the Lomb normalized periodogram:

<sup>n</sup>¼<sup>1</sup>½ � s tð Þ� <sup>n</sup> <sup>s</sup> cosð Þ <sup>ω</sup>kð Þ tn � <sup>τ</sup> h i<sup>2</sup>

<sup>n</sup>¼<sup>1</sup> cos <sup>2</sup>ð Þ <sup>ω</sup>kð Þ tn � <sup>τ</sup> <sup>þ</sup>

tan 2ð Þ¼ <sup>ω</sup>k<sup>τ</sup> <sup>∑</sup><sup>N</sup>

2.3 Nonuniform discrete Fourier transform (NUDFT)

the "nonuniform discrete Fourier transform" (NUDFT).

We start with the Euler summation formula:

1 2

fð Þτ dτ þ

Sk ¼ ∑ N n¼1

where s and σ<sup>2</sup> are the mean and variance of the data and the value of τ is

∑<sup>N</sup>

We present here a new method of treating unevenly spaced events which we call

Let us assume that sð Þ τ<sup>n</sup> are the exact values of the signal at the points given by

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

<sup>½</sup> <sup>f</sup>ð Þþ <sup>τ</sup><sup>1</sup> <sup>f</sup>ð Þ <sup>τ</sup><sup>N</sup> � þ <sup>Δ</sup><sup>τ</sup>

and make the following decomposition of the integral on the right hand side of

1 <sup>N</sup> <sup>∑</sup> N k¼1 S ∗

� � � � (4)

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

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 Lomb method [9] was extensively analyzed in Ref. [8]. We give here only

∑<sup>N</sup>

<sup>n</sup>¼<sup>1</sup> sin 2ð Þ <sup>ω</sup>ktn

<sup>n</sup>¼<sup>1</sup> cos 2ð Þ <sup>ω</sup>ktn

<sup>12</sup> <sup>f</sup> 0 ð Þ� τ<sup>N</sup> f 0 ð Þ τ<sup>1</sup> � � <sup>þ</sup> <sup>O</sup> <sup>Δ</sup>τ<sup>2</sup> � � (9)

∑<sup>N</sup>

sð Þ τ<sup>n</sup> exp ð Þ �iωkτ<sup>n</sup> : (8)

<sup>n</sup>¼<sup>1</sup>½ � s tð Þ� <sup>n</sup> <sup>s</sup> sin ð Þ <sup>ω</sup>kð Þ tn � <sup>τ</sup> h i<sup>2</sup>

9 >=

>;

(6)

<sup>n</sup>¼<sup>1</sup> sin <sup>2</sup>ð Þ <sup>ω</sup>kð Þ tn � <sup>τ</sup>

: (7)

<sup>k</sup> exp ð Þ �iωkτ<sup>n</sup>

s tð Þ<sup>n</sup> exp ð Þ �iωkτ<sup>n</sup> : (5)

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 tn, with the values s tð Þ<sup>n</sup> , where n is the beat number, n ¼ 1⋯N. Let us divide the interval ½ � t1; tN into equal subintervals:

$$
\Delta \tau = \frac{t\_N - t\_1}{N - 1},
\tag{1}
$$

and let us generate in the interval ½ � t1; tN evenly sampled times:

$$
\pi\_n = (n-1)\Delta\tau + t\_1. \tag{2}
$$

We will use the discrete Fourier transform (DFT) for a basis formed from the evenly sampled times τn. We will assume that

$$s(t\_n) = \frac{1}{N} \sum\_{k=1}^{N} \mathcal{S}\_k \exp\left(i\alpha\_k \tau\_n\right), \qquad \alpha\_k = 2\pi(k-1)/(N\Delta\tau). \tag{3}$$

The coefficient Sk will be determined by minimizing the expression

The RR Interval Spectrum, the ECG Signal, and Aliasing DOI: http://dx.doi.org/10.5772/intechopen.85327

$$\sigma = \sum\_{n=1}^{N} \left\{ \left[ s(t\_n) - \frac{1}{N} \sum\_{k=1}^{N} \mathbb{S}\_k \exp\left( i\alpha y\_k \tau\_n \right) \right] \left[ s(t\_n) - \frac{1}{N} \sum\_{k=1}^{N} \mathbb{S}\_k^\* \exp\left( -i\alpha y\_k \tau\_n \right) \right] \right\} \tag{4}$$

with the result

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

Fourier Transforms - Century of Digitalization and Increasing Expectations

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

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

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

Let us assume that the RR intervals are given at unevenly sampled times tn, with the values s tð Þ<sup>n</sup> , where n is the beat number, n ¼ 1⋯N. Let us divide the interval

<sup>N</sup> � <sup>1</sup> , (1)

τ<sup>n</sup> ¼ ð Þ n � 1 Δτ þ t1: (2)

Sk exp ð Þ iωkτ<sup>n</sup> , ω<sup>k</sup> ¼ 2πð Þ k � 1 =ð Þ NΔτ : (3)

<sup>Δ</sup><sup>τ</sup> <sup>¼</sup> tN � <sup>t</sup><sup>1</sup>

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

and let us generate in the interval ½ � t1; tN evenly sampled times:

The coefficient Sk will be determined by minimizing the expression

results. In Section 6, summary and conclusions are presented.

2. Spectral analysis of unevenly sampled data

2.1 Analysis according to beat numbers

evenly sampled times τn. We will assume that

1 <sup>N</sup> <sup>∑</sup> N k¼1

s tð Þ¼ <sup>n</sup>

½ � t1; tN into equal subintervals:

spectral analysis.

analysis.

frequency.

40

$$S\_k = \sum\_{n=1}^{N} s(t\_n) \exp\left(-i\alpha y\_k \tau\_n\right). \tag{5}$$

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:

$$P\_X(\alpha\_k) = \frac{1}{2\sigma^2} \left\{ \frac{\left[\sum\_{n=1}^N [s(t\_n) - \overline{s}] \cos\left(\alpha\_k(t\_n - \tau)\right)\right]^2}{\sum\_{n=1}^N \cos^2(\alpha\_k(t\_n - \tau))} + \frac{\left[\sum\_{n=1}^N [s(t\_n) - \overline{s}] \sin\left(\alpha\_k(t\_n - \tau)\right)\right]^2}{\sum\_{n=1}^N \sin^2(\alpha\_k(t\_n - \tau))} \right\} \tag{6}$$

where s and σ<sup>2</sup> are the mean and variance of the data and the value of τ is defined as

$$\tan\left(2o\_k\tau\right) = \frac{\sum\_{n=1}^{N} \sin\left(2o\_k t\_n\right)}{\sum\_{n=1}^{N} \cos\left(2o\_k t\_n\right)}.\tag{7}$$

#### 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 sð Þ τ<sup>n</sup> are the exact values of the signal at the points given by Eq. (2). The corresponding DFT is

$$S\_k = \sum\_{n=1}^{N} s(\tau\_n) \exp\left(-i\alpha \rho\_k \tau\_n\right). \tag{8}$$

Our aim is to find a good approximation to this expression in terms of the unevenly sampled signal s tð Þ<sup>n</sup> .

We start with the Euler summation formula:

$$\sum\_{n=1}^{N} f(\tau\_n) = \frac{1}{\Delta \tau} \int\_{\tau\_1}^{\tau\_N} f(\tau) d\tau + \frac{1}{2} [f(\tau\_1) + f(\tau\_N)] + \frac{\Delta \tau}{12} \left[ f'(\tau\_N) - f'(\tau\_1) \right] + O(\Delta \tau^2) \tag{9}$$

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

Fourier Transforms - Century of Digitalization and Increasing Expectations

$$\int\_{\tau\_1}^{\tau\_N} f(\tau)d\tau = \int\_{t\_1}^{t\_2} f(\tau)d\tau + \int\_{t\_2}^{t\_3} f(\tau)d\tau + \dots + \int\_{t\_{N-1}}^{t\_N} f(\tau)d\tau \tag{10}$$

an incomplete basis. It was found in Ref. [10] that for evenly sampled data with a sampling rate f <sup>S</sup>, the spectral amplitude S f ð Þ evaluated with FFT has the following

> ¼ S �f � f <sup>S</sup>

In order to avoid the aliasing symmetry of Eq. (17), the frequencies should be

<sup>f</sup> <sup>B</sup> <sup>¼</sup> <sup>f</sup> <sup>S</sup>

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

where f <sup>N</sup> is the Nyquist frequency for the RR intervals. As the ECG signal contains frequencies much greater than f <sup>N</sup> and the RR intervals are derived from the ECG signal, one cannot be sure that the spectral analysis of the RR intervals is free from aliasing. As a matter of fact, there are indications of aliasing in some rare cases [11–16]. One way to identify aliasing is to change the sampling rate and follow the changes in the spectrum. Unfortunately, for the RR intervals, one cannot speak about a definite sampling rate but rather can consider a distribution of sampling rates. The changes in sampling rate required to observe aliasing are of the same order as the fluctuations in the sampling rate. Therefore in practice it is almost impossible to observe consistent changes in the spectrum slightly changing the heart rate.

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

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

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

 ¼ S �f � nf <sup>S</sup> 

f <sup>S</sup> ¼ 1=Δτ ¼ 2f <sup>N</sup>, (19)

<sup>2</sup> : (18)

, (17)

symmetry properties:

rate from Eq. (1) in both cases:

basis to treat this problem.

4. Three experiments

43

spectrum and the spectrum of the ECG signal.

j j S f ð Þ ¼ S f � f <sup>S</sup>

The RR Interval Spectrum, the ECG Signal, and Aliasing

DOI: http://dx.doi.org/10.5772/intechopen.85327

where f is the frequency and n is an arbitrary integer.

bounded by the Nyquist frequency (denoted here by f <sup>B</sup>) according to

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

$$\int\_{\tau\_1}^{\tau\_N} f(\tau)d\tau = \frac{1}{2} [f(t\_1) + f(t\_2)](t\_2 - t\_1) + \dots + \frac{1}{2} [f(t\_{N-1}) + f(t\_N)](t\_N - t\_{N-1}) + O(\Delta \tau). \tag{11}$$

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

$$\begin{split} \sum\_{n=1}^{N} f(\mathbf{r}\_{n}) &= \frac{1}{2\Delta\tau} \{ [f(\mathbf{t}\_{1}) + f(\mathbf{t}\_{2})](\mathbf{t}\_{2} - \mathbf{t}\_{1}) + \dots + [f(\mathbf{t}\_{N-1}) + f(\mathbf{t}\_{N})](\mathbf{t}\_{N} - \mathbf{t}\_{N-1}) \} \\ &+ \frac{1}{2} [f(\mathbf{t}\_{1}) + f(\mathbf{t}\_{N})] + O(\Delta\tau). \end{split} \tag{12}$$

When tn are equally spaced, Eq. (12) becomes an identity with the Oð Þ¼ Δτ 0; therefore it seems to us that Eq. (12) is satisfied with an higher accuracy than just Oð Þ Δτ : Eq. (12) can be applied to approximate Eq. (8) with the substitution

$$f(t\_n) = \mathfrak{s}(t\_n) \exp\left(-i\alpha\_k t\_n\right),\tag{13}$$

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

$$S\_k = \sum\_{n=1}^{N} c\_n s(t\_n) \exp\left(-i\alpha\_k t\_n\right) + O(\Delta \tau),\tag{14}$$

where

$$\begin{aligned} c\_1 &= \frac{\Delta \tau + t\_2 - t\_1}{2\Delta \tau}, \\ c\_2 &= \frac{t\_3 - t\_1}{2\Delta \tau}, \\ \vdots \\ c\_{N-1} &= \frac{t\_N - t\_{N-2}}{2\Delta \tau}, \\ c\_N &= \frac{\Delta \tau + t\_N - t\_{N-1}}{2\Delta \tau}, \end{aligned} \tag{15}$$

with the inverse formula

$$s(\pi\_n) = \frac{1}{N} \sum\_{k=1}^{N} \mathcal{S}\_k \exp\left(i o y\_k \tau\_n\right) + \mathcal{O}(\Delta \tau), \quad o y\_k = 2\pi (k - 1)/(N\Delta \tau), \tag{16}$$

which is an interpolation formula for s tð Þ<sup>n</sup> at the evenly spaced points τ1⋯τN.

## 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 ð<sup>τ</sup><sup>N</sup> τ1

rule

ð<sup>τ</sup><sup>N</sup> τ1

> ∑ N n¼1

becomes

where

with the inverse formula

1 <sup>N</sup> <sup>∑</sup> N k¼1

sð Þ¼ τ<sup>n</sup>

3. Aliasing

42

<sup>f</sup>ð Þ<sup>τ</sup> <sup>d</sup><sup>τ</sup> <sup>¼</sup> <sup>1</sup>

fð Þ¼ τ<sup>n</sup>

2

1 2Δτ

þ 1 2 fð Þτ dτ ¼

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

ðt2 t1

Fourier Transforms - Century of Digitalization and Increasing Expectations

½ � f tð Þþ <sup>1</sup> f tð Þ<sup>2</sup> ð Þþ t<sup>2</sup> � t<sup>1</sup> ⋯ þ

½ f tð Þþ <sup>1</sup> f tð Þ <sup>N</sup> � þ Oð Þ Δτ :

Sk ¼ ∑ N n¼1 fð Þτ dτ þ

ðt3 t2

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

1 2

f g ½ � f tð Þþ <sup>1</sup> f tð Þ<sup>2</sup> ð Þþ t<sup>2</sup> � t<sup>1</sup> ⋯ þ ½ � f tð Þþ <sup>N</sup>�<sup>1</sup> f tð Þ <sup>N</sup> ð Þ tN � tN�<sup>1</sup>

f tð Þ¼ <sup>n</sup> s tð Þ<sup>n</sup> exp ð Þ �iωktn , (13)

cns tð Þ<sup>n</sup> exp ð Þþ �iωktn Oð Þ Δτ , (14)

When tn are equally spaced, Eq. (12) becomes an identity with the Oð Þ¼ Δτ 0; therefore it seems to us that Eq. (12) is satisfied with an higher accuracy than just Oð Þ Δτ : Eq. (12) can be applied to approximate Eq. (8) with the substitution

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

<sup>c</sup><sup>1</sup> <sup>¼</sup> <sup>Δ</sup><sup>τ</sup> <sup>þ</sup> <sup>t</sup><sup>2</sup> � <sup>t</sup><sup>1</sup>

cN�<sup>1</sup> <sup>¼</sup> tN � tN�<sup>2</sup>

cN <sup>¼</sup> <sup>Δ</sup><sup>τ</sup> <sup>þ</sup> tN � tN�<sup>1</sup>

which is an interpolation formula for s tð Þ<sup>n</sup> at the evenly spaced points τ1⋯τN.

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

<sup>c</sup><sup>2</sup> <sup>¼</sup> <sup>t</sup><sup>3</sup> � <sup>t</sup><sup>1</sup> <sup>2</sup>Δ<sup>τ</sup> ,

⋮

<sup>2</sup>Δ<sup>τ</sup> ,

<sup>2</sup>Δ<sup>τ</sup> ,

<sup>2</sup>Δ<sup>τ</sup> ,

Sk exp ð Þþ iωkτ<sup>n</sup> Oð Þ Δτ , ω<sup>k</sup> ¼ 2πð Þ k � 1 =ð Þ NΔτ , (16)

fð Þτ dτ þ ⋯ þ

ðtN tN�<sup>1</sup>

½ � f tð Þþ <sup>N</sup>�<sup>1</sup> f tð Þ <sup>N</sup> ð Þþ tN � tN�<sup>1</sup> Oð Þ Δτ :

fð Þτ dτ (10)

(11)

(12)

(15)

an incomplete basis. It was found in Ref. [10] that for evenly sampled data with a sampling rate f <sup>S</sup>, the spectral amplitude S f ð Þ evaluated with FFT has the following symmetry properties:

$$|\mathcal{S}(f)| = |\mathcal{S}(f \pm f\_{\mathcal{S}})| = |\mathcal{S}(-f \pm f\_{\mathcal{S}})| = |\mathcal{S}(\pm f \pm \eta f\_{\mathcal{S}})|,\tag{17}$$

where f is the frequency and n is an arbitrary integer.

In order to avoid the aliasing symmetry of Eq. (17), the frequencies should be bounded by the Nyquist frequency (denoted here by f <sup>B</sup>) according to

$$f\_B = \frac{f\_S}{2}.\tag{18}$$

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:

$$
\overline{f}\_{\mathcal{S}} = \mathbf{1}/\Delta \mathfrak{x} = \mathfrak{Y}\_{N^{\mathfrak{p}}} \tag{19}
$$

where f <sup>N</sup> is the Nyquist frequency for the RR intervals. As the ECG signal contains frequencies much greater than f <sup>N</sup> and the RR intervals are derived from the ECG signal, one cannot be sure that the spectral analysis of the RR intervals is free from aliasing. As a matter of fact, there are indications of aliasing in some rare cases [11–16]. One way to identify aliasing is to change the sampling rate and follow the changes in the spectrum. Unfortunately, for the RR intervals, one cannot speak about a definite sampling rate but rather can consider a distribution of sampling rates. The changes in sampling rate required to observe aliasing are of the same order as the fluctuations in the sampling rate. Therefore in practice it is almost impossible to observe consistent changes in the spectrum slightly changing the heart rate.

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.

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

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

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.

The relative power computed (from the ECG signal of Zahi) by four different methods, in the spectral range of 0– 0.5 Hz, (a) by DFT, (b) by spline interpolation of the RR data, (c) by NUDFT, and (d) from the ECG signal.

above range is presented in Figure 2d.

The RR Interval Spectrum, the ECG Signal, and Aliasing

DOI: http://dx.doi.org/10.5772/intechopen.85327

of aliasing in this region in the RR analysis.

Figure 2.

45

VLF is missing in J.C.'s ECG spectrum (see Figure 6a–6b).

## 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

Figure 1.

The relative power of the ECG signal of Zahi, (a) in the spectral range of 0–36 Hz and (b) in the spectral range of 0–12 Hz.

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 develop-

Fourier Transforms - Century of Digitalization and Increasing Expectations

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

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

The relative power of the ECG signal of Zahi, (a) in the spectral range of 0–36 Hz and (b) in the spectral range

intervals of 7-minute duration were chosen for analysis.

ment of aliasing.

Figure 1.

of 0–12 Hz.

44

4.1 The first experiment

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).

Figure 2.

The relative power computed (from the ECG signal of Zahi) by four different methods, in the spectral range of 0– 0.5 Hz, (a) by DFT, (b) by spline interpolation of the RR data, (c) by NUDFT, and (d) from the ECG signal.

Figure 3. The relative power of the ECG signal of Ori.

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

The relative power computed (from the ECG signal of Ofek) by two different methods, in the spectral range of

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

f <sup>1</sup> ¼ 1:234Hz; also its higher harmonics can be distinctively seen. The frequency f <sup>1</sup> is just the breathing frequency 74/min. In the same figure, one can also see the diffused peaks near the average heart rate frequency of 1.636 Hz and its higher

similar respiratory patterns exist in both the RR and in the ECG signals.

0–0.6 Hz, (a) by spline interpolation of the RR data, (b) from the ECG signal.

The RR Interval Spectrum, the ECG Signal, and Aliasing

DOI: http://dx.doi.org/10.5772/intechopen.85327

4.3 The third experiment

Figure 5.

47

Figure 4.

The relative power computed (from the ECG signal of Ori) by four different methods, in the spectral range of 0–0.52 Hz. (a) by DFT, (b) by spline interpolation of the RR data, (c) by NUDFT, (d) from the ECG signal.

#### 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

The RR Interval Spectrum, the ECG Signal, and Aliasing DOI: http://dx.doi.org/10.5772/intechopen.85327

#### Figure 5.

The relative power computed (from the ECG signal of Ofek) by two different methods, in the spectral range of 0–0.6 Hz, (a) by spline interpolation of the RR data, (b) from the ECG signal.

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 f <sup>1</sup> ¼ 1:234Hz; also its higher harmonics can be distinctively seen. The frequency f <sup>1</sup> is just the breathing frequency 74/min. In the same figure, one can also see the diffused peaks near the average heart rate frequency of 1.636 Hz and its higher

4.2 The second experiment

Figure 4.

46

Figure 3.

The relative power of the ECG signal of Ori.

Fourier Transforms - Century of Digitalization and Increasing Expectations

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

The relative power computed (from the ECG signal of Ori) by four different methods, in the spectral range of 0–0.52 Hz. (a) by DFT, (b) by spline interpolation of the RR data, (c) by NUDFT, (d) from the ECG signal.

Figure 6.

The relative power computed (from the ECG signal of J.C.) by two different methods, in the spectral range of 0– 0.46 Hz, (a) by spline interpolation of the RR data, (b) from the ECG signal.

Figure 7.

Figure 8.

49

The relative power computed (from the ECG signal of Ori with breathing rate of 0.04 Hz) by two different

methods, in the spectral range of 0–0.62 Hz, (a) by NUDFT, (b) from the ECG signal.

The RR Interval Spectrum, the ECG Signal, and Aliasing

DOI: http://dx.doi.org/10.5772/intechopen.85327

The relative power of the ECG signal of Alex with a breathing rate of 1.234 Hz.

harmonics. One should observe aliasing at about 1:636Hz� f <sup>1</sup> ¼ 0:402Hz. Indeed one can see diffused peaks around that frequency in Figure 9a, which displays the power estimation of the RR intervals using the NUDFT (which the Nyquist rate below is similar to the spline interpolation). The width of this region can be estimated by noting that the RR intervals have different instantaneous sampling rates which are equal to the inverse of the RR interval time. In Figure 10 we have calculated the distribution of the sampling rates by dividing the frequency region into 100 beans. We have shifted that distribution by subtracting f <sup>1</sup>. As one can see, the results are confined approximately to the region 0.32–0.47 Hz. Indeed the aliasing peaks of Figure 9a appear in this region. The pictures below the Nyquist frequency are very similar for the DFT, NUDFT, the spline interpolation, and the Lomb method (Figure 9b) with a similar aliasing behavior.

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 f <sup>1</sup>, as can be seen in Figure 9a and b. Both methods do not have the aliasing symmetry of the DFT as given by Eq. (17); therefore the results are not symmetric with respect to the Nyquist frequency (half the sampling rate), as it is satisfied, for example, in the case of the spline interpolation. We have found an

Figure 7.

harmonics. One should observe aliasing at about 1:636Hz� f <sup>1</sup> ¼ 0:402Hz. Indeed one can see diffused peaks around that frequency in Figure 9a, which displays the power estimation of the RR intervals using the NUDFT (which the Nyquist rate below is similar to the spline interpolation). The width of this region can be estimated by noting that the RR intervals have different instantaneous sampling rates which are equal to the inverse of the RR interval time. In Figure 10 we have calculated the distribution of the sampling rates by dividing the frequency region into 100 beans. We have shifted that distribution by subtracting f <sup>1</sup>. As one can see, the results are confined approximately to the region 0.32–0.47 Hz. Indeed the aliasing peaks of Figure 9a appear in this region. The pictures below the Nyquist frequency are very similar for the DFT, NUDFT, the spline interpolation, and the

The relative power computed (from the ECG signal of J.C.) by two different methods, in the spectral range of 0–

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 f <sup>1</sup>, as can be seen in Figure 9a and b. Both methods do not have the aliasing symmetry of the DFT as given by Eq. (17); therefore the results are not symmetric with respect to the Nyquist frequency (half the sampling rate), as it is satisfied, for example, in the case of the spline interpolation. We have found an

Lomb method (Figure 9b) with a similar aliasing behavior.

0.46 Hz, (a) by spline interpolation of the RR data, (b) from the ECG signal.

Fourier Transforms - Century of Digitalization and Increasing Expectations

Figure 6.

48

The relative power computed (from the ECG signal of Ori with breathing rate of 0.04 Hz) by two different methods, in the spectral range of 0–0.62 Hz, (a) by NUDFT, (b) from the ECG signal.

Figure 8. The relative power of the ECG signal of Alex with a breathing rate of 1.234 Hz.

Figure 11.

Figure 12.

51

The relative power of the ECG signal of Alex with a breathing rate of 1.404 Hz.

The RR Interval Spectrum, the ECG Signal, and Aliasing

DOI: http://dx.doi.org/10.5772/intechopen.85327

The relative power computed (from the ECG signal of Alex with a breathing rate of 1.404 Hz) by two different methods, in the spectral range of 0–1.6 Hz, (upper figure) by NUDFT, (lower figure) from the ECG signal.

Figure 9.

The relative power computed (from the ECG signal of Alex with a breathing rate of 1.234 Hz) by two different methods, in the spectral range of 0–1.5 Hz, (a) by NUDFT, (b) from the ECG signal.

Figure 10. A 100 bin histogram of the heart rates of Alex which are subtracted by the breathing rate of 1.234 Hz.

The RR Interval Spectrum, the ECG Signal, and Aliasing DOI: http://dx.doi.org/10.5772/intechopen.85327

Figure 11. The relative power of the ECG signal of Alex with a breathing rate of 1.404 Hz.

Figure 12.

The relative power computed (from the ECG signal of Alex with a breathing rate of 1.404 Hz) by two different methods, in the spectral range of 0–1.6 Hz, (upper figure) by NUDFT, (lower figure) from the ECG signal.

Figure 9.

Figure 10.

50

The relative power computed (from the ECG signal of Alex with a breathing rate of 1.234 Hz) by two different

A 100 bin histogram of the heart rates of Alex which are subtracted by the breathing rate of 1.234 Hz.

methods, in the spectral range of 0–1.5 Hz, (a) by NUDFT, (b) from the ECG signal.

Fourier Transforms - Century of Digitalization and Increasing Expectations

exact result at f <sup>1</sup> and a diffused aliasing around 0.4 Hz. It is interesting to note that both methods give almost the same result below and above the Nyquist frequency. One can interpret the appearance of the sharp peak at f <sup>1</sup> as a result of a partial destruction of aliasing symmetry due to uneven samplings.

4.In some cases the LF peak does not show up in the ECG spectrum. This led us

5. 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.

6.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.

7. 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

, Adi Ronen<sup>3</sup> and Yair Cassuto<sup>3</sup>

1 Department of Physics, Ben Gurion University of the Negev, Beer-Sheva, Israel

2 Department of Mathematics and Computer Sciences, Ben Gurion University of

3 Department of Life Sciences, Ben Gurion University of the Negev, Beer-Sheva,

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,

symmetry of the basis.

Author details

Israel

53

Alexander Gersten<sup>1</sup>

the Negev, Beer-Sheva, Israel

\*, Ori Gersten<sup>2</sup>

\*Address all correspondence to: alex.gersten@gmail.com

provided the original work is properly cited.

to suspect that part of the LF peak is of aliasing origin.

The RR Interval Spectrum, the ECG Signal, and Aliasing

DOI: http://dx.doi.org/10.5772/intechopen.85327

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–21].

## 6. Summary and conclusions

The ECG signal spectrum is bounded below the Nyquist frequency f<sup>B</sup> 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 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.

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:


exact result at f <sup>1</sup> and a diffused aliasing around 0.4 Hz. It is interesting to note that both methods give almost the same result below and above the Nyquist frequency. One can interpret the appearance of the sharp peak at f <sup>1</sup> as a result of a partial

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

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

The ECG signal spectrum is bounded below the Nyquist frequency f<sup>B</sup> 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 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

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

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

1. 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

2.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

3. 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

destruction of aliasing symmetry due to uneven samplings.

Fourier Transforms - Century of Digitalization and Increasing Expectations

speeded heart rate, and in constructing aliasing filters [18–21].

allowing to observe in detail the development of aliasing.

the external influences on the heart signal.

spectrum has its origin in aliasing.

sampled data.

52

5. Further developments

6. Summary and conclusions

analysis.


## Author details

Alexander Gersten<sup>1</sup> \*, Ori Gersten<sup>2</sup> , Adi Ronen<sup>3</sup> and Yair Cassuto<sup>3</sup>

1 Department of Physics, Ben Gurion University of the Negev, Beer-Sheva, Israel

2 Department of Mathematics and Computer Sciences, Ben Gurion University of the Negev, Beer-Sheva, Israel

3 Department of Life Sciences, Ben Gurion University of the Negev, Beer-Sheva, Israel

\*Address all correspondence to: alex.gersten@gmail.com

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## 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

[19] Keresnyei R. Selecting the optimal anti-aliasing filter for multichannel biosignal acquisition intended for intersignal phase shift analysis. Physiological

DOI: http://dx.doi.org/10.5772/intechopen.85327

The RR Interval Spectrum, the ECG Signal, and Aliasing

Measurement. 2015;36:N23-N24

2009

55

[20] J. Pučik, O. Ondraček Heart rate variability spectrum: Physiologic aliasing and nonstationarity

[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

considerations. In: Trends in Biomedical Engineering; Bratislava. 16–18 Sept.

[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

The RR Interval Spectrum, the ECG Signal, and Aliasing DOI: http://dx.doi.org/10.5772/intechopen.85327

[19] Keresnyei R. Selecting the optimal anti-aliasing filter for multichannel biosignal acquisition intended for intersignal phase shift analysis. Physiological Measurement. 2015;36:N23-N24

References

213-220

245-311

1995

1989

447-462

54

[1] Sayers BMA. Analysis of heart rate variability. Ergonomics. 1973;16:85-97

Fourier Transforms - Century of Digitalization and Increasing Expectations

[11] Witte H et al. Evidence of a previously undescribed form of

442-444

1989;20:1-12

physiologic signal.

Acta. 1990;49:59-68

563-572

410-419

Leuven; 2004

respiratory sinus arrhythmia (RSA)–the physiological manifestation of cardiac aliasing. Pflügers Archiv. 1988;412:

[12] Rother M et al. Cardiac aliasing–a possible cause for the misinterpretation

[13] Nilsson J, Panizza M, Hallett M. Principles of digital sampling of a

Electroencephalography and Clinical Neurophysiology. 1993;89:349-358

fluctuations in rabbits during different behavioral states. Biomedica Biochimica

physiological aliasing within the heart rate fluctuations by higher frequent respiratory movements. Journal of Physiology and Pharmacology. 1994;45:

[14] Zwiener U et al. Heart rate

[15] Zwiener U et al. Forms of

[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,

[18] Morren G. Advanced signal processing applied to in-vivo

spectroscopy and heart rate variability [thesis]. Katholieke Universiteit

Eggington S. A new method in applying power spectral statistics to examine cardio-respiratory interactions in fish. Journal of Theoretical Biology. 2006;241:

of cardiorespirographic data in neonates. Early Human Development.

[2] Akselrod SD et al. Power spectrum analysis of heart rate fluctuations: A quantitative probe of beat-to-beat cardiovascular control. Science. 1981;213:

[3] Kamath MV, Fallen EL. Power spectral analysis of HRV: A noninvasive

[4] Malik M et al. Heart rate variability,

physiological interpretation, and clinical use. Circulation. 1998;93:1043-1065

[5] Malik M, Camm AJ, editors. Heart Rate Variability. Armonk NY: Futura;

[6] Oppenheim AV, Schafer RW. Discrete-Time Signal Processing. Englewood Cliffs NJ: Prentice-Hall;

[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:

[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

signature of cardiac autonomic functions. Critical Reviews in Biomedical Engineering. 1993;21:

standards of measurements,

[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

Chapter 4

Abstract

Directional Denoising Using

Fourier filtering for image denoising consists in masking parts of the Fourier spectrum of an image and using inverse Fourier transform of the masked image to obtain the denoised one. In cases of directional noise, this process can induce artifacts, mainly because of the spatial coherence that exists in the theoretical noisefree image. Moreover, it can lead to loss of low-frequency content that is important in applications such as fringe projection technique, which aims at measuring 3D elevations of a surface. A method based on the principle of Fourier spectrum cloning for the denoising of images is proposed in this chapter. This method improves the PSNR and the SSIM ratio in comparison with spectrum masking denoising. The method will be detailed first, and then examples of image denoising

Fourier filtering is one of the main techniques used for the denoising of images corrupted by periodic noise. Most of image processing denoising algorithms tend to consider statistically defined noises, such as Gaussian, Poisson, or speckle noises [1]. However, in a relatively high number of cases, noises encountered in images are quasiperiodic and directional. These noises can be viewed as first-order, structured noises. Quasiperiodic noises are essentially due to AC perturbations or acquisition and reconstruction process errors in the case of three-dimensional reconstruction images. This type of noise is usually removed using filtering in the Fourier domain [2]. The Fourier transform is intrinsically well adapted, because it decomposes a signal on a basis of sine and cosine function which have an infinite support. The principle of Fourier filtering is usually the same: the Fourier spectrum exhibits some peaks that correspond to the frequencies of the noise, and the denoising operation consists in masking the part of the spectrum that contains the peaks after having detected them with the eye or simple or more complex algorithms. The removal of the Fourier spectrum peaks has a major drawback: the abrupt removal of all the Fourier coefficients induces artifacts and missing spatial frequencies in the

The idea developed in this chapter consists in cloning the values of the module of the spectrum around the removed part and to use a combination of these values to fill the removed part. The observed result is a reduced noise with fewer oscillations

Fourier Spectrum Cloning

Laurent Navarro and Jérôme Molimard

in two different applications will be presented.

spectrum cloning, periodic noise

1. Introduction

reconstructed images.

57

Keywords: Fourier transform, fringe projection, image denoising,

## Chapter 4
