基于DSP的数字滤波器的设计与仿真

MATLABCCS设定模拟信号的采样频率为48000Hz，。设计一种FIR低通滤波器，其参数为： FIR低通滤波器采样频率 阻带截止频率：16000Hz通带最大衰减 阻带最少衰减 滤波器系数 由MATLAB根据前述参数求得计／分析工具FDATOOL设计FIR滤波器图1HzhnznN1

ynhkxnxnkhkN是滤波器的阶数Yn表ynxn、hn

ynxnhnFDATOOL,根据滤波规定设立滤波器类型、通带截止频率、指定DesignFilter，生成滤波器系数。把生成的滤波器系数传到目的DSP。选择菜单Targets->ExporttoCodeComposerStudio(tm)IDEExporttoCHeaderFileCheaderfile32b，16bOK，保存该头文献，需指定文献名(filtercoeff．h)c：\ti\myprojects\fir工程中)。coeff.word,在每行的前面加上.word，比且把每行的最后的逗号去掉。4在本实验中使用MATLABMATLABDSPDSP中实现带通CCSDSP的滤波效果。MATLAB4.1生成fir图4.1MATLABCCS4.2是是是否有待滤波数据图4.2CCS5MATLABDSPFIRDSP中编写解决程序，MATLABFDATOOL能快速设计一种FIR,然后把滤波器系数以头文献形式导入CCS中,头文献中MATLAB辅助DSP实现FIR数字滤波器含滤波器阶数和系数数组,在MATLABDSP程序并显示、分析解决后的数据。使用该办法,便MATLAB中设计的滤波器系数变MATLAB中先编写程序，查看待滤波和已滤波的信号的时域和频域5.15.2所示。图5.1待滤波信号的时域 图5.2待滤波信号的频域FDATOOLFIR5.3图 FIR滤波器的参ExportNumeratorB5.45.4Export5.55.6图5.5已滤波信号时域波形 图5.6已滤波信号频域波形16bitOKfir.hC语言头文献。将系5.7MATLAB在CCS中编写汇编语言程序，进行调试，实现带通滤波的功效。在CCSIDELHM.pjtfir.asm。另外根据fir.cmdMATLAB生成的LHM.hinput1.datLHM.pjt工程文献夹下，进行编译、链接,生成LHM.out。data5.85.8dataMATLAB程序设计中，采样频率设立的过小，截止频率不小于采样FIR.mCCS在.hcoefffir.hinput.dat(7)在View的GraphTime/frequencygraphpropertydialog框，未修改抽样点数。显示的图形出现差别CCSViewGraphTime/frequencygraphpropertydialog框。5.9所示，输出滤波前5.9GraphpropertydialogOK5.105.105.115.115.125.135.125.13由上述的全部截图可知，输入参数中，通带截止频率和阻带截止频率是比较重要的两个参数，不适宜过小，过小滤掉的高频成分太多，很难还原出远波形的形状，也不适宜过大，过大滤波效果欠佳。在本课程4000Hz4500H显示图形时的，时域显示的抽样点数为1024，频域显示的带宽为0.5Hz5.105.125.11和5.13能够懂得，滤波前的时域信号的波形图的波形很杂乱，很难看出信FIR COEF： - - start: PORTRPA1,@x+1 PORTW*AR2+,PA0 PORTRFIR.title .usect .usect table: start: #- PORTRPA1,@xn MAC*AR2+0%,*AR3+0%,A PORTW@y,PA0 PORTRlhm.h-85,-64,-61,-36,8,62,110,136,96,42,-11,-44,-44,-11,39,85,88,37,-29,-83,-101,-73,-7,70,131,82,-5,-96,-151,-144,-72,37,191,162,61,-76,-191,-231,-174,-35,256,278,180,-7,-208,-337,-330,-176,316,445,392,156,-178,-475,-599,-470,-361,745,857,594,-6,-748,-1336,-1456,-386,2192,4154,5816,6768,6768,5816,4154,386,-893,-1456,-1336,-748,-6,594,857,361,-108,-470,-599,-475,-178,156,392,316,72,-176,-330,-337,-208,-7,180,256,132,-35,-174,-231,-191,-76,61,191,140,37,-72,-144,-151,-96,-5,131,124,70,-7,-73,-101,-83,-29,88,106,85,39,-11,-44,-44,-11,96,131,136,110,62,8,-36,-61,--Fir.asm .global_c_int00 .set154 .set ;TheLengthofInput.usect(ORDER-.usect(ORDER- output output AR0, AR2, AR3, AR5, AR6,;Copy #ORDER- #1,FIR_INDEX #xn,FIR_DATA ;Copyxn(data)toRPTZA,#ORDER- A, #(xn+ORDER-1),FIR_DATA #(a0+ORDER-1),FIR_COEFSTM#input,DATA_IN ;Getdatafrom"input"STM#output,DATA_OUT ;Writedatato"output"STM#D_LEN-1,BRCRPTBDnext-STM#ORDER,*D