Rezultati

 

Slika 1 prikazuje simulacijsko shemo. Za podane vrednosti koeficientov oziroma velicin je bilo potrebno optimmirati parameter tako, da je prenihaj sistemskega toka ostal znotraj 5% stacionarne vrednosti. Vpliv na lomno frekvenco je podan z zakonitostjo

Za nas primer je bil podan barov. Z rocnim poskusanjem je bilo ugotovljeno obmocje za , ki je bil uporabljen v funkciji fmin(sfuns, p1, p2).

Optimizacijska funkcija je naslednja:

 function yo=optim(xo)
 global om3;
 om3=xo;
 [t,x,y]=gear('nek4e',1,[],[1e-4,.001,.1]);
 ymax=max(y);
 stac=100.3345;
 yo=(ymax-1.05*stac)^2;

Glavni program, ki je poiskal minimalni in s tem tudi optimalni je:

F0=3000;Ke=1800;
omm=35;Kq=9;
Kce=.04;om3=1;
A=2;K1=0.05;
Bef=12000;

global om3;
p=fmin('optim',1,2);
Vt=Bef*Kce/p;

Program je vrnil optimalno vrednost , s katerim izracunamo kot

Slika 2 kaze optimaliziran odziv na katerem je razviden 5% prenihaj sistema.

 

Slika 2: Odziv pri