Adaptive Noise Canceller

Adaptive Noise Canceller
Project Report
In this report, an adaptive noise canceller will be presented

and some useful observations will be done with the canceling
process over the audio signal.
Group Members :
 Ahmad Naeem
 Mobeen Sohail
 Jibran Shahid
mee
30/12/2016
Noise canceling in audio signal with adaptive filter
Introduction :
All the physical systems when they operate, they produce physical noise, which is a mixing of an
infinitive number of sinusoidal harmonics with different frequencies. So, the initial signal
information is corrupted with this noise. This complex signal, may be very noisy, so much that
the human ear or other system which may follow it, cannot receive correct the initial signal.
So, an algorithm has to be invented which must be able to separate the physical noise from
the information signal and to output the information signal without noise. This is not possible,
as there is not a perfect system. So, this algorithm should have the ability to reduce the noise
level as much as it can.
A good noise-canceling algorithm is the algorithm, which is presented in this paper. This
algorithm is based on the Least Mean Square algorithm. This system, operates like an automatic
control system with feedback where the feedback performs an adaptation to the filter
coefficients, otherwise it adjust the filtering of the noise in each sample.
Problem formulation :
A simple audio signal is represented by the following mathematical equation
1 x(n) sin(2* pi* fo/ fs*n)
Where f 0 is the frequency of the audio signal, fs is the sampling frequency and nis the
discrete time base.

Fig.1 The simple audio signal without noise
When white noise is added to this signal then it becomes
2 c n sin 2* pi* fo/ fs*n noise n

So, the problem is to extract the useful signal from this complex signal.
Fig.2 The audio signal with noise
A simple adaptive filter noise canceller is presented in the following diagram.
0
Fig.3 Adaptive filter noise canceller
Brief adaptive filter theory :

One of the most successful adaptive algorithms is the LMS algorithm developed by Widrow and
his coworkers. Instead of computing Wopt which is the optimal solution of the wiener filter, the
LMS coefficients are adjusted from sample to sample in such a way as to minimize the MSE
(mean square error).
The LMS is based on the steepest descent algorithm where the weight vector is updated from
sample to sample as in the next equation.
3 Wk 1 Wk k
Where Wk and k are the weight or the impulse response w (n) of the filter and the true
gradient vectors respectively, at the kth sampling instant. Controls the stability and rate of
convergence.
The steepest descent algorithm in the previous equation still requires knowledge of R and P ,
since k is obtained by the next equation.
dJ
4 k 2P 2RW 0dW
The LMS algorithm is a practical method of obtaining estimations of the filter weights W k in real
time without the matrix inversion of the next equation Wopt R 1P or the direct computation of
the autocorrelation and cross correlation. The Widrow – Hopf LMS algorithm for updating the
weights from sample to sample is given by
5 Wk 1 Wk 2 ek X k
Where: ek yk WkTXk
Clearly, the LMS algorithm above does not require prior knowledge of the signal statistics (that
is the correlations R and P ), but instead uses their instantaneous estimations. The weights
obtained by the LMS algorithm are only estimations, but these estimations improve gradually
with time as the weights are adjusted and the filter learns the characteristics of the signal.
Eventually, the weights converge. The condition of convergence is:
6 0 1 max
Where max is the maximum eigenvalue of the input data covariance matrix In practice, Wk
never reaches the theoretical optimum (the Wiener solution), as the algorithm iteration step
can never be zero but it has a discrete value.
Furthermore, the LMS algorithm can be summarized in the next equation for the
implementation of the filter. N 1
7 y k wi k *u k i
i 0
8 h( 1) h( ) u( )e( )
In which: e( ) d( ) y k
In the above equation, T denotes the transpose of the vector hinv. The desired response d( )
can be found from the adaptive filter scheme to be
9 d( ) hxT( )w( g)
The notation u( ) denotes a vector of samples arranged as
T
10 u(k) u( ) u( 1)...............u( N 1)
Where N is the order of the filter.
Problem solution :
The problem solution may be based on standard adaptive techniques and especially on the
adaptive filter noise canceller.
This filter gets the physical noise and filters it in such a way that the algebraic addition with the
complex signal produces a new signal with a very low level of noise. The schematic diagram of
the adaptive filter noise canceller was given above and it is presented again in the following
diagram.
Fig.4 The used adaptive noise canceller
As it can be observed, the adaptive filter noise canceller is a system with feedback as a system
with automatic control. For this reason, there is stability problem. The problem arises when the
error signal’s amplitude increases in accordance with time or when no convergence is obtained.
The mathematical algorithm, which is used by the adaptive filter for the filter’s coefficients
update, is called LMS (Least mean square) algorithm. This algorithm tries to minimize the
square of the error signal so; the audio signal has been reconstructed quite well.
The system output is the error signal. The error signal is the algebraic subtraction between the
audio signal with noise and the noise. When the LMS algorithm performs a good convergence
then the only error signal is the audio signal. So, the error signal must be the noise canceller
output.
The LMS algorithm has also an initial step variable μ, which is the step of the convergence in
each iteration. In the matlab code where the adaptive filter noise canceller is built, the NLMS
algorithm used. The only difference is that the step variable μ, changes in each iteration with
the next equation:
11 ( ) 2
u k
where is a positive constant, usually less than 2 , and is a small positive constant. The
symbol u( ) 2 is the energy of the signal in a predefined buffer.
Implementation of the adaptive filter noise canceller in matlab and waveforms:

The adaptive noise canceller was implemented in matlab as it offers the opportunity for a
quick and very professional way to study the whole operation of the system.
The matlab code, which was written, is presented now.
function [x] = noice_canceller ;

clear all;
%********************** %****
Time Base *******
%**********************
N = 2000; n = 1:1:N;
%*********************** %****
System Vars
%***********************
fo = 500;
Fs = 44100;
db = 5;
%*********************** %*****
Signal Generation
%***********************
audio = sin(2*pi*fo/Fs*n); %+ 0.7*sin(2*pi*5*fo/Fs*n));
%**********************
%**** Noise Generation
%*********************
%noise_plus_audio = awgn(audio,db);
noise = wgn(1,N,db)/2;
%noise = 0.1*randn(1,N);
for i=1:1:N
noise_plus_audio(i) = (audio(i) + noise(i)/2);
end
%**************************************
%**** Initialization of Noise Canceller
%**************************************
w0 = zeros(1,32);
% Initial filter coefficients
%wi = zeros(1,255);
mu = 0.4;
% LMS step size
S = initnlms(w0,mu,[],1);
%*************************************** %*****
Apply Noise Canceller over signal
%***************************************
fori=1:1:N
[y,canceller_output(i),S] = adaptnlms(noise(i),noise_plus_audio(i),S);
%***************************************
%***** Calculate FFT of Signals
%***************************************
ff_mix = abs(fft(noise_plus_audio,1024));
ff_output = abs(fft(canceller_output,1024));
%***************************************
%***** Calculate Error
%***************************************
error(i) = audio(i) - canceller_output(i);
%*********************** %***
Plots *************
%***********************
subplot(6,1,1);
plot(noise_plus_audio,'r');
axis([0 N -2 2]);
title('Signal and Noise');
subplot(6,1,2);
plot(audio,'g');
axis([0 N -2 2]);
title('Audio Signal');
subplot(6,1,3);
plot(canceller_output,'b');
title('Canceller Output');
axis([0 N -2 2]);
subplot(6,1,4);
plot(ff_mix,'r');
title('FFT of MixedSignal');
subplot(6,1,5);
plot(ff_output,'r');
title(' FFT of Output');
subplot(6,1,6);
plot(error,'g');
axis([0 N -2 2]);
title('Error of FFT');
drawnow;
end
The output waveforms of the system are plotted in the next figure.
Conclusion :
To sum up, the adaptive noise canceller is a very efficient and useful system in many
applications with sound, video etc. The only disadvantage is that it needs digital signal
processor DSP for its operation
References:
 Simon Haykin, Adaptive Filter Theory, Prentice Hall 2002
 http://matlabcode.weebly.com/matlab-projects/noise-cancellation-from-audio-signalby-
using-adapter-filter
 Mathworks.com

Adaptive Noise Canceller

Uploaded by

Copyright:

Available Formats

Adaptive Noise Canceller

Uploaded by

Document Information

Original Description:

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Adaptive Noise Canceller

Uploaded by

Copyright:

Available Formats

Adaptive Noise Canceller

In this report, an adaptive noise canceller will be presented

1 x(n) sin(2* pi* fo/ fs*n)

discrete time base.

When white noise is added to this signal then it becomes

2 c n sin 2* pi* fo/ fs*n noise n

Fig.2 The audio signal with noise

A simple adaptive filter noise canceller is presented in the following diagram.

Brief adaptive filter theory :

The notation u( ) denotes a vector of samples arranged as

Fig.4 The used adaptive noise canceller

Implementation of the adaptive filter noise canceller in matlab and waveforms:

function [x] = noice_canceller ;

You might also like