实验十一IIR数字滤波器设计及应用.docxVIP

物理与电子信息学院----电子信息技术实验中心 数字信号处理及实验实验报告 实验题目 IIR数字滤波器设计与应用 姓名 组别 班级 光电14 学号 1443202000206 【实验目的】 加深理解IIR数字滤波器的特性，掌握IIR数字滤波器的设计原理与设计方法，以及IIR数字滤波器的应用。 【实验原理】 IIR数字滤波器一般为线性时不变的因果离散系统，N阶IIR数字滤波器的系统函数可以表达为z-1 H 系数{ai；i=1， 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 滤波器设计(直接数字滤波器设计法) 1.Butterworth滤波器设计 Butterworth滤波器是通带、阻带都单调衰减的滤波器。 (1) 调用buttord函数确定巴特沃斯滤波器的阶数，格式为 [N,Wc] = buttord(Wp,Ws,Ap,As) 输入参数： Ap,As为通带最大衰减和阻带最小衰减，以dB为单位。 Wp,Ws为归一化通带截频和阻带截频，归一化频率=实际数字频率π。因此0Wp，Ws1，1对应数字频率π 输出参数： N为滤波器的阶数； Wc为截频，0 Wc 1。 (2) 调用butter函数设计出巴特沃斯滤波器，格式为 [b,a] = butter(N,Wc,options) 输入参数： N和Wc是buttor函数返回的参数，含义见上。 Options=’low’, ’high’, ’bandpass’, ’stop’, 分别对应低通、高通、带通、带阻，默认情况下为低通或带通。 输出参数： b和a为设计出的IIR数字滤波器H(s)的分子多项式和分母多项式的系数矩阵。 2.Chebyshev I型滤波器设计 Chebyshev I型滤波器为通带纹波控制器：在通带呈现纹波特性，在阻带单调衰减。 [N,Wc] = cheb1ord(Wp, Ws, Ap, As) [b,a] = cheby1(N,Ap,Wc,options) 参数含义与butter中参数一致。 3. Chebyshev II 型滤波器设计 Chebyshev II 型滤波器为阻带纹波控制器:在阻带呈现纹波特性。在通带单调衰减。 [N,Wc] = cheb2ord(Wp, Ws, Ap, As) [b,a] = cheby2(N,As,Wc,options) 4. 椭圆滤波器设计 椭圆滤波器在通阻带都呈现纹波特性。 [N,Wc] = ellipord(Wp,Ws,Ap,As) [b,a] = ellip(N,Ap,As,Wc,options) 【实验结果与数据处理】 1、信号x[k]=1+cos(π4k)+cos(2π （1）设计IIR低通滤波器，滤除cos(2π MATLAB程序： clc,clear,close all [N,Wc]=b