Trapesintegrasjon

Trapesintegrasjon er en metode i numerisk analyse til beregning av verdien til et bestemt integral av en funksjon . Det gjøres ved å erstatte flaten under kurven som funksjonen beskriver, med et trapes eller flere med samme bredde. Nøyaktigheten av beregningen øker med antall slike trapeser man gjør bruk av.

 Bakgrunn 
Vanligvis kan ikke et integral beregnes eksakt. Man er derfor ofte henvist til å bruke numeriske metoder og elektroniske regnemaskiner. Disse tar vanligvis utgangspunkt i Riemann-summen som definerer integralet. Det er også hensiktsmessig å beregne funksjonsverdiene i diskrete punkt xi  med konstant avstand h.

Da integrasjon av en funksjon y = f(x) som beskriver en kurve i xy-planet, gir arealet under kurven, kan dette deles opp i mindre deler som kan beregnes mer nøyaktig. Trapesintegrasjon er basert på å dele arealet opp i mindre trapeser. Deles hele integrasjonen opp i n slike interval, har man for integralet



hvor x0 = a og xn = b er henholdsvis nedre og øvre grense i integralet, mens h = (b - a)/n. Denne fremgangsmåten kan gjøres mer nøyaktig i form av Simpson-integrasjon eller Romberg-integrasjon.

Trapesapproksimasjonen
Når funksjonen f(x) forandrer seg lite mellom integrasjonsgrensene a og b, er arealet under kurven tilnærmet lik med det til et trapes med to parallelle sider f(a) og f(b) med gjensidig avstand b - a. Derfor kan verdien til integralet skrives som I = T + E  hvor



er arealet til trapeset og E er feilen i tilnærmelsen. Jo mindre avstanden b - a er, desto mindre blir E. En mer nøyaktig analyse viser at



hvor z ligger i intervallet mellom a og b.K.E. Atkinson, An Introduction to Numerical Analysis, John Wiley & Sons, New York (1989). ISBN 0-471-50023-2.  Den viser det opplagte resultat at består funksjonen av rette linjer, så er den andrederiverte f " = 0 og tilnærmelsen er eksakt. For en konkav funksjon er f " < 0 slik at E > 0 og den eksakte verdien er større enn hva trapesapproksimasjonen gir. Det motsatte er tilfelle for en konveks funksjon.

Generelt tilfelle

I det generelle tilfellet når intervallet b - a  er vilkårlig stort, kan man dele hele arealet under kurvet opp i mindre trapeser som er så små at man kan bruke trapesapproksimasjonen for hver av dem. Hvis man velger en oppdeling med n slike trapeser, hver med samme høyde h = (b - a)/n, er verdien av  integralet tilnærmet lik



hvor nå unøyaktigheten er gitt ved



Den avtar derfor raskt ved økende antall n med trapesintervaller, det vil si desto mindre diskretiseringen h er. Hadde man i stedet valgt trapeser med forskjellige høyder, måtte man beregne mange flere funksjonsverdier i den tilsvarende summen. Også med like intervall vil en fordobling av deres antall medføre at halvparten av funksjonsverdiene forblir de samme. Det sparer også den tid som går med til en slik mer nøyaktig beregning.

Simpson-integrasjon

Da feilen ved trapesintegrasjon varierer med diskretiseringen som h 2, kan et mer presis resultat med nøyaktighet h 3 oppnås ved bruk av Richardson-ekstrapolasjon. Da T(h)  er resultatet av tilnærmelsen med n steg, vil en dobling av antall intervall gi resultatet T(h/2) i samme approksimasjon. Hvert punkt xk mellom a og b vil da gå over i punktet x2k, mens det nye punktet x2k+1 vil ligge mellom x2k og x2k+2 med x2n = b. Men dette kan nå gjøres enda mer nøyaktig ved bruk av kombinasjonen



hvor h1 = h/2  er den halverte intervallengden. Dette forbedrete resultatet for integralet går under navnet Simpson-integrasjon.

Da feilen i denne approksimasjonen viser seg å være



og ikke av orden h 3, er nøyaktigheten bedre enn forventet og bidrar til at denne integrasjonsmetoden er mye anvendt.R.L. Burden and D.J. Faires, Numerical Analysis,  Brooks/Cole, Boston (2005). ISBN 0-538-73351-9. Videre bruk av Richardson-ekstrapolasjon gjør den enda mer nøyaktig og blir omtalt som Romberg-integrasjon.

Referanser

 Eksterne lenker 

 «Trapezium formula», fra Encyclopedia of Mathematics
 «The Trapezoid Rule for Approximating Integrals», video fra YouTube (engelsk)
 «Trapezoidal Rule», video fra YouTube (engelsk)

Kategori:Matematisk fysikk
Kategori:Numerisk integrasjon