% 2. järku mudeli identifitseerimine
% Hüppelise sisendi 0...500 KATSEANDMED:
% ajavektor  t [h]
t=[0 0.5 1   2    3    4    5    6  10]'
% ja nendele ajahetketele vastavad väljundi väärtused
y=[8 9.0 11 13.6 14.8 15.2 15.4 15.45 15.5]'
plot(t,y)% katseandmete graafiku väljatrükk

% Staatiline ülekandeteguri Ko määramine:
deltaU=500 % sisendi muutus
deltaY=15.5 - 8 % väljundi muutus:
K0=deltaY/deltaU

% väljundi karakteristiku  normeerimine:
yy=y-8
yy=yy/7.5

plot(t,yy)% normeeritud karakteristiku väljatrükk

%abimuutujate leidmine:
t7=spline(yy,t,0.7)
t3=t7/3
y3=spline(t,yy,t3)
% vastavalt y3-le leiame tabelist hii_ruudu
y3n=[0.189 0.20 0.21 0.22 0.24 0.260 0.28 0.30 0.32];
hi2=[0.000 0.11 0.20 0.27 0.41 0.525 0.63 0.75 0.87];
hii_ruut=spline(y3n,hi2,y3) 
plot(hi2,y3n), grid, hold on, plot(hii_ruut,y3,'o')
hii=sqrt(hii_ruut)
% ajakonstantide T1 T2 leidmine
T= t7/2.4
T1= T*(1+hii)
T2= T*(1-hii)