16.01.2015 Views

OPTIM´ALN´I ROZHODOV´AN´I AˇR´IZEN´I

OPTIM´ALN´I ROZHODOV´AN´I AˇR´IZEN´I

OPTIM´ALN´I ROZHODOV´AN´I AˇR´IZEN´I

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

KAPITOLA 2. NELINEÁRNÍ<br />

PROGRAMOVÁNÍ 26<br />

5. Vypočtěte optimální konstantu r 0 proporcionálního regulátoru, který v regulačním<br />

obvodu reguluje systém s přenosem<br />

G(s) =<br />

5<br />

s(s + 1) 5 .<br />

Vstupní signál je skok řízení w(t) = 1(t). Kritérium kvality řízení je<br />

J(r 0 ) =<br />

∫ ∞<br />

0<br />

(<br />

e 2 (t) + 3u 2 (t) ) dt<br />

kde e(t) je regulační odchylka a u(t) je akční veličina.<br />

6. Mějme stejný regulační obvod jako v předchozím příkladě. Řízení je nulové, ale na<br />

vstupu systému působí porucha, kterou můžeme považovat za náhodný signál blízký<br />

bílému šumu. Navrhněte optimální konstantu proporcionálního regulátoru, při které<br />

je rozptyl regulační odchylky minimální.<br />

7. Je dán diskrétní systém popsaný stavovými rovnicemi<br />

x(k + 1) = Mx(k) + Nu(k)<br />

y(k) = Cx(k),<br />

který řídíme pomocí lineární zpětné vazby ve tvaru<br />

a) u(k) = −K 1 x(k)<br />

b) u(k) = −K 2 y(k).<br />

Vypočtěte optimální zpětnovazební matici K i v obou případech, která minimalizuje<br />

kritérium<br />

∞∑<br />

J(K) = y T (k)Qy(k) + u T (k)Ru(k).<br />

k=0<br />

Návod k řešení:<br />

Dosad’te rovnici regulátoru do stavové rovnice systému. Řešení této rovnice dosad ’te<br />

do kritéria. Použijte vztah x T Qx = tr(xx T Q). Optimální zpětnou vazbu určíme z<br />

podmínky ∂J = 0. Derivaci kritéria podle zpětnovazební matice K provedeme tak,<br />

∂K<br />

že stanovíme přírůstek kritéria J(K + ε∆K) a použijeme dvě následující tvrzení:<br />

1) Mějme maticovou funkci F(X) = (A + BX) k a ε je nekonečně malá veličina<br />

prvního řádu. Pak platí<br />

F (X + ε∆X) = (A + BX) k +<br />

k∑<br />

ε (A + BX) k−i−1 B∆X (A + BX) i + ∆(ε 2 )<br />

i=0<br />

2) Mějme maticovou funkci f(X) = tr (F(X)) a ε je nekonečně malá veličina prvního<br />

řádu. Platí-li<br />

f (X + ε∆X) = f (X) + ε tr (M(X)∆X) + ∆(ε 2 ),<br />

pak<br />

∂f<br />

∂X = MT (X)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!