3.6 Numerisk derivasjon med digitale verktøy

Fremover-, bakover- og sentraldifferanse, og bruk av regneark og CAS.

45 min

10 oppgaver

Numerisk derivasjonFremoverdifferanseSentraldifferanseRegnearkCAS

Du leser den tradisjonelle versjonen

Din fremgang i kapitlet

0 / 10 oppgaver

Kapitlets plass i kurset

Bygger på

3.1Grenseverdi

Numerisk derivasjon

Ikke alle funksjoner kan deriveres analytisk på en enkel måte. I praksis har vi ofte bare måledata (tabellverdier) uten et eksplisitt funksjonsuttrykk. Da må vi tilnærme den deriverte numerisk.

Ideen er enkel: vi bruker differansekvotienten $\displaystyle \frac{f(a+h) - f(a)}{h}$ med en liten, men endelig $h$ i stedet for å ta grenseverdien $h \to 0$ . Spørsmålet er: hvilken differanseformel gir best tilnærming?

Tre differanseformler

f

✏️Eksempel 1: Sammenligning av differanseformler

La $f(x) = x^3$ . Tilnaerm $f'(2) = 12$ med $h = 0{,}1$ ved hjelp av de tre differanseformlene.

Løsning:

$f(2) = 8$ , $f(2{,}1) = 9{,}261$ , $f(1{,}9) = 6{,}859$ .

Framoverdifferanse:
$\frac{9{,}261 - 8}{0{,}1} = 12{,}61 \quad \text{(feil: } 0{,}61\text{)}$

Bakoverdifferanse:
$\frac{8 - 6{,}859}{0{,}1} = 11{,}41 \quad \text{(feil: } 0{,}59\text{)}$

Sentraldifferanse:
$\frac{9{,}261 - 6{,}859}{0{,}2} = 12{,}01 \quad \text{(feil: } 0{,}01\text{)}$

Sentraldifferansen gir en mye bedre tilnærming!

📜Feilorden for differanseformlene

Feilen i tilnærmingene har følgende størrelseorden:

- Framoverdifferanse: Feil $\sim O(h)$ (proporsjonalt med $h$ )
- Bakoverdifferanse: Feil $\sim O(h)$ (proporsjonalt med $h$ )
- Sentraldifferanse: Feil $\sim O(h^2)$ (proporsjonalt med $h^2$ )

Det betyr at når vi halverer $h$ :
- Framover/bakover: feilen halveres omtrent
- Sentral: feilen reduseres til omtrent en fjerdedel

Sentraldifferansen konvergerer altså mye raskere mot den eksakte verdien.

📝Oppgave 1

La $f(x) = \sin x$ . Tilnaerm $f'(\pi/4)$ med sentraldifferanse.

Bruk $h = 0{,}1$ .

Bruk $h = 0{,}01$ .

Sammenlign feilene. Stemmer $O(h^2)$ -estimatet?

Numerisk derivasjon med regneark

I et regneark kan vi sette opp numerisk derivasjon systematisk:

A	B	C
$x$	$f(x)$	$f'(x) \approx$
$a - h$	$=f(A2)$
$a$	$=f(A3)$	$=(B4-B2)/(A4-A2)$
$a + h$	$=f(A4)$

Kolonne C bruker sentraldifferanse:

\displaystyle \frac{f(a+h) - f(a-h)}{2h}

.
Eksempel i regneark-syntaks:
- Celle A2: =A3-0.001
- Celle B2: =A2^3 (for

f(x) = x^3

)
- Celle C3: =(B4-B2)/(A4-A2)

📝Oppgave 2

Sett opp et regneark for å tilnærme den deriverte.

Lag en tabell med $x$ -verdier fra $0$ til $5$ med steglengde $0{,}5$ for $f(x) = x^2 - 3x + 1$ . Beregn sentraldifferansen i hvert indre punkt.

Gjenta med steglengde $0{,}1$ . Hvor mye bedre blir tilnærmingen?

Numerisk og symbolsk derivasjon med CAS

CAS-verktøy (Computer Algebra System) kan både:

1. Derivere symbolsk -- finne det eksakte uttrykket for $f'(x)$
2. Beregne numerisk -- gi en tilnærmet verdi av $f'(a)$

I GeoGebra CAS:
- Symbolsk: Derivative[x^3 + sin(x)] gir 3x^2 + cos(x)
- Numerisk: NDerivative[x^3 + sin(x), 2] gir en numerisk tilnærming av $f'(2)$

I Python med SymPy:

from sympy import *
x = symbols('x')
f = x**3 + sin(x)
diff(f, x)        # Symbolsk: 3*x**2 + cos(x)
diff(f, x).subs(x, 2).evalf()  # Numerisk verdi i x=2

📝Oppgave 3

Bruk CAS eller regneark til å løse oppgavene.

Finn $f'(x)$ symbolsk for $f(x) = x \cdot e^{-x^2}$ .

Bruk numerisk derivasjon (sentraldifferanse med $h = 0{,}001$ ) til å tilnærme $f'(1)$ .

Sammenlign det numeriske og analytiske svaret.

Løs oppgaven Tren

Avrundingsfeil ved veldig liten h

Man skulle tro at mindre

h

alltid gir bedre resultat. Men i praksis får vi problemer med avrundingsfeil når

h

er veldig liten:

$\frac{f(a + h) - f(a)}{h} \approx \frac{\text{nesten like tall}}{\text{lite tall}}$

Når $f(a+h) \approx f(a)$ , taper vi gjeldende siffer i subtraksjonen. For dobbeltpresisjon (float64) er $h \approx 10^{-8}$ typisk et bra valg for sentraldifferanse.

Tommelregel: Bruk $h \approx \sqrt{\varepsilon}$ for framover/bakover og $h \approx \varepsilon^{1/3}$ for sentral, der $\varepsilon \approx 10^{-16}$ er maskinpresisjonen.

📝Oppgave 4

Undersoesk effekten av steglengde på noeyyaktigheten.

La $f(x) = e^x$ . Beregn sentraldifferansen for $f'(1) = e$ med $h = 10^{-1}, 10^{-2}, \ldots, 10^{-8}$ . Lag en tabell over $h$ og feil.

For hvilken $h$ får du den minste feilen?

📝Oppgave 5

Gitt tabellverdier, tilnaerm den deriverte.

Tabellen viser $f(x)$ : $f(1) = 3{,}0$ , $f(1{,}5) = 4{,}2$ , $f(2) = 5{,}8$ , $f(2{,}5) = 7{,}9$ , $f(3) = 10{,}5$ . Tilnaerm $f'(2)$ .

Tilnaerm $f'(1)$ og $f'(3)$ .

Løs oppgaven Tren

Når bruker vi numerisk vs. analytisk derivasjon?

- Analytisk (for hånd eller CAS): Når vi har et funksjonsuttrykk og trenger et eksakt svar. Gir fullstendig informasjon om $f'(x)$ for alle $x$ .
- Numerisk: Når vi bare har tabelldata, når funksjonsuttrykket er for komplisert, eller når vi trenger en rask tilnærming i ett punkt.

I praksis brukes ofte en kombinasjon: CAS finner det analytiske uttrykket, som så evalueres numerisk i ønskede punkter.

Oppsummering

- Den deriverte kan tilnærmes numerisk med differansekvotienter
- Sentraldifferansen $\displaystyle \frac{f(a+h) - f(a-h)}{2h}$ er mest noeyyaktig ( $O(h^2)$ -feil)
- Regneark og CAS er nyttige verktøy for numerisk derivasjon
- For liten $h$ gir avrundingsfeil -- det finnes en optimal steglengde
- Numerisk derivasjon brukes når analytisk derivasjon er upraktisk eller umulig

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.