DSP中的FIR滤波器论文(2)

595*32768/10000,3113*32768/10000 4298*32768/10000 3113*32768/10000,595*32768/10000 -851*32768/10000,-531*32768/10000 322*32768/10000,423*32768/1000 0 -76*32768/10000,-289*32768/10000 -53*32768/10000,188*32768/10000 92*32768/10000,-77*32768/10000 -89*32768/10000,36*32768/10000 56*32768/10000,16*32768/10000 -32*32768/10000,-8*32768/10000

5.FIR滤波器的DSP实现

所选取的N=37，滤波器的算法为

y(n)=b0x(n)+b1x(n-1)+b2x(n-2)+···+b36x(n-36)

根据我们所选择使用的循环缓冲区法可以编写得到FIR滤波器的源程序如下：

.title "fir.asm" .mmregs .def start .bss yn,1

xn .usect "xn",37 hn .usect "hn",37 indata .usect "indata",300 outdata .usect "outdata",256 .data fir_coff:

.word -8*32768/10000,-32*32768/10000 .word 16*32768/10000,56*32768/10000 .word 36*32768/10000,-89*32768/10000 .word -77*32768/10000,92*32768/10000 .word 188*32768/10000,-53*32768/10000 .word -289*32768/10000,-76*32768/10000

介绍了FIR的滤波器的实现方法

.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/10000 .word -851*32768/10000,-531*32768/10000 .word 322*32768/10000,423*32768/10000 .word -76*32768/10000,-289*32768/10000 .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,-8*32768/10000 .text start:

ssbx frct ;使用小数运算 stm #hn,ar1 ;系数首地址 rpt #36 ;将系数移入循

mvpd fir_coff,*ar1+ ;环缓冲区 stm #-1,ar0 stm #outdata,ar5 stm #indata,ar4 stm #xn+36,ar2

rpt #36 ;将输入数据移入 mvdd *ar4+,*ar2+0% ;循环缓冲区 stm #xn+36,ar2 ;指向第一个输入 stm #hn+36,ar3

stm #37,bk ;缓冲区大小37 stm #255,brc ;块重复256次

介绍了FIR的滤波器的实现方法

stm #-1,ar0

rptb loop-1 ;块程序重复大小 rptz a,#36 ;计算一个输出 mac *ar2+0%,*ar3+0%,a sth a,*(yn) ;保存输出 mvkd *(yn),*ar5+

mvdd *ar4+,*ar2+0% ;读进一个输入 loop: nop

b loop .end

6.FIR滤波器的结果检验

(1)点击View—>Graph—>Time/Frequency...

如图设置好参数，如图5所示。点击OK，可得到输入信号的时域波形。

图5

（2）重复上述操作，改变Display Type，得到输入信号频域波形，输出信号时域波形，输出信号频域波形。如图6所示。

介绍了FIR的滤波器的实现方法

图6

可知能达到设计要求，将输入的较高频2500hz的信号滤除，让低频1000hz的信号通过。

7.调试问题

1）实验时，未做初始化设置直接进入，没有装入gel文件及相关参数，导致输入代码编译链接后不能载入程序。

解决方法：初始化设置软件仿真cpu和硬件仿真cpu，并设置其中所用的gel文件，保存并退出启动dsp选择忽略进入编辑编译初始化界面，由软件仿真cpu进入编辑编译工具。 2）在解决问题1后，载入程序，载入数据，并设置载入数据的初始地址和长度，运行程序，发现程序一直卡死在某个地方，无法显示正确的时域图和频域图。

解决方法：换了一台电脑运行，发现问题已被解决，可能原因是自己电脑安装时装少了一些文件。

介绍了FIR的滤波器的实现方法

8.心得体会

刚开始面对大作业的时候，我们都还是觉得很生疏，不知从哪里入手，尽管做过实

验。对于自己选的题目，具体需要做什么，能不能做出来，我们并不清楚。但当时间一天天过去后，我们开始认真起来了。从复习课本，再到复习实验，相互探讨编程思路，上网搜索资料，调试实验等，我们都在努力的自学着。通过这次FIR滤波器的设计，我们都有不同程度的收获。

首先，是学习上的巩固。一方面，通过复习课本和实验，对DSP的基础知识又有了很大的巩固。其次，通过对用Mathlab实现FIR滤波器的设计，熟悉了matlab软件的一些相关的窗口函数以及相关功能的调用，如怎样实现窗函数的调用，怎样实现滤波，等等。

第二，增进了与同学的交流。平常除了上课，大家也比较少接触。而每次做课程设计的时候，大家总会互相探讨，发表自己的看法，帮忙解决遇到的问题，分享自己的心得，于是交流也就不由自主的变得相对频繁些了，感觉也就更亲切了。

第三，对FIR的滤波器性能和作用有了更深一层的了解。FIR滤波器的应用十分广泛，当今许多信号处理系统和图像处理系统等都要求信号具有线性相位特性。在这方面，FIR滤波器有十分独特的优点，运用mathlab语言，我们能够很容易的设计出具有严格线性相位的FIR滤波系统，以及比较容易的实现。此外，我们还知道了FIR滤波器是永远稳定的，因为FIR滤波器的冲激响应是有限长序列，其系统函数为一个多项式，它所含的的极点多为原点。

介绍了FIR的滤波器的实现方法

参 考 文 献

【1】. 程佩青，数字信号处理，清华大学出版社 【2】.邹彦，唐东，DSP原理与应用，电子工艺出版社

【3】.王宏，MATLAB6.5在信号处理中的应用，清华大学出版社 【4】.谷萩隆嗣，数字滤波器与信号处理，科学出版社