COURSE: Advanced Digital Signal Processing (MET 1113)
LAB EXPERIMENT 1 (CLO 2, PLO 1, PLO 2, PLO 6, P7 ): STUDY OF ADAPTIVE
NOISE CANCELLATION USING RLS ADAPTIVE FILTERING
Rubric for the lab report
Excellent Very Good Good Poor Very Poor
Criteria 5 4 3 2 1
Analytical Skills: Fully able to Reasonably Not quite Poor ability to Very poor
Ability to use an use an able to use able to useuse an ability to use
recursive least recursive an recursive an recursive recursive least an recursive
square (RLS) filter least square least square least square square (RLS) least square
to extract useful (RLS) filter (RLS) filter (RLS) filter filter to (RLS) filter
information from a to extract to extract to extract extract useful to extract
noisy signal useful useful useful information useful
information information information from a noisy information
from a noisy from a noisy from a noisy signal from a noisy
signal signal signal signal
Proficiency in Using Fully able to Reasonably Not quite Poor ability to Very poor
Tools: run the code able to run able to run run the code ability to run
Ability to run the program the code the code program using the code
code program using using the program program the MATLAB program
the MATLAB MATLAB using the using the software. using the
software. software. MATLAB MATLAB MATLAB
software. software. software.
Observation: Fully able to Reasonably Not quite Poor ability to Very poor
Ability to show the show the able to show able to show show the ability to
output waveform. output the output the output output show the
waveform. waveform. waveform. waveform. output
waveform.
FACULTY OF ENGINEERING AND BUILT
ENVIRONMENT
MASTER OF SCIENCE IN ELECTRICALS AND
ELECTRONICS ENGINEERING
ADVANCED DIGITAL SIGNAL PROCESSING
MET 1113
LABORATORY MANUAL
STUDENT NAME :
STUDENT ID :
FACULTY OF ENGINEERING AND BUILT
ENVIRONMENT
EXPERIMENT 1 (CLO 2, PLO 1, PLO 2, PLO 6, P7):
STUDY OF ADAPTIVE NOISE CANCELLATION USING RLS ADAPTIVE
FILTERING
1.0 Objectives
1. To use an recursive least square (RLS) filter to extract useful information from a noisy
signal.
2. To show the output waveform using the MATLAB software.
2.0 Equipment
Software: MATLAB
3.0 Theory
The information bearing signal is a sine wave that is corrupted by additive white gaussian
noise. The adaptive noise cancellation system assumes the use of two microphones. A primary
microphone picks up the noisy input signal, while a secondary microphone receives noise that is
uncorrelated to the information bearing signal, but is correlated to the noise picked up by the
primary microphone.
The model illustrates in Figure 1 shows the ability of the Adaptive RLS filter to extract
useful information from a noisy signal.
Figure 1: Adaptive Noise Cancellation Demo
4.0 Procedure
1. The information bearing signal is a sine wave of 0.055 cycles/sample
2. Run the code program below:
>> signal = sin(2*pi*0.055*(0:1000-1)');
>> signalSource = dsp.SignalSource(signal,'SamplesPerFrame',100,...
'SignalEndAction','Cyclic repetition');
>> plot(0:199,signal(1:200));
>> grid; axis([0 200 -2 2]);
>> title('The information bearing signal')
Display the output results.
3. The noise picked up by the secondary microphone is the input for the RLS adaptive filter. The
noise that corrupts the sine wave is a lowpass filtered version of (correlated to) this noise. The
sum of the filtered noise and the information bearing signal is the desired signal for the adaptive
filter.
4. Run the code program below.
>> nvar = 1.0; % Noise variance
>> noise = randn(1000,1)*nvar; % White noise
>> noiseSource = dsp.SignalSource(noise,'SamplesPerFrame',100,...
'SignalEndAction','Cyclic repetition');
>> plot(0:999,noise);
>> title('Noise picked up by the secondary microphone');
>> grid; axis([0 1000 -4 4]);
Display the output results.
5. The noise corrupting the information bearing signal is a filtered version of 'noise'. Initialize the
filter that operates on the noise.
6. Run the code program below:
lp = dsp.FIRFilter('Numerator',fir1(31,0.5));% Low pass FIR filter
7. Set and initialize RLS adaptive filter parameters and values:
M = 32; % Filter order
delta = 0.1; % Initial input covariance estimate
P0 = (1/delta)*eye(M,M); % Initial setting for the P matrix
rlsfilt = dsp.RLSFilter(M,'InitialInverseCovariance',P0);
8. Running the RLS adaptive filter for 1000 iterations. As the adaptive filter converges, the
filtered noise should be completely subtracted from the "signal + noise". Also the error, 'e',
should contain only the original signal.
9. Run the code program below:
>> scope = timescope('TimeSpan',1000,'YLimits',[-2,2], ...
'TimeSpanOverrunAction','Scroll');
for k =1:10
n = noiseSource(); % Noise
s = signalSource();
d = lp(n) + s;
[y,e] = rlsfilt(n,d);
scope([s,e]);end
>> release(scope);
Display the output results.
10. Show the convergence plot of the adaptive filter response to the response of the FIR filter by
run the code program below.
>> H = abs(freqz(rlsfilt.Coefficients,1,64));
>> H1 = abs(freqz(lp.Numerator,1,64));
>> wf = linspace(0,1,64);
>> plot(wf,H,wf,H1);
>> xlabel('Normalized Frequency (\times\pi rad/sample)');
>> ylabel('Magnitude');
>> legend('Adaptive Filter Response','Required Filter Response');
>> grid;
>> axis([0 1 0 2]);
Display the output results.
5.0 Conclusion
__________________________________________________________________
__________________________________________________________________
__________________________________________________________________
__________________________________________________________________
