Simulation of Adaptive Noise Canceller for an ECG signal Analysis

ACEEE Int. J. on Signal & Image Processing, Vol. 03, No. 01, Jan 2012

Simulation of Adaptive Noise Canceller for an ECG
signal Analysis
Dr. D.C. Dhubkarya1, Aastha Katara2, and Raj Kumar Thenua2
1
Department of Electronics & Communication, BIET, Jhansi, India
Email: dcd3580@yahoo.com
2
Department of Electronics & Instrumentation, Anand Engineering College, Agra, India
Email: aastha.katara@gmail.com, kumarraj04in@gmail.com

Abstract— In numerous applications of signal processing, These are preferable because they are stable, and no special
communications and biomedical we are faced with the adjustments are needed for their implementation.Fig.1
necessity to remove noise and distortion from the signals. illustrates the general configuration for an Adaptive filter [4].
Adaptive filtering is one of the most important areas in digital
The adaptive filter has two inputs: the primary input d(n),
signal processing to remove background noise and distortion.
In last few years various adaptive algorithms are developed
which represents the desired signal corrupted with undesired
for noise cancellation. In this paper we have presented an noise, and the reference signal x(n), which is the undesired
implementation of LMS (Least M ean Square), NLM S noise to be filtered out of the system.
(Normalized Least Mean Square) and RLS (Recursive Least The goal of adaptive filtering systems is to reduce the
Square) algorithms on MATLAB platform with the intention noise portion, and to obtain the uncorrupted desired signal.
to compare their performance in noise cancellation application. In order to achieve this task, a reference of the noise signal is
We simulate the adaptive filter in MATLAB with a noisy ECG needed. That reference is fed to the system, and it is called a
signal and analyze the performance of algorithms in terms of reference signal x(n). However, the reference signal is typically
M SE (Mean Squared Error), SNR Improvement,
not the same signal as the noise portion of the primary signal
computational complexity and stability. The obtained results
shows that, the RLS algorithm eliminates more noise from
- it can vary in amplitude, phase or time delay. Therefore the
noisy ECG signal and has the best performance but at the cost reference signal cannot be simply subtracted from the primary
of large computational complexity and higher memory signal to obtain the desired portion at the output.
requirements. The basic idea for the adaptive filter is to predict the
amount of noise in the primary signal, and then subtract that
Index Terms— Adaptive filters, LMS, Mean Squared Error noise from it. The prediction is based on filtering the reference
(MSE), RLS signal x(n), which contains a solid reference of the noise
present in the primary signal. The noise in the reference signal
I. INTRODUCTION is filtered to compensate for the amplitude, phase and time
An adaptive filter has the property of self-modifying its delay, and then subtracted from the primary signal. This
frequency response to change the behavior in time, allowing filtered noise is the system’s prediction of the noise portion
the filter to adapt the response to the input signal of the primary signal, y(n). The resulting signal is called error
characteristics change. Due to this capability, the overall signal e(n), and it presents the output of the system. Ideally,
performance and the construction flexibility, the adaptive the resulting error signal would be only the desired portion
filters have been employed in many different applications, of the primary signal.
some of the most important are: telephonic echo cancellation, In this work we investigate the performance of various
radar signal processing, navigation systems, communications adaptive algorithms with the help of MATLAB simulation [7]
channel equalization and biomedical signals processing [1- and tested for an ECG signal. The paper is organized in four
3].The most common adaptive filters, which are used during sections; section 2 gives an idea of adaptive algorithms, in
the adaption process, are the finite impulse response (FIR) section 3 an Adaptive Noise Cancellation (ANC) model is
types. designed and finally the results are discussed in section 4.

II. ADAPTIVE ALGORITHMS
A. LEAST MEAN SQUARE A LGORITHM
The LMS algorithm [4], is a type of adaptive filter algorithm
that is also known as stochastic gradient-based algorithm as
it utilizes the gradient vector of the filter tap weights to
converge on the optimal wiener solution. With each iteration
of the LMS algorithm, the filter tap weights of the adaptive
filter are updated according to the following formula:
Figure 1. General Adaptive filter configuration

© 2012 ACEEE 1
DOI: 01.IJSIP.03.01. 47


Here x(n) is the input vector of time delayed input values,

In the RLS algorithm the estimate of previous samples of
The vector
output signal, error signal and filter weight is required that
represents the coefficients of the adaptive FIR filter tap weight
leads to higher memory requirements.
vector at time n.
The parameter µ is known as the step size parameter and
III. ADAPTIVE NOISE CANCELLATION
is a small positive constant. This step size parameter controls
the influence of the updating factor. Selection of a suitable Adaptive noise cancellation (ANC) is performed by
value for µ is imperative to the performance of the LMS subtracting noise from a received signal, and an operation
algorithm, if the value is too small the time the adaptive filter controlled in an adaptive manner is done during the
takes to converge on the optimal solution will be too long; if adaptation process to get an improved signal-to-noise ratio.
µ is too large the adaptive filter becomes unstable and its Noise subtraction from a received signal could generate
output diverges. disastrous results by causing an increase in the average
power of the output noise. However when filtering and
B.NORMALIZED LEAST MEAN SQUARED A LGORITHM
subtraction are controlled by an adaptive process, it is
In the standard LMS algorithm, when the convergence possible to achieve a superior system performance compared
factor µ is large, the algorithm experiences a gradient noise to direct filtering of the received signal. Fig.2 shows adaptive
amplification problem. In order to solve this difficulty, we can noise canceling system.
use the NLMS (Normalized Least Mean Square) algorithm. The ANC system composed of two separate inputs, a
The correction applied to the weight vector w(n) at iteration primary input or ECG signal source which is shown as s(n)
n+1 is “normalized” with respect to the squared Euclidian and a reference input that is the noise input shown as x(n) .
norm of the input vector x(n) at iteration n. The primary signal is corrupted by noise x1(n). The signal
We may view the NLMS algorithm as a time-varying step- x1(n) is highly correlated with noise signal or reference signal
size algorithm, calculating the convergence factor µ as in Eq. x(n). Desired signal d(n) results from addition of primary
(3)[5]. signal s(n) and correlated noise signal x1(n). The reference
signal x(n) is fed into adaptive filter and its output y(n) is
subtracted from desired signal d(n). Output of the summer
block is then fed back to adaptive filter to update filter
Where: α is the NLMS adaption constant, which optimize
coefficients. This process is run recursively to obtain the
the convergence rate of the algorithm and should satisfy the noise free signal which is supposed to be the same or very
condition 0< α<2, and c is the constant term for normalization similar to primary signal s(n) .
and is always less than 1.
In NLMS algorithm, the filter weights are updated by the
Eq. (4).

C. RECURSIVE LEAST SQUARE A LGORITHM
The RLS algorithm is known for its excellent performance
when working in time varying environments but at the cost
of an increased computational complexity and some stability
problems. In this algorithm the filter tap weight vector is
updated using Eq. (5) [7]. Figure 2. Adaptive Noise Cancellation system

IV. SIMULATION RESULTS
Eq. (6) and (7) are intermediate gain vector used to compute
The adaptive noise canceller was implemented in
tap weights.
MATLAB for three algorithms; LMS, NLMS and RLS [7]. In
the simulation the reference input signal x(n) was a white
Gaussian noise of power 1-dB generated using randn function
Where: λ is a small positive constant very close to, but in MATLAB, and source signal s(n) was a clean amplified
smaller than 1. ECG signal recorded with 12-lead configuration [6], the desired
The filter output is calculated using the filter tap weights signal d(n) ,obtained by adding a delayed version of x(n)
of previous iteration and the current input vector as in Eq. into clean signal s(n), d(n) = s(n) + x1(n) as shown in Fig.3.
(8).
© 2012 ACEEE 2
DOI: 01.IJSIP.03.01. 47


The simulation of the LMS, NLMS and RLS algorithms was algorithms.When input signal is non-stationary in nature,
carried out with the following specifications: the RLS algorithm proved to have the highest convergence
Filter order N=19, step size µ= 0.009 and iterations= 1000, speed, less MSE,and higher SNR Improvement but at the
c= 0.001 cost of large computational complexity and memory
The LMS filtered output is shown in Fig.4 (a), the mean requirement. The NLMS algorithm changes the step-size
squared error generated as per adaption of filter parameters according to the energy of input signals hence it is suitable
is shown in Fig.4 (b). The step size µ control the performance for both stationary as well as non-stationary environment
of the algorithm, if µ is too large the convergence speed is and its performance lies between LMS and RLS. Hence it
fast but filtering is not proper, if µ is too small the filter gives provides a trade-off in convergence speed and computational
slow response, hence the selection of proper value of step- complexity. The implementation of algorithms was
size for specific application is prominent to get good results. successfully achieved, with results that have a really good
Fig.5 and Fig.6 shows the output results for NLMS and RLS response.
algorithms respectively. If we investigate the filtered output
of all algorithms, LMS adopt the approximate correct output ACKNOWLEDGMENT
in 750 samples, NLMS adopt in 600 samples and RLS adopt
The authors gratefully acknowledge Dr. B.D. Gupta
in 250 samples. This shows that RLS has fast learning rate.
Director, Anand Engineering College, Agra, India and the
In Table-I performance analysis of all three algorithms is
department of Electronics and Communication, Bundelkhand
presented in term of Mean Squared Error (MSE),
Institute of Engineering and Technology (BIET), Jhansi, U.P.,
computational complexity and stability. It is clear from the
India for providing necessary support and research facilities.
Table-I, the computational complexity and stability problems
increases in an algorithm as we try to reduce the mean squared
REFERENCES
error.
[1] Bernard Widrow, John R. Glover, John M. Mccool, John
TABLE I. PERFORMANCE COMPARISON OF ADAPTIVE ALGORITHMS
Kaunitz, Charles S. Williams, Robert H. Hean, James R. Zeidler,
Eugene Dong, Jr. and Robert C. Goodlin, “Adaptive Noise
Cancelling: Principles and Applications”, Proceedings of the IEEE,
1975, Vol.63 , No. 12 , Page(s): 1692 – 1716.
[2] J. Benesty, F. Amand , A. Gilloire and Y. Grenier , “Adaptive
Filtering Algorithms for Stereophonic Acoustic Echo Cancellation”,
In Table-II SNR Improvement is presented for each International Conference on Acoustics, Speech, and Signal
algorithm. From Table-I & Table-II it is clear that the RLS Processing, 1995, vol.5, Page(s): 3099 – 3102.
[3] Simon Haykin, “Adaptive Filter Theory”, Prentice Hall, 4 th
algorithm has best performance but same time the
edition.
computational complexity is also increased. If we investigate [4] Paulo S.R. Diniz, “Adaptive Filtering: Algorithms and Practical
NLMS algorithm its performance is comparable with RLS Implemetations” ,ISBN 978-0-387-31274-3, Kluwer Academic
algorithm with slight additional complexity hence NLMS is Publisher © 2008 Springer Science+Business Media, LLC, pp.77-
the favorable choice for most of the industries. 195.
[5] Abhishek Tandon, M. Omair Ahmad, “An efficient, low-
TABLEII. COMPARISON OF SNR IMPROVEMENT
complexity, normalized LMS algorithm for echo cancellation” The
2nd Annual IEEE Northeast Workshop on Circuits and Systems,
2004. NEWCAS 2004, Page(s): 161 – 164.
[6] Ch. Renumadhavi, Dr. S.Madhava Kumar, Dr. A. G. Ananth,
Nirupama Srinivasan, “A New Approach for Evaluating SNR of
ECG Signals and Its Implementation”, Proceedings of the 6th
CONCLUSIONS WSEAS International Conference on Simulation, Modelling and
Optimization, Lisbon, Portugal, September 22-24, 2006.
The main objective of this paper was to implement an [7] Ying He Hong He, Yi Wu and Hongyan Pan, “The Applications
adaptive noise canceller for de-noising an ECG signal and and Simulation of Adaptive Filter in Noise Canceling”, 2008
test the performance of the system for various adaptive International Conference on Computer Science and Software
Engineering, 2008, Vol.4, Page(s): 1 – 4.

© 2012 ACEEE 3
DOI: 01.IJSIP.03.01. 47


Figure. 3(a) Clean ECG signal s(n);(b) Noise signal x(n);(c) desired
signal d(n)

Figure 4. MATLAB simulation for LMS algorithm; N=19, step
size=0.009

Figure 5. MATLAB simulation for NLMS algorithm; N=19, step
size=0.001

Figure 6. MATLAB simulation for RLS algorithm; N=19, λ=1

© 2012 ACEEE 4
DOI: 01.IJSIP.03.01. 47

Simulation of Adaptive Noise Canceller for an ECG signal Analysis

More Related Content

What's hot (19)

Viewers also liked (16)

Similar to Simulation of Adaptive Noise Canceller for an ECG signal Analysis (20)

More from IDES Editor (20)

Recently uploaded (20)

Simulation of Adaptive Noise Canceller for an ECG signal Analysis