filter函数 与filtfilt函数的效果区别
filter滤波器称为一维数字滤波器。filtfilt滤波器称为零相位数字滤波。其滤波算法是基于filter而来的。只是filtfilt实现了零相位。其基本实现过程为先让信号用filter滤波,再将信号时域反转再次通过filter滤波,这样两次滤波后相位实现了零相位。filter滤波有明显的延迟,filtfilt滤波延时比filter小很多。分析如下(程序参见matlab脚本,完整程序下载链接):
(1)设置三种频率(10Hz,20Hz,30Hz)构成的信号y1
f1=10;
f2=20;
f3=30;
Fs=150;
t=0:1/Fs:1;
y1=5*sin(2*pi*f1*t)+10*sin(2*pi*f2*t)+8*sin(2*pi*f3*t);
(2)设计滤波参数,分别用filter滤波器和filtfilt滤波器把30Hz的信号滤掉。
%滤波
fp=25;%滤掉频率超过25Hz的信号
fs=40;
wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; %DF指标(低通滤波器的通、阻带边界频)
[N,wp]=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆DF阶数N和通带截止频率wp
[B,A]=ellip(N,rp,rs,wp); %调用ellip计算椭圆带通DF系统函数系数向量B和A
y2=filter(B,A,y1); %滤波
y3=filtfilt(B,A,y1);
(3)绘制滤波效果图和FFT频谱分析图,对比滤波效果。
如图所示,第1个时域图可以发现经过filter滤波的信号比原始信号明显有延迟,而经过filtfilt滤波后的信号几乎没有延迟。第2频谱图可以发现,按照设计的参数,两个滤波器都能把频率为30Hz的信号滤除了。