基于MATLAB的FIR滤波器设计6.docVIP

EDA课程设计报告 滤波器设计参数：根据要求，要设计一个输入8位，输出8位的17阶线性相位FIR滤波器，所以采用图2(a)的方式，其中输入信号范围为：[±99，0，0，0， ±70，0，0，0， ±99，0，0，0， ±70，…]，此滤波器 Fs为44kHz,Fc为10.4kHz。 （一）FIR数字滤波器理论简述 有限冲激响应（FIR）数字滤波器和无限冲激响应（IIR）数字滤波器广泛应用于数字信号处理系统中。IIR数字滤波器方便简单，但它相位的非线性，要求采用全通网络进行相位校正，且稳定性难以保障。FIR滤波器具有很好的线性相位特性，使得它越来越受到广泛的重视。 有限冲击响应（FIR）滤波器的特点： 1 既具有严格的线性相位，又具有任意的幅度； 2 FIR滤波器的单位抽样响应是有限长的，因而滤波器性能稳定； 3只要经过一定的延时，任何非因果有限长序列都能变成因果的有限长序列，因而能用因果系统来实现； 4 FIR滤波器由于单位冲击响应是有限长的，因而可用快速傅里叶变换(FFT)算法来实现过滤信号，可大大提高运算效率。 5 FIR也有利于对数字信号的处理，便于编程，用于计算的时延也小，这对实时的信号处理很重要。 6 FIR滤波器比较大的缺点就是阶次相对于IIR滤波器来说要大很多。 FIR数字滤波器是一个线性时不变系统（LTI），N阶因果有限冲激响应滤波器可以用传输函数H（z）来描述， 在时域中，上述有限冲激响应滤波器的输入输出关系如下： 其中，x[n]和y[n]分别是输入和输出序列。 N阶有限冲激响应滤波器要用N＋1个系数描述，通常要用N+1个乘法器和N个两输入加法器来实现。乘法器的系数正好是传函数的系数，因此这种结构称为直接型结构，可通过式（1.2）来实现，如图1。 FIR滤波器具有对称结构，为线性相位滤波器： 。 (1.3) 这种对称性，可使得乘法器数量减半：对n价滤波器，当n为偶数时，乘法器的个数为n/2个；当n为奇数时，乘法器的个数为(n+1)/2个。在电路实现中，乘法器占用的逻辑单元数较多。乘法器的增加，意味着电路成本增加，另外对电路的工作速度也有影响。 N阶线性相位的因果FIR系统的单位冲激响应滤波器可用对称冲激响应 或者反对称冲激响应 来描述。 具有对称冲激响应的FIR传输函数的冲激响应可写成如下形式： 当N为偶数时 当N为奇数时 则FIR线性相位系统的结构可转化成如图2（a）和图2（b）所示。 图2（a） N为奇数 图2（b） N为偶数 （二）设计方案 随着数字技术日益广泛的应用，以现场可编程门阵列（FPGA）为代表的ASIC器件得到了迅速普及和发展，器件集成度和速度都在高速长。FPGA既具有门阵列的高逻辑密度和高可靠性，又具有可编码逻辑器件的用户可编程特性，可以分布式算法是一种以实现乘加运算为目的的运算方法。它与传统算法实现乘加运算的不同在于执行部分积运算的先后顺序不同。简单地说，分布式算法在完成乘加功能时是通过将各输入数据每一对应位产生的部分积预先进相加形成相应部分积，然后在对各部门积进行累加形成最终结果，而传统算法是等到所有乘积产生之后再进行相加来完成乘加运算的。与传统算法相比，分布式算法可极大FPGA有着规整的内部逻辑块阵列和丰富的连线资源，特别适合细粒度和高并行度结构特点的数字信号处理任务，如FIR、FFT等。对于线性相位因果FIR滤波器，它的系列具有中心对称特性，即h(i)=±h(N-1-i)。令s(i)=x(i) ±x(N-1-i)，对于偶对称，代入式（1）可得： 根据要求，要设计一个输入8位，输出8位的17阶线性相位FIR滤波器，所以采用图2(a)的方式，其中输入信号范围为：[±99，0，0，0， ±70，0，0，0， ±99，0，0，0， ±70，…]，此滤波器 Fs为44kHz,Fc为10.4kHz。MATLAB设计计算滤波器系数过程如下： FIR滤波器参数设置，因为是17阶，所以Specify order处填16，h(0)=0. FIR滤波器的幅频响应 FIR滤波器的相频响应 FIR滤波器的冲激响应 FIR滤波器系数 Matable 仿真图 仿真结果： 仿真结果： 对FIR滤波器的系数进行调整，整数化 可得FIR滤波器的参数为[-12 -19 13 30 -14 -54 14 170 255 170 14 -54 -14 30 13 -19 -12] 根据以上所说