滤波器分析和总结.docx

一、fir 数字滤波器 1、MATLAB 、fir 低通滤波器 b=fir1(79,0.1,hamming(80)); %80 point hamming window FIR filter, Wn=0.2 freqz(b,1) %fvtool(b,1) %Flter Visualization Tool 可心以查看幅度响应、冲激响应等 h=buffer(round(b*32768),8) %参数*2^15,把小数点移最高位后面后取整，然后排成 8 个一行便于输出。 csvwrite(table.txt,h) %参数表输出到文件中 %然后只要把参数表从 table.txt 中拷到代码中，每行前面加上 .word -3,-10,-17,-24,-30,-35,-37,-34 -26,-10,12,41,74,107,136,156 161,145,105,41,-46,-150,-263,-371 -461,-518,-524,-468,-339,-132,151,502 908,1347,1797,2228,2615,2931,3155,3271 3271,3155,2931,2615,2228,1797,1347,908 502,151,-132,-339,-468,-524,-518,-461 -371,-263,-150,-46,41,105,145,161 156,136,107,74,41,12,-10,-26 -34,-37,-35,-30,-24,-17,-10,-3 、fir 高通滤波器 clc wp=0.3*pi; ws=0.2*pi; Bt=wp-ws NO=ceil(6.6*pi/Bt); N=NO+mod(NO+1,2) wc=(wp+ws)/2/pi; b=fir1(N-1,wc,high,Hamming(N)); %freqz(b,1,512) freqz(b,1,512,15000) h=buffer(round(b*32768),8); csvwrite(tabl.txt,h) 、fir 带通滤波器wp1=0.2;wp2=0.3; ws1=0.1;ws2=0.4; wp=[wp1,wp2];ws=[ws1,ws2]; B=wp1-ws1;%计算过渡带宽N0=ceil(6.6*pi/B); N=N0+mod(N0+1,2) wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2; wc=[wc1,wc2]; b=fir1(N-1,wc,bandpass,hamming(N)); %80 point hamming window FIR filter, Wn=0.2 fvtool(b,1) %Flter Visualization Tool可心以查看幅度响应、冲激响应等h=buffer(round(b*32768),8) %参数*2^15,把小数点移最高位后面后取整，然后排成 8 个一行便于输出。 csvwrite(table.txt,h) %参数表输出到文件中 %然后只要把参数表从 table.txt 中拷到代码中，每 行前面加上 .word 15,9,0,-4,0,4,0,-11 -19,-15,0,13,15,6,0,6 18,19,0,-27,-38,-24,0,11 0,-12,0,35,61,48,0,-44 -47,-19,0,-21,-59,-60,0,83 117,74,0,-31,0,34,0,-99 -172,-134,0,118,127,50,0,54 152,155,0,-209,-293,-184,0,77 0,-84,0,239,415,323,0,-286 -307,-119,0,-131,-370,-378,0,516 731,464,0,-198,0,222,0,-660 -1175,-938,0,886,991,404,0,498 1512,1687,0,-2935,-4943,-3970,0,4557 6554,4557,0,-3970,-4943,-2935,0,1687 1512,498,0,404,991,886,0,-938 -1175,-660,0,222,0,-198,0,464 731,516,0,-378,-370,-131,0,-119 -307,-286,0,323,415,239,0,-84 0,77,0,-184,-293,-209,0,155 152,54,0,50,127,118,0,-134 -172,-99,0,34,0,-31,0,74 117,83,0,-60,-59,-21,0,-19 -47,-44,0,48,61,35,0,-12 0,11,0,-24,-38,-27,0,19 18,6,0,6,15,1