fir滤波器设计.doc

FIR滤波器设计报告学号：1020300148姓名：杨昆FIR滤波器的设计（1）滤波器的设计要求：利用MATLAB和CCS仿真环境设计一个FIR滤波器，设计滤波器的采样平率为600HZ，信号频率为100，220，280HZ的双正弦不等幅信号，要求滤去200HZ以上的信号成分。（2）设计思路：1. 低通滤波器，不需要很高的阶数。2. 在Matlab语言中确定滤波器的各种参数。3. 将Matlab语言转换成通用语言形式。4. 编写实际DSP应用程序。5. 在CCS中调试程序。（3）设计过程：1.使用Matlab中的bartlett窗函数。下面是在Matlab中设计的滤波器程序。clear all;clf;N=256;fs=600;dt=1/fs;for t=1:N f1=280; f2=100; f3=220; y(t)=(sin(2*pi*f1*t*dt)+sin(2*pi*f2*t*dt)+1.2*sin(2*pi*f3*t*dt); x(t)=round(10000*y(t);end figure(1)plot(y),title();fp=200fst=250NL=56 f1=200;w1=2*pi*f1/fswindow=bartlett(NL+1) n=1:NL+1hd=sin(w1*(n-NL/2)./(pi*(n-NL/2);hd(NL/2)=0.67;h=hd.*rot90(window)h1=h*32768;h2=fft(h,N);pyy=h2.*conj(h2);f=(0:(N/2-1);for i=1:N/2-1 f(i)=f(i)*fs/N;end figure(2)plot(f,pyy(1:N/2),title();yy1=filter(h,1,y);figure(3)plot(yy1),title();y=fft(y,N);pyy=y.*conj(y);f=(0:(N/2-1);for i=1:N/2-1 f(i)=f(i)*fs/N;end figure(4)subplot(1,2,1),plot(f,pyy(1:N/2),title();y=fft(yy1,N);pyy=y.*conj(y);f=(0:(N/2-1);for i=1:N/2-1 f(i)=f(i)*fs/N;end subplot(1,2,2),plot(f,pyy(1:N/2),title()2.仿真结果 图一图二 图三 图四 图一为滤波器的频谱图像，图二为滤波前后波形频谱对比，图三位滤波前信号波形，图四为滤波后信号波形。（4）CCS仿真：1.程序如下：.mmregs .ref filter_start .def _c_int00K_DATA_SIZE .set 256K_BUFFER_SIZE .set 64K_STACK_SIZE .set 256K_B .set 27K_CIR .set K_BUFFER_SIZE STACK .usect stack,K_STACK_SIZESYSTEM_STACK .set K_STACK_SIZE+STACKDATA_DP .usect filter_vars,0filterdata .usect filter_vars,K_DATA_SIZEbufferdatay .usect filter_vars,K_BUFFER_SIZE*2bufferdatax .usect filter_vars,K_BUFFER_SIZE*2 .data .global inputdatainputdata .WORD 19657,-7341,12931,-13597,-10392,1902,13725,-6166 .WORD 20923,-27713,1266,1176,793,15499,0,-15499 .WORD -793,-1176,-1266,27713,-20923,6166,-13725,-1902 .WORD 10392,13597,-12931,7341,-19657,0,19657,-7341 .WORD 12931,-13597,-10392,1902,13725,-6166,20923,-27713 .WORD 1266,1176,793,15499,0,-15499,-793,-1176 .WORD -1266,27713,-20923,6166,-13725,-1902,10392,13597 .WORD -12931,7341,-19657,0,19657,-7341,12931,-13597 .WORD -10392,1902,13725,-6166,20923,-27713,1266,1176 .WORD 793,15499,0,-15499,-793,-1176,-1266,27713 .WORD -20923,6166,-13725,-1902,10392,13597,-12931,7341 .WORD -19657,0,19657,-7341,12931,-13597,-10392,1902 .WORD 13725,-6166,20923,-27713,1266,1176,793,15499 .WORD 0,-15499,-793,-1176,-1266,27713,-20923,6166 .WORD -13725,-1902,10392,13597,-12931,7341,-19657,0 .WORD 19657,-7341,12931,-13597,-10392,1902,13725,-6166 .WORD 20923,-27713,1266,1176,793,15499,0,-15499 .WORD -793,-1176,-1266,27713,-20923,6166,-13725,-1902 .WORD 10392,13597,-12931,7341,-19657,0,19657,-7341 .WORD 12931,-13597,-10392,1902,13725,-6166,20923,-27713 .WORD 1266,1176,793,15499,0-15499,-793,-1176 .WORD -1266,27713, -20923, 6166, -13725, -1902, 10392,13597 .WORD -12931,7341,-19657,0,19657,-7341,12931,-13597 .WORD -10392,1902,13725,-6166,20923,-27713,1266,1176 .WORD 793,15499,0,-15499,-793,-1176,-1266,27713 .WORD -20923,6166,-13725,-1902,10392,13597,-12931,7341 .WORD -19657,0,19657,-7341,12931,-13597,-10392,1902 .WORD 13725,-6166,20923,-27713,1266,1176,793,15499 .WORD 0,-15499,-793,-1176,-1266,27713,-20923,6166 .WORD -13725,-1902,10392,13597,-12931,7341,-19657,0 .WORD 19657,-7341,12931,-13597,-10392,1902,13725,-6166 .WORD 20923,-27713,1266,1176,793,15499,0,-15499 .text_c_int00: .asg AR2,ORIGIN .asg AR3,INPUT .asg AR4,FILTER .asg AR5,OUTPUTSTART: SSBX FRCT SSBX INTM LD #DATA_DP,DP STM #STACK,SP CALL filter_start NOP NOP NOPLOOP B LOOP .def h0,h1,h2,h3,h4,h5,h6,h7,h8,h9,h10,h11,h12,h13,h14,h15,h16,h17,h18,h19,h20,h21,h22,h23,h24,h25,h26 .def filter_starth0.set 0000Hh1.set -003FHh2.set 008AHh3.set 0000Hh4.set -015BHh5.set 01F0Hh6.set 0000Hh7.set -03CDHh8.set 056EHh9.set 0000Hh10.set -0D92Hh11.set 1DDBHh12.set 4F2AHh13.set 2349Hh14.set -1049Hh15.set 0000Hh16.set 06C9Hh17.set -04E3Hh18.set 0000Hh19.set 02B7Hh20.set -0209Hh21.set 0000Hh22.set 0116Hh23.set -00BEHh24.set 0000Hh25.set 0035Hh26.set 0000H .textfilter_start: STM #K_CIR,BK STM #1,AR0 STM #inputdata,ORIGIN STM #bufferdatax,INPUT ; STM #bufferdatay,FILTER STM #filterdata,OUTPUT ;滤波 STM #K_DATA_SIZE-1,BRC RPTB filter_end-1 MVDD *ORIGIN+,*INPUT RPT #K_B-1-1 MAR *INPUT-0% MPY*INPUT+0%,#h26,BLDB,AMPY*INPUT+0%,#h25,BADDB,AMPY*INPUT+0%,#h24,BADDB,AMPY*INPUT+0%,#h23,BADDB,AMPY*INPUT+0%,#h22,BADDB,AMPY*INPUT+0%,#h21,BADDB,AMPY*INPUT+0%,#h20,BADDB,AMPY*INPUT+0%,#h19,BADDB,AMPY*INPUT+0%,#h18,BADDB,AMPY*INPUT+0%,#h17,BADDB,AMPY*INPUT+0%,#h16,BADDB,AMPY*INPUT+0%,#h15,BADDB,AMPY*INPUT+0%,#h14,BADDB,AMPY*INPUT+0%,#h13,BADDB,AMPY*INPUT+0%,#h12,BADDB,AMPY*INPUT+0%,#h11,BADDB,AMPY*INPUT+0%,#h10,BADDB,AMPY*INPUT+0%,#h9,BADDB,AMPY*INPUT+0%,#h8,BADDB,AMPY*INPUT+0%,#h7,BADDB,AMPY*INPUT+0%,#h6,BADDB,AMPY*INPUT+0%,#h5,BADDB,AMPY*