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

1、（2）旁瓣幅度下降速度要大，以利增加阻带衰减；（3）主瓣的宽度要窄，以获得较陡的过渡带。通常上述三点很难同时满足。当选用主瓣宽度较窄时，虽然得到较陡的过渡带，但通带和阻带的波动明显增加；当选用最小的旁瓣幅度时，虽能得到匀滑的幅度响应和较小的阻带波动，但过渡带加宽。因此，实际选用的窗函数往往是它们的折衷。在保证主瓣宽度达到一定要求的条件下，适当牺牲主瓣宽度来换取旁瓣波动的减少。表2.1几种常用的窗函数对比窗函数 旁瓣峰值衰减 过渡带 阻带最小衰减矩形窗 -13 -21三角形 -27 -25汉宁窗 -31 -44 海明窗 -41 -53布莱克曼窗 -57 -74凯瑟窗 -57 -80比较可得：凯塞

2、窗可提供变化的过渡带宽，通过改变B的值可达到最陡的过渡带：凯塞窗具有与海明窗相匹敌的特性，通过调整p的值，可将凯塞窗完全等价海明窗；凯塞窗最大旁瓣值比主瓣约低80dB，在所有的窗函数中旁瓣抑制度最高2。2.2.2用频率抽样法设计FIR滤波器所谓频率抽样法就是从频域出发，根据频域的采样定理，对给定的理想滤波器的频域响应进行等间隔采样。 （2.2）其中k=0,1,2N-1把当作待设计的滤波器频率响应的采样值,通过下式可以求出滤波器的系统函数和频率响应： （2.3）（2-3） （2.4）其中，是一个内插函数： （2.5）由于频谱的有限个采样值恢复出来的频率响应实际上是对理想频率响应的逼近，因此，这种

3、方法必然有一定的逼近误差。若被逼近的频率响应比较平滑，则各采样点之间的逼近误差较小；反之，则逼近误差较大。2.2 主控模块TMS320C5402 的主要特性:T MS320C5402 是T I 公司于1996 年推出的一种定点DSP 芯片, 采用先进的修正哈佛结构和8 总线结构, 使处理器的性能大大提高。其独立的程序和数据总线, 允许同时访问程序存储器和数据存储器, 实现高速并行操作。例如, 可以在一条指令中, 同时执行3 次读操作和1 次写操作。TMS320C5402 的运行速度为40MIPS （ 每秒执行40 百万条指令） 。指令周期为25ns。此外, 还可以在数据总线与程序总线之间相互传

4、送数据。从而使处理器具有在单个期内同时执行算术运算, 逻辑运算, 位移操作, 乘法累加运算以及访问程序和数据存储器的强大功能。包括: CPU, 总线, 存储器,在片外围电路, 串行口, 以及外部总线接口等部分组成。1 CPU40 位算数逻辑运算单元（ ALU ） , 包括1 个40 位桶形移位寄存器和2 个独立的40 位累加器。􀀂 17 位X 17 位并行乘法器, 与40 位专用加法器相连, 用于非流水线式单周期乘法/ 累加（ MAC）运算。比较, 选择, 存储单元（ CSSU） , 用于加法/比较选择。双地址生成器, 包括8 个辅助寄存器和2 个辅助寄存器算数运算单元（ A

5、RAU ） 。2、总线结构程序总线（ PB） 传送取自程序存储器的指令代码和立即数。3 条数据总线（ CB, DB 和EB） 将内部各单元连接在一起。4 条地址总线（ PAB, CAB, DAB 和EAB ）传送执行指令所需的地址。3 存储器192K 字可寻址存储空间（ 64K 字程序存储器, 64K 字数据存储器, 以及64K 字I/ O 空间） 。片内2K 字ROM。片内10K 字双存取RAM。4 在片外围电路 用来周期产生中断的定时器。8 位主机接口（ HPI） , 用来与主设备或主处理器接口。一个自动缓冲串口（ BSP） 和时分多路串行口（ TPM） , 都可用作标准同步串口。3 系统

6、软件设计3.1 软件整体设计FIR滤波器的实现结构有以下几种形式（1） 直接型结构 该种结构实现时需N个乘法器、N-1个延迟器和N-1个加法器。结构图如下:图3.1 直接型结构图（2） 级联型结构FIR系统函数H（z）是z-1的N-1次多项式，其分解为实系数二阶因子的乘积形式为：H（z）=（0k+1kz-1+2kz-2） 根据上式就可画出级联结构如下图。其中Nc为N/2的最大整数。该结构即是把H（z）的共轭零点或两个单个零点组成基本二阶节，H（z）为基本二阶因子的子系统函数积。图3.2 级联型结构图从图3.2中可看出，每个二阶因子用直接型实现，每一节控制一对零点，便于传输零点的控制。但系数1k

7、的个数比h（n）的个数多，其乘法次数比直接型多。（3） 线性相位型结构线性相位特性是指滤波器对不同频率的正弦波所产生的相移和正弦波的频率成直线关系。在滤波器通带内的信号通过滤波器后，除了由相频特性的斜率决定的延迟外，可以不失真地保留通带以内的全部信号。线性相位的因果的FIR系统的单位取样响应具有如下特性h（n）=h（N-1-n） 当N为偶数时，式变为h（n）z-n+z-（N-1-n） 当N为奇数时，此式为H（z）= h（n）z-n+z-（N-1-n）+h（）z-（） 由上面两个式子可知:实现直接形式的结构网络只需要N/2（N为偶数）或（N十1）/2次乘法（N为奇数），利用h（n）的对称比直接卷

8、积形式少用一半的乘法器。图3.3 N为偶数，线性相位结构图3.4 N为奇数，线性相位结构（4） 频率取样结构系统函数H（z）在单位圆上作N等分取样的值得到h（n）的离散傅立叶变换值H（k）:H（k）=H（）=|H（k）| = h（k） 根据内插公式，得（1-z-N） 上式即为频率取样结构的系统函数其中，H（k）为单位圆上的频率取样值。可见，该系统是由一个子FIR系统（1-z-N）和一个子IIR系统级联组成的。图3.5 频率取样结构图FIR滤波器的实现方法数字滤波器的实现方法一般有以下几种（1）用计算机软件实现 软件实现方法即是在通用的微型计算机上用软件实现。利用计算机的存储器、运算器和控制器把

9、滤波所要完成的运算编成程序通过计算机来执行，软件可由使用者自己编写，也可以使用现成的。但是这种方法速度慢，难以对信号进行实时处理，虽然可以用快速傅立叶变换算法来加快计算速度，但要达到实时处理要付出很高的代价。（2）采用DSP（Digital Signal Processing）处理器来实现 DSP处理器是专为数字信号处理而设计的，它的主要数字运算单元是一个乘累加器，能够在一个机器周期内完成一次乘累加运算，配有适合于信号处理的指令，具备独特的循环寻址和倒序寻址能力。这些特点都非常适合数字信号处理中的滤波器设计的有效实现，并且它速度快，成本低。用DSP芯片实现数字滤波除了具有稳定性好、精确度高、不

10、受环境影响外，还具有灵活性好的特点。（3）采用固定功能的专用信号处理器实现 专用信号处理器采用专用集成电路ASIC来实现，适用于过程固定而又追求高速的信号处理任务，是以指定的算法来确定它的结构，使用各种随机逻辑器件组成的信号处理器。它们体积小、保密性好，具有极高的性能，然而灵活性差。（4）用FPGA等可编程器件实现 使用相关开发工具和VHDL等硬件开发语言，通过软件编程用硬件实现特定的数字滤波算法。这一方法由于具有通用性的特点并可以实现算法的并行运算，无论是作为独立的数字信号处理，还是作为DSP芯片的协作处理器都是比较活跃的研究领域。可以采用DSP或FPGA来实现硬件电路。3.2 测量模块 C

11、CS算法流程图4 系统测试与总结4.1系统测试在View的Graph中单击Time/frequency出现graph property dialog框。将显示类型，图形名称，起始地址，抽样点数，数据类型等分别进行设置，输出各种波形。输出滤波前的信号时域波形。首先，设置如图4.1所示，输出滤波前的信号时域波形。图4.1 Graph property dialog设置框图单击OK后生成如下图4.2波形。图4.2滤波前信号波形图图4.2为滤波前的信号，波形很杂乱，从时域上很难看出信号的周期性。修改相应设置，输出如图4.3所示滤波前信号频谱。图4.3滤波前信号频谱波形图经过滤波后，时域波形频谱波形如图

12、4.4和4.5所示。图4.4滤波后信号时域波形图图4.5滤波后信号频谱图4.2总结体会通过本次课程设计，让我了解了数字滤波是信号处理技术中的重要部分，研究了数字滤波器的基本理论知识以及它实现方法。学习了数字滤波器的结构、设计理论，掌握了各种数字滤波器的原理和特性。并且复习了MATLAB的相关知识，并且利用MATLAB与CCS设计了低通滤波器以及高通滤波器，还研究了如何在定点DSP中实现数字滤波器的算法，掌握了CCS环境下的程序开发方法、调试工具的使用及优化级别的选择等。研究了MATLAB环境下FIR数字滤波器的设计方法，利用MATLAB软件编程实现FIR滤波器设计。DSP作为一门新兴学科，越来

13、越引起人们的关注，目前已广泛应用在各个领域。所以在今后的工作和学习生活中，综合运用模拟电子、数字电子和DSP基本原理等课程中所学的理论知识去独立完成一个项目的设计。不仅要完善以上的MATLAB设计方法，更要学会利用多种途径设计各种DSP最小系统，熟练C语言编程，用DSP芯片实现低通FIR数字滤波，除了具有准确度高、不受环境影响等优点外。通过对本课题的研究，关于数字滤波器理论知识得到了丰富和提高，熟悉了课题研究的整体思路和方法，认识到进行课题研究应具备的素质 参考文献1王念旭DSP基础与应用系统设计北京：北京航空航天大学出版社，20012张雄伟DSP集成开发与应用实例北京：电子工业出版社，200

14、23张雄伟DSP芯片的原理与开发应用（第二版）北京：电子工业出版社，20004彭宗启DSP技术原理及应用北京：电子工业出版社，19985朱铭铭基于模糊的DSP系统设计北京：电子工业出版社，20036李真芳，苏涛等.DSP程序开发一MATLAB调试及直接目标代码生成M.西安:西安电子科技大学出版社，2003附录用线性缓冲区实现FIR滤波器程序清单： .title “FIR1.ASM” .mmregs .def start x .usect “x”，8 PA0 .set 0 PA1 .set 1 .data COEF： .word 1*32768/10 .word 2*32768/10 .word

15、 -4*32768/10 .word 3*32768/10 .textstart: SSBX FRCT STM #x+7,AR2 STM #6,AR0 LD #x+1,DP PORTR PA1,x+1FIR1: RPTZ A,#6 MACD *AR2-,COEF,A STH A,*AR2 PORTW *AR2+,PA0 BD FIR1 PORTR PA1,*AR2+0 .end 用循环缓冲区实现FIR滤波器程序清单： .title “FIR2.ASM” .mmregs .bss y,1xn .usect “xn”,7b0 .usect “b0”,7PA0 .set 0PA1 .set 1tab

16、le: .word 4*32768/10 .word 5*32768/10 .word 6*32768/10 .word 7*32768/10 .text STM #b0,AR1 RPT #6 MVPD table,*AR1+ STM #xn+6,AR2 STM #b0+6,AR3 STM #7,BK STM #-1,AR0 LD #xn,DP PORTR PA1,xnFIR2: MAC *AR2+0%,*AR3+0%,A STH A,y PORTW y,PA0 BD FIR2 PORTR PA1,*AR2+0%汇编程序清单lhm.h文件内容：coeff .word -85, -64, -61

17、, -36, 8, 62, 110, 136, 131 .word 96, 42, -11, -44, -44, -11, 39, 85, 106 .word 88, 37, -29, -83, -101, -73, -7, 70, 124 .word 131, 82, -5, -96, -151, -144, -72, 37, 140 .word 191, 162,61, -76, -191, -231, -174, -35, 132 .word 256, 278, 180, -7, -208, -337, -330, -176, 72 .word 316, 445, 392, 156, -

18、178, -475, -599, -470, -108 .word 361, 745, 857, 594, -6, -748, -1336, -1456, -893 .word 386, 2192, 4154, 5816, 6768, 6768, 5816, 4154, 2192 .word 386, -893, -1456, -1336, -748, -6, 594, 857, 745 .word 361, -108, -470, -599, -475, -178, 156, 392, 445 .word 316, 72, -176, -330, -337, -208, -7, 180, 2

19、78 .word 256, 132, -35, -174, -231, -191, -76, 61, 162 .word 191, 140, 37, -72, -144, -151, -96, -5, 82 .word 131, 124, 70, -7, -73, -101, -83, -29, 37 .word 88, 106, 85, 39, -11, -44, -44, -11, 42 .word 96, 131, 136, 110, 62, 8, -36, -61, -64 .word -85Fir.asm内容： .title fir.asm .mmregs .global _c_in

20、t00ORDER .set 154 D_LEN .set 1024 ; The Length of Input Dataxn .usect xn, （ORDER-1）a0 .usect a0input .usect input, D_LEN output .usect output, D_LEN .copy HM.h ; Get coefs from coef.h .text .asg AR0, FIR_INDEX .asg AR2, FIR_DATA .asg AR3, FIR_COEF .asg AR5, DATA_IN .asg AR6, DATA_OUT_c_int00: SSBX F

21、RCT STM #a0,FIR_COEF ; Copy a0（coefs） to FIR_COEF（AR3） RPT #ORDER-1 MVPD #coeff,*FIR_COEF+ STM #1, FIR_INDEX STM #xn, FIR_DATA ; Copy xn（data） to FIR_DATA（AR2） RPTZ A,#ORDER-1 STL A, *FIR_DATA+ STM #（xn+ORDER-1）, FIR_DATA STM #（a0+ORDER-1）, FIR_COEF STM #input, DATA_IN ; Get data from STM #output, D

22、ATA_OUT ; Write data to STM #D_LEN-1, BRC RPTBD next-1 STM #ORDER, BK LD *DATA_IN+, AFIR: STL A, *FIR_DATA+% RPTZ A, （ORDER-1） MAC *FIR_DATA+0%, *FIR_COEF+0%, A ; FIRS STH A, *DATA_OUT+nextFIR_END: B FIR_END .endFIR.m的程序如下：fs=40000; %/采样HzN=1024 %数据个数T=1/fs; %采样周期n=0:N-1;df=n*（fs/N）%待滤波信号波形xin=randn

23、（1,1024）figure（1）plot（xin）%待滤波信号频谱xinff=abs（fft（xin）;figure（2）plot（df,xinff）%滤波后信号波形y_filter_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（input1.dat,w）;%打开文件fprintf（fid,1651 1 0 0 0n%输出文件头%dn,xto_ccs）;%输出fclose（fid）;