GV

GeezerVison

The World Wide Web, restored to 800 pixels of glory.

Home|About|Gaming cabinet|Archives|How it works|Starter sites

On this site

  • Sign in
  • ActionsCI/CD and automation
  • CodespacesCloud development environments
  • All Features
  • Documentation
  • Training
  • Certifications
  • Build with Copilot
  • Customize Copilot
  • Integrations
  • Trending repos
  • GitHub Sponsors
  • Developer articles
  • Enterprises

This view

  • Original page
  • New address

GeezerVison is showing github.com — python-control/examples/sisotool_example.py at main · ControlCorp/python-control

Expand file tree

/

sisotool_example.py

Copy path
More file actions
More file actions

Latest commit

 

History

History
History
46 lines (36 loc) · 802 Bytes
/

sisotool_example.py

Copy path

File metadata and controls

46 lines (36 loc) · 802 Bytes
Raw
Copy raw file
Download raw file
Open symbols panel
Edit and raw actions
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""sisotooldemo.py
Shows some different systems with sisotool.
All should produce smooth root-locus plots, also zoomable and clickable,
with proper branching
"""
#%%
import matplotlib.pyplot as plt
from control.matlab import *
# first example, aircraft attitude equation
s = tf([1,0],[1])
Kq = -24
T2 = 1.4
damping = 2/(13**.5)
omega = 13**.5
H = (Kq*(1+T2*s))/(s*(s**2+2*damping*omega*s+omega**2))
plt.close('all')
sisotool(-H)
#%%
# a simple RL, with multiple poles in the origin
plt.close('all')
H = (s+0.3)/(s**4 + 4*s**3 + 6.25*s**2)
sisotool(H)
#%%
# a branching and emanating example
b0 = 0.2
b1 = 0.1
b2 = 0.5
a0 = 2.3
a1 = 6.3
a2 = 3.6
a3 = 1.0
plt.close('all')
H = (b0 + b1*s + b2*s**2) / (a0 + a1*s + a2*s**2 + a3*s**3)
sisotool(H)
You are visitor 01402520 — GeezerVison does not store the pages you look up.
Best viewed with Netscape Navigator 4.7 or Internet Explorer 6 at 800×600. No cookies required.