Równanie stanu (teoria układów dynamicznych)

Równania stanu – sposób na reprezentację modelu matematycznego układu dynamicznego (zwłaszcza układu automatycznej regulacji). Znajomość stanu układu daje bardzo wiele, ale jeszcze więcej wiemy o układzie, gdy znamy związki zmiennej stanu z innymi ważnymi zmiennymi. Dlatego w opisie układu (w jego modelu matematycznym) kluczową rolę odgrywa związek rządzący zachowaniem się zmiennej stanu, czyli równania stanu. Opis układu za pomocą równań stanu nazywany jest też czasami opisem w przestrzeni stanów lub modelem zmiennych stanu.

Sformułowanie równań stanu

edytuj

Szeroką klasę układów dynamicznych, których miarą zmiany procesu w czasie jest pochodna wektora stanu, a stan procesu dla   zależy tylko od stanu w chwili początkowej   oraz od wymuszenia   dla   można opisać równaniem   z warunkiem początkowym   Jest to tak zwane równanie stanu.

Zwykle nie wszystkie zmienne stanu są dostępne (bezpośrednio mierzalne), czyli wektor stanu   nie jest zarazem wektorem odpowiedzi układu. Do pełnego opisu potrzebne jest jeszcze równanie wiążące   z   i z wektorem wymuszenia   o postaci:   Równanie to nazywamy równaniem wyjścia układu. Równanie wyjścia jest równaniem algebraicznym i dlatego nie zawiera się w nim opis dynamiki układu.

W ten sposób otrzymaliśmy równania stanu, czyli równania różniczkowe pierwszego rzędu, które stanowią matematyczny model układu fizycznego określonego przez powiązane między sobą:

  • zmienne: wejściowe   wyjściowe   i zmienne stanu   oraz
  • macierz stanu   macierz wyjść   macierz wejść   oraz macierz przenoszenia   (macierz przenoszenia   zwana też macierzą sterowania bezpośredniego pojawia się w równaniach stanu tylko dla układów właściwych – to znaczy mających transmitancję właściwą; dla układów ściśle właściwych – to znaczy mających transmitancję ściśle właściwą – macierz   nie występuje w równaniach stanu).
 
Schemat blokowy równań stanu dla układu ciągłego. Człon całkujący (integrator) zachowuje się jak element magazynujący (zob. stan układu)

Związek ten w przypadku modelu dla układu liniowego, ciągłego, stacjonarnego, o wielu wejściach i wielu wyjściach, deterministycznego, o parametrach stałych, można zapisać w postaci:

 
 

(jak już wyżej wskazano, drugie z powyższych równań nazywa się równaniem wyjścia), gdzie

 
 
 
 
 

Dla przypadku braku wymuszeń zewnętrznych   mamy do czynienia z tak zwanym układem swobodnym (opisanym jednorodnym równaniem stanu). W przypadku gdy takie wymuszenia pojawiają się w równaniu stanu, nazywane jest ono też czasami niejednorodnym równaniem stanu.

Dla przypadku układu o jednym wejściu   i jednym wyjściu   podane powyżej równania stanu przybierają postać:

 
 

gdzie   są odpowiednimi wektorami, a   jest czynnikiem skalarnym.

Dla przypadku układu niestacjonarnego, podane na wstępie równania stanu przybierają postać:

 
 

Dla przypadku układu nieliniowego, podane na wstępie równania stanu przybierają postać:

 
 

gdzie   i   są funkcjami nieliniowymi. Jeśli funkcja   jest liniową kombinacja stanów i wejść, to równania te można zapisać w notacji macierzowej jak w zapisach wyżej.

Dla przypadku układu stochastycznego, podane na wstępie równania stanu przybierają postać:

 
 

gdzie dodatkowo ujęte zostały wektory zakłóceń stochastycznych (wektory szumu):   – wektor szumu przetwarzania,   – wektor szumu pomiarowego.

 
Schemat blokowy równań stanu dla układu dyskretnego. Jednokrokowy operator opóźnienia   wskazuje, że mamy do czynienia z członem opóźniającym.

Dla przypadku modelu dyskretnego (z czasem dyskretnym), podane na wstępie równania stanu przybierają postać:

 
 

gdzie:   oznacza dyskretną chwilę czasu.

W przypadku tzw. układów uogólnionych, podane na wstępie równania stanu (zwane też w tym kontekście równaniami stanu układu normalnego) mają postać uogólnioną:

 
 

gdzie   jest macierzą o elementach rzeczywistych (tzn.  ) a   jest tu tzw. semiwektorem stanu (deskryptorem). Cechą układu uogólnionego jest to, że macierz   może nie mieć pełnego rzędu, to znaczy rząd   =  

Układ uogólniony nazywany jest układem singularnym, jeśli rząd   W przypadku szczególnym, gdy   wyżej podany układ uogólniony jest singularny, jeżeli   (tzn.   jest macierzą osobliwą).

Rozwiązanie równań stanu

edytuj

Dla równań przedstawionych na wstępie rozwiązanie ogólne dane jest równaniem (jest to tak zwany wzór Cauchy-Bellmana):

 

gdzie   jest macierzą przejścia.

Niejednoznaczność opisu równaniami stanu

edytuj

W przeciwieństwie do opisu typu wejście-wyjście (w postaci odpowiedniego równania różniczkowego, transmitancji lub całki splotowej) opis równaniami stanu nie jest jednoznaczny. Oznacza to, że:

  • różnym postaciom opisu w przestrzeni stanów może odpowiadać jeden opis transmitancyjny, a z drugiej strony
  • dla danej transmitancji istnieje nieskończenie wiele opisów w przestrzeni stanów (tym niemniej spośród różnych możliwych sposobów wyboru zmiennych stanu kilka jest szczególnie ciekawych).

Wprowadzenie przekształcenia nieosobliwego zmiennych stanu (to znaczy takiej zmiany współrzędnych w przestrzeni stanów, że przejście od jednych współrzędnych do drugich jest wzajemnie jednoznaczne) spowoduje, że ten sam układ będzie mógł być opisany innymi zmiennymi stanu i równaniami o innej postaci – przy zachowaniu tych samych własności.

Przekształcenie współrzędnych w przypadku układu liniowego oznacza, że współrzędne stanu   są kombinacjami liniowymi nowych współrzędnych z (i na odwrót), np:

 
 
 
 

Nieosobliwość przekształcenia oznacza, że wyznacznik macierzy przekształcenia   jest różny od zera, co zapewni jednoznaczność zamiany zmiennych   na   i na odwrót.

Równania można zapisać w skrócie jako   lub   Zastępując w poprzednich równaniach wektor stanu   przez   otrzymuje się nowe równania:

 
 

Przy użyciu nowych zmiennych stanu   otrzymuje się równania podobnej postaci, przy czym w miejscu macierzy   pojawia się macierz   zamiast wektora   – wektor   zamiast wektora   – wektor  

Najważniejsze jest spostrzeżenie, że macierz   charakteryzuje cecha podobieństwa do macierzy   Przesądza to o identyczności właściwości dynamicznych układu niezależnie od wyboru współrzędnych stanu.

Nieosobliwe przekształcenie współrzędnych stanu nie zmienia więc zasadniczych właściwości układu. Właściwości te są w istocie określone przez wartości własne, które są jednoznaczne w odróżnieniu od współrzędnych stanu i macierzy   układu.

Pojęciowo jest oczywiste, że właściwości układu nie mogą zależeć od wyboru zmiennych stanu. Istnieje nieskończenie wiele sposobów wyboru zmiennych stanu.

Powiązanie równań stanu z transmitancją

edytuj

Alternatywnym do równań stanu sposobem opisu układu dynamicznego (aczkolwiek zakładającym zerowy stan początkowy) jest transmitancja.

W przypadku układu o jednym wejściu i jednym wyjściu macierze stanu łączy z transmitancją następująca zależność:

 
 
Typowy model w przestrzeni stanów dla układu opisanego macierzami A, B, C, D

Istnieje nieskończenie wiele kombinacji macierzy A, B, C reprezentujących daną transmitancję G. Dla układu rzędu n = 2 macierze stanu można opisać następująco:

 
 
 

gdzie   jest dowolną liczbą rzeczywistą różną od 0.

W przypadku układu o wielu wejściach i wielu wyjściach w miejsce wspomnianej wyżej transmitancji określa się tak zwaną macierz transmitancji   Wówczas wyżej podana zależność przybiera postać:

 

(wytłuszczenia symboli we wzorach wskazują na to, że chodzi o wektory, a nie o wartości skalarne).

Postacie kanoniczne

edytuj

Jeśli wybiera się zmienne stanu dla układu o danej transmitancji, to spośród różnych możliwych sposobów wyboru tych zmiennych kilka jest szczególnie ciekawych.

Ze względów rachunkowych najbardziej pożądany jest sposób, w którym otrzymuje się macierz   układu w postaci diagonalnej, przy czym na przekątnej tej macierzy występują poszczególne wartości własne układu (czyli bieguny transmitancji). Model taki (ang. uncoupled form) można uzyskać, rozkładając transmitancję na ułamki proste. Jako zmienne wybiera się wtedy zmienne poszczególnych członów elementarnych połączonych równolegle. Zasadnicze wady tego opisu wynikają z braku ogólności metody (innymi słowy z braku jej kanoniczności), ponieważ przy wielokrotnych wartościach własnych układu nie jest możliwy rozkład na ułamki o mianownikach tylko pierwszego stopnia.

Cechy ogólności (kanoniczności) ma sposób wyboru zmiennych jako tzw. współrzędnych fazowych (ang. direct phase variable form) – w efekcie uzyskuje się model, który zawsze będzie odpowiadał modelowi sterowalnemu. Podobną strukturę ma model, w którym przyjmiemy   zmienne stanu ponumerujemy odwrotnie i uzależnimy wszystkie równania stanu od   (ang. direct feed-forward form) – uzyskuje się wówczas strukturę systematyczną i ogólną, a model taki odpowiada zawsze układowi obserwowalnemu. Modele sterowalne i obserwowalne są szczególnie pożądane z uwagi na to, że z niesterowalnością i nieobserwowalnością wiąże się pewna „niezręczność” w sposobie wprowadzenia sterowania lub wyprowadzenia wyjścia układu – w stosunku do układu swobodnego (zobacz też dekompozycja Kalmana).

Przykład

edytuj

Rozważmy przykład dla układu 4-wymiarowego o jednym wejściu i jednym wyjściu. Każda transmitancja ściśle właściwa może być zamieniona do przestrzeni stanów w następujący sposób. Transmitancję należy przekształcić tak, aby w mianowniku i liczniku pojawiły się odpowiednie współczynniki:

 

Sterowalna

edytuj
 
 

Powstały model jest sterowalny.

Obserwowalna

edytuj
 
 

Powstały model jest obserwowalny.

Równania dla układów ze sprzężeniem zwrotnym

edytuj

Sprzężenie zwrotne od wyjścia

edytuj
 
Typowy model w przestrzeni stanu dla układu ze sprzężeniem zwrotnym od wyjścia

W przypadku sprzężenia zwrotnego od wyjścia, wyjście mnożone jest przez macierz   i podawane jest na wejście systemu:

 

Ponieważ wartości   nie są ograniczone, można je w prosty sposób zanegować dla ujemnego sprzężenia zwrotnego. Obecność znaku minus (zwykle spotykana w notacji) ma charakter jedynie notacyjny – jego nieobecność nie ma wpływu na wyniki końcowe obliczeń. Wówczas równania:

 
 

można zapisać:

 
 

Przekształcając równanie wyjścia dla   i dokonując podstawienia w równaniu stanu, otrzymujemy:

 
 

Zaleta powyższego zapisu polega na tym, że na wartości własne macierzy   można wpływać poprzez odpowiednie określenie   przez dekompozycję własną wyrażenia  

Zakłada się przy tym, ze układ zamknięty jest sterowalny lub że niestabilne wartości własne   można uczynić stabilnymi poprzez odpowiedni dobór  

Często spotyka się uproszczenie dla takiego układu polegające na usunięciu   i określeniu   jako macierz jednostkową, co redukuje równania do postaci:

 
 

To znów sprowadza niezbędną dekompozycję własną jedynie do wyrażenia  

Zobacz też

edytuj

Bibliografia

edytuj
  • Tadeusz Kaczorek Teoria sterowania i systemów, Wydawnictwo Naukowe PWN, Warszawa 1993, ISBN 83-01-10936-X.
  • Tadeusz Kaczorek, A. Dzieliński, W. Dąbrowski, R. Łopatka, Podstawy teorii sterowania, Wydawnictwa Naukowo-Techniczne, Warszawa 2005, ISBN 83-204-2967-6.
  • Andrzej Markowski, Jerzy Kostro, Andrzej Lewandowski, Automatyka w pytaniach i odpowiedziach, Warszawa: Wydawnictwa Naukowo-Techniczne, 1979, ISBN 83-204-0110-0, OCLC 830948960.