-->

适用filtfilt与初始条件(以避免不连续性)连续块(Apply filtfilt on succ

2019-11-05 04:49发布

我们有一个不同的截止值的两个低通滤波器:

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初始条件参数?

文章来源: Apply filtfilt on successive blocks with initial conditions (to avoid discontinuity)