def sosfilter(sos, zi_in, x): y = x zi_out = zi_in for i in range(len(sos)): y, zi_out[i] = lfilter(sos[i,:3], sos[i,3:], y, zi = zi_in[i]) return y, zi_out