widebandcancelling.m

Wide-band noise cancellation

clear all;
clc;
% Wide-band noise cancellation
Fss=8000;
mu = 0.2;  % adaptation parameter
M = 11;     % filter length
delay = -10;    % delay may be + or - 
N = 40128;
load gong;y=y';s=y; % signal
noise = 0.2*randn(1,N);             % noise
y = s; %+ noise;
h =zeros(M,1);     % initial filter coefficients
for n = M+1+delay:N-M,  % adaptive algorithm
   x = y(n-delay-M+1:n-delay)';
   q(n) = h'*x;
   e(n) = y(n) - q(n);
   h = h + mu * e(n) * x;
end
hold off
plot(y,'g')
hold on
plot(q,'r')
hold on;
plot(e)
p8K = audioplayer(y/max(abs(y)),Fss);
playblocking(p8K);pause(1)
p8K = audioplayer(q/max(abs(q)),Fss);
playblocking(p8K);pause(1)
p8K = audioplayer(e/max(abs(e)),Fss);
playblocking(p8K);