基于DSP的FIR数字滤波器的实现
加入技术交流群
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个样本后的存储情况。本文引用地址：https://www.eepw.com.cn/article/150927.htm
评论