第六章 无限脉冲响应数字滤波器的的设计.ppt

第六章 无限长单位脉冲响应(IIR)  数字滤波器的设计方法 本章目录 数字滤波器的技术指标与设计方法 说明 数字低通滤波器技术指标 通带 阻带 过渡带 数字低通滤波器技术指标 通带内允许的最大衰减 阻带内允许的最小衰减 概念：片段常数特性 对于选频滤波器，一般对通带和阻带内的幅频响应曲线形状没有具体要求，只要求其波纹幅度小于某个常数，通常将这种要求称为“片段常数特性”。 所谓片段，是指“通带”和“阻带”，常数是指“通带波纹幅度δ1”和“阻带波纹幅度δ2”。 三、IIR数字滤波器的两类设计方法 6.2 模拟滤波器的设计 模拟低通滤波器的设计，即设计Ha(s)； 模拟滤波器的频率变换 （即模拟高通、带通、带阻滤波器的设计） 6.2.1 模拟低通滤波器的设计 常用的模拟滤波器 巴特沃斯（Butterworth）滤波器 具有单调下降的幅频特性 切比雪夫（Chebyshev）滤波器 幅频特性在通带或者阻带内有波动 椭圆（Ellipse）滤波器 在通带和阻带内都有纹波 贝塞尔（Bessel）滤波器等 通带内有较好的线性相位特性 理想模拟滤波器幅频特性 一、模拟滤波器的技术指标 幅度平方函数 模拟低通滤波器的技术指标 通带截止频率 通带最大衰减 阻带截止频率 阻带最小衰减 3dB截止频率 二、模拟滤波器的设计思路 Ha(s) 由给定的 、 、 和 求出 三、巴特沃斯低通滤波器的设计方法 巴特沃斯低通滤波器的幅度平方函数为： N为滤波器阶数 Ωc为3dB截止频率 巴特沃斯低通滤波器的特点 幅度特性随着Ω增加单调下降，下降的速度与阶数有关。 随着N增大，幅度下降的速度越快，过渡带越窄，在通带内更接近于1，在阻带内迅速接近于零，因而幅度特性更接近于理想的矩形频率特性。 不管N的取值是多少，都经过 点。 幅度平方函数的极点分布 幅度平方函数有2N个极点 这2N个极点等间隔分布在半径为Ωc的圆上（该圆称为巴特沃斯圆），间隔是?/Nrad。 这些极点以虚轴为对称轴，而且不会落在虚轴上。 Butterworth圆 求得Ha(s) 在S平面总共有2N个极点，为了保证所设计的滤波器是稳定的，将s平面左半平面的N个极点分配给Ha(s)，而将右半平面的N个极点分配给Ha(-s)，得： 对频率以Ωc归一化 Ha(p)的求得 Ha(s)的求得 滤波器阶数N的确定 滤波器Ωc的确定 四、模拟Butterworth低通滤波器设计步骤 由给定的设计指标Ωp、α p 、 Ωs和αs确定巴特沃斯滤波器的阶数N和频率Ωc 。 由阶数N，按所要求的形式，查表6.2.1,得到Ha(p)。 由Ha(p)去归一化得到Ha(s)。 已知通带截止频率fp=5kHz，通带最大衰减αp=2dB，阻带截止频率fs=12kHz，阻带最小衰减αs=30dB，按照以上技术指标设计巴特沃斯低通滤波器。 (3) 将Ha(p)去归一化，得到Ha(s) 巴特沃斯滤波器的MATLAB实现 IIR数字滤波器的Matlab仿真实现 IIR模拟滤波器设计 模拟滤波器到数字滤波器的转换 IIR数字滤波器设计 设数字滤波器系统函数为 模拟滤波器的系统函数为 函数buttord的调用格式 函数buttord的调用格式为 [n, Wc] = buttord(Wp, Ws, Rp, Rs) % 设计数字Butterworth滤波器 [n, Wc]= buttord(Wp, Ws, Rp, Rs,’s’) % 设计模拟Butterworth滤波器 其中，n为滤波器阶数，Wc为3dB截止频率。 Lowpass: Wp = .1, Ws = .2 Highpass: Wp = .2, Ws = .1 Bandpass: Wp = [.2 .7], Ws = [.1 .8] Bandstop: Wp = [.1 .8], Ws = [.2 .7] 函数butter的调用格式 函数butter的调用格式为 [b, a]= butter(n, Wc,‘ftype’,‘s’) % 设计模拟Butterworth滤波器 [b, a]= butter(n, Wc,‘ftype’) % 设计数字Butterworth滤波器 其中‘ftype’定义滤波器的类型； Wc为一个变量时，默认为低通，高通时‘high’ Wc为二元矢量时，默认为带通，带阻时‘stop’ 例6.2.1的MATLAB程序 wp=2*pi