TS Spraw, Mechatronika WAT, Semest IV, Teoria sterowania, Laboratorium, Sprawka, Obserwator stanu


Wojskowa Akademia Techniczna

Wydział Mechatroniki

INSTYTUT SYSTEMÓW MECHATRONICZNYCH

ĆWICZENIE LABORATORYJNE nr.2

Temat: METODY OPISU UKŁADÓW LINIOWYCH, STACJONARNYCH.

Badanie sterowalności i obserwowalności liniowego układu dynamicznego.”

Prowadzący: dr.inż. Krzysztof Motyl

Ocena:

Grupa

Podgrupa

Data

A6I1S1

1

22.04.08r.

Wykonawcy:

Układy liniowe , stacjonarne możemy opisać za pomocą:

Modelowanie układu dynamicznego:

0x01 graphic

W trakcie ćwiczeń laboratoryjnych badaliśmy powyższy liniowy układ, którego wymuszeniem, a zarazem wejściem jest siła F(t) , a wielkością wyjściową jest przemieszczenia masy y(t). Jest to układ o jednym wejściu i jednym wyjściu i możemy opisać nasz schemat za pomocą następującego równania: *

0x01 graphic

Układ ten rozwiązaliśmy w programie MATLAB-SIMULINK za pomocą kilku sposobów:

Sposób 1 Z wykorzystaniem równań różniczkowych

Nasze równanie* opisujące układ jest drugiego stopnia, ma więc 2 integratory. Podstawiając odpowiednio za 0x01 graphic
=>0x01 graphic
itd. Otrzymujemy następujące przekształcenia:

0x01 graphic

Równanie to przedstawiliśmy w postaci schematu blokowego następująco :

0x01 graphic

Po podstawieniu odpowiednich wartości parametrów:

0x01 graphic

Otrzymaliśmy następujące wykresy na naszym oscyloskopie:

0x01 graphic

Wnioski:

  1. Rysunek pierwszy od góry przedstawia nam wykres sygnału przed pierwszym całkowaniem. Można zauważyć bardzo krótki czas narastania oraz fakt iż sygnał ustabilizował się na tym samym poziomie co przed wymuszeniem.

  2. Drugi rysunek przedstawia zachowanie sygnału po 2 całkowaniu. Widać na nim, że czas narastania jest większy niż w pierwszym oraz, że jest także stabilny asymptotycznie.

  3. Ostatni wykres jest to sygnał po 3 całkowaniach. Możemy zauważyć, że nie jest już on stabilny asymptotycznie, a czas narastanie jest największy.

Sposób 2 Z wykorzystaniem funkcji przejścia

Nasze ogólne równanie * po przeprowadzeniu transformaty Laplace`a wygląda następująco:

0x01 graphic

Co mogliśmy przedstawić jako schemat blokowy za pomocą bloku Transfer Fcn w następujący sposób:

0x01 graphic

Podstawiając te same wartości parametrów otrzymaliśmy następujący wykres:

0x01 graphic

Wnioski:

Na podstawie przedstawionego wykresu bardzo łatwo możemy zauważyć iż jest to ten sam sygnał co w poprzednim przykładzie na ostatnim wykresie. Przykład ten pokazuje nam jak możemy uprościć sobie schemat budowy układu za pomocą transmitancji oraz że rozwiązanie to jest równoważne z wynikiem poprzedniego przykładu.

Sposób 3 Z wykorzystaniem przestrzeni stanów

Nasze ogólne równanie * możemy przekształcić do postaci macierzowych, które będą spełniać następujące równania przestrzeni stanu:

0x01 graphic

, gdzie :

- A- macierz stanu,

- B- macierz wejścia,

- C- macierz wyjścia,

- D- bezpośrednia macierz transmisji

Uzyskane przez nas przekształcenia wynoszą:

0x01 graphic

Schemat naszego układu mogliśmy zbudować używając bloku State-Space:

0x01 graphic

Następnie określając wartości naszych macierzy i parametrów:

0x01 graphic

Uzyskaliśmy następujące wykresy na oscyloskopie:

0x01 graphic

Wnioski:

Na podstawie wykresów możemy stwierdzić, że pierwszy w porównaniu do drugiego nie jest stabilny asymptotycznie, natomiast pierwszy osiąga dwukrotnie większą pierwszą amplitudę odchylenia od wartości ustalonej oraz ma dłuższy czas narastania. Czas regulacji jest podobny.

4. Obliczanie parametrów z punktu widzenia wej-wyj:

0x01 graphic

a20x01 graphic
(t) + a10x01 graphic
(t) + a00x01 graphic
(t) = b0y(t)

5.Analiza modelu matematycznego przy różnych wartościach współczynnika tłumienia:

Nasz schemat blokowy zbudowaliśmy w taki sposób, aby uzyskać jednocześnie wszystkie wykresy sygnałów o różnym tłumieniu.

0x01 graphic

Następnie uzupełniliśmy wartości naszych parametrów:

0x01 graphic

W wyniku czego uzyskaliśmy następujące wykresy sygnałów:

0x01 graphic

Wnioski:

Wykres żółty przedstawia nam sygnał bez tłumienia. Jest to człon niedotłumiony - oscylacyjny. Można zauważyć, że wykres nie dąży do stabilizacji.

Wykres różowy jest wykresem członu o tłumieniu ξ= 0,2. Jest to również człon oscylacyjny, jednak widać, że po pewnym czasie regulacji stabilizuje się on,

Wykres niebieski charakteryzuje się tłumieniem ξ=0.7 . Jest on także oscylacyjny, lecz widać że oscylacje te są bardzo małe.

Wykres czerwony jest już członem inercyjnym o dwóch jednakowych stałych czasowych o tłumieniu krytycznym równym ξ=1.

Wykres zielony jest członem inercyjnym o dwóch różnych stałych czasowych i tłumieniu nadkrytycznym ξ=2.

6.Sprawdzenie obserwowalności i sterowalności rozpatrywanego układu dynamicznego.

Po podstawieniu wartości parametrów:

m=2, k=200, c=5.

Nasze macierze wyglądają następująco:

0x01 graphic

A=0x01 graphic
, B=0x01 graphic
0x01 graphic
, C=0x01 graphic
0x01 graphic

  1. Aby układ dynamiczny był sterowalny: det[B,AB]≠0

det0x01 graphic
=-0,25 (≠ 0) => układ sterowalny,

  1. Aby układ był obserwowalny: det0x01 graphic
    ≠0

det0x01 graphic
=1 ( ≠ 0) => układ obserwowalny

Badany przez nas układ jest zarówno obserwowalny jak i sterowalny



Wyszukiwarka