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

发布时间 : 星期六 文章电子科技大学信号处理实验4 - 数字滤波器设计及实现更新完毕开始阅读9ac8e69e5122aaea998fcc22bcd126fff7055deb

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、 在设计滤波器时,不给出通带,阻带等参数,由学生自行根据需要设计各参数,在 设计过程中,加深对滤波器的理解。

报告评分:

指导教师签字: