Trapeziummethode bij integralen
In de toegepaste wiskunde komen we vaak het probleem tegen van het berekenen van de oppervlakte onder een curve. Theoretisch kan dit worden opgelost met behulp van bepaalde integralen. In de praktijk zijn echter niet alle functies even gemakkelijk analytisch te integreren. Sommige functies zijn complex, gegevens zijn alleen beschikbaar in tabelvorm, of modellen missen eenvoudige primitieve functies. In dergelijke situaties worden numerieke methoden essentieel. Een van de bekendste en meest gebruikte numerieke methoden is de trapeziumregel, een benaderingstechniek voor integralen die de oppervlakte onder een curve verdeelt in trapeziumvormige segmenten.
Basisbegrip van de bepaalde integraal
Geometrisch gezien kan de bepaalde integraal \(\int_a^bf(x)\,dx\) worden opgevat als het oppervlak begrensd door de curve \(y=f(x)\), de \(x\)-as en de verticale lijnen \(x=a\) en \(x=b\). Als \(f(x)\ge 0\), is de integraal een positief oppervlak. Als de functie negatieve waarden heeft op een bepaald interval, geeft de integraal een "oppervlak met teken".
Het voornaamste probleem ontstaat wanneer:
1. De functie \(f(x)\) heeft geen primitieve functie die kan worden uitgedrukt door elementaire functies.
2. De waarde van \(f(x)\) is alleen op sommige punten bekend (experimentele gegevens).
3. Symbolische berekeningen zijn inefficiënt of onmogelijk.
Hier wordt de trapeziummethode gebruikt om de integraalwaarde te benaderen door gebruik te maken van de functiewaarden op bepaalde punten.
Trapeziummethode idee
De trapeziummethode begint met een eenvoudig idee: benader een curve met een rechte lijn in elk subinterval. Waar de Riemann-sommethode gebruikmaakt van loodrechte zijden en een rechthoekig oppervlak, gebruikt de trapeziummethode twee eindpunten op het subinterval en verbindt deze met een rechte lijn. Het oppervlak onder de rechte lijn vormt een trapezium, geen rechthoek.
Stel dat we het volgende willen berekenen:
\[
\int_a^bf(x)\,dx
\]
We verdelen het interval \([a,b]\) in \(n\) subintervallen van gelijke lengte. De lengte van elk subinterval is:
\[
h=\frac{ba}{n}
\]
Scheidingspunten:
\[
x_0=a,\;x_1=a+h,\;x_2=a+2h,\;\dots,\;x_n=b
\]
Met functiewaarden:
\[
f(x_0), f(x_1), \dots, f(x_n)
\]
Bij elk subinterval \([x_{i-1}, x_i]\) wordt het oppervlak onder de curve benaderd door oppervlakte van een trapezium:
\[
A_i \approx \frac{h}{2}\left[f(x_{i-1}) + f(x_i)\right]
\]
Vervolgens wordt de totale oppervlakte (integrale benadering) verkregen door alle trapeziums bij elkaar op te tellen:
\[
\int_a^bf(x)\,dx \approx \sum_{i=1}^{n}\frac{h}{2}\left[f(x_{i-1}) + f(x_i)\right]
\]
Na het vereenvoudigen wordt de formule van de samengestelde trapeziumregel als volgt:
\[
\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]
\]
Let op het coëfficiëntenpatroon: de eindpunten \(x_0\) en \(x_n\) worden vermenigvuldigd met 1, terwijl de punten in het midden worden vermenigvuldigd met 2.
Eenvoudig rekenvoorbeeld
Stel dat we het volgende willen benaderen:
\[
\int_0^2 x^2\,dx
\]
Analytisch gezien is de integraal \(\left.\frac{x^3}{3}\right|_0^2=\frac{8}{3}\approx 2{,}6667\). We zullen echter de trapeziummethode gebruiken met \(n=4\).
1. Interval \([0,2]\), dan:
\[
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. Functiewaarde:
\(f(0)=0\)
\(f(0{,}5)=0{,}25\)
\(f(1)=1\)
\(f(1{,}5)=2{,}25\)
\(f(2)=4\)
Gebruik de formule voor een samengesteld trapezium:
\[
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
\]
Het geschatte resultaat van 2,75 ligt vrij dicht bij de exacte waarde van 2,6667, met een foutmarge van ongeveer 0,0833.
Fout in de trapeziummethode
De trapeziummethode is een benaderingsmethode, wat betekent dat er altijd een verschil is tussen de benaderde waarde en de werkelijke integraalwaarde. De omvang van de fout wordt beïnvloed door:
1. Aantal subintervallen \(n\): hoe groter \(n\), hoe kleiner \(h\), en over het algemeen hoe nauwkeuriger de resultaten.
2. Kromming van de functie: als de functie sterk gekromd is (een grote tweede afgeleide heeft), kan de trapeziummethode minder nauwkeurig zijn dan methoden van hogere orde.
Voor voldoende gladde functies (met continue tweede afgeleiden) kan de fout van de samengestelde trapeziummethode worden geschat:
\[
E_T = -\frac{(ba)}{12}h^2 f”(\xi)
\]
voor een of andere \(\xi\) tussen \(a\) en \(b\). Uit deze formule blijkt dat de fout evenredig is met \(h^2\). Dat wil zeggen, als we \(h\) halveren (d.w.z. het aantal subintervallen verdubbelen), wordt de fout met ongeveer een kwart verminderd.
Voordelen van de trapeziummethode
De trapeziummethode is om verschillende redenen populair:
1. Eenvoudig en snel
De berekening is relatief eenvoudig en kan zelfs handmatig worden uitgevoerd voor kleine gevallen.
2. Geschikt voor tabelgegevens
Als de waarde van \(f(x)\) alleen op discrete punten bekend is (bijvoorbeeld meetresultaten), kan deze methode direct worden toegepast.
3. Betere nauwkeurigheid dan de rechthoekige methode
Omdat deze methode twee eindpunten gebruikt en de functie benadert met een rechte lijn, is ze over het algemeen nauwkeuriger dan de rechthoekige benadering voor hetzelfde aantal partities.
4. Stabiel voor velen technische toepassingen
Veel natuurkundige en technische problemen (bijvoorbeeld het berekenen van arbeid, energie, dwarsdoorsnedeoppervlak, debiet, enzovoort) kunnen met deze methode worden aangepakt.
Beperkingen van de trapeziummethode
Hoewel deze methode nuttig is, kent ze wel beperkingen:
1. Minder nauwkeurig voor sterk gekromde functies
Omdat er op elk subinterval alleen een lineaire benadering wordt gebruikt, vereisen functies met snelle veranderingen of scherpe krommen een grote \(n\).
2. Geen hogere-ordemethode
Fouten van de orde van \(h^2\); andere methoden, zoals die van Simpson, kunnen kleinere fouten opleveren over hetzelfde aantal subintervallen (mits de functie voldoende glad is).
3. Gevoelig voor intervalselectie
Als het interval te breed is en \(n\) klein, kunnen de resultaten aanzienlijk afwijken.
Sluitend
De trapeziummethode voor integratie is een van de meest fundamentele technieken in numerieke integratie. Door een interval in kleine delen te verdelen en de curve in elk deel te vervangen door een rechte lijn, kunnen we de integraal benaderen zonder de primitieve functie analytisch te hoeven bepalen. De formule is eenvoudig, gemakkelijk te implementeren en vooral nuttig wanneer de gegevens discreet zijn of de functie moeilijk te integreren is.
Gebruikers van de trapeziummethode moeten echter wel letten op het aantal subintervallen en het gedrag van de te integreren functie. Voor sterk gekromde functies of wanneer een hoge precisie vereist is, kan de trapeziummethode nog steeds worden gebruikt door \(n\) te verhogen, of kunnen andere numerieke methoden zoals de methode van Simpson of Romberg worden overwogen. Met een goed begrip wordt de trapeziummethode een effectief en praktisch hulpmiddel voor het oplossen van diverse integraalproblemen in de wetenschap, techniek en economie.
