Pythonで二次遅れ系をシミュレーションをする。
\begin{align}
G(s)=\frac{\omega_{n}^2 }{s^2 + 2 \zeta \omega_{n} s + \omega_{n}^2}
\end{align}
コードは以下
from control.matlab import *
import matplotlib.pyplot as plt
s = tf('s')
zeta = 0.5
omega = 1
G = omega ** 2 / (s ** 2 + 2 * zeta * omega * s + omega ** 2)
y, t = step(G, 10)
plt.plot(t, y, 'k-')
plt.grid(color='k', linestyle='dotted', linewidth=1)
plt.show()
コメント