Dla inercji II-go rzędu z opóźnieniem.


te=ScopeData.time;
ue=ScopeData.signals(1).values;
ye=ScopeData.signals(2).values;
subplot(211)
plot(te,ue);grid
subplot(212)
plot(te,ye);grid

n=max(find(te<=step time));
t=te(n:end)-te(n);
u=ue(n:end)-ue(n-1);
y=ye(n:end)-ye(n);
plot(t,y);grid

v=diff(y)./diff(t);
tv=t(1:end-1);
[vV I]=max(v);
tpp=t(I);
ypp=y(I);

yu=max(y);

F%=ypp/yu*100%; (F%>26%)

U=max(u);

k=yu/U;

n10=max(find(y<=0.1*yu));
y10=y(n10);
t10=t(n10);

n90=max(find(y<=0.9*yu));
y90=y(n90);
t90=t(n90);

T=(t90-t10)/3.3;

tau=t10-0.53*T;

tr=7.2*tau;