fixpunkte und stabilitätsanalyse - itp.uni-hannover.de · alexander rachow, gregor efstradiadis...

Post on 30-Oct-2019

12 Views

Category:

Documents

0 Downloads

Preview:

Click to see full reader

TRANSCRIPT

Seite Alexander Rachow, Gregor Efstradiadis

Fixpunkte und Stabilitätsanalyse

1

Seite Alexander Rachow, Gregor Efstradiadis

Themenüberblick

Motivation 1D-Probleme Bifurkationen 2D-Probleme Fixpunkttypen Lotka-Volterra-Modelle

2

Seite Alexander Rachow, Gregor Efstradiadis

Motivation

Bisher: Lineare Dynamik Jetzt: Nichtlineare Systeme Ohne explizite Lösung Fixpunkte !!!

3

Seite Alexander Rachow, Gregor Efstradiadis

1D-Probleme

Betrachte Systeme der Form Einfaches Beispiel:

Explizite Lösung der DGL mit Separation

4

Seite Alexander Rachow, Gregor Efstradiadis

1D-Probleme Stattdessen: Graphische Analyse als Vektorfeld auf x-Achse

Graphische Darstellung:

i. Pfeil nach rechts

ii. Pfeil nach links Fixpunkt bei

5

Seite Alexander Rachow, Gregor Efstradiadis

1D-Probleme

6

Seite Alexander Rachow, Gregor Efstradiadis

1D-Probleme

Graphisches Stabilitätskriterium: i. Stabil Pfeile zeigen nach innen

ii. Instabil Pfeile zeigen nach außen

Quantitative Analyse auch möglich?

Lineare Stabilitätsanalyse

7

Seite Alexander Rachow, Gregor Efstradiadis

1D-Probleme

Ansatz: Taylor d

Lösung:

8

Seite Alexander Rachow, Gregor Efstradiadis

1D-Probleme Lineares Stabilitätskriterium:

i. Stabil

ii. Instabil

9

Seite Alexander Rachow, Gregor Efstradiadis

1D-Probleme

Zurück zum Beispiel Fixpunkt, wenn

10

Seite Alexander Rachow, Gregor Efstradiadis

1D-Probleme

Zurück zum Beispiel Fixpunkt, wenn

11

instabil

stabil

Seite Alexander Rachow, Gregor Efstradiadis

1D-Probleme

Erinnerung:

12

Seite Alexander Rachow, Gregor Efstradiadis

1D-Probleme

Bemerkung: Was passiert bei ?

13

stabil instabil

Seite Alexander Rachow, Gregor Efstradiadis

1D-Probleme

Hier ist halbstabiler Fixpunkt

14

Seite Alexander Rachow, Gregor Efstradiadis

Bifurkationen Übersicht: Sattel-Knoten-Bifurkationen Transkritische Bifurkationen Pitchfork-Bifurkationen

• Superkritische Bifurkationen • Subkritische Bifurkationen

Nicht perfekte Bifurkationen Beispiele

15

Seite Alexander Rachow, Gregor Efstradiadis

Was sind Bifurkationen

Bis jetzt: „simple“ 1D-Probleme • Dynamik der Vektorfelder sehr eingeschränkt

Nun: Abhängigkeit von Parametern • Struktur des Verlaufs ändert sich • Fixpunkte werden „zerstört“ oder „erschaffen“

Änderungen in der Dynamik heißen „Bifurkationen“ Parameterwerte an denen sie auftreten

„Bifurkationspunkt“ Bifurkationen liefern Modelle für Verstärkung und

Instabilität

16

Seite Alexander Rachow, Gregor Efstradiadis

Sattel-Knoten-Bifurkation

Fixpunkte für bei Bei Variation des Parameters bewegen sich

Fixpunkte aufeinander zu und annihilieren sich Normalform

unterschiedliche Fälle

17

Seite Alexander Rachow, Gregor Efstradiadis

Fall 1:

18

Seite Alexander Rachow, Gregor Efstradiadis

Fall 2:

19

Seite Alexander Rachow, Gregor Efstradiadis

Fall 3:

20

Seite Alexander Rachow, Gregor Efstradiadis

Bifurkationsdiagramm

21

Seite Alexander Rachow, Gregor Efstradiadis

Transkritische Bifurkationen

Es gibt Modelle in denen ein Fixpunkt für alle Parameterwerte existieren „muss“ Dieser Fixpunkt kann jedoch seine Stabilität

ändern Für solche Modelle werden transkritische

Bifurkationen verwendet Normalform

22

Seite Alexander Rachow, Gregor Efstradiadis

Fall 1:

23

Seite Alexander Rachow, Gregor Efstradiadis

Fall 2:

24

Seite Alexander Rachow, Gregor Efstradiadis

Fall 3:

25

Seite Alexander Rachow, Gregor Efstradiadis

Bifurkationsdiagramm

26

Seite Alexander Rachow, Gregor Efstradiadis

Pitchfork Bifurkation

Modell für Systeme mit Symmetrie Es gibt 2 Arten dieser Bifurkation

• Superkritische Pitchfork Bifurkation • Subkritische Pitchfork Bifurkation

27

Seite Alexander Rachow, Gregor Efstradiadis

Superkritische Pitchfork Bifurkation

Normalform invariant unter Kubischer Term wirkt „stabilisierend“

28

Seite Alexander Rachow, Gregor Efstradiadis

Fall 1:

29

Seite Alexander Rachow, Gregor Efstradiadis

Fall 2:

30

Seite Alexander Rachow, Gregor Efstradiadis

Fall 3:

31

Seite Alexander Rachow, Gregor Efstradiadis

Bifurkationsdiagramm

32

Seite Alexander Rachow, Gregor Efstradiadis

Subkritische Pitchfork Bifurkation

Normalform Kubischer Term ist hier nicht mehr

stabilisierend „treibt“ Funktion andere Normalform für „stabile Funktionen“

System invariant für erster stabilisierender Term ist neue Normalform

33

Seite Alexander Rachow, Gregor Efstradiadis

Bifurkationsdiagramm für „instabile Bifurkation“

34

Seite Alexander Rachow, Gregor Efstradiadis

Bifurkationsdiagramm für „stabile Bifurkation

35

Seite Alexander Rachow, Gregor Efstradiadis

Wichtiges zur „stabilen Bifurkation“

Für koexistieren zwei stabile Zustände • Welcher Fixpunkt für angenommen wird,

ergibt sich aus der Anfangsbed. • Daher ist der Fixpunkt im Ursprung gegen kleine

Störungen stabil, nicht jedoch gegen große der Ursprungsfixpunkt ist nur lokal stabil

36

Seite Alexander Rachow, Gregor Efstradiadis

Wichtiges zur „stabilen Bifurkation“

37

Die Existenz von versch. stabilen Zuständen ermöglicht Sprünge und Hysteresen

Seite Alexander Rachow, Gregor Efstradiadis

Nichtperfekte Bifurkationen

Normalerweise Systeme nicht perfekt symmetrisch, sondern nur näherungsweise Diese Syteme haben Unvollkommenheiten Betrachte hierzu :

• Für liegt normale Symmetrie vor • Für liegt Symmetriebrechung vor

Analyse hier deutlich schwieriger aufgrund zweier unabhängiger Parameter

38

Seite Alexander Rachow, Gregor Efstradiadis

Analyse der Fixpunkte

graphischer Ansatz: Plotte in ein Koordinatensyst.

Suche nach Schnittpunkten Schnittpunkte sind Fixpunkte

39

Seite Alexander Rachow, Gregor Efstradiadis

Analyse der Fixpunkte

Kritischer Fall: Gerade ist Tangente an Extrempunkt • Dann: Sattel-Knoten-Bifurkation

Suche Werte für h an denen Bifurkation auftritt: Für den Wert am lokalen Maximum gilt:

Sattel-Knoten-Bifurkation tritt auf , wenn

40

Seite Alexander Rachow, Gregor Efstradiadis

Analyse der Fixpunkte

41

Seite Alexander Rachow, Gregor Efstradiadis

Analyse der Fixpunkte

Bifurkationsdiagramme von für festes

42

Seite Alexander Rachow, Gregor Efstradiadis

Analyse der Fixpunkte

Bifurkationsdiagramme von für festes

43

Seite Alexander Rachow, Gregor Efstradiadis

Kurze Anmerkung zu Katastrophen

Spitzenkatastrophe

44

Seite Alexander Rachow, Gregor Efstradiadis

Beispiel : Populationsentwicklung

Hier: Falterart aus Kanada Modell bedient sich der „Trennung von

Zeiträumen“ • Population wächst schnell (char. Zeitraum:

Monate) • Bäume wachsen langsam (char. Zeitraum Jahre)

Bei Betrachtung der Populationsentwicklung können also die „Waldvariablen“ konstant angenommen werden

45

Seite Alexander Rachow, Gregor Efstradiadis

Populationsentwicklung

Für die Population gilt: wächst logistisch mit Wachstumsrate und

Tragfähigkeit (hierbei fest) (Todesrate aufgr. von Raubtieren) wächst

logistisch

46

Seite Alexander Rachow, Gregor Efstradiadis

Populationsentwicklung

Hier Untersuche Modell nun auf einen sog.Ausbruch

47

Seite Alexander Rachow, Gregor Efstradiadis

Dimensionlose Formulierung

48

Seite Alexander Rachow, Gregor Efstradiadis

Analyse der Fixpunkte

Erster Fixpunkt bei (immer instabil) Weitere Fixpunkte durch Lösung von

49

Seite Alexander Rachow, Gregor Efstradiadis

Analyse der Fixpunkte

50

Seite Alexander Rachow, Gregor Efstradiadis

Bifurkationskurven errechnen

Wir untersuchen Kurven im für die eine Sattel-Knoten-B. vorliegt Wir können jedoch die Parameter nicht als

Funktion voneinander ausdrücken Wähle Parametrisierung

51

Seite Alexander Rachow, Gregor Efstradiadis

Bifurkationskurven errechnen

Für S-K-B müssen erfüllt sein: Einsetzen von (3) in (1) liefert:

Einsetzen von (4) in (3) liefert:

Bifurkationskurven sind (5) und (4)

52

Seite Alexander Rachow, Gregor Efstradiadis

Bifurkationskurven

53

Seite Alexander Rachow, Gregor Efstradiadis

2D-Probleme

Betrachte nun nichtlineare Probleme vom Typ Fixpunkte dann gegeben durch

54

Seite Alexander Rachow, Gregor Efstradiadis

2D-Probleme

baut Phasenebene auf ist dann Vektorfeld auf der Phasenebene

55

Seite Alexander Rachow, Gregor Efstradiadis

2D-Probleme

Erst lineare Systeme untersuchen! Allgemeine Lösung:

Schreibe

Mit und

56

Seite Alexander Rachow, Gregor Efstradiadis

Sattelpunkt

57

Seite Alexander Rachow, Gregor Efstradiadis

Knoten

: Stabil : Instabil

58

Seite Alexander Rachow, Gregor Efstradiadis

Periodische Lösung

59

Seite Alexander Rachow, Gregor Efstradiadis

Spirale

: stabil : instabil

60

Seite Alexander Rachow, Gregor Efstradiadis

Sonderfall: Stern-Knoten

61

Seite Alexander Rachow, Gregor Efstradiadis

Sonderfall: Entarteter-Knoten

62

Seite Alexander Rachow, Gregor Efstradiadis

Zusammenfassung

Was haben wir bisher gelernt? Fixpunkte zur Betrachtung der Dynamik Stabilität graphisch oder analytisch

bestimmen Graphisch: Vektorflusses visualisieren Analytisch: Lineare Stabilitätsanalyse In 2D viel mehr Möglichkeiten als in 1D

63

Seite Alexander Rachow, Gregor Efstradiadis

Lotka-Volterra-Modelle

Beschreiben Wechselwirkung von mehreren Populationen Beispiel: Räuber-Beute-Problem Populationszahlen: (Beute), (Räuber) Ungestörte Wachstumsrate:

Begegnung Räuber-Beute:

64

Seite Alexander Rachow, Gregor Efstradiadis

Lotka-Volterra-Modelle

Insgesamt also: Betrachte die Dynamik mit „Mathematica“!

65

Seite Alexander Rachow, Gregor Efstradiadis

Quellenverzeichnis

Strogatz, Steven H. (1994): Nonlinear Dynamics and Chaos. Massachusetts

Suter, Dieter (2010): Skript zur Vorlesung „Analytische Mechanik“. [cited 02.06.2015]. https://e3.physik.uni-dortmund.de/~suter/Vorlesung/Physik_III_WS10/2.8_Chaos.pdf

[cited 29.05.2015] http://de.wikipedia.org/wiki/Kritischer_Punkt_(Dynamik)

66

top related