%使用Filtfilt函數(shù)進(jìn)行濾波
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];            %濾波器分母系數(shù)
B=[0.06745 0.1348 0.06745];     %濾波器分子系數(shù)
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)');