b) Når veks aurebestanden mest, og kor stor er veksten då?
Løysing
Vi kopierer funksjonen frå regresjonsanalyseverktøyet til grafikkfeltet. Så bruker vi CAS til å analysere funksjonen. Her må vi finne det punktet der grafen er brattast. Det vil vere i vendepunktet til grafen til .
I linje 3 og 4 kontrollerer vi at nullpunktet til den dobbeltderiverte faktisk er eit vendepunkt ved å sjekke om forteiknet skiftar ved nullpunktet, noko det gjer. Aurebestanden vaks raskast nesten 5 år etter at kalkinga starta, og då vaks bestanden med 2 740 aurar per år.
I oppgåve FM-23 viser vi at vi alltid har eitt vendepunkt for ein logistisk funksjon, og at funksjonen er brattast der, så vi treng ikkje å sjekke dette fleire gonger.
FM-21
I 1970 dukka det for første gong opp løvetann på Fureøya. I nokre år framover vart talet på løvetannplantar kartlagt. Resultatet viser vi i tabellen nedanfor.
Talet på løvetannplantar
År
1970
1971
1971
1973
1974
1975
1976
1977
1978
1979
1980
Talet på løvetann-plantar
3
15
80
201
300
501
731
915
1 131
1 350
1 490
a) Finn den logistiske modellen nx som passar best til utviklinga av talet på løvetannplantar der x står for talet på år etter 1970.
Løysing
Når x betyr talet på år etter 1970, får vi at x=0 betyr året 1970, x=1 betyr 1971 og så vidare. Vi legg derfor tala frå og med 0 til og med 10 i éin kolonne i reknearket i GeoGebra og talet på løvetannplantar i neste kolonne. Så vel vi "Regresjonsanalyse" og regresjonsmodellen logistisk.
Biletet viser at ein logistisk modell passar bra med tala. Den modellen som passar best, er
nx=1675,391+46,3e-0,58x
b) Kor mange løvetannplantar blir det på lang sikt etter denne modellen?
Løysing
Når x→∞, vil e-0,58x→0 og nemnaren gå mot 1. Funksjonen går derfor mot talet i teljaren.
Talet på løvetannplantar på lang sikt, eller det vi kallar bereevna, vil derfor vere 1 675 etter modellen.
c) Når vaks talet på løvetannplantar raskast, og kor raskt vaks det då?
Løysing
Vi kopierer funksjonen frå regresjonsanalyseverktøyet til grafikkfeltet. Så bruker vi CAS til å analysere funksjonen.
Etter modellen auka talet på løvetannplantar mest i det sjette året etter 1970, det vil seie i 1976. Då auka talet på løvetannplantar med 244 plantar per år.
d) Bruk regresjon i Python til å lage ein modell som passar best mogleg med tala i tabellen. Teikn funksjonen.
Løysing
Forslag til kode:
Koden gir denne utskrifta:
"Funksjonen blir n(x) = 1675/(1+46.30e^(-0.58x))."
Kommentar til koden: I linje 12 har vi lagt til p0 = [max(y_verdiar),1,1] i kommandoen curve_fit. Dersom vi ikkje gjer det, feiler regresjonen. Tala max(y_verdiar), 1 og 1 er startverdiar for dei tre variablane B, a og k, som programmet bruker i dei numeriske berekningane for å komme fram til dei endelege verdiane for dei tre variablane. Jo rettare startverdiane er, jo større er sjansen for at dei numeriske berekningane blir vellykka.
max(y_verdiar) er eit godt val for startverdien til B fordi han plukkar ut den største y-verdien i målingane, og vi veit at i ein logistisk modell er B den øvre grenseverdien for funksjonen. For dei to andre parametrane a og k set vi berre startverdien lik 1, som er den automatiske startverdien for parametrar ved bruk av curve_fit.
FM-22
Talet på hjort har auka kraftig i eit område. Tabellen viser talet på hjort for enkelte år i perioden 1973–2003.
Talet på hjort
Årstal
1973
1976
1979
1983
1988
1993
1998
2003
Talet på hjort,N
500
840
1 350
2 250
3 500
4 300
4 700
4 900
a) Finn ein logistisk funksjon Nt som passar med tala i tabellen. La t stå for talet på år etter 1973.
Løysing
Vi skriv inn tala i tabellen inn i reknearkdelen i GeoGebra og reknar ut verdiane for t. Så vel vi verktøyet "Regresjonsanalyse" og regresjonsmodellen "Logistisk".
Vi ser at grafen som kjem opp, passar godt til punkta. Den logistiske funksjonen som passar best med tala i tabellen, er
Nt=49881+9,10e-0,203t
b) Teikn grafen til den logistiske funksjonen du fann i a), og punkta han er basert på, i det same koordinatsystemet.
Løysing
Vi kopierer grafen (inkludert punkta) frå regresjonsverktøyet over i grafikkfeltet.
c) Kva er bereevna til hjortebestanden i dette området?
Løysing
Her må vi basere oss på modellen. Når t→∞, vil e-0,203t→0, nemnaren i funksjonen, gå mot 1, og talet på hjort vil gå mot talet i teljaren. Det er dette vi kallar bereevna for hjortebestanden.
Bereevna er 4 988 hjort etter modellen, noko vi rundar av til 5 000.
d) Når vaks hjortebestanden raskast, og kor raskt vaks han då?
Løysing
Vi må finne vendepunktet til funksjonen, og vi bruker CAS.
Hjortebestanden vaks mest omtrent 11 år etter 1973, det vil seie i 1984, og då vaks bestanden med 253 dyr per år.
e) Kor mange hjort var det då hjortebestanden vaks raskast i forhold til bereevna?
Løysing
Vi får at hjortebestanden var på halvparten av bereevna då han vaks raskast.
Vil det vere slik for alle logistiske modellar, trur du?
FM-23
Vi skal sjå litt meir på den generelle logistiske funksjonen
fx=B1+a·e-kx
Ved å analysere han kan vi få eit godt bilete på korleis ein logistisk modell fungerer. Vi føreset at konstantane a,k og B alle er større enn null slik dei er når det er snakk om vekst i ein bestand.
a) Bruk CAS og finn ut for kva x-verdi funksjonen veks raskast. Kor stor er veksten då?
Løysing
Vi får at grafen til f veks mest når x=lnak, og då veks han med 14Bk per eining på x-aksen. Det betyr at for alle logistiske modellar kan vi finne når veksten er størst, og kor stor han er, ved å bruke resultata i linje 2 og linje 4.
For å sjekke at nullpunktet til den dobbeltderiverte faktisk er eit vendepunkt, kan vi ikkje løyse ulikskapen f''x>0 slik vi vanlegvis ville ha gjort fordi GeoGebra ikkje veit at konstantane a,k og B er positive slik føresetnaden er. Det vi kan gjere, er å bruke den trippelderiverte funksjonen. Eit vendepunkt er eit ekstremalpunkt for den deriverte. Då kan vi bruke dobbeltderiverttesten på den deriverte (altså den trippelderiverte til f) for å sjekke om det er eit toppunkt eller eit botnpunkt for den deriverte. Vi ser av linje 3 at den trippelderiverte er negativ, som betyr at grafen til den deriverte vender den hole sida ned. Dette viser at den deriverte har eit toppunkt der, og vi har vist at veksten er størst der.
b) Dersom fx er talet på individ i ein bestand, kor mange individ er det når bestanden veks raskast?
Løysing
Sidan B er den verdien funksjonen går mot når x→∞, det vi kallar bereevna, er talet på individ når funksjonen veks raskast, lik halve bereevna.
c) Er det mogleg å løyse oppgåve a) og b) ved rekning for hand? Gjer eit forsøk!
I løysinga har vi enkelt kvitta oss med nemnaren og uttrykket framfor parentesen i teljaren sidan begge desse ikkje kan vere null. Vi må no sjekke at løysinga faktisk er x-koordinaten til eit vendepunkt. Det er slitsamt å skulle derivere éin gong til for hand slik vi gjorde enkelt med CAS. Vi prøver å løyse ulikskapen f''0>0 slik vi vanlegvis gjer i staden.
Kommentar til løysinga: Sidan dei tre konstantane a,k og B alle er større enn null, veit vi at vi kan dele med dei utan å måtte snu ulikskapsteiknet i overgangen mellom tredje og fjerde linje.
Vi får altså at grafen vender den hole sida si opp på venstresida av nullpunktet og motsett på høgre side, som vi forventa. Då veit vi at løysinga er x-koordinaten til eit vendepunkt, og vidare at grafen vil stige raskast her.