基于DSP的FIR滤波器的设计报告.docx

1、基于DSP的FIR滤波器的设计报告1概述 -0-1数字滤波器的研究现状 -0-2系统硬件设计 -1 -2.1系统设计方案 -1-2.2主控模块 -3-3系统软件设计 -4-3软件整体设计 -4-3.2测量模块 -7 -4系统测试与总结 -8-4系统测试 -8-4.2总结体会 -10 -.参考文献 -11 -附录 -12 -1概述1.1数字滤波器的研究现状20世纪60年代起山于计算机技术、集成工艺和材料工业的发展，滤波器发 展上了一个新台阶，并且朝着低功耗、高精度、小体积、多功能、稳定可靠和价 廉方向努力，其中小体积、多功能、高精度、稳定可靠成为70年代以后的主攻 方向。导致RC有源滤波器、数字

2、滤波器、开关电容滤波器和电荷转移器等各 种滤波器的飞速发展，到70年代后期，上述儿种滤波器的单片集成已被研制出 来并得到应用。80年代，致力于各类新型滤波器的研究，努力提高性能并逐渐 扩大应用范围。90年代至现在主要致力于把各类滤波器应用于各类产品的开发 和研制。当然，对滤波器本身的研究仍在不断进，在国内外的研究中，设计FIR 滤波器所涉及的乘法运算方式有：并行乘法、位串行乘法和采用分布式算法的乘 法。并行乘法运行速度快，但占用的硬件资源极大。如果滤波器的阶数增加，乘 法器位数也将变大，硬件规模将变得十分庞大。位串行乘法器的实现方法主要是 通过对乘法运算进行分解，用加法器来完成乘法的功能，也即

3、无乘法操作的乘法 器。位串行乘法器使得乘法器的硬件规模达到了最省，但是山于是审行运算，使 得它的运算周期过长，运算速度与硬件规模综合考虑时不是最优的。1.2论文主要完成的工作本文主要应用软件设汁FIR数字滤波器，并对所设计的滤波器进行仿 真。应用DSP集成开发环境一CCS调试汇编程序，用TMS320C5402来实现了 FIR数字滤波。具体内容有：(1)研究数字滤波器的基本理论和实现方法，从原理上理解和分析数字 滤波器，并通过对数字滤波器的结构和设讣方法的学习，为设计数字滤波器奠定 理论基础。(2)研究FIR数字滤波器的设计方法，及如何用滤波器设计工具设计各 种类型的数字滤波器，并进行仿真和分析

4、。(3)对数字滤波器的DSP实现进行研究。2系统硬件设计2.1系统设计方案FIR滤波器的设讣方法主要有窗函数设计法和频率抽样设计法等，其中窗函 数设讣法是最基本的设讣方法。在设IIFIR滤波器中，一个最重要的讣算就是加 窗，釆用矩形窗是最直接和简便的方法，但采用矩形窗存在较大的Gibbis效应, 且矩形窗的笫一旁瓣与主瓣相比仅衰减l3dB,因此实际设计中一般采用其他窗 函数。主要介绍儿种常用的窗函数和频率抽样设汁法等。利用窗函数设计FIR（一）窗函数法的基本思想窗函数设讣的基本思想是要选取某一种合适的理想频率选择性滤波器，然后 将它的脉冲响应截断以得到一个线性相位和因果的FIR滤波器。因此这种

5、方法的 重点在于选择某种合适的窗函数和一种理想滤波器。对于给定的滤波器技术指 标，选择滤波器长度和具有最窄主瓣宽度和尽可能小的旁瓣衰减的某个窗函数。（二）儿种常用的窗函数工程中比较常用的窗函数有：矩形窗函数、三角窗函数、汉宇窗函数、海明 窗函数、布莱克曼窗函数和凯塞窗函数。窗函数的选择原则是：（1） 具有较低的旁瓣幅度，尤其是第一旁瓣幅度：（2） 旁瓣幅度下降速度要大，以利增加阻带衰减；（3） 主瓣的宽度要窄，以获得较陡的过渡带。通常上述三点很难同时满足。当选用主瓣宽度较窄时，虽然得到较陡的过渡 带，但通带和阻带的波动明显增加；当选用最小的旁瓣幅度时，虽能得到匀滑的 幅度响应和较小的阻带波动，

6、但过渡带加宽。因此，实际选用的窗函数往往是它 们的折衷。在保证主瓣宽度达到一定要求的条件下，适当牺牲主瓣宽度来换取旁 瓣波动的减少。表2.1儿种常用的窗函数对比减窗函数旁瓣峰值衰减过渡带阻带最小衰矩形窗-134%-21三角形-27%-25汉宁窗-31%-44海明窗-418%-53布莱克曼窗-57,2%-74凯瑟窗-57,0%-80比较可得：凯塞窗可提供变化的过渡带宽，通过改变B的值可达到最陡的过渡带:凯塞窗具有与海明窗相匹敌的特性，通过调整p的值，可将凯塞窗完全等价 海明窗；凯塞窗最大旁瓣值比主瓣约低80dB,在所有的窗函数中旁瓣抑制度最高2。2.2.2用频率抽样法设计FIR滤波器所谓频率抽样

7、法就是从频域出发，根据频域的釆样定理，对给定的理想滤波 器的频域响应进行等间隔采样。乞(討)|“ W (22)其中k=0丄2 N-1把乩当作待设计的滤波器频率响应的采样值H伙)，通过下式可以求岀滤 波器的系统函数H(z)和频率响应H(ejw):(2.3)(2.4)心加(5)冷黔HU =H伙)(w = k) (2-3)山于频谱的有限个采样值恢复出来的频率响应实际上是对理想频率响应的 逼近，因此，这种方法必然有一定的逼近误差。若被逼近的频率响应比较平滑, 则各采样点之间的逼近误差较小；反之，则逼近误差较大。2.2主控模块TMS320C5402的主要特性：TMS320C5402是TI公司于1996年

8、推出的一种定点DSP芯片，采用先进 的修正哈佛结构和8总线结构，使处理器的性能大大提高。其独立的程序和数据 总线，允许同时访问程序存储器和数据存储器，实现高速并行操作。例如，可以 在一条指令中，同时执行3次读操作和1次写操作o TMS320C5402的运行速度 为40MIPS (每秒执行40百万条指令)。指令周期为25nso此外，还可以在数 据总线与程序总线之间相互传送数据。从而使处理器具有在单个期内同时执行算 术运算，逻辑运算，位移操作，乘法累加运算以及访问程序和数据存储器的强大 功能。包括:CPU,总线，存储器，在片外围电路，串行口，以及外部总线接口等部 分组成。1 CPU40位算数逻辑运

9、算单元(ALU),包括1个40位桶形移位寄存器和2个独 立的40位累加器。口 17位X 17位并行乘法器，与40位专用加法器相连，用 于非流水线式单周期乘法/累加(MAC)运算。比较，选择，存储单元(CSSU),用 于加法/比较选择。双地址生成器，包括8个辅助寄存器和2个辅助寄存器算数 运算单元(ARAU ) o2、总线结构程序总线(PB)传送取自程序存储器的指令代码和立即数。3条数据总线 (CB, DB和EB)将内部各单元连接在一起。4条地址总线(PAB, CAB, DAB和 EAB )传送执行指令所需的地址。3存储器192K字可寻址存储空间(64K字程序存储器，64K字数据存储器，以及 6

10、4K字I/O空间)。片内2K字ROMo片内10K字双存取RANL4在片外围电路用来周期产生中断的定时器。8位主机接口( HPI),用来与主设备或主处理 器接口。一个自动缓冲串口( BSP)和时分多路串行口( TPM),都可用作标准同 步串口。3系统软件设计3.1软件整体设计FIR滤波器的实现结构有以下儿种形式(1)直接型结构该种结构实现时需N个乘法器、N-1个延迟器和N-1个加法器。结构图如下：图3.1直接型结构图(2)级联型结构FIR系统函数H(z)是Z1的N-1次多项式，其分解为实系数二阶因子的乘积 形式为：严CH(z)= ( P Ok+P IkZ1+ 3 2kZ2)k=根据上式就可画出级

11、联结构如下图。其中Nc为N/2的最大整数。该结构 即是把H(z)的共辄零点或两个单个零点组成基本二阶节，H(z)为基本二阶因子的 子系统函数积。图3. 2级联型结构图从图3.2中可看出，每个二阶因子用直接型实现，每一节控制一对零点，便 于传输零点的控制。但系数B lk的个数比h(n)的个数多，其乘法次数比直接型 多。(3)线性相位型结构线性相位特性是指滤波器对不同频率的正弦波所产生的相移和正弦波的频 率成直线关系。在滤波器通带内的信号通过滤波器后，除了山相频特性的斜率决 定的延迟外，可以不失真地保留通带以内的全部信号。线性相位的因果的FIR系统的单位取样响应具有如下特性h(n)二土 h(NTn

12、)当N为偶数时，式变为H二乞 h(n) z+z 1 n，71 = 0当N为奇数时，此式为匕 n-H(z)= 土 肛门才+严于勺+吐匚！)异”=() 2由上面两个式子可知:实现直接形式的结构网络只需要N/2(N为偶数)或(N十1)/2次乘法(N为奇数)，利用h(n)的对称比直接卷积形式少用一半的乘法器。图3.3 N为偶数，线性相位结构图3. 4 N为奇数，线性相位结构(4)频率取样结构系统函数H(z)在单位圆上作N等分取样的值得到h(n)的离散傅立叶变换值H(k): nt 伽H(k)=H(|y/ )=IH(k)ld， =h(k) Jn=0根据内插公式，得1 z H 伙)H(z)(l/) -7 -

13、N A=0 WN Z-2n.上式即为频率取样结构的系统函数其中，H(k)为单位圆上 的频率取样值。可见，该系统是由一个子FIR系统(1-Z-N)和一个子IIR系统图3. 5频率取样结构图FIR滤波器的实现方法数字滤波器的实现方法一般有以下儿种用计算机软件实现软件实现方法即是在通用的微型计算机上用软件实现。利用计算机的存储 器、运算器和控制器把滤波所要完成的运算编成程序通过计算机来执行，软件可 山使用者自己编写，也可以使用现成的。但是这种方法速度慢，难以对信号进行 实时处理，虽然可以用快速傅立叶变换算法来加快计算速度，但要达到实时处理 要付出很高的代价。(2)采用DSP (Digital Sig

14、nal Processing)处理器来实现DSP处理器是专为数字信号处理而设计的，它的主要数字运算单元是一个乘 累加器，能够在一个机器周期内完成一次乘累加运算，配有适合于信号处理的指 令，具备独特的循环寻址和倒序寻址能力。这些特点都非常适合数字信号处理中 的滤波器设计的有效实现，并且它速度快，成本低。用DSP芯片实现数字滤波除 了具有稳定性好、精确度高、不受环境影响外，还具有灵活性好的特点。(3)采用固定功能的专用信号处理器实现专用信号处理器采用专用集成电路ASIC来实现，适用于过程固定而乂追求高 速的信号处理任务，是以指定的算法来确定它的结构，使用各种随机逻辑器件组 成的信号处理器。它们体积

15、小、保密性好，具有极高的性能，然而灵活性差。(4)用FPGA等可编程器件实现使用相关开发工具和VHDL等硬件开发语言，通过软件编程用硬件实现特定 的数字滤波算法。这一方法山于具有通用性的特点并可以实现算法的并行运算， 无论是作为独立的数字信号处理，还是作为DSP芯片的协作处理器都是比较活跃 的研究领域。可以采用DSP或FPGA来实现硬件电路。3.2测量模块同出迪瞥伫勺频i制CCS算法流程图4系统测试与总结4.1系统测试在 View 的 Graph 中单击 Tiine/frequency 出现 graph property dialog 框。将显示类型，图形名称，起始地址，抽样点数，数据类型等分

16、别进行设置， 输出各种波形。输出滤波前的信号时域波形。首先，设置如图4. 1所示，输出 滤波前的信号时域波形。图 4. 1 Graph property dialog 设置框图单击OK后生成如下图4.2波形。图4. 2滤波前信号波形图图4.2为滤波前的信号，波形很杂乱，从时域上很难看出信号的周期性。修 改相应设置，输出如图4.3所示滤波前信号频谱。图4. 3滤波询信号频谱波形图经过滤波后，时域波形频谱波形如图4.4和4.5所示。图4. 4滤波后信号时域波形图图4. 5滤波后信号频谱图4.2总结体会通过本次课程设计，让我了解了数字滤波是信号处理技术中的重要部分，研 究了数字滤波器的基本理论知识以

17、及它实现方法。学习了数字滤波器的结构、设 计理论，掌握了各种数字滤波器的原理和特性。并且复习了 MATLAB的相关知识, 并且利用MATLAB与CCS设计了低通滤波器以及高通滤波器，还研究了如何在定 点DSP中实现数字滤波器的算法，掌握了 CCS环境下的程序开发方法、调试工具 的使用及优化级别的选择等。研究了 MATLAB环境下FIR数字滤波器的设计方法, 利用MATLAB软件编程实现FIR滤波器设计。DSP作为一门新兴学科，越来越引 起人们的关注，LI前已广泛应用在各个领域。所以在今后的工作和学习生活中， 综合运用模拟电子、数字电子和DSP基本原理等课程中所学的理论知识去独立完 成一个项LI

18、的设计。不仅要完善以上的MATLAB设计方法，更要学会利用多种途 径设讣各种DSP最小系统，熟练C语言编程，用DSP芯片实现低通FIR数字滤波, 除了具有准确度高、不受环境影响等优点外。通过对本课题的研究，关于数字滤 波器理论知识得到了丰富和提高，熟悉了课题研究的整体思路和方法，认识到进 行课题研究应具备的素质参考文献1王念旭.DSP基础与应用系统设计.北京：北京航空航天大学出版社, 20012张雄伟.DSP集成开发与应用实例.北京：电子工业出版社，20023张雄伟.DSP芯片的原理与开发应用（第二版）.北京：电子工业岀版社,20004彭宗启.DSP技术原理及应用.北京：电子工业出版社，199

19、85朱铭铭.基于模糊的DSP系统设计.北京：电子工业出版社，20036李真芳，苏涛等.DSP程序开发一 MATLAB调试及直接标代码生成M. 西安:西安电子科技大学出版社，2003附录用线性缓冲区实现FIR滤波器程序清单: title“FIR1. ASM” mmregs.defstartX usect“X”，8PAOset0PAIset1 dataCOEF： word1*32768/10 word 2*32768/10.word -4*32768/10 word 3*32768/10word-4*32768/10word2*32768/10word1*32768/10 textstart:SS

20、BXFRCTSTM#x+7, AR2STM#6, AROLD#x+l,DPPORTRPAI, x+lFIRl:RPTZA, #6MACD*AR2-, COEF, ASTHA, *AR2PORTW*AR2+, PAOBDFIR1PORTR PAI, *AR2+0end用循环缓冲区实现FIR滤波器程序清单:.title “FIR2.ASM” mmregsdefbssstarty, 1xn usect“xn” ,7bO usect“bO” ,7PAOset0PAIset1 datatable: word1*32768/10 word2*32768/10 word3*32768/10 word4*32

21、768/10 word5*32768/10 word6*32768/10 word7*32768/10 textstart:SSBX FRCTSTM#b0, ARIRPT#6MVPDtable, *AR1+STM#xn+6, AR2STM#b0+6, AR3STM#7,BKSTM#-l, AROLDttxn, DPPORTR PAI, xnFIR2:RPTZ A, #6MAC *AR2+0%, *AR3+0%, ASTH A, yPORTW y, PAOBD FIR2PORTR PAI, *AR2+0%end汇编程序清单lhmh文件内容： datacoeff word-85, -64, -61

22、, -36, & 62, 110, 136, 131 word96, 42, -11, -44, -44, -11, 39, 85, 106 word8& 37, -29, -83, -101, -73, -7, 70, 124wordword131, 82, 5, -96, -131, -144, 72, 37, 140191, 162,61, -76, -191, -231, -174, -35, 132 word256, 27& 180, -7, -20& -337, -330, -176, 72wordword316, 445, 392, 156, -17& -475, -599, -

23、470, -108361, 745, 857, 594, -6, -748, -1336, -1456, -893 word386, 2192, 4154, 5816, 676& 676& 5816, 4154, 2192 word386, -893, -1456, -1336, -748, -6, 594, 857, 745 word361, -10& -470, -599, -475, -17& 156, 392, 445 word316, 72, -176, -330, -337, -20& -7, 180, 278 word256, 132, -35, -174, -231, -191

24、, -76, 61, 162 word191, 140, 37, 一72, -144, 151, -96, 5, 82 word131, 124, 70, 7, -73, -101, -83, -29, 37 word8& 106, 85, 39, -11, -44, -44, -11, 42 word96, 131, 136, 110, 62, & -36, -61, -64 wordFir. asm 内容:-85 title fir.asmmmregsglobal _c_int00ORDER set 154;The Leng th of Input;Get coefs from coef

25、hD_LEN . set 1024Dataxn usect xn, (ORDER-1)a0 .usect d0, (ORDER-1) input usect input, D_LEN output usect output, D_LENcopy HM. h textasgasgARO, FIRINDEXAR2, FIR.DATA.asgAR3, FIR_COEF.asgAR5, DATA.IN.asg_c_int00：AR6, DATA_OUTSSBXFRCTSTMFIR_C0EF(AR3)RPT#aO, FIR_COEF ; Copy aO(coefs) toSORDER-lMVPD#coe

26、ff,*FIR_COEF+STM#1, FIR_INDEXSTMFIR_DATA(AR2)RPTZ#xn, FIR_DATA ; Copy xn(data) toA,#ORDER-1STLA, *FIR_DATA+STM#(xn+ORDER-1), FIR_DATASTM# (aO+ORDER-1), FIR_COEFSTM#input, DATA_IN ; Get data from inputSTM#output, DATA_OUT ; Write data to outputSTM#D_LEN-1, BRCRPTBDnex tTSTMSORDER, BKLD*DATA_I+, AFIR:

27、 STLA, *FIR_DATA+%RPTZA, (ORDER-1)MAC*FIR_DATA+0%, *FIR_COEF+O%, A ; FIRSSTHA, *DATA_OUT+nextFIR END:BendFIR.ENDFIR. m的程序如下： fsMOOOO; %/采样 HzN二1024 %数据个数T二1/fs; %采样周期n二 O:NT;df=n*(fs/N)%待滤波信号波形xin=randn(l,1024)figure (1)plot (xin)%待滤波信号频谱xinff=abs(fft(xin);figure (2)plot (df, xinff)%滤波后信号波形y_f订ter_out=filter (B, 1, xin)% Y = FILTER(B, A, X) figure (3)plot (y_filter_out)%滤波后信号频谱yff=fft(y_filter_out);figure (4)plot (df,yff)xin=xin/max(xin) ;%归一化xto_ccs 二 round (32767*xin) % 取整fid=fopen(, inputl. dat,， w) ;%打开文件fprintf (fid/ 1651 1 0 0 0n ) ;%输出文件头fprintf (fid, %dn，, xto_ccs) 输出 fclose (fid);