fir_b.m

数字信号课程设计

x1=wavread('音乐.wav');
fs=22050;
wp1=2*pi*1200/fs;wp2=2*pi*3000/fs;
ws1=2*pi*1000/fs;ws2=2*pi*3200/fs;
wp=(wp1+ws1)/2;
ws=(wp2+ws2)/2;
wdelta=wp1-ws1;

N=ceil(8*pi/wdelta);       %取整
wn=[wp ws];

beta=3.39;
window=kaiser(N+1,beta); 
[b,a]=fir1(N,wn/pi,'high',window); 

figure(1)
freqz(b,a,512);
title('FIR带通滤波器');

f2=filter(b,a,x1);

figure(2)
subplot(2,1,1)
plot(x1)
title('FIR带通滤波器滤波前的时域波形');
subplot(2,1,2)
plot(f2);
title('FIR带通滤波器滤波后的时域波形');
F0=fft(f2,1024);
f=fs*(0:511)/1024;

figure(3)
y2=fft(x1,1024);
subplot(2,1,1);
plot(f,abs(y2(1:512)));
title('FIR带通滤波器滤波前的频谱')
xlabel('频率/Hz');
ylabel('幅值');
subplot(2,1,2)
F2=plot(f,abs(F0(1:512)));
title('FIR带通滤波器滤波后的频谱')
xlabel('频率/Hz');
ylabel('幅值'); 
sound(f2,44100);    %滤波后的语音信号