Trapesmetoden i integraler
I anvendt matematikk støter vi ofte på problemet med å beregne arealet under en kurve. Teoretisk sett kan dette løses ved hjelp av bestemte integraler. I praksis er imidlertid ikke alle funksjoner enkle å integrere analytisk. Noen funksjoner er komplekse, data er bare tilgjengelige i tabellform, eller modeller mangler enkle antidiriverte. I slike situasjoner blir numeriske metoder essensielle. En av de mest kjente og hyppig brukte numeriske metodene er trapesregelen, en tilnærmingsteknikk for integraler som deler arealet under en kurve inn i trapesformede former.
Grunnleggende konsept med bestemt integral
Geometrisk kan det bestemte integralet \(\int_a^bf(x)\,dx\) forstås som arealet avgrenset av kurven \(y=f(x)\), \(x\)-aksen og de vertikale linjene \(x=a\) og \(x=b\). Hvis \(f(x)\ge 0\), er integralet et positivt areal. Hvis funksjonen har negative verdier på et visst intervall, gir integralet et "fortegnsareal".
Hovedproblemet oppstår når:
1. Funksjonen \(f(x)\) har ingen antiderivasjon som kan uttrykkes av elementærfunksjoner.
2. Verdien av \(f(x)\) er bare kjent på noen punkter (eksperimentelle data).
3. Symbolsk beregning er ineffektiv eller umulig.
Det er her trapesmetoden brukes til å tilnærme integralverdien ved å bruke funksjonsverdiene på bestemte punkter.
Idé for trapesformet metode
Trapesmetoden starter med en enkel idé: tilnærme en kurve med en rett linje i hvert delintervall. Mens Riemann-summetoden bruker vinkelrette sider og et rektangulært areal, bruker trapesmetoden to endepunkter på delintervallet og forbinder dem med en rett linje. Arealet under den rette linjen danner en trapes, ikke et rektangel.
Anta at vi ønsker å beregne:
\[
\int_a^bf(x)\,dx
\]
Vi deler intervallet \([a,b]\) inn i \(n\) delintervaller av samme lengde. Lengden på hvert delintervall er:
\[
h = \frac{ba}{n}
\]
Delepunkter:
\[
x_0=a,\;x_1=a+h,\;x_2=a+2h,\;\punkter,\;x_n=b
\]
Med funksjonsverdier:
\[
f(x_0), f(x_1), ∫, f(x_n)
\]
På hvert delintervall \([x_{i-1}, x_i]\) er arealet under kurven tilnærmet beregnet med trapesarealet:
\[
A_i \approx \frac{h}{2}\venstre[f(x_{i-1}) + f(x_i)\høyre]
\]
Så får man det totale arealet (integraltilnærmingen) ved å legge sammen alle trapesene:
\[
\int_a^bf(x)\,dx \approx \sum_{i=1}^{n}\frac{h}{2}\left[f(x_{i-1}) + f(x_i)\right]
\]
Hvis det ryddes opp, blir formelen for den sammensatte trapesformede regelen:
\[
\int_a^b f(x)\,dx \approx \frac{h}{2}\left[f(x_0)+2f(x_1)+2f(x_2)+\cdots+2f(x_{n-1})+f(x_n)\right]
\]
Legg merke til koeffisientmønsteret: endepunktene \(x_0\) og \(x_n\) multipliseres med 1, mens punktene i midten multipliseres med 2.
Enkelt beregningseksempel
Anta at vi ønsker å nærme oss:
\[
\int_0^2 x^2\,dx
\]
Analytisk er integralet \(\left.\frac{x^3}{3}\right|_0^2=\frac{8}{3}\approx 2{,}6667\). Vi vil imidlertid bruke trapesmetoden med \(n=4\).
1. Intervall \([0,2]\), så:
\[
h=\frac{2-0}{4}=0{,}5
\]
2. Titik-titik: \(x_0=0\), \(x_1=0{,}5\), \(x_2=1\), \(x_3=1{,}5\), \(x_4=2\)
3. Funksjonsverdi:
\(f(0)=0\)
\(f(0{,}5)=0{,}25\)
\(f(1)=1\)
\(f(1{,}5)=2{,}25\)
\(f(2)=4\)
Bruk formelen for den sammensatte trapesformen:
\[
T=\frac{h}{2}\left[f(x_0)+2f(x_1)+2f(x_2)+2f(x_3)+f(x_4)\right]
\]
\[
T=\frac{0{,}5}{2}\left[0+2(0{,}25)+2(1)+2(2{,}25)+4\right]
\]
\[
T=0{,}25\left[0{,}5+2+4{,}5+4\right]=0{,}25(11)=2{,}75
\]
Det omtrentlige resultatet på 2,75 er ganske nær den nøyaktige verdien på 2,6667, med en feil på omtrent 0,0833.
Feil i trapesformet metode
Trapesmetoden er en tilnærmingsmetode, som betyr at det alltid er en forskjell mellom den omtrentlige verdien og den faktiske integralverdien. Størrelsen på feilen påvirkes av:
1. Antall delintervaller n: jo større n, desto mindre h, og generelt desto mer nøyaktige er resultatene.
2. Funksjonens krumning: hvis funksjonen er svært krum (har en stor andrederivert), kan trapesmetoden være mindre nøyaktig enn metoder av høyere orden.
For tilstrekkelig glatte funksjoner (som har kontinuerlige andrederiverte) kan feilen til den sammensatte trapesmetoden estimeres:
\[
E_T = -\frac{(ba)}{12}h^2 f"(\xi)
\]
for noen ξ mellom ξ og ξ. Fra denne formelen kan man se at feilen er proporsjonal med ξ. Det vil si at hvis vi halverer ξ (dvs. dobler antallet delintervaller), reduseres feilen med omtrent en fjerdedel.
Fordeler med trapesmetoden
Trapesmetoden er populær av flere grunner:
1. Enkelt og raskt
Beregningen er relativt enkel, og kan til og med gjøres manuelt for små tilfeller.
2. Egnet for tabelldata
Hvis verdien av \(f(x)\) bare er kjent i diskrete punkter (f.eks. måleresultater), kan denne metoden anvendes direkte.
3. Bedre nøyaktighet enn den rektangulære metoden
Fordi den bruker to endepunkter og tilnærmer funksjonen med en rett linje, er den generelt mer nøyaktig enn den rektangulære tilnærmingen for samme antall partisjoner.
4. Stabil for mange tekniske applikasjoner
Mange fysikk- og ingeniørproblemer (f.eks. beregning av arbeid, energi, tverrsnittsareal, utladning og så videre) kan løses med denne metoden.
Begrensninger ved trapesmetoden
Selv om denne metoden er nyttig, har den begrensninger:
1. Mindre nøyaktig for svært buede funksjoner
Fordi den bare bruker en lineær tilnærming på hvert delintervall, krever funksjoner med raske endringer eller skarpe kurver en stor \(n\).
2. Ikke en metode av høy kvalitet
Feil i størrelsesorden \(h^2\); andre metoder som Simpsons kan gi mindre feil over samme antall delintervaller (forutsatt at funksjonen er tilstrekkelig jevn).
3. Følsom for intervallvalg
Hvis intervallet er for bredt og \(n\) er lite, kan resultatene avvike betydelig.
Lukking
Trapesmetoden for integrasjon er en av de mest grunnleggende teknikkene innen numerisk integrasjon. Ved å dele et intervall inn i små deler og erstatte kurven med en rett linje i hver del, kan vi tilnærme integralet uten å måtte finne den antideriverte analytisk. Formelen er enkel, lett å implementere og spesielt nyttig når dataene er diskrete eller funksjonen er vanskelig å integrere.
Brukere av trapesmetoden må imidlertid være oppmerksomme på antall delintervaller og hvordan funksjonen som integreres oppfører seg. For svært krumme funksjoner eller når høy presisjon er nødvendig, kan trapesmetoden fortsatt brukes ved å øke \(n\), eller andre numeriske metoder som Simpsons eller Rombergs metoder kan vurderes. Med god forståelse blir trapesmetoden et effektivt og praktisk verktøy for å løse ulike integralproblemer innen vitenskap, ingeniørfag og økonomi.
