电子科技大学信号处理实验4 - 数字滤波器设计及实现

xlabel('频率/Hz');ylabel('幅值/A');

Ffs=20000; Wp=2200*2/Ffs; Ws=3500*2/Ffs; Rp=1;

Rs=20;%最小阻带衰减改为>20dB [N,Wn]=buttord(Wp,Ws,Rp,Rs); [num,den]=butter(N,Wn);

y=filter(num,den,x); subplot(2,2,3); plot(t,y);

title('IIR滤波器滤波后的时域图'); xlabel('时间/s');ylabel('y'); r_f=fft(y)/n; r_f=fftshift(r_f); subplot(2,2,4); plot(f,abs(r_f));

title('IIR滤波器滤波后的频域图'); xlabel('频率/Hz');ylabel('幅值/A');

（6） 利用“fdatool”设计IIR级联和并联实现得到的滤波器系数，滤波代码

和结果图。（拓展要求）

九、实验结论：说明：红色字在写报告时去掉！

1、 该实验中，实现同样功能的IIR滤波器需要6阶，FIR滤波器需要19阶，可以明显的看出，实现同样的功能FIR需要更多的阶数。

2、由FIR、IIR滤波器的图可以看出，在通带范围内，FIR具有良好的线性 相位。

3、降低截止频率后，由于旁瓣数量减小，导致滤波后的旁瓣幅值下降。

十、总结及心得体会：

1、 butter既可以用于模拟滤波器的设计，也可以用于数字滤波器的设计。用 于模拟滤波器的设计时，在函数的参数中加‘s’。

2、 freqz函数用于，当已知频响函数的系数时画图。当已知单位冲击响应时， 可以通过在时域的值后加0，再DFT变换，求得其频域响应。

十一、对本实验过程及方法、手段的改进建议：

1、 增加带通，带阻，高通滤波器的设计，进一步比较FIR与IIR滤波器的差异。

2、 在设计滤波器时，不给出通带，阻带等参数，由学生自行根据需要设计各参数，在 设计过程中，加深对滤波器的理解。

报告评分：

指导教师签字：