Tilbake
1.5
Integrasjon ved substitusjon

1.5 Integrasjon ved substitusjon

Bruke substitusjon som metode for å løse integraler.

55 min
19 oppgaver
SubstitusjonVariabelskifteSammensatte funksjonerKjerneregelen baklengs
Du leser den tradisjonelle versjonen
Din fremgang i kapitlet
0 / 19 oppgaver
Kapitlets plass i kurset

Hvorfor numerisk integrasjon?

Mange funksjoner lar seg ikke integrere analytisk. For eksempel har ex2dx\int e^{-x^2} \, dx ingen elementær antiderivert. I slike tilfeller bruker vi numeriske metoder for å beregne tilnærmede verdier av bestemte integraler.

Numerisk integrasjon er også svært nyttig når vi bare har datapunkter (målinger) og ikke kjenner funksjonsuttrykket.

Numerisk integrasjon
Numerisk integrasjon er metoder for å beregne tilnærmede verdier av bestemte integraler abf(x)dx\int_a^b f(x) \, dx.

Grunnideen er å dele intervallet [a,b][a, b] inn i nn delintervaller, hver med bredde h=ban\displaystyle h = \frac{b - a}{n}, og tilnærme arealet under grafen med enklere geometriske figurer.

Delingspunktene er x0=a, x1=a+h, x2=a+2h, , xn=bx_0 = a, \ x_1 = a + h, \ x_2 = a + 2h, \ \ldots, \ x_n = b.

📜Trapesformelen
Med nn delintervaller og h=ban\displaystyle h = \frac{b-a}{n} er trapesformelen:

Tn=h2[f(x0)+2f(x1)+2f(x2)++2f(xn1)+f(xn)]T_n = \frac{h}{2}\Big[f(x_0) + 2f(x_1) + 2f(x_2) + \cdots + 2f(x_{n-1}) + f(x_n)\Big]

=h2[f(x0)+f(xn)+2k=1n1f(xk)]= \frac{h}{2}\left[f(x_0) + f(x_n) + 2\sum_{k=1}^{n-1} f(x_k)\right]

Trapesformelen tilnærmer arealet ved å erstatte kurvebuen mellom to nabopunkter med en rett linje, slik at hvert delområde blir et trapes.

✏️Trapesformelen med 4 delintervaller

Bruk trapesformelen med n=4n = 4 til å tilnærme 02ex2dx\int_0^2 e^{-x^2} \, dx.

Vi har a=0a = 0, b=2b = 2, n=4n = 4 og h=204=0,5\displaystyle h = \frac{2-0}{4} = 0{,}5.

Delingspunktene og funksjonsverdiene:

kkxkx_kf(xk)=exk2f(x_k) = e^{-x_k^2}
00e0=1e^0 = 1
10,5e0,250,7788e^{-0{,}25} \approx 0{,}7788
21,0e10,3679e^{-1} \approx 0{,}3679
31,5e2,250,1054e^{-2{,}25} \approx 0{,}1054
42,0e40,0183e^{-4} \approx 0{,}0183

T4=0,52[1+2(0,7788)+2(0,3679)+2(0,1054)+0,0183]T_4 = \frac{0{,}5}{2}\Big[1 + 2(0{,}7788) + 2(0{,}3679) + 2(0{,}1054) + 0{,}0183\Big]
=0,25[1+1,5576+0,7358+0,2108+0,0183]= 0{,}25 \cdot [1 + 1{,}5576 + 0{,}7358 + 0{,}2108 + 0{,}0183]
=0,253,52250,8806= 0{,}25 \cdot 3{,}5225 \approx 0{,}8806
Den eksakte verdien er 0,8821\approx 0{,}8821, så tilnærmingen er ganske god.
📜Simpsons regel
Med et partall nn delintervaller og h=ban\displaystyle h = \frac{b-a}{n} er Simpsons regel:

Sn=h3[f(x0)+4f(x1)+2f(x2)+4f(x3)+2f(x4)++4f(xn1)+f(xn)]S_n = \frac{h}{3}\Big[f(x_0) + 4f(x_1) + 2f(x_2) + 4f(x_3) + 2f(x_4) + \cdots + 4f(x_{n-1}) + f(x_n)\Big]

Koeffisientene følger mønsteret 1,4,2,4,2,,4,2,4,11, 4, 2, 4, 2, \ldots, 4, 2, 4, 1.

Simpsons regel tilnærmer kurvebuen mellom tre nabopunkter med en parabel (andregradsfunksjon), noe som gir langt bedre nøyaktighet enn trapesformelen.

✏️Simpsons regel med 4 delintervaller

Bruk Simpsons regel med n=4n = 4 til å tilnærme 02ex2dx\int_0^2 e^{-x^2} \, dx.

Vi bruker de samme verdiene som i forrige eksempel: h=0,5h = 0{,}5.

S4=0,53[f(0)+4f(0,5)+2f(1)+4f(1,5)+f(2)]S_4 = \frac{0{,}5}{3}\Big[f(0) + 4f(0{,}5) + 2f(1) + 4f(1{,}5) + f(2)\Big]

=0,53[1+4(0,7788)+2(0,3679)+4(0,1054)+0,0183]= \frac{0{,}5}{3}\Big[1 + 4(0{,}7788) + 2(0{,}3679) + 4(0{,}1054) + 0{,}0183\Big]

=0,53[1+3,1152+0,7358+0,4216+0,0183]= \frac{0{,}5}{3}\Big[1 + 3{,}1152 + 0{,}7358 + 0{,}4216 + 0{,}0183\Big]

=0,535,29090,8818= \frac{0{,}5}{3} \cdot 5{,}2909 \approx 0{,}8818

Den eksakte verdien er 0,8821\approx 0{,}8821. Simpsons regel gir her en feil på bare 0,0003\approx 0{,}0003, mye bedre enn trapesformelen (0,0015\approx 0{,}0015).

📊Utforsk: Trapesmetoden

Trapesmetoden med 4 trapeser tilnærmer 02(x2+1)dx\int_0^2 (x^2 + 1)\, dx. Øk antall trapeser (siste tall i kommandoen) og se hvordan tilnærmingen nærmer seg den eksakte verdien 143\displaystyle \frac{14}{3}.

📜Feilestimering
Feil for trapesformelen:
ET(ba)312n2M2|E_T| \leq \frac{(b-a)^3}{12n^2} \cdot M_2
der M2=maxaxbf(x)M_2 = \max_{a \leq x \leq b} |f''(x)|.

Feil for Simpsons regel:
ES(ba)5180n4M4|E_S| \leq \frac{(b-a)^5}{180n^4} \cdot M_4
der M4=maxaxbf(4)(x)M_4 = \max_{a \leq x \leq b} |f^{(4)}(x)|.

Simpsons regel konvergerer altså som 1n4\displaystyle \frac{1}{n^4}, mens trapesformelen konvergerer som 1n2\displaystyle \frac{1}{n^2}. Dobling av nn gir ca. 4 ganger bedre nøyaktighet for trapes, men 16 ganger bedre for Simpson.

✏️Feilestimering og valg av n

Hvor mange delintervaller nn trenger vi i trapesformelen for å beregne 0111+x2dx\displaystyle \int_0^1 \frac{1}{1+x^2} \, dx med feil mindre enn 0,0010{,}001?

Vi har f(x)=11+x2\displaystyle f(x) = \frac{1}{1+x^2}.

f(x)=2x(1+x2)2\displaystyle f'(x) = \frac{-2x}{(1+x^2)^2}

f(x)=2(1+x2)2+2x2(1+x2)2x(1+x2)4=6x22(1+x2)3\displaystyle f''(x) = \frac{-2(1+x^2)^2 + 2x \cdot 2(1+x^2) \cdot 2x}{(1+x^2)^4} = \frac{6x^2 - 2}{(1+x^2)^3}

For x[0,1]x \in [0, 1]: f(x)|f''(x)| har maksimum i x=0x = 0 der f(0)=2|f''(0)| = 2. Altså M2=2M_2 = 2.

Feilgrensen: ET(10)312n22=16n2\displaystyle |E_T| \leq \frac{(1-0)^3}{12n^2} \cdot 2 = \frac{1}{6n^2}

Vi krever 16n2<0,001\displaystyle \frac{1}{6n^2} < 0{,}001, altså n2>10,006166,7\displaystyle n^2 > \frac{1}{0{,}006} \approx 166{,}7, som gir n>12,9n > 12{,}9.

Vi trenger altså minst n=13n = 13 delintervaller.

Praktisk bruk

Når bruker vi numerisk integrasjon?

1. Funksjonen har ingen elementær antiderivert: ex2dx\int e^{-x^2}\,dx, sinxxdx\displaystyle \int \frac{\sin x}{x}\,dx, 1+x3dx\int \sqrt{1 + x^3}\,dx
2. Funksjonen er gitt som datapunkter: Målinger fra eksperimenter, fartsmålinger, temperaturdatalogger
3. Kontrollregning: Sjekke analytiske svar med numerisk tilnærming

Tommelregel for valg av metode:
- Trapesformelen er enklest å huske og bruke
- Simpsons regel gir bedre nøyaktighet med samme antall beregninger
- Dobling av nn og sammenligning gir praktisk feilkontroll

✏️Numerisk integrasjon fra tabell

Farten til en bil (i m/s) er målt hvert 2. sekund:

tt (s)0246810
v(t)v(t) (m/s)0512182225

Beregn tilnærmet tilbakelagt strekning s=010v(t)dts = \int_0^{10} v(t)\,dt med trapesformelen.

Her er n=5n = 5, h=2h = 2 s, a=0a = 0 og b=10b = 10.

T5=22[v(0)+2v(2)+2v(4)+2v(6)+2v(8)+v(10)]T_5 = \frac{2}{2}\Big[v(0) + 2v(2) + 2v(4) + 2v(6) + 2v(8) + v(10)\Big]

=1[0+2(5)+2(12)+2(18)+2(22)+25]= 1 \cdot \Big[0 + 2(5) + 2(12) + 2(18) + 2(22) + 25\Big]

=0+10+24+36+44+25=139 m= 0 + 10 + 24 + 36 + 44 + 25 = 139 \text{ m}

Bilen har tilbakelagt omtrent 139139 meter i løpet av 1010 sekunder.

📝Oppgave 1

Løs oppgavene:

a

Beregn 01x2dx\int_0^1 x^2\,dx med trapesformelen og n=4n = 4. Sammenlign med eksakt verdi.

b

Beregn samme integral med Simpsons regel og n=4n = 4. Sammenlign feilen.

📝Oppgave 2

Løs oppgavene:

a

Beregn 131xdx\displaystyle \int_1^3 \frac{1}{x}\,dx med trapesformelen og n=4n = 4.

b

Finn eksakt verdi av integralet og beregn feilen.

📝Oppgave 3

Løs oppgavene:

a

Skriv opp formelen for S6S_6 (Simpsons regel med 6 delintervaller) for 03f(x)dx\int_0^3 f(x)\,dx.

b

Hva er spesielt med Simpsons regel når funksjonen er et polynom av grad 3\leq 3?

📝Oppgave 4

Løs oppgavene:

a

Vannføringen i en elv (i m3/s\text{m}^3/\text{s}) er målt hver time:

tt (timer)01234
Q(t)Q(t)1014201812

Beregn total vannmengde (i m3\text{m}^3) som har passert i løpet av 4 timer, med trapesformelen.

📝Oppgave 5

Løs oppgavene:

a

Beregn 01ex2dx\int_0^1 e^{-x^2}\,dx med trapesformelen for n=2n = 2 og n=4n = 4. Kommenter konvergensen.

b

Beregn samme integral med Simpsons regel for n=4n = 4.

📝Oppgave 6

Løs oppgavene:

a

Bruk feilformelen for trapesmetoden til å vise at n10n \geq 10 delintervaller er nok for å beregne 01exdx\int_0^1 e^x\,dx med feil <0,001< 0{,}001.

📝Oppgave 7

Løs oppgavene:

a

Beregn 0πsinxdx\int_0^{\pi} \sin x \, dx med Simpsons regel og n=6n = 6. Sammenlign med eksakt verdi.

📝Oppgave 8

Løs oppgavene:

a

En funksjon ff er kjent i følgende punkter:

xx00,511,52
f(x)f(x)1,001,221,652,323,50

Beregn 02f(x)dx\int_0^2 f(x)\,dx med både trapesformelen og Simpsons regel.

📝Oppgave 9

Løs oppgavene:

a

Vis at trapesformelen med n=1n = 1 gir T1=ba2[f(a)+f(b)]\displaystyle T_1 = \frac{b-a}{2}[f(a) + f(b)], og tegn en figur som illustrerer dette.

📝Oppgave 10

Løs oppgavene:

a

Finn det minste nn slik at Simpsons regel gir 0111+x2dx\displaystyle \int_0^1 \frac{1}{1+x^2}\,dx med feil <106< 10^{-6}.

📝Oppgave 11

Løs oppgavene:

a

Beregn 011+x3dx\int_0^1 \sqrt{1 + x^3}\,dx med Simpsons regel og n=4n = 4.

📝Oppgave 12

Løs oppgavene:

a

Forklar hvorfor trapesformelen systematisk overestimerer integralet dersom ff er konveks (f(x)>0f''(x) > 0) på hele [a,b][a, b].

b

Gjelder tilsvarende for Simpsons regel? Begrunn svaret.

📝Oppgave 13

Løs oppgavene:

a

Akselerasjonen a(t)a(t) til en rakett de første 10 sekundene er målt:

tt (s)0246810
a(t)a(t) (m/s2^2)0815202325

Raketten starter fra ro. Beregn farten v(10)=010a(t)dtv(10) = \int_0^{10} a(t)\,dt med Simpsons regel.

📝Oppgave 14

Løs oppgavene:

a

Integralet 01sinxxdx\displaystyle \int_0^1 \frac{\sin x}{x}\,dx kan ikke beregnes analytisk. Bruk Simpsons regel med n=6n = 6 til å finne en tilnærmet verdi. (Hint: sinxx1\displaystyle \frac{\sin x}{x} \to 1 når x0x \to 0.)

Oppsummering

Trapesformelen: Tn=h2[f(x0)+2f(x1)++2f(xn1)+f(xn)]\displaystyle T_n = \frac{h}{2}[f(x_0) + 2f(x_1) + \cdots + 2f(x_{n-1}) + f(x_n)]
- Koeffisienter: 1,2,2,,2,11, 2, 2, \ldots, 2, 1
- Feil: O(h2)O(h^2) - konvergerer som 1n2\displaystyle \frac{1}{n^2}

Simpsons regel: Sn=h3[f(x0)+4f(x1)+2f(x2)++4f(xn1)+f(xn)]\displaystyle S_n = \frac{h}{3}[f(x_0) + 4f(x_1) + 2f(x_2) + \cdots + 4f(x_{n-1}) + f(x_n)]
- Koeffisienter: 1,4,2,4,2,,4,11, 4, 2, 4, 2, \ldots, 4, 1 (krever partall nn)
- Feil: O(h4)O(h^4) - konvergerer som 1n4\displaystyle \frac{1}{n^4}
- Eksakt for polynomer av grad 3\leq 3

Praktisk: Doble nn og sammenlign - dersom svaret endrer seg lite, har vi god nøyaktighet.

Repetisjonsoppgaver
Din fremgang
0deloppgaver0 / 5 oppgaver

Dette kapitlet er skrevet av Anthropics toppmodeller (Claude Opus og Claude Fable) og er foreløpig ikke manuelt gjennomgått — kvalitetskontrollen gjøres av uavhengige KI-agenter, og innmeldte feil rettes fortløpende. Funnet en feil? Meld fra, så retter vi den. Les mer om hvordan innholdet lages.