yyqz.txt

用汉明窗作的语音去噪低通滤波器

x1=wavread('D:\Program Files\MATLAB71\work\原信号.wav');
t=0:1/(8*1024):86496/(8*1024);%86497B
A=0.5;
d=[A*sin(2*pi*4000*t)]';
x2=x1+d;
wavwrite(x2,8192,8,'D:\Program Files\MATLAB71\work\加噪后信号.wav');
figure(1)
subplot(3,1,1);
plot(t,x1);
title('原信号时域波形');
subplot(3,1,2);
plot(t,x2);
title('加噪后时域波形');
fs=8*1024;
f=fs*(0:1023)/1024;
y1=fft(x1,1024);
y2=fft(x2,1024);
%%%%%以下为滤波器
wp=2*pi*3.0/8
ws=2*pi*3.5/8
width=ws-wp
N=ceil(6.2*pi/width)
wn=(wp+ws)/pi/2
b=fir1(N,wn,hanning(N+1));
figure(2)
freqz(b,1)
%%%%%%滤波
x3=filter(b,1,x2);
y3=fft(x3,1024);
figure(3)
subplot(3,1,1);
plot(f,abs(y1));
%axis([0 9000 0 320]);
title('原信号频域波形');
subplot(3,1,2);
plot(f,abs(y2));
axis([0 9000 0 100]);
title('加噪后频域波形');
subplot(3,1,3);
plot(f,abs(y3));
title('滤波后频域波形');
wavwrite(x3,8192,8,'D:\Program Files\MATLAB71\work\滤波后信号.wav');
figure(1)
subplot(3,1,3);
plot(t,x3);
title('滤波后信号时域波形');