FIR滤波器的原理及设计.docxVIP

选题?2?实验讲义 实验名称：基于分布式算法的?FIR?滤波器设计 1．数字滤波器基础知识 （（ （ （ （ 无线通信等领域都有着非常广泛的应用。通过滤波运算，将一组输入数据序列转变为另一组 输出数据序列，从而达到修正时域或频域中信号属性的目的。数字滤波器就是用于完成这种 信号滤波功能，用有限精度算法来实现的一种离散时间线性时不变（LTI）系统。相比于模 拟滤波器，数字滤波器具有以下优点：?1）数字滤波器的频域特性容易控制，性能指标优良； （2）数字滤波器可以工作在极低的频率，可以方便地实现模拟滤波器难以实现的线性相位 系统；?3）数字滤波器工作稳定，一般不会受到外部环境的影响；?4）数字滤波器的灵活性 和可重用性高，只需要简单编程就可以修改滤波器的特性，设计周期短。数字滤波器的实现 可以采用专用?DSP?芯片，通过编写程序，利用软、硬件结合完成滤波器设计，也可以采用 市面上通用的数字滤波器集成电路来实现，但这两种方法无法适应高速应用场合。随着集成 电路技术的高速发展，FPGA?应用越来越普及，FPGA?器件具有芯片密度大、执行效率高， 速度快，集成度高等优点，用?FPGA?芯片作为滤波器的设计载体，可以实现高速信号滤波 功能。 1.1?FIR?数字滤波器特点 数字滤波器通常分为?IIR（无限冲激响应）和?FIR(有限冲激响应)两种。FIR?滤波器具有 以下特点：（1）可以做成严格的线性相位，同时又可以具有任意的幅度特性（2）单位冲激 响应是有限长的，所以一定是稳定的，因此在实际中得到广泛的应用。 1.2?FIR?滤波器结构 设?FIR?滤波器的单位冲激响应为?h(n)?，?0???n???N???1?， 系统函数??H?( 系统函数??H?(Z?)??????h(n)?z ?n 差分方程形式为：?y( 差分方程形式为：?y(n)?????h(k?)?x(n???k?)? （1） N??1 k??0 基本结构（直接型）： s(??Ny(n)?? s(??N y(n)??????h(k?)s(n???k?)?? （2） 1.3?结构简化 当?FIR?滤?波?器?具?有?严?格?线?性?相?位?时?，?h(n)?满?足?h(n)???h(?N???1???n)?或 h(n)????h(?N???1???n)?设?N?为偶数，令 s(0)???x(n)???x[n???(?N???1)] s(1)???x(n???1)???x[n???(?N???2)] ? N N ??1)???x[n???( ??1)]???x(n?? ) 2 2 2 则（1）式可以简化为： N ?1 2 k??0 与（1）相比，所需乘法器数量降为原来的一半。 y(n)?? y(n)?????h(k?)?x(n???k?)??????h(?L???l???m)?x[n???(?L???l???m)] 1.4?设计规模的改进 当滤波器的长度?N?增加时，可将求和分配到几个独立的?M?阶并行?DA?的?LUT?中。以（1） 式为例，设?N???L???M?， N??1 L?1?M??1 k??0 l??0?m?0 这样，长度为?N?的滤波器实现可以分解为?L?个长度为?M?的滤波器，可运用流水线加法器累加 结果。 即?s(n? k)???[s??(n 即?s(n? k)???[s??(n k) s?(n k)s?(n k)] 在很多?DSP?应用场合中，滤波器系数一般为常数，在这种情况下，可以利用分布式 （Distributed?Algorithm，DA）算法原理将求乘积和运算转变为移位和加法运算来实现。 在滤波器规模比较小的情况下，采用?DA?算法可以减少电路规模，更容易实现流水处理，从 而有更高的执行效率。以（2）式为例，DA?算法的基本思想为： 在（2）式中，设h(k)为常数，将?s(n k)表示为?B+1?位二进制补码形式，其中，最高位为 符号位，用?s?(n k)表示，下标?B?表示第?B?位，数值位共有?B?位，用?s?(n k)?表示，其中 B b 0 b B 1?。 B 1 0 补 将式中各二进制位按权展开，则有下式成立， s(n k) 2?B?s?(n k) B B?1 s?(n?k)?2b b b?0 代入（2）式，并改变式中求和的次序，可得 y(n) N 2 1  h(k)s(n?k) h(k)[?2? h(k)[?2?B?s??(n k) N 1 2 B?1 B  s?(n?k)?2b?] b 2?B?[?? h(k)s 2?B?[?? h(k)s??(n k)]??? h(k)?? s?(n k)?2b N N 1 1 2 2 B?1 B b k?0 k?0 b?0 2?B? 2?B?[ h(k)s??(n k)]