我们有一个不同的截止值的两个低通滤波器:
b, a = signal.butter(2, 0.125)
b2, a2 = signal.butter(2, 0.140)
当应用第一过滤器来x[0:10000]
和第二到x[10000:20000]
与lfilter
,我们必须使用初始条件的输出为“连续的”,如在的答案这里看到连续性问题,当施加在连续时间帧的IIR滤波器 :
zi = lfilter_zi(b, a)
x[0:10000], zi = lfilter(b, a, x[0:10000], zi=zi)
x[10000:20000], zi = lfilter(b2, a2, x[10000:20000], zi=zi)
问题:如何做同样的施加时filtfilt
(向前和向后滤波),以在连续的块使用过滤器时,为了确保连续性,因为没有zi
初始条件参数?