第5章--IIR-滤波器设计1.ppt

第4章IIR滤波器设计 葛青 4 1滤波器的概念及分类 数字滤波器 通过对数字信号的运算处理 改变信号频谱 完成滤波作用的算法或装置 数字滤波器具有稳定 适应性强 线性相位等优点 避免了模拟滤波器固有的电压漂移 温度漂移和噪声等问题 2 滤波原理 线性时不变系统输入输出关系可用如下数学函数来描述时域 频域 幅频特性 表示信号通过该滤波器后各频率成分的衰减情况相频特性 反映各频率成分通过滤波器后在时间上的延时情况 联系无失真传输条件 2 滤波原理 滤波原理 2 滤波原理 滤波器的设计就是用一个因果稳定的线性移不变系统的系统函数去逼近理想滤波器的性能 数字滤波器设计 确定传输函数H z 的过程 滤波器的设计原理 表征滤波器频率响应的特征参量 3个参量 1 幅度平方响应 幅度的平方与频率的关系 2 相位响应 3 群延迟响应 物理意义 相位对角频率的导数的负值 若滤波器通带内群延迟响应 常数 则为线性相位滤波器 联系无失真传输问题 无失真传输系统 如果系统响应仅是输入信号在时间上的延时和幅度上的放大或缩小 就认为信号在传输过程中没有失真 这种系统称为无失真传输系统 无失真传输系统 3 滤波器的分类 经典滤波器 用于提取有用的信号 抑制不需要的干扰 或将信号分成不同的频率成分 现代滤波 信号检测和信号估计 经典滤波器和现代滤波器 模拟和数字滤波器的理想幅频特性曲线 模拟滤波器幅频特性曲线 数字滤波器幅频特性曲线 低通滤波器 LPAF LPDF Lowpassfilter带通滤波器 BPAF BPDF Bandpassfilter高通滤波器 HPAF HPDF Highpassfilter带阻滤波器 BSAF BSDF Bandstopfilter模拟 Analog 滤波器 是用来处理模拟信号的处理系统 数字 Digital 滤波器 是用来处理数字信号的处理系统 IIR infiniteimpulseresponse 有限冲激响应滤波器 单位冲激响应h n 是无限长序列FIR finiteimpulseresponse 无限冲激响应滤波器 单位冲激响应h n 是有限长序列 数字滤波器的分类 4 IIR系统与FIR系统 从系统的单位采样响应来区分无限长单位冲激响应 IIR 系统 单位冲激响应h n 是无限长序列有限长单位冲激响应 FIR 系统 单位冲激响应h n 是有限长序列 五 IIR系统与FIR系统 从系统函数的构造来区分IIR系统 至少有一个极点 包括全极点系统 分子只有常数项 和零极点系统 分子不止常数项 有反馈环路 采用递归型结构 FIR系统 收敛域内无极点 是全零点系统 无反馈环路 多采用非递归结构 IIR系统与FIR系统 从系统函数的构造来区分 IIR系统与FIR系统 当0 n N 1时 h n 0具有这种形式的滤波器称为FIR滤波器 FiniteImpulseResponse 其差分方程为 2 bi中至少有一项不为零 并且分母中至少有一个根不为分子所抵消 具有这种形式的滤波器称为IIR滤波器 InfiniteImpulseResponse 其差分方程为 IIR系统与FIR系统 IIR系统 至少有一个极点 包括全极点系统 分子只有常数项 和零极点系统 分子不止常数项 有反馈环路 采用递归型结构 FIR系统 收敛域内无极点 是全零点系统 无反馈环路 多采用非递归结构 4 数字滤波器的性能要求 一个理想滤波器 要求所在通频带内幅频响应是一常数 相位频率相应为零或是频率的线性函数 但一个实际的滤波器要是不可能得到上述幅频和相频响应 以低通滤波器为例 频率响应有通带 过渡带及阻带三个范围 4 数字滤波器的性能要求 p 通带衰减 s 阻带衰减 p 通带截止频率 s 阻带截止频率 p s 过渡带 通带允许的最大衰减阻带最小衰减 4 数字滤波器的性能要求 式中假定被归一化为1 按设计任务 确定滤波器性能要求 制定技术指标一个因果稳定的离散LSI系统的系统函数H z 逼近此性能指标利用有限精度算法实现此系统函数 如运算结构 字长的选择等实际技术实现 软件法 硬件法或DSP芯片法 数字滤波器的设计过程 数字滤波器的设计 IIR滤波器设计主要内容包括 巴特沃思 切比雪夫模拟低通滤波器设计 脉冲响应不变法和双线性变换法的数字化变换方法 数字高通 带通和带阻滤波器的设计 而FIR滤波器是直接采用的数字式设计方法 针对FIR滤波器特征 首先介绍了其线性相位的实现条件 然后介绍了窗函数法和频率抽样法的设计方法 设计IIR数字滤波器一般有以下两种方法 1 模拟滤波器 首先设计一个合适的模拟滤波器 然后将它转换成满足给定指标的数字滤波器 这种方法适合于设计幅频特性比较规则的滤波器 例如低通 高通 带通 带阻等 2 直接在频域或者时域中进行数字滤波器设计 由于要联立方程 设计时需要计算机作辅助设计 4 2IIR数字滤波器的设计 IIR滤波器以模拟低通滤波器为基础的设计方法 为了设计其他的选频滤波器 高通 带通 带阻等 需要对低通滤波器进行频率转换 在设计过程中有两种不同的变换 频带变换和模拟 数字变换 根据这两种变换的先后次序 引出两种设计方法 4 2IIR数字滤波器的设计 其它非低通数字滤波器的设计过程 IIR的设计方法 借用模拟滤波器 先设计一个合适的模拟滤波器 然后变换成满足给定指标的数字滤波器 设计步骤将给定的数字滤波器的技术指标转换为模拟滤波器的技术指标 根据转换后的技术指标设计模拟原型滤波器 按照一定规则将模拟滤波器转换为数字滤波器 4 2 1模拟低通滤波器的设计 常用模拟低通滤波器及其特点1 巴特沃思 Butterworth 滤波器 特点是从通带中心向两边幅频特性单调下降 2 切比雪夫型I Chebyshev 滤波器 通带内有等波纹幅频特性函数 阻带内有单调下降的幅频特性函数 3 切比雪夫型II Chebyshev 滤波器 通带内有单调下降的幅频特性函数 阻带内有等波纹幅频特性函数 4 椭圆滤波器 又称考尔Cauer滤波器或双切比雪夫滤波器 通带和阻带内均有等波纹幅频特性 5 贝塞尔 Bessel 滤波器 又称最大平坦群时延滤波器 在中心频率附近逼近线性相位特性 模拟滤波器的设计过程 1 根据要求确定设计指标 2 选择滤波器类型 3 计算滤波器阶数 4 通过查表或计算确定滤波器的系统函数Ha s 5 综合实现并装配调试 4 2 1模拟低通滤波器的设计 1 巴特沃思低通滤波器 1 基本性质巴特沃思滤波器以巴特沃思函数来近似滤波器的系统函数 巴特沃斯滤波器是根据幅频特性在通频带内具有最平坦特性定义的滤波器 N为滤波器阶数 c为3dB截止频率 巴特沃思滤波器的低通模平方函数表示 巴特沃思低通滤波平方幅频特性函数 巴特沃斯滤波器的主要特征a 对所有的N 即 b 是的单调下降函数 c 随着阶次N的增大而更接近于理想低通滤波器 归一化原型滤波器是指截止频率已经归一化成的低通滤波器 对于截止频率为某个低通滤波器 则令替归一化原型滤波器系统函数中的 对于其他高通 带通 带阻滤波器 可应用后面讨论到的频带变换法 由其变换得出 巴特沃斯低通滤波器的幅度平方函数为 N为滤波器的阶次 为3dB截频 2 巴特沃思低通滤波器设计过程 通带最大衰减和阻带最小衰减可由定义式求出 设 0处幅度归一化到1 则 2 巴特沃思低通滤波器设计过程 1 求出N Nmin应取向上取整 2 如技术指标未给出 则可由下式计算 或采用上式求 阻带指标有富裕量 下式则通带指标有富裕量 查表求得归一化传输函数 令代替归一化原型滤波器系统函数中的 即得到实际滤波器传输函数Ha s 3 求系统函数Ha s 低通巴特沃斯滤波器设计步骤 1 由求滤波器阶次N 2 由N查表 求出归一化传输函数 3 令代入 得实际滤波器传输函数Ha s Nmin应取向上取整 N及计算 或 通带边频 通带的最大衰减为 阻带边频为 阻带的最小衰减为 例 设计一巴特沃思低通滤波器 使其满足以下指标 解 根据滤波器技术指标 取N 4 查表得四阶巴特沃思多项式 得归一化系统函数表达式 确定阶次N 巴特沃斯归一化低通滤波器参数 巴特沃斯归一化低通滤波器参数 巴特沃斯归一化低通滤波器参数 代入 1 巴特沃思滤波器设计步骤 1 根据技术指标求求出滤波器的阶数N和3dB截止频率 Nmin应取向上取整 或 1 巴特沃思滤波器设计步骤 2 求系统函数Ha s 1 查表求Ha s 查表求得归一化传输函数 令代替归一化原型滤波器系统函数中的 得到实际滤波器传输函数 求滤波器的极点 并由s平面左半平面的极点构成 幅度平方函数极点为 将s平面左半平面的N个极点分配给 2 极点法求系统函数Ha s 例 已知fp 5kHz p 2dB fs 12kHz s 30dB 设计巴特沃斯低通滤波器 解 1 确定阶次N 取N 5 其极点为 归一化极点 2 查表求极点 P1 2 0 3090 j0 9511 p3 4 0 8090 j0 5857 p0 1 000所以 3 求 c 所以 s比题目中给出的指标小 过渡带符合要求 4 去归一化 求Ha S 将p s c代入Ha P 得 切比雪夫滤波器引入原因 Butterworth滤波器频率特性 无论在通带与阻带都随频率单调变化 因此如果在通带边缘满足指标 则在阻带内肯定会有富裕量 也就是会超过指标的要求 因此并不经济 更有效的方法是将指标的精度要求均匀的分布在通带内 或均匀分布在阻带内 或同时均匀分布在通带与阻带内 这时就可设计出阶数较低的滤波器 这种精度均匀分布的办法可通过选择具有等波纹特性的逼近函数来完成 4 2 2切比雪夫滤波器设计 通带内等波纹 阻带内单调 切比雪夫1型通带内单调 阻带内等波纹 切比雪夫2型切比雪夫1型的幅度平方函数 是小于1的正数 称为纹波参数 是表示通带内纹波大小的一个参数 愈大 纹波也愈大 p称为有效通带截止频率 频率通常对 p归一化 在通带内 在1和之间起伏变化 在阻带内是单调下降的 当N为奇数时 滤波器在处的幅度响应为1 当N为偶数时 滤波器在处的幅度响应为 当时 4 2 2切比雪夫滤波器设计 基本性质 的确定 设允许的通带纹波为 那么 阶数N的确定 1 查表得 2 查表得幅度平方函数的极点 其中 的确定 极点分布 2N个等角度间隔 间隔为 N 分布在为长半轴 为短半轴的椭圆上的点 切比雪夫滤波器的传输函数 系数A由s 0时滤波器幅度响应的值确定 当N为奇数时 当N为偶数时 去归一化得到滤波器的传输函数Ha s 切比雪夫低通滤波器的设计步骤 1 由给定的设计指标确定切比雪夫滤波器的参数 N和 p 2 确定Ha s 1 查表求得归一化传输函数 令代替归一化原型滤波器系统函数中的 即得到实际滤波器传输函数 切比雪夫低通滤波器的设计步骤 归一化Chebyshev低通滤波器的传输函数为 切比雪夫低通滤波器的设计步骤 归一化Chebyshev低通滤波器的传输函数为 为了得到单位峰值增益 一般选 切比雪夫低通滤波器的设计步骤 2 确定Ha s 2 滤波器的极点 并由s平面左半平面的极点构成Ha s 切比雪夫低通滤波器的设计步骤 系数A由s 0时滤波器幅度响应的值确定 当N为奇数时 当N为偶数时 切比雪夫低通滤波器的设计步骤 去归一化得到滤波器的传输函数Ha s 切比雪夫滤波器是一种全极型滤波器 它是由切比雪夫多项式的正交函数推导出来的 采用了在通带内等波动 在通带外单调衰减的准则去逼近理想特性 它比巴特沃斯滤波器在通带内有更均匀的特性 是在所有全极型滤波器中过渡带最窄的最