发泡机 东莞网站建设,17173游戏排行榜,计算机应用网站建设与维护是做什么,优秀网站设计 pdf最近学习自动控制原理#xff0c;关于控制系统的一些#xff0c;老师用布置了一些作业说要用matlab画#xff0c;我试试python首先介绍一下所使用的库#xff1a;control matplotlib sympy1.control库#xff1a;用来计算脉冲响应与阶跃响应Paste_Image.png2.sympy#x…最近学习自动控制原理关于控制系统的一些老师用布置了一些作业说要用matlab画我试试python首先介绍一下所使用的库control matplotlib sympy1.control库用来计算脉冲响应与阶跃响应Paste_Image.png2.sympy用以化简多项式为和的形式方便写参数比如我们使用 sympy 来验证等式x2y2(x(2xy)1/2y)(x−(2xy)1/2y)from sympy import Symbol, expandfrom fractions import Fractionx 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/499794113.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 ctlimport matplotlib.pyplot as plt#画图包from sympy import Symbol, expanddef step_plot(s,color):#y,tctl.step(s)plt.plot(t,y,color,linewidth0.6)plt.xlabel(Time(seconds),fontsize9)plt.ylabel(Amplitude,fontsize9)return pltdef impulse_plot(s,color):#y,tctl.impulse(s)plt.plot(t,y,color,linewidth0.6)plt.xlabel(Time(seconds),fontsize9)plt.ylabel(Amplitude,fontsize9)return pltsSymbol(s)ySymbol(y)listexpand((s**20.6*s1)*(s**23*s9)*(s5))sctl.tf([45],[1,8.6,29.8,67.4,51,45])bctl.tf([1],[1,0.6,1])ax1plt.subplot(311)#在图表2中创建子图1ax1impulse_plot(s,red)ax2plt.subplot(312)#在图表2中创建子图2ax2impulse_plot(b,blue)ax3plt.subplot(313)ax1impulse_plot(s,red)ax2impulse_plot(b,-.)plt.show()
