aspec wyklad234, PWR - Automatyka i Robotyka W10, semestr 4, Podstawy automatyki, Podstawy Automatyki


METODY PRZESTRZENI STANÓW

Sposoby przejścia między różnymi opisami

Opis obiektu w przestrzeni stanów nie jest jednoznaczny, tj. różne zestawy macierzy A, b, c mogą opisywać obiekt o tych samych własnościach. Powstaje pytanie jak wybrać wektor stanu mając opis obiektu w postaci transmitancji. Istnieje kilka standardowych sposobów przejścia, które pozwalają na uzyskanie tzw. kanonicznych postaci macierzy A, b, c (metoda bezpośrednia, równoległa, iteracyjna). Weźmy obiekt opisany transmitancją :

0x01 graphic

0x01 graphic

0x01 graphic
, gdzie

0x01 graphic
, czyli

0x01 graphic
.

Rysujemy schemat blokowy wewnętrznej struktury obiektu pamiętając, że mnożenie przez s-1 oznacza całkowanie.

0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
rn-1

0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
rn-2

...... ....

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
U(s) E(s) Y(s)

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
∫ ∫ ... ∫ r0

xn xn-1 x2 x1

0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
-pn-1

0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
-pn-2

.... . . . . .

0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
-p0

Na poszczególne zmienne stanu xi (t) wybieramy wielkości :

0x01 graphic

Po przejściu do dziedziny czasu otrzymuje się układ równań różniczkowych :

0x01 graphic

równanie wyjścia :

0x01 graphic

Wykorzystując zapis macierzowy :

0x01 graphic

0x01 graphic

W ten sposób otrzymaliśmy jedną z postaci kanonicznych Brunowsky'ego - Luenbergera macierzy A, b, c. Jest to tzw. postać kanoniczna dla sterowania. Stosując inny tok postępowania można otrzymać tzw. postać kanoniczną dla obserwowania.

0x01 graphic
; 0x01 graphic

0x01 graphic

Przejście odwrotne, tzn. z opisu przestrzeni stanów do transmitancji wykonuje się następująco :

0x01 graphic

Przejście pomiędzy równaniem różniczkowym a transmitancją wykonuje się za pomocą przekształcenia Laplace'a.

Sterowalność i obserwowalność

W przypadku opisu wektorem stanu, może zdarzyć się tak, że przy dostępnym sterowaniu nie jesteśmy w stanie wpływać na niektóre zmienne. Mówimy wtedy, że te zmienne nie są sterowalne. Analogicznie niektóre zmienne stanu mogą w ogóle nie wpływać na wyjście. Mówimy wtedy, że te zmienne nie są obserwowalne. Ogólnie zmienne stanu można podzielić na cztery kategorie.

x - zmienne stanu

0x08 graphic

sterowalne

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
i obserwowalne

0x08 graphic

niesterowalne y(t)

0x08 graphic
i obserwowalne

0x08 graphic

u(t) sterowalne

0x08 graphic
i nieobserwowalne

0x08 graphic

niesterowalne

i nieobserwowalne

Pojęcie sterowalności i obserwowalności prześledźmy na przykładzie obiektu jakim jest myślowy model zegara wskazówkowego. Wielkości wejściowe są wprowadzane do obiektu za pomocą pokręteł umieszczonych z tyłu. Wielkościami wyjściowymi jest położenie wskazówek na cyferblacie. Cały mechanizm zegara przetwarza położenie pokręteł na położenie wskazówek. Zmiennymi stanu są sprężyny znajdujące się wewnątrz zegara. Zmienne stanu reprezentują zmagazynowaną w obiekcie energię. Kółka zębate są członami proporcjonalnymi i stanowią fizyczne połączenie pomiędzy wejściem, zmiennymi stanu i wyjściem. Można je interpretować jako współczynniki macierzy. Zmienne stanu są sterowalne, jeżeli kręcąc pokrętłami mamy wpływ na daną sprężynę. Zmienne stanu są obserwowalne, gdy stan danej sprężyny ma wpływ na położenie wskazówek. Sterowalność i obserwowalność zależy od wewnętrznych połączeń w strukturze obiektu. Dla obiektów liniowych wykryto [Kalman] odpowiedni warunek matematyczny, pozwalający testować wewnętrzną strukturę obiektu.

Twierdzenia o sterowalności i obserwowalności

Układ liniowy stacjonarny o jednym wejściu jest sterowalny wtedy i tylko wtedy, gdy rząd macierzy sterowalności Q jest pełny, gdzie

0x01 graphic
,

czyli rz Q = n.

W praktyce sprawdza się warunek : czy det Q ≠ 0.

Układ liniowy stacjonarny o jednym wyjściu jest obserwowalny wtedy i tylko wtedy, gdy rząd macierzy obserwowalności N jest pełny.

0x01 graphic
, czyli rz N = n.

W praktyce sprawdza się warunek det N ≠ 0.

Obiekt zapisany w postaci kanonicznej do sterowania jest w pełni sterowalny. Obiekt zapisany w postaci kanonicznej do obserwowania jest w pełni obserwowalny.

Liniowa transformacja wektora stanu.

Obiekt jest opisany przy pomocy następujących równań stanu :

0x01 graphic

Chcemy by obiekt był opisany przy pomocy innego wektora stanu, zdefiniowanego w sposób następujący :

0x01 graphic

gdzie T - macierz transformacji, spełniająca następujące warunki :

dim T = n, co oznacza, że wektory x i v mają ten sam wymiar,

rz T = ndet T ≠ 0 - macierz jest nieosobliwa

W nowych współrzędnych obiekt jest opisany równaniami :

0x01 graphic

Transformacja T wektora stanu działa w sposób następujący na trójkę macierzy A, b, c.

0x01 graphic

Opis za pomocą wektora stanu nie jest jednoznaczny. Dwa obiekty mają taką samą dynamikę wtedy i tylko wtedy, gdy ich równania charakterystyczne maję te same pierwiastki. Liniowa transformacja wektora stanu nie zmienia położenia pierwiastków równania charakterystycznego.

Cele konstruowania układu regulacji automatycznej

0. Zapewnić stabilność

1. Spełnić wymogi, co do kształtu przebiegu przejściowego.

2. Spełnienie wymogów co do odchyłki statycznej.

3. Zapewnić sterowanie optymalne wg. określonego kryterium.

W punkcie 1. wymogi są określone za pomocą czasu regulacji, odchyłki dynamicznej, konieczności wyeliminowania zakłóceń.

Koncepcja realizacji celu.

Stosujemy proporcjonalne sprzężenie zwrotne od poszczególnych zmiennych stanu. Zakładamy, że wszystkie stany są dostępne pomiarowo.

0x08 graphic
0x08 graphic
x1

0x08 graphic
0x08 graphic
z(t) x2 y(t)

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x01 graphic
. . . . . . c

0x08 graphic
_ xn

0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
g0

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
g1

. . . ...

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
gn

Po wprowadzeniu sprzężenia zwrotnego, nowe sterowanie ma postać :

0x01 graphic
, gdzie

0x01 graphic
.

Idąc tym tropem otrzymujemy (dla uproszczenia zapisu pomijamy w zapisie zależność od czasu):

0x01 graphic

Wprowadzenie sprzężenia zwrotnego doprowadziło do transformacji T trójki macierzy opisujących obiekt sterowania.

0x01 graphic

Sformułowanie zagadnienia

Dane:

1. Własności obiektu sterowania w postaci wektora współczynników wielomianu charakterystycznego :

0x01 graphic

2. Pożądane własności układu regulacji (obiekt z odpowiednim sprzężeniem zwrotnym) - również w postaci wektora współczynników wielomianu charakterystycznego :

0x01 graphic

Pożądany wielomian charakterystyczny dobieramy tak, aby były spełnione warunki 0 i 1. Uzyskujemy to przez odpowiednie umiejscowienie biegunów.

3. Macierze A, b są w postaci kanonicznej do sterowania.

Szukane :

Znaleźć taki wektor g , aby układ sterowania miał pożądane własności określone wektorem 0x01 graphic
.

Rozwiązanie :

Z jednej strony zachodzi :

0x01 graphic

0x01 graphic

Z drugiej strony mamy :

0x01 graphic
,

zatem otrzymujemy układ równań :

0x01 graphic

i ostatecznie

0x01 graphic
.

Wynik ten jest słuszny, jeżeli macierze A,b,c są w postaci kanonicznej dla sterowania. W przeciwnym przypadku musimy jeszcze dokonać odpowiedniej transformacji liniowej - o czym potem.

Obserwator stanu

Rozważania przeprowadzone wcześniej są słuszne, gdy wszystkie zmienne stanu są dostępne pomiarowo. Często zdarza się jednak tak, że poszczególne zmienne nie są dostępne dla układu regulacji. Wynika to z faktu , że nie jesteśmy w stanie w sposób zadowalający wyznaczać pochodnych wyższych rzędów ze zmierzonych sygnałów. Obserwatory stanu konstruuje się właśnie w celu udostępnienia zmiennych stanu.

Wydawałoby się, że nie ma tu żadnego problemu. Wystarczyłoby działać sygnałem wejściowym u(t) na układ o takich samych własnościach dynamicznych jak obiekt sterowania. Jednakże w praktyce nie znamy warunków początkowych, w jakich znajduje się obiekt. Jak również zakłócenia mogą powodować rozbieżność między rzeczywistymi zmiennymi stanu a ich estymatami generowanymi przez obserwator. Stąd też zachodzi konieczność konstruowania takiego obserwatora, dla którego ta rozbieżność dążyłaby asymptotycznie do zera.

Na początku lat 70-tych Luenberger zaproponował następującą strukturę obserwatora:

0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
u(t) 0x01 graphic
c 0x01 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

x(t) y(t)

0x01 graphic
0x01 graphic
0x01 graphic
. . . . . . . 0x01 graphic

Naszym celem jest takie dobranie parametrów obserwatora (tzn. macierzy 0x01 graphic
) tak, aby : 0x01 graphic
, gdzie 0x01 graphic

Odejmując stronami równania dynamiki obiektu i obserwatora otrzymuje się :

0x01 graphic

0x01 graphic

Przyjmujemy, że macierze obserwatora są następujące :

0x01 graphic
0x01 graphic

czyli : 0x01 graphic

Otrzymaliśmy autonomiczny układ dynamiczny ze względu na e(t), którego rozwiązaniem jest (rozwiązanie jednorodnego układu równań różniczkowych) :

0x01 graphic

Jeżeli wartości własne macierzy 0x01 graphic
leżą w lewej półpłaszczyźnie, to spełniony jest warunek e(t) → 0. Wektor k dobieramy tak, by otrzymać jak najszybszą zbieżność.

-Sformułowanie zagadnienia

Dane :

1. wektor współczynników wielomianu charakterystycznego obiektu sterowania - p,

2. wektor współczynników wielomianu charakterystycznego obserwatora - 0x01 graphic
,

3. macierze A, c są w postaci kanonicznej dla obserwowania.

Szukane :

Tak dobrać wektor k, aby macierz 0x01 graphic
miała pożądane własności dynamiczne.

Rozwiązanie :

Z jednej strony mamy :

0x01 graphic

Prowadzi to do równości :

0x01 graphic

Dualność

Zagadnienia konstruowania regulatora i obserwatora są dualne. Podsumujmy uzyskane dotychczas rezultaty.

DANE :

Opis obiektu sterowania w przestrzeni stanów za pomocą macierzy A,b,c. Własności dynamiczne obiektu są określone przez współczynniki wielomianu charakterystycznego :

0x01 graphic
: 0x01 graphic

REGULATOR

OBSERWATOR

CEL

Tak dobrać wektor g , aby układ sterowanie miał pożądane własności dynamiczne :

0x01 graphic

Tak dobrać wektor k, tak aby obserwator miał pożądane własności dynamiczne :

0x01 graphic

TWIERDZENIE 1

1. Jeżeli obiekt jest w pełni sterowalny oraz

2. Macierze A, b są w postaci kanonicznej dla sterowania,

TO : 0x01 graphic

1. Jeżeli obiekt jest w pełni obserwowalny oraz

2. Macierze A, c są w postaci kanonicznej dla obserwowania,

TO : 0x01 graphic

TWIERDZENIE 2

Jeżeli obiekt jest w pełni sterowalny, ale macierze A, b są w dowolnej postaci

TO : 0x01 graphic

0x01 graphic

Jeżeli obiekt jest w pełni obserwowalny, ale macierze A, c są w dowolnej postaci,

TO : 0x01 graphic

0x01 graphic

Macierz W jest utworzona ze współczynników wielomianu charakterystycznego obiektu sterowania :

0x01 graphic

Algorytm konstruowania układu regulacji

1. Wyznacz współczynniki wielomianu charakterystycznego dla obiektu sterowania, pamiętając by pn = 1.

2. Przedstaw obiekt w postaci kanonicznej dla obserwowania, sprawdź sterowalność (obserwowalność wynika z przyjętej postaci).

Jeżeli niedostępne są zmienne stanu to :

3. Określ pożądane własności dynamiczne obserwatora podając położenia biegunów :0x01 graphic
.

4. Wyznacz współczynniki wielomianu charakterystycznego obserwatora, korzystając z zależności :

0x01 graphic
.

5. Wyznacz : 0x01 graphic
i przyjmij 0x01 graphic
, 0x01 graphic

6. Określ pożądane własności dynamiczne układu sterowania podając położenie biegunów : 0x01 graphic
.

7. Wyznacz współczynniki wielomianu charakterystycznego układu sterowania :

0x01 graphic

8. Wyznacz 0x01 graphic



Wyszukiwarka