4.1. Trapezfunktion#
In dieser Aufgabe sollen Sie die in der Vorlesung vorgestellte Trapezfunktion selbst programmieren.
Aufgabenstellung#
Aufgabenteil A#
Schreiben Sie eine Funktion, die die in der Vorlesung vorgestellte Trapezregel zum numerischen Integrieren von Funktionen implementiert. Die Funktion soll ein Numpy Array an Funktionswerten und eine Schrittweite \(\Delta x\) entgegen nehmen und daraus den Flächeninhalt bestimmen. Gehen Sie davon aus, dass für die zu den Funktionswerten gehörigen x-Werte gilt: \(x_0, x_0+\Delta x, x_0 + 2\Delta x, ...\).
Aufgabenteil B#
Wählen sie eine beliebige Funktion und überprüfen Sie die von Ihnen erstellte Funktion mit der Trapezfunktion des Scipy Moduls.
Lösungsvorschlag#
Show code cell content
import matplotlib.pyplot as plt
import numpy as np
import scipy as sp
from scipy import integrate
def fkt(xi):
yi = np.sin(xi)
return yi
def trapez(werte, schrittweite):
flaeche = 0
for i in range(len(werte)-1):
flaeche += werte[i] * schrittweite
flaeche += (werte[i+1] - werte[i]) * schrittweite / 2
return flaeche
xi = np.arange(0,2*np.pi,0.2)
yi = fkt(xi)
flaeche_sp = sp.integrate.trapz(yi, xi)
flaeche_selbst = trapez(yi, 0.2)
print('Eigene Trapezfunktion:', flaeche_selbst)
print('Scipy:', flaeche_sp)
Eigene Trapezfunktion: 0.0034463689419717627
Scipy: 0.0034463689419721374
/tmp/ipykernel_4902/2441456229.py:20: DeprecationWarning: 'scipy.integrate.trapz' is deprecated in favour of 'scipy.integrate.trapezoid' and will be removed in SciPy 1.14.0
flaeche_sp = sp.integrate.trapz(yi, xi)