fir_h.m

来自「数字信号课程设计」· M 代码 · 共 41 行

fs=22050;
x1=wavread('音乐.wav');
wp=2*pi*5000/fs;ws=2*pi*4800/fs;

wdelta=wp-ws;
N=ceil(8*pi/wdelta);       %取整
wn=(wp+ws)/2;

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);    %滤波后的语音信号