butterworth滤波例子

liushuangq05

查阅了下本版块,进一步增加一个butterworth带通滤波的例子和程序给大家,希望对新手有用.
采用3阶带通滤波,数据采样率为200点每秒,所以下面的程序段设置为[0.4/100 20/100].
[b,a]=butter(3,[0.4/100 20/100]);  %0.4,20是hz单位


滤波器幅频特性如下图:

liushuangq05

clear;clc; nam=dir('data\*.txt'); aa='data\'; [k1,k2]=size(nam); [b,a]=butter(3,[0.4/100 20/100]); %0.4,20是hz单位 if k1~=0 for k=1:1:length(nam(:,1)) %for k=1:1:1 namname{k,1}=nam(k,1).name; %end [k1,k2]=size(namname); gg=strcat(aa,namname{k,1}); c=strtok(gg,'.'); c=strcat(c,'.flt'); fp=fopen(gg,'r'); fp1=fopen(c,'w'); data=fscanf(fp,'%f %f %f',[3,inf]); d1=data(1,:); d2=data(2,:); d3=data(3,:); %clear data; dd1=filter(b,a,d1); dd2=filter(b,a,d2); dd3=filter(b,a,d3); data=[dd1; dd2; dd3]; fprintf(fp1,'%f %f %f\n',data); fclose(fp1); fclose(fp); end end plot(d1);hold on;plot(dd1+5,'r');axis tight;hold off 以上为完整程序段. 幅频特性图如下(上图发错了,它是滤波前后对比图):

liushuangq05

测试数据见附件.放在程序目录下的data(自己建立)文件夹中即可. 程序是批处理形式的.

liushuangq05

另外补充IIR型滤波器不是最小相位滤波器,所以存在一定的相移.本算例有3个采样点的延迟.