Wiskunde · Begrippenlijst

Wat is stabiliteit (Lyapunov)?

Definitie 19.10 Universitaire wiskunde — Bachelor jaar 3 · Hoofdstuk 19 — Gewone differentiaalvergelijkingen

Een evenwicht is een punt xˉ\bar x met F(xˉ)=0F(\bar x) = 0 (zodat φt(xˉ)=xˉ\varphi_t(\bar x) = \bar x). Het heet stabiel als er voor elke ε>0\varepsilon > 0 een δ>0\delta > 0 is zodat x0xˉ<δ\norm{x_0 - \bar x} < \delta impliceert dat de oplossing voor alle t0t \geq 0 bestaat met φt(x0)xˉ<ε\norm{\varphi_t(x_0) - \bar x} < \varepsilon; en asymptotisch stabiel als bovendien φt(x0)xˉ\varphi_t(x_0) \to \bar x voor alle x0x_0 nabij xˉ\bar x.

Voorbeelden

Voorbeeld 19.9 (Het vlak, geklasseerd)

Voor x=Axx' = Ax met AM2(R)A \in M_2(\R) inverteerbaar wordt het faseportret nabij 00 beslist door τ=trA\tau = \operatorname{tr}A en δ=detA\delta = \det A, via de eigenwaarden λ±=τ±τ24δ2\lambda_\pm = \frac{\tau \pm \sqrt{\tau^2 - 4\delta}}2:

  • δ<0\delta < 0: reële eigenwaarden met tegengesteld teken — een zadel; twee banen komen binnen, twee gaan weg, alle andere vliegen voorbij. Altijd instabiel.
  • δ>0\delta > 0, τ24δ\tau^2 \geq 4\delta: reële eigenwaarden met hetzelfde teken (=signτ= \operatorname{sign}\tau) — een knoop, stabiel dan en slechts dan als τ<0\tau < 0; de banen raken de trage eigenrichting.
  • δ>0\delta > 0, τ2<4δ\tau^2 < 4\delta, τ0\tau \neq 0: toegevoegd complexe eigenwaarden τ2±iω\frac\tau2 \pm \iu\omega — een spiraal (focus), stabiel dan en slechts dan als τ<0\tau < 0; de oplossingen zijn eτt/2\eu^{\tau t/2} maal rotaties met periode 2πω\frac{2\pi}\omega.
  • τ=0\tau = 0, δ>0\delta > 0: zuiver imaginaire eigenwaarden — een centrum: gesloten banen (ellipsen), stabiliteit zonder asymptotische stabiliteit, precies het grensgeval dat Stelling 19.12 voor niet-lineaire stelsels niet kan beslissen (het onderste evenwicht van de slinger, Probleem 19.1, zit hier).

De grensparabool τ2=4δ\tau^2 = 4\delta draagt de ontaarde knopen (jordanblokken: banen met één enkele raakrichting). Alles wordt van twee getallen afgelezen — daarom is de eerste reflex vóór elk vlak faseportret tr\operatorname{tr} en det\det te berekenen; zo geeft A=(011c)A = \bigl(\begin{smallmatrix}0 & 1\\ -1 & -c\end{smallmatrix}\bigr) (de gedempte oscillator): δ=1>0\delta = 1 > 0, τ=c\tau = -c: een stabiele spiraal voor 0<c<20 < c < 2, een stabiele knoop voor c2c \geq 2 — onderdemping tegenover overdemping, in één oogopslag.

Lees in het hoofdstuk →