最近学习自动控制原理,关于控制系统的一些,老师用布置了一些作业说要用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 已知响应函数求单位阶跃响应或脉冲响应...相关推荐

  1. python函数库 阶跃 信号函数 调用_python 已知响应函数求单位阶跃响应或脉冲响应...

    最近学习自动控制原理,关于控制系统的一些,老师用布置了一些作业说要用matlab画,我试试python 首先介绍一下所使用的库:control matplotlib sympy 1.control库: ...

  2. python已知y求x_当使用python知道X坐标值未知时,如何根据Y坐标获取Path的X坐标

    当我只有一个已知的Y坐标方程,即P = a * b(其中a& b定义的值为0.8,150)和x坐标完全未知时,如何在曲线绘图上得到一个点没有连接x和y的方程(例如:y = mx b; #i没有 ...

  3. python已知两边求第三边_探究“已知一个三角形两边及其夹角,求第三边”的问题...

    探究"已知一个三角形两边及其夹角求第三边"的问题 [知识点] 余弦定理. [对应版本章节] 本节课是人民教育出版社出版的高中数学 A 版数学必修 5 第一章"解三角 形& ...

  4. python已知两边求第三边_已知两边求第三边公式

    如果是三角形是直角三角形,知道两边,可以用勾股定理求出第三边.如果是三角形是普通三角形(锐角.钝角三角形),那这个条件下只能求出第三边的范围:两边之和大于第三边,两边之差小于第三边. 求边公式 只知道 ...

  5. 二阶系统单位阶跃响应 c语言,已知二阶系统的单位阶跃响应为

    习 题 3-1 已知二阶系统的单位阶跃响应为 )1.536.1sin(5.1210)(2.1°+?=?t e t c t 试求系统的超调量σ%.峰值时间p t 和调节时间s t %)2(±=Δ. 3- ...

  6. 已知函数和y值求自变量用matlab,matlab已知函数值求自变量,求助

    思路是这样的:你得到最小的函数值后,把它带入你作图的方程反解自变量的值.我不知道你画的这个函数图像的方程什么,我用一个二次方程y=x.^2+x给你做示范:先做图:x=-10:0.1:10;y=x.^2 ...

  7. python怎么根据点来拟合曲线_2019_nCoV_利用python根据已知点求拟合曲线及简单预测(无实际意义)...

    前言 本文仅做根据已知点求拟合曲线的几种方法的python实现,无任何实际意义 数据来源(另一篇博文) 利用Python爬取新冠肺炎疫情实时数据,Pyecharts画2019-nCoV疫情地图 参考 ...

  8. 已知序列求蝶形运算_在N=32的基 2 时间抽取发 FFT 运算流图中,从 到 需______级蝶形运算过程。 (4.0分)_学小易找答案...

    [单选题]互相关函数性质正确的是: [计算题]两个有限长的复序列 x[n] 和 h[n] ,其长度分别为 N 和 M ,设两序列的线性卷积为 y[n]=x[n]*h[n] ,回答下列问题: (1) 序 ...

  9. C语言--已知三边求三角形面积

    已知三边求面积用海伦公式S=sqrt[p(p-a)(p-b)(p-c)] p=(a+b+c)/2 #include<stdio.h> #include<math.h> int ...

  10. 已知分布函数求概率密度例题_高中数学必修一函数知识点总结

    同学们,今天开始讲解函数章节学习,函数这章极其重要,因为函数是高中数学重要的枢纽章节,高中数学除了立体几何和概率统计和函数没有关系之外,所有章节多多少少和函数有关系,所以函数学不好高中数学很难突破10 ...

最新文章

  1. Spring-ApplicationContext解读
  2. 学委收作业神器-查未交同学
  3. python怎么输入文件数据库_python学习日记——文件及数据库
  4. Skywalking部署常见问题以及注意事项
  5. hbuid 集成svn_HBuilder如何配置SVN的步骤详解
  6. BUS HOUND调试USB驱动遇到的错误代码解析
  7. CentOS 6.4 安装 media wiki 1.23.6(转)
  8. postgre管理员 无法访问表_PostgreSQL常见问题处理方法
  9. 什么是微服务(通俗易懂)
  10. 手机怎么设置企业邮箱的服务器,如何在手机上设置网易企业邮箱
  11. Java 行业的相关就业前景分析
  12. 新年寄语 —— 奋斗2020
  13. linux 端口耗尽,短连接以致端口耗尽
  14. 从MapX到MapXtreme
  15. 技校计算机英语,技校计算机英语的一体化教学论文
  16. 上标及下标 Unicode
  17. October 2006
  18. 单片机 李广弟_单片机基础(第3版)
  19. PCB和软件设计经验(1)
  20. 为什么安监控需要公网ip_三种方法告诉你项目超过255个摄像机怎么设置IP?

热门文章

  1. 无插件使用Eclipse和Resin调试WEB应用(Debug Web App In Ecli...
  2. 英文名称:DSPE-PEG10-Mal的试剂分子式是C68H127N2O21P
  3. DedeCMS 标签使用方法汇总!
  4. 算法专题训练(1)股票问题
  5. 大规模网络拓扑测量与分析研究1-论文
  6. python中sklearn.datasets.make_blobs()函数用法
  7. [PTA] 第四章 题解
  8. 翻译2Momentum, business cycle, and time varying expected returns(后附读后感)(图表上传不了)
  9. 酷博社区卫生系统激活健康档案的生命活力
  10. 重力传感器、加速度传感器以及陀螺仪的区别