www.pudn.com > speechenhancement.rar > whpro.m
clear;
close all;
load('mtlb','mtlb');
M=length(mtlb);
s=mtlb';
>M=input('input the length of signal>>>>>M=');
n=1:M;
>s=exp(-0.002*n).*sin(pi*n/50);
w=0.5*randn(1,M);>try 0.1 0.4,1
x=s+w;
Rss=xcorr(s,s);
Rww=xcorr(w,w);
[h,e]=WH(Rss,Rww,M)
ss=filter(h,2,x);
subplot 221
plot(n,s);
title('Source signal');
subplot 222
plot(n,w);
title('noise');
subplot 223
plot(n,x);
title('Observated signal');
subplot 224
plot(n,ss);
title('Filtered signal');
figure;
plot(n,ss-s);
title('Error between filtered signal and Source signal')