基于TMS320F2812截止頻率为2kHz低通滤波的设计与实现.docVIP

截止频率为2KHz的低通滤波器设计实现 滤波器是一种用来消除干扰杂讯的器件，将输入或输出经过过滤而得到纯净的直流电。对特定频率的频点或该频点以外的频率进行有效滤除的电路，就是滤波器，其功能就是得到一个特定频率或消除一个特定频率。 中心频率（Center Frequency）：滤波器通带的中心频率f0，一般取f0=（f1+f2）/2，f1、f2为带通或带阻滤波器左、右相对下降1dB或3dB边频点。窄带滤波器常以插损最小点为中心频率计算通带带宽。 截止频率（Cutoff Frequency）：指低通滤波器的通带右边频点及高通滤波器的通带左边频点。通常以1dB或3dB相对损耗点来标准定义。相对损耗的参考基准为：低通以DC处插损为基准，高通则以未出现寄生阻带的足够高通带频率处插损为基准。 通带带宽（BWxdB）：（下图）指需要通过的频谱宽度，BWxdB=（f2-f1）。f1、f2为以中心频率f0处插入损耗为基准，下降X（dB）处对应的左、右边频点。通常用X=3、1、0.5 即BW3dB、BW1dB、BW0.5dB 表征滤波器通带带宽参数。分数带宽（fractional bandwidth）=BW3dB/f0×100[%]，也常用来表征滤波器通带带宽。 （1） 式(1)的系统的差分方程表达式为: （2） 式中，y( n) 为输出信号，即经过滤波之后的信号; N为滤波器阶数，式( 2) 就是线性移不变系统的卷积和公式，即为FIR 滤波器要实现的差分方程。 提取滤波器参数： CCS环境下的FIR滤波器的设计及软件仿真： 流程图： C语言实现：从上述分析可知，FIR 滤波器的算法实际上是一个乘累加运算过程。在这个累加运算中，主要涉及到三个要素: 输入信号x( n) 、冲击响应即滤波器系数h( n) 和输出信号y( n) 。 基于Matlab 的FIR 数字滤波器的设计方法方便快捷，其主要有基于FDATool 的FIR 滤波器的设计和基于Matlab 内建函数的FIR 滤波器的设计。 以下设计一个低通滤波器，其通带截止频率fp =2000Hz，阻带起始频率fst = 3000Hz，抽样频fs =15000Hz，阻带衰减不少于－50dB。这里采用窗口函数设计法设计FIR 低通滤波器，用Matlab 内建函数实现h( n) 的计算。由于海明窗的阻带最小衰减为－53dB，且可以提供较小的过渡带宽，所以选择海明窗［3］。编写Matlab 程序，求出过度带宽，由过度带宽确定窗口长度N，求出低通滤波器的截止频率，进而求出低通滤波器的系数h(n) 。Matlab 程序如下: fp = 2000; fst = 3000; fs = 15000; wp = 2* fp /fs; ws = 2* fst /fs; delt = ws － wp; N0 = ceil( 6.6 /delt) ; N = N0 + mod( N0 + 1,2) ; n = N－1; wn = ( ws + wp) /2; h = fir1( n,wn) ; n = 0: N－1,hn = round( h* 32768),stem( n,h,. ) , hw = fft( h,512) ,w =[0: 511]/512; plot( w,20* log10( abs( hw) ) ) ,grid; 将得到的33 个系数转化为Q15 格式，即分别乘以32768，转化后的低通滤波器系数如下:h[]= { 31，62，50，－39，－177，－218，0，421，668，289，－730，－1655，－1323，988，4773，8341，9804，8341，4773，988，－1323，－1655，－730，289，668，421，0，－218，－177，－39，50，62，31} CCS 是供用户开发和调试DSP和MCU程序的集成开发软件。Code Composer StudioTM IDE 提供强健、成熟的核心功能与简便易用的配置和图形可视化工具，使系统设计更快。实现FIR数字低通滤波器的程序流程如图。 由于CCS 是采用工程( projiect) 来集中管理应用文档的，所以必须将它们全部添加到工程中。需要注意的是，采用C 代码写的源程序需要将C 语言的标准支持库文件rts．lib 添加到工程中，该库放在CCS 的c \ti \c2000 \cgtools \lib 下。 基于CCS 的图形可视化工具，可以分别查看输入信号x( n) 、输出信号y( n) 的时域波形和频域波形，混频输入信号波形图。通过图形对比，所设计的数字低通FIR 滤波器滤掉了输入信号的高频分量而只