%minimalizace kvadraticke chyby pri identifikaci kmitavych soustav
%funkce minim([k,T,ksi],h,t) vypocita hodnoty prechodove charakteristiky
%systemu F2=tf(k,[T*T T*ksi*2 1]) v casech danych vektorem t a tyto hodnoty porovna
%s hodnotami ve vektoru h. Vraci kvadratickou chybu, kterou algoritmus fminsearch 
%minimalizuje a tim identifikuje parametry soustavy F2.

par = fminsearch('minim',[6,6,1],[],h,t) %[k,T,ksi]=[6,6,1]
k = par(1); %ziskani hodnot identifikovanych parametru
T=par(2);
ksi=par(3);
F2=tf(k,[T*T T*ksi*2 1]) % identifikovana soustava