Digital Garden
Maths
Calculus
Integrals
Numerical Integration

Numerical Integration

Die numerische Integration wird genutzt, wenn sich eine Stammfunktion nicht durch elementare Funktionen ausdrücken lässt, die numerische Auswertung der Stammfunktion zu komplex oder zeitaufwendig ist. In diesen Fällen verwendet man spezielle Näherungsverfahren, die sog. numerische Integration.

Trapezformel

Wir setzen voraus, dass die stetige Funktion f(x)f(x) im Integrationsintervall axba\leq x \leq b oberhalb der x-Achse verläuft damit wir es uns besser vorstellen können jedoch muss dies nicht der Fall sein um die Trapezformel zu verwenden.

Wir zerlegen zuerst das Integrationsintervall in nn Teilintervalle gleicher Länge

h=banh = \frac{b-a}{n}

Die Punkte der Teilintervalle auf der x-Achse nennen wir Stützstellen.

x0=a,x1=x0+h,...,xk=x0+kh,...,xn=bx_0=a,\,x_1=x_0+h,...,x_k=x_0+k\cdot h,...,x_n=b

Die Funktionswerte der Stützstellen nennen wir Stützwerte.

yk=f(xk)=f(x0+kh)=f(a+kh) mit k=0,1,...,ny_k=f(x_k)=f(x_0+k\cdot h)=f(a+k\cdot h)\text{ mit } k=0,1,...,n

Mit der Formel für den Flächeninhalt eines geometrischen Trapez A=a+b2hA=\frac{a+b}{2}h können wir die Fläche des ersten Trapezes berechnen A1=y0+y12hA_1=\frac{y_0+y_1}{2}h und genau so dann die restlichen Flächen der Trapeze

An=yn1+yn2hA_n=\frac{y_{n-1}+y_n}{2}h

trapezFlaecheFormel trapezFormel

Die Summe aller Trapezflächen geben dann eine gute Näherung für den gesuchten Flächenwert. Welches zu einer ziemlich kurzen Formel zusammen gezogen werden kann indem man die Brüche gleichnamig macht und ein wenig ausklammert und faktorisiert.

abf(x)dxA1+...+An=y0+y12h+y1+y22h+...+yn1+yn2h=(y0+2y1+2y2+...+2yn1+yn)h2=(121+2)h\begin{align*} \int_{a}^{b}{f(x)\,dx}&\approx A_1+...+A_n \\ &= \frac{y_0+y_1}{2}h + \frac{y_1+y_2}{2}h +... + \frac{y_{n-1}+y_n}{2}h \\ &= (y_0+2y_1+2y_2+...+2y_{n-1}+y_n)\frac{h}{2} \\ &= (\frac{1}{2}\sum{}_1 + \sum{}_2)h \end{align*}

wobei 1=\sum{}_1= Die Summe der beiden äusseren Stützwerte (y1,yn)(y_1,y_n) und 2=\sum{}_2= Die Summe der inneren Stützwerte (y2,...yn1)(y_2,...y_{n-1}).

Wenn nn \to \infty dann liefert die Trapezformel den exakten Integralwert.

Mehr dazu kann man auch hier (opens in a new tab) finden.

Simpsonsche Formel

link (opens in a new tab) link (opens in a new tab)