DSP的FIR设计(低通滤波)C语言编写要点.pdfVIP

一、设计目的 低通滤波器设计。 本设计中使用的信号为 信息信号： signal=sin(2*pi*sl*n*T) 高频噪声 1：noise1=0.7*sin(2*pi*ns1*n*T) 高频噪声 2：noise2=0.4*sin(2*pi*ns2*n*T) 混合信号： x=(signal+noise1+noise2) 其中 sl=500Hz，ns1=3000Hz，ns2=8000Hz，T=1/20000。混合信号波形为滤波器 输入信号波形， 信息信号波形为输出信号波形， 滤波器的效果为滤除两个高频噪 声。 二、FIR 滤波器基本理论 （1）FIR 滤波器的特点 数字滤波器的功能， 就是把输入序列通过一定的运算变换成输出序列。 它的 实现方法有很多，其中比较常用到的是无限长脉冲响应滤波器 IIR 和有限长脉冲 响应滤波器 FIR 两种。 在计算量相等的情况下， IIR 数字滤波器比 FIR 滤波器的幅频特性优越，频 率选择性也好。但是，它有着致命的缺点，其相位特性不好控制。它的相位特性 f( ) argH( ej ) 是使频率产生严重的非线性的原因。但是在图像处理、数据传 输等波形传递系统中都越来越多的要求信道具有线性的相位特性。 在这方面 FIR 滤波器具有它独特的优点， 设 FIR 滤波器单位脉冲响应 h(n)长度为 N ，其系统函 数 H(z)为 N 1 n H (z) h(n)z n 0 1 H(z)是 z 的(N-1) 次多项式， 它在 z 平面上有 (N-1)个零点， 原点 z=0 是(N-1) 阶重 极点。因此， H(z)永远稳定，它可以在幅度特性随意设计的同时，保证精确、严 格的线性相位。 （2）FIR 滤波器的基本结构 数字滤波是将输入的信号序列，按规定的算法进行处理，从而得到所期望的 输出序列， FIR滤波器的差分方程为： N 1 y(n) a x(n k ) k k 0 对上式进行 Z变换得到 FIR滤波器的传递函数为： N 1 Y z k H z b z k X z i 0 由上式可以看出， H(z)是 z 1 的N-1次多项式，它在 z平面内有 N-1个零点，同 时在原点处有 N-1个重极点。N阶滤波器通常采用 N个延迟单元、N个加法器与 N+1 个乘法器，取图中 (a)、 (b)两种结构。 图 FIR滤波器的一般结构 因为 FIR滤波器的单位抽样响应是有限长的，所以它永远是稳定的。另外， 若对 h(n)提出一些约束条件，那么可以很容易地使 H(z)具有线性相位，这在信 号处理的很多领域是非常重要的。 FIR滤波