DSP原理及应用课程设计报告FIR数字滤波器应用

《DSP原理及应用课程设计报告FIR数字滤波器应用》由会员分享，可在线阅读，更多相关《DSP原理及应用课程设计报告FIR数字滤波器应用（16页珍藏版）》请在装配图网上搜索。

1、徐州师范大学物理与电子工程学院课程设计报告- I -*大大 学学 * 学学 院院本 科 生 课 程 设 计课课 程程 名名 称：称： dsp 原理及应用原理及应用 题题 目：目： FIR 数字滤波器应用 专专 业业 班班 级：级： 08 电信 学学 生生 姓姓 名：名： * 学学 生生 学学 号：号： * 日日 期：期： 2011.6.17 指指 导导 教教 师：师： * *教务部印制教务部印制- II -目目 录录目 录.II摘 要.IIIABSTRACT.III1 引言.41.1 研究背景.41.2 研究目的.41.3 研究意义.42 设计方案简述.521 系统方案设计流程图.52.2 设

2、计原理.52.2.1 FIR 滤波器的基本结构.63. FIR 滤波器的详细设计.73.1 方案选择.73.2 低通滤波器的 MATLAB 实现 .73.3 FIR 滤波器的 DSP 实现 .93.3.1 汇编语言设计.93.3.2 具体的操作步骤.114 实验结果及分析.124.1 实验步骤与内容.124.2 实验过程中出现的错误及解决的办法.134.3 CCS 程序运行后的各种输出结果 .145 总结.156 参考文献.16- III -摘 要数字信号处理器，也称 DSP 芯片，是针对数字信号处理需要而设计的一种具有特殊结构的处理器。它是现代电子技术的、相结合的产物。目前的 DSP 芯片以

3、其强大的数据处理功能在通信和其他信号处理领域得到广泛注意并已成为开发应用的热点技术。DSP 芯片是一种特别适合于进行数字信号处理运算的微处理器。主要应用是实时快速的实现各种数字信号处理算法，如卷积及各种变换等。其中利用 DSP 来实现数字滤波器就是很重要的一种应用，本文深入研究基于美国德州仪器公司(TI)TMS320C54XX DSP 芯片的滤波器系统软件实现方法，用窗口设计法实现 FIR 滤波器。 关键词：DSP；TMS320C5410；FIR 滤波器;自适应滤波器AbstractDigital signal processor, also known as DSP chips, digit

4、al signal processing required for the design of a special structure of the processor. It is a modern electronic technology, the product of the combination. The current DSP chip with its powerful data processing capabilities in the areas of communications and other signal processing has become widesp

5、read attention and development and application of a hot technology. DSP chip is particularly suitable for digital signal processing algorithms of the microprocessor. Main application is to achieve a variety of real-time fast digital signal processing algorithms such as convolution and a variety of t

6、ransformation. One advantage of the DSP implementation of digital filters is a very important application of this in-depth research is based on Texas Instruments (TI) TMS320C54XX DSP chip filter system software implementation, with the window design method of FIR filter.KeywordsKeywords：DSP; TMS320C

7、5410; FIR filter; adaptive filter - 4 -1 引言数字滤波器在数字通信、语音图象处理、谱分析、模式识别、自动控制等领域得到了广泛的应用。相对于模拟滤波器，数字滤波器没有漂移，能够处理低频信号，频率响应特性可做成非常接近于理想的特性，且精度可以达到很高，容易集成等。和 IIR 滤波器相比，FIR 滤波器可以得到严格的线性相位。而且FIR 滤波器可以用非递归方法实现，有限精度的计算不会产生振荡，同时由于量化舍入以及系数的不准确所引起的误差的影响比 IIR 滤波器要小得多。再者FIR 滤波器可采用 FFT 算法，在相同阶数下，运算速度可以快得多。1.1 研究背景数

8、字滤波器是 DSP 最基本的应用之一。在系统设计中，滤波器的好坏将直接影响系统的性能。使用 DSP 进行数字信号处理可以从噪声中提取信号，即对一个具有噪声和信号的混合源进行采样，然后经过数字滤波器来滤除噪声，从而提取有用信号。用 DSP 芯片实现数字滤波除了具有稳定性好、精确度高、不受外界影响等优点外，还具有灵活性好的特点，因而可通过可编程 DSP 芯片来实现数字滤波器的特性。1.2 研究目的通过课程设计，加深对 DSP 芯片 TMS320C54x 的结构、工作原理的理解，获得 DSP 应用技术的实际训练，掌握设计较复杂 DSP 系统的基本方法。通过使用汇编（或者 C）语言编写具有完整功能的

9、FIR 系统，使学生加深对所学知识的理解，进一步巩固汇编（或者 C）语言讲法规则。学会编制结构清晰、风格良好、数据结构适当的汇编（或者 C）语言程序，从而具备解决综合性实际问题的能力。1.3 研究意义数字滤波器分为有限冲激响应滤波器 FIR 和无限冲激响应滤波器 IIR。FIR滤波器属于经典滤波器，优点就是不存在系统极点，FIR 滤波器是绝对的稳定的系统，FIR 滤波器还确保了线性相位，在信号处理中占有极其重要的地位。数字滤波器一直以来就是数字信号处理器 DSP 最广为认知的应用，FIR 滤波器的单位响应 b（n）为有限长序列，若 b（n）为实数，且满足偶对称：b（n）=b(N-1-n)的条件

10、，称为对称 FIR 滤波器。系数对称 FIR 滤波器在数字信号处理中应用十分广泛，常用于相位失真要求较高的场合。如：通信系统：调制解调器综合业务数据网；希尔伯特变换器：要求输入输出信号正交；高保真音响系统音乐的相位失真必须减到最小，尽可能逼真地重现原来的声音等。- 5 -2 设计方案简述21 系统方案设计流程图 图 2.1 系统方案设计流程图2.2 设计原理数字滤波器的设计问题就是寻找一组系数 ai 和 bi，使得其性能在某种意义上逼近所要求的特性。如果在 s 平面上去逼近，就得到模拟滤波器，如果在 z平面上去逼近，则得到数字滤波器。数字滤波是将输入的信号序列，按规定的算法进行处理，从而得到所

11、期望的输出序列。一个线性位移不变系统的输出序列 y(n)和输入序列 x(n)之间的关系，应满足常系数线性差分方程： （2.3.1） x(n)为输入序列，y(n)为输出序列，ai、bi 为滤波器系数，N 为滤波器的阶数。使用 MATLAB 编写程序产生待滤波数据使用 MATLAB 滤波器的设计分析工具FDATOOL 设计 FIR 滤波器使用 CCS 编写程序实现滤波功能查看程序运行后的输入输出波形并比较滤波效果0 )()()(110ninyainxbnyMiiNii- 6 -2.2.1 FIR 滤波器的基本结构在式 2.1.1 中，若所有的 ai 均为 0，则得 FIR 滤波器的差分方程： （2

12、.3.2）对式（2.3.1）进行 z 变换，可得 FIR 滤波器的传递函数： （2.3.4）由此可得到 FIR 滤波器的结构如图 1 所示。FIR 滤波器的单位冲击响应 h(n)是一个有限长序列。若 h(n)为实数，且满足偶对称或奇对称的条件，即 h(n)=h(N-1-n)或 h(n)=-h(N-1-n)，则 FIR 滤波器具有线性相位特性。 )()(10Niiinxbny )()()(10NiiizbzXzYzH图 2.3.1 FIR 滤波器的结构如图- 7 -3. FIR 滤波器的详细设计3.1 方案选择在设计数字滤波器时通常采用 MATLAB 来进行辅助设计和仿真用 MATLAB设计 F

13、IR 数字滤波器的方法有很多种!限于篇幅!这里只介绍比较常用的窗函数设计方法最优化设计法及滤波器设计与分析工具 FDATOOL 设计法.(1)窗函数设计法窗函数法的设计核心是从给定的频率特性,通过加窗确定有限长单位脉冲响应序列 h（n）根据给定的滤波器技术指标，选择滤波器长度 N 和窗函数 w（n）使其具有最窄宽度的主瓣和最小的旁瓣。工程中常用的窗函数共有 6 种，即矩形窗巴特利特窗汉宁窗汉明窗布莱克曼窗和凯塞窗#以基于窗函数设计法的FIR1 和 FIR2 函数命令可以方便地得到滤波器的系数,并可以通过使用freqz（）和 plot（）函数绘制滤波器的幅频响应和相频响应。（2）最优化设计法MA

14、TLAB 中的 remez 函数可采用 Park-McClenllan 算法设计线性相位 FIR 滤波器,这种滤波器使期望频率响应与实际频率响应之间的最大误差最小，使其频率响应呈现出等波纹特性，因此有时也称为等波纹滤波器。这种方法的缺点是无法直接设定滤波器阻带衰减参数，只能通过调整参数进行多次试验来得到所需的滤波器。（3）FDATool 设计法FDATool 是 MATLAB 信号处理工具箱里专用的滤波器设计与分析工具。FDATool 界面为滤波器设计提供了一个交互式的设计环境，用户可以根据对幅值和零极点图的设置，设计几乎所有基本的常规滤波器。采用 FDATool 设计法是其界面的上半部分为特

15、性区，用来显示滤波器的各种特性；下半部分为参数设定区，用来设定滤波器的各种参数，包括滤波器的类型设计方法阶数频率设定等。通过对比以上 3 种方法,就会发现利用窗函数设计滤波器更为直观和方便。本文便是利用这种方法首先设计出 FIR 滤波器的系数,然后根据其结果再在 DSP上实现。3.2 低通滤波器的 MATLAB 实现用窗函数法设计 FIR 低通滤波器，其技术指标如下：kHz, 250fs20pfkHz，kHz，通带最大衰减，阻带最小衰减。30sf3dBpA 50dBsA 将其换算成数字域的性能指标如下：通带截止频率，通带最大衰0.16p减 ；阻带截止频率，阻带最小衰减 。3dBpA 0.238

16、s50dBsA 根据窗函数法的设计原则，海明窗（hamming）可提供大于 50dB 的衰减。- 8 -要求滤波器的过渡带为：,利用海明窗0.240.160.08sp设计的滤波器的过渡带 ，所以低通滤波器单位脉冲响应的长度为：8 / N ，取 N=101。881000.08N通过 FIR1 函数设计出滤波器的单位冲激响应序列，用它来作为在下一章中DSP 设计程序中的系数。 (a)幅度响应 (b)相位响应图 3.2.1 hamming 窗滤波器的幅度响应与相位响应- 9 -3.3 FIR 滤波器的 DSP 实现3.3.1 汇编语言设计所选取的 N=101，滤波器的算法为 y(n)= x(n)+

17、x(n-1)+ x(n-2)+ x(n-100)根据我们所选择使用的循环缓冲区法可以编写得到 FIR 滤波器的源程序如下： .title fir.asm .mmregs .def start .bss yn,1xn .usect xn,101hn .usect hn,101indata .usect indata,300outdata .usect outdata,256 .datafir_coff: .word -8*32768/10000,-32*32768/10000 .word 16*32768/10000,56*32768/10000 .word 36*32768/10000,-89

18、*32768/10000 .word -77*32768/10000,92*32768/10000 .word 188*32768/10000,-53*32768/10000 .word -289*32768/10000,-76*32768/10000 .word 423*32768/10000,322*32768/10000 .word -531*32768/10000,-851*32768/10000 .word 595*32768/10000,3113*32768/10000 .word 4298*32768/10000 .word 3113*32768/10000,595*32768/

19、10000 .word -851*32768/10000,-531*32768/10000 .word 322*32768/10000,423*32768/10000 .word -76*32768/10000,-289*32768/10000- 10 - .word -53*32768/10000,188*32768/10000 .word 92*32768/10000,-77*32768/10000 .word -89*32768/10000,36*32768/10000 .word 56*32768/10000,16*32768/10000 .word -32*32768/10000,-

20、8*32768/10000start: ssbx frct ;使用小数运算 stm #hn,ar1 ;系数首地址 rpt #101 ;将系数移入循 mvpd fir_coff,*ar1+ ;环缓冲区 stm #-1,ar0 stm #outdata,ar5 stm #indata,ar4 stm #xn+101,ar2 rpt #101 ;将输入数据移入 mvdd *ar4+,*ar2+0% ;循环缓冲区 stm #xn+101,ar2 ;指向第一个输入 stm #hn+101,ar3 stm #101,bk ;缓冲区大小 37 stm #255,brc ;块重复 256 次 stm #-1,

21、ar0 rptb loop-1 ;块程序重复大小 rptz a,# 101 ;计算一个输出 mac *ar2+0%,*ar3+0%,a sth a,*(yn) ;保存输出 mvkd *(yn),*ar5+ mvdd *ar4+,*ar2+0% ;读进一个输入loop: nop b loop .end - 11 -3.3.2 具体的操作步骤（1）打开 FDATOOL,根据滤波要求设置滤波器的类型、通带截止频率、指定阶数、采样频率等。设置完设计参数后再单击 Design Filter 按钮，生成滤波器系数。（2）把生成的滤波器系数传到目标 DSP。选择菜单 Targets-Export to Co

22、de Composer Studio(tm)IDE，打开 Export to C Header File 对话框，选择 C header file，然后指定变量名(滤波器阶数和系数向量)，输出数据类型可以选择浮点型或 32 b，16 b 整型等，再根据自己安装选择目标板板号和处理器号，单击OK，保存该头文件，同时需指定文件名(filtercoeffh)和路径(保存在c：timyprojectsfir 工程中)。（3）修改 CCS 汇编程序，删除掉数据前的所有文字，在其开头加上.data，第二行加上 coeff .word,在每行的前面都加上.word，并且把每行的最后的逗号去除掉。（4）编译汇

23、编程序，若有错误，则根据错误提示进行修改；若没有错误，则继续往下执行。（5）加载初始化 DATA 数据。运行程序，查看输入输出的波形，修改相应的参数进行调试- 12 -4 实验结果及分析4.1 实验步骤与内容MATLAB辅助DSP 实现FIR ,其总体过程为在DSP 中编写处理程序，在MATLAB中利用滤波器设计、分析工具( FDATOOL) ,根据指定的滤波器性能快速设计一个FIR ,再把滤波器系数以头文件形式导入CCS 中,头文件中MATLAB 辅助DSP 实现FIR 数字滤波器含滤波器阶数和系数数组,在MATLAB中调试、运行DSP 程序并显示、分析处理后的数据。使用该方法,便于采用汇编

24、语言来实现程序。头文件名不变,当MATLAB中设计的滤波器系数改变时,相应头文件中系数也改变,方便程序调试、仿真。利用 FDA TOOL 设计 FIR 滤波器的参数，操作步骤如图 4.1.1 所示。图 4.1.1 FDA TOOL 设计 FIR 滤波器的参数操作步骤图- 13 -图 4.1.2 已滤波的信号的时域波形 图 4.1.3 已滤波的信号的频域波形MATLAB 滤波器设计工具在完成 FIR 设计后,在 Export as 中选择 16bit 符号整数输出，然后单击 OK 按钮。将系数稍作修改，调整后拷贝到程序的系数初始化空间即可。2.在 CCS 中编写汇编程序，进行调试，实现带通滤波功

25、能。在 CCS IDE 中建立 LHM.pjt 工程,用汇编语言编写处理主程序 fir.asm。另根据板上的存储器配置方式,编写存储器配置文件 fir.cmd 文件,将 MATLAB 生成的 LHM.h 和input1.dat 文件拷贝到 LHM.pjt 工程文件夹下，进行编译、链接,生成可执行文件 LHM.out。 4.2 实验过程中出现的错误及解决的办法1.在 MATLAB 程序设计中，采样频率设置的过小，截止频率大于采样频率的一半，运行的时候图形出现错误。2. FIR.m 中的采样频率要和开始设置时的采样频率一致，否则结果会出现偏差。3.在 CCS 中未定义标号，程序运行错误。4.在.h

26、 文件中未把 coeff 顶格写，编译时出现错误。5.未将 fir.h 和 input.dat 文件考入工程文件中。6.间接寻址过程中，程序中丢失。7.在 View 的 Graph 中单击 Time/frequency 出现 graph property dialog框，未修改抽样点数。显示的图形出现差异- 14 -4.3 CCS 程序运行后的各种输出结果在 View 的 Graph 中单击 Time/frequency 出现 graph property dialog 框。将显示类型，图形名称，起始地址，抽样点数，数据类型等分别进行设置，输出各种波形。输出滤波前的信号时域波形。首先，设置如图

27、 4.1.4 所示，输出滤波前的信号时域波形。图 4.1.4 滤波前数据类型设置单击 OK 后生成如下图 4.1.5 波形。图 4.1.5 为滤波前的信号，波形很杂乱，从时域上很难看出信号的周期性。修改相应设置，输出如图 4.1.6 所示滤波前信号频谱。经过滤波后，时域波形和信号频谱分别如图 4.1.7 和图 4.1.8 所示。- 15 -图 4.1.5 滤波前时域波形 图 4.1.6 滤波前信号频谱图 4.1.7 滤波后时域波形 图 4.1.8 滤波后信号频5 总结在这次课程设计中，我了解到了数字滤波是信号处理技术中的重要部分，研究了数字滤波器的基本理论知识以及它实现方法。学习了数字滤波器的

28、结构、设计理论，掌握了各种数字滤波器的原理和特性，取得了很多收获。 首先，通过复习课本和实验，对 DSP 的基础知识又有了很大的巩固。其次，通过对用 Mathlab 实现 FIR 滤波器的设计，熟悉了 Mathlab 软件的一些相关的窗口函数以及相关功能的调用，如怎样实现窗函数的调用，怎样实现滤波等。 第二，增进了与同学的交流。平常除了上课，大家也比较少接触。而每次做课程设计的时候，大家总会互相探讨，发表自己的看法，帮忙解决遇到的问题，分享自己的心得，于是交流也就不由自主的变得相对频繁些了，感觉也就更亲切了。 第三，对 FIR 的滤波器性能和作用有了更深一层的了解。FIR 滤波器的应用十分广泛

29、，当今许多信号处理系统和图像处理系统等都要求信号具有线性相位特性。在这方面，FIR 滤波器有十分独特的优点，运用 MATLAB 语言，我们能够很容易的设计出具有严格线性相位的 FIR 滤波系统，以及比较容易的实现。此外，我们还知道了 FIR 滤波器是永远稳定的，因为 FIR 滤波器的冲激响应是有限长序列，其系统函数为一个多项式，它所含的的极点多为原点。在今后的工作和学习生活中，综合运用模拟电子、数字电子和 DSP 基本原理等课程中所学的理论知识去独立完成一个项目的设计。不仅要完善以上的MATLAB 设计方法，更要学会利用多种途径设计各种 DSP 最小系统，熟练汇编语言编程。- 16 -在整个设

30、计过程中，首先我要感谢的是老师在课程设计上给予我的指导、提供给我的支持和帮助，这是我能顺利完成这次报告的主要原因。其次，我要感谢帮助过我的同学，他们也为我解决了不少我不太明白的程序编写的难题。同时也感谢学院为我提供良好的做课程设计的环境。6 参考文献1（日）谷秋隆嗣著. 数字滤波器与信号处理. 科学出版社.2003.2 潘松，黄继业，王国栋. 现代 DSP 技术. 西安电子科技大学出版社. 20033 陈金鹰.DSP 技术及应用. 机械工业出版社.2004.64 孙宗瀛. DSP 原理设计与应用M. 北京：清华大学出版社，2002.5 汪安民. TMS320C54XX DSP 实用技术. 清华

31、大学出版社. 2002:5-156 TMS320C54x Code Composer Studio Tutorial. Texas Instruments Incorporated. 20007 薛年喜.Matlab 在数字信号处理中的应用M.北京：清华大学出版社，2003.8王立宁. Matlab 与通信仿真M. 北京: 人民邮电出版社, 2000. 276281.9余成波，杨菁，杨如民编著. 数字信号处理及 MATLAB 实现. 清华大学出版社.2005.10程佩青.数字滤波与快速傅里叶变换.北京:清华大学出版社,198811胡广书. 数字信号处理理论、算法与实现. 清华大学出版社. 199712彭启 .DSP 与实时数字信号处理.成都:电子科技大学出版社,1995