《IIR和FIR数字滤波器的设计及其结构研究文章研究报告》由会员分享，可在线阅读，更多相关《IIR和FIR数字滤波器的设计及其结构研究文章研究报告（16页珍藏版）》请在金锄头文库上搜索。
1、1 目录目录 一 数字滤波器设计的意义 1 二 IIR 数字滤波器 1 2 1 IIR 数字滤波器的设计与方法 1 2 2IIR 数字滤波器的两种设计方法 3 2 2 1 脉冲响应不变法 3 2 2 2 双线性变换法简介 4 2 2 3 脉冲响应不变法与双线性变换法的优缺点比较 6 三 线性相位 FIR 数字滤波器 6 3 1 线性相位 FIR 数字滤波器的基本特性 7 3 2 利用窗函数法设计线性相位 FIR 数字滤波器 7 3 3 利用频率取样法设计线性相位 FIR 数字滤波器 9 3 4 窗函数法与频率取样法的比较 10 四 IIR 和 FIR 数字滤波器的基本结构 10 4 1IIR
2、数字滤波器基本结构 10 4 1 2IIR 系统的直接实现形式 11 4 1 2IIR 系统的级联实现形式 11 4 1 3IIR 系统的并联实现形式 12 4 2FIR 数字滤波器基本结构 12 4 2 1FIR 系统的直接实现形式 12 4 2 2FIR 系统的级联实现形式 12 五 IIR 和 FIR 数字滤波器的比较 13 5 1IIR 数字滤波器的主要优点 13 5 2IIR 数字滤波器的主要缺点 13 5 3FIR 数字滤波器的主要优点 13 5 4FIR 数字滤波器的主要缺点 13 六 数字滤波器比较概括性总结 14 参考文献 14 1 摘要 数字滤波器是指输入 输出均为数字信号
3、 通过数值运算处理改变输入信 号所含频率成分的相对比例 或者滤除某些频率成分的数字器件和程序 经典 数字滤波器从滤波特性上分类 可以分成低通 高通 带通和带阻等滤波器 根据数字滤波器冲激响应的时域特性 可以分成无限脉冲响应数字滤波器 简称 IIR 和有限脉冲响应数字滤波器 简称 FIR IIR 和 FIR 数字滤波器的设计方法 及其结构各不相同 关键字 数字滤波器 IIR FIR 一 数字滤波器设计的意义 在数字信号处理中 FIR 数字滤波器是最常用的单元之一 它用于将输入信 号 x n 的频率特性进行特定的修改 转换成另外的输出序列 y n 与 IIR 滤波 器相比较 在设计和实现上 FIR
4、 滤波器具有如下优越性 1 相位响应可为严格的线性 因此它不存在延迟失真 只有固定的时间延迟 2 由于不存在稳定性问题 所以设计相对简单 3 只包含实数算法 不涉及复数算法 不需要递推运算 另外 也应看到 IIR 滤波器虽然设计简单 但主要是用于设计具有分段 常数特性的滤波器 如低通 高通 带通和带阻等 往往脱离不了模拟滤波器 的格局 而 FIR 滤波器则要灵活的多 尤其是他易于适应某些特殊应用 如构 成数字微分器或希尔伯特变换器等 因而有更大的适应性和广阔的应用领域 二 IIR 数字滤波器 2 1 IIR 数字滤波器的设计与方法 数字滤波器是对数字信号实现滤波的线性时不变系统 数字滤波实质上
5、是一 种运算过程 实现对信号的运算处理 输入数字信号 数字序列 通过特定的 运算转变为输出的数字序列 因此 数字滤波器本质上是一个完成特定运算的 数字计算过程 也可以理解为是一台计算机 描述离散系统输出与输入关系的 1 卷积和差分方程只是给数字信号滤波器提供运算规则 使其按照这个规则完成 对输入数据的处理 时域离散系统的频域特性为 2 1 jejwjw eHeXeY 其中 分别是数字滤波器的输出序列和输入序列的频域特性 jw eY jw eX 或称为频谱特性 是数字滤波器的单位取样响应的频谱 又称为数 je eH 字滤波器的频域响应 输入序列的频谱经过滤波后 因此 jw eX jejw eHe
6、X 只要按照输入信号频谱的特点和处理信号的目的 适当选择 使得滤波 je eH 后的满足设计的要求 这就是数字滤波器的滤波原理 jejw eHeX 数字滤波器根据其冲激响应函数的时域特性 可分为两种 即无限长冲激 响应 IIR 数字滤波器和有限长冲激响应 FIR 数字滤波器 IIR 数字滤波器的 特征是 具有无限持续时间冲激响应 需要用递归模型 来实现 其差分方程 为 2 2 1 1 10 nybnxany N i i N i i 系统函数为 2 3 N k k k M r r r za zb ZH 1 0 1 设计IIR滤波器的任务就是寻求一个物理上可实现的系统函数 使其频 zH 率响应满足
7、所希望得到的频域指标 即符合给定的通带截止频率 阻带截 zH 止频率 通带衰减系数和阻带衰减系数 IIR数字滤波器是一种离散时间系统 其系统函数为 2 4 1 1 0 zx zy ZH N k k k M k k k za zb 假设M N 当M N时 系统函数可以看作一个IIR的子系统和一个 M N 的FIR 子系统的级联 IIR数字滤波器的设计实际上是求解滤波器的系数和 akbk 1 它是数学上的一种逼近问题 即在规定意义上 通常采用最小均方误差准则 去逼近系统的特性 如果在S平面上去逼近 就得到模拟滤波器 如果在z平面 上去逼近 就得到数字滤波器 2 2IIR 数字滤波器的两种设计方法
8、2 2 1 脉冲响应不变法 利用模拟滤波器来设计数字滤波器 也就是使数字滤波器能模仿模拟滤波器 的特性 这种模仿可以从不同的角度出发 脉冲响应不变法是从滤波器的脉冲 响应出发 使数字滤波器的单位脉冲响应序列模仿模拟滤波器的冲激响应 nh 即将进行等间隔采样 使正好等于ha t 的采样值 满足 ths ths nh nThnh s 如果令是的拉普拉斯变换 为的Z变换 利用采样序 sHa tha zH zh 列的Z变换与模拟信号的拉普拉斯变换的关系得 2 5 2 1 X 1 k k a ez st k T jsX T jks T zX as 可看出 脉冲响应不变法将模拟滤波器的 S 平面变换成数字
9、滤波器的 Z 平面 这个从s到z的变换是从 S 平面变换到 Z 平面的标准变换关系式 脉冲 st ez 响应不变法的映射关系如图 1 1 所示 图 1 1 脉冲响应不变法的映射关系 由 2 5 式 数字滤波器的频率响应和模拟滤波器的频率响应间的关系为 2 6 2 H 1 a T kw j T eH k jw 这就是说 数字滤波器的频率响应是模拟滤波器频率响应的周期延拓 正如 采样定理所讨论的 只有当模拟滤波器的频率响应是限带的 且带限于折叠频 率以内时 即 j 3 T T 3 T T o o 11 jIm z Re z Z S 1 2 7 0 jHa 2T s 才能使数字滤波器的频率响应在折叠
10、频率以内重现模拟滤波器的频率响应 而不产生混叠失真 即 2 8 1 T w jH T eH a jw w 但是 任何一个实际的模拟滤波器频率响应都不是严格限带的 变换后就会 产生周期延拓分量的频谱交叠 即产生频率响应的混叠失真 这时数字滤波器 的频响就不同于原模拟滤波器的频响 而带有一定的失真 当模拟滤波器的频 率响应在折叠频率以上处衰减越大 越快时 变换后频率响应混叠失真就越小 这时 采用脉冲响应不变法设计的数字滤波器才能得到良好的效果 脉冲响应 不变法中的频响混叠现象如图 1 2 所示 图 1 2 脉冲响应不变法中的频响混叠现象 对某一模拟滤波器的单位冲激响应ha t 进行采样 采样频率为
11、fs 若使fs 增加 即令采样时间间隔 T 1 fs 减小 则系统频率响应各周期延拓分量之 间相距更远 因而可减小频率响应的混叠效应 2 2 2 双线性变换法简介 脉冲响应不变法的主要缺点是产生频率响应的混叠失真 这是因为从 S 平面 到 平面是多值的映射关系所造成的 为了克服这一缺点 可以采用非线性频 率压缩方法 将整个频率轴上的频率范围压缩到 T T之间 再用z esT 转换到 Z 平面上 也就是说 第一步先将整个 S 平面压缩映射到 S1平面的 T T一条横带里 第二步再通过标准变换关系z es1T将此横带变换到整 3 2 j a H o o 2 3 T ej H T 2 T T T 2
12、 1 个Z平面上去 这样就使 S 平面与Z平面建立了一一对应的单值关系 消除了 多值变换性 也就消除了频谱混叠现象 映射关系如图 1 3 所示 图 1 3 双线性变换的映射关系 为了将 S 平面的整个虚轴 j 压缩到 S1 平面 j 1 轴上的到段T T 上 可以通过以下的正切变换实现 2 2 tan 2 1T T 9 当 1 由经过 0 变化到时 由 经过 0 变化到 也即映T T 射了整个 j 轴 将式 2 9 写成 2 2 1 1 1 1 2 2 2 T j j T j j e T e e T e T j 将此关系解析延拓到整个 S 平面和 S1 平面 令 则得SJ 1SJ 2 10 T
13、s e e e T s sT 1 1 12 1 再将 S1 平面通过以下标准变换关系映射到 Z 平面 1 ts e Z 从而得到 S 平面和 Z 平面的单值映射关系为 2 11 1 1 1 12 z z T S o 11 Z jIm z Re z T j 1 1 T S1 S j oo 1 2 12 s T s T z 2 1 2 1 式 2 11 与式 2 12 是 S 平面与 Z 平面之间的单值映射关系 这种变换 都是两个线性函数之比 因此称为双线性变换 式 2 9 与式 2 10 的双线性变换符合映射变换应满足的两点要求 首先 把 可得 jw ez 2 13 js 即 S 平面的虚轴映射
14、到 Z 平面的单位圆 其次 将代入式 2 13 得 js j T j T z 2 2 因此 2 14 22 22 2 T 2 z T 由此看出 当 0 时 z 0 时 z 1 也就是说 S 平面的左 半平面映射到 Z 平面的单位圆内 S 平面的右半平面映射到 Z 平面的单位圆外 S 平面的虚轴映射到 Z 平面的单位圆上 因此 稳定的模拟滤波器经双线性变 换后所得的数字滤波器也一定是稳定的 2 2 3 脉冲响应不变法与双线性变换法的优缺点比较 从以上讨论可以看出 脉冲响应不变法使得数字滤波器的单位脉冲响应完全 模仿模拟滤波器的单位冲激响应 也就是时域逼近良好 而且模拟频率 和数 字频率 之间呈线
15、性关系 T 最大缺点是有频率响应的混叠效应 双线性变换法与脉冲响应不变法相比 其主要的优点是避免了频率响应的混叠 现象 这是因为 S 平面与 Z 平面是单值的一一对应关系 S 平面整个 j 轴单 值地对应于 Z 平面单位圆一周 即频率轴是单值变换关系 但是双线性变换的这个特点是靠频率的严重非线性关系而得到的 由于这种 频率之间的非线性变换关系 导致一个线性相位的模拟滤波器经双线性变换后 得到非线性相位的数字滤波器 不再保持原有的线性相位了 而且这种非线性 1 关系要求模拟滤波器的幅频响应必须是分段常数型的 即某一频率段的幅频响 应近似等于某一常数 三 线性相位 FIR 数字滤波器 3 1 线性
16、相位 FIR 数字滤波器的基本特性 单位抽样响应 h n 是有限长的 因此 FIR 数字滤波器一定是稳定的 经延 时 h n 总可变成因果序列 所以 FIR 数字滤波器总可以由因果系统实现 h n 为有限长 可以用 FFT 实现 FIR 数字滤波器 FIR 的系统函数是 Z 1 的多 项式 故 IIR 的方法不适用 FIR 的相位特性可以是线性的 因此 它有更广 泛的应用 3 2 利用窗函数法设计线性相位 FIR 数字滤波器 用窗函数设计滤波器首先要对滤波器提出性能指标 一般是给定一个理想的 频率响应 使所设计的 FIR 滤波器的频率响应去逼近所要求的 理想的滤波器的相应 窗函数法设计的任务在于寻找一个可实现 有 限长单位脉冲响应 的传递函数 3 1 0 1 去逼近 我们知道 一个理想的频率响应的傅理叶变换 3 2 1 2 2 0 所得到的理想单位脉冲响应往往是一个无限长序列 对经过适 当的加权 截断处理才得到一个所需要的有限长脉冲响应序列 对应不同的加 权 截断 就有不同的窗函数 所要寻找的滤波器脉冲响应就等于理想脉冲响 应和窗函数的乘积 即 由此可见 窗函数的性质就决定 了滤波器的
