传递函数

创建传递函数有两种方式

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 sympyclass Time_Response:''' 时域响应'''s, t = sympy.symbols('s, t')def __init__(self, fun, doprint=False):''' fun: 返回关于s的传递函数的 functiondoprint: 输出公式'''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_Responseorange = '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 求解的响应曲线：一毛一样

Python 自动控制原理 control相关推荐

python自动控制原理_自动控制原理（山东联盟-中国石油大学（华东））知到答案2020年MOOCPython语言程...
自动控制原理(山东联盟-中国石油大学(华东))知到答案2020年MOOCPython语言程序设计测试答案更多相关问题设有一半径为a的介质圆球,置于一均匀磁场B之中,且绕通过其球心的某一固定轴,以角 ...
Auto Control 001 自动控制原理自动控制的一般概念
一 . 自动控制系统的组成自动控制装置:自动控制装置的组成当中涉及到了这样一些东西: 需要有(1)被控对象,那么这些被控对象需要有谁来控制呢?一定要有控制器,这些控制器,我们也把它叫做(2)自动控制 ...
使用python-control库实现MATLAB自动控制原理常用函数：Bode图 Nyquist图根轨迹
最近在学习自动控制原理,电脑装的windows/ubuntu双系统,matlab安装在windows上,日常用ubuntu,写作业切换系统过于繁琐,因此想要找到python中可以替换MATLAB Co ...
python安装插件报错原因_Sublime Text3 python自动补全问题——Sublime Text3安装Anaconda插件...
学习python的时候在编辑器的选择上会有很多选择,我最终还是选择了sublime text3. 相对于其他编辑器,sublime text有以下特性: 插件多,类似GoSublime,Emmet信 ...
python自动化办公能做什么-用Python自动办公，做职场高手（完结)
教程目录: ┣━07.S2 Word自动化处理,又快又好做文档 ┃ ┣━36 本章介绍 ┣━08.[Word]S2-1 轻松用Python快速生成Word文档 ┃ ┣━45.[真实案例]S2-1-3 ...
《狂飙》壁纸太帅，Python自动切换太酷(8)
小朋友们好,大朋友们好! 我是猫妹! 要说最近什么电视剧最火? 非<狂飙>莫属. <狂飙>剧名来自毛主席诗词"国际悲歌歌一曲,狂飙为我从天落". 导演借用& ...
自动控制原理matlab pdf下载,第五章-matlab在自动控制原理中的应用.pdf
第五章-matlab在自动控制原理中的应用.pdf 第第五五章章 MATLAB MATLAB MATLAB MATLAB在自动控制原理的应用在自动控制原理的应用 5 5 5 5 1 1 1 1 控制系 ...
Python自动生成代码 - 通过tkinter图形化操作生成代码框架
Python自动生成代码 - 通过tkinter图形化操作生成代码框架背景脚本代码 Demo_CodeGenerator.py display.py FileHandler.py: 脚本运行结果: ...
自动控制原理-＞绪论
绪论历史发展英文名词课程目标研究历史经典控制的相关概念历史发展基于国防科技大学学堂在线课程<自动控制原理>. 英文名词反馈feedback 动态过程dynamic proc ...

Python 自动控制原理 control

传递函数

输入响应

系统绘图

Laplace 逆变换

Python 自动控制原理 control相关推荐

最新文章

热门文章