用python建模_用Python建模一个线性系统

我知道这有点旧，但搜索引起了我的这个问题。当我找不到一个好的模块时，我把它放在一起。这并不多，但如果有人发现他们在这里，这是一个好的开始。

import matplotlib.pylab as plt

import numpy as np

import scipy.signal

def bode(G,f=np.arange(.01,100,.01)):

plt.figure()

jw = 2*np.pi*f*1j

y = np.polyval(G.num, jw)/np.polyval(G.den, jw)

mag = 20.0*np.log10(abs(y))

phase = np.arctan2(y.imag, y.real)*180.0/np.pi % 360

plt.subplot(211)

#plt.semilogx(jw.imag, mag)

plt.semilogx(f,mag)

plt.grid()

plt.gca().xaxis.grid(True, which='minor')

plt.ylabel(r'Magnitude (db)')

plt.subplot(212)

#plt.semilogx(jw.imag, phase)

plt.semilogx(f,phase)

plt.grid()

plt.gca().xaxis.grid(True, which='minor')

plt.ylabel(r'Phase (deg)')

plt.yticks(np.arange(0, phase.min()-30, -30))

return mag, phase

f=scipy.signal.lti([1],[1,1])

bode(f)

编辑：我回到这里，是因为有人upvoted这个答案，你应该尝试Control Systems Library。他们已经使用匹配的语法和一切来实现大部分Matlab控制系统工具箱。