数字滤波器设计PPT.ppt

第六章　数字滤波器设计 ;6.1 引 言 ;　　同样，与模拟滤波器类似，数字滤波器按频率特性也有低通、高通、带通和带阻等之分，滤波器的性能指标通常也习惯在频域给出。常用数字滤波器的幅度特性示意图如图6.1所示。与模拟滤波器不同的是，由于序列的傅里叶变换具有以2π为周期的周期性，因此，数字滤波器的频率响应也有这种周期性。 低通滤波器的通带处于0或2π的整数倍频率附近，高通滤波器的通带则处于π的奇数倍频率附近。 ;图6.1 各种数字滤波器的幅度特性 ;图6.1 各种数字滤波器的幅度特性 ;图6.2 实际低通滤波器的幅度特性;　　通带内衰减（波动）和阻带衰减（波动）通常用分贝表示， 对于图6.2， 我们令 ;　　滤波器的频率特性除了幅度特性外，还有相位特性θ(ω)。 一般对θ(ω)并无过多要求，只要保证滤波器稳定就可以了。 但在有些场合要求θ(ω)具有一定的性质，如线性相位特性， 即要求θ(ω)=-τω（τ为延时常数）等。 　　数字滤波器可用N阶差分方程 ;　　数字滤波器按其单位采样响应长度可分为无限冲激响应滤波器（IIR）和有限冲激响应滤波器（FIR）两类。按照滤波器的实现方式则又可以分为递归滤波器和非递归滤波器两类。式(6-5)中若ak=0(k=1,…,N)， 该滤波器为FIR滤波器；若ak≠0(k=1, …,N)，该滤波器为IIR滤波器。一般情况下递归滤波器对应于IIR滤波器， 而非递归滤波器对应于FIR滤波器。 ;6.2 IIR滤波器设计方法 ;　　(2) 滤波器系统函数的零点和极点位置完全决定了滤波器的幅度和相位响应。所以，通过合理设置数字滤波器系统函数的零、极点，即可得到符合要求的滤波特性。 这种方法往往需要多次调整零、极点位置，称为零、极点累试法。 　 ;　　　　(3) 计算机辅助设计法。 这是一种最优化设计方法。 它先确定一种最优化准则， 例如设计出的实际频率响应的幅度与理想频率响应的幅度的均方误差最小准则，或它们的最大误差最小准则等， 然后确定满足该最佳准则的滤波器系数ak、bi。 这种设计一般不易得到滤波器系数的显式表达式， 而是需要进行大量的迭代运算，需用计算机辅助设计完成。 ;6.2.1 根据模拟滤波器设计IIR数字滤波器 　　目前IIR数字滤波器设计中用得较多的是借助于模拟滤波器的设计方法。借用模拟滤波器设计有一套相当成熟的方法， 可以给数字滤波器的设计带来很大方便。该方法的设计步骤是： 　　（1） 若所需设计的数字滤波器是低通的，按一定规则先将给出的数字低通滤波器的技术指标转换为模拟低通滤波器的技术指标。 　　（2） 根据转换后的技术指标设计模拟低通滤波器的传递函数Ha(s)。  ;　　（3） 再按一定规则将Ha(s)转化成H(z)，完成低通滤波器的设计。  　　(4) 若所设计的是高通、带通或带阻滤波器，那么还需将高通、带通或带阻数字滤波器的技术指标转化为低通模拟滤波器的技术指标，然后按步骤（2）设计低通模拟滤波器 Ha(s)，再将Ha(s)最终转换为所需的H(z)。 ;　　1. 选定模拟低通滤波器原型 　　由滤波器理论，高通、带通、 带阻滤波器均可以利用变量变换方法，分别由低通滤波器变换得到。 所以要先根据要求选定模拟低通滤波器的设计方法。 模拟低通滤波器有巴特沃兹（Butterworth）、切比雪夫（Chebyshev）和椭圆（Elliptic）滤波器等。  　　;　　2. 由模拟滤波器完成IIR数字滤波器设计 　　利用模拟滤波器来完成数字滤波器的设计， 就是先确定模拟滤波器的Ha(s)进而确定数字滤波器的系统函数H(z)。 它实际上是由S平面到Z平面间的一种映射转换，此时必须满足两种基本要求：  　　（1）H(z)的频率响应要能模仿Ha(s)的频率响应，即S平面的虚轴必须映射到Z平面的单位圆上。 ;　　（2） 因果稳定的模拟滤波器Ha(s)转换成数字滤波器H(z)，仍是因果稳定的。也就是S平面左半平面（Re［s］0）应该映射到Z平面的单位圆以内（|z|1）。  　　将系统函数Ha(s)从S平面映射到Z平面可以有多种方法， 工程上常用的是冲激响应不变法（或叫脉冲响应不变法）和双线性变换法， 下面我们对这两种方法略作介绍。 ;　　1) 冲激响应不变法 　　冲激响应不变法是使数字滤波器的单位采样响应序列h(n)模仿模拟滤波器冲激响应ha(t)。将模拟滤波器的冲激响应加以等间隔（间隔为T）采样，并使h(n)正好等于ha(t)的采样值， 即满足 ;　　下面我们分析采用冲激响应不变法时，S平面和Z平面之间的映射关系。若令Hz(s)是ha(t)的拉普拉斯变换，H(z)为h(n)的Ｚ变换，则利用已知的序列的Ｚ变换与模拟信号的拉普拉斯变换的关系， 可得 ;　　将模拟滤波器的Ha(