1
Kryterium stabilności
2
Stabilność liniowych układów sterowania
Układ zamknięty liniowy i stacjonarny opisany równaniem (1) jest stabilny, jeżeli dla
skończonej wartości zakłócenia przy dowolnych wartościach początkowych jego odpowiedź
ustalona przyjmuje skończone wartości.
)
(
)
(
)
(
)
(
)
(
)
(
1
1
1
1
1
1
t
u
b
dt
t
u
d
b
dt
t
u
d
b
t
x
a
dt
t
x
d
a
dt
t
x
d
a
o
m
m
m
m
m
m
o
n
n
n
n
n
n
+
+
+
=
+
+
+
−
−
−
−
−
−
⋯
⋯
(1)
Transmitancja operatorowa tego układu ma postać:
)
(
)
(
)
(
0
1
1
1
0
1
1
1
s
M
s
L
a
s
a
s
a
s
a
b
s
b
s
b
s
b
s
G
n
n
n
n
m
m
m
m
=
+
+
+
+
+
+
+
+
=
−
−
−
−
⋯
⋯
(2)
Stąd jego równanie charakterystyczne
0
)
(
0
1
1
1
=
+
+
+
+
=
−
−
a
s
a
s
a
s
a
s
M
n
n
n
n
⋯
(3)
W ujęciu matematycznym warunkiem koniecznym i dostatecznym na to, ażeby układ
zamknięty był stabilny jest, aby pierwiastki równania charakterystycznego (3) miały ujemne
części rzeczywiste. Rozwiązanie tego równania wystarczy, więc dla stwierdzenia czy dany
układ liniowy jest stabilny. Jednak w praktyce ta metoda nie zawsze jest dogodna i
wystarczająca.
Z tego względu zostały opracowane metody pozwalające na badanie stabilności bez
rozwiązywania równania charakterystycznego są to tzw. kryteria stabilności. Kryteria te
dzielą się na: algebraiczne do, których należą kryteria Routha i Hurwitza oraz
częstotliwościowe Michajłowa i Nyquista.
Jeżeli pierwiastki równania charakterystycznego mogą być wyrażone jako
I ich części rzeczywiste są ujemne (leżą na lewej półpłaszczyźnie) to układ jest stabilny.
Metoda wymaga obliczenie pierwiastków r.char , stąd pojawiły się jej odpowiedniki
(wymagające mniej obliczeń arytmetycznych)
Warunki stabilności:
Re(s
i
)<0 dla i=1,2,3,4........ – stabilny asymptotycznie
Re(s
i
)=0 dla dowolnego i, pozostałe Re(s
i
)<0 - granica stabilności
Re(s
i
)>0 dla dowolnego i – układ niestabilny
Dla pierwiastków rzeczywistych
s=σ
Dla pierwiastków zespolonych:
s=σ±jω
Przykłady funkcji:
2.1
1.1
Kryterium Hurwitza
Warunkiem koniecznym i dostatecznym stabilności układu liniowego i stacjonarnego jest,
ażeby wszystkie współczynniki wielomianu charakterystycznego transmitancji tego układu
istniały i były dodatnie, a ponadto ażeby wyznacznik
n
∆
(1) zwany wyznacznikiem Hurwitza
oraz jego podwyznaczniki
1
3
2
,
,
,
−
∆
∆
∆
n
⋯
były dodatnie.
Jeżeli którykolwiek współczynnik jest ujemny lub równy zeru albo którykolwiek
podwyznacznik jest ujemny to układ jest niestabilny.
Jeśli dowolny z podwyznaczników jest równy zeru to oznacza, że równanie charakterystyczne
układu ma między innymi pierwiastki urojone i wtedy układ jest na granicy stabilności. Na
jego wyjściu występują drgania o ustalonej amplitudzie.
0
2
4
1
3
0
2
3
1
4
2
5
3
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
a
a
a
a
a
a
a
a
a
a
a
a
a
a
a
n
n
n
n
n
n
n
n
n
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
−
−
−
−
−
−
−
=
∆
(4)
1
1
−
=
∆
n
a
2
3
1
2
−
−
−
=
∆
n
n
n
n
a
a
a
a
3
1
4
2
5
3
1
3
0
−
−
−
−
−
−
−
=
∆
n
n
n
n
n
n
n
n
a
a
a
a
a
a
a
a
⋮
1
3
5
0
2
4
1
3
3
1
4
2
5
3
1
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
a
a
a
a
a
a
a
a
a
a
a
a
a
a
a
a
n
n
n
n
n
n
n
n
n
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
−
−
−
−
−
−
−
−
=
∆
Prz ykł ad 1
Za pomocą kryterium Hurwitza zbadać stabilność układu zamkniętego, którego równanie
charakterystyczne ma postać:
0
2
11
7
4
6
)
(
2
3
4
5
=
+
+
+
+
+
=
s
s
s
s
s
s
M
Zauważmy, że spełniony jest warunek konieczny stabilności, ponieważ wszystkie
współczynniki równania charakterystycznego są dodatnie.
Wyznacznik Hurwitza utworzony ze współczynników wielomianu tego równania ma postać:
2
7
6
0
0
0
11
4
1
0
0
2
7
6
0
0
0
11
4
1
0
0
2
7
6
5
=
∆
Obliczamy wartość wyznacznika za pomocą polecenia det, które wprowadzamy w oknie
poleceń MATLAB-a według poniższej składni. Argumentem polecenia det zapisanym w
nawiasach okrągłych jest macierz współczynników wyznacznika Hurwitza, która z kolei
zapisana jest w nawiasach kwadratowych. Poszczególne elementy wierszy tej macierzy
oddzielone są odstępami, natomiast wiersze – oddzielone są średnikami.
» delta_5=det([6 7 2 0 0;1 4 11 0 0;0 6 7 2 0;0 1 4 11 0;0 0 6 7 2]),
delta_5 =
-5846
Ujemna wartość wyznacznika Hurwitza wskazuje na to, że badany układ jest niestabilny.
Kryterium Routha
Pierwszym krokiem w uproszczeniu kryterium Hurwitza, nazywanym kryterium Routha, jest
umieszczenie współczynników równania
w dwóch wierszach. Pierwszy wiersz składa się z nieparzystych współczynników, natomiast
drugi wiersz z parzystych współczynników licząc od najwyższej potęgi wielomianu
charakterystycznego. Dla równania tego pierwsze dwa wiersze tablicy są
następujące:
Następnym krokiem jest wypełnienie następnych wierszy tablicy Routha w następujący
sposób:
i tak dalej. Kolumna z lewej strony tablicy Routha jest kolumną odniesienia i służy do
identyfikacji obliczeń. Ostatni wiersz tablicy Routha ma zawsze w tej kolumnie element s
0
.
Po skompletowaniu tablicy Routha ostatnim krokiem jest określenie znaków spółczynników
pierwszej kolumny tablicy, która zawiera informacje o pierwiastkach równania. Przyjęte
zostało następujące założenie:
Wszystkie pierwiastki równania charakterystycznego znajdują się w lewej półpłaszczyźnie
jeśli wszystkie elementy pierwszej kolumny tablicy Routha mają ten sam znak. Liczba zmian
znaków w elementach pierwszej kolumny równa jest liczbie pierwiastków w prawej
półpłaszczyźnie.
Przykład
Rozważmy równanie
w którym nie brakuje elementów i wszystkie współczynniki są tego samego znaku. Spełniony
jest warunek konieczny dotyczący współczynników, jednak warunek dostateczny musi zostać
jeszcze sprawdzony. Pierwszą czynnością jest zainicjowanie tablicy, w kolumnie z lewej
strony znajdują się potęgi s, natomiast współczynniki wielomianu rozdziela się pomiędzy
pierwszy drugi wiersz w sposób pokazany poniżej. W pierwszym wierszu znajdują się w
kolejności współczynniki nieparzyste, natomiast w drugim parzyste
Tablica kompletowana jest poczynając od góry wiersz po wierszu, obliczając elementy
następnego wiersza. Każdy obliczany element wyprowadzany jest na podstawie czterech
elementów znajdujących się w dwóch wyższych wierszach, dwa z nich są w lewej kolumnie
i dwa w kolumnie znajdującej się na prawo od obliczanego elementu. W każdym przypadku,
obliczany element ma ujemny wyznacznik z czterech znajdujących się wyżej elementów,
podzielony jest przez lewy dolny element wyznacznika. Dla przykładu, pierwszy element
wiersza s
2
drugi element wiersza s
2
pierwszy element wiersza s
1
i tak dalej.
Skompletowana tablica Routha pokazana jest poniżej. Liczba pierwiastków wielomianu M(s)
znajdująca się w prawej półpłaszczyźnie jest równa liczbie zmian znaków lewej kolumny
tablicy, przesuwając się z góry na dół.
W tym przykładzie są dwie zmiany znaków w lewej kolumnie co oznacza, że wielomian M(s)
ma dwa pierwiastki w prawej półpłaszczyźnie.
Tablica Routha – przypadki szczególne
Zero w pierwszej kolumnie tablicy Routha
W pierwszym przypadku, jeśli zero pojawia się w pierwszym elemencie wiersza, wówczas
wszystkie elementy w następnym wierszu mają wartości równe nieskończoności i dalsze
wypełnianie tablicy nie jest możliwe. Aby poradzić sobie z tą sytuacją zastępuje się pierwszy
element w pierwszej kolumnie przez bardzo mała liczbę dodatnią ε i kontynuuje się
obliczanie pozostałych elementów. Przypadek ten zostanie zilustrowany przez następujący
przykład.
Rozważmy następujące równanie charakterystyczne układu liniowego
Nie wszystkie współczynniki mają ten sam znak, czyli na pewno występują pierwiastki
w prawej półpłaszczyźnie. Sprawdźmy przy użyciu kryterium Routha ile pierwiastków
znajduje się w prawej półpłaszczyźnie. Przy kompletowaniu tablicy Routha dwa pierwsze
wiersze
uzyskuje się bezpośrednio ze współczynników wielomianu. Brakujące współczynniki
uzupełnia się zerami.
Współczynniki drugiego wiersza można podzielić przez 3, co pozwoli na uproszczenie
obliczeń.
Tablica ta nie może być dalej kompletowana w zwykły sposób, ponieważ nie można dzielić
przez zero. W pierwszej kolumnie pojawiło się zero, przy czym nie cały wiersz jest zerowy.
Sytuacja z zerem w pierwszej kolumnie rozwiązywana jest w ten sposób, że zamiast zera
wprowadza się bardzo małą liczbę dodatnią ε. Dla powyższego wielomianu zastępując zero
w pierwszej kolumnie przez ε i po wyznaczeniu kolejnych elementów tablicy w zależności od
ε
otrzymuje się
Dla wszystkich wyrażeń w pierwszej kolumnie zawierających ε wyznacza się granicę ε → 0,
przy założeniu dodatniej wartości ε, na przykład dla
Uzyskane znaki elementów pierwszej kolumny tablicy Routha
W tym przypadku są trzy zmiany znaku w pierwszej kolumnie, więc badany wielomian ma
trzy pierwiastki w prawej półpłaszczyźnie.
Zerowy wiersz w tablicy Routha
- świadczy o jednym z 3 przypadków
1.Równanie ma przynajmniej jedną parę pierwiastków o przeciwnych znakach
2. Równanie ma jedną lub więcej par pierwiastków sprzężonych na osi urojonych
3 Równanie ma pary pierwiastków tworzących symetrie wokół początku układu
W sytuacji gdy pojawia się cały wiersz zerowy w tablicy Routha, tworzy się równanie
pomocnicze p(s) = 0, które formuje się ze współczynników wiersza znajdującego się powyżej
wiersza zerowego w tablicy Routha. Rozwiązując równanie pomocnicze otrzymuje się
również pierwiastki równania oryginalnego. Aby dalej wypełniać tablicę Routha wykonuje się
następujące kroki:
1. Tworzy się równanie pomocnicze p(s) = 0 przez użycie współczynników z wiersza
znajdującego się powyżej wiersza zerowego.
2. Wyznacza się pochodną równania pomocniczego względem s;
3. Zastępuje się wiersz zerowy współczynnikami wielomianu
4. Kontynuuje się wypełnianie tablicy Routha z użyciem nowo utworzonego wiersza
współczynnikami zastępującymi wiersz zerowy.
5.
Interpretuje się zmianę znaków współczynników w pierwszej kolumnie tablicy Routha
w zwykły sposób.
Przykład
Rozważmy następujące równanie charakterystyczne układu liniowego
z badania współczynników widać, że w wielomianie występują pierwiastki z prawej
półpłaszczyzny.Tablica Routha zaczyna się następująco:
współczynniki wiersza s
4
można podzielić przez 4:
Pojawił się wiersz zerowy. Wprowadzamy równanie pomocnicze ze współczynników
znajdujących się nad wierszem zerowym w wierszu s
4
Różniczkując wielomian p(s) względem s otrzymuje się
Otrzymanymi współczynnikami 4 oraz 2 zastępujemy wiersz zerowy.
Współczynniki równania dp(s)/ds.
Współczynniki uzyskanego wiersza można podzielić przez 2. Pozostała część tablicy Routh
jest następująca
Wielomian pomocniczy jest czwartego rzędu, czyli w równaniu M(s) występują dwie pary
pierwiastków. Z badania pierwszej kolumny tablicy Routha uzyskanej dla równania M(s)
widać,że występuje jedna zmiana znaku, czyli jedna para pierwiastków jest o przeciwnych
znakach. Dwa pozostałe pierwiastki muszą znajdować się na osi urojonej . Z analizy tablicy
Routha dla tego przykładowego wielomianu wyznaczyliśmy następujące typy pierwiastków:
Lewa półpłaszczyzna LP = 3
Prawa półpłaszczyzna PP = 1
Oś urojona IA = 2
Układ ze strojonym parametrem
równanie charakterystyczne:
Tablica Routha dla równania
Na podstawie powyższej tablicy uzyskuje się dwa warunki stabilności: z wiersza s
2
,
otrzymanym warunkiem stabilności jest 34 − K > 0, z wiersza s
1
warunek
(−K
2
+ 612)/(34−K) > 0, natomiast dla s
0
warunek K > 0.
Z rozważenia tych dwóch warunków otrzymany zakres stabilności dla parametru K
0 < K < 24.739
Wyszukiwarka
Podobne podstrony:
instrukcja 7 kryt Nyquista id 215675
wykład 6 instrukcje i informacje zwrotne
Instrumenty rynku kapitałowego VIII
05 Instrukcje warunkoweid 5533 ppt
Instrukcja Konwojowa
2 Instrumenty marketingu mix
Promocja jako instrument marketingowy 1
Promocja jako instrument marketingowy
Instrukcja do zad proj 13 Uklad sterowania schodow ruchom
Instrukca 6 2
więcej podobnych podstron