(完整版)基于DSP的数字滤波器的设计与仿真

1、2.1 系统功能介绍一个实际的应用系统中，总存在各种干扰。数字滤波器在语音信号处理、信号频谱估计、信号去噪、无线通信中的数字变频以及图像信号等各种信号处理中都有广泛的应用，数字滤波器也是使用最为广泛的信号处理算法之一。在本设计中，使用MATLAB模拟产生合成信号，然后利用CCS进行滤波。设定模拟信号的采样频率为48000Hz,。设计一个FIR低通滤波器，其参数为：滤波器名称：FIR低通滤波器采样频率：Fs=48000Hz通带截止频率：15000Hz阻带截止频率：16000Hz通带最大衰减：0.1dB阻带最少衰减：80dB滤波器系数：由MATLAB根据前述参数求得。2.2 总体设计方案流程图图

2、1 总体设计方案13.1 滤波器原理对于一个FIR滤波器系统，它的冲击响应总是又限长的，其系统函数可记为:H(z)=h(n)z-n其中N-1是FIR的滤波器的阶数,z-n为延时结，h(n)为端口信号函数。最基本的FIR滤波器可用下式表示：y(n)=h(k)x(n-k)其中 x(n-k)输入采样序列，h(k)是滤波器系数，N是滤波器的阶数 Y(n)表示滤波器的输出序列，也可以用卷积来表示输出序列 y(n)与x(n)、h(n)的关系，如y(n)=x(n)*h(n)3.2 操作步骤(1) 打开FDATOOL，根据滤波要求设置滤波器类型、通带截止频率、指定阶数、采样频率等。指定完设计参数后单击按钮De

3、signFilter，生成滤波器系数。(2) 把 生 成 的 滤 波 器 系 数 传 到 目 标DSP。 选 择 菜 单Targets-ExporttoCodeComposerStudio(tm)IDE,打开ExporttoCHeaderFile对话框，选择Cheaderfile，指定变量名(滤波器阶数和系数向量)，输出数据类型可选浮点型或32b，16b整型等，根据自己安装选择目标板板号和处理器号，单击OK,保存该头文件，需指定文件名(filtercoeff.h)和路径(保存在c：timyprojectsfir工程中)。(3) 修改CCS汇编程序，删掉数据前的所有文字，在开头加上.data，第

4、二行加coeff.word,在每行的前面加上.word，比且把每行的最后的逗号去掉。(4) 编译汇编程序，如果有错误，按错误进行修改；没错误，则往下执行。(5) 加载初始化DATA数据。运行程序，查看输入输出波形，修改相应参数进行调试主要内容和步骤下：2第 4 章详细设计在本实验中使用MATLAB模拟产生信号， 观察滤波前的时域波形和频域波形。MATLAB仿真后，使用得到的滤波器参数，进行DSP编程，在DSP中实现带通滤波，并使用CCS的频谱分析功能，查看DSP的滤波效果。MATLAB程序流程图如图4.1所示图 4.1MATLAB 程序流程图CCS汇编程序流程图如图4.2所示图 4.2CCS

5、汇编程序流程图第 5 章实验过程35.1 汇编语言实验步骤与内容MATLAB辅助DSP实现FIR，其总体过程为在DSP中编写处理程序，在MATLAB中利用滤波器设计、分析工具（FDATOOL），根据指定的滤波器性能快速设计一个FIR，然后把滤波器系数以头文件形式导入CCS中，头文件中MATLAB辅助DSP实现FIR数字滤波器含滤波器阶数和系数数组， 在MATLAB中调试、运行DSP程序并显示、分析处理后的数据。使用该方法，便于采用汇编语言来实现程序。头文件名不变，当MATLAB中设计的滤波器系数改变时，相应头文件中系数也改变，方便了程序调试、仿真。（1） 在MATLAB中先编写程序， 查看待滤

6、波和已滤波的信号的时域和频域波形，如图5.1和图5.2所示。图 5.1 待滤波信号的时域图图 5.2 待滤波信号的频域图利用FDATOOL设计FIR滤波器的参数，操作步骤如图5.3所示。4打开Export，把Numerator改为B,如图5.4所示。1X1LJ图 5.4Export 设置已滤波的信号的时域和频域波形，如图5.5和图5.6所示图 5.3FIR 滤波器的参数5图 5.5 已滤波信号时域波形图图 5.6 已滤波信号频域波形图MATLAB滤波器设计工具在完成FIR设计后，如图5.7所示，在Exportas中选择16bit符号整数输出，然后单击OK按钮。生成fir.hC语言头文件。将系数

7、稍作修改，调整后拷贝到程序的系数初始化空间即可。(2)在CCS中编写汇编语言程序， 进行调试， 实现带通滤波的功能。 在CCSIDE中建立LHM.pjt工程,用汇编语言编写处理主程序fir.asm。另外根据板上的存储器配置方式,编写存储器配置文件fir.cmd文件，将MATLAB生成的LHM.h和inputl.dat文件拷贝到LHM.pjt工程文件夹下，进行编译、链接,生成可执行文件LHM.out。运行程序，查看输入输出波形，修改相应参数进行调试5.2 实验过程中出现的错误及解决的办法(1)在MATLAB程序设计中，采样频率设置的过小，截止频率大于采样频率的一半，运行的时候图形出现错误。(2)

8、FIR.m中的采样频率要和开始设置时的采样频率一致，否则结果会出现图 5.7MATLAB 输出数值转换图(3)加载初始化data数据，图5.8所示。图 5.8data 初始化图6(3)在CCS中未定义标号，程序运行错误。(4)在.h文件中未把coeff顶格写，编译时出现错误。(5)未将fir.h和input.dat文件考入工程文件中。6)间接寻址过程中，程序中丢失。(7)在View的Graph中单击Time/frequency出现graphpropertydialog框，未修改抽样点数。显示的图形出现差异5.3 CCS 程序运行后的各种输出结果在View的Graph中单击Time/freque

9、ncy出现graphpropertydialog框。将显示类型，图形名称，起始地址，抽样点数，数据类型等分别进行设置，输出各种波形。输出滤波前的信号时域波形。首先，设置如图5.9所示，输出滤波前图5.10为滤波前的信号，波形很杂乱，从时域上很难看出信号的周期性。修改相应设置，输出如图5.11所示滤波前信号频谱。偏差。的信号时域波形。7图 5.11 滤波前信号频谱波形图经过滤波后，时域波形频谱波形如图5.12和 5.13所示图 5.12 滤波后信号时域波形图图 5.13 滤波后信号频谱图由上述的所有截图可知，输入参数中，通带截止频率和阻带截止频率是比较重要的两个参数，不宜过小，过小滤掉的高频成分

10、太多，很难还原出远波形的形状，也不宜过大，过大滤波效果欠佳。在本课程设计中，通带截止频率选用4000Hz，阻带截止频率选用4500Hz，因此，显示图形时的，时域显示的抽样点数为1024，频域显示的带宽为0.5Hz。对比滤波前后信号的时域图5.10和5.12、频谱波形图5.11和5.13可以知道，滤波前的时域信号的波形图的波形很杂乱，很难看出信号的周期性，滤波后周期性也相对的明显起来了；滤波前信号分布在整个频带上，滤波后阻带范围内的基本被滤掉，这些都可以从上述截图看出。8附件：源程序清单用线性缓冲区实现 FIR 滤波器程序清单.title“FIR1.ASM”.mmregs.defstartx.u

11、sect“x”，8PA0.set0PA1.set1.dataCOEF：.word1*32768/10.word2*32768/10.word-4*32768/10.word3*32768/10.word-4*32768/10.word2*32768/10.word1*32768/10.textstart:SSBXFRCTSTM#x+7,AR2STM#6,AR0LD#x+1,DPPORTRPA1,x+1FIR1:RPTZA,#6MACD*AR2-,COEF,ASTHA,*AR2PORTW*AR2+,PA0BDFIR1PORTRPA1,*AR2+0.end用循环缓冲区实现 FIR 滤波器程序清单:

12、.title“FIR2.ASM”9.mmregs.def.bssstarty,1xn.usect“xn”,7b0.usect“b0”,7PA0.set0PA1.set1.datatable:.word1*32768/10.word2*32768/10.word3*32768/10.word4*32768/10.word5*32768/10.word6*32768/10.word7*32768/10.textstart:SSBXFRCTSTM#b0,ARlRPT#6MVPDtable,*ARl+STM#xn+6,AR2STM#b0+6,AR3STM#7,BKSTM#-1,AR0LD#xn,DPP

13、ORTRPAl,xnFIR2:RPTZA,#6MAC*AR2+0%,*AR3+0%,ASTHA,yPORTWy,PA0BDFIR210PORTRPAl,*AR2+0%.end汇编程序清单lhm.h 文件内容.datacoeff.word-85,-64,-61,-36,8,62,110,136,131.word96,42,-11,-44,-44,-11,39,85,106.word88,37,-29,-83,-101,-73,-7,70,124.word131,82,-5,-96,-151,-144,-72,37,140.word191,162,61,-76,-191,-231,-174,-35

14、,132.word256,278,180,-7,-208,-337,-330,-176,72.word316,445,392,156,-178,-475,-599,-470,-108.word361,745,857,594,-6,-748,-1336,-1456,-893.word386,2192,4154,5816,6768,6768,5816,4154,2192.word386,-893,-1456,-1336,-748,-6,594,857,745.word361,-108,-470,-599,-475,-178,156,392,445.word316,72,-176,-330,-337

15、,-208,-7,180,278.word256,132,-35,-174,-231,-191,-76,61,162.word191,140,37,-72,-144,-151,-96,-5,82.word131,124,70,-7,-73,-101,-83,-29,37.word88,106,85,39,-11,-44,-44,-11,42.word96,131,136,110,62,8,-36,-61,-64.word-85Fir.asm 内容：.titlefir.asm.mmregs.global_c_int00ORDER.set154D_LEN.set1024xn.usectxn,(OR

16、DER-1).usecta0,(ORDER-1);TheLengthofInputDataa011STM#a0,FIR_COEF;Copya0(coefs)toFIR_COEF(AR3)RPT#ORDER-1MVPD#coeff,*FIR_COEF+STM#1,FIR_INDEXSTM#xn,FIR_DATA;Copyxn(data)toFIR_DATA(AR2)RPTZA,#ORDER-1STLA,*FIR_DATA+STM#(xn+ORDER-1),FIR_DATASTM#(a0+ORDER-1),FIR_COEFSTM#input,DATA_IN;GetdatafrominputSTM#

17、output,DATA_OUT;WritedatatooutputSTM#D_LEN-1,BRCRPTBDnext-1STM#ORDER,BKLD*DATA_IN+,AFIR:STLA,*FIR_DATA+%RPTZA,(ORDER-1)MAC*FIR_DATA+0%,*FIR_COEF+0%,A;FIRSSTHA,*DATA_OUT+nextinput.usectinput,D_LENoutput.usectoutput,D_LEN.copyHM.h.text.asgAR0,FIR_INDEX.asgAR2,FIR_DATA.asgAR3,FIR_COEF.asgAR5,DATA_IN.asgAR6,DATA_OUT_c_int00:SSBXFRCT;Getcoefsfromcoef.h12FIR_END:BFIR_END.endFIR.m 的程序如下：fs=40000;%/采样 HzN=1024%数据个数T=1/fs;%采样周期n=0:N-1;df=n*(fs/N)%待滤波信号波形xin=randn(1,1024)figure(1)plot(xin)%待滤波信号频谱xinff=abs(fft(xin);figure(2)plot(df,xi