Download
1 / 58
580 likes | 720 Views
IIR DF 的设计方法. 一、滤波器的设计方法. DF 是一种具有 频率选择性 的离散线性系统，即选频滤波器。它在确定信号与随机信号的数字处理中有着广泛的应用。 数字滤波器的设计是确定其 系统函数 并实现的过程。. 二、滤波器设计的步骤. 1.根据任务，确定性能指标。 2.用因果系统的线性时不变系统函数去逼近。 3.用有限精度算法实现这个系统函数。（包括选择运算结构、选择合适的字长、有效数字处理方法。） 4.用适当的软、硬件技术实现 包括采用：通用计算机软件、数字滤波器硬件、或者二者结合。. 三、性能指标. 我们在进行滤波器设计时，需要确定其性能指标。
E N D
一、滤波器的设计方法 • DF是一种具有频率选择性的离散线性系统，即选频滤波器。它在确定信号与随机信号的数字处理中有着广泛的应用。 • 数字滤波器的设计是确定其系统函数并实现的过程。
二、滤波器设计的步骤 1.根据任务，确定性能指标。 2.用因果系统的线性时不变系统函数去逼近。 3.用有限精度算法实现这个系统函数。（包括选择运算结构、选择合适的字长、有效数字处理方法。） 4.用适当的软、硬件技术实现 包括采用：通用计算机软件、数字滤波器硬件、或者二者结合。
三、性能指标 • 我们在进行滤波器设计时，需要确定其性能指标。 • 理想滤波器物理不可实现的。（由于从一个频带到另一个频带之间的突变） • 要物理可实现：应从一个带到另一个带之间设置一个过渡带且在通带和止带内也不应该严格为1或零。应给以较小容限。
1、低通滤波器的性能指标 |H(ejw)|或|H(f) δ1:通带的容限 1 δ2：阻带容限 Ap 1-δ1 通带截止频率：fp(wp)又称为通带上限频率。 通带衰减：Ap 阻带截止频率：fs(ws)又称阻带下限截止频率。 阻带衰减：As δ2 As fs ws fp wp f w
2、高通滤波器的性能指标 |H(ejw)|或|H(f) Ap 通带截止频率：fp(wp)又称为通带下限频率。 通带衰减：Ap 阻带截止频率：fs(ws)又称阻带上限截止频率。 阻带衰减：As 1 As fs ws fp wp f w
3、带通滤波器的性能指标 |H(ejw)|或|H(f) Ap 通带截止频率：上限截止频率fp2(wp2)，下限截止频率fp1(wp1)。 通带衰减：Ap 阻带截止频率：上限截止频率fs2(ws2)，下限截止频率fs1(ws1)。 阻带衰减：As 1 As f w fs1 ws1 fp1 wp1 fs2 ws2 fp2 wp2
4、带阻滤波器的性能指标 |H(ejw)|或|H(f) Ap 通带截止频率：上限截止频率fp2(wp2)，下限截止频率fp1(wp1)。 通带衰减：Ap 阻带截止频率：上限截止频率fs2(ws2)，下限截止频率fs1(ws1)。 阻带衰减：As 1 As fs1 ws1 f w fp1 wp1 fs2 ws2 fp2 wp2
四、H(z)如何推导出 • 有时根据提出对滤波器的性能要求、频率特性（低、高、带通、带阻）来设计系统H(z). • 有时根据时域波形提出要求来设计-->单位冲激响应h(n)的形状。 • 有时也直接给出H(z).(但要求因果稳定).
五、确定DF的采用的结构及运算结构的好坏 • 确定DF的采用的结构将会影响DF的精度、稳定性、经济性及运算速度等很多重要性质。 • 1.计算复杂性 一个运算结构应含有最少的乘法器和最少的延时器。乘法器最费时间， 延时器最费存储单元。 2. 存储器长度的有限和运算结构将会有影响。即有时会希望使用一种运算结构，虽然它的乘法器和延时器并不是最少的，但它对存储器的有限字长效应是最不敏感的。
三、IIR DF 的设计方法 • 设 计 IIR 数 字 滤 波 器 系 统 函 数 方 法： • 1、简单滤波器的零、极点累试法 • 2、间 接 方 法 • 3、直 接 方 法
1、简单滤波器的零、极点累试法 在z平面上直接设计IIR数字滤波器，以滤波器响应作为依据，直接在z平面上，通过多次选定极点和零点位置逼近该响应。 即在单位园内设置一对共轭极点 ，频响在w0处就有一峰值。 r越近于1，极点位置越接近单位园，则峰值就越尖锐。同理，若在单位园内设置一对零点 ，频响就会在w1处出现零值，即可实现陷波。如特性还达不到要求，可再移动零、极点，这样作二、三次调整后，就可以获得一些简单的DF.这种方法，可以设计一些简单阶数很低（1~2阶)的DF。 Im[z] * Re[z] *
例子 • 从图中清楚表明，极点愈靠近单位园（r愈靠近1），带通特性愈尖锐。
2、间 接 方 法 • 由于模拟滤波器设计技术是非常成熟的，归一化各种模拟低通滤波器的系统函数已有表可查，利用成熟的设计技术，可得到一个间接设计IIR DF的方法，即间接设计方法。 • 这 种 方 法 通 常 要 先 设 计 一 中 间 滤 波 器 ， 然 后 通 过 映 射 或 频 率 变 换 完 成 最 终 IIR 数 字 滤 波 器 的 设 计。 这 种 间 接 设 计 方 法中 包 括： (1) 由模 拟滤波器设计数字滤波器 (2) 频 率 变 换 法（分为模拟频率变换法和数字频率变换法）来设计数字滤波器
3、直 接 方 法 • 直 接 方 法( 计 算 机 辅 助 设 计 法 ) （1）在频域利用幅度平方误差最小法直接设计IIR数字滤波器。 （2）在时域直接设计IIR数字滤波器 此 法 根 据 性 能 指 标 和 一 定 的 逼 近 准 则， 直 接 利 用 计 算 机 完 成 设 计。
一、由模拟滤波器设计数字滤波器步骤 • 这 种 方 法 设 计 经 以 下 三 步： • 1、数字滤波器的技术指标转换成模拟低通滤波器指标 • 2、模拟低通滤波器设计 • 3、映射实现：从模拟低通滤波器再转换成数字滤波器
1、数字滤波器的技术指标转换成模拟滤波器指标1、数字滤波器的技术指标转换成模拟滤波器指标 • 根 据 给 定 设 计 要 求， 把 数 字 滤 波 器 的 性 能 指 标 变 成 模 拟 滤 波 器低通 的 性 能 指 标。 • 这是因为模拟滤波器设计方法已经很成熟，它不仅有完整的设计公式，还有完善的图表供查阅；另外，还有一些典型的滤波器类型可供使用。（采用冲激不变法和双线性变换法）
2、模拟滤波器设计 • 设 计 出 符 合 要 求 的 模 拟 滤 波 器 的 系 统 函 数Ha(s)。可以选择多种类型的滤波器。如Butterworth,Chebyshev,Ellipse,Bessel等。
3、 映 射 实 现 • 利 用 一 定 的 映 射 方 法， 把 模 拟 滤 波 器 系 统 函 数 数 字 化， 完 成 IIR 数 字 滤 波 器 系 统 函 数 的 设 计。（采用冲激不变法和双线性变换法）
三、由模拟滤波器设计数字滤波器这种方法的局限性三、由模拟滤波器设计数字滤波器这种方法的局限性 • 用 这 方 法 设 计IIR 数 字 滤 波 器 存 在 一 定 的 局 限 性： • (1) 这 种 方 法 只 适 于 设 计 幅度 响 应 分 段 恒 定 的 滤 波 器， 不 能 解 决 多 带 或 任 意 幅 度 的 滤 波 器 的 设 计 问 题。 • (2) 它 只 是 一 种 可 行 的 方 法， 而 不 是 最 优 的 方 法。
一、设计模拟低通滤波器 • 首先将要设计的数字滤波器的指标，转变成模拟低通原型滤波器的指标后，设计“模拟低通原型”滤波器。 • 在此介绍常用的Butterworth巴特渥斯滤波器模拟低通滤波器的设计。 • 它有严格的设计公式，现成的曲线和图表供设计。
典型模拟滤波器的特点 1、Butterworth巴特渥斯滤波器 它具有单调下降的幅频特性；即最平幅度。 2、Chebyshev切比雪夫滤波器 在通带或阻带有波纹，可提高选择性。 3.Bessel贝塞尔滤波器 在通带内有较好的线性相位特性。 4.Ellipse椭圆滤波器 其选择性相对前三种是最好的。
二、模拟滤波器设计思想 • 根据模拟滤波器设计要求，求出相应的模拟系统函数. • 使其逼近某个理想滤波器的特性。(滤波器的特性包括有：幅度特性、相位特性/群时延特性)，在此我们采用幅度平方函数特性来设计。
三、根据幅度平方函数确定系统函数1、求滤波器的幅度平方函数三、根据幅度平方函数确定系统函数1、求滤波器的幅度平方函数 • 设 计 模 拟 滤 波 器 经 常 要 借 助 其 幅 度 平 方 函 数 其中:Ha(s)是模拟滤波器的系统函数。 • 假 设p1, z1为Ha(s) 的一个零点和一个极点，则-p1, -z1必为Ha(-s)的一个零点和极点，Ha(s)、Ha(-s)的零极点成象限对称分布。所以必然有如下形式： z1 p1 -z1 -p1 ＊ ＊ ＊ ＊
2、根据幅度平方函数设计模拟滤波器的系统函数的步骤2、根据幅度平方函数设计模拟滤波器的系统函数的步骤 • 我 们 知 道， 实 际 滤 波 器 都 是 稳 定 的， 因 此 其 极 点 一 定 位 于S 平 面 左 半 平 面， 这 样 可 根 据 幅 度 平 方 函 数 通 过 如 下 步 骤 分 配 零、 极 点 来 设 计 出 模 拟 滤 波 器 的 系 统 函 数 。 • (1)由 来 确 定 象 限 对 称的 S 平 面 函 数。 • (2)将 因 式 分 解， 得 到 各 零 点 和 极 点。 • (3) 按 照 与Ha(s) 的 低 频 特 性 或 高 频 特 性 的 对 比 就 可 确 定 出 增 益 常 数。
(1)由 来 确 定 象 限 对 称的 S 平 面 函 数。 • 将 • 代入 中即得到s平面函数。
(2)将 因 式 分 解， 得 到 各 零 点 和 极 点。 • 将左半平面的极点归于Ha(s)。 • 如无特殊要求，可取 的对称零点的任一半作为 Ha(s) 的零点。 • 如要求是最小相位延时滤波器，则应取左半平面零点作为Ha(s) 的零点。且 轴上的零点或极点都是偶次的，其中一半属于Ha(s)。
(3) 按 照 与Ha(s) 的 低 频 特 性 或 高 频 特 性 , 确 定 出 增 益 常 数。 • 由 的条件，代入可求得增益常数。
例子 • 根据以下幅度平方函数 确定系统函数Ha(s).
Butterworth巴特渥斯低通滤波器1、幅度平方函数 • Butterworth 低 通 滤 波 器 具 有 通 带 最 平 幅 度逼近特 性， 是 一 全 极 点 型 滤 波 器，且极点均匀分布上Ωc的园上，并且与虚轴对称。 • 其最主要特点：在通带内，幅频最平坦，随着频率的升高而单调下降。其幅 度 平 方 函 数为 • 其 中N 为 整 数， 表 示 滤 波 器 的 阶 次，Ωc 定 义 为 截 止 频 率， 为 振 幅 响 应 衰 减 到 - 3dB 处 的 频 率。
2、Butterworth滤波器的极点分布 由 可知Butterworth的零点全部在S=∞处，它是全极点型滤波器，且分布在半径为Ωc 的圆上,呈象限对称分布。 为了得到稳定的滤波器，s左半平面的极点必须分配给Ha(s)，s右半平面的极点分配给Ha(-s)。 取其分布在左平面的极点, 设计出巴特沃斯低通滤波器.
3、Butterworth的幅度响应及极点分布 其中左半平面构成Butterworth滤波器的系统函数 极点不会落在S平面上的虚轴上
4、Butterworth滤波器阶数N与幅度响应的关系 当N增大时，滤波器的特性曲线变得陡峭，则更接近理想矩形幅度特性。
6、Butterworth滤波器的特点（总结） (1)当Ω=0时， 即Ω=0处无衰减。 (2)当Ω=Ωc时， 或 通带最大衰减 在止带内的逼近是单调变化的，不管N为多少，所有幅频特性曲线都经过-3dB点，或说衰减3dB,这就是3dB不变性。
(3) 在Ω<Ωc的通带内： 有最大平坦的幅度特性，即N阶Butterworth低通滤波器在Ω=0处： 前(2N-1)阶导数为零，因而Butterworth又称最平幅度特性滤波器。 随着Ω由0变到Ωc ，|Ha(jΩ)|2单调减小，N越大，减小越慢，也就是通带内特性越平坦。
(4) 在Ω>Ωc，即在过渡带及阻带中，|Ha(jΩ)|2 也随Ω增加而单调减小，但是Ω/Ωc >1，故比通带内衰减的速度要快得多，N越大，衰减速度越大。 当Ω=Ωs，即频率为阻带截止频率时，衰减为： 2为阻带最小衰减。 (5)滤波器的特性完全由其阶数N决定。N越大，则通带内在更大范围内更接近于1，在止带内迅速地接近于零，因而幅频特性更接近于理想的矩形频率特性。
7、归一化的Butterworth滤波器的系统函数 设实际所给的频率为 ，归一化后的频率为，对低通模拟滤波器，令 所以， , ，并令归一化复数变量为p, 则有 定义一个衰减函数 ，所以， 使滤波器的指标与幅频特性联系。
8、Butterworth滤波器设计步骤 • 1，将实际频率归一化，的归一化幅平方特性 • 2，求C和N • 则 • 即， • 解得： （1）
由于 当C=1时， 由得 k=1,2,,N （2） 所以 （3）
9、例子 试设计一模拟低通巴特沃思滤波器，要求截止频率 ， 通带最大衰减 ， 阻带起始频率 ， 阻带最小衰减 。
首先，将频率归一化 ， ， ， 因为 ，所以C=1 由（1）有 取N=5
由（2），有5个极点分别是： 再由（3），可求出
总的归一化转移函数是 的级连，即 最后，
四、映射实现的方法 • 由模拟滤波器映射成数字滤波器的方法，也即，数字滤波器能模仿滤波器的特性。 主要有以下几种映射方法： • 冲 激 响 应 不 变 法 • 阶 跃 响 应 不 变 法 • 双 线 性 变 换 法
冲激不变法（脉冲响应不变法） 由雷道（Rader),戈尔登(Golden)提出
一、变 换 原 理1、什么是冲激不变法 • 冲 激 响 应 不 变 法 是 从 时 域 出 发， 要 求 数 字 滤 波 器 的冲 激 响 应 h(n) 对 应 于 模 拟 滤 波 器ha(t) 的 等 间 隔 抽样。 h(n)=ha(nT) ， 其 中 T 是 抽 样 周 期。 因 此 时 域 逼 近 良 好。
