Trapezregel

Die Trapezregel beschreibt ein mathematisches Verfahren, wie man das Integral einer Funktion im Intervall [a,b][a,b] numerisch annähert. Das entspricht der Fläche unter der Kurve f(x)f(x) bei kartesischer Darstellung.
Dazu ersetzt man die Fläche unter der Kurve durch ein Trapez, oder bei Stückelung des Intervalls durch mehrere Trapeze.
Man kann die Kurve f(x)f(x) näherungsweise durch eine Sehne zwischen den Funktionswerten an den Stellen aa und bb ersetzen. Dies führt zur Sehnentrapezformel. Man kann aber auch in der Mitte von [a,b][a,b] die Tangente an f(x)f(x) legen und erhält dann die Tangententrapezformel.

Beispiel

J(f)=∫0233x−1 dx J(f) = \int\limits_{0}^{2}3^{3x-1} \, \mathrm{d}x =(7289ln⁡(3))=73,6282396649… = \over{728 }{ {9 \ln(3)}} = 73{,}6282396649\dots
Nun gilt es mit Hilfe der Trapezformel als Näherungsverfahren dieses Integral zu bestimmen.

Sehnentrapezformel

Sehnentrapezformel.png
Sehnentrapez
J(f)=∫abf(x) dx=Q(f)+E(f) J(f) = \int\limits_{a}^{b}f(x)\, \mathrm{d}x = Q(f) + E(f)
Das Trapez wird gebildet aus der Grundlinie [a,b][a,b] (dem Intervall auf der xx-Achse), den senkrechten Geraden [a,f(a)][a,f(a)] und [b,f(b)][b,f(b)] sowie der Sehne als Verbindungsgerade zwischen f(a)f(a) und f(b)f(b). Diese Sehne ersetzt die Kurve f(x)f(x).
Die Sehnentrapezformel ergibt sich aus dem Flächeninhalt des beschriebenen Trapezes:
Q(f)=b−a2 (f(a)+f(b)) Q(f) = \dfrac{b-a}{2}\ (f(a)+f(b))
Diese Formel - und auch die folgenden - kann man herleiten aus der "Allgemeinen Quadraturformel für eine Teilfläche" (siehe Numerische Quadratur).
Ist f(x)f(x) zweimal stetig differenzierbar in [a,b][a,b], dann gilt für das Restglied E(f)E(f) folgende Abschätzung (siehe Numerische Quadratur):
∣E(f)∣≤((b−a)312)max⁡a≤x≤b∣f′′(x)∣ \left| E(f) \right| \le \over{(b-a)^3 }{ 12} \max_{a\le x \le b} {\left| f''(x) \right|}
Ist f(x)f(x) zusätzlich noch reellwertig, dann gilt mit einer Zwischenstelle ζ\zeta aus [a,b][a,b]:
E(f)=−((b−a)312)⋅f′′(ζ) E(f) = - \over{(b-a)^3 }{ 12} \cdot {{f''(\zeta)}}
Anschließend an das obige Beispiel:
Q(f)=((2−0)2)⋅(f(0)+f(2)) Q(f) = \over{(2-0) }{ 2} \cdot \left( f(0) + f(2) \right) =(7303)=243,33333… = \over{730 }{ 3} = 243,33333\dots

Zusammengesetzte Sehnentrapezformel

J(f)=∫abf(x) dx=Q(f)(n)+E(n)(f) J(f) = \int\limits_{a}^{b}f(x)\, \mathrm{d}x = Q(f)^{(n)} + E^{(n)}(f)
Um das Integral noch besser annähern zu können, unterteilt man das Intervall [a,b][a,b] in nn nebeneinanderliegende gleich große Teilintervalle der Länge hh (mit h=(b−a)/nh = (b - a) / n). In jedem Teilintervall wendet man die Sehnentrapezformel für die einzelnen Teilflächen an und addiert danach die entstandenen Näherungen. Damit erhält man die summierte (bzw. zusammengesetzte) Sehnentrapezformel:
Q(f)=h ((12)⋅f(a) + ∑i=1n−1f(a+i⋅h) + (12)⋅f(b)) Q(f)=h \ \left(\over{1}{ 2} \cdot f(a) \ + \ \sum\limits_{i=1}^{n-1} f(a + i \cdot h) \ + \ \over{1}{ 2} \cdot f(b) \right)
mit
h=((b−a)n) h = \over{(b - a) }{ n}
Die Fehlerabschätzung für das Restglied lautet:
∣E(f)∣≤((b−a)12)h2max⁡a≤x≤b∣f′′(x)∣\left| E(f) \right| \le \over{(b-a) }{ 12}h^2 \max_{a\le x \le b} {\left| f''(x) \right|}
bzw. für reellwertige Funktionen mit einer Zwischenstelle ζ\zeta aus dem Intervall [a,b][a,b]:
E(f)=−((b−a)12)⋅h2⋅f′′(ζ) E(f)= - \over{(b-a) }{ 12} \cdot h^2 \cdot f''(\zeta)
Anschließend an das obige Beispiel: Sei die Schrittweite h=(13)≥n=6h = \over{1 }{ 3} \ge n = 6
Q(6)(f)=(13)⋅((12) f(0)+f((13))+f((23))+f(1)+f((43))+f((53))+(12) f(2)) Q^{(6)}(f) = \over{1 }{ 3} \cdot \left( \over{1 }{ 2} \ f(0) + f\left(\over{1 }{ 3}\right) + f\left(\over{2 }{ 3}\right) + f(1) + f\left(\over{4 }{ 3}\right) + f\left(\over{5 }{ 3}\right) + \over{1 }{ 2} \ f(2) \right) =(7289)=80,88888… = \over{728 }{ 9} = 80,88888\dots

Tangententrapezformel

Tangententrapezformel.png
Tangententrapez
Die obere Seite des Trapezes wird hier gebildet, indem man in der Mitte des Intervalls [a,b][a,b] eine Tangente an f(x)f(x) legt. Die restlichen Seiten sind die Grundlinie [a,b][a,b] (das Intervall auf der xx-Achse) und die senkrechten Geraden an den Stellen aa und bb bis zur Tangente.
Die Tangententrapezformel ergibt sich aus dem Flächeninhalt des beschriebenen Trapezes:
Q(f)=(b−a) f((a+b2)) Q(f) = (b - a) \ f\left(\over{a + b }{ 2} \right)
Diese Formel - und auch die folgenden - kann man herleiten aus der "Allgemeinen Quadraturformel für eine Teilfläche" (siehe Numerische Quadratur).
Damit lässt sich das Integral darstellen als
J(f)=∫abf(x) dx=Q(f)+E(f) J(f) = \int\limits_{a}^{b}f(x)\, \mathrm{d}x = Q(f) + E(f)
Ist f(x)f(x) zweimal stetig differenzierbar in [a,b][a,b], dann gilt für das Restglied E(f)E(f) folgende Abschätzung (siehe Numerische Quadratur):
∣E(f)∣≤((b−a)324)max⁡a≤x≤b∣f′′(x)∣\left| E(f) \right| \le \over{(b-a)^3 }{ 24} \max_{a\le x \le b} {\left| f''(x) \right|}
Ist f(x)f(x) zusätzlich noch reellwertig, dann gilt mit einer Zwischenstelle ζ\zeta aus [a,b][a,b]:
E(f)=−((b−a)324)⋅f′′(ζ) E(f) = - \over{(b-a)^3 }{ 24} \cdot {{f''(\zeta)}}
Anschließend an das obige Beispiel:
 Q(f)=(2−0)⋅f(1)=18 \ Q(f) = (2-0) \cdot f(1) = 18

Zusammengesetzte Tangententrapezformel

Um das Integral noch besser annähern zu können unterteilt man das Intervall [a,b][a,b] in nn nebeneinanderliegende gleich große Teilintervalle der Länge hh. In jedem Teilintervall wendet man die Tangententrapezformel für die einzelnen Teilflächen an und addiert danach die entstandenen Näherungen. Damit erhält man die summierte Tangententrapezformel:
Q(f)=h⋅∑i=1nf(a + h⋅(2i−12)) Q(f)=h \cdot \sum\limits_{i=1}^n f(a \ + \ h \cdot \over{2i - 1 }{ 2})
mit
h=((b−a)n) h = \over{(b - a) }{ n}
Die Fehlerabschätzung für das Restglied lautet:
∣E(f)∣≤((b−a)24) h2max⁡a≤x≤b∣f′′(x)∣\left| E(f) \right| \le \over{(b - a) }{ 24} \ h^2 \max_{a\le x \le b} {\left| f''(x) \right|}
bzw. für reellwertige Funktionen mit einer Zwischenstelle ζ\zeta aus dem Intervall [a,b][a,b]:
E(f)=((b−a)24)⋅h2⋅f′′(ζ) E(f)=\over{(b - a) }{ 24} \cdot h^2 \cdot f''(\zeta)
Anschließend an das obige Beispiel: Sei die Schrittweite h=(13)≥n=6h = \over{1 }{ 3} \ge n = 6
Q(6)(f)=(13)⋅(f((16))+f((36))+f((56))+f((76))+f((96))+f((116))) Q^{(6)}(f) = \over{1 }{ 3} \cdot \left( f\left(\over{1 }{ 6}\right) + f\left(\over{3 }{ 6}\right) + f\left(\over{5 }{ 6}\right) + f\left(\over{7 }{ 6}\right) + f\left(\over{9 }{ 6}\right) + f\left(\over{11 }{ 6}\right) \right) =(36439) = \over{364 \sqrt 3 }{ 9} =70,05183266… = 70,05183266\dots

Siehe auch

 
 

Religion und Mathematik sind nur verschiedene Ausdrucksformen derselben göttlichen Exaktheit.

Kardinal Michael Faulhaber

Copyright- und Lizenzinformationen: Diese Seite basiert dem Artikel Trapezregel aus der frеiеn Enzyklοpädιe Wιkιpеdιa und stеht unter der Dοppellizеnz GNU-Lιzenz für freie Dokumentation und Crеative Commons CC-BY-SA 3.0 Unportеd (Kurzfassung). In der Wιkιpеdιa ist eine Listе dеr Autorеn des Originalartikels verfügbar. Da der Artikel geändert wurde, reicht die Angabe dieser Liste für eine lizenzkonforme Weiternutzung nicht aus!
Anbieterkеnnzeichnung: Mathеpеdιa von Тhοmas Stеιnfеld  • Dοrfplatz 25  •  17237 Blankеnsее  • Tel.: 01734332309 (Vodafone/D2)  •  Email: cο@maτhepedιa.dе
Datenschutzerklärung

 

Numerische Integration