Logistisk vekst

Når veksten møter en grense. S-kurven, bæreevnen og vendepunktet, med det klassiske gjærforsøket fra 1913.

R1
Utforsk

Start her: utforsk

Slipp noen gjærceller ned i en kolbe med sukkerløsning. Først er det rikelig med mat og plass. Men hva skjer etter hvert, når det blir trangt og maten tar slutt?

Punktene er ekte målinger av gjærmengden i en kultur, time for time. Dataene ble publisert av T. Carlson i 1913. Se først på punktene, og prøv så modellene.

5101520−5050100150200250300350400450500550600650700750
Modell (tt = time): f(t)=28,47⋅e0,2239tf(t) = 28{,}47 \cdot e^{0{,}2239t}
Den deriverte: f′(t)=6,373⋅e0,2239tf'(t) = 6{,}373 \cdot e^{0{,}2239t}
I time 9: modellen gir 213,5, og vekstfarten er 47,8 per time.
Forklaringsgrad R2=−0,3421R^2 = −0{,}3421. Jo nærmere 1, jo bedre passer modellen til dataene, men sjekk alltid om modellen gir mening utenfor dataene.

Lim gjerne inn egne målinger.

Prøv dette:

  1. Beskriv formen på punktene med egne ord. Hva skjer i starten, i midten og på slutten?
  2. Når vokste gjæren raskest? Hvor mye økte den da per time?
  3. Velg eksponentiell modell. Når begynner den å bomme? Hvorfor?
  4. Velg logistisk modell. Hvilket tall nærmer gjærmengden seg?
Forklaring

Forklaring

Den logistiske modellen

Gjæren vokste nesten eksponentielt i starten. Etter hvert ble det trangt, maten tok slutt og avfallsstoffene hopet seg opp. Veksten bremset, og til slutt flatet mengden ut rundt 660. Grafen blir formet som en S: sakte start, rask vekst i midten og metning til slutt. Dette kalles logistisk vekst, og det beskriver alt fra bakterier og dyrebestander til spredning av sykdommer, nye teknologier og rykter.

N(t)=K1+ae−btN(t) = \frac{K}{1 + a e^{-bt}}

Vendepunktet

Logistisk vekstLogistisk vekstVekst som starter nesten eksponentielt, men flater ut mot en øvre grense, bæreevnen. Grafen er S-formet. har ett vendepunktVendepunktEt punkt der grafen skifter krumning, fra hul side opp til hul side ned eller omvendt. Der skifter f′′(x)f''(x) fortegn., og det ligger alltid midt mellom 0 og K:

N=K2na˚rt=ln⁡abN = \frac{K}{2} \quad \text{når} \quad t = \frac{\ln a}{b}

I vendepunktet er vekstfarten størst:

Nmaks′=K⋅b4N'_{\text{maks}} = \frac{K \cdot b}{4}

Før vendepunktet vokser veksten (N′′>0N'' > 0), etter vendepunktet avtar den (N′′<0N'' < 0).

Modell: Form S-kurven

Endre KK, aa og bb. Den oransje grafen er vekstfarten N′(t)N'(t). Legg merke til at toppen av den oransje grafen alltid ligger rett under vendepunktet.

246810121416182020406080100120140160180200220vendepunkt
N(t)N(t)N′(t)N'(t)
N(0)=4,8N(0) = 4{,}8, bæreevne K=150K = 150
Vendepunkt ved t=5,67t = 5{,}67, største vekstfart Kb4=22,5\frac{Kb}{4} = 22{,}5

Sammenligning med eksponentiell vekst

Helt i starten, når NN er liten sammenlignet med KK, ligner logistisk vekst på eksponentiell vekst. Forskjellen blir tydelig først når bestanden nærmer seg bæreevnen.

Eksempler

Gjennomgåtte eksempler

Eksempel 1: Les av modellen

En bestand er N(t)=10001+49e−0,5tN(t) = \dfrac{1000}{1 + 49e^{-0{,}5t}}, der tt er år. Finn startverdien, bæreevnen, vendepunktet og den største vekstfarten.

Startverdi: N(0)=10001+49=20N(0) = \dfrac{1000}{1 + 49} = 20.
Bæreevne: K=1000K = 1000.
Vendepunkt: t=ln⁡490,5≈7,78t = \dfrac{\ln 49}{0{,}5} \approx 7{,}78 år, da er N=500N = 500.
Største vekstfart: 1000⋅0,54=125\dfrac{1000 \cdot 0{,}5}{4} = 125 individer per år.
Eksempel 2: Tolk gjærmodellen

Regresjonen gir omtrent N(t)=6641+67e−0,539tN(t) = \dfrac{664}{1 + 67e^{-0{,}539t}}. Tolk modellen.

Bæreevnen er omtrent 664. Kulturen kan ikke vokse mer enn det med den maten og plassen som finnes.
Vendepunktet er ved t=ln⁡670,539≈7,8t = \dfrac{\ln 67}{0{,}539} \approx 7{,}8 timer. Da var veksten raskest.
Største vekstfart er 664⋅0,5394≈89\dfrac{664 \cdot 0{,}539}{4} \approx 89 enheter per time.
Eksempel 3: Når er bestanden 900?

Med N(t)=10001+49e−0,5tN(t) = \dfrac{1000}{1 + 49e^{-0{,}5t}}: Når er N(t)=900N(t) = 900?

1+49e−0,5t=1000900=1091 + 49e^{-0{,}5t} = \dfrac{1000}{900} = \dfrac{10}{9}, så 49e−0,5t=1949e^{-0{,}5t} = \dfrac{1}{9}.
e−0,5t=1441e^{-0{,}5t} = \dfrac{1}{441}, så −0,5t=ln⁡1441-0{,}5t = \ln\dfrac{1}{441}.
t=2ln⁡441≈12,2t = 2\ln 441 \approx 12{,}2 år.
Vanlige feil

Vanlige feil

  • Å tro at bestanden kan passere KK. I modellen nærmer N(t)N(t) seg KK nedenfra, men når den aldri.
  • Å bruke eksponentiell modell på data som flater ut. Den eksponentielle modellen vil spå ubegrenset vekst.
  • Å tro at vekstfarten er størst i starten. Den er størst i vendepunktet, når bestanden er halvparten av bæreevnen.
  • Å lese aa som startverdien. Startverdien er K1+a\frac{K}{1 + a}.
Huskelapp

Huskelapp

  • N(t)=K1+ae−btN(t) = \dfrac{K}{1 + a e^{-bt}}, med bæreevne KK og N(0)=K1+aN(0) = \dfrac{K}{1 + a}.
  • lim⁡t→∞N(t)=K\lim\limits_{t \to \infty} N(t) = K (horisontal asymptote).
  • Vendepunkt: N=K2N = \frac{K}{2} ved t=ln⁡abt = \frac{\ln a}{b}.
  • Største vekstfart: Kb4\frac{Kb}{4}.
  • Starten ligner eksponentiell vekst.
Oppgaver

Øv selv