Trapezregel

Trapezregel

Studium
∫abf(x) dx≈h2(f(x0)+2f(x1)+…+2f(xn−1)+f(xn))\int_a^b f(x)\,dx \approx \frac{h}{2} \left( f(x_0) + 2 f(x_1) + \ldots + 2 f(x_{n-1}) + f(x_n) \right)

Voraussetzung: h=b−an,h = \frac{b - a}{n}, xk=a+kh,x_k = a + k h, n∈N,n \in \mathbb{N}, n≥1n \geq 1

a,ba, b
untere und obere Integrationsgrenze
ff
Funktion, die integriert wird
nn
Anzahl der Teilintervalle
hh
Schrittweite, Breite eines Teilintervalls
xkx_k
Stützstellen, x₀ = a bis xₙ = b

Die Trapezregel nähert ein Integral durch Trapeze unter den Sehnen zwischen benachbarten Stützstellen an.

Erklärung

Die Trapezregel berechnet ein bestimmtes Integral näherungsweise, wenn keine Stammfunktion bekannt ist oder ff nur an einzelnen Stellen gegeben ist. Man teilt [a,b][a, b] in nn gleich breite Teilintervalle der Breite hh und ersetzt den Graphen in jedem Teilintervall durch die Sehne. Die Fläche darunter ist ein Trapez mit dem Inhalt h2(f(xk−1)+f(xk))\frac{h}{2}\left(f(x_{k-1}) + f(x_k)\right).

Addiert man alle Trapeze, kommen die inneren Stützstellen zweimal vor, deshalb der Faktor 22. Für n=1n = 1 ergibt sich die einfache Trapezregel b−a2(f(a)+f(b))\frac{b - a}{2}\left(f(a) + f(b)\right).

Ist ff zweimal stetig differenzierbar, gilt für den Fehler

∫abf(x) dx−T(h)=−(b−a) h212 f′′(ξ)\int_a^b f(x)\,dx - T(h) = -\frac{(b - a) \, h^2}{12} \, f''(\xi)

für ein ξ∈[a,b]\xi \in [a, b]. Halbiert man hh, wird der Fehler etwa viermal kleiner (Ordnung 2). Lineare Funktionen werden exakt integriert. Bei konvexen Funktionen (f′′>0f'' > 0) liegt die Sehne über dem Graphen, die Trapezregel überschätzt dann das Integral.

Beispiele

∫01x2 dx=13\int_0^1 x^2\,dx = \frac{1}{3} mit n=4n = 4, also h=0,25h = 0{,}25:

T=0,252(0+2⋅(0,0625+0,25+0,5625)+1)=0,34375T = \frac{0{,}25}{2} \left( 0 + 2 \cdot (0{,}0625 + 0{,}25 + 0{,}5625) + 1 \right) = 0{,}34375

Der Fehler ist 0,34375−13=196≈0,01040{,}34375 - \frac{1}{3} = \frac{1}{96} \approx 0{,}0104. Wegen f′′=2f'' = 2 stimmt das genau mit der Fehlerformel 1⋅0,25212⋅2\frac{1 \cdot 0{,}25^2}{12} \cdot 2 überein.

∫121x dx=ln⁡2≈0,6931\int_1^2 \frac{1}{x}\,dx = \ln 2 \approx 0{,}6931. Mit n=2n = 2, h=0,5h = 0{,}5:

T=0,52(1+2⋅11,5+12)=1724≈0,7083T = \frac{0{,}5}{2} \left( 1 + 2 \cdot \frac{1}{1{,}5} + \frac{1}{2} \right) = \frac{17}{24} \approx 0{,}7083

Mit n=4n = 4 ergibt sich 0,69700{,}6970. Der Fehler sinkt von etwa 0,01520{,}0152 auf 0,00390{,}0039, also ungefähr auf ein Viertel.

Typische Fehler

  • Die inneren Funktionswerte nicht verdoppeln oder die Randwerte verdoppeln.
  • hh und nn verwechseln: Bei nn Teilintervallen gibt es n+1n + 1 Stützstellen.
  • Den Faktor 12\frac{1}{2} vor der Klammer vergessen.
  • Bei unstetigen oder stark schwankenden Funktionen der Näherung trauen, ohne nn zu erhöhen und die Ergebnisse zu vergleichen.

Verwandte Formeln