1 from scipy import signal 2 import matplotlib.pyplot as plt 3 4 s1 = signal.lti([1, 2], [1, 4, 4, 6, 0]) 5 w, mag, phase = signal.bode(s1) 6 7 plt.figure() 8 plt.semilogx(w, mag) 9 plt.figure() 10 plt.semilogx(w, phase) 11 plt.show() 12