Kurs:Mathematische Modellbildung/Themen/Räuber-Beute-Modelle/Implementation - Uni

Modellierungszyklus 3 - Niveau Uni[Bearbeiten]

Auswahl der Software[Bearbeiten]

Ziel: die gegenseitige Beeinflussung zweier Populationen und Regelmäßigkeiten in der Beziehung der Populationen darstellen → Anwendung von Lotka-Volterra-Gleichungen

Tabellenkalkulation: Diskretisierung Lotka-Volterra (rekursive, schrittweise Berechnung und Plot)
Maxima: Ableitungsberechnungen, Hessematrix, Eigenwerte (Gleichgewichtszustand als lokales Minimum)
Octave: Vektorfeld zu gekoppelten Differentialgleichungen, Diskretisierung (Darstellung als Orbit), Höhenlinien

Diskretisierung Lotka-Volterra[Bearbeiten]

Fichtenpopulation

Für die Änderungsrate der Fichtenpopulation wurde sowohl die Vermehrung der Fichten als auch die Schädigung der Fichten durch die Borkenkäfer betrachtet. Es ergibt sich folgender Zusammenhang:

$F'(t)=a\cdot F(t)-b\cdot F(t)\cdot B(t)$

F'(t): Änderungsrate Fichten
a: Reproduktionsrate Fichten
b: Sterberate Fichten pro Borkenkäfer
F(t): Anzahl Fichten zum Zeitpunkt t
B(t): Anzahl Borkenkäfer zum Zeitpunkt t

Änderungsrate der Fichten[Bearbeiten]

Betrachtung der Vermehrung der Fichten, die sich aus der Multiplikation der Reproduktionsrate mit der aktuellen Fichtenanzahl ergibt.
Der Subtrahend berücksichtigt, wie tödlich ein Zusammentreffen von Fichte und Borkenkäfer ist. Die getöteten Fichten werden abgezogen.

Borkenkäferpopulation[Bearbeiten]

Für die Änderungsrate der Borkenkäferpopulation wurden sowohl die Sterberate der Borkenkäfer als auch die Reproduktionsrate der Borkenkäfer pro Fichte berücksichtigt. Es ergibt sich folgender Zusammenhang:

$B'(t)=-c\cdot B(t)+d\cdot F(t)\cdot B(t)$

B'(t): Änderungsrate Borkenkäfer
c: Sterberate Borkenkäfer
d: Reproduktionsrate Borkenkäfer pro Fichte
F(t): Anzahl Fichten zum Zeitpunkt t
B(t): Anzahl Borkenkäfer zum Zeitpunkt t

Änderungsrate der Borkenkäfer[Bearbeiten]

Sterberate der Borkenkäfer: Multiplikation der Sterberate mit der aktuellen Borkenkäferanzahl. Da die Population durch Sterbefälle zurückgeht, wird die Konstante c mit einem negativen Vorzeichen versehen.
zweiter Summand: wie nahrhaft ist ein Zusammentreffen von Fichte und Borkenkäfer und wie wirkt sich dies auf die Reproduktion der Borkenkäfer aus.
Zusammentreffen von Borkenkäfer und Fichte: Reproduktion der Borkenkäfer steigt → Term hat positiven Einfluss auf die Änderung

Wahl der Parameter[Bearbeiten]

Parameter a (Reproduktionsrate Fichten)[Bearbeiten]

$a={1 \over 30}\cdot 10={1 \over 3}$
Eine Fichte braucht im Durchschnitt 30 Jahre, um geschlechtsreif zu werden
pro Fichte entwickeln sich 10 der Samen auch zu einem neuen Baum

Parameter c (Sterberate der Borkenkäfer)[Bearbeiten]

$c={1 \over 2}$
Lebenserwartung eines Borkenkäfers beträgt 2 Jahre

Parameter b und d[Bearbeiten]

schwer zu schätzen, deshalb ist eine weitere Annahme für die Berechnung dieser Parameter nötig :

In einem Ökosystem gibt es immer einen Gleichgewichtszustand (genauso viele Individuen sterben wie geboren werden) → keine Änderung der Populationsgröße (Änderungsrate=0)

Annahmen für Gleichgewichtszustand[Bearbeiten]

$B_{Gleichgewicht}=200000$

$F_{Gleichgewicht}=100000$ (in Tausend)

Durch Nullsetzen und Umstellen der Änderungsraten erhält man:

$b={a \over B_{Gleichgewicht}}={{1 \over 3} \over 200000}\approx 1,66\cdot 10^{-6}$

$d={c \over F_{Gleichgewicht}}={{1 \over 2} \over 100000}=0,000005$

Festlegung der Schrittweite[Bearbeiten]

Es wurde eine Schrittweite von 0,1 festgelegt, um möglichst viele Punkte in kleinen Abständen zu erhalten und einen ersten Plot durchführen zu können.

Dabei ergibt sich:

$F(t)=F(t-1)+0,1\cdot (a\cdot F(t-1)-b\cdot F(t-1)\cdot B(t-1))$

$B(t)=B(t-1)+0,1\cdot (-c\cdot B(t-1)+d\cdot F(t-1)\cdot B(t-1))$

Aus Zyklus 1 wird weiterhin verwendet:[Bearbeiten]

$F(0)=65520000$
$B(0)=102080$

Lotka-Volterra-Gleichung[Bearbeiten]

Die Anzahl der Fichten wurde durch 1 000 geteilt, um beide Graphen in einem Koordinatensystem darstellen zu können.

Vektorfeld und Orbitdarstellung des zyklischen Verlaufs mit Octave[Bearbeiten]

Darstellung der Differentialgleichung als Vektorfeld mit Octave
Hierbei stellen die Vektoren mit ihrem Anteil in x- und y-Richtung die Änderung an dem jeweiligen Punkt dar und der zyklische Verlauf der Populationsentwicklung lässt sich an den sich andeutenden Orbits bereits erkennen.