from scipy import signal
import matplotlib.pyplot as plt

s1 = signal.lti([1, 2], [1, 4, 4, 6, 0])
w, mag, phase = signal.bode(s1)

plt.figure()
plt.semilogx(w, mag)
plt.figure()
plt.semilogx(w, phase)
plt.show()