# A hybrid approach to the simultaneous eliminating of power-line interference and associated ringing artifacts in electrocardiograms

- Xiaolin Zhou
^{1}Email author and - Yuanting Zhang
^{1, 2}

**12**:42

https://doi.org/10.1186/1475-925X-12-42

© Zhou and Zhang; licensee BioMed Central Ltd. 2013

**Received: **19 October 2012

**Accepted: **1 April 2013

**Published: **14 May 2013

## Abstract

### Background

The second-order, infinite impulse response notch filter is widely used to remove electrical power line noise in electrocardiograms (ECGs). However this filtering process often introduces spurious ringing artifacts in the vicinity of raw signal with sharp transitions. It is challenging to simultaneously remove these two types of noise without losing vital information about cardiac activities.

### Objective

Our objective is to devise a method to remove the power-line interference without introducing artifacts nor losing vital information. To this end we have developed the "hybrid approach" involving two-sided filtration and multi-iterative approximation techniques. The two-sided filtration technique can suppress the interference but some cardiac components are lost. The lost information can be restored using multi-iterative approximation technique.

### Results

For evaluation, four artificial data sets, each including 91 ECGs of different heart rates, were generated by a dynamical model. Four publicly-accessible sets of clinical data (MIT-BIH Arrhythmia, QT, PTB Diagnostic ECG, and T-Wave Alternans Challenge Databases) were also selected. Our new hybrid approach and the existing method were tested with these two types of signal under various pre-determined conditions. In contrast with the existing method, the hybrid approach can provide more than 27.40 dB and 37.77 dB reduction in signal distortion for 95% and 60% of artificial ECGs respectively; it can provide in excess of 11.78 dB and 17.48 dB reduction in distortion for 95% and 60% of these real records respectively.

### Conclusions

Overall, a significant reduction in signal distortion is demonstrated. These test results indicate that the newly proposed approach outperforms the traditional method assessed on both the artificial and clinical ECGs and suggest it could be of practical use for clinicians in the future.

## Background

Biopotential signals, such as electrocardiogram (ECG), often suffer from power-line interference (PLI, 50 or 60 Hz) since the recorded signal is an output of the electric fields of coupling states surrounding main power lines (PLs) and the power of the body. PLI is probably the most common problem encountered in the processing of biopotential signals. Essentially, a notch filter is adapted for minimizing PLI because of its ability to reject narrow band noise. Indeed, the second-order, infinite impulse response (IIR) notch filter is routinely applied for this purpose [1–3]. Because of the transient response effect of the notch filter, the impulse response of this type filter generally has an oscillatory behavior, which may cause microvolt-level ringing artifacts (RAs, typically ranging between 0 and 40 μV) in the immediate regions of input signal with sharp transitions. Besides, it will cause undesirable attenuation in signal components at frequencies close to the center frequency (50 or 60 Hz). Tolerable signal distortion needs a narrow stopband bandwidth (SBW); however, a narrower SBW results in a longer transient response time (TRT); whilst a longer TRT often incurs more serious RAs. It is an inherent contradiction. When an ECG signal is being processed, the RAs occur in the right side of QRS complexes, and consequently, this implies that many cardiac components are lost in ST-T regions. Serious distortion (signal distortion caused by the SBW itself and the appreciable RAs) may make the ECG signal more difficult to interpret, particularly for the ST-T segment analysis, QT interval estimation, the detection of Ventricular Late Potentials (VLPs) and so on [4–6]. Removal of PLI however should be done with utmost stringent efforts not to eliminate or distort the raw signals without introducing artifacts nor losing vital information [7–9]. Many sophisticated digital methods have been investigated to cope with either 50 or 60 Hz interference [10–12], and they satisfy the requirement for suppression and even elimination of PLI during ECG signals acquisition. However, it is impossible to design an IIR notch filter to remove PLI without causing distortion [13–15], and this problem is still unsolved in practice. In this paper we address the challenges to simultaneously remove the PLI and RAs without losing critical cardiac components by developing a new method which we call the "hybrid approach".

Being motivated by early pioneering work on investigating the RAs phenomena caused by the suppression of PLI, in this paper, the hybrid approach comprising two-sided filtration and multi-iterative approximation techniques, is proposed to simultaneously minimizing the PLI and associated RAs. In the first instance, the two-sided filtration technique is partitioned into four steps, which are applied to eliminate the PLI, localize the RAs and remove them afterward, whilst handling the boundary effects which are caused by the practical causal filter. To deal with the inherent contradictions, next the multi-iterative approximation technique is accomplished in three steps, which are adopted to sequentially reconstruct the lost cardiac information. A combination may thus prove to be more effective in eliminating these two types of noise. The elaborate scheme of the hybrid approach is stated in Sect. 3.

From the practical viewpoints of industrial and clinical applications, a bio-model-oriented diagnostic signal processing technique should be evaluated on the clinical data that are acquired from numerous subjects, as well as the artificial data that cover a variety of specific pre-determined conditions. Using artificial ECGs has a typical advantage, in that the signal distortion can be precisely calculated, since the ideal signal (i.e., "true" signal) can be reset to any desired case. In this study, the performances of the proposed and existing methods are evaluated in detail with artificial ECGs which are generated by an open-source program [16, 17], as well as four well-used, clinical ECG databases. All of these real data are accessible to the public at Physionet [18]. We compared these two methods with respect to signal distortion in the absence and presence of artificial PLs, respectively. We also examined noise reduction in the presence of artificial PLs. Specifically, these data sets are classified into four groups: with and without the addition of artificial electrical PLs under the notch filter with the center frequency of 50 and 60 Hz, respectively. The first two groups that without PLs, are designed for quantitatively investigating the signal distortion. The other two groups, which mixed with the artificial PLs, are chosen for evaluating the distortion in an environment with interference, together with examining the capacities of the newly presented and old methods for removing PLI. The relationship of SBW and TRT of notch filters has not been well delineated, to provide more insight in the next section we also detail the key properties of finite impulse response (FIR) and IIR notch filters, so that their properties can be compared with each other. Additionally, related RAs are quantitatively examined, and the challenges are outlined.

### Problem statement

*f*

_{0}. The magnitude response of an ideal notch filter may be given as below,

where *ω* = 2*πf*/*f*_{
s
} is the normalized digital frequency, *ω*_{0} = 2*πf*_{0}/*f*_{
s
} is the normalized center frequency at *f*_{0}. *f*_{
s
} is the sampling rate, and *f* is the specified frequency. In practice, the notch filter has a SBW at *f*_{0}, that is, Δ*ω* = 2*π* Δ*f*/*f*_{
s
}. Δ*ω* and Δ*f* are normalized and digital SBWs, respectively.

#### FIR and IIR notch filters

*H*

_{ f }(

*z*) denote the transfer function of a second-order, FIR notch filter,

*γ*= cos

*ω*

_{0}.

*H*

_{ f }(

*z*) is simple and easy to implement. However, a disadvantage in using this kind filter is that the SBW of

*H*

_{ f }(

*z*) is relatively large, which could not meet the specifications [19]. In order to be applicable at narrow SBW situations, a

*H*

_{ f }(

*z*) based second-order, IIR notch filter is then commonly used [3, 20],

*β*

_{0}= 1/(1 +

*λ*) and

*λ*= tan(Δ

*ω*/2). We regulate 0 <

*λ*< 1 in this study. Eq. (3.a) can be rearranged as follows,

*a*

_{0}= 1,

*b*

_{0}=

*b*

_{2}=

*β*

_{0},

*a*

_{1}=

*b*

_{1}= − 2

*γβ*

_{0}, and

*a*

_{2}= (1 −

*λ*)

*β*

_{0}. First provided

*γ*≠ − 1,

*H*

_{ i }(

*z*) contains two poles (

*α*

_{1}and

*α*

_{2}) inside the unit circle |

*z*| = 1 at

*z*complex plane,

where *α*_{1} + *α*_{2} = − *a*_{1}, *α*_{1} ⋅ *α*_{2} = *a*_{2}. The stability and settling time of *H*_{
i
}(*z*) are characterized by these two poles.

*x*[

*n*] and

*y*[

*n*] be the input and output signals at discrete time

*n*, respectively. This filter can be implemented by the following difference equation,

*k*refers to the

*k*th-order index of

*H*

_{ i }(

*z*). By deduction, in essence, Eq. (5.a) can be identified as follows,

*h*[

*k*] represents the impulse response of

*H*

_{ i }(

*z*) (see Appendix A),

*H*_{
i
}(*z*) is a stable system, since |*a*_{2}| < 1 and |*a*_{1}| < 1 + *a*_{2}. When *H*_{
i
}(*z*) contains a pair of complex valued poles or a single negative pole, *h*[*k*] will cycle back and forth between negative and positive during the transient state [21], which indicates that *h*[*k*] is associated with an oscillatory behavior.

*K*is a positive integer. One interpretation of Eq. (7) is that it represents a FIR notch system. Figure 1(a)-(d) display FIR and IIR notch filters calculated by Eqs. (6) and (3.b), respectively. Notably, the higher order

*K*of the filter, the weaker intensity the pass-band ripples and greater the selectivity. Consulting Figure 1(a)-(d), it is clear that this kind FIR filter has the following limitations: (A) The order

*K*is considerably higher than that of an equivalent second-order IIR filter meeting the same requirements. It thus has far more computational complexity. (B) Because of many pass-band ripples, signals that include the information of interest inside the relevant frequency bands will be grossly distorted. This is an issue related to the pseudo-Gibbs phenomena.

#### Inherent contradictions

*H*

_{ i }(

*z*) is often used for removing PLI. Figure 2(a)-(c) display an input clinical ECG that is corrupted by real PLI, the output of the input signal after passing it to

*H*

_{ i }(

*z*) and the relevant residue portion, respectively. In Figure 2(b) we see that PLI is well canceled. Ideally, PLI should be eliminated without any undesirable effects to distort the raw signals. Unfortunately, this cannot be achieved completely. RAs can sometimes interfere with the ST-T regions. The

*H*

_{ i }(

*z*) operates with IIR of oscillation, which probably distorts the input as well as attenuates the amplitude by producing RAs. Consequently, important cardiac components will be lost. In contrast to the ECG signal, RAs are not remarkable. In other words, they may pollute the ST segments (1 to 20 μV) [8], but they are not distinguishable by the naked eye unless the intensity of which is greater than a threshold, such as 1 μV. Just as in Figure 2(b), it seems as if there were no RAs. However, if we carefully check Figure 2(c) (as indicated by arrows), we find that they still exist.

Sharp transitions of signal may generate noticeable and intolerable RAs in the immediate vicinities of these abrupt changes. Figure 2(d)-(f) is another example, in which we clearly see the spurious effect of RAs: RAs with the amplitude up to 30 μV, as shown in Figure 2(f). Therefore, the RAs should be carefully removed to prevent the distortion of input signal. The main cause of RAs is due to the abrupt bandstop of *H*_{
i
}(*z*), spectral components that lie within the Δ*f*, as well as those close to *f*_{0} ± Δ*f*/2, will be attenuated; this is the frequency-domain description. In the time domain, the cause of RAs is *H*_{
i
}(*z*) itself: infinite impulse and oscillatory responses.

*n*

_{0}is a specific time. For simplifying the mathematics that follows the processing, by definition,

*δ*[

*n*] starts at 0 (

*n*

_{0}= 0), and goes to

*∞*. By making

*δ*[

*n*] pass through the system

*H*

_{ i }(

*z*), and letting

*y*

_{ f }[

*n*] be the output,

*ν*[

*n*] =

*δ*[

*n*] −

*y*

_{ f }[

*n*] the difference. The variance of

*ν*[

*n*] is then calculated by,

*λ*,

*γ*

^{2}+

*λ*

^{2}− 1 ≤ 0 (i.e., tan(Δ

*ω*/2) ≤ sin

*ω*

_{0}, commonly, Δ

*ω*/2 ≪

*ω*

_{0}), then the pole radius

*ρ*for

*H*

_{ i }(

*z*) is given by,

Referring to Eqs. (9.b) and (10), we come to the overall conclusions: (*i*) If *λ* ≪ 1, ${\sigma}_{v}^{2}\phantom{\rule{0.25em}{0ex}}$ can be loosely interpreted as a linear function of Δ*f*. With Δ*f* increasing, more and more signal components in the stop and pass bands will be modified in both the amplitude ("ripple") and the phase. (*ii*) Eq. (10) expresses that when Δ*f* → 0, then *ρ* → 1. In other words, the wider the Δ*f*, the closer the location of poles to the origin in the *z* plane, meaning the system *H*_{
i
}(*z*) settles more rapidly (also meaning the system has a shorter duration of TRT) [22]. Therefore, it indicates that the duration of RAs tends to decrease as Δ*f* increases. It is an inherent contradiction of notch filters.

*n*

_{0}= 1 s (amplitude of the impulse is 1 mV). By making Δ

*f*with an increment of 0.1 Hz from 0.5 to 10.0 Hz, we calculated the outputs

*y*

_{ f }[

*n*] by Eq. (5.a) at

*f*

_{0}= 50 Hz with each Δ

*f*. We thus obtained 96 outputs. For each output, at a specified Δ

*f*, the duration of RAs is defined as the time from impulse to the point

*w*, at which,

Figure 2(h) shows the output at Δ*f* = 3.0 Hz, and Figure 2(i) displays the corresponding residual components. Figures 2(j) and (k) plot all the calculated results. Observe that these results agree with the aforementioned conclusions.

### The principle of this algorithm

As we have seen, a cardinal implication is how the QRS complex affects the output of system *H*_{
i
}(*z*); to the system it always poses a big impulse. Because *H*_{
i
}(*z*) is a causal filter, noticeably time-decaying RAs can only occur on the right side of QRS complexes, see Figures 2(e), (f), (h) and (i). This implies that RAs only depend on the input waveforms regardless of whether the output waveforms which are of steep transitions or not. In addition, Eqs. (9.b) and (10) offer insights on the determinants of *H*_{
i
}(*z*), which are uniquely controlled by its Δ*f*. These key intrinsic properties of *H*_{
i
}(*z*) would be applied in this newly developed approach.

*x*[

*n*] denotes the input, the number of samples is

*L*and

*x*

^{ T }[

*n*] denotes the counterpart of the signal

*x*[

*n*], that is,

*x*

^{ T }[

*n*] =

*x*[

*L*− 1 −

*n*]. Then we construct a mirror extended signal,

We explore two-sided filtration and multi-iterative approximation techniques to eliminate the probable PLI and RAs which are contained in *x*_{
me
}[*n*]. For the various filter parmeters, let ${y}_{\mathit{me}}^{r}\left[n\right]\phantom{\rule{0.25em}{0ex}}$ denote outputs of the signal *x*_{
me
}[*n*] after passing it to the systems *H*_{
i
}(*z*), ${y}_{\mathit{me}}^{n}\left[n\right]\phantom{\rule{0.25em}{0ex}}$ denote outputs of this new method when applied to *x*_{
me
}[*n*].

#### Two-sided filtration technique

The design procedure can be summarized as follows:

Step 1 - Initialization

Given *f*_{
s
} and *f*_{0}, at a specified Δ*f*, we use Eq. (3.b) to calculate filter coefficients *a*_{1}, *a*_{2}, *b*_{0}, *b*_{1} and *b*_{2}. *bCons* = ⌊ *f*_{
s
}/*f*_{
b
} ⌋, ⌊·⌋ represents a round operator, *f*_{
b
} = 125 Hz is a constant. *bCons* is an integer that is chosen to accommodate different *f*_{
s
}. If *bCons* < 2, let*bCons* = 2.

Step 2 - Twice notch filtering

*i*) First notch filter suppressing. We pass the

*x*

_{ me }[

*n*] to the notch filter, that is the Eq. (5.a) with filter coefficients calculated in the previous step, and get the output ${y}_{\mathit{me}}^{r}\left[n\right]$ as the system output. Then we can obtain the differential components

*dy*

_{ me }[

*n*] using the derivative filter followed,

Figure 3(a) illustrates *y*_{
me
}^{
r
}[*n*] with the lost cardiac components, as indicated by the arrows. By means of signal-mirror extension, for a single filtering of *x*_{
me
}[*n*], it includes two operations: one is conducted in *x*[*n*] in the forward direction, see Figure 3(b), which shows the first half of *dy*_{
me
}[*n*] (0 *≤ n ≤ L* – 1 ); the other is conducted in the counterpart *x*^{T}[*n*], which is equivalent to be operated in *x*[*n*] in the backward direction, as displayed in Figure 3(c), which shows the next half of *dy*_{
me
}[*n*] (*L ≤ n ≤* 2*L* – 1) at relevant positions of the first half. The "relevant positions" means a sample at *n* is corresponded to the sample at *n** (*n** = 2*L* – 1 – *n*), hereinafter the same. Because *H*_{
i
} (*z*) is a causal system, detectable RAs in both halves of *dy*_{
me
}[*n*] lie on right and left arms of the same QRS complexes, respectively.

(*ii*) Second notch filter suppressing. Sequentially, let *dy*_{
me
}[*n*] be filtered by Eq. (5.a) with the same filter, we obtain the output $d{y}_{\mathit{me}}^{\text{'}}\left[n\right]$. *dy*_{
me
}[*n*] may involove the PLI, RAs and distorted components which are caused by the boundary effect [3]. However, $d{y}_{\mathit{me}}^{\text{'}}\left[n\right]$ can only contains the major information of RAs and distorted portions. Figure 3(d) shows the first half of $d{y}_{\mathit{me}}^{\text{'}}\left[n\right]$ and Figure 3(e) illuminates the second half of $d{y}_{\mathit{me}}^{\text{'}}\left[n\right]$. In order to localize RAs in the next step, we need the second notch filtering operation. In fact, the possible PLI is directly suppressed in this step.

Step 3 - RAs localization

*cs*

_{ me }[

*n*] implemented by the derivative operation as below,

*bCons*/2 samples. Next we let

*cs*

_{ me }[

*n*] pass through a low-pass filter shown as follows, and let

*ls*

_{ me }[

*n*] denote as the relevant output,

where *lOrder* = 4 · *bCons*, intrinsic delay of *H*_{l 0}(*z*) is (*lOrder* − 1)/2 samples and the gain is *lOrder*. For extracting the information of RAs, this low-pass filter is introduced to filter out fluctuations around the input sample and suppress undesirable outliers. Figure 3(f) displays the first half of *ls*_{
me
}[*n*] and Figure 3(g) shows the second half. From Figures 3(f) and (g), it is clear that triangle-like spikes localize the RAs.

Step 4 - Determination of the output

*ds*

_{ me }[

*n*], which results from the following operation,

*ds*

_{ me }[

*n*]. Due to the same consideration in

*Step 3*, we let

*ds*

_{ me }[

*n*] be processed by another low-pass filter defined as below,

*jOrder*= 16 ·

*bCons*, the phase delay of

*H*

_{l 1}(

*z*) is (

*jOrder*− 1)/2 samples and gain is

*jOrder*, and let

*js*

_{ me }[

*n*] denote the output. Figure 3(h) shows the first half of

*js*

_{ me }[

*n*], and we see Figures 3(b) and (c), it is obvious that to each sample at position

*n*: (1) if

*js*

_{ me }[

*n*] > 0, it represents that it gathers the information of RAs that lie within one arm of QRS complexes (left or right); (2) if

*js*

_{ me }[

*n*] < 0, it means that it gather the information of RAs that lie within the other arm of QRS complexes (right or left), see Figures 3(b), (c) and (h). Therefore, to each sample, we eliminate RAs based upon the criteria as below,

- (i)
If

*js*_{ me }[*n*] < 0, ${y}_{\mathit{me}}^{n}\left[n\right]={y}_{\mathit{me}}^{r}\left[n\right]+d{y}_{\mathit{me}}^{\text{'}}\left[n\right]$; - (ii)
If

*js*_{ me }[*n*] > 0, ${y}_{\mathit{me}}^{n}\left[n\right]={y}_{\mathit{me}}^{r}\left[{n}^{*}\right]+d{y}_{\mathit{me}}^{\text{'}}\left[{n}^{*}\right]$; - (iii)
If

*js*_{ me }[*n*] = 0 and*ls*_{ me }[*n*] <*ls*_{ me }[*n**],*y*^{ n }_{ me }[*n*] =*y*^{ r }_{ me }[*n*] +*dy*^{ ' }_{ me }[*n*]; - (iv)
Otherwise, ${y}_{\mathit{me}}^{n}\left[n\right]={y}_{\mathit{me}}^{r}\left[{n}^{*}\right]+d{y}_{\mathit{me}}^{\text{'}}\left[{n}^{*}\right].$

Furthermore, another benefit of criteria (*i*)-(*iv*) is avoiding the boundary effect. Figure 3(j) shows the output ${y}_{\mathit{me}}^{n}\left[n\right]\phantom{\rule{0.25em}{0ex}}$ processed by this new method, and Figure 3(i) shows the residue portion ${x}_{\mathit{me}}\left[n\right]-{y}_{\mathit{me}}^{n}\left[n\right]$. Pay special attention to the details of Figures 3(a) and (j), and see the regions with arrows in Figure 3(j), then consult Figures 3(b), (c) and (i), we easily find that RAs are mostly eliminated. Note worthily, variables *bCons*, *lOrder* and *jOrder* are self-adaptive with respect to *f*_{
s
}. Although the introduction of both-sided filtration increases the memory requirements and the proposed method also increases the complexity by introducing logical tests in *Step 4*, two-sided filtration restores the phase distortions caused by the one-sided operation.

#### Multi-iterative approximation technique

The single implementation of previous technique has two minor drawbacks, both of which are presented in Figure 3(i): (1) it is based upon the assumption that RAs consist of non-overlapping portions at positions of both directions of *x*[*n*]. In fact, it would not be able to distinguish overlapped RAs. However, *H*_{
i
}(*z*) of a small Δ*f* may incur some appreciable parts (greater than 1 μV) overlapped in the middle of adjacent QRS complexes, as shown by the second and forth arrows in Figure 3(i). (2) Another major limitation encountered in using such a practical filter is that the system will attenuate waveforms constituted by high frequencies surrounding *f*_{0}, as we see from positions with first and third arrows in Figure 3(i). As in the former disadvantage, intuitively, we might apply the *H*_{
i
}(*z*) with a larger Δ*f* to overcome this. From Eq (9.b), however, such a filter may cause more attenuation in signal components at frequencies close to *f*_{0}, since the duration of RAs and the amount of lost components are mutually dependent upon each other.

A novel method that may overcome the fundamental problems is to repeat the processing of two-sided filtration several times with various Δ*f*; we name it as multi-iterative approximation technique. Specifically, it can be partitioned into ternary steps as below,

Step I - Denoising with a large reference Δ*f*. In order to make distinctive RAs do not overlap each other, the *H*_{
i
}(*z*) with a large enough Δ*f* should be adapted. Whereas, it does not denote that the larger the Δ*f*, the better the performance. Consulting Eqs. (9.b) and (10) and Figures 2(j), (k), it is clear that duration of RAs tends to decrease drastically as Δ*f* increases at small Δ*f* cases; but it decreases slightly at large Δ*f* cases. Furthermore, ${\sigma}_{v}^{2}$ is roughly linear with respect to the Δ*f*. Thus, an appropriate compromise should be made to fit the task: an empirical and experimental Δ*f* = 6.0 Hz is employed here.

Input *x*_{
me
}[*n*] is depicted in Figure 2(d). Figure 4(b) illustrates the residual part ${x}_{\mathit{me}}\left[n\right]-{y}_{\mathit{me}}^{n}\phantom{\rule{0.12em}{0ex}}\left[n\right]$, which is obtained at Δ*f* = 6.0 Hz, and Figure 4(f) shows the related PSD spectrum. In the time domain, certain amounts of components of QRS complexes have been obviously lost, as we can see from Figure 4(b). It can likewise be seen in Figure 4(f) that parts of the signal of interest, which are close to *f*_{0}, have been suppressed because of this large Δ*f* in the frequency domain.

Step II - Reconstruction of details with a small target Δ*f*. In *Step I*, we achieve a short duration of RAs, and thus RAs can be eliminated. However, because of the large SBW Δ*f*, many cardiac components are lost. Thereby, in addition to possible PLs, residual part ${x}_{\mathit{me}}\left[n\right]-{y}_{\mathit{me}}^{n}\left[n\right]$ also contains lost components of original signal. In order to extract these useful details, ${x}_{\mathit{me}}\left[n\right]-{y}_{\mathit{me}}^{n}\left[n\right]$ is then to be processed by the two-sided filtration technique with a small Δ*f* (Δ*f* < 6.0 Hz). Let ${x}_{\mathit{me}}^{d}\left[n\right]={x}_{\mathit{me}}\left[n\right]-{y}_{\mathit{me}}^{n}\left[n\right]$ denote the input, ${y}_{\mathit{me}}^{d}\left[n\right]$ the output.

*Step I* and *Step II* are complementary with respect to Eqs (9.b) and (10). To illustrate, Figure 4(c) shows the ${y}_{\mathit{me}}^{d}\left[n\right]$ at Δ*f* = 2.0 Hz and Figure 4(g) displays the relevant PSD spectrum of ${y}_{\mathit{me}}^{d}\left[n\right]$. It can be observed from Figures 4(b), (c), (f) and (g) that most details have been reconstructed.

Step III - Reconstruction of slight details with the same target Δ*f*. We can see the vicinities indicated by arrows in Figure 4(g), ${y}_{\mathit{me}}^{d}\left[n\right]$ may still contain a little valuable information, just as the results shown in Figure 4(c): the amplitude of which is greater than 1 μV. To further reduce the distortion effects, similarly, let ${x}_{\mathit{me}}^{s}\left[n\right]={x}_{\mathit{me}}^{d}\left[n\right]-{y}_{\mathit{me}}^{d}\left[n\right]$, ${x}_{\mathit{me}}^{s}\left[n\right]$ is then be processed by the two-sided filtration technique with the same Δ*f* in *Step II*. Let ${y}_{\mathit{me}}^{s}\left[n\right]$ as the output. Figures 4(d) and (h) show ${y}_{\mathit{me}}^{s}\left[n\right]$ and the relevant spectrum, respectively. We see that spectrum content in pass bands of Figure 4(h) is flatter than that of Figure 4(g), this represents that most lost details have been reconstructed, since the amplitude of ${y}_{\mathit{me}}^{s}\left[n\right]$ is substantially less than 1 μV as shown in Figure 4(d).

*f*(i.e., [

*f*

_{0}− Δ

*f*/2,

*f*

_{0}+ Δ

*f*/2]); consequentially, these components can never be reconstructed due to the so-called "frequency overlap", as we can see from Figure 4(d). It is an inherent weakness of

*H*

_{ i }(

*z*). Figure 5 shows the flowcharts of the newly proposed algorithm; Figure 5(a) depicts the flowchart of two-sided filtration technique; Figure 5(b) represents the flowchart of multi-iterative approximation technique. Notably, Δ

*f*is the only parameter that needs to be specified in this method.

### Materials and evaluation

#### Artificial and clinical ECG data sets

We simulated four artificial ECG data sets with *f*_{
s
} of 250, 360, 500 and 1000 Hz, respectively. The generator generates realistic ECGs with user-settable parameters, such as "sampling frequency", "internal sampling frequency" and so on. This generator can be accessed from [16, 17]. For this study, we used the generator to simulate 10-second-duration data at "internal sampling frequency" of 2000 Hz (but 720 Hz for ECGs sampled at 360 Hz) with the specified "mean heart rates", and regulated other parameters with default values [17]. To assess the performance of these two methods in an environment with various sharp transitions, to each data set, we simulated ECGs with the "mean heart rates" ranging from 50 to 140 beats per minute (BPM), in increments of 1 BPM. We thus obtained 91 data for each data set.

**Four publicly-accessible sets of clinical data are selected for evaluation**

Databases |
f
(Hz) | Data numbers/Channels (Durations) | Brief description |
---|---|---|---|

QTDB | 250 | 105/2-lead | It was chosen to represent a wide variety of QRS and ST-T morphologies. |

(15 min.) | |||

MITDB | 360 | 48/2-lead | It was obtained from 47 subjects and contains affluent arrhythmia information. |

(30 min.) | |||

TWADB | 500 | 100/multi-lead | Including subjects with risk factors, such as myocardial infarctins, transient ischemia, ventricular and so on, as well as healthy subjects. |

(2 min.(appro.)) | |||

PTBDB | 1000 | 549/15-lead | It was collected from 290 healthy volunteers and patients. |

(1.92 min.(appro.)) |

#### Performance metrics

The main power supply is not perfectly stable. In some countries, tolerance of the frequency variation of PLs is 1% [23]. In poor electrical environments, however, the variation may up to 3% [24]. Because the bandwidth of real PLs is varying, in order to assess the performance of this method under different situations, to both clinical and artificial data, we conducted these tests by changing Δ*f* of the *H*_{
i
}(*z*) from 1.0 to 4.0 Hz with an increment of 0.1 Hz, respectively. In addition, to determine the capacity of eliminating PLI, artificial PLs that were simulated by sinusoidal functions (with the amplitude of 0.1 mV) at industrial frequencies (50 and 60 Hz, respectively), were added to the raw data. To each lead, at a specified Δ*f*, different situations can be categorized into four groups: ❶ At *f*_{0} = 50 Hz, implement *H*_{
i
}(*z*) without PLI; ❷ At *f*_{0} = 60 Hz, implement *H*_{
i
}(*z*) without PLI; ❸ At *f*_{0} = 50 Hz, implement *H*_{
i
}(*z*) with the addition of 50 Hz PLI; ➍ At *f*_{0} = 60 Hz, implement *H*_{
i
}(*z*) with the addition of 60 Hz PLI.

*H*

_{ i }(

*z*). To each input signal

*x*[

*n*], we obtained the outputs

*y*

^{ r }[

*n*] and

*y*

^{ n }[

*n*] by processing input signal

*x*[

*n*] with the IIR filter, that is Eq. (5.a) and this new method, respectively. We calculated the relevant ratio of percentage root-mean-square difference (

*rPRD*, in units of decibels [

*dB*]), and it is given by,

where $\mathit{PR}{D}_{r}^{2}=\frac{{\sigma}_{r}^{2}}{A{m}^{2}}$, $\mathit{PR}{D}_{n}^{2}=\frac{{\sigma}_{n}^{2}}{A{m}^{2}}$, $A{m}^{2}={\displaystyle \sum}_{n=0}^{\mathit{L}-1}x{\left[n\right]}^{2}$, ${\sigma}_{r}^{2}={\displaystyle \sum}_{n=0}^{\mathit{L}-1}{\left(x\left[n\right]-{y}^{r}\left[n\right]\right)}^{2}$ and ${\sigma}_{n}^{2}={\displaystyle \sum}_{n=0}^{\mathit{L}-1}{\left(x\left[n\right]-{y}^{n}\left[n\right]\right)}^{2}$. *rPRD* illustrates, for the new method, how distortion of the input *x*[*n*] is quantitatively lessened in comparison with the old method. It is not feasible to access the purely clinical signals, since these real signals might already contain PLs and other broadband noises. We hence let these signals be processed via the new approach first, and let the resulting signals be the inputs *x*[*n*] here.

To each data set with a specified group, wherein first, for the results *rPRD* s calculated from all of the various SBWs, we can figure out a threshold, at which, those results that, in excess of a certain percentage of total *rPRD* s, are greater than this threshold. Let *rPRD*|_{95%} and *rPRD*|_{60%} denote the thresholds with percentages of 95% and 60%, respectively. To facilitate comparisons with the old method in examining the capability of minimization of distortion at a specific SBW, then, for those results *rPRD* s outputted by the same SBW Δ*f*, at which (i.e., this single Δ*f*), we can define $\mathit{rPRD}{|}_{50\%}^{\mathit{\Delta f}}$ denote the threshold with a percentage of 50%. Therefore, *rPRD*|_{95%}, *rPRD*|_{60%} and $\mathit{rPRD}{|}_{50\mathit{\%}}^{\mathit{\Delta f}}$ represent the overall performances of this newly presented method.

## Results and discussion

*f*, each curve is associated to a histogram in the same column. Tables 2 and 3 give the statistical results (

*rPRD*|

_{95%}and

*rPRD*|

_{60%}) corresponding to the histograms in Figures 6 and 7, respectively.

**Statistical results of all of the four artificial ECG data sets**

f
(Hz) | Without PLI | Additive PLI (0.1 mV) | ||||||
---|---|---|---|---|---|---|---|---|

f | f | f | f | |||||

rPRD| | rPRD| | rPRD| | rPRD| | rPRD| | rPRD| | rPRD| | rPRD| | |

250 | 28.82 | 38.82 | 33.20 | 42.53 | 29.49 | 40.25 | 35.93 | 45.29 |

360 | 28.91 | 38.75 | 34.76 | 42.60 | 29.67 | 40.48 | 36.86 | 45.60 |

500 | 28.09 | 38.60 | 33.01 | 41.05 | 27.88 | 39.20 | 34.66 | 43.68 |

1000 | 27.40 | 37.77 | 32.70 | 41.19 | 27.62 | 38.12 | 33.78 | 42.69 |

**Statistical results of all of the four clinical ECG data sets**

f | Without PLI | Additive PLI (0.1 mV) | ||||||
---|---|---|---|---|---|---|---|---|

f | f | f | f | |||||

rPRD| | rPRD| | rPRD| | rPRD| | rPRD| | rPRD| | rPRD| | rPRD| | |

250 | 17.05 | 21.86 | 15.07 | 20.58 | 17.46 | 22.14 | 16.02 | 21.10 |

360 | 15.25 | 19.78 | 11.78 | 17.48 | 15.29 | 19.90 | 12.24 | 17.71 |

500 | 14.66 | 26.16 | 14.71 | 24.86 | 15.53 | 26.66 | 15.76 | 26.71 |

1000 | 14.67 | 24.20 | 15.88 | 23.85 | 16.58 | 25.38 | 18.07 | 26.07 |

Figure 6(a)-(e) and Figure 7(a)-(e) plot the results using the artificial data set and QTDB (*f*_{
s
} = 250 Hz), respectively. To groups ❶ and ❷, as Δ*f* increases, it tends to produce larger $\mathit{rPRD}{|}_{50\%}^{\mathit{\Delta f}}$ for both artificial and clinical data sets as a whole. By contrast, to both groups ❸ and ➍, as Δ*f* increases, it tends to produce $\mathit{rPRD}{|}_{50\%}^{\mathit{\Delta f}}$ with relationships that look like the exponential decay for artificial data sets, while relationships with $\mathit{rPRD}{|}_{50\%}^{\mathit{\Delta f}}$ decreasing first and then increasing for QTDB data as a whole. With regard to different notch frequencies, the performances are significantly different for both artificial and clinical ECG data sets, as we see from Figures 6(e) and 7(e). From Table 2, to four groups, the minimum *rPRD*|_{95%} and *rPRD*|_{60%} were 28.82 dB and 38.82 dB for artificial ECGs, respectively. Similarly, in all of four groups, the minimum *rPRD*|_{95%} and *rPRD*|_{60%} were 15.07 dB and 20.58 dB for QTDB ECGs as shown in Table 3, respectively. It is worth noting that QTDB data was chosen specifically to contain a broad variety of QRS morphologies (i.e., various kinds of abrupt discontinuities) [25]. At this point, the results of QTDB ECGs are relatively more objective in the present study.

Figures 6(f)-(i) and 7(f)-(i) display the results using the artificial data set and MITDB (*f*_{
s
} = 360 Hz), respectively. Figures 6(j) and 7(j) are the relevant $\mathit{rPRD}{|}_{50\%}^{\mathit{\Delta f}}$ relations of Figures 6(f)-(i) and 7(f)-(i), respectively. MITDB includes in excess of 53 leads with abnormal, wide size or low intensity QRS waveforms, which are of low frequency components (i.e., slow transitions or smooth variations). Therefore, this turns out that less distortion results from the old method, since RAs always occur surrounding the notch frequency *f*_{0}, that is the high-frequency end in most ECG spectra. Hence, the performance of the new method is not very good for MITDB data, as illustrated in Figure 7(j). Likewise, with respect to the different notch frequencies, the performances differ markedly for both artificial and clinical ECG data sets, as we can see from Figures 6(j) and 7(j). However, for the same notch frequency (50 or 60 Hz, respectively), the performances are not significantly different for MITDB data without or with the addition of artificial PLs. Of these statistical results, as displayed in Tables 2 and 3, for four groups, the minimum *rPRD*|_{95%} and *rPRD*|_{60%} were 28.91 dB and 38.75 dB for artificial ECGs, respectively; for the QTDB recordings with four groups, the minimum *rPRD*|_{95%} and *rPRD*|_{60%} were 11.78 dB and 17.48 dB, respectively.

The statistical results of artificial data set and TWADB (*f*_{
s
} = 500 Hz) are shown in Figures 6(k)-(p) and 7(k)-(p), respectively. In terms of artificial ECG data with four groups, as Δ*f* increases, they exhibit consistent tendencies with these of artificial ECGs sampled at 250 and 360 Hz. However, for the TWADB data with four groups, as Δ*f* increases, the performances depict no significant consistency as those of QTDB and MITDB, since the $\mathit{rPRD}{|}_{50\%}^{\mathrm{\Delta f}}$ are divergent at sides of low and large SBWs, while relatively convergent in the middle of SBWs. In addition, for the QTDB and MITDB data sets, results with high probabilities lie in low sides of SBWs (i.e., less than 22 dB), but for the TWADB data the results of high probabilities lie within high sides (i.e., greater than 30 dB). In Tables 2 and 3, for the four groups, we see that the minimum *rPRD*|_{95%} and *rPRD*|_{60%} were 27.88 dB and 38.60 dB, 14.66 dB and 24.86 dB, for artificial and TWADB data sets, respectively. Furthermore, from Figures 6(e), (j), (p) and Table 3, for this clinical data set, the performance of the new method is better than that of the QTDB and MITDB data sets as a whole.

It can be clearly seen from Figures 6(q)-(u) and 7(q)-(u) that, for the artificial data set and PTBDB (*f*_{
s
} = 1000 Hz), respectively, the results of the artificial ECGs reveal no significant difference compared with three earlier used artificial data sets. However, the clinical ECGs have significant differences in the histograms as well as the $\mathit{rPRD}{|}_{50\%}^{\mathit{\Delta f}}$ relations to these of three previous clinical ECGs. As far as the histograms of the four groups are concerned, they exhibit Gaussian probability distributions, since a huge number of results was obtained (a total of 204228 results) for each group, as shown in Figure 7(q)-(t). In terms of $\mathit{rPRD}{|}_{50\%}^{\mathit{\Delta f}}$ relations, as shown in Figure 7(u), four curves show tendencies toward convergence as Δ*f* increases. Again from Tables 2 and 3, the minimum *rPRD*|_{95%} and *rPRD*|_{60%} were 27.40 dB and 37.77 dB for artificial ECGs of four groups, respectively; and for the PTBDB recordings, the minimum *rPRD*|_{95%} and *rPRD*|_{60%} were 14.67 dB and 23.85 dB, respectively.

Generally, for the artificial ECGs of a specific group, differences of results which were obtained from different data sets, show no special significance, since the only differences of each data set are the sampling rates *f*_{
s
} and the duration *L*. Regarding the four clinical ECG data sets (each including four groups), however, they showed obviously different performances, since these data sets were obtained from various subjects and with different emphases. In particular, the MITDB recordings, comprising many and a broad variety of wide size QRS complexes, show relatively poor performance. From all the results of artificial and clinical ECGs sampled at each *f*_{
s
}, we can find that the performance exhibited by the artificial data set is significantly better than that of the relevant clinical data set. The reason is that the clinical ECGs may contain many low frequency leads. In summary, all of these test results indicate that the proposed method has the capacity to well reduce the PLI, and simultaneously, greatly minimize the RAs for both artificial and clinical ECGs.

### Benefits and limitations

As previously mentioned, we always want the SBW of a notch filter to be very narrow for suppressing PLI. The FIR notch filter, which is defined by Eq. (7), does not encounter the "infinite impulse" problem, but it requires a large degree *K* to meet the specification, and this often comes at the cost of high computational complexity. Additionally, another problem with the FIR notch filter is that it may still produce RAs since it is oscillatory in the range of the "finite response", because *K* is a large number. In contrast, our hybrid approach provides primary benefits in eliminating specific interferences. It is not limited to reducing fixed fundamental PLI (50 or 60 Hz) but also applicable to removing the high-frequency harmonics for some worse cases, since the *f*_{0} can be tailored to any desired frequency in this approach.

Thus far, all of our discussions are based upon the tan(Δ*ω*/2) ≤ sin *ω*_{0}. In certain applications, however, some rare cases may still exist where tan(Δ*ω*/2) > sin *ω*_{0}, especially those ECG monitor systems with low sampling rates *f*_{
s
}. In general, such situations occur at *f*_{
s
} = 2*f*_{0} + *δ*_{
f
} and 0 ≤ *δ*_{
f
} ≪ *f*_{
s
}. By Eq. (4), consider the following two possible cases,

- (i)
*γ*= − 1 That is,*f*_{0}=*f*_{ s }/2. According to Eq. (3.a),*H*_{ i }(*z*) represents a first-order, IIR notch filter that has only one pole*α*inside the unit circle, for this case,${\rho}^{\text{'}}=\left|\alpha \right|=-\alpha ={a}_{2}=\frac{1-\lambda}{1+\lambda},\left(\alpha <0,{\rho}^{\text{'}}\propto 1/\mathrm{\Delta}f\right)$(20)

*δ*[

*n*], we can obtain the output variance,

*f*

_{0}will pollute valuable information of low frequencies. It is a limitation not caused by the method but

*f*

_{ s }itself. Therefore, we recommend

*δ*

_{ f }≥ 4.0 Hz for practical application based upon the current study and literature [23, 24].

- (ii)− 1 <
*γ*< 0 and*γ*^{2}+*λ*^{2}− 1 > 0 That is, |*α*_{1}| ≠ |*α*_{2}|. This yields (see Appendix C),$\begin{array}{c}\hfill \mathrm{Max}\left(\left|{\alpha}_{1}\right|,\left|{\alpha}_{2}\right|\right)=\left|\gamma \right|\cdot \frac{1+\sqrt{1-{\gamma}^{-2}+{\gamma}^{-2}{\lambda}^{2}}}{1+\lambda},\hfill \\ \hfill \left({\alpha}_{1}<0,{\alpha}_{2}<0,\mathrm{Max}\left(\left|{\alpha}_{1}\right|,\left|{\alpha}_{2}\right|\right)\propto \mathrm{\Delta}f\right)\hfill \end{array}$(22)

Max(|*α*_{1}|, |*α*_{2}|) denotes the larger one of |*α*_{1}| and |*α*_{2}|. For |*α*_{1}| ≠ |*α*_{2}|, the settling time of *H*_{
i
}(*z*) is mainly determined by Max(|*α*_{1}|, |*α*_{2}|) [22]. Within this situation, using the *H*_{
i
}(*z*) of a larger Δ*f* would not be able to achieve shorter durations of RAs but with more cardiac components lost. The "larger" reference Δ*f* is then set to the specified target Δ*f* in *Step I* to meet the uniformity within the approach. Thus, it is worth emphasizing that more precautions should be taken when Δ*f* is adjusted to avoiding overlapping RAs in the middle of adjacent QRS complexes. Indeed, it is essential to remember this limitation as well, since Eq. (10) is not universal but with conditions.

## Conclusions

Instead of frequency response, specifications of second-order, IIR notch filter may be given in terms of the impulse response. In this study, we started with the impulse response of *H*_{
i
}(*z*), the output ${\sigma}_{v}^{2}$ and the settling time that, concerning its behaviour in the time domain, have been quantitatively investigated. Minimizing the RAs involved in the signal of sharp transitions by removing PLI is a great concern in the processing of biopotential signal. The detection and analysis of the VLPs in the ECG signal is highly sensitive to the residual PLI and the RAs after the QRS complexes as a result of the filtering technique applied. The artifacts may become considerable in cases of high and steep complexes. Owing to the intrinsic properties of *H*_{
i
}(*z*), we proposed a hybrid approach, which utilizes *H*_{
i
}(*z*), to reliably suppress PLI as well as to aviod the generation of RAs. It is applicable to different *f*_{
s
} and easy to implement. In fact, Δ*f* is the only parameter that needs to be specified for this approach. Problems are greatly mitigated via these techniques. Sufficient results and performance statistics are provided to validate the reliability of this method in the test environment with a variety of conditions (e.g., artificial and clinical ECGs, of various *f*_{
s
}, with altering Δ*f*, etc.). An eventual consideration related to practice is *f*_{0} of PLs. To the artificial PLs used in this study, *f*_{0} was set to 50 or 60 Hz without varying. However, *f*_{0} of real PLs, similar to its bandwidth, may fluctuate over a small range. Even so, we can refer to many previous studies for how to adaptive tracking of the *f*_{0} with serious drift.

## Appendix A: Proof of the Equation (6)

*x*[

*n*] can be expressed as the sum of linear superposition of unit impulse functions

*δ*[

*n*−

*k*],

*k*=

*n*,

*δ*[

*n*−

*k*] = 1; otherwise, if

*k*≠

*n*,

*δ*[

*n*−

*k*] = 0. By definition, to each impulse function

*δ*[

*n*−

*k*], the output is denoted as the impulse response

*h*[

*n*−

*k*]. Consulting Eq. (5.a), we obtain,

*H*

_{ i }(

*z*) is a causal system, then we conclude that,

- (i)If
*n*= 0,*h*[*n*− 1] =*h*[*n*− 2] = 0. Ignore*δ*[*n*− 1] and*δ*[*n*− 2], since*δ*[*n*− 1] =*δ*[*n*− 2] = 0, then,$\phantom{\rule{4em}{0ex}}h\left[n\right]={b}_{0}$(a.3) - (ii)If
*n*= 1, because*h*[*n*− 2] = 0, and*δ*[*n*] =*δ*[*n*− 2] = 0, which results in,$\phantom{\rule{4em}{0ex}}h\left[n\right]=\left(-{a}_{1}\right)h\left[n-1\right]+{b}_{1}$(a.4) - (iii)If
*n*= 2,*δ*[*n*] =*δ*[*n*− 1] = 0, we have,$\phantom{\rule{4em}{0ex}}h\left[n\right]=\left(-{a}_{1}\right)h\left[n-1\right]+\left(-{a}_{2}\right)h\left[n-2\right]+{b}_{2}$(a.5) - (iv)If
*n*≥ 3,*δ*[*n*] =*δ*[*n*− 1] =*δ*[*n*− 2] = 0, this yields,$\phantom{\rule{4em}{0ex}}h\left[n\right]=\left(-{a}_{1}\right)h\left[n-1\right]+\left(-{a}_{2}\right)h\left[n-2\right]$(a.6)

We thus obtain Eqs (5.b) and (6) in terms of the filter coefficients *a*_{0}, *a*_{1}, *a*_{2}, *b*_{0}, *b*_{1} and *b*_{2}.

## Appendix B: Proof of the Equation (9.b)

*v*[

*n*] =

*δ*[

*n*] −

*y*

_{ f }[

*n*] is equivalent to the input

*δ*[

*n*] to be processed by the following comb filter,

where *c*_{0} = *λβ*_{0}. If *γ* ≠ − 1, *H*_{
c
}(*z*) also contains two poles at *α*_{1} and *α*_{2} inside the unit circle |*z*| at *z* plane.

_{c}represents the integral taken around the unit circle in counter-clockwise direction, and let

## Appendix C: Proof related to the Equation (22)

*λ*, which results in,

*λ*) ≤ 0 is true, from Eq. (c.5), this yields,

*λ*)

^{2}≤ 0. Due to

*λ*> 0, thus the assumption ℱ (

*λ*) ≤ 0 is false, then ℱ (

*λ*) > 0. Because

*G*(

*λ*) > 0, we therefore have,

## Declarations

### Acknowledgements

This work was supported in part by the National Basic Research Program 973 (2010CB732606), the Guangdong Innovation Research Team Fund for Low-Cost Healthcare Technologies in China, the External Cooperation Program of the Chinese Academy of Sciences (GJHZ1212), the Key Lab for Health Informatics of Chinese Academy of Sciences, the Enhancing Program of Key Laboratories of Shenzhen City (ZDSY20120617113021359) and the Young Scientists Fund of the National Science Foundation of China (61002002). The authors would like to thank Prof. Emma Pickwell-MacPherson and Dr. Yan Chen for their assistance, advice, and encouragement.

## Authors’ Affiliations

## References

- McManus CD, Neubert KD, Cramer E: Characterization and elimination of AC noise in electrocardiograms: a comparison of digital filtering methods.
*Comput Biomed Res*1993, 26: 48–67. 10.1006/cbmr.1993.1003View ArticleGoogle Scholar - Dotsinsky I, Stoyanov T: Power-line interference cancellation in ECG signals.
*Biomed Instrum Technol*2005, 39(2):155–162.Google Scholar - Pei SC, Tseng CC: Elimination of AC interference in electrocardiogram using IIR notch filter with transient suppression.
*IEEE Trans Biomed Eng*1995, 42(11):1128–1132. 10.1109/10.469385View ArticleGoogle Scholar - Breithardt G, Cain ME, El-Sherif N, Flowers NC, Hombach V, Janse M, Simson MB, Steinbeck G: Standards for analysis of ventricular late potentials using high-resolution or signal-averaged electrocardiography: a statement by a task force committee of the European Society of Cardiology, the American Heart Association, and the American College of Cardiology.
*J Am Coll Cardiol*1991, 17(5):999–1006. 10.1016/0735-1097(91)90822-QView ArticleGoogle Scholar - Santangeli P, Infusino F, Sgueglia GA, Sestito A, Lanza GA: Ventricular late potentials: a critical overview and current applications.
*J Electrocardiol*2008, 41(4):318–324. 10.1016/j.jelectrocard.2008.03.001View ArticleGoogle Scholar - Subramanian AS, Gurusamy G, Selvakumar G: Detection of Ventricular Late Potentials using Wavelet-Neural Approach.
*Eur J Sci Res*2011, 58: 11–20.Google Scholar - Kuchar DL, Thorburn CW, Sammel NL: Late potentials detected after myocardial infarction: natural history and prognostic significance.
*Circulation*1986, 74(6):1280–1289. 10.1161/01.CIR.74.6.1280View ArticleGoogle Scholar - Prutchi D, Norris M:
*Design and development of medical electronic instrumentation*. New Jersey: John Wiley E Sons, Inc; 2005.Google Scholar - Mousa A, Yilmaz A: Comparative analysis on wavelet-based detection of finite duration low-amplitude signals related to ventricular late potentials.
*Physiol Meas*2004, 25: 1443–1457. 10.1088/0967-3334/25/6/010View ArticleGoogle Scholar - Van Alstk JA, Scbilder TS: Removal of base-line wander and power-he interference from the ECG by an efficient FIR filter with a reduced number of taps.
*IEEE Trans Biomed Eng*1985, 32(12):1052–1060.View ArticleGoogle Scholar - Mitov IP: A method for reduction of power line interference in the ECG.
*Med Eng Phys*2004, 26: 879–887. 10.1016/j.medengphy.2004.08.014View ArticleGoogle Scholar - Mahesh SC, Aggarwala RA, Uplane MD: Interference reduction in ECG using digital FIR filters based on Rectangular window.
*WSEAS Trans Sig Proc*2008, 4(5):340–349.Google Scholar - Rangayyan RM:
*Biomedical signal analysis: a case-study approach*. New York: IEEE press; 2002.Google Scholar - Kaur M, Singh B: Powerline Interference Reduction in ECG Using Combination of MA Method and IIR Notch.
*IJRTE*2009, 2(6):125–129.Google Scholar - Luo S, Johnston P: A review of electrocardiogram filtering.
*J Electrocardiol*2010, 43(6):486–496. 10.1016/j.jelectrocard.2010.07.007View ArticleGoogle Scholar - McSharry PE, Clifford GD, Tarassenko L, Smith LA: A dynamical model for generating synthetic electrocardiogram signals.
*IEEE Trans Biomed Eng*2003, 50(3):289–294. 10.1109/TBME.2003.808805View ArticleGoogle Scholar - ECGSYN: A realistic ECG waveform generator. http://physionet.cps.unizar.es/physiotools/ecgsyn/Matlab/
- The Physionet ECG Database http://physionet.org/physiobank/database/
- Olguin DO, Lara FB, Chapa SOM: Adaptive notch filter for EEG signals based on the LMS algorithm with variable step-size parameter. In
*Proceedings of the 39th International Conference on Information Sciences and Systems*. Baltimore; 2005.Google Scholar - Hirano K, Nishimura S, Mitra S: Design of digital notch filters.
*IEEE Trans Commun*1974, 22(7):964–970. 10.1109/TCOM.1974.1092311View ArticleGoogle Scholar - Hellerstein J, Diao Y, Parekh S, Tilbury DM:
*Feedback control of computing systems*. Chichester: Wiley-IEEE Press; 2004.View ArticleGoogle Scholar - Vegte JVD:
*Fundamentals of Digital Signal Processing*. New York: Prentice Hall; 2003.Google Scholar - Levkov C, Mihov G, Ivanov R, Daskalov I, Christov I, Dotsinsky I: Removal of power-line interference from the ECG: a review of the subtraction procedure.
*Biomed Eng Online*2005, 4: 1–18. 10.1186/1475-925X-4-1View ArticleGoogle Scholar - Kumaravel N, Senthil A, Sridhar KS, Nithiyanandam N: Integrating the ECG power-line interference removal methods with rule-based system.
*Biomed Sci Instrum*1995, 31: 115–135.Google Scholar - Laguna P, Mark RG, Goldberg A, Moody GB: A database for evaluation of algorithms for measurement of QT and other waveform intervals in the ECG.
*Comput Cardiol*1997, 24: 673–676.Google Scholar - Press WH, Flannery BP, Teukolsky SA, Vetterling WT:
*Numerical recipes in C*. Melbourne: Cambridge University Press; 1988.Google Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.