dgl1 08
dgl1 08
dgl1 08
8. Stabilität
A. Allgemeines.
Wir betrachten ein allgemeines DGL-System erster Ordnung:
y0(t) = f (t, y(t)) (8.1)
mit y(t) ∈ Rn, hinreichend glatter rechter Seite f : I × D → Rn,
I ⊂ R offenes Intervall und D Gebiet in Rn. y∗ sei eine spezielle
Lösung von (8.1). Wir fragen nach dem Verhalten benachbarter“
”
Lösungen y(t; t0, y0).
Definition (8.2)
a) y∗ heißt stabil (auch Ljapunov–stabil) auf I, falls y∗ auf I
definiert ist und es zu t0 ∈ I und ε > 0 stets ein δ > 0 gibt, mit
∀ y0 : ky0 − y∗(t0)k < δ ⇒ ky(t; t0, y0) − y∗(t)k < ε (∀ t ∈ I) .
138
b) Kann man in a) die Größe δ > 0 unabhängig von t0 wählen,
so heißt y∗ gleichmäßig stabil auf I.
c) Ist y∗ auf [ a, ∞[ ⊂ I erklärt, so heißt y∗ dort asymptotisch
stabil, falls y∗ stabil ist und es zu jedem t0 ≥ a ein δ(t0) > 0 gibt
mit
Bemerkungen (8.3)
a) In (8.2) a) wird implizit vorausgesetzt, dass y(·; t0, y0) für
ky0 − y∗(t0)k < δ (zumindest) auf I erklärt ist. Zumeist wird auch
bei der Definition der Stabilität ein Intervall der Form I = [a, ∞[
zugrunde gelegt. Für autonome DGL lässt sich darüber hinaus
139
(zeitliche Verschiebungsinvarianz!) a := t0 := 0 setzen.
Stabilität bedeutet dann, dass die Trajektorien mit Anfangswer-
ten in der Nähe von y∗(0) für alle Zeiten t ≥ 0 in einem ε–
Schlauch um y∗(t) verbleiben, während sie sich im Falle der
asymtotischen Stabilität y∗(t) für t → ∞ annähern.
b) Mit Hilfe der Transformation
z(t) := y(t) − y∗(t)
z0(t) = f (t, z(t) + y∗(t)) − f (t, y∗(t)) (8.4)
=: f ∗(t, z(t))
kann man sich o.B.d.A. auf die Untersuchung der Stabilität der
Nulllösung z∗ = 0 zurück ziehen. In der Tat ist z∗ := 0 ein
Gleichgewichtspunkt des DGL-Systems z0 = f ∗(t, z), d.h., es
gilt f ∗(t, z∗(t)) ≡ 0.
140
B. Stabilität bei linearen DGL.
Gegeben sei ein homogenes, lineares DGL-System
∀ t ≥ t0 ≥ a : Y(t) Y(t0)−1 ≤ M .
141
c) y∗ = 0 ist genau dann asymptotisch stabil, falls gilt:
lim kY(t)k = 0 .
t→∞
Beweis:
Nach Satz (6.8) lautet die allg. Lösung von (8.5) y(t) = Y(t) c,
c ∈ Rn. Entsprechend ist Lösung der AWA gegeben durch
y(t; t0, y0) = Y(t) Y(t0)−1 y0 .
zu a): Ist y∗ = 0 stabil, so ist jede Lösung von y0 = A y be-
schränkt. Damit sind auch die Spaltenvektoren der Fundamen-
talmatrix Y(t) beschränkt und damit auch kY(t)k selbst.
Umgekehrt folgt aus kY(t)k ≤ M (∀t), M > 0:
142
zu b): Ist y∗ = 0 gleichm. stabil, so folgt Y(t) Y(t0)−1 y0 ≤ ε
für alle t ≥ t0 ≥ a und ky0k ≤ δ. Damit ergibt sich
Y(t) Y(t0)−1 ≤ ε/δ. Die Umkehrung erhält man wie in a).
zu c): Dies folgt unmittelbar aus obiger Darstellung der Lösun-
gen.
Satz (8.7)
Sei λ(t) der größten EW der symmetrischen Matrix A(t) + A(t)T.
Ist dann
Z∞
λ(t) dt = −∞ , a ≤ t0 < ∞,
t0
so ist y∗ = 0 asymptotisch stabil.
143
Beweis:
d d
ky(t)k2 = y(t)T y(t)
dt dt
= (A(t) y(t))T y(t) + y(t)T (A(t) y(t))
= y(t)T A(t)T + A(t) y(t)
≤ λ(t) y(t)T y(t)
= λ(t) · ky(t)k2 .
t0
144
Satz (8.8) (Stabilitätssatz II)
Gegeben sei ein lineares DGL-System y0 = Ay, mit konstanter(!)
Matrix A ∈ R(n,n). λj bezeichne die EWe von A, j = 1, . . . , n,
und g(λj ) bzw. a(λj ) deren geom. bzw. algebr. Vielfachheiten.
Dann gelten:
a) y∗ = 0 ist strikt stabil ⇔ ∀ j : Re (λj ) < 0,
b) y∗ = 0 ist gleichmäßig stabil ⇔ ∀ j : Re (λj ) ≤ 0 und
Re (λj ) = 0 ⇒ g(λj ) = a(λj ) .
c) In allen anderen Fällen ist y∗ = 0 instabil.
Beweis:
Man erhält die obigen Aussagen unmittelbar aus der Darstellung
des Fundamentalsystems in Abschnitt 6 D, insbesondere aus den
Darstellungen (6.21), (6.22) und (6.28).
145
Beispiele (8.9)
a) Der Gleichgewichtspunkt y∗ = 0 des folgenden DGL-Systems
ist instabil
4 8 4
y0(t) = −1 −2
A y(t) = 1 y(t)
−2 −4 −6
Die EWe λ1 = −4 und λ2 = 0 der Koeffizientenmatrix erfüllen
zwar Re (λj ) ≤ 0, jedoch ist λ2 = 0 doppelter Eigenwert mit
der geometrischen Vielfachheit g(λ2) = 1 < a(λ2) = 2.
146
b) Das DGL-System
y10 (t)
! ! ! !
−1 1 y1(t) 5
= A y(t) + b = +
y20 (t) −1 −1 y2(t) 1
besitzt den Gleichgewichtspunkt y∗ = (3, −2)T.
Die Transformation z := y − y∗ liefert das homogene DGL-
System z0 = A z mit der obigen Matrix A.
Die Eigenwerte von A sind λ1,2 = −1 ± i. Damit ist der Gleich-
gewichtspunkt y∗ nach (8.8) a) strikt stabil.
Ohne Beweis zitieren wir noch ein Kriterium für strikte Stabi-
lität, welches mit dem charakteristischen Polynom der Matrix A
arbeitet und ohne die konkrete Berechnung der Eigenwerte aus-
kommt! Dieses Kriterium findet in der Regelungstheorie vielfach
Anwendung.
147
Satz (8.10) (Kriterium von Routh und Hurwitz∗)
Pn k
Gegeben sei ein reelles Polynom p(z) = k=0 ak z , an > 0 .
Es sind äquivalent:
↑ w2
w →
1
↑ w2
w →
1
153
c) Asymptotisch stabiler Knotenpunkt 2. Art:
λ1, λ2 < 0, λ1 6= λ2, wj (t) = wj0 eλj t
↑ w2
w →
1
154
d) Instabiler Knotenpunkt 2. Art:
λ1, λ2 > 0, λ1 6= λ2, wj (t) = wj0 eλj t
↑ w2
w →
1
155
e) Instabiler Sattelpunkt:
λ1 < 0 < λ 2
↑ w2
w →
1
156
f) Parallele Geraden:
λ1 = 0, λ2 6= 0;
Beachte: 0 ist kein isolierter stationärer Punkt!
↑ w2
w1 →
↑
w2
w1 →
↑ w2
w →
1
159
i) Instabiler Strudelpunkt:
λ1 = λ2 = λ = α + i β, β 6= 0, α>0
↑ w2
w →
1
160
j) Stabiler Wirbelpunkt:
λ1 = λ2 = λ = i β, β 6= 0
↑ w2
w1 →
161
D. Der nichtlineare Fall.
Wir betrachten wieder ein autonomes (nichtlineares) DGL-
System
y0(t) = f (y(t)) . (8.17)
y∗ = 0 ∈ Rn sei Gleichgewichtspunkt des Systems, also f (0) = 0.
Häufig gelingt es, die Stabilität der DGL (8.17) durch Lineari-
sierung der rechten Seite nachzuweisen. Ist f hinreichend glatt,
so liefert die Taylor–Entwicklung von f um y = 0:
y0(t) = A y(t) + g(y(t)), A := Jf (0),
(8.18)
g (0) = 0, g(y) = O(kyk2) .
Ist die verkürzte“ DGL y0 = A y nun strikt stabil, bzw.
”
instabil, so überträgt sich dies auch auf das nichtlineare DGL-
System.
162
Satz (8.19) (Stabilitätssatz III)
Unter den obigen Voraussetzungen lässt sich zeigen:
a) Gilt für alle EWe λj von A = Jf (0) : Re (λj ) < 0, so ist
y∗ = 0 ein strikt stabiler Gleichgewichtspunkt von y0 = f (y).
Warnung (8.20)
Die Stabilität bzw. Instabilität des linearisierten DGL-Systems
überträgt sich nicht notwendig auf das nichtlineare DGL-System,
falls für alle EWe Re (λj ) ≤ 0 gilt, es jedoch zumindest einen
Eigenwert λk gibt mit Re (λk ) = 0.
163
Beispiel (8.21)
Wir betrachten die DGL des physikalischen Pendels
g
ϕ00(t) = − sin ϕ(t) = − ω 2 sin ϕ(t).
`
Mit y1 := ϕ, y2 := ϕ0 ergibt sich das DGL-System:
y10 (t) = y2(t)
y20 (t) = −ω 2 sin(y1(t)) .
Die Gleichgewichtspunkte sind demnach y1k = k π, y2k = 0, k ∈ Z.
Die Linearisierung um diese Punkte ergibt:
! ! !
d y1 0 1 y1 − y1k
= + O(ky − yk k2) .
dt y2 −ω 2(−1)k 0 y2 − y2k
164
(
± iω, falls k gerade
λ1,2 =
±ω, falls k ungerade .
Für ungerade k ist ϕk also ein instabiler Sattelpunkt. Für
gerade k ist ϕk ein stabiler, jedoch kein asymptotisch stabiler
Gleichgewichtspunkt. Letzteres folgt nicht aus (8.19)!
↑φ
−π π
φ→
166
Satz (8.23) (Stabilitätssatz IV)
a) Ist V (y) eine Ljapunov–Funktion von f , so ist y∗ = 0 ein
gleichmäßig stabiler Gleichgewichtspunkt.
b) Ist V (y) sogar eine strenge Ljapunov–Funktion, so ist y∗ = 0
ein asymptotisch stabiler Gleichgewichtspunkt.
Beweis:
zu a) Zu 0 < ε ≤ R sei Mε := min{V (y) : kyk = ε} > 0.
Wegen der Stetigkeit von V existiert 0 < δ < ε mit
kyk ≤ δ ⇒ V (y ) < M ε .
Ist y0 ∈ Kδ (0) Anfangswert, so folgt für die zugehörige Trajek-
torie
d
V (y(t)) = ∇ V (y(t))T f (y(t)) ≤ 0 ,
dt
167
also ist V (y(t)) < Mε für alle t ≥ t0, und somit auch ky(t)k < ε
für t ≥ t0.
zu b) Ist V (y) strenge Ljapunov–Funktion, so fällt V (y(t)) für
t ≥ t0. Daher existiert E := limt→∞ V (y(t)). Wäre E > 0,
so gäbe es ein α, 0 < α < ε, mit V (y) < E für alle kyk < α.
n Trajektorie ganz in oα ≤ kyk ≤ ε.
Daher verläuft die Sei
nun M1 := max ∇ V (y)Tf (y) | α ≤ kyk ≤ ε , also nach Voraus-
setzung M1 < 0. Dann folgt:
Zt
d
V (y(t)) − V (y0) = V (y(τ )) dτ ≤ M1(t − t0)
dτ
t0
und damit lim V (y(t)) = −∞, im Widerspruch zu den Eigen-
t→∞
schaften von V . Damit ist lim V (y(t)) = 0, und damit auch
t→∞
lim y(t) = 0.
t→∞
168
Bemerkung (8.24)
Analog lässt sich zeigen: Gilt für eine C 1–Funktion V = V (y)
die Eigenschaft (8.22) a) sowie
h∇ V (y) , f (y)i > 0, ∀ y : 0 < kyk ≤ R ,
so ist y∗ = 0 ein instabiler Gleichgewichtspunkt.
170