Matlab 单位阶跃响应曲线三维图、上升时间、峰值时间、最大过调量和调整时间

168 阅读 0 评论 111 点赞

我是靠谱客的博主秀丽白猫，这篇文章主要介绍Matlab 单位阶跃响应曲线三维图、上升时间、峰值时间、最大过调量和调整时间，现在分享给大家，希望可以做个参考。

1.Matlab作单位阶跃响应曲线的三维图

定义闭环系统传递函数如下：
H(s)=1/(s^2+2*zeta*s+1)

复制代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
%Matlab作单位阶跃响应曲线的三维图
%标准二阶系统响应曲线 wn=1
t=0:0.2:10;
zeta=[0 0.2 0.4 0.8 1];
for n=1:6
num=1;
den=[1 2*zeta(n) 1];
[y(1:51,n),x,t]=step(num,den,t) ;
end


plot(t,y)


mesh(t,zeta,y')

2.用Matlab求上升时间、峰值时间、最大过调量和调整时间

定义闭环系统传递函数如下：
H(s)=25/(s^2+6*s+25)

复制代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
%求上升时间、峰值时间、最大过调量和调整时间
%可以应用到高阶系统
num=[0 0 25];
den=[1 6 25];
t=0:0.005:5;
[y,x,t]=step(num,den,t) ;
r=1;
while y(r)<1.0001;r=r+1;end;
rise_time=(r-1)*0.005


[ymax,tp]=max(y);
peak_time=(tp-1)*0.005


max_overshoot=ymax-1


s=1001;while y(s)>0.98&y(s)<1.02;s=s-1;end;
settle_time=(s-1)*0.005

运行后结果：
rise_time = 0.5550
peak_time = 0.7850
max_overshoot = 0.0948
settle_time = 1.1850

求单位脉冲响应的两种方法（求斜坡响应也可用此方法）

在初始条件为零时，H(s)的单位脉冲响应也就是sH(s)的单位阶跃响应
H(s)=1/（s^2+0.2*s+1）

num=[1];
den=[1 0.2 1]；
impulse(num,den)

或

num=[0 1 0];
den=[1 0.2 1]；
step(num,den)