Powierzchnia ślizgowa, warunek osiągania, dowód Lapunowa, chattering jako konsekwencja sgn.
Zanim wejdziemy w matematykę, trzy obrazy które pokazują o co chodzi.
Wyobraź sobie skater'a startującego z krawędzi półrury (half-pipe). Najpierw spada w kierunku dna rury — to faza reaching. Po dotarciu na dno jedzie po nim ku środkowi — to faza sliding. Dno rury to powierzchnia ślizgowa: zaprojektowana ścieżka, po której chcemy, żeby układ poruszał się stabilnie.
Klasyczny termostat: gdy temperatura jest niższa od zadanej — grzeje pełną mocą. Gdy wyższa — wyłącza grzanie. Bang-bang, dwa stany. SMC robi identyczną rzecz, tylko w dowolnej wymiarowości: daje +1 lub −1 zależnie od tego, po której stronie powierzchni jesteśmy. Jak termostat — szybkie przełączanie wokół zadanej wartości daje średnio efekt ciągły.
Hamulec elektryczny w pociągu daje siłę proporcjonalną do prędkości — im wolniej, tym słabiej. Nigdy nie zatrzymasz się dokładnie. Hamulec szczękowy (mechaniczny) hamuje tą samą siłą niezależnie od prędkości — i zatrzymuje cię w skończonym czasie. PID jest jak elektryczny hamulec, SMC — jak mechaniczny.
Wszystkie trzy mówią to samo. SMC: definiujemy zaprojektowaną „dolinę" w przestrzeni stanu ( — analogia 1). Sterowanie zachowuje się jak przełącznik (analogia 2): pcha mocno w jedną stronę gdy , w drugą gdy . To daje stałą siłę napędzającą trajektorię do — stąd zbieżność w czasie skończonym (analogia 3).
Najważniejsza różnica między SMC a PID widoczna jest w jednym wykresie: jak zbliżają się do celu w czasie. Najprostszy układ: (jakaś zmienna , której chcemy zerową wartość). Dwa kontrolery:
Analogia hamulca:
Dolny wykres (log |s|) pokazuje to formalnie. PID: — w skali log to prosta opadająca w nieskończoność. SMC: — w skali log krzywa się rozbiega w dół (logarytm zera = −∞), bo s dochodzi dokładnie do 0.
Kompromis: SMC „kupuje" dokładność i predicable time za cenę nieciągłości w sterowaniu (chattering). PID ma gładkie sterowanie, ale nigdy nie domyka błędu.
Mamy układ drugiego rzędu — np. masa-sprężyna-tłumik z modułu 1. Stan to para . Wprowadźmy uchyb (jak daleko jesteśmy od celu) i jego pochodną .
Definiujemy powierzchnię ślizgową jako prostą formułę:
Co to znaczy geometrycznie: w przestrzeni równanie wyznacza linię prostą przechodzącą przez początek układu z nachyleniem . Czyli .
Co to znaczy dynamicznie: gdy trajektoria jest na tej linii, dynamika błędu redukuje się do równania pierwszego rzędu . Rozwiązaniem jest — eksponencjalny zanik z czasem charakterystycznym . Wymiar dynamiki został zredukowany z 2 do 1.
Tylko jeden parametr: . Większe — szybsza zbieżność na manifoldzie, ale wymagająca większej początkowej siły. Praktycznie dla planta z dynamiką w skali sekund.
Powierzchnia ślizgowa to manifold w przestrzeni stanu, na którym chcemy znaleźć trajektorię. Definiujemy ją jako kombinację błędu i jego pochodnej:
Jeśli pętla utrzyma , to , a stąd — eksponencjalna zbieżność do zera z czasem charakterystycznym . To dynamika pierwszego rzędu uzyskana z planta drugiego rzędu — wymiar został zredukowany.
Sterowanie ślizgowe to dwa zadania na raz: (1) doprowadzić trajektorię z dowolnego stanu początkowego do w czasie skończonym; (2) utrzymać mimo zaburzeń. Pierwsze zadanie nazywa się reaching phase, drugie sliding phase.
Spróbuj: λ = 0.5 daje wolną zbieżność (τ = 2 s). λ = 5 daje τ = 0.2 s — dziesięciokrotnie szybciej. Cena: większe i większy peak w fazie reaching.
Chcemy, żeby malało z gwarantowaną prędkością. Nie wolniej niż jakaś stała :
Co to znaczy w trzech zdaniach: (1) Gdy , chcemy (s maleje). (2) Gdy , chcemy (s rośnie ku zeru). (3) W obu przypadkach maleje co najmniej o na sekundę. Więc czas dotarcia do zera jest ograniczony:
To jest cała magia SMC w jednej linijce. Reaching condition + nierówność = górne ograniczenie czasu dotarcia, niezależne od kształtu zaburzenia.
Mając plant z ograniczonym zaburzeniem , projektujemy:
Pierwsza część () — nominalna kompensacja: gdyby zaburzenia nie było, ta wartość trzymałaby trajektorię na manifoldzie. Druga część () — pchacz: stała siła w stronę manifoldu, jak magnes z analogii.
Definiujemy „energię" błędu: . Liczymy jak się zmienia w czasie:
Stosujemy nierówność :
Dla , prawa strona jest ujemna gdy . Czyli energia maleje z gwarantowaną szybkością. Ze wstępu o lemacie porównawczym dochodzimy do:
Wniosek: kontroler dotrze do manifoldu w najdłuższym przypadku w sekund. Nie 7 sekund, nie 10 — dokładnie ta liczba lub mniej. To jest jakościowo różne od PID, gdzie zbieżność jest tylko w granicy .
Predykcja teoretyczna: . Dla mamy , więc liniowo maleje z prędkością 2.5. Czas osiągania: s (4 kroki).
| k | t | x | s = −x | sgn(s) | u | ẋ = u + d |
|---|---|---|---|---|---|---|
| 0 | 0.0 | −1.0000 | 1.0000 | +1 | +2.0 | +2.5 |
| 1 | 0.1 | −0.7500 | 0.7500 | +1 | +2.0 | +2.5 |
| 2 | 0.2 | −0.5000 | 0.5000 | +1 | +2.0 | +2.5 |
| 3 | 0.3 | −0.2500 | 0.2500 | +1 | +2.0 | +2.5 |
| 4 | 0.4 | 0.0000 | 0.0000 | ← | reaching: s = 0 | |
| 5 | 0.5 | +0.2500 | −0.2500 | −1 | −2.0 | −1.5 |
| 6 | 0.6 | +0.1000 | −0.1000 | −1 | −2.0 | −1.5 |
Od zaczyna się chattering — przy idealnej dyskretyzacji trajektoria zostałaby na . W praktyce powoduje przelot przez 0 i oscylacje wokół manifoldu. Ekwiwalentne sterowanie utrzymujące to — moduł 4 (super-twisting, boundary layer) pokazuje jak to wygładzić.
Gdy układ jest dokładnie na manifoldzie (), to również — bo ma się tam utrzymać, nie tylko być. Z warunku wyznaczamy, jakie utrzymuje to równanie:
To jest źródło pierwszego członu sterowania z §4 — nie ad-hoc, tylko wprost z warunku pozostawania na manifoldzie.
W formule jest nieciągłość: nie jest gładką funkcją w punkcie . Z punktu widzenia teorii równań różniczkowych pojawia się problem: co jest w ogóle rozwiązaniem takiego równania?
Filippov (1960) podał definicję: w punkcie nieciągłości bierzemy powłokę wypukłą wszystkich pochodnych z otoczenia. Trajektoria musi spełniać , gdzie to ten zbiór. Praktycznie: gdy trajektoria jest na , wybieramy z dopuszczalnych pochodnych tę, która utrzymuje . Tę pochodną daje właśnie sterowanie .
W praktyce inżynierskiej można nie wnikać. Ważne jest to: średnia filippovska wartość w fazie sliding to dokładnie — czyli kompensuje zaburzenie. Stąd niewrażliwość na matched zaburzenia.
Plant z zaburzeniem:
Mówimy, że zaburzenie jest matched, jeśli istnieje skalarne :
Co znaczy: zaburzenie wchodzi do układu tym samym kanałem co sterowanie. Można je traktować jako efektywne zaburzenie sterowania .
Dla planta MSD ze wstrzyknięciem matched:
Pod warunkiem , suma pod sgn ma znak — czyli zawsze pcha do manifoldu. Po wejściu w sliding, średnia filippovska wartość przełącznika sgn jest dokładnie , więc . Człon przełączający kompensuje zaburzenie aktywnie.
Plant ze zaburzeniem unmatched :
wchodzi w równaniu pierwszej pochodnej, gdzie sterowania nie ma. SMC tego nie skompensuje — sgn(s) trafia tylko do , nie do . Trzeba rozszerzyć (integral SMC, hierarchical SMC) albo użyć innej techniki (H∞).
Phase portrait to wykres trajektorii w 2D przestrzeni . Najsilniejsza wizualizacja SMC, bo pokazuje na raz: powierzchnię (zielona linia), wektor pole (strzałki — czerwone gdy , niebieskie gdy ), trajektorię (niebieska linia). Faza reaching to ruch z dowolnego punktu w stronę zielonej linii. Faza sliding to ruch po niej do origin.
Co tu widać.Trajektoria startuje z punktu (e₀, ė₀) (pomarańczowa kropka), płynie w kierunku przeciwnym do wektorów pola (w prawej półpłaszczyźnie s > 0 — czerwone strzałki — w lewej s < 0 — niebieskie). Po dotarciu do prostej (zielona) ślizga się po niej do origin. Im większe λ, tym bardziej stroma prosta — szybsza zbieżność po osiągnięciu manifoldu, ale dłuższa faza reaching.
Kandydat Lapunowa jest niedodatnio monotoniczny pod warunkiem , gdzie . Czas dotarcia do ograniczony jest przez — skończony, nie asymptotyczny. To jest jakościowa różnica względem PID, gdzie zbieżność jest tylko asymptotyczna.
Spróbuj: zwiększ η do 25 — faza reaching skraca się dramatycznie. Zmniejsz η do 2 — przy zaburzeniu D = 1 (które plant nominalnie ma) trajektoria może w ogóle nie dotrzeć do manifoldu, lub dryfować poza nim.
Reguła ogólna: dla planta n-D z m wejściami, projektujemy m powierzchni ślizgowych. Manifold jest -wymiarowy. Na nim dynamika redukuje się z n do .
Stan: . Wejście: siła na wózek. Wybór powierzchni:
Wagi projektuje się tak, żeby na dynamika redukowana (3D) była stabilna. Standardowo: rozwiązuje się redukowane DARE.
6 osobnych powierzchni , po jednej dla każdego przegubu. Mogą być sprzężone (w niektórych konstrukcjach zawiera pochodne sąsiednich przegubów). Klasyczne podejście: dekoupling przez transformację do postaci kanonicznej.
Pełny formalizm: Edwards & Spurgeon rozdz. 4–5. Praktyka robotyczna: Slotine & Li „Applied Nonlinear Control" rozdz. 7.
Moduł 4 pokazuje 3 sposoby na to: boundary layer (saturacja zamiast sgn), super-twisting (drugiego rzędu SMC z gładkim u), terminal SMC.
Chattering to nieuchronna konsekwencja : w idealnej teorii ciągłej, gdy ślizga się po manifoldzie, sterowanie przełącza się nieskończenie szybko z na , żeby utrzymać . W rzeczywistości:
Konsekwencja: oscyluje w pasie wokół zera o szerokości proporcjonalnej do . Sterowanie przełącza się prostokątnie z amplitudą — z punktu widzenia silnika to katastrofalna jakość mocowania, dla zaworu — szybkie zużycie. Zoom-in pokazuje wprost: dla ms widać wyraźne oscylacje pozycji rzędu kilku procent referencji.
Skala chatteringu jest liniowa w i w . Stąd dwa typowe ataki na problem: szybsze próbkowanie (drogie sprzętowo), albo zmiana algorytmu (boundary layer, super-twisting — następny moduł).
Wybierz SMC, gdy spełnione wszystkie warunki:
Wybierz coś innego, gdy: