control system using anaconda

 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 signal
import matplotlib.pyplot as plt
s1 = signal.lti([1], [1, 1])
w, mag, phase = signal.bode(s1)
plt.figure()
plt.semilogx(w, mag)    # bode magnitude plot
plt.figure()
plt.semilogx(w, phase)  # bode phase plot
plt.show()