我们有两个截止值不同的低通滤波器:
b, a = signal.butter(2, 0.125)
b2, a2 = signal.butter(2, 0.140)
在将第一个过滤器应用于x[0:10000]
并将第二个过滤器应用于x[10000:20000]
时,我们必须使用初始条件使输出“连续”,如Continuity issue when applying an IIR filter on successive time-frames的答案:
lfilter
问题:由于没有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
(正向和反向过滤)时如何做以确保连续块上使用过滤器时的连续性?