zfir1.m

用MATLAB实现数字信号处理的程序

%Fir1函数：设计具有标准频率响应的FIR滤波器
Wn=[0.4 0.7];
n=38;
b1=fir1(n,Wn,'stop');               %使用默认的海明窗时设计的带阻滤波器
window=chebwin(n+1,30);
b2=fir1(n,Wn,'stop',window);        %使用切比雪夫窗时设计的带阻滤波器
[H1,W1]=freqz(b1,1,512,2);
[H2,W2]=freqz(b2,1,512,2);
subplot(2,1,1),plot(W1,20*log10(abs(H1)));
xlabel('归一化频率'),ylabel('幅度/dB ');
grid;
subplot(2,1,2),plot(W2,20*log10(abs(H2)));
xlabel('归一化频率'),ylabel('幅度/dB');
grid;