超低截止频率高通滤波器基线快速恢复

直接滤波

t=0:1/500:50;
x=sin(2*pi*t)+10;
[b,a]=butter(2,0.05/500,'high');
y1=filter(b,a,x);
figure
subplot(211)
plot(x)
title('ori')
subplot(212)
plot(y1)
title('filted')

计算初始状态

z=filtic(b,a,[0,0],[1,1]);
y1=filter(b,a,x,z*x(1));
figure
subplot(211)
plot(x)
title('ori')
subplot(212)
plot(y1)
title('filted')

fiftic函数是用于计算初始化状态，对于高通滤波器而言，初始状态应该是输入是1，输出是0