04.03.2013 Aufrufe

Blatt 12 - Institut für Mathematik - TU Berlin

Blatt 12 - Institut für Mathematik - TU Berlin

Blatt 12 - Institut für Mathematik - TU Berlin

MEHR ANZEIGEN
WENIGER ANZEIGEN

Sie wollen auch ein ePaper? Erhöhen Sie die Reichweite Ihrer Titel.

YUMPU macht aus Druck-PDFs automatisch weboptimierte ePaper, die Google liebt.

Technische Universität <strong>Berlin</strong> Wintersemester 2009/2010<br />

<strong>Institut</strong> <strong>für</strong> <strong>Mathematik</strong><br />

Prof. Dr. P. Wittbold<br />

Dr. Frank Jochmann<br />

http://www.math.tu-berlin.de/∼wittbold/<br />

<strong>12</strong>. Übungsblatt zu Differentialgleichungen I<br />

(Stabilität, ω-Limes Menge, Lyapunov-Funktion, Abgabe: bis Fr. 29.01. in den Tutorien)<br />

Übung<br />

1. Aufgabe<br />

Seien F : R N → R N lokal Lipschitz-stetig und x : [0, ∞) → R N eine beschränkte Lösung der<br />

Dgl. x ′ = F(x).<br />

Zeigen Sie, dass die ω-Limes Menge zusammenhängend ist.<br />

Definition: Seien F : D ⊂ RN → RN lokal Lipschitz-stetig und <strong>für</strong> z ∈ D sei<br />

φ(z, ·) : [0, T + (z)) → D die maximale Lösung des rechtsseitigen AWP x ′ = F(x), x(0) = z.<br />

Eine Menge M ⊂ D heisst Attraktor, wenn M invariant unter dem Fluss φ ist und jeder<br />

Punkt in M innerer Punkt des Einzugsbereiches E(M) von M ist. (Der Einzugsbereich ist die<br />

Menge aller Punkte z ∈ D mit T + (z) = ∞ und dist (φ(z, t), M) t→∞<br />

−→ 0.)<br />

2. Aufgabe<br />

Zeigen Sie, dass der Einzugsbereich E(M) eines Attraktors M ⊂ D offen ist (und somit E(M)<br />

eine offene Umgebung von M ist).<br />

3. Aufgabe<br />

Seien A ∈ R n×n schiefadjungiert, d.h. A ∗ = −A, r > 0 und g ∈ C 1 (R, R) mit g(r 2 ) = 0 und<br />

g(u) > 0 <strong>für</strong> alle u ∈ (0, ∞) \ {r 2 }.<br />

Bestimmen Sie den Einzugsbereich von<br />

Ist Kr ein Attraktor?<br />

F(x) def<br />

= Ax − g(|x| 2 )x<br />

Kr = {x ∈ R n : |x| = r}.<br />

4. Aufgabe<br />

Suppose that f : R N → R N is a locally Lipschitz-continuous function, such that 0 is an equilibrium<br />

point of the differential equation<br />

x ′ = F(x).<br />

Furthermore, let V : R n → R be a Lyapunov-function for this system with<br />

V (y) → ∞ <strong>für</strong> |y| → ∞.<br />

1. Show that, for each y ∈ R n , the initial value problem<br />

x ′ = F(x), x(0) = y, t ∈ [0, ∞)<br />

has a unique global solution φ(y, ·) : [0, ∞) → R n .


2. Suppose that, in addition,<br />

Show that, for each y ∈ R n ,<br />

V (0) = 0, V (y) > 0 and ∇V (y) · f(y) < 0 ∀y = 0.<br />

φ(y, t) t→∞<br />

−→ 0.<br />

Tutoriumsaufgaben<br />

1. Aufgabe<br />

Für α > 0, β > 0, γ > 0, δ > 0 betrachten wir das Volterra-Lotka-System<br />

x ′ = F(x) def<br />

<br />

x1(α − βx2)<br />

=<br />

.<br />

(δx1 − γ)x2<br />

1. Zeigen Sie, dass dieses System die stationären Punkte (0, 0) und z = (γ/δ, α/β) hat.<br />

2. Untersuchen Sie diese stationären Punkte auf Stabilität.<br />

Hinweis : Betrachten Sie die Lyapunov-Funktion<br />

2. Aufgabe<br />

V (x) def<br />

=<br />

x1<br />

z1<br />

u −1 (δu − γ)du +<br />

x2<br />

z2<br />

u −1 (βu − α)du.<br />

Sei D ⊂ R N offen und f : D → R N lokal Lipschitz stetig. Sei K ⊂ D eine konvexe, kompakte<br />

und positiv invariante Teilmenge unter dem von dem System x ′ = F(x) erzeugten Fluss. Zeigen<br />

Sie, dass K mindestens eine Ruhelage enthält.<br />

3. Aufgabe<br />

Seien<br />

F(x) def<br />

=<br />

Bestimmen Sie den Einzugsbereich von<br />

und zeigen Sie, dass K1 ein Attraktor ist.<br />

x2<br />

<br />

+ (1 − |x|<br />

−x1<br />

2 )x.<br />

K1 = {x ∈ R 2 : |x| = 1}<br />

4. Aufgabe<br />

Seien F : D ⊂ R N → R N lokal Lipschitz-stetig und <strong>für</strong> z ∈ D sei φ(z, ·) : [0, T + (z)) → D die<br />

maximale Lösung des rechtsseitigen AWP x ′ = F(x), x(0) = z.<br />

Beweisen oder widerlegen Sie durch ein Gegenbeispiel:<br />

• Für jede nichtleere, positiv invariante Menge M ⊂ D ist auch M positiv invariant.<br />

• Sei x0 ∈ D stabiler Gleichgewichtspunkt. Dann ist E({x0}) offen.<br />

• Sei x0 ∈ D asymptotisch stabiler Gleichgewichtspunkt. Dann ist {x0} ein Attraktor.<br />

• Für jede nichtleere Menge M ⊂ D und alle z ∈ E(M) gilt: ω(z) ⊂ M. Hier ist ω(z) die<br />

ω-Limes Menge.


Hausaufgaben<br />

1. Aufgabe 7 Punkte<br />

Seien<br />

F(x) def<br />

<br />

x2<br />

= + |x|<br />

−x1<br />

2<br />

<br />

π<br />

sin x<br />

|x|<br />

falls x ∈ R N \ {0} und F(0) def<br />

= 0.<br />

• Zeigen Sie, dass F auf R N stetig differenzierbar ist.<br />

• Untersuchen Sie die Nulllösung von x ′ = F(x) auf Stabilität.<br />

• Untersuchen Sie die Nulllösung von x ′ = F(x) auf asymptotische Stabilität.<br />

2. Aufgabe 8 Punkte<br />

Für α ∈ R sei<br />

F(x) def<br />

<br />

<br />

=<br />

.<br />

x1x2<br />

α(x1 − 1) 3 (x2 − 1) 3<br />

1. Untersuchen Sie den stationären Punkt y = (0, 1) des Systems x ′ = F(x) auf Stabilität in<br />

Abhängigkeit von α.<br />

2. Untersuchen Sie den stationären Punkt z = (1, 0) auf Stabilität in Abhängigkeit von α im<br />

Fall α = 0.<br />

Hinweis zu (2): Betrachten Sie die Lyapunov-Funktion<br />

Vα(x) def<br />

= α<br />

x1<br />

u<br />

1<br />

−1 (u − 1) 3 x2<br />

du +<br />

0<br />

u(1 − u) −3 du.<br />

3. Aufgabe 5 Punkte<br />

Bestimmen Sie zu gegebenem α ∈ R die ω-Limesmenge <strong>für</strong> die Lösung des Anfangswertproblems<br />

wobei<br />

x ′ = Aαx, x(0) = (1, 0)<br />

Aα<br />

def<br />

=<br />

<br />

−α 1<br />

.<br />

−2 −α

Hurra! Ihre Datei wurde hochgeladen und ist bereit für die Veröffentlichung.

Erfolgreich gespeichert!

Leider ist etwas schief gelaufen!