Slika 1 prikazuje simulacijsko shemo. Za podane vrednosti koeficientov oziroma velicin je bilo potrebno optimmirati parameterre tako, da je prenihaj cim manjsi.
Optimizacijska funkcija je naslednja:
function [out]=itae(in) global ystat; global cl; global Kq; global Av; Av=in out=0; [t,x,y]=gear('sim5',100); myt=0; for q=1:size(t,1) if myt==0 if y(q)>ystat myt=t(q); end end end out=(max(y)-ystat)*(myt);
Glavni program, ki je poiskal optimalne parametre je imel naslednje konstante:
clear Fo=105; Kk=80; wnk=30; Kq=10; Av=15; Ak=0.7; wn=50; cn=0.25; n=25; Vth=1; cl=0.08; Vt=35; Bef=14000; Qb=1; %popravki wn=1; Kq=200; n=10; cl=1; Av=32.8; cn=1; ystat=150;
Slika 2 kaze optimaliziran odziv na katerem je razviden odziv sistema.