IIR数字滤波器及设计方法.ppt

取 N = 2 6.6 设计IIR滤波器的频率变换法 设计各种IIR数字滤被器，通常是先设计一个归一化截止频率的模拟低通原型滤波器，然后用如下两种方法转换为数字滤波器： 模拟频带变换法 数字频带变换法 将归一化的模拟低通原型滤波器经模拟频带变换成所需类型的模拟滤波器，然后用冲激响应不变法或双线性变换法转换为所需类型的数字滤波器 将归一化的模拟低通原型滤波器通过冲激响应不变法或双线性变换法转换为数字低通滤波器，然后用数字频带变换法，转换为所需类型数字滤波器 6.6 设计IIR滤波器的频率变换法 模拟频带变换法 两种等效的设计方法 模拟域频率变换关系可通过查表得到 截止频率为Ωc的模拟低通滤波器到其它频率选择性滤波器的转换公式 一般直接用归一化原型转换，取Ωc=1, 可使设计过程简化 数字化可采用脉冲响应不变法与双线性变换法 把第一种方法中的两步合成一步来实现 把数字滤波器的性能要求转换为对应的模拟滤波器的性能要求 根据性能要求设计模拟滤波器 通过脉冲响应不变法或双线性变换法，转化为数字滤波器 6.6.1 模拟低通滤波器变换成数字低通滤波器 设计步骤 由于采用脉冲响应不变法设计，因此模拟低通的 fc=1 kHz 例 用脉冲响应不变法设计三阶巴特沃思数字低通滤波器，采样频率为fs=4 kHz，其3 dB截止频率为fc=1 kHz 已知所求数字巴特沃思低通滤波器 fs=4 kHz，fc=1 kHz 设计模拟低通滤波器 查表得归一化三阶模拟巴特沃思原型低通滤波器的传递函数 去归一化 ， s=s/Ωc = s/(2πfc )代入上式 通过脉冲响应不变法转化为数字滤波器 将上式因式分解 将部分分式系数代入下式 该映射变换须满足以下要求：  S平面虚轴 jΩ必须映射到Z平面的单位圆ejω上 S平面的左半平面必须映射到Z平面单位圆的内部|z|1 利用模拟滤波器来设计数字滤波器 常用的变换方法有：脉冲响应不变法和双线性变换法 使稳定的模拟系统映射为具有相同幅频特性的稳定的数字系统 6.4 用脉冲响应不变法设计IIR数字滤波器 使数字滤波器的单位脉冲响应序列h(n)与相应的模拟滤波器的冲激响应ha(t)，在采样点处的量值相等，即h(n)=ha(nT) 6.4.1 变换原理 第1章中采样序列的Z变换与模拟信号拉氏变换的关系（P48） 脉冲响应不变法将模拟滤波器的S平面变换成数字滤波器的Z平面 脉冲响应不变法的映射关系 数字滤波器与模拟滤波器的频率响应间的关系为 (P49) 数字滤波器的频率响应是模拟滤波器频率响应的周期延拓 只有当模拟滤波器的频率响应是限带的，且满足采样定理时，数字滤波器的频率响应才不产生混叠失真 提高采样频率为fs，可减小频率响应的混叠效应 6.4.2 混叠失真 脉冲响应不变法中的频响混叠现象 脉冲响应不变法设计中，由一个较为复杂的模拟系统函数求出数字滤波器系统函数是一个很麻烦的变换过程 因为乘积的z变换并不等于各部分变换的乘积，所以不宜采用级联分解 但各项和的z变换是线性关系，因而用部分分式表达系统函数，特别适合于对复杂模拟系统函数的变换 6.4.3 模拟滤波器的数字化方法 设模拟滤波器的系统函数Ha(s)只有单阶极点，且分母的阶次大于分子的阶次，将Ha(s)用部分分式表示 数字滤波器的单位脉冲响应等于对ha(t)的采样，因此 将Ha(s) 拉氏反变换 u(t)是单位阶跃函数 模拟滤波器与数字滤波器系统函数的关系 对h(n)求Z变换，得数字滤波器的系统函数 S平面的单极点s=sk对应于Z平面上z=eskT处的单极点 Ha(s)与H(z)的部分分式的系数相同，都是Ak 模拟滤波器因果稳定，则数字滤波器也因果稳定  脉冲响应不变法只能保证极点有这种代数对应关系 与模拟滤波器的系统函数Ha(s) 相比 Ha(s)与H(z)间的变换关系为 脉冲响应不变法的设计步骤，可直接将 Ha(s)写成单极点的部分分式和的形式，然后将各部分分式用上式关系替代，即得到H(z) 由于数字滤波器频率响应幅度与采样间隔T成反比(p49) |ω|π 为使数字滤波器增益不随采样频率变化，通常令 h(n)=Tha(nT) 此时: 将上式代入到 例 设模拟滤波器的系统函数为 试利用脉冲响应不变法将Ha(s)转换成IIR数字滤波器的系统函数H(z) 解 数字滤波器的系统函数为 设T=1，则有 模拟滤波器的频响Ha(jΩ)以及数字滤波器的频响H(ejω)分别为: 由于Ha(jΩ)不是限带的，所以H(ejω)产生了频谱混叠失真 T取不同值时的幅频特性图 数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应，因此时域特性逼近好