conda install numpy scipy matplotlib # if not yet installed
conda install -c conda-forge control slycot
-----------------------------------------------------
# secord.py - demonstrate some standard MATLAB commands
# RMM, 25 May 09
import os
import matplotlib.pyplot as plt # MATLAB plotting functions
from control.matlab import * # MATLAB-like functions
# Parameters defining the system
m = 250.0 # system mass
k = 40.0 # spring constant
b = 60.0 # damping constant
# System matrices
A = [[0, 1.], [-k/m, -b/m]]
B = [[0], [1/m]]
C = [[1., 0]]
sys = ss(A, B, C, 0)
# Step response for the system
plt.figure(1)
yout, T = step(sys)
plt.plot(T.T, yout.T)
plt.show(block=False)
# Bode plot for the system
plt.figure(2)
mag, phase, om = bode(sys, logspace(-2, 2), plot=True)
plt.show(block=False)
# Nyquist plot for the system
plt.figure(3)
nyquist(sys)
plt.show(block=False)
# Root lcous plot for the system
rlocus(sys)
if 'PYCONTROL_TEST_EXAMPLES' not in os.environ:
plt.show()
output:
src::
https://github.com/python-control/python-control
https://github.com/python-control/python-control/tree/master/examples
https://python-control.readthedocs.io/en/0.9.1/intro.html
https://anaconda.org/conda-forge/control
create transfer function
https://python-control.readthedocs.io/en/0.9.1/generated/control.tf.html
>>> # Create a MIMO transfer function object >>> # The transfer function from the 2nd input to the 1st output is >>> # (3s + 4) / (6s^2 + 5s + 4). >>> num = [[[1., 2.], [3., 4.]], [[5., 6.], [7., 8.]]] >>> den = [[[9., 8., 7.], [6., 5., 4.]], [[3., 2., 1.], [-1., -2., -3.]]] >>> sys1 = tf(num, den)
ploting
mag, phase, omega = bode(sys1)
from scipy import signalimport matplotlib.pyplot as plt
s1 = signal.lti([1], [1, 1])w, mag, phase = signal.bode(s1)
plt.figure()plt.semilogx(w, mag) # bode magnitude plotplt.figure()plt.semilogx(w, phase) # bode phase plotplt.show()
Comments
Post a Comment