指标

分析方法

峰值时间tp

[y,t]=step(sys);%y为输出变量响应值，t为仿真时间向量

[Y,k]=max(y);%Y为y的峰值，k为响应时间

tp=t(k); %获取峰值时间

超调量σ%

C=dcgain(sys);%求取系统的终值

[Y,k]=max(y);%Y为y的峰值，k为响应时间

sita=100*(Y-C)/C;%计算超调量

上升时间tr

C=dcgain(sys);%求取系统终值

  n=1;

while y(n)<C

  n=n+1

  end%通过循环求取第一次到达终值时的时间

tr=t(n);%获得上升时间

C=dcgain(sys);%求取系统终值

  n=1;

while y(n)<0.1*C

  n=n+1

  end

  m=1;

while y(m)<0.9*C

  m=m+1

  end

tr=t(m)-t(n);%获得上升时间

调节时间ts

[y,t]=step(sys);%y为输出变量响应值，t为仿真时间向量

C=dcgain(sys);%求取系统终值

i=length(t);%求仿真时间向量长度

while(y(i)>0.98*C)&(y(i)<1.02*C)%

i=i-1

end%循环求取终值在±2%内的时间

ts=t(i);%调节时间