Scipy.signal.iirfilter equivalent in DSP.jl and more

One problem of the current implementation can be shown here (input signal is rolling)

image

orange signal (filtered signal) starts from 0, it should start from it’s previous value (so nearer from -1).

filter should “keep” an internal state