《dddIIR数字滤波器》PPT课件.ppt

数字滤波器 DF 的定义输入和输出均是数字信号 通过一定运算关系 数值运算 改变输入数字信号所含频率成份的相对比例或滤除某些频率成份的器件 如何用数字滤波器处理模拟信号 通过A DC和D AC 用数字滤波器对模拟信号进行处理 第6章无限脉冲响应数字滤波器的设计 数字滤波器的数学描述 差分方程 系统函数 1 数字滤波器的分类 1 一般分类经典滤波器 选频滤波器 适用于信号中有用的频率成分和希望滤除的频率成分占用不同的频带的情况 通过选频滤波器达到滤波的目的 现代滤波器 适用于信号和干扰的频带相互重叠的情况 通过现代滤波器利用信号的统计分布规律 从干扰中最佳提取信号 如 维纳滤波器 卡尔曼滤波器和自适应滤波器等 本书只讨论选频滤波器的设计 2 从滤波器的功能上来分类分为低通 高通 带通 带阻滤波器 理想滤波器幅度特性 特点 1 h n 非因果且无限长 物理不可实现 只能尽可能逼近 2 DF的频率响应以2 为周期 低频区域处于 的偶数倍附近 高频区域处于 的奇数倍附近 3 从滤波器的实现网络结构或从单位脉冲响应来分类无限脉冲响应 IIR 数字滤波器有限脉冲响应 FIR 数字滤波器 在数字滤波器中 一般考察其半个周期 0 的频域特性 在模拟滤波器中 通常考察其 0 范围内频率域的特性 N阶IIR滤波器系统函数 N 1阶FIR滤波器系统函数 2 数字滤波器的技术要求数字滤波器的传输函数H ejw H ejw 系统的幅频特性 表示信号通过该滤波器后各频率成分衰减情况 系统的相频特性 反映各频率成分通过滤波器后在时间上的延时情况 p 通带截止频率 通带频率范围 0 p S 阻带截止频率 阻带频率范围 s C 3dB截止频率 P 通带最大衰减 S 阻带最小衰减 1 通带内幅度响应误差范围 2 阻带内幅度响应误差范围 如将 H ej0 归一化为1 上两式则表示成 3 数字滤波器设计方法IIR滤波器设计方法 1 先设计模拟滤波器 AF 的系统函数Ha s 然后按某种变换 将Ha s 转换成数字滤波器的系统函数H z 2 借助CAD 计算机辅助设计 在频域或时域直接设计 FIR滤波器设计方法 1 经常采用的是窗函数设计法和频率采样法 2 用计算机辅助的切比雪夫最佳一致逼近法设计 6 3模拟滤波器 AF 的设计 模拟滤波器的理论和设计方法已发展得相当成熟 且有若干典型的模拟滤波器可以选择 如 巴特沃斯 Butterworth 滤波器 切比雪夫 Chebyshev 滤波器 椭圆 Kllipse 滤波器 贝塞尔 Bessel 滤波器等 这些滤波器都有严格的设计公式 现成的曲线和图表供设计人员使用 AF的表示 AF的设计 找到ha t 或Ha s 更多时候 即找到典型滤波器公式中的阶和未知参数 各种理想模拟滤波器的幅度特性 可实现模拟滤波器的幅度特性 滤波器 无源滤波器 由R L C等元件组成 有源滤波器 由运放 R C等元件组成 有源滤波器的工作频率难以做得很高 而无源滤波器的工作频率很高 设为某一频率正弦量 有源低通滤波器 为改善滤波效果 常采用二阶低通滤波器 有源高通滤波器 只需将低通滤波器的R C互换 1 模拟低通滤波器的设计指标及逼近方法 1 设计指标有 p s p s其中 p和 s分别称为通带截止频率和阻带截止频率 p是通带 0 p 内允许的最大衰减系数 s是阻带 s 内允许的最小衰减系数 图中 c称为3dB截止频率 因 如果 0处幅度已归一化到1 即 Ha j0 1 技术指标高指 p越小 s越大 过渡带越小越好 技术指标可归纳为 2 用模拟滤波器逼近方法设计滤波器步骤 给出模拟滤波器的技术指标 设计传输函数Ha s 使其幅度平方函数满足给定指标ap和as Ha j 2 Ha j Ha j Ha s Ha s S j 确定Ha s 系统Ha s 应是稳定的系统 因此 极点应位于S左半平面内 特点 单调变化 当 0时 H j 1 3dB点必通过 当 c时 H j N越大 曲线越陡峭 当N 就是理想滤波器 2 Butterworth低通滤波器的设计方法巴特沃斯低通滤波器的模方函数 Ha j 2用下式表示 N 滤波器阶数 c是3dB截止频率 Butterworth低通滤波器的幅度函数只由阶数N控制 巴特沃斯滤波器的设计步骤 1 根据给出的技术指标 P S p S 求N c N可能有小数部分 取大于等于N的最小整数 k 0 1 2N 1 2N个极点等间隔分布在半径为 c的圆上 间隔是 Nrad 其中 一半是Ha s 的极点 另一半是Ha s 的极点 则左半平面的N个极点属于Ha s 稳定系统 右半平面N个极点属于Ha s N 3 设N 3 极点有6个 它们分别为 取s平面左半平面的极点s0 s1 s2组成Ha s 极点在s平面呈象限对称 分布在Buttterworth圆上 共2N点 极点间的角度间隔为 极点不落在虚轴上 N为奇数 实轴上有极点 N为偶数 实轴上无极点 一半极点在左半平面 一半极点在右半平面 c 称为归一化频率 可知 只要知道阶数N 就可得到极点pk 查表P120表6 1可直接得到极点pk和归一化系统函数Ha p 图表法 为使滤波器能统一设计 将所有的频率对3dB截止频率 c归一化 归一化后的Ha s 表示为 将Ha p 去归一化 得到系统函数Ha s 表6 1巴特沃斯归一化低通滤波器参数 表示两极点P1 PN 2 阶数N的确定 阶数N的大小主要影响幅度特性下降的速度 它应该由技术指标确定 将 p代入幅度平方函数中 将 s代入幅度平方函数中 用上式求出的N可能有小数部分 应取大于等于N的最小整数 关于3dB截止频率 c 如果技术指标中没有给出 可以按照下面两式求出 通常是用一个算出 c 然后用另一个 反过来 来检验 总结以上 低通巴特沃斯滤波器的设计步骤如下 根据技术指标 p p s s 求出滤波器的阶数N 2 求出归一化极点pk 得到归一化传输函数Ha p 3 将Ha p 去归一化 将p s c代入Ha p 得到实际的滤波器传输函数Ha s 解 1 设计模拟滤波器的指标为 p 2 fp 104 rad s ap 2dB s 2 fs 2 4 104 rad s as 30dB 2 确定阶数N和 c 例 通带截止频率fp 5kHz 最大衰减 p 2dB 阻带截止频率fs 12kHz 最小衰减 s 30dB 设计巴特沃斯低通滤波器 3 由N 查表确定Ha p b0 1 0000 b1 3 2361 b2 5 2361 b3 5 2361 b4 3 2361 N 5 4 将p s c代入Ha p 去归一化 得Ha s 3 为将Ha p 去归一化 先求3dB截止频率 c 检验 可以看出 满足 s 30dB的真实fs在10 525kHz处 与12kHz比 还有富裕量 6 3 4Chebyshev低通滤波器的设计方法 提出的背景巴特沃斯滤波器的频率特性曲线 无论在通带和阻带都是频率的单调函数 因此 当通带边界处满足指标要求时 通带内肯定会有余量 因此 更有效的设计方法应该是将精确度均匀地分布在整个通带内 或者均匀分布在整个阻带内 或者同时分布在两者之内 这样 就可用阶数较低的系统满足要求 这可通过选择具有等波纹特性的逼近函数来达到 切比雪夫多项式 当N 0时 C0 x 1当N 1时 C1 x x当N 2时 C2 x 2x2 1当N 3时 C3 x 4x3 3x当N 4时 C3 x 8x4 8x2 1切比雪夫多项式的递推公式为CN 1 x 2xCN x CN 1 x 切比雪夫 型与巴特沃斯低通的幅度函数平方曲线 通带内 在1和间等波纹起伏 N为奇数 N为偶数 通带外 迅速单调下降趋向0 Chebyshev滤波器的三个参量 Wp 给定 e N 滤波器阶数N的确定 阻带衰减越大所需阶数越高 滤波器的设计步骤 归一化 2 根据技术指标求出滤波器阶数N及 计算传函 极点法 计算传函 查表法 2 通带波纹1dB 6 4用模拟滤波器设计IIR数字滤波器 利用模拟滤波器设计IIR数字滤波器设计过程 设计技术成熟有相当简便的公式和图表 要求DF特性模仿AF的特性 实际上是个映射问题Mapping 若转换后的H z 稳定且满足技术要求 则设计过程要满足 1 因果稳定的模拟滤波器转换成数字滤波器 仍是因果稳定的 2 数字滤波器的频率响应模仿模拟滤波器的频响 s平面的虚轴映射到z平面的单位圆上 满足上述转换关系的映射方法有 脉冲响应不变法和双线性变换法 一 脉冲响应不变法的基本思想使h n 模仿ha t 让h n 正好等于ha t 的采样值二 变换方法设模拟滤波器Ha s 只有单阶极点 且分母多项式的阶次高于分子多项式的阶次 将Ha s 用部分分式表示 Ha s LT 1 Ha s ha t 时域采样 h n ZT h n H z 所以说脉冲响应不变法是一种时域上的变换方法 si为Ha s 的单阶极点 LT 1 Ha s U t 为单位阶跃函数 对ha t 进行等间隔采样 采样间隔为T 得到 对上式进行Z变换 得到数字滤波器的系统函数H z 结论 1 S平面的单极点s si映射到Z平面的极点z esiT 2 Ha s 部分分式的系数与H z 部分分式的系数相同 3 左半平面映射到单位圆内 模拟滤波器稳定则数字滤波器稳定 三 S平面和Z平面之间的映射关系1 采样信号的拉氏变换与相应的序列的Z变换之间的映射关系 1 设ha t 的采样信号表示为 2 对进行拉氏变换 得到 2 模拟信号的拉氏变换与相应的序列的Z变换之间的映射关系 模拟信号ha t 的傅里叶变换Ha j 和其采样信号的傅里叶变换之间的关系满足 结论 采样信号的拉氏变换是原模拟信号的拉氏变换在S平面沿虚轴以 s 2 T为周期进行的周期延拓 将s j 代入上式 得 说明 采用脉冲响应不变法将AF变换为数字DF时Ha s 沿虚轴以 s 2 T为周期进行周期延拓 再经过Z eST的映射关系映射到Z平面上 从而得到H z 上面的这种转换是否满足对转换关系提出的2点要求 设 S j Z r ejw脉冲响应不变法标准映射关系 Z eSTr ejw e j T e T ej T得到 r e T w T 频率域的坐标变换是线性的 因果稳定模仿频响 由采样信号的拉氏变换与相应的序列的Z变换之间的映射关系 0时 S平面的右半平面映射到Z平面的单位圆外 r z 1 以上分析结论 若Ha s 是因果稳定的 则转换后的H z 也是因果稳定的 r e T 因果稳定的分析 当 不变 模拟角频率 变化2 T整数倍 映射值不变 S平面上每一条宽度为2 T的水平横带都重迭地映射到Z平面的整个全平面上每条水平横带的左半部分映射到Z平面单位圆内 水平横带右半部分映射到Z平面的单位圆外j 虚轴上每2 T段都对应着单位圆一周 z esT是周期函数 由上面分析结果 S平面与Z平面的映射关系满足转换条件 但存在着多值 s 单值 z 映射关系 3 DF的频响是AF频响的周期延拓 频率混叠 只有AF频响限于 T之内 DF频响才不失真地复现AF频响 否则 设计出来的DF在w 附近产生频率混叠 2 采样信号的拉氏变换是其模拟信号的拉氏变换以2 T为周期 沿虚轴进行周期化 脉冲响应不变法的应用受限只适合设计带限滤波器 如 低通 带通滤波器的设计 不适合高通 带阻滤波器的设计 假设没有频率混叠现象 即满足 按照上式 并将关系式s j 代入 T 代入得到 数字滤波器的频响可以很好模仿模拟滤波器的频响 四 脉冲响应不变法的优缺点优点 1 频率变换是线性关系 w T 模数字滤波器可以很好重现模拟滤波器的频响特性 2 数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应 时域特性逼近好 缺点 1 有频谱混迭失真现象 S平面到 Z平面有多值映射关系 2 由于频谱混迭 使应用受到限制 T 失真 但运算量 实现困难 冲激响应不变法将模拟滤波器转换为数字滤波器 采样间隔为0 1s S平面极点 Z平面极点 例 已知模拟低通滤波器的系统函数为 用脉冲响应不变法将转换成数字滤波器的系统函数 解 首先将展为部分分式的形式 极点为 数字滤波器的系统函数为 设T 1s 得 设T 0 1s 得 将三者的幅度特性用它们的最大值归一化后的幅度特性如图所示 AF的幅度特性 DF的幅度特性 例 已知模拟滤波器的传输函数Ha s 为用脉冲响应不变法将Ha s 转换成数字滤波器的系统函数H z 解 首先将Ha s 写成部分分式 极点为 根据 H z 的极点为 按照 经过整理 得到当 T 1s时用H1 z 表示 T 0 1s时用H2 z 表示 则 将Ha j H1 ejw H2 ejw 的幅度特性用它们最大值归一化后 得到它们的幅度特性曲线 如下图所示 很轻的混叠现象 严重的混叠现象 例 用脉冲响应不变法设计巴特沃思数字低通滤波器 要求在频率小于的通带内 幅度特性下降小于等于1dB 在频率到之间的阻带衰减大于等于15dB 解 1 待设计的数字低通滤波器的技术指标为 查表得到归一化系统函数为 求3dB截止频率 去归一化 4 求相应的数字滤波器的系统函数 5 检查设计的数字滤波器是否满足给定的指标要求 脉冲响应不变法的主要缺点 产生频率谱混迭现象 原因 模拟低通的最高频率超过了折叠频率 T 数字化后在w 形成频谱混叠现象 解决方法 采用非线性压缩方法 将整个频率轴上的频率范围压缩到 T之间 而后再用Z eST转换到Z平面上 6 4 4用双线性变换法设计IIR数字低通滤波器 一 双线性变换法消除频谱混迭的原理1 非线性压缩 S平面 S1平面映射 双线性变换法用正切变换实现非线性频率压缩 设Ha s s j 经过非线性频率压缩后用Ha s1 s1 j 1表示 则 上式表明 当 1从 T经过0变化到 T时 则由 经过0变化到 这样实现了s平面上整个虚轴完全压缩到s1平面上虚轴的 T之间的转换 T 时域采样间隔 由上面可得 1 T T 2 S1平面到Z平面的映射 将S1平面映射到Z平面上 用标准映射Z eS1T 代入上式 双曲正切 th x sh x ch x 3 双线性不变法的映射关系映射过程 从s平面映射到s1平面 再从s1平面映射到z平面 S平面与Z平面是一一对应的单值映射关系 消除了脉冲响应不变法的多值映射关系 消除了频谱混迭现象 只要模拟滤波器Ha s 因果稳定 其极点应位于S左半平面 转换成的H z 也是因果稳定的 位于单位圆内 0 0 1 j 1 j T T 2 Ttan 1T Z es1T S平面 S1平面 Re z jIm z 0 z平面 二 模拟频率 和数字频率 之间的关系将 Z ej S j 代入S Z平面映射关系式 说明 s平面上 与z平面的 成非线性正切关系 当 增加时 增加得很快 当 趋于 时 趋于 由于这种非线性关系 消除了频率混叠现象 代价 影响数字滤波器频响逼真模拟滤波器的频响的逼真度 存在幅度失真和相位失真 频率关系 非线性 三 双线性变换法特点1 优点消除了频谱混迭失真 频率映射表明S平面与Z平面是一一对应的单值映射关系 避免了脉冲响应不变法的频谱 混迭 现象 2 缺点以频率变换的非线性为代价 模拟域和数字域进行非线性映射 其瞬时响应不如脉冲响应不变法好 四 双线性变换法的幅度失真和相位失真情况如果 的刻度是均匀的 通过非线性正切关系 映射到z平面的 刻度不均匀 随 增加越来越密 即边界频率发生畸变 如果模拟滤波器具有片段常数特性 则转换到z平面仍具有片段常数特性 适于片段常系数滤波器的设计 幅度特性失真 相位特性失真 例 已知模拟滤波器的系统函数 试用双线性变换法将转换成数字滤波器的系统函数 设T 1s 则 解 双线性变换法设计的数字滤波器的频响限制在0和 之间 不存在频率响应的混叠现象 例 已知Ha s a a s a 1 RC 求H z 解 1 用脉冲响应不变法时 先确定极点 s a 则2 用双线性不变法时 T为采样间隔 其中 设 a 1000 T 0 001和T 0 002 H1 z 和H2 z 的归一化幅频特性如图所示 双线性变换法 脉冲响应不变法 模拟滤波器幅频特性 250Hz 500Hz 存在混叠 不存在混叠 曲线形状差别大 五 用模拟滤波器设计IIR数字低通滤波器步骤 1 确定数字低通滤波器的性能指标 Wp ap Ws aS 2 确定相应模拟低通滤波器的性能指标 脉冲响应不变法 W T双线性变换法 2 T tan W 3 设计模拟低通滤波器的传输函数Ha S 4 将Ha S 转换成数字滤波器系统函数H z 利用脉冲响应不变法设计时 应先确定Ha s 极点 按Z eST映射极点 求出H z 若利用双线性变换法 设采样间隔为T 则 例6 5 要求设计一数字低通滤波器 指标为 在通带内频率低于0 2 rad s时 容许幅度误差在1dB以内 在频率0 3 之间 阻带衰减大于15dB 指定模拟滤波器采用巴特沃斯低通滤波器 分别用脉冲响应不变法和双线性变换法设计H z 解 采用脉冲响应不变法 1 数字域指标 P 0 2 rad p 1dB S 0 3 rad S 15dB 2 模拟域指标 取T 1 则 p Wp T 0 2 rad s s 0 3 rad s P 1dB s 15dB 3 设计模拟低通AF 要求采用巴特沃斯低通滤滤器 取N 6 查表 得归一化传输函数Ha p 去归一化 c p 100 1ap 1 1 2N 0 7032rad s将p s c代入Ha p 得到去归一化传输函数 4 求DF的H z 求出Ha S 的极点Si由映射Z esit确定出相应数字滤波器的传输函数H z 用脉冲响应不变法设计的数字低通滤波器的幅度特性 用双线性变换法设计数字低通IIR滤波器 1 数字域指标 P 0 2 rad p 1dB S 0 3 rad S 15dB 2 模拟域指标 取T 1 P 1dB s 15dB 3 设计模拟低通AF 要求采用巴特沃斯低通滤滤器 取N 6 查表 得归一化传输函数Ha p 去归一化 c p 100 1ap 1 1 2N 0 7662rad s去归一化传输函数Ha s Ha p p s c 4 求数字Filter的H z 例 用双线性变换法设计数字低通滤波器 要求在频率小于等于100Hz的通带内 幅度的衰减特性不大于2dB 在频率大于等于300Hz的阻带内 衰减不小于15dB 采样频率为1000Hz 采用巴特沃思型滤波器 采样频率为1000Hz 所以T 0 001s 待设计的数字低通滤波器的技术指标为 解 1 确定待设计的数字低通滤波器的技术指标 2 确定模拟低通滤波器的技术指标 设T 1s 3 设计模拟巴特沃思低通滤波器 取N 2 4 由求其相应数字滤波器的系统函数 5 校核所设计的数字滤波器是否满足给定的指标要求 所设计的数字滤波器满足给定的指标要求 6 5模拟高通 带通 带阻滤波器的设计 1 高通的设计 为了简单 通带指标通常为 归一化以为基准 模拟高通滤波器的设计步骤如下 A 确定高通滤波器的技术指标 通带下限频率 p 阻带上限频率 s 通带最大衰减 p 阻带最小衰减 s B 确定相应低通滤波器的设计指标 将高通滤波器的边界频率转换成低通滤波器的边界频率 各项设计指标为 低通滤波器通带截止频率 p 1 p 低通滤波器阻带截止频率 s 1 s 通带最大衰减仍为 p 阻带最小衰减仍为 s C 设计归一化低通滤波器G p D 求模拟高通的H s 将G p 先转换成归一化高通H q 再去归一化 得H s 补充 高通的设计步骤 1 确定归一化的高通滤波器技术指标 2 频率转换成归一化的低通技术指标 3 由低通指标设计归一化的巴特沃斯低通滤波器 G p 4 频率反变换为高通H q 再去归一化求模拟高通的H s 失效 例 设计高通滤波器 fp 200Hz fs 100Hz 幅度特性单调下降 fp处最大衰减为3dB 阻带最小衰减 s 15dB 解 高通技术要求 fp 200Hz p 3dBfs 100Hz s 15dB 归一化频率 低通技术要求 设计归一化低通G p 采用巴特沃斯滤波器 故 求模拟高通H s 2 模拟带通的设计 归一化基准 表6 2 2 与 的对应关系 频率变换关系 带通的设计步骤 1 确定带通滤波器的归一化技术指标 2 频率转换成归一化低通指标 归一化基准 频率变换关系 s与 s的绝对值可能不相等 一般取绝对值小的 s 这样保证在较大的 s处更能满足要求 3 设计归一化低通G p 求N G p 失效 4 频率反变换为带通H q 再去归一化求模拟带通的H s 2 低通到带通的频率变换 低通与带通滤波器的幅度特性如图 频率变换公式 上式称为低通到带通的频率变换公式 利用该式将带通的边界频率转换成低通的边界频率 下面推导由归一化低通到带通的转换公式 将q j 代入上式 得到 为去归一化 将q s B代入 得到 模拟带通的设计步骤 A 确定模拟带通滤波器的技术指标 即 带通上限频率 u 带通下限频率 l 下阻带上限频率 s1 上阻带下限频率 s2 通带中心频率 02 l u 通带宽度B u l 与以上边界频率对应的归一化边界频率如下 B 确定归一化低通技术要求 s与 s的绝对值可能不相等 一般取绝对值小的 s 这样保证在较大的 s处更能满足要求 通带最大衰减仍为 p 阻带最小衰减亦为 s C 设计归一化低通G p D 将G p 转换成带通H s 例 设计模拟带通滤波器 通带带宽B 2 200rad s 中心频率 0 2 1000rad s 通带内最大衰减 p 3dB 阻带 s1 2 830rad s s2 2 1200rad s 阻带最小衰减 s 15dB 解 1 模拟带通的技术要求 0 2 1000rad s p 3dB s1 2 830rad s s2 2 1200rad s s 15dBB 2 200rad s 0 0 B 5 s1 s1 B 4 15 s2 s2 B 6 2 模拟归一化低通技术要求 取 s 1 833 p 3dB s 15dB 3 设计模拟归一化低通滤波器G p 取N 3 查表得 4 求模拟带通H s 模拟带阻的设计 归一化基准 带阻的设计步骤 1 确定带阻滤波器的归一化技术指标 2 频率转换成归一化低通指标 归一化基准 频率变换关系 s与 s的绝对值可能不相等 一般取绝对值小的 s 这样保证在较大的 s处更能满足要求 3 设计归一化低通G p 求N G p 失效 4 频率反变换为带阻H q 再去归一化求模拟带阻的H s 3 低通到带阻的频率变换 低通与带阻滤波器的幅频特性如图 l和 u分别是下通带截止频率和上通带截止频率 s1和 s2分别为阻带的下限频率和上限频率 0为阻带中心频率 02 u l 阻带带宽B u l B作为归一化参考频率 相应的归一化边界频率为 u u B l l B s1 s1 B s2 s2 B 02 u l 根据 与 的对应关系 可得到 且 u l 1 p 1 上式称为低通到带阻的频率变换公式 将上式代入p j 并去归一化 可得 上式就是直接由归一化低通转换成带阻的频率变换公式 设计带阻滤波器的步骤 A 确定模拟带阻滤波器的技术要求 下通带截止频率 l 上通带截止频率 u阻带下限频率 s1 阻带上限频率 s2阻带中心频率 02 u l 阻带宽度B u l它们相应的归一化边界频率为 l l B u u B s1 s1 B s2 s2 B 20 u l以及通带最大衰减 p和阻带最小衰减 s B 确定归一化模拟低通技术要求 取 s和 s的绝对值较小的 s 通带最大衰减为 p 阻带最小衰减为 s C 设计归一化模拟低通G p D 将G p 转换成带阻滤波器H s 例 设计模拟带阻滤波器 其技术要求为 l 2 905rad s s1 2 980rad s s2 2 1020rad s u 2 1105rad s p 3dB s 25dB 试设计巴特沃斯带阻滤波器 解 1 模拟带阻滤波器的技术要求 l 2 905 u 2 1105 s1 2 980 s2 2 1020 02 l u 4 2 1000025 B u l 2 200 l l B 4 525 u u B 5 525 s1 s1 B 4 9 s2 5 1 20 l u 25 2 归一化低通的技术要求 3 设计归一化低通滤波器G p 4 带阻滤波器的H s 为 数字高通 带通和带阻滤波器的设计 例如高通数字滤波器等 具体设计步骤如下 1 确定所需类型数字滤波器的技术指标 2 将所需类型数字滤波器的技术指标转换成所需类型模拟滤波器的技术指标 转换公式为 3 将所需类型模拟滤波器技术指标转换成模拟低通滤波器技术指标 4 设计模拟低通滤波器 5 将模拟低通通过频率变换 转换成所需类型的模拟滤波器 6 采用双线性变换法 将所需类型的模拟滤波器转换成所需类型的数字滤波器 例6 6设计数字高通滤波器 通带截止频率 p 0 8 rad 通带衰减不大于3dB 阻带截止频率 s 0 44 rad 阻带衰减不小于15dB 采用巴特沃斯型滤波器解 1 数字高通的技术指标为 p 0 8 rad p 3dB s 0 44 rad s 15dB 2 模拟高通的技术指标计算如下 令T 1 则有 3 模拟低通滤波器的技术指标计算如下 将 p和 s对3dB截止频率 c归一化 这里 c p 4 设计归一化模拟低通滤波器G p 模拟低通滤波器的阶数N计算如下 查表 得到归一化模拟低通传输函数G p 为 为去归一化 将p s c代入上式得到 5 将模拟低通转换成模拟高通 将上式中G s 的变量换成1 s 得到模拟高通Ha s 6 用双线性变换法将模拟高通H s 转换成数字高通H z 实际上 5 6 两步可合并成一步 即 例6 7设计一个数字带通滤波器 通带范围为0 3 rad到0 4 rad 通带内最大衰减为3dB 0 2 rad以下和0 5 rad以上为阻带 阻带内最小衰减为18dB 采用巴特沃斯型模拟低通滤波器 解 1 数字带通滤波器技术指标为通带上截止频率 u 0 4 rad通带下截止频率 l 0 3 rad阻带上截止频率 s2 0 5 rad阻带下截止频率 s1 0 2 rad通带内最大衰减 p 3dB 阻带内最小衰减 s 18dB 2 模拟带通滤波器技术指标如下 设T 1 则有 通带中心频率 带宽 将以上边界频率对带宽B归一化 得到 u 3 348 l 2 348 s2 4 608 s1 1 498 0 2 804 3 模拟归一化低通滤波器技术指标 归一化通带截止频率 p 1 p 3dB s 18dB 4 设计模拟低通滤波器 5 将归一化模拟低通转换成模拟带通 6 通过双线性变换法将Ha s 转换成数字带通滤波器H z 下面将 5 6 两步合成一步计算 将上式代入 5 中的转换公式 得 将上面的p等式代入G p 中 得 例6 8设计一个数字带阻滤波器 通带下限频率 l 0 19 阻带下截止频率 s1 0 198 阻带上截止频率 s2 0 202 通带上限频率 u 0 21 阻带最小衰减 s 13dB l和 u处衰减 p 3dB 采用巴特沃斯型 解 1 数字带阻滤波器技术指标 l 0 19 rad u 0 21 rad p 3dB s1 0 198 rad s2 0 202 rad s 13dB 2 模拟带阻滤波器的技术指标 设T 1 则有 阻带中心频率平方为 20 l u 0 421阻带带宽为B u l 0 07rad s 将以上边界频率对B归一化 l 8 786 u 9 786 s1 9 186 s2 9 386 20 l u 85 98 3 模拟归一化低通滤波器的技术指标 p 1 p 3dB 4 设计模拟低通滤波器 5 将G p 转换成模拟阻带滤波器Ha s 6 将Ha s 通过双线性变换 得到数字阻带滤波器H z 6 6IIR数字滤波器的直接设计法 1 零极点累试法称为零极点累试法 在确定零极点位置时要注意 1 极点必须位于z平面单位圆内 保证数字滤波器因果稳定 2 复数零极点必须共轭成对 保证系统函数有理式的系数是实的 例6 9图 a 零极点分布 b 幅度特性 2 在频域利用幅度平方误差最小法直接设计IIR数字滤波器设IIR滤波器由K个二阶网络级联而成 系统函数用H z 表示 式中 A是常数 ai bi ci di是待求的系数 Hd ej 是希望设计的滤波器频响 如果在 0 区间取N点数字频率 i i 1 2 N 在这N点频率上 比较 Hd ej 和 H ej 写出两者的幅度平方误差E为 式中共有 4K 1 个待定的系数 求它们的原则是使E最小 下面我们研究如何求出 4K 1 系数 E是 4K 1 个未知数的函数 用下式表示 上式 表示4K个系数组成的系数向量 为推导公式方便 令 为选择A使E最小 令 设 k是 的第k个分量 ak或bk或ck或dk 因为 式中H i表示对Hi函数共轭 将上式具体写成对ak bk ck dk的偏导 得到 式中 k 1 2 3 K i 1 2 3 N 同理求得 由于系统函数是一个有理函数 极 零点均以共轭成对的形式存在 对于极点z1 一定有下面关系 a 要求的幅度特性 b k 1 2时的幅度特性 例设计低通数字滤波器 其幅度特性