数字滤波器设计之思路或vhdl.docVIP

§ 11.2 FIR数字滤波器设计 FIR数字滤波器主要用来实现信号预处理、防混叠、带选、抽选/插补、滤波和视频卷积等功能的，以下简要介绍其结构原理和设计方案。 11.2.1 FIR滤波器结构与原理简要 N阶的FIR滤波器系统的传递函数为 （11-1） 它有N-1阶极点在z=0处，有N-1个零点位于有限z平面的任何位置。式（11-1）的系统差分方程表达式为 （11-2） 式(11-2)就是输入序列x(n)与单位冲击响应h(n)的线性卷积。由式(11-2)可知n时刻的输入y(n)仅于n时刻的输入以及过去N-1个输入值有关。可以直接画出其网络结构，如图11-2所示，这种结构称直接型。 如果FIR滤波器的单位冲击响应h(n)对称，例如单位冲击响应h(n)满足： 0(1(…( (11-3) 则该因果系统具有严格的线性相位。当为偶数时，式(11-2)可化为 （11-4） 把式（11-3）代入式（11-4）得 （11-5） 把式（11-5）画成网络结构，如图11-3。同理，当为奇数时，式(11-2)可化为 （11-6） 把式（11-6）画成网络结构，如图11-4。 设计FIR滤波器的方法有窗函数法、频率取样法和等波纹优化设计法等。窗函数法是最简单的设计方法，也称为傅里叶级数法。先给定所要求的理想滤波器频率响应，要求设计一个FIR滤波器频率响应来逼近。但设计是在时域进行的，因而先对进行傅里叶反变换得到： （11-7） 由于是矩形频率特性，故一定是无限长的序列，且是非因果的，而要 图 11-4 设计的是FIR滤波器，其必然是有限长的序列，所以要用有限长的来逼近无限长的，最有效的方法是截断，或者用一个有限长的窗口函数序列来截取，即： （11-8） 因而窗口函数序列的形状及长度的选择就很关键。窗口函数主要有以下几种： 矩形窗： 汉宁（Hanning）窗： （11-9） 海明（Hamming）窗，又称为升余弦窗： （11-10） 它是汉宁窗的改进，可以得到旁瓣更小的效果。 布拉克曼（Blackman）窗，又称二阶升余弦窗： （11-11） 在这个窗函数中加上了二次谐波分量，可以进一步抑制旁瓣。另外还有凯泽窗和三角窗等。 11.2.2 FIR滤波器设计方案确定 FIR滤波器的实现方案可以有多种，有的在速度上有优势，有的在资源上占优势，有的则在结构上占优势，以下介绍几种方案，以作比较选择。 1. 方案一 图11-5 FIR滤波器设计方案一 这个方案利用了具有严格线性FIR滤波器的单位冲击响应的对称特性，对图11-2所示的直接型结构进行简化，得到如图11-3和图11-4所示的结构。这样就可以减少一半的乘法运算，节省相当一部分资源。具体结构如图11-5所示（由图11-4的结构而来），运算采用并行方式进行。但是这种结构只能被采用在M为偶数时，当M为奇数时，则必须换用类似于图11-3的结构。可见这种结构对于固定参数的FIR滤波器比较适合，对于参数可变的FIR滤波器就不大适用了。 2．方案二 此方案就是引用如图11-2所示的直接型，它没有经过任何化简，运算也是采用并行的方式进行，设计结构如图11-6所示。这种方案速度快，通用性强，适合于参数可变的FIR滤波器；但是对于实现N阶8位的FIR滤波器，要有N个8×8个乘法器，从而占用大量的资源，对于资源相对紧张的FPGA来说，采用这种方案不够实际。 图11-6 FIR滤波器设计方案二 3．方案三 为了能使FPGA实现可变参数的FIR滤波器，要求其网络结构必须是通用型的，而且不能使用过多的运算，以利节省逻辑资源。基于这点，方案三在方案二的基础上进行改进，网络结果仍采用直接型，在进行乘法运算时，用串行乘法来代替并行乘法，这样以牺牲时间为代价换得尽可能少的占用资源。其中的×的串行乘法器是用次移位相加的方式实现的。 以下是串行乘法的推导： 图11-7 串行乘法器模块图 假设数据、是两个位的乘数，是乘积，，其中表示乘数的第位数据，则有： 乘法器结构如图11-7所示。这一方案总体上是可行的，只是其中的个乘法器仍占用了较多的资源。 4．方案四 方案四与前3个方案有较大的不同，它不采用乘法器进行乘法运算，而是采用乘法表，用查表来得到两数相乘的积，即把两乘数作为ROM的地址，对应单元的数据就是该两乘数的乘积（RO