ste.m

<MATLAB 7辅助控制系统设计与仿真>书籍内的源代码.

function [sigma,tp,ts]=ste(y,t) %函数定义
 [mp,tf]=max(y);
cs=length(t);
yss=y(cs);
sigma=(mp-yss)/yss    %计算超调量
tp=t(tf)        %计算峰值时间
%计算调节时间
i=cs+1;n=0;
while n==0,
      i=i-1;
        if i==1,
            n=1;
        elseif y(i)>1.05*yss,
            n=1;
        end;
end;
t1=t(i);cs=length(t);j=cs+1;
n=0;
while n==0,
    j=j-1;
    if j==1,
        n=1;
    elseif y(j)<0.95*yss,   %选择5%的误差带
        n=1;
    end;
end;
t2=t(j);
if t2<tp
    if t1>t2
        ts=t1
    end
elseif t2>tp,
    if t2<t1,
        ts=t2
    else
        ts=t1
    end
end