Konstante Lösungen

Die konstanten Lösungen (auch Gleichgewichtspunkte oder kritische Punkte genannt) erhält man, indem man die rechten Seiten der Lotka-Volterra-Gleichungen gleich Null setzt:

${\displaystyle N(a-bP)=0,\qquad P(cN-d)=0.}$

Es gibt also genau zwei konstante Lösungen, nämlich den trivialen Gleichgewichtspunkt ${\displaystyle (0,0)}$ und den inneren Gleichgewichtspunkt

${\displaystyle (N^{*},P^{*})=\left({\frac {d}{c}},{\frac {a}{b}}\right).}$

Ein erstes Integral

Eine Methode zum Auffinden nicht-konstanter Lösungen besteht darin, ein erstes Integral, also eine Invariante der Bewegung, zu suchen. Volterra findet eine solche auf folgendem Weg: Multipliziert man die erste Grundgleichung mit ${\displaystyle c}$ und die zweite mit ${\displaystyle b}$, und addiert anschließend die beiden Gleichungen, so verschwinden die Terme mit dem Produkt ${\displaystyle NP}$, und man erhält

${\displaystyle c\,{\frac {\mathrm {d} N}{\mathrm {d} t}}+b\,{\frac {\mathrm {d} P}{\mathrm {d} t}}=ac\,N-bd\,P}$.

Durch Multiplikation der ersten Grundgleichung mit ${\displaystyle d/N}$ und der zweiten mit ${\displaystyle a/P}$ und anschließender Addition kommt man zu

${\displaystyle {\frac {d}{N}}\,{\frac {\mathrm {d} N}{\mathrm {d} t}}+{\frac {a}{P}}\,{\frac {\mathrm {d} P}{\mathrm {d} t}}=ac\,N-bd\,P}$.

Subtraktion dieser beiden Gleichungen ergibt

${\displaystyle c\,{\frac {\mathrm {d} N}{\mathrm {d} t}}-{\frac {d}{N}}\,{\frac {\mathrm {d} N}{\mathrm {d} t}}+b\,{\frac {\mathrm {d} P}{\mathrm {d} t}}-{\frac {a}{P}}\,{\frac {\mathrm {d} P}{\mathrm {d} t}}=0}$.

Durch Integration dieser letzten Gleichung erreicht man schließlich die Beziehung

${\displaystyle V(N,P):=c\,N-d\,\ln N+b\,P-a\,\ln P=}$ const.

Umgekehrt kann man die totale Ableitung der so definierten Funktion ${\displaystyle V}$ nach ${\displaystyle t}$ berechnen:

${\displaystyle {\begin{aligned}{\frac {\mathrm {d} }{\mathrm {d} t}}V(N,P)&={\frac {\partial V(N,P)}{\partial N}}\cdot {\frac {\mathrm {d} N}{\mathrm {d} t}}+{\frac {\partial V(N,P)}{\partial P}}\cdot {\frac {\mathrm {d} P}{\mathrm {d} t}}\\&=\left(c-{\frac {d}{N}}\right)\,N(a-bP)+\left(b-{\frac {a}{P}}\right)\,P(cN-d)\\&=0;\end{aligned}}}$

so gelangt man ebenfalls zu der Aussage, dass ${\displaystyle V}$ auf den Lösungen der Grundgleichungen konstant (invariant) ist; eine Lösung der Lotka-Volterra-Gleichung kann also ihre Niveaulinien von ${\displaystyle V}$ nicht verlassen.

Ein anderer Weg zum Auffinden einer Invarianten der Bewegung besteht darin, die Lotka-Volterra-Gleichungen mit Hilfe eines eulerschen Multiplikators in eine exakte Differentialgleichung umzuformen und diese dann zu integrieren.

Stabilität

Da ${\displaystyle V}$ als erstes Integral auch eine Ljapunow-Funktion ist, und da ${\displaystyle V}$ am inneren Gleichgewichtspunkt ein striktes lokales Minimum besitzt, folgt aus dem ersten Kriterium von Ljapunow, dass dieser Gleichgewichtspunkt stabil ist.

Periodizität

Da die Funktion ${\displaystyle V}$ im Quadranten ${\displaystyle Q_{+}}$ strikt konvex ist und ihr Minimum im inneren Gleichgewichtspunkt annimmt, bilden die Niveaulinien von ${\displaystyle V}$ geschlossene Kurven im Phasenraum. Da jede Lösung in einer Niveauline von ${\displaystyle V}$ enthalten sein muss, folgt aus der Eindeutigkeit und einer Betrachtung des lokalen Richtungsfeldes die Periodizität der Lösungen.

Loi du cycle périodique. – Les fluctuations des deux espèces sont périodique.

(Gesetz der periodischen Zykel: Die Schwankungen der beiden Klassen sind periodisch.)

Erhaltung der Mittelwerte

Aus der Periodizität der Lösungen folgt mit ein paar Zeilen Rechnung das

Loi de la conservation des moyennes. – Les moyennes pendant une période des nombres des individus des deux espèces sont independent des conditions initiales, et égales aux nombres qui correspondent à l'ètat stationnaire, pour les valeurs données des «coefficients d’accroissement» ${\displaystyle \epsilon _{1},\epsilon _{2}}$, et des «coefficients de voracité» ${\displaystyle \gamma _{1},\gamma _{2}}$.

(Gesetz der Erhaltung der Mittelwerte: Die zeitlichen Mittelwerte über eine Periode der Individuenzahlen der beiden Klassen hängen nicht von den Anfangsbedingungen ab, und sind gleich den Anzahlen des Gleichgewichtszustands, für die gegebenen „Wachstumskoeffizienten“ ${\displaystyle \epsilon _{1},\epsilon _{2}}$, und die „Fresskoeffizienten“ ${\displaystyle \gamma _{1},\gamma _{2}}$.)

Das bedeutet, die zeitlichen Mittelwerte erfüllen die Gleichungen

${\displaystyle {\overline {N_{1}}}={\frac {\epsilon _{2}}{\gamma _{2}}}\quad {\mbox{und}}\quad {\overline {N_{2}}}={\frac {\epsilon _{1}}{\gamma _{1}}}.}$

Auf den ersten Blick verwirrend ist hier, dass der Mittelwert der Beutetierpopulation ${\displaystyle N_{1}(t)}$ nur von Sterbe- und Fressrate der Raubtierpopulation und nicht von der Reproduktionsrate der Beutetiere abhängt. Dagegen ist der Mittelwert der Raubtierpopulation ${\displaystyle N_{2}(t)}$ nur von Reproduktions- und Sterberate der Beutetierpopulation und nicht von Fress- und Sterberate der Raubtiere abhängig. Dabei ist die gleichgewichtige Anzahl der Beute umso höher, je ungünstiger die Parameter für die Raubtiere sind. Die gleichgewichtige Anzahl der Räuber ist dagegen umso höher, je günstiger die Parameter für die Beutetiere sind.

Verständlich wird diese Eigenschaft des Lotka-Volterra-Modells, wenn man sich die in der mathematischen Begründung zur Anwendung kommende Modellbildung ansieht: die Kontrolle über die Populationszahl der einen Tierklasse obliegt hier ausschließlich der jeweils anderen Klasse.

Störung der Mittelwerte

Das wegen seiner biologischen Interpretation interessanteste dieser Gesetze ist das

Loi de la perturbation des moyennes. – Si l’on détruit les deux espèces uniformément et proportionnellement aux nombres de leurs individus (assez peu pour que les fluctuations subsistent), la moyenne du nombres des individus de l’espèce dévorée croît et celles de l’espèce devorant diminue.

(Gesetz der Störung der Mittelwerte: Werden die beiden Tierklassen gleichmäßig und proportional zur Populationsgröße dezimiert, und ist die Rate der Dezimierung der Beute kleiner als deren Reproduktionsrate, so wächst der Mittelwert der Beutepopulation, und der Mittelwert der Räuberpopulation sinkt.)

Tatsächlich beweist Volterra eine quantitative Version: Ist ${\displaystyle \alpha }$ die Zerstörungsrate der Beutelebewesen, und ${\displaystyle \beta }$ die Zerstörungsrate der Räuber, so sind die Mittelwerte für die Lösungen der gestörten Lotka-Volterra-Gleichungen

${\displaystyle {\overline {N_{1}}}={\frac {\epsilon _{2}+\beta }{\gamma _{2}}}\quad {\mbox{und}}\quad {\overline {N_{2}}}={\frac {\epsilon _{1}-\alpha }{\gamma _{1}}}.}$

Das bedeutet: die über eine Lotka-Volterra-Periode gemittelte Anzahl der Beutelebewesen steigt genau dann an, wenn die Räuber dezimiert werden – ziemlich unabhängig von einer Dezimierung der Beute, solange diese nicht ausgerottet wird. Umgekehrt sinkt die mittlere Anzahl der Räuber immer dann, wenn die Beutelebewesen dezimiert werden, und dieses Absinken hängt nicht davon ab, wie stark die Räuber zusätzlich dezimiert werden (solange diese nicht ausgerottet werden).

Intraspezifische Konkurrenzterme

Eine erste Erweiterung der Lotka-Volterra-Gleichungen entsteht durch Subtraktion von Termen proportional zu ${\displaystyle N_{1}^{2}}$ bzw. ${\displaystyle N_{2}^{2}}$, die die intraspezifische Konkurrenz modellieren. Es gibt verschiedene Möglichkeiten, die Form ${\displaystyle -\alpha _{1}N_{1}^{2}}$ und ${\displaystyle -\alpha _{2}N_{2}^{2}}$ der neu hinzugekommenen Terme zu begründen:

• Mit den empirischen Untersuchungen zur Bevölkerungsentwicklung nach Pierre-François Verhulst, siehe logistische Gleichung.
• Durch die Annahme, die (ungestörte) Wachstumsrate einer Population sei proportional zur Differenz zwischen einer Kapazitätsgrenze und der tatsächlichen Populationszahl.
• Durch eine Analyse des Einflusses von intraspezifischen Begegnungen auf die Populationszahl ähnlich Volterras Begründung des Terms zur Modellierungen der interspezifischen Konkurrenz: Eine intraspezifische Begegnung ist mit einer gewissen Wahrscheinlichkeit eine Konkurrenz um eine Ressource, bei der ein Individuum den Kürzeren zieht.

Die daraus entstehenden Lotka-Volterra-Konkurrenzgleichungen der Theoretischen Biologie sind ein klassischer Ansatz zur Beschreibung der Dynamik einer stark vereinfachten Biozönose, bestehend aus einer nachwachsenden Ressource ${\displaystyle r}$ und mindestens 2 darum konkurrierender Arten:

${\displaystyle {\frac {\mathrm {d} N_{1}}{\mathrm {d} t}}=N_{1}ar-m_{1}N_{1}}$

${\displaystyle {\frac {\mathrm {d} N_{2}}{\mathrm {d} t}}=N_{2}br-m_{2}N_{2}}$,

wobei a,b exponentielle Wachstumsraten sind und die m (mortality rate) Sterberaten darstellen. Die zu einem Zeitpunkt verfügbare Ressourcenmenge wird angenommen als: ${\displaystyle r_{\mathrm {max} }-(pN_{1}+qN_{2})}$

Damit ergibt sich:

${\displaystyle {\frac {\mathrm {d} N_{1}}{\mathrm {d} t}}=N_{1}a(r_{\mathrm {max} }-(pN_{1}+qN_{2}))-m_{1}N_{1}}$

${\displaystyle {\frac {\mathrm {d} N_{2}}{\mathrm {d} t}}=N_{2}b(r_{\mathrm {max} }-(pN_{1}+qN_{2}))-m_{2}N_{2}}$

durch Ausmultiplizieren und Ersetzung der Koeffizienten

${\displaystyle \epsilon _{1}=ar_{\mathrm {max} }-m_{1}}$

${\displaystyle \gamma _{1}=aq}$

${\displaystyle \alpha _{1}=ap}$

${\displaystyle \epsilon _{2}=m_{2}-br_{\mathrm {max} }}$

${\displaystyle \gamma _{2}=-bp}$

${\displaystyle \alpha _{2}=-bq}$

kommt man zu Gleichungen der Form

${\displaystyle {\begin{aligned}{\frac {\mathrm {d} N_{1}}{\mathrm {d} t}}&=N_{1}(\epsilon _{1}-\gamma _{1}N_{2}-\alpha _{1}N_{1}),\\{\frac {\mathrm {d} N_{2}}{\mathrm {d} t}}&=N_{2}(-\epsilon _{2}+\gamma _{2}N_{1}-\alpha _{2}N_{2}),\end{aligned}}}$

die wiederum zwei Gleichgewichtslagen zulassen: den trivialen Gleichgewichtspunkt ${\displaystyle (0,0)}$, und den inneren Gleichgewichtspunkt ${\displaystyle (N_{1}^{*},N_{2}^{*})}$, der durch ein lineares Gleichungssystem gegeben ist:

${\displaystyle {\begin{pmatrix}\alpha _{1}&\gamma _{1}\\-\gamma _{2}&\alpha _{2}\end{pmatrix}}{\begin{pmatrix}N_{1}^{*}\\N_{2}^{*}\end{pmatrix}}={\begin{pmatrix}\epsilon _{1}\\-\epsilon _{2}\end{pmatrix}}}$

Durch Lösen dieses Gleichungssystems findet man den Gleichgewichtspunkt

${\displaystyle {\begin{pmatrix}N_{1}^{*}\\N_{2}^{*}\end{pmatrix}}={\frac {1}{\alpha _{1}\alpha _{2}+\gamma _{1}\gamma _{2}}}{\begin{pmatrix}\alpha _{2}\epsilon _{1}+\gamma _{1}\epsilon _{2}\\\gamma _{2}\epsilon _{1}-\alpha _{1}\epsilon _{2}\end{pmatrix}},}$

der unter der Bedingung ${\displaystyle \gamma _{2}\epsilon _{1}>\alpha _{1}\epsilon _{2}}$ im ersten Quadranten liegt.

Zu diesem erweiterten Lotka-Volterra-Systems gibt es auch eine Ljapunow-Funktion:

${\displaystyle V(N_{1},N_{2})=\gamma _{2}(N_{1}-N_{1}^{*}\ln N_{1})+\gamma _{1}(N_{2}-N_{2}^{*}\ln N_{2}),}$

mit der die Voraussetzungen des Zweiten Kriteriums von Ljapunow für den Gleichgewichtspunkt ${\displaystyle (N_{1}^{*},N_{2}^{*})}$ erfüllt sind. Daraus folgt, dass dieser Gleichgewichtspunkt jetzt asymptotisch stabil ist.

Mehr als zwei Klassen von Lebewesen

Ein großer Teil von Volterras Buch bezieht sich auf Erweiterungen seines Systems auf mehr als zwei Klassen von Lebewesen, die in unterschiedlichen Weisen miteinander interagieren.

Fischereidaten

In der Einleitung zu Volterras Buch findet sich eine Tabelle, die zu den Jahren 1905 und 1910–1923 und zu drei Fischereihäfen jeweils den prozentualen Anteil der Knorpelfische (Sélaciens), also insbesondere der Haie, am gesamten Fischfang des Fischereihafens enthält:

«Cela prouve pendant la période 1915–1920, où la pêche était moins intense à cause de la guerre, un accroissement relatif de la classe des Sélaciens qui, particulièrement voraces, se nourissent d’autres poissons. Les statistique inclinent donc à penser qu’une diminution dans l’intensité de la destruction favorise les espèces les plus voraces.»

Diese Statistiken zeigen in den Jahren 1915 bis 1920, als der Fischfang im Mittelmeer wegen des Ersten Weltkriegs weniger intensiv war, einen erhöhten Anteil an Raubfischen, der dann mit der Intensivierung der Fischerei nach 1920 wieder zurückgeht. Das dritte Lotka-Volterra-Gesetz, die Verschiebung der Mittelwerte, bietet hierfür eine plausible Erklärung.

Medizinische Epidemiologie

In der Theoretischen Biologie sowie in der medizinischen Epidemiologie finden Modelle vom Lotka-Volterra-Typ zur Beschreibung der Ausbreitungsprozesse von Krankheiten Verwendung. Einige Beispiele finden sich in SI-Modell, SIR-Modell und SIS-Modell.

Wirtschaftswissenschaften

Dem Goodwin-Modell zur Erklärung von Konjunkturschwankungen liegen Lotka-Volterra-Gleichungen zugrunde, wobei der Lohnquote die Rolle des Räubers und der Beschäftigungsquote die Rolle der Beute zukommt.

Gerold Blümle entwickelte ein Konjunkturmodell, in dem (mathematisch) der Investitionsquote die Rolle der Raubtiere zukommt, und der Streuung oder Varianz der Gewinne die Rolle der Beutetiere. Bei Frank Schohl kommt der Varianz der Renditenänderungen der Unternehmen die Rolle der Raubtiere, der Varianz der Angebotsänderungen der Unternehmen die Rolle der Beutetiere zu.