最近学习自动控制原理，关于控制系统的一些，老师用布置了一些作业说要用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()
