Mange integraler lar seg ikke løse med penn og papir på en praktisk måte. Da kan vi i stedet lage en algoritme som beregner en god tilnærming ved å dele opp arealet under kurven i mange små biter og summere dem. I R2 skal du både forstå hvordan algoritmene fungerer matematisk og kunne skrive dem som faktisk kjørende kode.
Numerisk integrasjon betyr at du bruker en algoritme til å beregne tilnærmede verdier av bestemte integraler når eksakt regning er vanskelig eller umulig. Du deler arealet under en kurve inn i mange små figurer, summerer dem, og kjører beregningen i et program som Python.
Et bestemt integral gir oss arealet under en kurve mellom to grenser. For enkle funksjoner kan vi bruke antideriverte og regne eksakt. Men mange funksjoner, for eksempel de som dukker opp i fysikk, biologi og økonomi, har ingen enkel antiderivert. Da er numerisk integrasjon løsningen.
Kompetansemålet handler om to ting: å forstå den matematiske ideen bak algoritmen, og å omsette den til et program. Det kobler saman abstrakt matematikk og praktisk problemløsning, noe som er sentralt i kjerneelementene modellering og abstraksjon i læreplanen.
Rektangelmetoden deler intervallet fra a til b inn i n like store delintervaller med bredde h = (b - a) / n. For hvert delintervall beregner du funksjonsverdien i et punkt, for eksempel venstre kant eller midtpunktet, og ganger med bredden h. Summen av alle disse rektanglene er tilnærmingen til integralet.
Trapesmetoden er litt mer nøyaktig. I stedet for et rektangel bruker du et trapes i hvert delintervall, altså du trekker en rett linje mellom de to endepunktene på kurven. Formelen blir: integral tilnærmet lik h/2 ganger (f(x_0) + 2*f(x_1) + 2*f(x_2) + ... + 2*f(x_{n-1}) + f(x_n)). Jo flere delintervaller du bruker, jo mindre blir feilen.
Midtpunktmetoden er en variant der du bruker funksjonsverdien i midten av hvert delintervall. Den gir ofte bedre nøyaktighet enn venstre- eller høyre-rektangelmetoden med samme antall intervaller.
Algoritmen kan skrives som en funksjon som tar inn funksjonsuttrykket, grensene a og b, og antall delintervaller n. Inne i funksjonen beregner du h, starter en sum på null, og bruker en løkke som går gjennom alle delintervallene og legger til bidraget fra hvert trapes eller rektangel.
Et eksempel på trapesmetoden i Python kan se slik ut: du definerer en funksjon trapes(f, a, b, n), beregner h = (b - a) / n, setter s = f(a) + f(b), og bruker en for-løkke fra 1 til n-1 der du legger til 2 * f(a + i * h). Til slutt returnerer du h / 2 * s.
Når du tester koden, bør du prøve den på et integral du kan beregne eksakt, for eksempel integralet av x^2 fra 0 til 1 som er 1/3. Da kan du se hvor stor feilen er, og observere at feilen minker når du øker n.
En vanlig feil er å blande indekser i løkken, slik at du enten hopper over et punkt eller teller et punkt dobbelt. Skriv ut noen mellomverdier mens du tester, så ser du raskt om noe er galt. En annen feil er å glemme at trapesmetoden dobler alle indre punkter men ikke endepunktene, noe som lett glipper.
Pass også på at du skiller mellom tilnærmingen og den eksakte verdien. Numerisk integrasjon gir alltid en feil, og det er en del av oppgaven å diskutere hvor stor feilen er og hva som påvirker den.
For å øve best: start med rektangelmetoden, test den på et enkelt integral, og sjekk feilen. Gå deretter videre til trapesmetoden og sammenlign. Prøv å øke n fra 10 til 100 til 1000 og se hva som skjer med nøyaktigheten. Til slutt kan du prøve å generalisere koden slik at den fungerer for hvilken som helst funksjon du sender inn som argument.
Bruk trapesmetoden med n = 4 til å tilnærme integralet av f(x) = x^2 fra 0 til 1.
Fasit: Tilnærmet verdi: 0,34375. Eksakt verdi: 1/3 ≈ 0,33333. Feil ≈ 0,01042.
Steg 1: Beregn h = (1 - 0) / 4 = 0,25. Steg 2: Finn x-verdiene: x0 = 0, x1 = 0,25, x2 = 0,5, x3 = 0,75, x4 = 1. Steg 3: Beregn funksjonsverdiene: f(0) = 0, f(0,25) = 0,0625, f(0,5) = 0,25, f(0,75) = 0,5625, f(1) = 1. Steg 4: Sett inn i trapesformelen: (0,25 / 2) * (0 + 2*0,0625 + 2*0,25 + 2*0,5625 + 1) = 0,125 * 2,75 = 0,34375.
Prøv Lekselos gratis i 14 dager for å lese hele forklaringen, se eksemplet steg for steg og øve på matematikk for realfag tilpasset deg.
Lekselos loser barnet gjennom oppgaven ett steg om gangen, så det lærer i stedet for å kopiere. Du får ukeplan, prøver og gloser samlet, og ser hva som er gjort. Forankret i LK20, laget i Norge. 14 dager gratis, uten kort.
Analytisk integrasjon gir en eksakt verdi ved hjelp av antideriverte og regneregler. Numerisk integrasjon gir en tilnærmet verdi ved hjelp av en algoritme, og brukes når eksakt løsning er vanskelig eller umulig å finne.
Trapesmetoden gir vanligvis bedre nøyaktighet enn rektangelmetoden for samme antall delintervaller, fordi den tar høyde for at kurven endrer seg mellom punktene. For glatte funksjoner er trapesmetoden et godt valg.
Feilen minker når du øker antall delintervaller n. For trapesmetoden er feilen omtrent proporsjonal med h^2, der h = (b - a) / n. Dobler du n, blir feilen omtrent fire ganger mindre.
Læreplanen sier ikke hvilket språk du må bruke, men Python er det vanligste valget i norsk skole. Du kan bruke vanlige Python-konstruksjoner som for-løkker og funksjoner uten å trenge eksterne biblioteker.
Du trenger ikke et formelt bevis, men du bør kunne forklare ideen bak metoden, sette opp formelen, bruke den på konkrete tall, og diskutere nøyaktighet. I tillegg kan du bli bedt om å skrive eller tolke kode som utfører algoritmen.
Forankring i læreplanen: Kompetansemål (LK20): Fag: Læreplan i matematikk for realfag (matematikk R) Trinn: Matematikk R2 Kompetansemål: Eleven skal kunne utvikle algoritmer for å beregne integraler numerisk, og bruke programmering til å utføre algoritmene. Nynorsk: Eleven skal kunne utvikle algoritmar for å berekne integral numerisk, og bruke programmering til å utføre algoritmane. Kjerneelement: Modellering og anvendelser, Abstraksjon og generalisering (Kilde: Utdanningsdirektoratet)
Lekselos loser deg gjennom oppgaver med hint, ikke fasit, og lager øving tilpasset deg. Prøv gratis.