IIR数字滤波器设计及应用.docxVIP

本科学生验证性实验报告 学号 114090315 姓名 李开斌 学院 物电学院 专业、班级 11 电子 实验课程名称 数字信号处理 教师及职称 李宏宁 开课学期 2013 至 2014 学年 下 学期 填报时间 2014 年 5 月 14 日 云南师范大学教务处编印 bj z j bj z j j0 N H(z) i ai z * i 实验序号 09 实验名称 IIR 数字滤波器设计及应用 实验时间 2014 年 5 月 14 日 实验室 同析 3 栋 313 1. 实验目的 2. 加深理解 IIR 数字滤波器的特性， 掌握 IIR 数字滤波器的设计原理与设计方法， 以及 IIR 数字滤波器的应用。 调用 butter 函数设计出巴特沃斯滤波器， 格式为 [b,a] = butter(N,Wc,options) 输入参数： N 和 Wc 是 buttord 函数返回的参数，含义见上。 Options= low , high , bandp分a别ss对应, 低通s、top高通,、 带通、带阻，默认情况下为低通或带通。 输出参数： b 和 a 为设计出的 IIR 数字滤波器 H(s)的分子多项式和分母多项式的 系数矩阵。 ( 1) Chebyshev I型滤波器设计 Chebyshev I 型滤波器为通带纹波控制器：在通带呈现纹波特性，在阻带单调衰 减。 [N,Wc] = cheb1ord(Wp, Ws, Ap, As) [b,a] = cheby1(N,Ap,Wc,options) 参数含义与 butter 中参数一致。 2.(2) Chebyshev II 型滤波器设计 Chebyshev II 型滤波器为阻带纹波控制器： 在阻带呈现纹波特性。 [N,Wc] = cheb2ord(Wp, Ws, Ap, As) [b,a] = cheby2(N,As,Wc,options) 椭圆滤波器设计 椭圆滤波器在通阻带都呈现纹波特性。 [N,Wc] = ellipord(Wp,Ws,Ap,As) [b,a] = ellip(N,Ap,As,Wc,options) 利用 Matlab 中的函数分析并绘出常用基本信号的波形。 注意事项 ： （ 1）在使用 MATLAB 时应注意中英输入法的切换， 在中文输入法输入程序时得 到的程序是错误的； （2）MATLAB 中两个信号相乘表示为 x.*u,中间有个‘ .，同样两个信号相除也 是如此； （ 3）使用 MATLAB 编写程序时，应新建一个 m 文件，而不是直接在 Comandante 窗口下编写程序； 在使用ＭＡＴＬＡＢ编程时，应该养成良好的编写习惯。 5．实验数据处理方法 比较法 画图法 1．实验现象与结果 2 1.信号 x[k] 1 cos( k) cos(2 k) ,确定设计指标，实现各种 IIR 数字滤波器 43 以实现以下信号处理。 2 ( 1)设计 IIR 低通滤波器，滤除 cos(2 k) 的成分。 3 % 低通滤波器 [N,Wc]=buttord(0.250,0.677,3,60) [b,a]=butter(N,Wc,low) freqz(b,a); axis([0,1,-120,0]); grid on title( 巴特沃斯低通数字滤波器 ) 巴特沃斯低通数字滤波器 0 )(edB eud -50 itdu n nga M -100 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Normalized Frequency ( rad/sample) )seerged(esah )seerged(esah 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1Normalized Frequency ( rad/sample) -200-400 -200 -400 -200-400 -200 -400 ( 2)设计 IIR 高通滤波器，滤除 1 cos( k) 的成分 4 % 高通滤波器 [N,Wc]=buttord(0.250,0.677,3,60) [b,a]=butter(N,Wc,high) freqz(b,a); axis([0,1,-120,0]); grid on title( 巴特沃斯高通数字滤波器 ) 巴特沃斯高通数字滤波器0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Normalized Frequency ( rad/sample)200)seerged(esah0 巴特沃斯高通数字滤波器 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Normalized Frequency ( rad/sample) 200 )seerged