基于DSP的 FIR滤波器的设计Word文件下载.docx

1、关键词：数字滤波器；DSP；TMS320C54x；系统开发环境CCS；FIR滤波器ABSTRACTThe digital filter are one kind of equipment which by the digital multiplier unit, the accumulator and the time delay unit is composed. Its function is to inputs the discrete signal the numeric code to carry on operation processing, achieves the chang

2、e signal frequency spectrum the goal. Because the electronic accounting machine technology and large scale integrated circuits development, the digital filter available computer software has realized, also the available large scale integration digit hardware real-time realizes. Along with the DSP wi

3、despread application, also gradually enhances to the DSP design requirements, filters design has the absolute quality function to the DSP performances enhancement, realizes the principle in TMS320C54x under system development environment CCS to FIR filters DSP to carry on the discussion. Using the C

4、 language design corresponding filter, through the experiment simulation, the FIR filter which realizes on DSP completes the predetermined filter task from the input signal and output signals time domain and the frequency range curve observation.key word: Digital filter; DSP；System development envir

5、onment CCS；FIR filter第1章绪论1.1设计背景数字信号处理就是用数字信号处理器（DSP）来实现各种算法。由于具有精度高、灵活性强等优点，已广泛应用在数字图像处理、数字通信、数字音响、声纳、雷达等领域。而数字滤波技术又是进行数字信号处理的最基本手段之一。它是对数字输入信号进行运算，产生数字输出信号，以改善信号品质，提取有用信息，或者把组合在一起的多个信号分量分离开来为目的。数字滤波与模拟滤波相比有很多优点，它除了可避免模拟滤波器固有的电压漂移、温度漂移和噪声等问题外，还能满足滤波器对幅度和相位的严格要求。低通有限冲激响应滤波器（低通FIR滤波器）有其独特的优点，因为FIR系统

6、只有零点，因此，系统总是稳定的，而且容易实现线性相位和允许实现多通道滤波器。根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应（IIR）滤波器和有限冲激响应（FIR）滤波器。与IIR滤波器相比，FIR的实现是非递归的，总是稳定的；更重要的是，FIR滤波器在满足幅频响应要求的同时，可以获得严格的线性相位特性。因此，它在高保真的信号处理，如数字音频、图像处理、数据传输、生物医学等领域得到广泛应用。FIR滤波器的设计是数字信号处理技术的基础，也是DSP芯片的重要组成部分。FIR滤波器性能的好坏直接影响着DSP的运行速度和精度，对现代电子技术的发展起决定性作用，FIR滤波器的设计在以后的一段相当

7、长的时间里将持续主导DSP，而DSP随着电子技术的不断发展，应用领域愈来愈广泛1.2设计要求利用C语言在CCS环境中编写一个FIR滤波器程序，并能利用已设计好的滤波器对常用信号进行滤波处理。1.3设计思路简介在TMS320C54x系统开发环境CCS（Code Composer Studio）下对FIR滤波器的DSP实现原理进行讨论。通过MATLAB设计确定FIR滤波器系数。利用C语言设计相应的滤波器，通过实验仿真，从输入信号和输出信号的时域和频域曲线可看出在DSP上实现的FIR滤波器能完成预定的滤波任务。第2章 系统开发平台与环境2.1 CCS开发环境CCS提供了配置、建立、调试、跟踪和分析程

8、序的工具，它便于实时、嵌入式信号处理程序的编制和测试，它能够加速开发进程，提高工作效率。CCS提供了基本的代码生成工具，它们具有一系列的调试、分析能力。CCS支持如下图1.1所示的开发周期的所有阶段。图 2.12.2 SEED-DEC2812开发实验箱 SEED-DECxxxx系列嵌入式DSP开发板本着模块化、总线型、开放式、系列化的设计思想，采用统一的系统结构、模块结构和机械结构，以多种典型DSP处理器构成具有标准总线和相同物理尺寸的高性能嵌入式DSP开发板。SEED-DEC2812 嵌入式DSP开发板原理框图如图1.2所示：图 2.2第3 章FIR滤波器设计过程3.1 FIR滤波器设计总框

9、图FIR滤波器设计总框图如图3.1所示图3.1 FIR滤波器设计总框图3.2 FIR滤波器设计的原理设a i（i=0，1，2，N-1）为滤波器的冲激响应,输入信号为x（n）,则FIR滤波器的输入输出关系为：在数字信号处理应用中往往需要设计线性相位的滤波器，FIR滤波器在保证幅度特性满足技术要求的同时，很容易做到严格的线性相位特性。为了使滤波器满足线性相位条件，要求其单位脉冲响应为实序列，且满足偶对称或奇对称条件即或者 。这样，当N为偶数时，偶对称线性相位FIR滤波器的差分方程表达式为 由上可见，FIR滤波器不断对输入样本x（n）延时后，再做乘法累加算法，将滤波器结果y（n）输出。因此，FIR实

10、际上是一种乘法累加运算。对于线性相位FIR而言，利用线性相位FIR滤波器系数的对称特性，可以采用结构精简的FIR结构将乘法器数目减少一半。FIR滤波器的结构如图3.2所示：图3.2 FIR滤波器的结构 FIR滤波器的设计方法循环缓冲算法:对于N级的FIR滤波器，在数据存储器中开辟一个称之为滑窗的N个单元的缓冲区，滑窗中存放最新的N个输入样本。每次输入新的样本时，一新样本改写滑窗中的最老的数据，而滑窗中的其他数据不需要移动。利用片内BK（循环缓冲区长度）寄存器对滑窗进行间接寻址，环缓冲区地址首位相邻。下面，以N=5的FIR滤波器循环缓冲区为例，说明循环缓冲区中数据是如何寻址的。5级循环缓冲区的结

11、构如图3.3所示，顶部为低地址。图3.3 循环缓冲区的结构当第一次执行完之后，间接寻址的辅助寄存器指向x（n-4）。然后，从I/O口输入数据x（n+1），将原来存放x（n-4）的数据存储单元改写为x（n+1）。接着，进行第二次乘法累加运算，最后指向x（n-3）。然后从I/O口输入数据x（n+2） ，将原来存放 x（n-4）的数据存储器单元改写为x（n+2）。第4章 系统软件设计4.1 用MATLAB确定参数Matlab确定滤波器系数程序如下：Clear allb=1,1;a=1,0.2,-0.24; % 滤波器函数系数impulse=1,zeros（1,9）;h=filter（b,a,impu

12、lse）;subplot（3,2,1）;stem（h）; %绘制脉冲响应subplot（3,2,2）;impz（b,a,10）;h,w=freqz（b,a）; %绘制频率响应subplot（3,2,3）;plot（w,20*log10（abs（h）; %幅频曲线subplot（3,2,5）;plot（w,angle（h）*180/pi）;%相频曲线subplot（3,2,4）;zplane（b,a）;%零极点分布图pd=-unwrap（angle（h）./w;gd=grpdelay（b,a）;subplot（3,2,6）;plot（w,gd,-,w,pd,:）;图4.1Matlab确定滤波器系

13、数4.2 程序流程图程序流程图如图4.2所示图4.2 程序流程图4.3程序源代码4.3.1 FILTER.c代码：#include DSP28_Device.hfilter.hext_inf.hunsigned int DataLong=1024;unsigned int i,j;unsigned int k= 0;unsigned short Ad_data2048;short Ad_data_signed2048;short filter_result1024;unsigned int convcount = 0;volatile unsigned int adconvover =0;in

14、terrupt void ISRTimer2（void）;interrupt void ad（void）;void main（void） /*初始化系统*/ InitSysCtrl（）; /*关中断*/ DINT; IER = 0x0000; IFR = 0x0000; /*初始化PIE中断*/ InitPieCtrl（）; /*初始化PIE中断矢量表*/ InitPieVectTable（）; /初始化cputimer InitCpuTimers（）; /*设置中断服务程序入口地址*/ EALLOW; PieVectTable.TINT2 = &ISRTimer2; PieVectTable

15、.ADCINT = &ad; EDIS; /*开中断*/ IER |= M_INT1; /ADC中断 EINT; ERTM; /*设置CPU*/ ConfigCpuTimer（&CpuTimer2, 150, 22）; StartCpuTimer2（）; /*开中断*/ IER |= M_INT14; InitAdc（）; for（;） if （adconvover=1） /*fir滤波处理*/ for（i=0;i=16）; _nassert（n for（j=0;j=m） break; else y0=（long）xi+j*（long）hi; acc=acc+y0; *y+=（short）（a

16、ccs）;第5章 系统仿真5.1仿真设置1.fir_filter（Ad_data_signed,hfir,filter_result,ORDER_FIR,DataLong,ROUND_FIR）和adconvover=0处设置断点2.在Graph Property Dialog中设置Start Address和Display Type如图5.1、图5.2、图5.3、图5.4图 5.1图 5.2图 5.3图 5.43.设置实验板输入信号的参数，再运行程序得到仿真图。5.2 仿真图（1）标准矩形波，频率为300Hz，振幅为800mV。滤波前的波形如图5.5和图5.6，滤波后的波形如图5.7和图5.8

17、。图 5.5图 5.6图 5.7图 5.8（2）噪声三角波，频率为400Hz，振幅为800mV。滤波前的波形如图5.9和图5.10，滤波后的波形如图5.11和图5.12。图 5.9图 5.10图 5.11图 5.12总结通过两个星期的努力，完成了本设计的工作，在设计过程中，充分的利用图书馆的资料，利用网络查看了大量相关的论文，认真地阅读、理解，做下笔记，并将书中内容加以充分利用。通过这段时间的学习，发现了自己的很多不足，自己知识的很多漏洞，看到了自己的实践经验还是比较缺乏，理论联系实际的能力还急需提高。这次课程设计实现了一个简单的FIR滤波器的设计，在罗老师的精心指导和同组同学的帮助下较好的完

18、成了课程设计的任务并达到了理想的目的。在这里我对罗老师表示感谢，感谢她的耐心和热心。通过这次课程设计使我对DSP芯片有了更深一步的了解，既巩固了理论知识，又提高了实践能力。对我而言，知识上的收获重要，精神上的丰收更加可喜。这次设计让我知道了学无止境的道理。我们每一个人永远不能满足于现有的成就，人生就像在爬山，一座山峰的后面还有更高的山峰在等着你。挫折是一份财富，经历是一份拥有。这次课程设计必将成为我人生旅途上一个非常美好的回忆！参考文献1高西全,丁玉美.数字信号处理教程M.西安：西安电子出版社，2008年.2孙宗瀛,谢鸿林.TMS320C5xDSP原理设计与应用M.北京：清华大学出版社，2002年.3乔瑞萍,崔涛,张芳娟.TMS320C54xDSP原理及应用M.西安：西安电子科技大学出版社，2005年.4张雄伟. DSP芯片的原理与开发应用（第四版）M.北京：电子工业出版社，2008年.5郑红.TMS320C54XDSP应用系统设计M. 北京：北京航空航天大学出版社，2002年.