% % KritIdent konstrukce minimalizovane ucelove funkce"KritIdent") % vektor.parametr x obsahuje koeficienty a2,a1 prenos.fce % verze Matlab 7 a vyssi function f=KritIdent(x) global a1 a2 Tmax k a1v a2v Jv a1=x(1); a2=x(2); % naplneni vektoru x koeficienty prenosove funkce figure(2); plot(a1,a2,'b x','LineWidth',2,'MarkerSize',10); hold on; sim('Ident',Tmax); % spousti "Ident"(simulaci s aktualnimi koef.) % nutno volat z aktualniho adresare!) figure(1) plot(Time,ys,'b','LineWidth',1); hold on;% vykresli prechod.fci v aktual.kroku k=k+1; a1v(k)=a1; a2v(k)=a2; Jv(k)=J; pause % stiskni "Enter" pro pokracovani vypoctu f=J; % vraci hodnotu kriterialniho ohodnoceni vypoctenou simulaci