IIR数字滤波器-精简.ppt

第一节数字滤波器的基本概念 数字滤波器 输入输出均为数字信号 经过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件 优点 精度高 稳定 体积小 重量轻 灵活 不要求阻抗匹配 能实现模拟滤波器 AF 无法实现的特殊滤波功能 第6章IIR数字滤波器设计 一 数字滤波器的分类 滤波器的种类很多 分类方法也不同 1 从功能上分 低通 带通 高通 带阻 2 从实现方法上分 FIR IIR 3 从设计方法上来分 Butterworth 巴特沃斯 Chebyshev 切比雪夫 Ellipse 椭圆 等 4 从处理信号分 经典滤波器 现代滤波器 1 经典滤波器 假定输入信号x n 中的有用成分和希望去除的成分 各自占有不同的频带 当x n 经过一个线性系统 即滤波器 后即可将欲去除的成分有效地去除 如果信号和噪声的频谱相互重叠 那么经典滤波器将无能为力 此时可以设计现代滤波器来解决 2 模拟滤波器和数字滤波器 经典滤波器从功能上分又可分为 1 低通滤波器 LPAF LPDF Lowpassanalogfilter Lowpassdigitalfilter 2 高通滤波器 HPAF HPDF Highpassanalogfilter Highpassdigitalfilter 3 带通滤波器 BPAF BPDF Bandpassanalogfilter Bandpassdigitalfilter 4 带阻滤波器 BSAF BSDF Bandstopanalogfilter Bandstopdigitalfilter 3 模拟滤波器的理想幅频特性 4 数字滤波器的理想幅频特性 用一个因果稳定的离散LSI系统的系统函数H z 逼近此性能指标 按设计任务 确定滤波器性能要求 制定技术指标 利用有限精度算法实现此系统函数 如运算结构 字长的选择等 实际技术实现 软件法 硬件法或DSP芯片法 1 数字滤波器的设计过程 二 数字滤波器的设计和指标 2 数字滤波器的技术要求 我们通常用的数字滤波器一般属于选频滤波器 假设数字滤波器的传输函数H ej 用下式表示 幅频特性 H ej 信号通过滤波器后的各频率成分衰减情况 相频特性 各频率成分通过滤波器后在时间上的延时情况 理想滤波器不可实现 只能以实际滤波器逼近 通带 阻带 过渡带 常用的技术指标是通带内和阻带内允许的衰减 一般用dB数表示 通带内允许的最大衰减用 p表示 阻带内允许的最小衰减用 s表示 p和 s分别定义为 上式中 归一化为1 H z H z 1 的极点既是共轭的 又是以单位圆成镜像对称的 H z 的极点 单位圆内的极点 3 表征滤波器频率响应的特征参量 幅度平方响应 定义 相位对角频率的导数的负值 群延迟响应 若滤波器通带内为常数 则为线性相位滤波器 相位响应 相位响应 用一因果稳定的离散LSI系统逼近给定的性能要求 s平面逼近 模拟滤波器的设计 z平面逼近 数字滤波器的设计 4 IIR数字滤波器的设计方法 思路 先设计模拟滤波器 再转换为数字滤波器 即为求滤波器的各系数 第二节模拟滤波器的设计 模拟滤波器的理论和设计方法已发展得相当成熟 且有若干典型的模拟滤波器供我们选择 如 1 巴特沃斯 Butterworth 滤波器2 切比雪夫 Chebyshev 滤波器3 椭圆 Ellipse 滤波器4 贝塞尔 Bessel 滤波器这些滤波器都有严格的设计公式 现成的曲线和图表供设计人员使用 1 模拟低通滤波器的设计指标及逼近方法 模拟低通滤波器的设计指标有 p p s和 s p和 s分别称为通带截止频率和阻带截止频率 p是通带 0 p 中的最大衰减系数 s是阻带 s的最小衰减系数 c称为3dB截止频率 滤波器的技术指标给定后 需要设计一个传输函数Ha s 希望其幅度平方函数满足给定的指标 p和 s 一般滤波器的单位冲激响应为实数 有 2 巴特沃斯低通滤波器的设计方法 1 巴特沃斯低通滤波器的幅度平方函数 Ha j 2用下式表示 N称为滤波器的阶数N越大 越接近理想滤波器 N越大 滤波器的实现也越复杂 低通巴特沃斯滤波器的设计步骤如下 根据技术指标 p p s s 求出滤波器的阶数N N的大小主要影响幅度特性下降的速度 2 查表得到归一化的巴特沃斯低通滤波器传输函数Ha p 3 将Ha p 去归一化 由巴特沃斯低通滤波器的幅度平方函数求出 c 将p s c代入Ha p 得到实际滤波器传输函数Ha s 求出的N可能有小数部分 取大于等于N的最小整数 例 已知通带截止频率fp 5kHz 通带最大衰减 p 2dB 阻带截止频率fs 12kHz 阻带最小衰减 s 30dB 按照以上技术指标设计巴特沃斯低通滤波器 解 1 确定阶数N 由N 5 直接查表 p261 表6 4 得到 极点 0 3090 j0 9511 0 8090 j0 5878 1 0000 其中 b0 1 0000 b1 3 2361 b2 5 2361 b3 5 2361 b4 3 2361 3 为将Ha p 去归一化 先求3dB截止频率 c 将p s c代入Ha p 中得到 检验 可以看出 满足 s 30dB的真实fs在10 525kHz处 与12kHz比 还有富裕量 第三节用脉冲 冲激 响应不变法设计IIR数字低通滤波器 利用模拟滤波器来设计数字滤波器 就是从已知的模拟滤波器传递函数Ha s 设计数字滤波器的系统函数H z 它归根结底是一个由S平面映射到Z平面的变换 这个变换通常是复变函数的映射变换 为了保证转换后的H z 稳定且满足技术要求 这个映射变换必须满足以下两条基本要求 2 数字滤波器的频率响应模仿模拟滤波器的频响 s平面的虚轴映射为z平面的单位圆 相应的频率之间成线性关系 因果稳定的模拟滤波器转换成数字滤波器 仍是因果稳定的 即S平面的左半平面必须映射到Z平面单位圆的内部 一 脉冲响应不变法的转换原理 核心原理 通过对连续函数ha t 等间隔采样得到离散序列ha nT 令h n ha nT T为采样间隔 它是一种时域上的转换方法 转换步骤 设模拟滤波器的传输函数为Ha s 相应的单位冲激响应是ha t 设模拟滤波器Ha s 只有单阶极点 且分母多项式的阶次高于分子多项式的阶次 将Ha s 用部分分式表示 式中si为Ha s 的单阶极点 将Ha s 进行逆拉氏变换得到ha t 式中u t 是单位阶跃函数 对ha t 进行等间隔采样 采样间隔为T 得到 对上式进行Z变换 得到数字滤波器的系统函数H z 设ha t 的采样信号用表示 二 存在的问题 其拉氏变换 采样信号的拉氏变换与相应的序列的Z变换之间的映射关系 多值映射 存在的问题 频率响应混叠失真 数字滤波器的频率响应是模拟滤波器频率响应的周期延拓 脉冲响应不变法的频率混叠现象 只有当模拟滤波器的频率响应是限带的 且带限于折叠频率以内时 才能使数字滤波器的频率响应在折叠频率以内重现模拟滤波器的频率响应 而不产生混叠失真 即 但实际的模拟滤波器频率响应都不是严格限带的 非理想 变换后就会产生周期延拓分量的频谱交叠 即产生频率响应的混叠失真 这时数字滤波器的频响就不同于原模拟滤波器的频响 有一定的失真 当模拟滤波器的频率响应在折叠频率以上处衰减越大 越快时 变换后频率响应混叠失真就越小 这时 采用脉冲响应不变法设计的数字滤波器才能得到良好的效果 数字滤波器频率响应幅度还与采样间隔T成反比 如果采样频率很高 即T很小 数字滤波器具有太高的增益 这是不希望的 为了使数字滤波器增益不随采样频率而变化 可以作以下简单的修正 令 h n Tha nT 则有 很小的数在分子比在分母安全 例 已知模拟滤波器的传输函数Ha s 为 用脉冲响应不变法将Ha s 转换成数字滤波器的系统函数H z 解 首先将Ha s 写成部分分式 极点为 那么H z 的极点为 设T 1s时用H1 z 表示 则 T 0 1s时用H2 z 表示 则 模拟滤波器 通带很窄 阻带拖了很长的尾巴 并非严格限带数字滤波器 T 0 1s 幅度特性更近似模拟的 只在折叠频率pi处有轻微混叠 T 1s 频谱混叠严重 脉冲响应不变法优点 数字滤波器的单位脉冲响应h n 完全模仿模拟滤波器的单位冲激响应h t 也就是时域逼近良好 模拟频率 和数字频率 之间呈线性关系 T 因而 一个线性相位的模拟滤波器 例如贝塞尔滤波器 通过脉冲响应不变法得到的仍然是一个线性相位的数字滤波器 脉冲响应不变法缺点 有频率响应的混叠效应 所以脉冲响应不变法只适用于限带的模拟滤波器 低通或带通滤波器 而且高频衰减越快 混叠效应越小 对于高通和带阻滤波器 由于它们在高频部分不衰减 因此会产生混叠现象 第四节用双线性变换法设计IIR数字低通滤波器 一 变换原理 了解决上一节的混叠问题 我们应想办法让S到Z域的映射为单值映射 方法 先采用线性压缩的方法 将模拟角频率 压缩至折叠频率以内 即 若 0 压缩后得到 1 1 T 0 T 找Z与S的关系 一般来说 为使AF与DF的某一频率有对应关系 可引入常数C 变换常数C的选取 常数C用来调节频率间的对应关系 若希望AF与DF在低频处有较为确切的对应关系 可以选择 解释 在低频处有较为确切的对应关系 即要求低频处 当 1较小时 处于低频处 此时有 而 和 1的对应公式为 二 优点 问题及其解决办法 一 优点 解决了混叠问题 二 问题 三 解决 通过频率预畸来解决 利用模拟滤波器设计IIR数字低通滤波器的步骤 确定数字低通滤波器的技术指标 通带截止频率 p 通带衰减 p 阻带截止频率 s 阻带衰减 s 将数字低通滤波器的技术指标转换成模拟低通滤波器的技术指标 如果采用脉冲响应不变法 边界频率的转换关系为 如果采用双线性变换法 边界频率的转换关系为 3 按照模拟低通滤波器的技术指标设计模拟低通滤波器 4 将模拟滤波器Ha s 从s平面转换到z平面 得到数字低通滤波器系统函数H z 例 设计低通数字滤波器 要求在通带内频率低于0 2 rad时 容许幅度误差在1dB以内 在频率0 3 到 之间的阻带衰减大于15dB 指定模拟滤波器采用巴特沃斯低通滤波器 T 1s 分别用脉冲响应不变法和双线性变换法设计滤波器 解 1 用脉冲响应不变法设计数字低通滤波器 数字低通的技术指标为 p 0 2 rad p 1dB s 0 3 rad s 15dB 模拟低通的技术指标为T 1s p 0 2 rad s p 1dB s 0 3 rad s s 15dB 设计巴特沃斯低通滤波器 先计算阶数N及3dB截止频率 c 根据阶数N 6 查表得到归一化传输函数为 由式 得到 c 0 7032rad s 最后 去归一化 将p s c代入Ha p 中 得到实际的传输函数Ha s 用脉冲响应不变法将Ha s 转换成H z 首先将Ha s 进行部分分式分解 然后再转换z esT 2 用双线性变换法设计数字低通滤波器 数字低通技术指标仍为 p 0 2 rad p 1dB s 0 3 rad s 15dB 模拟低通的技术指标为 设计巴特沃斯低通滤波器 阶数N计算如下 取N 6 并求得 c 0 7662rad s 根据N 6 查表6 2 1得到的归一化传输函数Ha p 与脉冲响应不变法得到的相同 为去归