二阶系统不同参数的响应（Python）

最新推荐文章于 2024-03-12 23:01:27 发布

君王的羔羊

最新推荐文章于 2024-03-12 23:01:27 发布

阅读量1.6k

点赞数 2

本文链接：https://blog.csdn.net/GQ_Sonofgod/article/details/109161107

版权

自动控制原理专栏收录该内容

2 篇文章 0 订阅

订阅专栏

python计算系统的阶跃响应和脉冲响应
在这里插入图片描述

from scipy.signal import lti,step2,impulse2
import matplotlib.pyplot as plt

f,((ax1,ax2,ax3),(ax4,ax5,ax6),(ax7,ax8,ax9),(ax10,ax11,ax12)) = plt.subplots(4,3,sharex='col',sharey='row') # 开启subplots模式
for i in range(1,13,1):
    exec("s%s=lti([%d],[1,2,5])"%(i,i))
    exec("t%s,y%s=step2(s%s)"%(i,i,i))
    exec("ax%s.plot(t%s,y%s,'r',label='s%s Step Response',linewidth=0.5)"%(i,i,i,i))
    print (i)
plt.show()

$$
G(s)={4\over s^3+2s^2+10s+8}
$$

$G(s)={4\over s^3+2s^2+10s+8}$

import numpy as np
import control as ctl
import matplotlib.pyplot as plt

def step_plot(s):
    t,y=ctl.step_response(s)
    plt.plot(t,y,'b',linewidth=0.6)
    plt.title('Step Response',fontsize=9)
    plt.xlabel('Time(seconds)',fontsize=9)
    plt.ylabel('Amplitude',fontsize=9)
    plt.show()

def impulse_plot(s):
    t,y=ctl.impulse_response(s)
    plt.plot(t,y,'b',linewidth=0.6)
    plt.title('Impulse Response',fontsize=9)
    plt.xlabel('Time(seconds)',fontsize=9)
    plt.ylabel('Amplitude',fontsize=9)
    plt.show()


s=ctl.tf([4],[1,2,10,8])
step_plot(s)

在这里插入图片描述
`
伯德图绘制（幅频曲线magnitude和相频曲线phase）

from scipy import signal
import matplotlib.pyplot as plt

s1 = signal.lti([2.557e-9,1], [4.999e-15,1.268e-8, 1])
w, mag, phase = s1.bode()
print(w)
plt.figure()
plt.semilogx(w, mag)    # Bode magnitude plot
plt.figure()
plt.semilogx(w, phase)  # Bode phase plot
plt.show()