Python 自动控制原理 control的详细解说
作者:荷碧·TongZJ 时间:2023-04-17 09:53:45
传递函数
创建传递函数有两种方式:
import control as ctrl
# 方式 1
s = ctrl.tf('s')
sys = 100 / (s ** 2 + 10 * s + 100)
# 方式 2
sys = ctrl.tf([100], [1, 10, 100])
对 tf 这个类,它内置的方法可求解 零点、极点、特征参数、特征根
# 求零点
sys.zero()
# 求极点
sys.pole()
# 特征参数、特征根
sys.damp()
输入响应
# 阶跃动态指标
step_info(sys)
# 阶跃响应
t, response = step_response(sys, T)
# 脉冲响应
t, response = impulse_response(sys, T)
t, response = initial_response(sys, T)
t, response = forced_response(sys, T)
T 是响应的时间,可以是 float (即时间上限),也可以是数组
阶跃动态指标是 dict 类型,包括:'RiseTime', 'SettlingTime', 'SettlingMin', 'SettlingMax', 'Overshoot', 'Undershoot', 'Peak', 'PeakTime', 'SteadyStateValue'
系统绘图
# Nyquist图, 可传入列表
nyquist_plot(sys)
# Bode图, 可传入列表
bode_plot(sys)
# 根轨迹图
root_locus(sys)
绘图使用的是 matplotlib.pyplot,所以执行完函数后,要加上 plt.show() 才会显示图像
Laplace 逆变换
可能是我太弱找不到这个包的 Laplace 逆变换函数,也可能是这个包真的没有这个函数
于是我利用 sympy 这个包求解:定义时域响应这个类,__call__ 使其可以计算时间数组 (np.array) 的响应
import sympy
class Time_Response:
''' 时域响应'''
s, t = sympy.symbols('s, t')
def __init__(self, fun, doprint=False):
''' fun: 返回关于s的传递函数的 function
doprint: 输出公式'''
sys = fun(self.s)
self.f_t = sympy.integrals.inverse_laplace_transform(sys, s=self.s, t=self.t)
if doprint:
sympy.pprint(self.f_t)
def __call__(self, time):
''' 使自身可作为函数被调用'''
response = list(map(lambda i: float(self.f_t.subs({self.t: i})), time))
return np.array(response)
设置 doprint 为 True,则可以输出时域响应的方程 —— 但是问题在于,自动控制原理里面的 Laplace 变换是默认 F(s) 各阶导数的初始值均为 0 的,这个条件我没有办法加入到 sympy 的求解过程里,所以结果看起来就有些奇怪
import control as ctrl
import matplotlib.pyplot as plt
import numpy as np
# 自定义类所在的模块
from mod.math_model import Time_Response
orange = 'orange'
blue = 'deepskyblue'
plt.rcParams['font.sans-serif'] = ['Microsoft YaHei']
f_t = Time_Response(lambda s: 100 / (s ** 2 + 10 * s + 100) / s, doprint=True)
t = np.linspace(0, 1, 100)
plt.subplot(1, 2, 1)
plt.title('sympy 计算')
plt.plot(t, f_t(t), c=orange)
s = ctrl.tf('s')
sys = 100 / (s ** 2 + 10 * s + 100)
t, response = ctrl.step_response(sys, T=1)
plt.subplot(1, 2, 2)
plt.title('control 计算')
plt.plot(t, response, c=blue)
plt.show()
对比 sympy 和 control 求解的响应曲线:一毛一样
来源:https://hebitzj.blog.csdn.net/article/details/124029533
标签:Python,自动,控制,control
![](/images/zang.png)
![](/images/jiucuo.png)
猜你喜欢
Python代码需要缩进吗
2022-05-07 18:21:15
用Dreamweaver设计限时自动关闭网页
2009-07-10 13:13:00
![](https://img.aspxhome.com/file/UploadPic/20072/200723113358457s.jpg)
关于keras中卷积层Conv2D的学习记录
2022-07-16 17:33:16
![](https://img.aspxhome.com/file/2023/0/77890_0s.png)
Python中类和对象的绑定及非绑定方法详解
2023-03-13 12:27:32
MySQL多表查询与7种JOINS的实现举例
2024-01-12 23:34:32
![](https://img.aspxhome.com/file/2023/1/68581_0s.png)
Python实现对百度云的文件上传(实例讲解)
2022-04-16 03:15:46
![](https://img.aspxhome.com/file/2023/1/107681_0s.jpg)
重温Javascript继承机制
2011-07-04 12:17:23
python numpy--数组的组合和分割实例
2023-02-17 03:58:39
![](https://img.aspxhome.com/file/2023/2/87942_0s.jpg)
ubuntu lamp(apache+mysql+php) 环境搭建及相关扩展更新
2023-10-16 03:26:55
对numpy中数组转置的求解以及向量内积计算方法
2023-01-23 16:39:52
python获取引用对象的个数方式
2023-06-19 23:10:30
如何利用pygame实现简单的五子棋游戏
2023-04-12 17:48:03
![](https://img.aspxhome.com/file/2023/8/121018_0s.jpg)
python Canny边缘检测算法的实现
2023-04-23 20:41:55
![](https://img.aspxhome.com/file/2023/4/105724_0s.png)
python神经网络VGG16模型复现及其如何预测详解
2022-10-16 06:59:15
![](https://img.aspxhome.com/file/2023/0/123460_0s.png)
抽取oracle数据到mysql数据库的实现过程
2024-01-14 03:07:12
JS实现键值对遍历json数组功能示例
2024-04-10 10:52:40
![](https://img.aspxhome.com/file/2023/6/136886_0s.png)
用SQL语句添加删除修改字段、一些表与字段的基本操作、数据库备份等
2011-12-01 07:53:11
python lambda 表达式形式分析
2021-05-06 06:05:58
![](https://img.aspxhome.com/file/2023/0/71320_0s.png)
Python可视化神器pyecharts绘制仪表盘
2021-12-29 15:57:53
![](https://img.aspxhome.com/file/2023/7/102957_0s.png)
MySql8.0以上版本正确修改ROOT密码的方法
2024-01-22 12:24:40
![](https://img.aspxhome.com/file/2023/5/103545_0s.jpg)