基于DSP的FIR数字滤波器的实现

liyf

摘要：在数字信号处理应用中，滤波占有十分重要的地位，如对信号的过滤、检测、预测等，都要广泛地用到滤波器。文中研究了FIR滤波器窗函数算法的基本思想给出了在定点DSP芯片上实现FIR数字滤波器的设计方法，并给出了仿真结果。
关键词：DSP；FIR滤波器；循环寻址；汉宁窗

0 引言
    数字信号处理现已在通信与信息系统、信号与信息系统、自动控制、需达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理应用中，滤波占有十分重要的地位，如对信号的过滤、检测、预测等，都要广泛地用到滤波器。IIR数字滤波器的设计保留了一些典型模拟滤波器优良的幅度特性，但所涉及的滤波器相位特性一般是非线性的，而FIR滤波器则可在保证幅度特性并满足技术要求的同时，也很容易做到严格的线性相位特性。

1 基于窗函数法的FIR滤波器设计
1．1 单位冲激响应
    首先应根据技术要求确定待求滤波器的单位冲激响应hd(n)。如果给出待求滤波器的频率为

，那么单位取样响应则可用下式求出：
   


    当

较复杂或不能用封闭公式

表示时，就不能用上式求出hd(n)。此时可以对

从ω=0到ω=2π采样M个点， 采样值为

，k=0，1，2，…，M-1，并用2π/M代替(1)式中的dω，此时(1)式可近似写成：
   


    这样，根据频率采样定理，hM(n)与hd(n)应满足如下关系：
   


    因此，如果M选得较大，就可以保证在窗口内，hM(n)有效逼近hd(n)。实际计算(2)式时，可以用

的M点采样值，来进行M点IDFT
(IFFT)，从而得到(2)式的值。
    如果给出通带阻带衰减和边界频率的要求，则可选用理想滤波器作为逼近函数，从而用理想滤波器的特性作傅立叶逆变换，以求出hd(n)。若理想低通滤波器为：
   


    那么，求出的单位取样响应hd(n)则为：
   


    为保证线性相位，通常应取α=(N-1)／2。

                          
                       
                          
                               
1．2 过渡带及阻带衰减
    根据对过渡带及阻带衰减的要求，设计时可选择窗函数的形状，并估计窗口长度N。设待求滤波器的过渡带用△ω表示，它近似等于窗函数的主瓣宽度。由于过渡带△ω近似与窗口长度N成反比。即N=A／△ω，其中A决定于窗口形式，例如，矩形窗A=4π，哈明窗A=8π等。按照过渡带及阻带衰减情况，选择窗函数形式。其设计原则是在保证阻带衰减的情况下，尽量选择主瓣比较窄的窗函数。
1．3 单位取样响应的计算
    计算滤波器的单位取样响应h(n)时，可按正式进行：
   


    式中，ω(n)是上面选择好的窗函数。如果要求线性相位，则要求hd(n)和ω(n)均对(N-1)／2对称，前面介绍的几种窗函数已保证对(N-1)
／2偶对称。而如要求h(n)对(N-1)／2奇对称，只要保证hd(n)对(N-1)／2奇对称即可。
    验算技术指标是否满足要求时其设计出的滤波器频率响应可采用下式进行计算：
   


    计算上式时，可使用FFT算法。如果

不能满足要求，那么，根据具体情况，可重复上述设计，直到满足要求为止。

2 基于DSP的FIR数字滤波器实现方案
2．1 滤波系统的差分方程
    若输入信号x(n)和滤波器的单位冲激响应h(n)在频域分别为

和

，则其输出信号的频率响应为

。根据离散傅氏变换的性质，可以得到滤波系统的差分方程

。其中，h(n)的计算可依据滤波器的参数，并借助MATLAB仿真工具计算得到。
2．2 实现方案
    由于在卷积运算之前，h(n)的N个数值是已知的，因此，可先在程序存储器中开辟N个单元来存放h(n)。由于输入序列x(n)是不断变化的，因此，在数据存储器中可开辟N个存储单元并对其进行移位寄存，其初始值分别为x(n)、x(n-1)……x(n-N-1)，然后采用循环寻址的方法对其进行访问。每次输入新的样本时，可以以新样本改写滑窗中的最老的数据，而滑窗中的其他数据则不需要移动。利用片内8 kB(循环缓冲区长度)寄存器可对滑窗进行间接寻址，循环缓冲区地址首尾相邻。8级循环缓冲区的结构如图1所示，其中顶部为低地址，图1中的(a)(b)(c)分别为初始状态、输入1个和2个样本后的存储情况。




                          
                       
                          
                               

3 仿真结果
    下面是一组信号采样序列样本x(n)，其中存在有高频干扰。现以x(n)作为输入序列，然后滤除其中的干扰成分。



    本设计的线性相位低通FIR数字滤波器的截止频率ωc为0．2πrad，窗口长度N为11。根据上述原理及实现方案，若采用汉宁窗来实现，则可依据算式计算出用汉宁窗设计时的各h(n)的系数：



    在CCS2．0软件中观察x(n)的输入曲线如图2所示，图3所示是其y(n)输出曲线。




                          
                       
                          
                               
    由图2和图3两图对比可见，经过滤波后，其输入曲线变平滑了。并且，根据计算所得：







4 结束语
    通过实验结果可见，本设计用DSP设计的FIR数字滤波器的输出曲线对应值与计算所得值完全相符，由此证明，本设计的程序是正确的。