数字滤波器设计演示幻灯片.pptVIP

第六章　数字滤波器设计 ;6.1 引 言 ;　　同样，与模拟滤波器类似，数字滤波器按频率特性也有低通、高通、带通和带阻等之分，滤波器的性能指标通常也习惯在频域给出。常用数字滤波器的幅度特性示意图如图6.1所示。与模拟滤波器不同的是，由于序列的傅里叶变换具有以2π为周期的周期性，因此，数字滤波器的频率响应也有这种周期性。 低通滤波器的通带处于0或2π的整数倍频率附近，高通滤波器的通带则处于π的奇数倍频率附近。 ;图6.1 各种数字滤波器的幅度特性 ;图6.1 各种数字滤波器的幅度特性 ;　　图6.1所示的理想滤波器的幅度特性有理想、陡截止的通带和无穷大衰减的阻带两个范围，这显然是无法实现的，因为它们的单位取样响应均是非因果和无限长的。实践中只能用一种因果可实现的滤波器去与之逼近，使其满足给定的误差容限。一个实际滤波器的幅度特性在通带中允许有一定的波动，阻带衰减则应大于给定的衰减要求，且在通带与阻带之间允许有一定宽度的过渡带。 ;　　图6.2示出了一个实际低通滤波器的幅度特性， 特性曲线中有通带、过渡带和阻带三个区间。通带范围是0≤ω≤ωp， 在通带内，幅度特性以误差δ1逼近于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滤波器设计方法 ;　　（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(s)从S平面变换成Z平面的数字滤波器H(z)。若令s=σ+jΩ，z=rejω，代入式（6-29)得到rejω=eσTejΩT， 则有 ;　　将式（6-28）示于图6.5，S平面上每一条宽度为2π/T的横向条带区都将重叠地映射到整个Z平面上，而第一条横向条带的左半边映射到Z平面单位圆内，右半边映射到Z平面的单位圆外，S平面的虚轴映射到Z平面的单位圆上。当模拟频率Ω从-π/T变化到π/T 时，数字频率ω则从-π变化到π，且由式（6-30b）可知，ω与Ω之间成线性关系。由