Logistisches Wachstum am Beispiel einer Infektionskrankheit Eine Infektionskrankheit (z.B. Corona) breitet sich in einem Land aus.
Zu Beginn des Infektionsgeschehens sind nur wenige Personen infiziert. Die Zahl der Infizierten B ( t ) B(t) B ( t ) wächst angenähert exponentiell (lila \textcolor{660099}{\text{lila} } lila Teil des Graphens). Je mehr Personen infiziert sind, desto mehr Personen werden neu infiziert.
(Ein rein exponentielles Wachstum zeigt die gr u ¨ ne \textcolor{006400}{\text{grüne}} gr u ¨ ne Kurve.)
Nach einer gewissen Zeit wird es immer schwerer, Nichtinfizierte (Gesunde) zu finden, die angesteckt werden können. Das Wachstum nimmt ab und wird zu einem begrenzten Wachstum , mit einer Schranke S S S ( roter (\textcolor{cc0000}{\text{roter} } ( roter Teil des Graphens in der Abbildung).
Funktionsgraph für logistisches Wachstum
▸ Herleitung der Differenzialgleichung für logistisches Wachstum
MerkeLogistische Wachstumsfunktion Die Berechnung des Bestandes B ( t ) B(t) B ( t ) der logistischen Wachstumsfunktion ist wie folgt gegeben:
B ( t ) = B 0 ⋅ S B 0 + ( S − B 0 ) ⋅ e − S ⋅ k ⋅ t \displaystyle B(t)=\frac{B_0\cdot S}{B_0+(S-B_0)\cdot e^{-S\cdot k\cdot t}} B ( t ) = B 0 + ( S − B 0 ) ⋅ e − S ⋅ k ⋅ t B 0 ⋅ S
Nach Kürzen mit B 0 B_0 B 0 erhält man:
B ( t ) = S 1 + ( S B 0 − 1 ) ⋅ e − S ⋅ k ⋅ t \displaystyle B(t)=\frac{S}{1+\left(\dfrac{S}{B_0}-1\right)\cdot e^{-S\cdot k\cdot t}} B ( t ) = 1 + ( B 0 S − 1 ) ⋅ e − S ⋅ k ⋅ t S Dabei gilt:
B ( t ) B(t) B ( t ) : ist der Bestand zur Zeit t t t ,
k k k : ist die Wachstumskonstante ,
S S S : ist die Sättigungsgrenze (Schranke),
B 0 B_0 B 0 : ist der Anfangsbestand zur Zeit t = 0 t=0 t = 0 , also der Startwert.
Beispiel Insektenpopulation In einem Labor werden Insekten in Glasbehältern gezüchtet. Anfangs besteht die Population aus 25 25 25 Insekten. Am nächsten Tag wurden 30 30 30 Insekten gezählt. Man schätzt, dass es aufgrund des beschränkten Platzes maximal 400 400 400 Insekten werden können.
a) Wie lautet die Wachstumsfunktion für die Zahl der Insekten?
b) Wie sieht der zugehörige Graph aus?
Lösung zu a) Gegeben sind der Anfangswert der Insektenanzahl B 0 = 25 B_0=25 B 0 = 25 und die Schranke S = 400 S=400 S = 400 . Die Werte werden in die Gleichung für die logistische Funktion eingesetzt:
B ( t ) \displaystyle B(t) B ( t ) = = = S 1 + ( S B 0 − 1 ) ⋅ e − S ⋅ k ⋅ t \displaystyle \dfrac{S}{1+\left(\dfrac{S}{B_0}-1\right)\cdot e^{-S\cdot k\cdot t}} 1 + ( B 0 S − 1 ) ⋅ e − S ⋅ k ⋅ t S ↓ Setze B 0 = 25 B_0=25 B 0 = 25 und S = 400 S=400 S = 400 ein.
= = = 400 1 + ( 400 25 − 1 ) ⋅ e − 400 ⋅ k ⋅ t \displaystyle \dfrac{400}{1+\left(\dfrac{400}{25}-1\right)\cdot e^{-400\cdot k\cdot t}} 1 + ( 25 400 − 1 ) ⋅ e − 400 ⋅ k ⋅ t 400 ↓ Vereinfache.
= = = 400 1 + ( 16 − 1 ) ⋅ e − 400 ⋅ k ⋅ t \displaystyle \dfrac{400}{1+\left(16-1\right)\cdot e^{-400\cdot k\cdot t}} 1 + ( 16 − 1 ) ⋅ e − 400 ⋅ k ⋅ t 400 ↓ Vereinfache.
= = = 400 1 + 15 ⋅ e − 400 ⋅ k ⋅ t \displaystyle \dfrac{400}{1+15\cdot e^{-400\cdot k\cdot t}} 1 + 15 ⋅ e − 400 ⋅ k ⋅ t 400
Bisher gilt: B ( t ) = 400 1 + 15 ⋅ e − 400 ⋅ k ⋅ t B(t)=\dfrac{400}{1+15\cdot e^{-400\cdot k\cdot t}} B ( t ) = 1 + 15 ⋅ e − 400 ⋅ k ⋅ t 400
Die unbekannte Wachstumskonstante k k k kann über die Aussage, dass nach 1 1 1 Tag 30 30 30 Insekten vorhanden sind, berechnet werden.
Setze B ( 1 ) = 30 B(1)=30 B ( 1 ) = 30 :
B ( t ) \displaystyle B\left(t\right) B ( t ) = = = 400 1 + 15 ⋅ e − 400 ⋅ k ⋅ t \displaystyle \dfrac{400}{1+15\cdot e^{-400\cdot k\cdot t}} 1 + 15 ⋅ e − 400 ⋅ k ⋅ t 400 ↓ Setze t = 1 t=1 t = 1 und B ( 1 ) = 30 B(1)=30 B ( 1 ) = 30 ein.
30 \displaystyle 30 30 = = = 400 1 + 15 ⋅ e − 400 ⋅ k ⋅ 1 \displaystyle \dfrac{400}{1+15\cdot e^{-400\cdot k\cdot 1}} 1 + 15 ⋅ e − 400 ⋅ k ⋅ 1 400 ⋅ ( 1 + 15 ⋅ e − 400 ⋅ k ⋅ 1 ) \displaystyle \cdot(1+15\cdot e^{-400\cdot k\cdot 1}) ⋅ ( 1 + 15 ⋅ e − 400 ⋅ k ⋅ 1 ) ↓ Löse nach k k k auf.
30 ⋅ ( 1 + 15 ⋅ e − 400 ⋅ k ⋅ 1 ) \displaystyle 30\cdot(1+15\cdot e^{-400\cdot k\cdot 1}) 30 ⋅ ( 1 + 15 ⋅ e − 400 ⋅ k ⋅ 1 ) = = = 400 \displaystyle 400 400 : 30 \displaystyle :30 : 30 1 + 15 ⋅ e − 400 ⋅ k \displaystyle 1+15\cdot e^{-400\cdot k} 1 + 15 ⋅ e − 400 ⋅ k = = = 400 30 \displaystyle \dfrac{400}{30} 30 400 ↓ Kürze den Bruch.
1 + 15 ⋅ e − 400 ⋅ k \displaystyle 1+15\cdot e^{-400\cdot k} 1 + 15 ⋅ e − 400 ⋅ k = = = 40 3 \displaystyle \dfrac{40}{3} 3 40 − 1 \displaystyle -1 − 1 15 ⋅ e − 400 ⋅ k \displaystyle 15\cdot e^{-400\cdot k} 15 ⋅ e − 400 ⋅ k = = = 40 3 − 1 \displaystyle \dfrac{40}{3}-1 3 40 − 1 ↓ Vereinfache die rechte Seite.
15 ⋅ e − 400 ⋅ k \displaystyle 15\cdot e^{-400\cdot k} 15 ⋅ e − 400 ⋅ k = = = 40 3 − 3 3 \displaystyle \dfrac{40}{3}-\dfrac{3}{3} 3 40 − 3 3 ↓ Fasse zusammen.
15 ⋅ e − 400 ⋅ k \displaystyle 15\cdot e^{-400\cdot k} 15 ⋅ e − 400 ⋅ k = = = 37 3 \displaystyle \dfrac{37}{3} 3 37 : 15 \displaystyle :15 : 15 e − 400 ⋅ k \displaystyle e^{-400\cdot k} e − 400 ⋅ k = = = 37 45 \displaystyle \dfrac{37}{45} 45 37 ln \displaystyle \ln ln ↓ Löse die e-Funktion auf.
− 400 ⋅ k \displaystyle -400\cdot k − 400 ⋅ k = = = ln ( 37 45 ) \displaystyle \ln\left(\dfrac{37}{45}\right) ln ( 45 37 ) : ( − 400 ) \displaystyle :(-400) : ( − 400 ) k \displaystyle k k = = = ln ( 37 45 ) ( − 400 ) \displaystyle \dfrac{\ln\left(\dfrac{37}{45}\right)}{(-400)} ( − 400 ) ln ( 45 37 ) k \displaystyle k k ≈ ≈ ≈ 0 , 00049 \displaystyle 0{,}00049 0 , 00049
Somit folgt: B ( t ) = 400 1 + 15 ⋅ e − 400 ⋅ 0 , 00049 ⋅ t = 400 1 + 15 ⋅ e − 0 , 196 ⋅ t B(t)=\dfrac{400}{1+15\cdot e^{-400\cdot0{,}00049\cdot t}}=\dfrac{400}{1+15\cdot e^{-0{,}196\cdot t}} B ( t ) = 1 + 15 ⋅ e − 400 ⋅ 0 , 00049 ⋅ t 400 = 1 + 15 ⋅ e − 0 , 196 ⋅ t 400
Antwort: Die Wachstumsfunktion für die Zahl der Insekten lautet:
B ( t ) = 400 1 + 15 ⋅ e − 0 , 196 ⋅ t \displaystyle B(t)=\dfrac{400}{1+15\cdot e^{-0{,}196\cdot t}} B ( t ) = 1 + 15 ⋅ e − 0 , 196 ⋅ t 400 Lösung zu b) Für die Zeichnung des Graphen ist es sinnvoll, einige Funktionswerte zu berechnen:
B(t): Anzahl Insekten
(gerundete Werte)
Graphische Darstellung der logistischen Wachstumsfunktion "Zahl der Insekten"
Weitere Eigenschaften der logistischen Funktion 1. Schnittpunkt mit der y y y -Achse: Zur Bestimmung des Schnittpunktes mit der y y y -Achse setze t = 0 t=0 t = 0 :
B ( t ) \displaystyle B(t) B ( t ) = = = S 1 + ( S B 0 − 1 ) ⋅ e − S ⋅ k ⋅ t \displaystyle \dfrac{S}{1+\left(\dfrac{S}{B_0}-1\right)\cdot e^{-S\cdot k\cdot t}} 1 + ( B 0 S − 1 ) ⋅ e − S ⋅ k ⋅ t S ↓ Setze t = 0 t=0 t = 0 ein.
B ( 0 ) \displaystyle B(0) B ( 0 ) = = = S 1 + ( S B 0 − 1 ) ⋅ e 0 \displaystyle \dfrac{S}{1+\left(\dfrac{S}{B_0}-1\right)\cdot e^0} 1 + ( B 0 S − 1 ) ⋅ e 0 S ↓ e 0 = 1 e^0=1 e 0 = 1
B ( 0 ) \displaystyle B(0) B ( 0 ) = = = S 1 + S B 0 − 1 \displaystyle \dfrac{S}{1+\dfrac{S}{B_0}-1} 1 + B 0 S − 1 S ↓ Vereinfache den Nenner.
B ( 0 ) \displaystyle B(0) B ( 0 ) = = = S S B 0 \displaystyle \dfrac{S}{\dfrac{S}{B_0}} B 0 S S ↓ Vereinfache.
B ( 0 ) \displaystyle B(0) B ( 0 ) = = = B 0 \displaystyle B_0 B 0
2. Asymptotisches Verhalten: Zur Bestimmung des asymptotischen Verhaltens der logistischen Wachstumsfunktion berechne den Grenzwert lim t → ∞ B ( t ) \displaystyle\lim_{t\rightarrow\infty} B(t) t → ∞ lim B ( t ) :
B ( t ) = S 1 + ( S B 0 − 1 ) ⋅ e − S ⋅ k ⋅ t B(t)=\dfrac{S}{1+\left(\dfrac{S}{B_0}-1\right)\cdot e^{-S\cdot k\cdot t}} B ( t ) = 1 + ( B 0 S − 1 ) ⋅ e − S ⋅ k ⋅ t S
Betrachte zunächst nur den Nenner des Bruches:
1 + ( S B 0 − 1 ) ⋅ e − S ⋅ k ⋅ t 1+\left(\dfrac{S}{B_0}-1\right)\cdot e^{-S\cdot k\cdot t} 1 + ( B 0 S − 1 ) ⋅ e − S ⋅ k ⋅ t
Berechne den Grenzwert für t → ∞ t\rightarrow\infty t → ∞ :
lim t → ∞ 1 + ( S B 0 − 1 ) ⋅ e − S ⋅ k ⋅ t ⏟ → 0 = 1 + ( S B 0 − 1 ) ⋅ 0 = 1 \displaystyle\lim_{t\rightarrow\infty}1+\left(\frac{S}{B_0}-1\right)\cdot\underbrace{ e^{-S\cdot k\cdot t}}_{\rightarrow 0}=1+\left(\frac{S}{B_0}-1\right)\cdot 0=1
t → ∞ lim 1 + ( B 0 S − 1 ) ⋅ → 0 e − S ⋅ k ⋅ t = 1 + ( B 0 S − 1 ) ⋅ 0 = 1
▸ Warum geht die obige e-Funktion gegen null?
Der ganze Bruch S 1 + ( S B 0 − 1 ) ⋅ e − S ⋅ k ⋅ t \dfrac{S}{1+\left(\dfrac{S}{B_0}-1\right)\cdot e^{-S\cdot k\cdot t}} 1 + ( B 0 S − 1 ) ⋅ e − S ⋅ k ⋅ t S geht gegen S S S , denn lim t → ∞ S 1 = S \displaystyle\lim_{t\rightarrow\infty}\dfrac{S}{1}=S
t → ∞ lim 1 S = S
Der Funktionsgraph nähert sich asymptotisch der Schranke S S S .
3. Wendepunkt Für die Berechnung des Wendepunktes müssen folgende Bedingungen erfüllt sein:
B ′ ′ ( t w ) = 0 B''(t_w)=0 B ′′ ( t w ) = 0 und Vorzeichenwechsel von B ′ ′ ( t ) B''(t) B ′′ ( t ) an der Stelle t w t_w t w
Die 2. 2. 2. Ableitung von B ( t ) B(t) B ( t ) muss berechnet werden.
Die 1. 1. 1. Ableitung ist bekannt: B ′ ( t ) = k ⋅ B ( t ) ⋅ ( S − B ( t ) ) B'(t)=k\cdot B(t)\cdot (S-B(t)) B ′ ( t ) = k ⋅ B ( t ) ⋅ ( S − B ( t ))
Mit der Produktregel wird B ′ ( t ) B'(t) B ′ ( t ) abgeleitet:
B ′ ′ ( t ) \displaystyle B''(t) B ′′ ( t ) = = = k ⋅ ( B ′ ( t ) ⋅ ( S − B ( t ) ) + B ( t ) ⋅ ( − B ′ ( t ) ) ) \displaystyle k\cdot\Big (B'(t)\cdot(S-B(t))+B(t)\cdot(-B'(t))\Big) k ⋅ ( B ′ ( t ) ⋅ ( S − B ( t )) + B ( t ) ⋅ ( − B ′ ( t )) ) ↓ Vereinfache.
B ′ ′ ( t ) \displaystyle B''(t) B ′′ ( t ) = = = k ⋅ ( B ′ ( t ) ⋅ S − B ′ ( t ) ⋅ B ( t ) − B ′ ( t ) ⋅ B ( t ) ) \displaystyle k\cdot\Big(B'(t)\cdot S-B'(t)\cdot B(t)-B'(t)\cdot B(t)\Big) k ⋅ ( B ′ ( t ) ⋅ S − B ′ ( t ) ⋅ B ( t ) − B ′ ( t ) ⋅ B ( t ) ) ↓ Fasse zusammen.
B ′ ′ ( t ) \displaystyle B''(t) B ′′ ( t ) = = = k ⋅ ( B ′ ( t ) ⋅ S − 2 ⋅ B ′ ( t ) ⋅ B ( t ) ) \displaystyle k\cdot\Big(B'(t)\cdot S-2\cdot B'(t)\cdot B(t)\Big) k ⋅ ( B ′ ( t ) ⋅ S − 2 ⋅ B ′ ( t ) ⋅ B ( t ) )
Beim Wendepunkt muss B ′ ′ ( t W ) = 0 B''(t_W)=0 B ′′ ( t W ) = 0 sein:
0 \displaystyle 0 0 = = = k ⋅ ( B ′ ( t W ) ⋅ S − 2 ⋅ B ′ ( t W ) ⋅ B ( t W ) ) \displaystyle k\cdot\Big(B'(t_W)\cdot S-2\cdot B'(t_W)\cdot B(t_W)\Big) k ⋅ ( B ′ ( t W ) ⋅ S − 2 ⋅ B ′ ( t W ) ⋅ B ( t W ) ) ↓ Wegen k ≠ 0 k\neq 0 k = 0 muss die Klammer null werden.
0 \displaystyle 0 0 = = = B ′ ( t W ) ⋅ S − 2 ⋅ B ′ ( t W ) ⋅ B ( t W ) \displaystyle B'(t_W)\cdot S-2\cdot B'(t_W)\cdot B(t_W) B ′ ( t W ) ⋅ S − 2 ⋅ B ′ ( t W ) ⋅ B ( t W ) ↓ Klammere B ′ ( t W ) B'(t_W) B ′ ( t W ) aus.
0 \displaystyle 0 0 = = = B ′ ( t W ) ⋅ ( S − 2 ⋅ B ( t W ) ) \displaystyle B'(t_W)\cdot (S-2\cdot B(t_W)) B ′ ( t W ) ⋅ ( S − 2 ⋅ B ( t W )) ↓ Die logistische Funktion ist streng monoton steigend .
Daher ist B ′ ( t ) > 0 B'(t) >0 B ′ ( t ) > 0 und die Klammer muss null werden.
0 \displaystyle 0 0 = = = S − 2 ⋅ B ( t W ) \displaystyle S-2\cdot B(t_W) S − 2 ⋅ B ( t W ) + 2 ⋅ B ( t W ) \displaystyle +2\cdot B(t_W) + 2 ⋅ B ( t W ) ↓ Löse nach B ( t W ) B(t_W) B ( t W ) auf.
2 ⋅ B ( t W ) \displaystyle 2\cdot B(t_W) 2 ⋅ B ( t W ) = = = S \displaystyle S S : 2 \displaystyle :2 : 2 B ( t W ) \displaystyle B(t_W) B ( t W ) = = = S 2 \displaystyle \dfrac{S}{2} 2 S
Nun muss der Vorzeichenwechsel von B ′ ′ ( t ) B''(t) B ′′ ( t ) an der Stelle t w t_w t w überprüft werden.
B ′ ′ ( t ) = B ′ ( t ) ⋅ ( S − 2 ⋅ B ( t ) ) B''(t)=B'(t)\cdot (S-2\cdot B(t)) B ′′ ( t ) = B ′ ( t ) ⋅ ( S − 2 ⋅ B ( t ))
Die logistische Funktion ist streng monoton steigend .
Daher ist B ′ ( t ) > 0 B'(t) >0 B ′ ( t ) > 0 und die Klammer ( S − 2 ⋅ B ( t ) ) (S-2\cdot B(t)) ( S − 2 ⋅ B ( t )) bestimmt das Vorzeichen von B ′ ′ ( t ) B''(t) B ′′ ( t ) .
Es gilt:
für t < t w t<t_w t < t w ist B ( t ) < S 2 B(t)<\dfrac{S}{2} B ( t ) < 2 S ⇒ S − 2 ⋅ B ( t ) > 0 \;\Rightarrow\;S-2\cdot B(t)>0 ⇒ S − 2 ⋅ B ( t ) > 0
und
für t > t w t>t_w t > t w ist B ( t ) > S 2 B(t)>\dfrac{S}{2} B ( t ) > 2 S ⇒ S − 2 ⋅ B ( t ) < 0 \;\Rightarrow\;S-2\cdot B(t)<0 ⇒ S − 2 ⋅ B ( t ) < 0
B ′ ′ ( t ) B''(t) B ′′ ( t ) hat an der Stelle t w t_w t w einen Vorzeichenwechsel.
Die beiden Bedingungen sind somit erfüllt.
4. Maximale Wachstumsgeschwindigkeit An der Stelle der größten Steigung ist die Wachstumsgeschwindigkeit maximal.
Die größte Steigung tritt im Wendepunkt auf, d.h. v max = B ′ ( t W ) v_{\max}=B'(t_W) v m a x = B ′ ( t W ) .
Setze in die erste Ableitung von B ( t ) B(t) B ( t ) die Koordinaten des Wendepunktes W P ( t W ∣ S 2 ) WP\left(t_W\Big\vert\dfrac{S}{2}\right) W P ( t W 2 S ) ein:
v max \displaystyle v_{\max} v m a x = = = B ′ ( t W ) \displaystyle B'(t_W) B ′ ( t W ) ↓ Setze die erste Ableitung von B ( t ) B(t) B ( t ) ein.
= = = k ⋅ B ( t W ) ⋅ ( S − B ( t W ) ) \displaystyle k\cdot B(t_W)\cdot (S-B(t_W)) k ⋅ B ( t W ) ⋅ ( S − B ( t W )) ↓ Zur Zeit t W t_W t W ist B ( t W ) = S 2 B(t_W)=\dfrac{S}{2} B ( t W ) = 2 S . Setze für B ( t w ) B(t_w) B ( t w ) den Wert S 2 \dfrac{S}{2} 2 S ein.
= = = k ⋅ S 2 ⋅ ( S − S 2 ) \displaystyle k\cdot\dfrac{S}{2}\cdot\left(S-\dfrac{S}{2}\right) k ⋅ 2 S ⋅ ( S − 2 S ) ↓ Berechne die Klammer.
= = = k ⋅ S 2 ⋅ S 2 \displaystyle k\cdot \dfrac{S}{2}\cdot \dfrac{S}{2} k ⋅ 2 S ⋅ 2 S ↓ Fasse zusammen.
= = = k ⋅ S 2 4 \displaystyle k\cdot \dfrac{S^2}{4} k ⋅ 4 S 2
▸ Herleitung der Funktionsgleichung durch Integration
Übungsaufgaben: Logistisches Wachstum Weitere Aufgaben zum Thema findest du im folgenden Aufgabenordner:Aufgaben zu Wachstums- und Zerfallsprozessen