Replace bilinear_zpk with the scipy implementation
Created by: asouthgate
Added an extra test to make sure complex filter dtypes are preserved
Addresses https://github.com/gwpy/gwpy/issues/1630
Using scipy's implementation we also avoid problems with old implementation casting to float as in https://github.com/gwpy/gwpy/issues/1630#issuecomment-1721674653