stepchar.m

来自「求系统的时域性能指标的函数」· M 代码 · 共 19 行

function [pos,tr,ts,tp]=stepchar(bktfg0,delta)
[y,t]=step(bktfg0);
[mp,ind]=max(y);
dimt=length(t);
yss=y(dimt);
pos=100*(mp-yss)/yss;
tp=t(ind);
for i=1:dimt
   if y(i)<1,
   else
      tr=t(i);
      break;
   end;
end;
for i=1:length(y)
   if y(i)<=(1-delta)*yss|y(i)>=(1+delta)*yss
      ts=t(i);
   end;
end;