10.01.2015 Views

Szabályozástechnika Matlab Gyakorlatok, Villamosmérnöki - Index of

Szabályozástechnika Matlab Gyakorlatok, Villamosmérnöki - Index of

Szabályozástechnika Matlab Gyakorlatok, Villamosmérnöki - Index of

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.

Állapotteres Leírás, Irányíthatóság, Megfigyelhetőség Hetthéssy Jenő, Bars Ruth, Barta András, 2005<br />

tartalmaz információt az x 3 állapotváltozóról (az x 3 állapotváltozó nem megfigyelhető). Az y kimenőjel<br />

irányítható, mivel az u bemenőjel az x 1 állapotváltozón keresztül befolyásolja.<br />

Ezeket a tulajdonságokat az eredeti alakból kiindulva is megállapíthatjuk a Kalman-féle irányíthatósági és<br />

megfigyelhetőségi mátrixok rangjának vizsgálatával.<br />

n−1<br />

A rendszer irányíthatósági mátrixa C = ⎡ ⎣ b Ab ... A b⎤<br />

⎦<br />

számítható a ctrb utasítással.<br />

» Co=ctrb(A,b)<br />

vagy az LTI struktúrával:<br />

» Co =ctrb(H)<br />

A rendszer teljesen állapotirányítható, ha<br />

o<br />

Co<br />

rangja megegyezik az állapotváltozók n számával.<br />

» rank(Co)<br />

Esetünkben a rendszer nem teljesen állapotirányítható, mivel rank( Co<br />

) = 2< n .<br />

A rendszer kimeneti irányítható, ha a coy<br />

= c*<br />

C mátrix rangja megegyezik a kimenőjelek számával.<br />

o<br />

» rank(c*Co)<br />

Mivel ez az érték 1, megegyezik a kimenőjelek számával, a rendszer kimeneti irányítható.<br />

⎡ c ⎤<br />

⎢ cA ⎥<br />

A rendszer megfigyelhetőségi mátrixa Ob<br />

= ⎢ ⎥<br />

⎢ ... ⎥<br />

⎢ n 1 ⎥<br />

⎣cA −<br />

⎦<br />

A MATLAB obsv utasításával számolva:<br />

» Ob=obsv(A,c)<br />

» Ob =obsv(H)<br />

» rank(Ob)<br />

A rendszer megfigyelhető, ha a megfigyelhetőségi mátrix rangja n. Mivel rank( Ob<br />

) = 2< n , a rendszer<br />

nem teljesen megfigyelhető.<br />

Határozzuk meg a rendszer átviteli függvényét.<br />

» [num,den]=ss2tf(A,b,c,d)<br />

Az átviteli függvény zérus-pólus alakban:<br />

» Hzpk=zpk(H)<br />

Figyeljük meg, hogy két zérus kiejt két pólust.<br />

» Hzpk=minreal(Hzpk)<br />

Az átviteli függvény alakban információt vesztünk a rendszer állapotairól. Az átviteli függvény csak az<br />

irányítható és megfigyelhető állapotváltozókról tartalmaz információt.<br />

A párhuzamos reprezentációt megkaphatjuk a részlettörtes felbontással is:<br />

» [num,den]=tfdata (H,’v’)<br />

» [r,p,k]=residue(num,den)<br />

0 0 1<br />

H()<br />

s = + +<br />

s+ 3 s+ 2 s+<br />

1<br />

63

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

Saved successfully!

Ooh no, something went wrong!