matlab单位阶跃响应与单位脉冲响应,python 已知响应函数求单位阶跃响应或脉冲响应...
最近学习自动控制原理,关于控制系统的一些,老师用布置了一些作业说要用matlab画,我试试python
首先介绍一下所使用的库:control matplotlib sympy
1.control库:用来计算脉冲响应与阶跃响应
Paste_Image.png
2.sympy:用以化简多项式为和的形式方便写参数
比如我们使用 sympy 来验证等式:x2+y2=(x+(2xy)1/2+y)(x−(2xy)1/2+y)
from sympy import Symbol, expand
from fractions import Fraction
x = Symbol('x')
y = Symbol('y')
expand((x+(2xy)Fraction(1, 2)+y)(x-(2xy)Fraction(1, 2)+y))
x2 + y*2
参照详细例子:
http://blog.csdn.net/lanchunhui/article/details/49979411
3.matplotlib 用以画图
分别是单位脉冲响应和单位节约响应的图像
单位脉冲响应.png
单位阶跃响应.png
![)QYB1_(CM2XKGM6}$)19R]P.png](http://upload-images.jianshu.io/upload_images/5854570-07199bf87612f8e6.png?imageMogr2/auto-orient/strip%7CimageView2/2/w/1240)
源码如下:
import control as ctl
import matplotlib.pyplot as plt#画图包
from sympy import Symbol, expand
def step_plot(s,color):#
y,t=ctl.step(s)
plt.plot(t,y,color,linewidth=0.6)
plt.xlabel('Time(seconds)',fontsize=9)
plt.ylabel('Amplitude',fontsize=9)
return plt
def impulse_plot(s,color):#
y,t=ctl.impulse(s)
plt.plot(t,y,color,linewidth=0.6)
plt.xlabel('Time(seconds)',fontsize=9)
plt.ylabel('Amplitude',fontsize=9)
return plt
s=Symbol("s")
y=Symbol("y")
list=expand((s**2+0.6*s+1)*(s**2+3*s+9)*(s+5))
s=ctl.tf([45],[1,8.6,29.8,67.4,51,45])
b=ctl.tf([1],[1,0.6,1])
ax1=plt.subplot(311)#在图表2中创建子图1
ax1=impulse_plot(s,'red')
ax2=plt.subplot(312)#在图表2中创建子图2
ax2=impulse_plot(b,'blue')
ax3=plt.subplot(313)
ax1=impulse_plot(s,'red')
ax2=impulse_plot(b,'-.')
plt.show()
matlab单位阶跃响应与单位脉冲响应,python 已知响应函数求单位阶跃响应或脉冲响应...相关推荐
- python函数库 阶跃 信号函数 调用_python 已知响应函数求单位阶跃响应或脉冲响应...
最近学习自动控制原理,关于控制系统的一些,老师用布置了一些作业说要用matlab画,我试试python 首先介绍一下所使用的库:control matplotlib sympy 1.control库: ...
- python已知y求x_当使用python知道X坐标值未知时,如何根据Y坐标获取Path的X坐标
当我只有一个已知的Y坐标方程,即P = a * b(其中a& b定义的值为0.8,150)和x坐标完全未知时,如何在曲线绘图上得到一个点没有连接x和y的方程(例如:y = mx b; #i没有 ...
- python已知两边求第三边_探究“已知一个三角形两边及其夹角,求第三边”的问题...
探究"已知一个三角形两边及其夹角求第三边"的问题 [知识点] 余弦定理. [对应版本章节] 本节课是人民教育出版社出版的高中数学 A 版数学必修 5 第一章"解三角 形& ...
- python已知两边求第三边_已知两边求第三边公式
如果是三角形是直角三角形,知道两边,可以用勾股定理求出第三边.如果是三角形是普通三角形(锐角.钝角三角形),那这个条件下只能求出第三边的范围:两边之和大于第三边,两边之差小于第三边. 求边公式 只知道 ...
- 二阶系统单位阶跃响应 c语言,已知二阶系统的单位阶跃响应为
习 题 3-1 已知二阶系统的单位阶跃响应为 )1.536.1sin(5.1210)(2.1°+?=?t e t c t 试求系统的超调量σ%.峰值时间p t 和调节时间s t %)2(±=Δ. 3- ...
- 已知函数和y值求自变量用matlab,matlab已知函数值求自变量,求助
思路是这样的:你得到最小的函数值后,把它带入你作图的方程反解自变量的值.我不知道你画的这个函数图像的方程什么,我用一个二次方程y=x.^2+x给你做示范:先做图:x=-10:0.1:10;y=x.^2 ...
- python怎么根据点来拟合曲线_2019_nCoV_利用python根据已知点求拟合曲线及简单预测(无实际意义)...
前言 本文仅做根据已知点求拟合曲线的几种方法的python实现,无任何实际意义 数据来源(另一篇博文) 利用Python爬取新冠肺炎疫情实时数据,Pyecharts画2019-nCoV疫情地图 参考 ...
- 已知序列求蝶形运算_在N=32的基 2 时间抽取发 FFT 运算流图中,从 到 需______级蝶形运算过程。 (4.0分)_学小易找答案...
[单选题]互相关函数性质正确的是: [计算题]两个有限长的复序列 x[n] 和 h[n] ,其长度分别为 N 和 M ,设两序列的线性卷积为 y[n]=x[n]*h[n] ,回答下列问题: (1) 序 ...
- C语言--已知三边求三角形面积
已知三边求面积用海伦公式S=sqrt[p(p-a)(p-b)(p-c)] p=(a+b+c)/2 #include<stdio.h> #include<math.h> int ...
- 已知分布函数求概率密度例题_高中数学必修一函数知识点总结
同学们,今天开始讲解函数章节学习,函数这章极其重要,因为函数是高中数学重要的枢纽章节,高中数学除了立体几何和概率统计和函数没有关系之外,所有章节多多少少和函数有关系,所以函数学不好高中数学很难突破10 ...
最新文章
- Spring-ApplicationContext解读
- 学委收作业神器-查未交同学
- python怎么输入文件数据库_python学习日记——文件及数据库
- Skywalking部署常见问题以及注意事项
- hbuid 集成svn_HBuilder如何配置SVN的步骤详解
- BUS HOUND调试USB驱动遇到的错误代码解析
- CentOS 6.4 安装 media wiki 1.23.6(转)
- postgre管理员 无法访问表_PostgreSQL常见问题处理方法
- 什么是微服务(通俗易懂)
- 手机怎么设置企业邮箱的服务器,如何在手机上设置网易企业邮箱
- Java 行业的相关就业前景分析
- 新年寄语 —— 奋斗2020
- linux 端口耗尽,短连接以致端口耗尽
- 从MapX到MapXtreme
- 技校计算机英语,技校计算机英语的一体化教学论文
- 上标及下标 Unicode
- October 2006
- 单片机 李广弟_单片机基础(第3版)
- PCB和软件设计经验(1)
- 为什么安监控需要公网ip_三种方法告诉你项目超过255个摄像机怎么设置IP?
热门文章
- 无插件使用Eclipse和Resin调试WEB应用(Debug Web App In Ecli...
- 英文名称:DSPE-PEG10-Mal的试剂分子式是C68H127N2O21P
- DedeCMS 标签使用方法汇总!
- 算法专题训练(1)股票问题
- 大规模网络拓扑测量与分析研究1-论文
- python中sklearn.datasets.make_blobs()函数用法
- [PTA] 第四章 题解
- 翻译2Momentum, business cycle, and time varying expected returns(后附读后感)(图表上传不了)
- 酷博社区卫生系统激活健康档案的生命活力
- 重力传感器、加速度传感器以及陀螺仪的区别