Logistische Funktion

Logistische Funktion für den Fall ${\displaystyle G=1\,;\,k=1\,;\,f(0)=0{,}5}$

Die logistische Funktion charakterisiert eine stetige eindimensionale Wahrscheinlichkeitsverteilung (die logistische Verteilung) und ist eine funktionelle Darstellung von Sättigungsprozessen aus der Klasse der sogenannten Sigmoidfunktionen mit unbegrenzter zeitlicher Ausdehnung.

Der Graph der Funktion beschreibt eine S-förmige Kurve, ein Sigmoid. Heute ist der Name logistische Kurve eindeutig der S-Funktion zugeordnet, wohingegen noch bis ins 20. Jahrhundert gelegentlich auch der Logarithmus mit dem italienischen Namen der logistischen Kurve (curva logistica) belegt wurde.

Die Funktion wird manchmal auch mit Expit bezeichnet, da die Umkehrfunktion der logistischen Funktion die Logit-Funktion ist.

Beschreibung

Die logistische Funktion, wie sie sich aus der diskreten logistischen Gleichung ergibt, beschreibt den Zusammenhang zwischen der verstreichenden Zeit und einem Wachstum. Hierzu wird das Modell des exponentiellen Wachstums modifiziert durch eine sich mit dem Wachstum verbrauchende Ressource, die eine obere Schranke darstellt.

Zur Anfangszeit ist der Funktionswert nicht 0, sondern es gilt ${\displaystyle f(0)>0}$.

Es gilt:

• Die obere Schranke ${\displaystyle G}$ bildet eine Grenze für den Funktionswert ${\displaystyle f(t)}$.
• Das Wachstum ${\displaystyle f'(t)}$ ist proportional zu:
  • dem aktuellen Bestand ${\displaystyle f(t)}$,
  • der noch vorhandenen Kapazität ${\displaystyle G-f(t)}$
  • und einer Wachstumskonstanten ${\displaystyle k}$.

Diese Entwicklung wird daher durch eine Bernoullische Differentialgleichung der Form

${\displaystyle f'(t)=k\cdot f(t)\cdot \left(G-f(t)\right)}$

mit einer Proportionalitätskonstanten ${\displaystyle k}$ beschrieben. Das Lösen dieser Differentialgleichung ergibt:

${\displaystyle f(t)=G\cdot {\frac {1}{1+\mathrm {e} ^{-k\cdot G\cdot t}\left({\frac {G}{f(0)}}-1\right)}}.}$

Am Anfang ist das Wachstum klein, da die Population und somit die Zahl der sich vermehrenden Individuen gering ist. In der Mitte der Entwicklung (genauer: im Wendepunkt) wächst die Population am stärksten, bis sie durch die sich erschöpfenden Ressourcen gebremst wird.

Anwendungen

Beispiel einer Epidemie: Krankheits- und Todesfälle (schwarz) im Verlauf der Ebolafieber-Epidemie in Westafrika bis Juli 2014 (annähernd logistische Funktionen)

Die logistische Gleichung beschreibt einen sehr häufig auftretenden Zusammenhang, wie der Beschreibung einer Population von Lebewesen, beispielsweise einer idealen Bakterien­population, die auf einem Bakteriennährboden begrenzter Größe wächst. Ein weiteres Beispiel ist (annähernd) die Verbreitung einer Infektionskrankheit mit anschließender permanenter Immunität, bei der mit der Zeit eine abnehmende Anzahl für die Infektionskrankheit anfällige Individuen übrig bleiben. Eine Anwendung findet die logistische Funktion auch im SI-Modell der mathematischen Epidemiologie.

Die Funktion findet weit über den Modellen aus der Biologie hinaus Anwendung. Auch der Lebenszyklus eines Produktes im Markt kann mit der logistischen Funktion nachgebildet werden. Weitere Anwendungsbereiche sind Wachstums- und Zerfallsprozesse in der Sprache (Sprachwandelgesetz, Piotrowski-Gesetz) sowie die Entwicklung im Erwerb der Muttersprache (Spracherwerbsgesetz).

Lösung der Differentialgleichung

Sei ${\displaystyle F\colon \mathbb {R} \to \mathbb {R} ,t\mapsto kt(G-t)}$. ${\displaystyle F}$ ist stetig. Es gilt, die Differentialgleichung

${\displaystyle {\frac {\mathrm {d} f}{\mathrm {d} t}}(t)=F(f(t))}$

zu lösen.

Die Differentialgleichung lässt sich mit dem Verfahren „Trennung der Variablen“ lösen. Es gilt ${\displaystyle F(t)\neq 0}$ für alle ${\displaystyle t\in \mathbb {R} }$, also ist die Abbildung

${\displaystyle \Phi (t):=\int _{0}^{t}{\frac {1}{F(s)}}\,{\mathrm {d} }s}$

auf den reellen Zahlen wohldefiniert.

Nach der Trennung der Variablen ist die Lösung ${\displaystyle f}$ der obigen Differentialgleichung also identisch mit der Lösung der Differentialgleichung

${\displaystyle \int _{f(0)}^{f(t)}{\frac {1}{F(s)}}\,\mathrm {d} s=t.}$

Durch Partialbruchzerlegung ergibt sich

${\displaystyle \int _{f(0)}^{f(t)}{\frac {1}{F(s)}}\,\mathrm {d} s=\int _{f(0)}^{f(t)}{\frac {1}{ks(G-s)}}\,\mathrm {d} s={\frac {1}{kG}}\int _{f(0)}^{f(t)}\left({\frac {1}{s}}+{\frac {1}{G-s}}\right)\,\mathrm {d} s.}$

Nach dem Hauptsatz der Differential- und Integralrechnung ist das obige Integral

${\displaystyle =\left.{\frac {1}{kG}}{\bigg (}\ln(s)-\ln(G-s){\bigg )}\right|_{s=f(0)}^{s=f(t)}={\frac {1}{kG}}\left(\ln \left({\frac {f(t)}{G-f(t)}}\right)+c\right)}$

wobei

${\displaystyle c:=-\ln {\frac {f(0)}{G-f(0)}}=\ln \left({\frac {G}{f(0)}}-1\right).}$

Es gilt also, die Funktionsgleichung

${\displaystyle \ln {\frac {f(t)}{G-f(t)}}=kGt-c}$

zu lösen, solange die ${\displaystyle f(t)}$ zwischen ${\displaystyle 0}$ und ${\displaystyle G}$ liegen, was wegen der Voraussetzung ${\displaystyle 0<f(0)<G}$ angenommen werden kann. Dabei ist ${\displaystyle \ln }$ der natürliche Logarithmus. Die Anwendung der Exponentialfunktion auf beiden Seiten führt zu

${\displaystyle \mathrm {e} ^{kGt-c}\,=\,{\frac {f(t)}{G-f(t)}}}$

und anschließende Kehrwertbildung zu

${\displaystyle (*)\quad \quad \quad \mathrm {e} ^{-kGt+c}\,=\,{\frac {G-f(t)}{f(t)}}\,=\,{\frac {G}{f(t)}}-1.}$

Wir bringen nun die ${\displaystyle 1}$ auf die linke Seite, bilden dann erneut den Kehrwert, und erhalten schließlich

${\displaystyle {\frac {f(t)}{G}}\,=\,{\frac {1}{1+\mathrm {e} ^{-kGt+c}}}}$

und daraus

${\displaystyle (**)\quad \quad \quad f(t)\,=\,G\cdot {\frac {1}{1+\mathrm {e} ^{-kGt+c}}}.}$

Setzen wir die Definition von ${\displaystyle c}$ in die gefundene Lösung (**) ein, so kommen wir zur oben behaupteten Lösung der logistischen Differentialgleichung:

${\displaystyle f(t)\,=\,G\cdot {\frac {1}{1+\mathrm {e} ^{-kGt+c}}}\,=\,G\cdot {\frac {1}{1+\mathrm {e} ^{-kGt}\mathrm {e} ^{c}}}\,=\,G\cdot {\frac {1}{1+\mathrm {e} ^{-kGt}\left({\frac {G}{f(0)}}-1\right)}}.}$

An dieser Funktionsgleichung liest man leicht ab, dass die Werte immer zwischen ${\displaystyle 0}$ und ${\displaystyle G}$ liegen, weshalb die Lösung für alle ${\displaystyle -\infty <t<\infty }$ gilt. Das kann man im Nachhinein natürlich auch durch Einsetzen in die Differentialgleichung bestätigen.

Berechnung des Wendepunkts

Zur Bestimmung des Wendepunktes der Lösungsfunktion ${\displaystyle f}$ bestimmen wir zunächst mittels Produktregel die Ableitungen

${\displaystyle {\begin{aligned}f'(t)&=k\cdot f(t)\cdot (G-f(t))\\f''(t)&=k\cdot f'(t)\cdot (G-f(t))+k\cdot f(t)\cdot (-f'(t))\\&=k\cdot f'(t)\cdot (G-f(t)-f(t))=k\cdot f'(t)\cdot (G-2\cdot f(t))\end{aligned}}}$

und bestimmen die Nullstelle ${\displaystyle t_{\mathrm {W} }}$ der zweiten Ableitung:

${\displaystyle f''(t_{\mathrm {W} })=k\cdot f'(t_{\mathrm {W} })\cdot (G-2\cdot f(t_{\mathrm {W} }))=0}$

${\displaystyle G-2\cdot f(t_{\mathrm {W} })=0}$

${\displaystyle G=2\cdot f(t_{\mathrm {W} })}$

${\displaystyle f(t_{\mathrm {W} })={\tfrac {G}{2}}.}$

Damit kennen wir den Funktionswert im Wendepunkt und stellen fest, dass die Population im Wendepunkt gerade die halbe Sättigungsgrenze überschreitet. Zur Bestimmung von ${\displaystyle t_{\mathrm {W} }}$ verwenden wir für ${\displaystyle f(t_{\mathrm {W} })={\tfrac {G}{2}}}$ die Lösungsformel und rechnen wie folgt:

${\displaystyle G\cdot {\frac {1}{1+\mathrm {e} ^{-k\cdot G\cdot t_{\mathrm {W} }}\cdot \mathrm {e} ^{c}}}={\frac {G}{2}}}$

${\displaystyle 1+\mathrm {e} ^{-k\cdot G\cdot t_{\mathrm {W} }}\cdot \mathrm {e} ^{c}=2}$

${\displaystyle \mathrm {e} ^{-k\cdot G\cdot t_{\mathrm {W} }}\cdot \mathrm {e} ^{c}=1=\mathrm {e} ^{0}}$

${\displaystyle -k\cdot G\cdot t_{\mathrm {W} }+c=0}$

${\displaystyle t_{\mathrm {W} }={\frac {c}{k\cdot G}}.}$

Für ${\displaystyle \textstyle G>f(0)}$ folgt mit ${\displaystyle \textstyle \mathrm {e} ^{c}={\frac {G}{f(0)}}-1\Rightarrow c=\ln \left({\frac {G}{f(0)}}-1\right)}$ weiter:

${\displaystyle t_{\mathrm {W} }={\frac {\ln \left({\frac {G}{f(0)}}-1\right)}{k\cdot G}}.}$

Damit ist der Wendepunkt vollständig bestimmt und es gibt nur diesen einen. Durch Einsetzen von ${\displaystyle f(t_{\mathrm {W} })={\tfrac {G}{2}}}$ in die erste Ableitung erhält man die maximale Wachstumsgeschwindigkeit:

${\displaystyle f'(t_{\mathrm {W} })=k\cdot {\frac {G}{2}}\cdot \left(G-{\frac {G}{2}}\right)=k\cdot {\frac {G}{2}}\cdot {\frac {G}{2}}}$

${\displaystyle f'(t_{\mathrm {W} })={\frac {k\cdot G^{2}}{4}}.}$

Weitere Darstellungen

Aus

${\displaystyle f(t)=G\cdot {\frac {1}{1+\mathrm {e} ^{-k\cdot G\cdot t}\left({\frac {G}{f(0)}}-1\right)}}}$

folgt:

${\displaystyle f(t)=G\cdot {\frac {1}{1+\mathrm {e} ^{-k\cdot G\cdot t}\cdot {\frac {G}{f(0)}}-\mathrm {e} ^{-k\cdot G\cdot t}}}\cdot {\frac {f(0)}{f(0)}}={\frac {G\cdot f(0)}{f(0)+\mathrm {e} ^{-k\cdot G\cdot t}\cdot G-\mathrm {e} ^{-k\cdot G\cdot t}\cdot f(0)}}}$

${\displaystyle ={\frac {G\cdot f(0)}{f(0)+\left(G-f(0)\right)\cdot \mathrm {e} ^{-k\cdot G\cdot t}}}}$

oder auch:

${\displaystyle f(t)={\frac {G}{2}}\cdot \left(\tanh \left({\frac {kG}{2}}(t-t_{\mathrm {W} })\right)+1\right)}$,

wobei ${\displaystyle t_{\mathrm {W} }}$ die oben berechnete Wendestelle ist:

${\displaystyle t_{\mathrm {W} }={\frac {\ln \left({\frac {G}{f(0)}}-1\right)}{k\cdot G}}.}$

Siehe auch

• Logistische Regression
• SI-Modell

Literatur

• Volker Oppitz, Volker Nollau: Taschenbuch Wirtschaftlichkeitsrechnung. Quantitative Methoden der ökonomischen Analyse. Carl Hanser Verlag, München u. a. 2004, ISBN 3-446-22463-7.
• Volker Oppitz: Gabler-Lexikon Wirtschaftlichkeitsrechnung. Mit Anwendersoftware für Praxis und Studium. Gabler-Verlag Wiesbaden 1995, ISBN 3-409-19951-9
• Nicholas F. Britton: Essential Mathematical Biology. 3. printing. Springer, London u. a. 2005, ISBN 1-85233-536-X, (Springer undergraduate mathematics series).
• Peter Schönfeld: Methoden der Ökonometrie. 2 Bände. Vahlen, Berlin u. a. 1969–1971.
• Norman R. Draper, Harry Smith: Applied Regression Analysis. 3rd Edition. Wiley-Interscience, New York NY u. a. 1998, ISBN 0-471-17082-8, (Wiley Series in Probability and Statistics. Texts and References Section).