云南师范大学dsp(数字信号处理)试验报告9资料.doc

本科学生实验报告 学号 114090395 姓名 李华福 学院 物电学院 专业、班级 11电子 实验课程名称 数字信号处理（实验） 教师及职称 李宏宁 开课学期 第三 至 第四 学年 下 学期 填报时间 2014 年 5 月 14 日 云南师范大学教务处编印 实验设计方案 实验序号 09 实验名称 IIR数字滤波器设计及应用 实验时间 2014年5月14日 实验室 数字信号处理实验室（313） 实验目的 加深理解IIRIIR数字滤波器的设计原理与设计方法，以及IIR数字滤波器的应用。 2． 实验原理、实验流程或装置示意图 IIR数字滤波器一般为线性时不变的因果离散系统，N阶IIR数字滤波器的系统函数可表达为Z-1 的有理项式，即： 系数{ai；i1,2,…,N}至少有一个非零。对于因果IIR数字滤波器，满足M=N。 IIR数字滤波器的设计主要通过成熟的模拟滤波器设计方法来实现：将数字滤波器设计指标转换为模拟滤波器设计指标，设计出相应的模拟滤波器H(s)，再经过脉冲响应不变法或双线性变换法得到所需的IIR数字滤波器H(z)。 IIR数字滤波器设计的重要环节是模拟原型低通滤波器的设计，主要包括Butterworth、Chebyshev和椭圆等滤波器。 MATLAB 信号处理工具箱中提供了IIR滤波器设计的函数。 IIR 滤波器阶数选择 buttord -巴特沃斯(Butterworth)滤波器阶数选择。 cheb1ord -切比雪夫(Chebyshev)I 型滤波器阶数选择。 cheb2ord -切比雪夫(Chebyshev)II 型滤波器阶数选择。 ellipord -椭圆(Elliptic)滤波器阶数选择。 IIR 滤波器设计 butter -巴特沃斯(Butterworth)滤波器设计 cheby1 -切比雪夫(Chebyshev)I 型滤波器设计 cheby2 -切比雪夫(Chebyshev)II 型滤波器设计 ellip -椭圆(Elliptic)滤波器设计 maxflat -通用的巴特沃斯(Butterworth)低通滤波器设计 yulewalk -Yule-Walker 滤波器设计(直接数字滤波器设计法) Butterworth滤波器设计 Butterworth滤波器是通带、阻带都单调衰减的滤波器。 调用buttord函数确定巴特沃斯滤波器的阶数，格式为 [N,Wc] = buttord(Wp,Ws,Ap,As) 输入参数：Ap,As为通带最大衰减和阻带最小衰减，以dB为单位。Wp,Ws为归一化通带截频和阻带截频，0Wp,Ws1 。 输出参数：N为滤波器的阶数；Wc为截频，0 Wc 1。 调用butter函数设计出巴特沃斯滤波器，格式为[b,a] = butter(N,Wc,options) 输入参数：N和Wc是buttord函数返回的参数，含义见上。Options=’low’,’high’,’bandpass’,’stop’,分别对应低通、高通、带通、带阻，默认情况下为低通或带通。 输出参数：b和a为设计出的IIR数字滤波器H(s)的分子多项式和分母多项式的系数矩阵。 【例2.3.1】 设计一个数字滤波器，要求在0~0.2（rad）内衰耗不大于3dB，在0.6~ （rad）内衰耗不小于60dB。 [解] 若模拟低通滤波器原型采用巴特沃思滤波器： [N,Wc]=buttord(0.2,0.6,3,60); [b,a]=butter(N,Wc); reqz(b,a);axis([0,1,-120,0]); 设计结果如图2.3.1所示。可见IIR是数字滤波器相位为非线性。 图2.3.1 Buttertworth低通滤波器 (2) Chebyshev I型滤波器设计 Chebyshev I型滤波器为通带纹波控制器：在通带呈现纹波特性，在阻带单调衰减。 [N,Wc] = cheb1ord(Wp, Ws, Ap, As) [b,a] = cheby1(N,Ap,Wc,options) 参数含义与butter中参数一致。 【例2.3.2】 设计一个数字带通滤波器，要求在100~200Hz通带内纹波不大于3dB，