Tuner Upper Removal Demonstration: Difference between revisions

From Class Wiki
Jump to navigation Jump to search
← Older editNewer edit →
No edit summary
No edit summary
Line 1: Line 1:
  % Demonstration of LMS algorithm for noise cancellation.
  % Demonstration of LMS algorithm for noise cancellation.  Tuner upper problem.
% Rich Kozick, Spring 1997
% Rich Kozick, Spring 1997
% Rob Frohne's  modifications for Macintosh 2000.  
% Rob Frohne's  modifications for Macintosh 2000
% Desired signal
% and Linux 2006 or so.
clear all
 
Totaltime=1;
 
  speak('Hit a key and speak the signal.');  
clear all
pause;
Fs=8000; % 8 Khz sampling for Linux.
[st, Fs] = recordsound(Totaltime, 22050, 1);
T0 = 2; % 2 seconds
s = st';
 
Ls = length(s);
system("espeak 'After hitting enter in the command window, speak the signal.'");
% Interference + random noise
st = record(T0,Fs);
speak('Hit a key and make the noise!');  
%st=wavread('Hello.wav');
pause;
st=st(:,1);
[nt,Fs] = recordsound(Totaltime, 22050, 1);
 
n = nt';
system("espeak 'Here is the signal.'")
%Sign = 0.01;
soundsc(st);
  %Dn=20;  % Delay of the noise that appears in y.
 
%n = n(1:Ls) + Sign*randn(Ls,1);
Ls = length(st);
%an = [0 .01 -.5 1 -.5 .1 .01 0];
T0=length(st)/Fs;
an = 4*[0 0 0 0 0 0 0 0 0 0 0 0 0 0 .5 1 .5];
%  Make the tuner upper noise or even several at the same time.
bn = [1];
t=0:1/Fs:T0-1/Fs;
%sys = tf(an,bn,1/Fs);
n =  10*(sin(2*pi*100*pi*t) + cos(2*pi*600*t) + 1.5*sin(2*pi*850*t));
%bode(sys);
 
%figure(1);
% Add him or them to the desired signal.
nf = filter(an,bn,n);
x = st' + n;
  %nf = [nf(Dn:(Ls)); zeros(Dn-1,1)] + Sign*randn(Ls,1); 
 
y = s + nf;  
system("espeak 'Here is the noisy signal.'")
Sigx = 0.01;
soundsc(x,Fs);
bx = [1];      % bx and ax are filtering on n to produce x
 
ax = [1];
N = 64;        % Length of adaptive filter
Dx = 1;        % Delay of n that appears in x
 
x = filter(bx, ax, n);%x = [x(Dx:(Ls)); zeros(Dx-1,1)] + Sigx*randn(Ls,1);
% LMS algorithm for adaptive noise cancellation
%x = n + Sigx*randn(Ls,1); 
 
speak('Here is the noisy signal.')
h = ones(N,1);
soundsc(y,Fs);
mu = 0.1;
N = 20;        % Length of adaptive filter  
for k=N:Ls
% LMS algorithm for adaptive noise cancellation  
  xk = x(k:-1:(k-N+1));
h = zeros(N,1);
  y(k) = h'*xk';
mu = 1/(10*N*var(x));
  e(k) = y(k);
%mu = 1.0;
  h = h - 2*mu*e(k)*xk'/(xk*xk');
for k=N:Ls    
end
    xk = x(k:-1:(k-N+1));    
 
  nhat(k) = h'*xk;    
% The signal estimate is in the vector e
  e(k) = - y(k) + nhat(k);    
system("espeak 'Here is a scaled version of the tail of the cleaned signal.'");
  h = h - mu*e(k)*xk;%/(xk'*xk);
skip =1000;
end  
soundsc(e(skip:length(e)),Fs);
% The signal estimate is in the vector e
speak('Here is the cleaned signal.');
soundsc(e,Fs);

Revision as of 14:43, 6 December 2015

% Demonstration of LMS algorithm for noise cancellation.  Tuner upper problem.

% Rich Kozick, Spring 1997 % Rob Frohne's modifications for Macintosh 2000 % and Linux 2006 or so.


clear all Fs=8000; % 8 Khz sampling for Linux. T0 = 2;  % 2 seconds

system("espeak 'After hitting enter in the command window, speak the signal.'"); st = record(T0,Fs); %st=wavread('Hello.wav'); st=st(:,1);

system("espeak 'Here is the signal.'") soundsc(st);

Ls = length(st); T0=length(st)/Fs; % Make the tuner upper noise or even several at the same time. t=0:1/Fs:T0-1/Fs; n = 10*(sin(2*pi*100*pi*t) + cos(2*pi*600*t) + 1.5*sin(2*pi*850*t));

% Add him or them to the desired signal. x = st' + n;

system("espeak 'Here is the noisy signal.'") soundsc(x,Fs);

N = 64;  % Length of adaptive filter

% LMS algorithm for adaptive noise cancellation

h = ones(N,1); mu = 0.1; for k=N:Ls 

 xk = x(k:-1:(k-N+1));
 y(k) = h'*xk';
 e(k) = y(k);
 h = h - 2*mu*e(k)*xk'/(xk*xk');

end

% The signal estimate is in the vector e system("espeak 'Here is a scaled version of the tail of the cleaned signal.'"); skip =1000; soundsc(e(skip:length(e)),Fs);

Retrieved from "https://fweb.wallawalla.edu/class-wiki/index.php?title=Tuner_Upper_Removal_Demonstration&oldid=11521"