Gew¨ohnliche Differentialgleichungen und Dynamische Systeme
Gew¨ohnliche Differentialgleichungen und Dynamische Systeme Gew¨ohnliche Differentialgleichungen und Dynamische Systeme
nicht leer ist. Wenn x ∈ ∩ n≥0 I s0 ...s n , dann gilt x ∈ I s0 , f µ (x) ∈ I s1 , usw. und somit h(x) = (s 0 s 1 . . .), woraus die Surjektivität folgt. Um die Stetigkeit von h zu zeigen, wählen wir x ∈ Λ mit h(x) = s 0 s 1 . . .. Sei ɛ > 0 und wähle n ∈ N mit 1/2 n < ɛ. Die Menge der Intervalle I t0 ...t n für alle möglichen Kombinationen von t 0 . . . t n sind disjunkt und ihre Vereinigung enthält Λ. Es gibt 2 n+1 solche Intervalle und I s0 ...s n ist eines von ihnen. Wir wählem nun δ > 0, so daß aus |x − y| < δ und y ∈ Λ die Beziehung y ∈ I s0 ...s n folgt. Damit gilt aber, daß h(x) und h(y) in der ersten n + 1 Stellen übereinstimmen und somit d(h(x), h(y)) < 1 2 < ɛ. n Damit ist die Stetigkeit von h gezeigt. Analog folgt die Stetigkeit von h −1 . Es bleibt die Konjugiertheit der Flüsse zu zeigen. □ Theorem 4.28 Es gilt h ◦ f µ = σ ◦ h. Beweis: Ein Punkt x ∈ Λ ist eindeutig durch ∩ n≥0 I s0 ...s n festgelegt. Es ist so daß da f µ (I s0 ) = I. Es ergibt sich I s0 ...s n = I s0 ∩ f −1 µ (I s1 ) ∩ . . . ∩ f −n µ (I sn ), f µ (I s0 ...s n ) = I s1 ∩ f −1 µ (I s 2 ) ∩ . . . ∩ f −(n−1) µ (I sn ) = I s1 ...s n , h ◦ f µ (x) = h ◦ f µ (∩ ∞ n=0 I s 0 ...s n ) = h(∩ ∞ n=1I s1 ...s n ) = s 1 s 2 . . . = σ ◦ h(x). Zusammenfassend haben wir gezeigt: □ Theorem 4.29 Sei f µ (x) = µx(1 − x) mit µ > 2 + √ 5. Dann gilt: a) Es gibt 2 n periodische Punkte der Periode n. b) Die Menge der periodischen Punkte liegt dicht in Λ. c) f µ besitzt einen dichten Orbit in Λ. Bemerkung: Der Nachweis der periodischen Punkte kann auch direkt erfolgen, denn der Graph von f n µ schneidet die Gerade y = x mindestens 2n -mal. Damit hat f n µ mindestens 2n Fixpunkte, bzw. f µ besitzt mindestens 2 n periodische Punkte der Periode n. Bemerkung: Die Konjugiertheit zur Shiftdynamik ist die strengste Definition von Chaos. Für praktische Zwecke gibt es weichere Definitionen, siehe [SU03]. In zwei Raumdimensionen erfolgt der Nachweis chaotischen Verhaltens of über ein Zwischenmodell, das Smale’sche Hufeisen, siehe Abschnitt 6.2. 52
Einschub: Periodenverdopplung und die Feigenbaumkaskade Für µ ≤ 3 existiert nach Lemma 4.20 als einzige periodische Lösung ein stabiler Fixpunkt. Für µ > 3 erhalten wir eine Kaskade von Periodenverdopplungen. Übung 4.30 Untersuche f µ mittels xppaut, insbesondere für µ ∈ (3, 4). Was passiert für µ > 4, was für 3.825 ≤ µ ≤ 3.86 ? Lösung. Eine Möglichkeit Übung 1.4 zu bearbeiten ist die folgende: kopiere das folgende xppaut- ODE file qf.ode 2 v m x # m=mu ist dummy-variable, d.h. parameter o m o m*x*(1-x) d in ein Verzeichnis und starte dort xppaut, d.h. xppaut qf.ode. Dann führe in xppaut die folgenden Schritte aus: 1) unter Viewaxes wähle m ∈ [2.5, 4.5] für x-axes und x ∈ [0.1] für y-axes. 2) unter Numerics, Method wähle discrete, sowie Total=200 und tRansient=150; hiermit (wie mit allem anderen) kann man spielen. 3) unter Graphic stuff wähle (E)dit und dann Line type: 0. 4) unter Initialconds wähle Mice und klick im Fenster rum. 4b) unter Initialconds wähle New und dann x = 0.1 (z.B.). Anschließend wähle Initialconds, Range und dann m ∈ [2.5, 4.5], 200 Steps. Im Graphikfenster von xppaut sollte nun in etwa das Bild in Abb. 19 links erscheinen. x 8 4 2 1 µ Abbildung 19: Die Feigenbaumkaskade f µ , numerisch und schematisch. Die numerische Simulation zeigt mit wachsendem µ eine Reihe von Periodenverdopplungen, siehe Abb.19. Die Werte µ j , bei denen diese stattfinden, häufen sich bei µ = µ ∞ . Überraschen- 53
- Seite 1 und 2: Gewöhnliche Differentialgleichunge
- Seite 3 und 4: Inhaltsverzeichnis 1 Einleitung 4 2
- Seite 5 und 6: Im nichtautonomen Fall lautet dies
- Seite 7 und 8: Geschichte: Die Beschreibung der Me
- Seite 9 und 10: ii) ‖λx‖ = |λ|‖x‖ iii)
- Seite 11 und 12: Ist g ≡ 0, so gilt: x(t) = c 1 x
- Seite 13 und 14: Die Exponentialreihe e At kann wie
- Seite 15 und 16: 2 1.5 1 0.5 0 -0.5 -1 -1.5 -2 -2 -1
- Seite 17 und 18: ε δ x=x(t) x 0 x 1 Abbildung 4: S
- Seite 19 und 20: Betrachte das autonome System ẋ =
- Seite 21 und 22: Dieses Lemma gilt allgemein für T
- Seite 23 und 24: Unteraum. Wie bei Differentialgleic
- Seite 25 und 26: 3 Nichtlineare Systeme Wir betracht
- Seite 27 und 28: nur aufhören zu existieren, falls
- Seite 29 und 30: Ist f differenzierbar, so hängt di
- Seite 31 und 32: Beweis: Die Differenz r(t) = x(t)
- Seite 33 und 34: Theorem 3.16 Sei f : R d × [t 0 ,
- Seite 35 und 36: und für die gesuchte Näherung spe
- Seite 37 und 38: und somit die Konvergenz der Lösun
- Seite 39 und 40: Beweis: Zum Nachweis dieses Satzes
- Seite 41 und 42: x µ Abbildung 11: Transkritische B
- Seite 43 und 44: kann für α = 0 explizit gelöst w
- Seite 45 und 46: O − (x) eines Punktes x ∈ R dur
- Seite 47 und 48: 1 0.8 0.6 0.4 0.2 0 0 0.2 0.4 0.6 0
- Seite 49 und 50: zu a) In diesem Fall ist p ∈ Λ n
- Seite 51: iii) Sei a ∈ Σ 2 und δ > 0 gege
- Seite 55 und 56: Theorem 4.33 Perioden-Verdopplungs-
- Seite 57 und 58: esitzt n-periodische Lösungen y n
- Seite 59 und 60: ii) Die periodische Lösung heißt
- Seite 61 und 62: Beispiel 4.43 (Pitchfork-Bifurkatio
- Seite 63 und 64: 5 Dynamik in der Nähe eines Fixpun
- Seite 65 und 66: und Damit ist T (·, f) : Cb 0 →
- Seite 67 und 68: Offensichtlich ist W s = E s . Für
- Seite 69 und 70: erhalten wir W c = {(µ, x, y) ∈
- Seite 71 und 72: Lemma 5.20 Es sei ˜f ∈ C k für
- Seite 73 und 74: Lemma 5.25 Für jedes η ∈ (0, β
- Seite 75 und 76: Beispiel 5.30 Wir betrachten das di
- Seite 77 und 78: Wir suchen nun Koordinatentransform
- Seite 79 und 80: Zu α 30 a 3 : Die Nichtresonanzbed
- Seite 81 und 82: 6 Homokline und heterokline Lösung
- Seite 83 und 84: Abbildung 24: Schnitt der Lösung m
- Seite 85 und 86: mit 0 ≤ ɛ ≪ 1 ein kleiner Para
- Seite 87 und 88: q ε p ε Abbildung 27: Unendlich v
- Seite 89 und 90: i) Λ enthält eine abzählbare Men
- Seite 91 und 92: wobei q u ɛ (t 0 ) = q u ɛ (t 0 ,
- Seite 93 und 94: Bemerkung 6.21 Durch die Koordinate
- Seite 95 und 96: (x,y,z ) 1 q q (r 0 θ ,z) p p Abbi
- Seite 97 und 98: Wir wollen nun eine kurze Zusammenf
- Seite 99 und 100: [Ise96] [Kel67] [KK74] [Kre85] Arie
Einschub: Periodenverdopplung <strong>und</strong> die Feigenbaumkaskade<br />
Für µ ≤ 3 existiert nach Lemma 4.20 als einzige periodische Lösung ein stabiler Fixpunkt. Für<br />
µ > 3 erhalten wir eine Kaskade von Periodenverdopplungen.<br />
Übung 4.30 Untersuche f µ mittels xppaut, insbesondere für µ ∈ (3, 4). Was passiert für<br />
µ > 4, was für 3.825 ≤ µ ≤ 3.86 ?<br />
Lösung. Eine Möglichkeit Übung 1.4 zu bearbeiten ist die folgende: kopiere das folgende xppaut-<br />
ODE file qf.ode<br />
2<br />
v m x<br />
# m=mu ist dummy-variable, d.h. parameter<br />
o m<br />
o m*x*(1-x)<br />
d<br />
in ein Verzeichnis <strong>und</strong> starte dort xppaut, d.h. xppaut qf.ode. Dann führe in xppaut die folgenden<br />
Schritte aus:<br />
1) unter Viewaxes wähle m ∈ [2.5, 4.5] für x-axes <strong>und</strong> x ∈ [0.1] für y-axes.<br />
2) unter Numerics, Method wähle discrete, sowie Total=200 <strong>und</strong> tRansient=150; hiermit<br />
(wie mit allem anderen) kann man spielen.<br />
3) unter Graphic stuff wähle (E)dit <strong>und</strong> dann Line type: 0.<br />
4) unter Initialconds wähle Mice <strong>und</strong> klick im Fenster rum.<br />
4b) unter Initialconds wähle New <strong>und</strong> dann x = 0.1 (z.B.). Anschließend wähle<br />
Initialconds, Range <strong>und</strong> dann m ∈ [2.5, 4.5], 200 Steps. Im Graphikfenster von xppaut sollte<br />
nun in etwa das Bild in Abb. 19 links erscheinen.<br />
x<br />
8<br />
4<br />
2<br />
1<br />
µ<br />
Abbildung 19: Die Feigenbaumkaskade f µ , numerisch <strong>und</strong> schematisch.<br />
Die numerische Simulation zeigt mit wachsendem µ eine Reihe von Periodenverdopplungen,<br />
siehe Abb.19. Die Werte µ j , bei denen diese stattfinden, häufen sich bei µ = µ ∞ . Überraschen-<br />
53