zfiltfilt.m

来自「用MATLAB实现数字信号处理的程序」· M 代码 · 共 12 行

%使用Filtfilt函数进行滤波
t=0:0.001:0.1;
x1=sin(2*pi*40*t);
x2=0.5*rand(size(t));           %添加随机噪声
x=x1+x2;
A=[1 -1.143 0.4128];            %滤波器分母系数
B=[0.06745 0.1348 0.06745];     %滤波器分子系数
y=filter(B,A,x);                %普通滤波器
z=filtfilt(B,A,x);              %零相位滤波
plot(t,x,'g:',t,y,'r--',t,z,'b*');
legend('原始信号','filter(X)','filtfilt(X)');