[NumCS] Finish up adaptive quadrature

This commit is contained in:
2025-10-28 11:29:31 +01:00
parent 0a17657d2c
commit 15ce3a1d3b
5 changed files with 18 additions and 2 deletions

Binary file not shown.

View File

@@ -1,4 +1,17 @@
\newsection \newsection
\subsection{Adaptive Quadratur} \subsection{Adaptive Quadratur}
Der lokale Fehler einer zusammengesetzten Quadraturformel auf dem Gitter $\mathcal{M} := \{ a = x_0 < x_1 < \dots < x_m = b \}$ ist (für $f \in C^2([a, b])$). Der lokale Fehler einer zusammengesetzten Quadraturformel auf dem Gitter $\mathcal{M} := \{ a = x_0 < x_1 < \dots < x_m = b \}$ ist (für $f \in C^2([a, b])$):
Test \begin{align*}
\left| \int_{x_k}^{x_{k + 1}} f(t) \dx t - \frac{f(x_k) + f(x_{k + 1})}{2}(x_{k + 1} - x_k) \right| \leq (x_{k + 1} - x_k)^3 ||f''||_{L^\infty([x_k, x_{k + 1}])}
\end{align*}
Also ist es nur sinnvoll, das Gitter zu verfeinern wo $|f''|$ gross ist.
Auf Seiten 150 - 151 im Skript findet sich Code, um eine adaptive Quadratur durchzuführen.
\setLabelNumber{all}{3}
\fancyremark{Adaptive Quadratur in Python} Mit \texttt{scipy.integrate.quad} können wir einfach eine adaptive Quadratur durchführen und benutzt \texttt{QUADPACK}.
Mit \texttt{scipy.integrate.quadrature} können wir die Gauss-Quadratur verwenden.
Für $x \in \R^d$, also eine mehrdimensionale Funktion der Dimension $d$ können wir \texttt{scipy.integrate.nquad} verwenden. Mehr dazu im nächsten Kapitel
% TODO: Possibly explain the graphs and / or add code for computation using scipy

View File

@@ -1 +1,2 @@
\newsection
\subsection{Quadratur in $\R^d$ und dünne Gitter} \subsection{Quadratur in $\R^d$ und dünne Gitter}

View File

@@ -1 +1,2 @@
\newsection
\subsection{Monte-Carlo Quadratur} \subsection{Monte-Carlo Quadratur}

View File

@@ -1 +1,2 @@
\newsection
\subsection{Methoden zur Reduktion der Varianz} \subsection{Methoden zur Reduktion der Varianz}