数字滤波器的设计

nullnull第七章 数字滤波器的 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 滤波器在实际信号处理中起到了非常重要的作用。任何检测的信号都含有噪声，而滤波是去除噪声的基本手段。本章介绍IIR滤波器和FIR滤波器的设计。IIR滤波器设计主要内容包括：巴特沃思、切比雪夫模拟低通滤波器设计；脉冲响应不变法和双线性变换法的数字化变换 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 ；数字高通、带通和带阻滤波器的设计。而FIR滤波器是直接采用的数字式设计方法。针对FIR滤波器特征，首先介绍了其线性相位的实现条件，然后介绍了窗函数法和频率抽样法的设计方法。null§ 7.1 概述7.1.1 滤波器的分类 这里主要讨论经典滤波器的设计。 按功能划分经典滤波器又可分为低通、高通、带通、带阻四种滤波器 图7-1 理想低通、高通、带通和带阻滤波器幅度特性 null 经典滤波器设计从实现方法上分为IIR滤波器和FIR滤波器。它是一个线性时不变离散时间系统，如果滤波器用单位脉冲响应序列 表示，其输入 与输出 之间的关系可以表示为： 的Z变换称为系统函数。IIR滤波器和FIR滤波器的系统函数分别是：nullδ1 ：通带衰减 δ2 ：阻带衰减 ωc ：通带截止频率 ωst：阻带截止频率 ωc －ωst ：过渡带 7.1.2 数字滤波器的性能要求 一个理想滤波器，要求所在通频带内幅频响应是一常数；相位频率相应为零或是频率的线性函数。但一个实际的滤波器要是不可能得到上述幅频和相频响应。以低通滤波器为例，频率响应有通带、过渡带及阻带三个范围。null图7-2 逼近理想低通滤波器得容限图通带最大衰减。null当幅度下降到 ， ，此时 ，称 为3dB通带截止频率。阻带最小衰减7.1.3 数字滤波器设计方法概述7.1.3 数字滤波器设计方法概述 设计IIR数字滤波器一般有以下两种方法： 1. 模拟滤波器:首先设计一个合适的模拟滤波器，然后将它转换成满足给定指标的数字滤波器，这种方法适合于设计幅频特性比较规则的滤波器，例如低通、高通、带通、带阻等。 2.直接在频域或者时域中进行数字滤波器设计，由于要联立方程，设计时需要计算机作辅助设计。 null7.2无限脉冲响应数字滤波器的设计1.设计的一般方法 IIR滤波器以模拟低通滤波器为基础的设计方法，为了设计其他的选频滤波器（高通，带通，带阻等），需要对低通滤波器进行频率转换，在设计过程中有两种不同的变换,频带变换和模拟/数字变换。根据这两种变换的先后次序，引出两种设计方法。 null图7-3 IIR滤波器的设计 流程 快递问题件怎么处理流程河南自建厂房流程下载关于规范招聘需求审批流程制作流程表下载邮件下载流程设计 2.巴特沃思低通滤波器 （1）基本性质 巴特沃思滤波器以巴特沃思函数来近似滤波器的系统函数。巴特沃斯滤波器是根据幅频特性在通频带内具有最平坦特性定义的滤波器。 null(7.2.1).下面归纳了巴特沃斯滤波器的主要特征 a.对所有的N， 。 巴特沃思滤波器的低通模平方函数表示b.  对所有的N， 即c. 是 的单调下降函数。d. 随着阶次N的增大而更接近于理想低通滤波器。null 如图7-4，可以看出滤波器的幅频特性随着滤波器阶次N的增加而变得越来越好，在截止频率 处的函数值始终为1/2的情况下，通带内有更多的频带区的值接近于1；在阻带内更迅速的趋近于零。 图7-4 巴特沃思低通滤波平方幅频特性函数null 在以后的设计和 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 时，经常用归一化巴特沃思低通滤波器为原型滤波器，一旦归一化低通滤波器的系统函数确定后，其它巴特沃思低通滤波、高通、带通、带阻滤波器的传递函数都可以通过变换法从归一化低通原型的传递函数 得到。归一化原型滤波器是指截止频率 已经归一化成 的低通滤波器。对于截止频率为某个 的低通滤波器，则令 代替归一化原型滤波器系统函数中的 ，即 对于其他高通、带通、带阻滤波器，可应用后面讨论到的频带变换法，由其变换得出。 null(2)系统函数和极点分布 设巴特沃斯的系统函数 ，则频率响应是 令上式分母为零可以得到 的2N个极点Sk 并解得当N为偶数时则 nullN为奇数： 的极点均匀分布在s平面的单位圆上，共有2N个角度间隔为/N的极点，极点关于j轴对称，不会落在虚轴上。将左半平面的极点构成 。 (3)设计过程巴特沃思低通滤波技术指标关系式为 ； 为通带边界频率 ； 为阻带边界频率null代入（7.2.1） 化简后得 两式相比消去后得 由此得 取满足上式的最小整数N作为滤波器的阶数。 null将N 带入式（7.2.4）或式(7.2.5)可得截止频率 或 查表求得归一化传输函数 ,令 代替归一化原 型滤波器系统函数中的 , 即 代入 ，即得到实际滤波器传输函数。例7.1 设计一巴特沃思低通滤波器，使其满足以下指标: 通带边频 ,通带的最大衰减为 ，阻带边频为 ,阻带的最小衰减为 。 null解:滤波器技术指标为 ， ， , 确定阶次N，代入式（7.2.6） 取 N=4 查表得四阶巴特沃思多项式，得归一化系统函数表   达式 由式（7.2.7）得 用 替换式（7.2.8）中的s，构成巴特沃思滤波器传输函数H(s)为 null3. 切比雪夫I滤波器切比雪夫I型滤波器的幅度平方函数为 （1）基本性质 是N阶切比雪夫多项式，定义为 nullN=0，C0(x)=1 N=1, C1(x)=x N=2, C2(x)=2x2-1=2xC1(x)-C0(x)迭代公式：CN(x)=2xCN-1(x)-CN-2(x) N>1N=偶数，CN(x)为偶函数 N=奇数，CN(x)为奇函数null切比雪夫多项式曲线1) x:[0, 1] 设：=Ncos-1x , x: 01N=0 ,  =0 C0(x):1N=1 ,  :/2 0, C1(x):01N=2 ,  :    /2 0, C2(x):-1  01N=3 ,  : 3 /2 0, C2(x):0  -1  01N=4 ,  : 2 0, C2(x):1  0  -1  012) x:[1,  ) nullx=0, N=even, |CN(0)|=1, N=odd, CN(0)=0 x=1, CN(1)=1 |x|≤1, CN(x)在[-1, +1]之间波动， N增加，波动次数增加 |x|>1, CN(x)单调上升，N增加，上升速度增加。nullnull切比雪夫滤波器的幅频响应有如下特点： 1)0< < 时， 在1与 之间等幅波动， 越 大，波动幅度越大。 2) ＝0时，当N为奇数，则 ＝1 ，当N为偶 数，则 。 3) ＝ 时，对所有的N值， 都取同样的值 。 4) > 时，曲线单调下降, 越大，N越大，曲线衰 减越快。 null(2)设计过程 1)根据要求的滤波器指标确定波纹参数 和阶数N。 由允许的通带波纹 确定。则 滤波器的阶数N由阻带允许的衰减确定。 2)查表求得归一化传输函数 ,令 代替归 一化原型滤波器系统函数中的 ，即得到实际滤波 器传输函数 。 null1.变换原理 数字滤波器的单位脉冲响应序列h(n)正好等于模拟滤波器的冲激响应ha(t)的采样值，即 h(n)=ha(nT), T为采样周期.如以 Ha(s) 及 H（z）分别表示 ha(t) 的拉氏变换及 h(n) 的Z变换，即数字滤波器的系统函数 便是的Z变换 。 7.2.2脉冲响应不变法( Impulse Invariance ) 下面我们分析从模拟滤波器到数字滤波器S平面和Z平面之间的映射关系。 null设抽样信号抽样信号的拉氏变换序列h(n)的z变换nullnullnullnull(2) 抽样信号与原信号关系null2.混叠失真 利用抽样序列的Z变换与模拟信号的拉普拉斯变换的关系，得 由采样定律可知，如果模拟滤波器的频率响应带限于折叠频率 以内 这时数字滤波器的频响才能不失真地重现模拟滤波器的频响（存在于折叠频率 以内）null 但任何一个实际的模拟滤波器，其频率响应都不可能是严格带限的，因此不可避免地存在频谱的交叠，即产生频率响应的混叠失真。 原模拟信号的频带不是限于之间，则会在的奇数倍附近产生频率混叠，从而映射到Z平面上， 附近产生频率混叠。这种频率混叠现象会使设计出的数字滤波器在 附近的频率特性，程度不同的偏离模拟滤波器在 附近的频率特性，严重时使数字滤波器不满足给定的技术指标。 null解决混叠的方法： 1)滤波器指标以模拟域形式给出，此时已确定，采样频率T增加，混叠减小。未确定，但已定，采样频率增加，为保证不变，必有增加，增加滤波器的阶数N，混叠2)滤波器指标以数字域形式给出，此时减小。null 设模拟滤波器的系统函数若只有单阶极点，且分母的阶数高于分子阶数 N＞M，则可表达为部分分式形式 其拉氏反变换为 是单位阶跃函数，对ha(t)采样得到数字滤波器的单位脉冲响应序列 3.数字化设计null再对h(n)取Z 变换，得到数字滤波器的传递函数 第二个求和为等比级数之和，要收敛的话 必有所以 null例7.2 IIR低通滤波器的设计指标如下： 　　1）通带截止频率ωP＝0.1πrad， 阻带起始频率ωs＝0.25πrad 2）通带最大衰减αp=3dB， 阻带最小衰减αs=15dB4.设计举例解：1) 根据DF指标, 将DF指标转换为归一化LPF指标nullnull4) 求H(z)null7.2.3双线性变换法 （Bilinear Transform)1.变换原理 s平面到z平面的映射关系 二次映射法 为了将S平面的jΩ轴压缩到S1平面jΩ1轴上的 到 一段上，可通过以下的正切变换实现： （7.2.11） 这里C是待定常数，下面会讲到用不同的方法确定C 。 null 当 由 时, 由 经过０变化到 ，即S平面的整个 轴被压缩到S1平面的 一段. 式（7.2.11）又可以写成 将这一关系解析扩展至整个S平面， 则得到S平面到S1平面的映射关系： 再将 S1 平面通过 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 变换关系映射到Z平面，即令null从而得到s平面与z平面的单值映射关系 (7.2.12)null2.变换常数的选择 式（7.2.11）中常数的选择可以使模拟滤波器的频响特性和数字滤波器的频响特性在不同的频率范围有对应的关系，起到调节二者频带间关系的作用。选择的方法有两种。 1)使模拟滤波器和数字滤波器的频响特性在低频部分有较确切的对应关系，即当 较小时，有由此得 2)使数字滤波器的某一特定频率(例如截止频率 与模拟原型滤波器的特定频率 严格对应，即null由于在待定的模拟频率和待定的数字频率处频率响应应严格相等，因而可以较准确的控制截止频率位置。 nullnullS/2-S/2null3. 模拟角频率Ω和数字角频率ω的映射关系nullnull总结计算H(Z)步骤如下： 设给定数字低通滤波器的通带截止频率 、阻带截止频率 、通带波动 和阻带波动 。 (1)利用公式 对通带和阻带截止频率 和 进行预畸变，求出模拟低通滤波器的通带和阻带截止频率 和 。预畸变函数式为 (2)求满足指标 、 、 和 要求的模拟低通滤波器的传输函数 。 (3)利用双线性变换公式 将 映射成 4.数字化设计 null例7.3 设计IIR滤波器代替如下性能的模拟LPF: 1)fp=50Hz fs=125Hz 2)αp≤3dB, αs≥15dB 3)采样频率f=1kHz解: 1)用脉冲响应不变法 ; 2)双线性变换法nullnullnullnull(3)数字带通的设计a.思路:归一化LPF→模拟BPF→数字BPFnull将代入上式nullnullnullnull归一化LPF指标为：null 前面我们学习了模拟低通滤波器，数字低通滤波器的设计，对于数字高通和带阻的设计，可以借助于模拟滤波器的频率变换设计一个所需类型的模拟滤波器，再通过双线性变换将其换算成所需类型的数字滤波器。1.高通滤波器 和 分别是低通的归一化通带截止频率和归一化阻带截止频率， 和 分别是高通的归一化通带下限频率和归一化阻带上限频率。 7.2.4数字高通、带通和带阻滤波器设计null低通的 从 经过 和 到0时，高通的 则从0经过 和 到 ，因此 和 之间关系为 ，即是低通到高通的频率变换公式。 null总结步骤为： (1)确定数字高通滤波器的技术指标 、 。 (2)将数字高通滤波器的技术指标转换成高通模拟滤 波器的技术指标 、 ，转换公式为 。 (3)利用频率变换 将模拟高通滤波器技术指标转 换成归一化模拟低通滤波器 的技术指标。 (4)设计模拟低通滤波器 。 (5)将模拟低通滤波器 通过频率转换转换成模拟高 通滤波器 ，并去归一化后得 。 null(6)采用双线性变换，将所需类型的模拟滤波器转换成所需类型的数字滤波器. 2.带通滤波器 、 分别是模拟带通滤波器通带的下限和上限频率， 是下阻带的上限频率， 是上阻带的下限频率，令 为通带带宽，用B做为归一化参考频率， null令 为通带的中心频率，归一化 可以找到 和 的转换关系。 null , ，所以有 由此可以得到低通滤波器的技术指标 、 ，可以设计低通滤波器的转移函数 。 ，将 代入得到 去归一化，将 代入上式得： 可得模拟带通滤波器的传递函数是 null数字带通滤波器的传递函数3.带阻滤波器 、分别是模拟带阻滤波器通带的截止频率， 是下限频率， 是上限频率，与带通模拟滤波器一样，令 为通带带宽，用B做为归一化参考频率， null 和 的转换关系： , ，所以有null由带通滤波器的推导过程易知 可得模拟带阻滤波器的传递函数是 数字带阻滤波器的传递函数 null例7.6 设计一个数字高通滤波器，要求通带截止频率ωp=0.8πrad，通带衰减不大于3dB，阻带截止频率ωs=0.5πrad,阻带衰减不小于18dB。采用巴特沃斯型滤波器。 解：(1)数字高通的技术指标为 ωp=0.8πrad,αp=3dB; ωs=0. 5πrad,αs=18dB (2) 模拟高通的技术指标计算如下： 令T=2，则有 null(3)设计归一化模拟低通滤波器G(p)。模拟低通滤波器的阶数N计算如下。 取N＝2 null(4)查表得到归一化模拟低通传输函数G(p)为 为去归一化，将p=s/Ωc代入上式得到： (5)将模拟低通转换成模拟高通。将上式中G(s) 的变量换成1/s，得到模拟高通Ha(s)： (6)用双线性变换法将模拟高通H (s)转换成数字高通H(z)：