利用该案例理论见前文
Ahlers:FIR滤波器窗函数设计法案例2之理论准备zhuanlan.zhihu.comMatlab实现效果:
可以看出阻带的频率衰减是大于40dB,通带的最大起伏只有0.05dB,小于0.2dB, 满足题意。
Maltab实现代码:
Fs=1200;
M=48;
n=0:M;
om=linspace(0,pi,200); %for displaying frequency response
wc=0.415*pi; %cut-off frequency
h=sinc(n-M/2)-wc/pi*sinc(wc/pi*(n-M/2));
hn=h.*(hann(M+1))';
H=freqz(hn,1,om);
subplot(2,1,1);
plot(om*1200/(2*pi),20*log10(abs(H)));
xlabel('Frequency[Hz]');
ylabel('|H(w)|[dB]');
title('Magnitude res