timepro.m

来自「文件是编辑好的m文件」· M 代码 · 共 26 行

function [pos,tr,ts,TP1]=timepro(y,t);
% 计算阶跃响应曲线的性能指标

% Copyright 2001-2002 ECUST. 
%   $Revision: 1.2 $  $Date: 2001/05/11 17:26:43 $

tend=t(length(t));tsp=t(2)-t(1);
[mp,ind1]=max(y);yss=mean(y(round(tend*0.9):tend));
pos=(mp-yss)*100/yss;tp=ind1;
[minp,ind2]=min(y(ind1+1:tend));
[mp1,ind3]=max(y(ind1+ind2+1:tend));
Rs=(mp-yss)/(mp1-yss);
TP=ind2+ind3;TP1=TP+ind1;
inda=find(y(1:ind1)>0.95*yss & y(1:ind1)<1.05*yss);
a=inda(max(find(y(inda)<yss)));b=inda(min(find(y(inda)>yss)));
tr=a+(b-a)*(yss-y(a))/(y(b)-y(a))-1;
tss=0.05;%0.02~0.05
i=length(t);n=0;
while n==0,i=i-1;
    if (i==1| y(i)>=(1+tss)*yss),n=1;end;
end;t1=i;
i=length(t);n=0;
while n==0,i=i-1;
    if (i==1 | y(i)<=(1-tss)*yss),n=1;break;end;
end;t2=i;
if t2>t1,ts=t2;else ts=t1;end;