实验六(IIR滤波器设计)分析.ppt

实验目的 加深理解IIR数字滤波器的特性，掌握IIR数字滤波器的设计原理与设计方法，以及IIR数字滤波器的应用。 实验六 IIR数字滤波器设计及应用 N阶IIR数字滤波器的系统函数为: IIR数字滤波器的设计主要通过成熟的模拟滤波器设计方法来实现，将数字滤波器设计指标转换为模拟滤波器设计指标，设计出相应的模拟滤波器Ha(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滤波器是通带、阻带均单调衰减的滤波器。 ? 调用buttord函数确定巴特沃斯滤波器的阶数，格式为 [N,Wc] = buttord(Wp, Ws, δp, δs) 输入参数：δp, δs分别为通带最大衰减和阻带最小衰减，以dB为单位。 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(z)的分子多项式和分母多项式的系数矩阵。 二. 实验原理 设计一个数字低通滤波器，要求在0~0.2?(rad)内衰减不大于3dB，在0.6?~?(rad)内衰减不小于60dB。 若模拟低通原型采用Butterworth滤波器： [N,Wc]=buttord(0.2,0.6,3,60); [b,a]=butter(N,Wc); freqz(b,a); axis([0,1,-120,0]); grid on 二. 实验原理 设计结果： 二. 实验原理 2. Chebyshev I型滤波器设计 Chebyshev I型滤波器为通带波纹滤波器：在通带呈现纹波特性，在阻带单调衰减。 [N,Wc] = cheb1ord(Wp, Ws, δp, δs) [b,a] = cheby1(N, δp, Wc, options) 参数含义与butter中参数一致。 二. 实验原理 设计一个数字带通滤波器，在100Hz ~ 200Hz通带内纹波不大于3dB，通带两边各50 Hz外是阻带，衰耗不小于40dB。抽样频率为1000 Hz 。 Wp =[100 200]/500; Ws =[100-50 200+50]/500; Rp = 3; Rs = 40; [N,Wn] = cheb1ord(Wp,Ws,Rp,Rs) [b,a] = cheby1(N,Rp,Wn); freqz(b,a,512,1000); title(Chebyshev Type I Bandpass Filter) axis([0,500,-80,0]); grid on 二. 实验原理 二. 实验原理 3. Chebyshev II 型滤波器设计 Chebyshev II 型滤波器为阻带波纹