Return a digital filter from an analog one using a bilinear transform.
The bilinear transform substitutes (z-1) / (z+1) for s.
scipy.signal.upfirdn
scipy.signal.findfreqs