plt.semilogx(f, mag)
plt.grid()
plt.gca().xaxis.grid(True, which='minor')
plt.ylabel('Magnitude (db)')
plt.subplot(212)
plt.semilogx(f, phase)
plt.grid()
plt.gca().xaxis.grid(True, which='minor')
plt.ylabel('Phase (deg)')
plt.yticks(np.arange(0, phase.min() - 30, -30))
return mag, phase
f = scipy.signal.lti([1], [1, 1])
bode(f)
Моделирование линейной системы с помощью Python