diff options
-rw-r--r-- | buch/chapters/010-potenzen/polynome.tex | 239 | ||||
-rw-r--r-- | buch/chapters/030-geometrie/chapter.tex | 3 | ||||
-rw-r--r-- | buch/chapters/030-geometrie/images/einheitskreis.pdf | bin | 19706 -> 20005 bytes | |||
-rw-r--r-- | buch/chapters/030-geometrie/images/einheitskreis.tex | 4 | ||||
-rw-r--r-- | buch/chapters/030-geometrie/uebungsaufgaben/3.tex | 169 | ||||
-rw-r--r-- | buch/papers/dreieck/main.tex | 20 | ||||
-rw-r--r-- | buch/papers/dreieck/references.bib | 36 | ||||
-rw-r--r-- | buch/papers/dreieck/teil0.tex | 45 | ||||
-rw-r--r-- | buch/papers/dreieck/teil1.tex | 89 | ||||
-rw-r--r-- | buch/papers/dreieck/teil2.tex | 112 | ||||
-rw-r--r-- | buch/papers/dreieck/teil3.tex | 73 | ||||
-rw-r--r-- | vorlesungen/slides/hermite/hermiteentwicklung.tex | 15 | ||||
-rw-r--r-- | vorlesungen/slides/hermite/loesung.tex | 19 | ||||
-rw-r--r-- | vorlesungen/slides/hermite/normalhermite.tex | 29 | ||||
-rw-r--r-- | vorlesungen/slides/hermite/normalintegrale.tex | 9 | ||||
-rw-r--r-- | vorlesungen/slides/hermite/skalarprodukt.tex | 22 |
16 files changed, 793 insertions, 91 deletions
diff --git a/buch/chapters/010-potenzen/polynome.tex b/buch/chapters/010-potenzen/polynome.tex index 5f119e5..981e444 100644 --- a/buch/chapters/010-potenzen/polynome.tex +++ b/buch/chapters/010-potenzen/polynome.tex @@ -13,20 +13,30 @@ Operationen konstruieren lassen, sind die Polynome. \index{Polynom}% Ein {\em Polynome} vom Grad $n$ ist die Funktion \[ -p(x) = a_nx^2n + a_{n-1}x^{n-1} + \dots + a_2x^2 + a_1x + a_0, +p(x) = a_nx^n + a_{n-1}x^{n-1} + \dots + a_2x^2 + a_1x + a_0, \] wobei $a_n\ne 0$ sein muss. Das Polynom heisst {\em normiert}, wenn $a_n=1$ ist. \index{normiert}% +\index{Grad eines Polynoms}% Die Menge aller Polynome mit Koeffizienten in der Menge $K$ wird mit $K[x]$ bezeichnet. \end{definition} Die Menge $K[x]$ ist heisst auch der {\em Polynomring}, weil $K[x]$ -mit der Addition, Subtraktion und Multiplikation von Polynomen ein -Ring mit $1$ ist. -Im Folgenden werden wir uns auf die Fälle $K=\mathbb{R}$ und $K=\mathbb{C}$ -beschränken. +mit der Addition, Subtraktion und Multiplikation von Polynomen eine +algebraische Struktur bildet, die man einen Ring mit $1$ nennt. +\index{Ring}% +Im Folgenden werden wir uns auf die Fälle $K=\mathbb{Q}$, $K=\mathbb{R}$ +und $K=\mathbb{C}$ beschränken. + +Für den Grad eines Polynoms gelten die bekannten Rechenregeln +\begin{align*} +\deg (a(x) + b(x)) &\le \operatorname{max}(\deg a(x), \deg b(x)) +\\ +\deg (a(x)\cdot b(x)) &=\deg a(x) + \deg b(x) +\end{align*} +für beliebige Polynome $a(x),b(x)\in K[x]$. In Abschnitt~\ref{buch:orthogonalitaet:section:orthogonale-funktionen} werden Familien von Polynomen konstruiert werden, die sich durch eine @@ -35,12 +45,14 @@ Diese Polynome lassen sich typischerweise auch als Lösungen von Differentialgleichungen finden. Ausserdem werden hypergeometrische Funktionen \[ -\mathstrut_pF_q\biggl(\begin{matrix}a_1,\dots,a_p\\b_1,\dots,b_q\end{matrix};z\biggr), +\mathstrut_pF_q\biggl( +\begin{matrix}a_1,\dots,a_p\\b_1,\dots,b_q\end{matrix};z +\biggr), \] die in Abschnitt~\ref{buch:rekursion:section:hypergeometrische-funktion} definiert werden, zu Polynomen, wenn mindestens einer der Parameter $a_k$ negativ ganzzahlig ist. -Polynome sind also bereits eine Vielfältige Quelle von speziellen +Polynome sind also bereits eine vielfältige Quelle von speziellen Funktionen. Viele spezielle Funktionen werden aber komplizierter sein und @@ -53,6 +65,7 @@ Dank des folgenden Satzes kann dies immer mit Polynomen geschehen. \begin{satz}[Weierstrass] \label{buch:potenzen:satz:weierstrass} +\index{Weierstrass, Satz von}% Eine auf einem kompakten Intervall $[a,b]$ stetige Funktion $f(x)$ lässt sich durch eine Folge $p_n(x)$ von Polynomen gleichmässig approximieren. @@ -69,6 +82,189 @@ ebenfalls als Approximationen dienen können. Weitere Möglichkeiten liefern Interpolationsmethoden der numerischen Mathematik. +\subsection{Polynomdivision, Teilbarkeit und grösster gemeinsamer Teiler} +Der schriftliche Divisionsalgorithmus für Zahlen funktioniert +auch für die Division von Polynomen. +Zu zwei beliebigen Polynomen $p(x)$ und $q(x)$ lassen sich also +immer zwei Polynome $a(x)$ und $r(x)$ finden derart, dass +$p(x) = a(x) q(x) + r(x)$. +Das Polynom $a(x)$ heisst der {\em Quotient}, $r(x)$ der {\em Rest} +der Division. +Das Polynom $p(x)$ heisst {\em teilbar} durch $q(x)$, geschrieben +$q(x)\mid p(x)$, wenn $r(x)=0$ ist. + +\subsubsection{Grösster gemeinsamer Teiler} +Mit dem Begriff der Teilbarkeit geht auch die Idee des grössten +gemeinsamen Teilers einher. +Ein gemeinsamer Teiler zweier Polynome $a(x)$ und $b(x)$ +ist ein Polynom $g(x)$, welches beide Polynome teilt, also +$g(x)\mid a(x)$ und $g(x)\mid b(x)$. +\index{grösster gemeinsamer Teiler}% +Ein Polynome $g(x)$ heisst grösster gemeinsamer Teiler von $a(x)$ +und $b(x)$, wenn jeder andere gemeinsame Teiler $f(x)$ von $a(x)$ +und $b(x)$ auch ein Teiler von $g(x)$ ist. +Man beachte, dass die skalaren Vielfachen eines grössten gemeinsamen +Teilers ebenfalls grösste gemeinsame Teiler sind, der grösste gemeinsame +Teiler ist also nicht eindeutig bestimmt. + +\subsubsection{Der euklidische Algorithmus} +Zur Berechnung eines grössten gemeinsamen Teilers steht wie bei den +ganzen Zahlen der euklidische Algorithmus zur Verfügung. +Dazu bildet man die Folgen von Polynomen +\[ +\begin{aligned} +a_0(x)&=a(x) & b_0(x) &= b(x) +& +&\Rightarrow& +a_0(x)&=b_0(x) q_0(x) + r_0(x) && +\\ +a_1(x)&=b_0(x) & b_1(x) &= r_0(x) +& +&\Rightarrow& +a_1(x)&=b_1(x) q_1(x) + r_1(x) && +\\ +a_2(x)&=b_1(x) & b_2(x) &= r_1(x) +& +&\Rightarrow& +a_2(x)&=b_2(x) q_2(x) + r_2(x) && +\\ +&&&&&\hspace*{2mm}\vdots&& +\\ +a_{m-1}(x)&=b_{m-2}(x) & b_{m-1}(x) &= r_{m-2}(x) +& +&\Rightarrow& +a_{m-1}(x)&=b_{m-1}(x)q_{m-1}(x) + r_{m-1}(x) &\text{mit }r_{m-1}(x)&\ne 0 +\\ +a_m(x)&=b_{m-1}(x) & b_m(x)&=r_{m-1}(x) +& +&\Rightarrow& +a_m(x)&=b_m(x)q_m(x).&& +\end{aligned} +\] +Der Index $m$ ist der Index, bei dem zum ersten Mal $r_m(x)=0$ ist. +Dann ist $g(x)=r_{m-1}(x)$ ein grösster gemeinsamer Teiler. + +\subsubsection{Der erweiterte euklidische Algorithmus} +Die Konstruktion der Folgen $a_n(x)$ und $b_n(x)$ kann in Matrixform +kompakter geschrieben werden als +\[ +\begin{pmatrix} +a_k(x)\\ +b_k(x) +\end{pmatrix} += +\begin{pmatrix} +b_{k-1}(x)\\ +r_{k-1}(x) +\end{pmatrix} += +\begin{pmatrix} +0 & 1\\ +1 & -q_{k-1}(x) +\end{pmatrix} +\begin{pmatrix} +a_{k-1}(x)\\ +b_{k-1}(x) +\end{pmatrix}. +\] +Kürzen wir die $2\times 2$-Matrix als +\[ +Q_k(x) = \begin{pmatrix} 0&1\\1&-q_k(x)\end{pmatrix} +\] +ab, dann ergibt das Produkt der Matrizen $Q_0(x)$ bis $Q_{m}(x)$ +\[ +\begin{pmatrix} +g(x)\\ +0 +\end{pmatrix} += +\begin{pmatrix} +r_{m-1}(x)\\ +r_{m}(x) +\end{pmatrix} += +Q_{m}(x) +Q_{m-1}(x) +\cdots +Q_1(x) +Q_0(x) +\begin{pmatrix} +a(x)\\ +b(x) +\end{pmatrix}. +\] +Zur Berechnung des Produktes der Matrizen $Q_k(x)$ kann man rekursiv +vorgehen mit der Rekursionsformel +\[ +S_{k}(x) = Q_{k}(x) S_{k-1}(x) +\qquad\text{mit}\qquad +S_{-1}(x) += +\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}. +\] +Ausgeschrieben bedeutet dies Matrixrekursionsformel +\[ +S_{k-1}(x) += +\begin{pmatrix} +c_{k-1} & d_{k-1} \\ +c_k & d_k +\end{pmatrix} +\qquad\Rightarrow\qquad +Q_{k}(x) S_{k-1}(x) += +\begin{pmatrix} +0&1\\1&-q_k(x) +\end{pmatrix} +\begin{pmatrix} +c_{k-1} & d_{k-1} \\ +c_k & d_k +\end{pmatrix} += +\begin{pmatrix} +c_k&d_k\\ +c_{k+1}&d_{k+1} +\end{pmatrix}. +\] +Daraus lässt sich für die Matrixelemente die Rekursionsformel +\[ +\begin{aligned} +c_{k+1} &= c_{k-1} - q_k(x) c_k(x) \\ +d_{k+1} &= d_{k-1} - q_k(x) d_k(x) +\end{aligned} +\quad +\bigg\} +\qquad +\text{mit Startwerten} +\qquad +\bigg\{ +\begin{aligned} +\quad +c_{-1} &= 1, & c_0 &= 0 \\ +d_{-1} &= 0, & d_0 &= 1. +\end{aligned} +\] +Wendet man die Matrix $S_m(x)$ auf den Vektor mit den Komponenten +$a(x)$ und $b(x)$, erhält man die Beziehungen +\[ +g(x) = c_{k-1}(x) a(x) + d_{k-1}(x) b(x) +\qquad\text{und}\qquad +0 = c_k(x) a(x) + d_k(x) b(x). +\] +Dieser Algorithmus heisst der erweiterte euklidische Algorithmus. +Wir fassen die Resultate zusammen im folgenden Satz. + +\begin{satz} +Zu zwei Polynomen $a(x),b(x) \in K[x]$ gibt es Polynome +$g(x),c(x),d(x)\in K[x]$ +derart, dass $g(x)$ ein grösster gemeinsamer Teiler von $a(x)$ und $b(x)$ +ist, und ausserdem +\[ +g(x) = c(x)a(x)+d(x)b(x) +\] +gilt. +\end{satz} + \subsection{Faktorisierung und Nullstellen} % wird später gebraucht um bei der Definition der hypergeometrischen Reihe % die Zaehler- und Nenner-Polynome als Pochhammer-Symbole zu entwickeln @@ -77,11 +273,24 @@ numerischen Mathematik. % Wird gebraucht für die Potenzreihen-Methode % Muss später ausgedehnt werden auf Potenzreihen -\subsection{Polynom-Berechnung} -Die naive Berechnung der Werte eines Polynoms beginnt mit der Berechnung -der Potenzen. -Die Anzahl nötiger Multiplikationen kann minimiert werden, indem man -das Polynom als +\subsection{Berechnung von Polynom-Werten} +Die naive Berechnung der Werte eines Polynoms $p(x)$ vom Grad $n$ +beginnt mit der Berechnung der Potenzen von $x$. +Da alle Potenzen benötigt werden, wird man dazu $n-1$ Multiplikationen +benötigen. +Die Potenzen müssen anschliessend mit den Koeffizienten multipliziert +werden, dazu sind weitere $n$ Multiplikationen nötig. +Der Wert des Polynoms kann also erhalten werden mit $2n-1$ Multiplikationen +und $n$ Additionen. + +Die Anzahl nötiger Multiplikationen kann mit dem folgenden Vorgehen +reduziert werden, welches auch als das {\em Horner-Schema} bekannt ist. +\index{Horner-Schema}% +Statt erst am Schluss alle Terme zu addieren, addiert man so früh +wie möglich. +Zum Beispiel multipliziert man $(a_nx+a_{n-1})$ mit $x$, was auf +die Multiplikationen beider Terme mit $x$ hinausläuft. +Mit dieser Idee kann man das Polynom als \[ a_nx^n + @@ -95,10 +304,10 @@ a_0 = ((\dots((a_nx+a_{n-1})x+a_{n-2})x+\dots )x+a_1)x+a_0 \] -schreibt. +schreiben. Beginnend bei der innersten Klammer sind genau $n$ Multiplikationen -und $n+1$ Additionen nötig, im Gegensatz zu $2n$ Multiplikationen -und $n$ Additionen bei der naiven Vorgehensweise. +und $n$ Additionen nötig, man spart mit diesem Vorgehen also +$n-1$ Multiplikationen. diff --git a/buch/chapters/030-geometrie/chapter.tex b/buch/chapters/030-geometrie/chapter.tex index f3f1d39..0b2842b 100644 --- a/buch/chapters/030-geometrie/chapter.tex +++ b/buch/chapters/030-geometrie/chapter.tex @@ -42,7 +42,7 @@ wie die Berechnung der Länge von Ellipsen- oder Hyperbelbögen auf die Notwendigkeit führt, neue spezielle Funktionen zu definieren. \input{chapters/030-geometrie/trigonometrisch.tex} -\input{chapters/030-geometrie/sphaerisch.tex} +%\input{chapters/030-geometrie/sphaerisch.tex} \input{chapters/030-geometrie/hyperbolisch.tex} \input{chapters/030-geometrie/laenge.tex} \input{chapters/030-geometrie/flaeche.tex} @@ -54,5 +54,6 @@ die Notwendigkeit führt, neue spezielle Funktionen zu definieren. %\uebungsaufgabe{0} \uebungsaufgabe{1} \uebungsaufgabe{2} +\uebungsaufgabe{3} \end{uebungsaufgaben} diff --git a/buch/chapters/030-geometrie/images/einheitskreis.pdf b/buch/chapters/030-geometrie/images/einheitskreis.pdf Binary files differindex 0b514eb..d708377 100644 --- a/buch/chapters/030-geometrie/images/einheitskreis.pdf +++ b/buch/chapters/030-geometrie/images/einheitskreis.pdf diff --git a/buch/chapters/030-geometrie/images/einheitskreis.tex b/buch/chapters/030-geometrie/images/einheitskreis.tex index c38dc19..a194190 100644 --- a/buch/chapters/030-geometrie/images/einheitskreis.tex +++ b/buch/chapters/030-geometrie/images/einheitskreis.tex @@ -41,6 +41,7 @@ \fill[color=blue] (\a:\r) circle[radius=0.05]; \draw[color=blue,line width=1.4pt] (\r,0) -- (\r,{\r*tan(\a)}); +\fill[color=blue] (\r,{\r*tan(\a)}) circle[radius=1.0pt]; \node[color=blue] at (\r,{0.5*\r*tan(\a)}) [right] {$\tan\alpha$}; \draw[color=blue,line width=0.4pt] ({\r*cos(\a)},0) -- (\a:\r); @@ -53,6 +54,7 @@ \draw[color=blue] (-0.1,{\r*sin(\a)}) -- (0.1,{\r*sin(\a)}); \draw[color=blue,line width=1.4pt] (0,\r) -- ({\r/tan(\a)},\r); +\fill[color=blue] ({\r/tan(\a)},\r) circle[radius=1.0pt]; \node[color=blue] at ({0.5*\r/tan(\a)},\r) [above] {$\cot\alpha$}; \draw[color=darkgreen,line width=1pt] (0,0) -- (\b:\r); @@ -61,9 +63,11 @@ \fill[color=darkgreen] (\b:\r) circle[radius=0.05]; \draw[color=darkgreen,line width=1.4pt] (0,\r) -- ({\r/tan(\b)},\r); +\fill[color=darkgreen] ({\r/tan(\b)},\r) circle[radius=1.0pt]; \node[color=darkgreen] at ({0.5*\r/tan(\b)},\r) [above] {$\cot\beta$}; \draw[color=darkgreen,line width=1.4pt] (\r,0) -- (\r,{\r*tan(\b)}); +\fill[color=darkgreen] (\r,{\r*tan(\b)}) circle[radius=1.0pt]; \node[color=darkgreen] at (\r,{0.5*\r*tan(\b)}) [right] {$\tan\beta$}; \draw[color=darkgreen,line width=0.4pt] (\b:\r) -- (0,{\r*sin(\b)}); diff --git a/buch/chapters/030-geometrie/uebungsaufgaben/3.tex b/buch/chapters/030-geometrie/uebungsaufgaben/3.tex new file mode 100644 index 0000000..6a501fb --- /dev/null +++ b/buch/chapters/030-geometrie/uebungsaufgaben/3.tex @@ -0,0 +1,169 @@ +\def\cas{\operatorname{cas}} +Die Funktion $\cas$ definiert durch +$\cas x = \cos x + \sin x$ hat einige interessante Eigenschaften. +Wie die gewöhnlichen trigonometrischen Funktionen $\sin x$ und $\cos x$ +ist $\cas x$ $2\pi$-periodisch. +Die Ableitung und das Additionstheorem benötigen bei den gewöhnlichen +trigonometrischen Funktionen aber beide Funktionen, im Gegensatz zu den +im folgenden hergeleiteten Formeln, die nur die Funktion $\cas x$ brauchen. +\begin{teilaufgaben} +\item +Drücken Sie die Ableitung von $\cas x$ allein durch Werte der +$\cas$-Funktion aus. +\item +Zeigen Sie, dass +\[ +\cas x += +\sqrt{2} \sin\biggl(x+\frac{\pi}4\biggr) += +\sqrt{2} \cos\biggl(x-\frac{\pi}4\biggr). +\] +\item +Beweisen Sie das Additionstheorem für die $\cas$-Funktion +\begin{equation} +\cas(x+y) += +\frac12\bigl( +\cas(x)\cas(y) + \cas x\cas (-y) + \cas(-x)\cas(y) -\cas(-x)\cas(-y) +\bigr) +\label{buch:geometrie:uebung3:eqn:addition} +\end{equation} +\end{teilaufgaben} +Youtuber Dr Barker hat die Funktion $\cas$ im Video +{\small\url{https://www.youtube.com/watch?v=bn38o3u0lDc}} vorgestellt. + +\begin{loesung} +\begin{teilaufgaben} +\item +Die Ableitung ist +\[ +\frac{d}{dx}\cas x += +\frac{d}{dx}(\cos x + \sin x) += +-\sin x + \cos x += +\sin(-x) + \cos(-x) += +\cas(x). +\] +\item +Die Additionstheoreme angewendet auf die trigonometrischen Funktionen +auf der rechten Seite ergibt +\begin{align*} +\sin\biggl(x+\frac{\pi}4\biggr) +&= +\sin x \cos\frac{\pi}4 + \cos x \sin\frac{\pi}4 +&&& +\cos\biggl(x-\frac{\pi}4\biggr) +&= +\cos(x)\cos\frac{\pi}4 -\sin x \sin\biggl(-\frac{\pi}4\biggr) +\\ +&= +\frac{1}{\sqrt{2}} \sin x ++ +\frac{1}{\sqrt{2}} \cos x +&&& +&= +\frac{1}{\sqrt{2}} \cos x ++ +\frac{1}{\sqrt{2}} \sin x +\\ +&=\frac{1}{\sqrt{2}} \cas x +&&& +&= +\frac{1}{\sqrt{2}} \cas x. +\end{align*} +Multiplikation mit $\sqrt{2}$ ergibt die behaupteten Relationen. +\item +Substituiert man die Definition von $\cas(x)$ auf der rechten Seite von +\eqref{buch:geometrie:uebung3:eqn:addition} und multipliziert aus, +erhält man +\begin{align*} +\eqref{buch:geometrie:uebung3:eqn:addition} +&= +{\textstyle\frac12}\bigl( +(\cos x + \sin x) +(\cos y + \sin y) ++ +(\cos x + \sin x) +(\cos y - \sin y) +\\ +&\qquad ++ +(\cos x - \sin x) +(\cos y + \sin y) +- +(\cos x - \sin x) +(\cos y - \sin y) +\bigr) +\\ +&= +\phantom{-\mathstrut} +{\textstyle\frac12}\bigl( +\cos x\cos y ++ +\cos x\sin y ++ +\sin x\cos y ++ +\sin x\sin y +\\ +& +\phantom{=-\mathstrut{\textstyle\frac12}\bigl(}\llap{$\mathstrut +\mathstrut$} +\cos x\cos y +- +\cos x\sin y ++ +\sin x\cos y +- +\sin x\sin y +\\ +& +\phantom{=-\mathstrut{\textstyle\frac12}\bigl(}\llap{$\mathstrut +\mathstrut$} +\cos x\cos y ++ +\cos x\sin y +- +\sin x\cos y +- +\sin x\sin y +\bigr) +\\ +& +\phantom{=} +-\mathstrut{\textstyle\frac12}\bigl( +\cos x\cos y +- +\cos x\sin y +- +\sin x\cos y ++ +\sin x\sin y +\bigr) +\\ +&= \cos x \cos y ++ +\cos x \sin y ++ +\sin x \cos y +- +\sin x \sin y. +\intertext{Die äussersten zwei Terme passen zum Additionstheorem für den +Kosinus, die beiden inneren Terme dagegen zum Sinus. +Fasst man sie zusammen, erhält man} +&= +(\sin x\cos y + \cos x \sin y) ++ +(\cos x\cos y - \sin x \sin y) +\\ +&= +\sin (x+y) + \cos(x+y) += +\cas(x+y). +\end{align*} +Damit ist das Additionstheorem für die Funktion $\cas$ bewiesen. +\qedhere +\end{teilaufgaben} +\end{loesung} diff --git a/buch/papers/dreieck/main.tex b/buch/papers/dreieck/main.tex index 75ba410..d7bc769 100644 --- a/buch/papers/dreieck/main.tex +++ b/buch/papers/dreieck/main.tex @@ -3,19 +3,21 @@ % % (c) 2020 Hochschule Rapperswil % -\chapter{Dreieckstest und Beta-Funktion\label{chapter:dreieck}} -\lhead{Dreieckstest und Beta-Funktion} +\chapter{$\int P(t) e^{-t^2} \,dt$ in geschlossener Form? +\label{chapter:dreieck}} +\lhead{Integrierbarkeit in geschlossener Form} \begin{refsection} \chapterauthor{Andreas Müller} \noindent -Mit dem Dreieckstest kann man feststellen, wie gut ein Geruchs- -oder Geschmackstester verschiedene Gerüche oder Geschmäcker -unterscheiden kann. -Seine wahrscheinlichkeitstheoretische Erklärung benötigt die Beta-Funktion, -man kann die Beta-Funktion als durchaus als die mathematische Grundlage -der Weindegustation -bezeichnen. +Der Risch-Algorithmus erlaubt, eine definitive Antwort darauf zu geben, +\index{Risch-Algorithmus}% +\index{elementare Stammfunktion}% +ob eine elementare Funktion eine Stammfunktion in geschlossener Form hat. +Der Algorithmus ist jedoch ziemlich kompliziert. +In diesem Kapitel soll ein spezieller Fall mit Hilfe der Theorie der +orthogonale Polynome, speziell der Hermite-Polynome, behandelt werden, +wie er in der Arbeit \cite{dreieck:polint} untersucht wurde. \input{papers/dreieck/teil0.tex} \input{papers/dreieck/teil1.tex} diff --git a/buch/papers/dreieck/references.bib b/buch/papers/dreieck/references.bib index d2bbe08..47bd865 100644 --- a/buch/papers/dreieck/references.bib +++ b/buch/papers/dreieck/references.bib @@ -4,32 +4,12 @@ % (c) 2020 Autor, Hochschule Rapperswil % -@online{dreieck:bibtex, - title = {BibTeX}, - url = {https://de.wikipedia.org/wiki/BibTeX}, - date = {2020-02-06}, - year = {2020}, - month = {2}, - day = {6} +@article{dreieck:polint, + author = { George Stoica }, + title = { Polynomials and Integration in Finite Terms }, + journal = { Amer. Math. Monthly }, + volume = 129, + year = 2022, + number = 1, + pages = {80--81} } - -@book{dreieck:numerical-analysis, - title = {Numerical Analysis}, - author = {David Kincaid and Ward Cheney}, - publisher = {American Mathematical Society}, - year = {2002}, - isbn = {978-8-8218-4788-6}, - inseries = {Pure and applied undegraduate texts}, - volume = {2} -} - -@article{dreieck:mendezmueller, - author = { Tabea Méndez and Andreas Müller }, - title = { Noncommutative harmonic analysis and image registration }, - journal = { Appl. Comput. Harmon. Anal.}, - year = 2019, - volume = 47, - pages = {607--627}, - url = {https://doi.org/10.1016/j.acha.2017.11.004} -} - diff --git a/buch/papers/dreieck/teil0.tex b/buch/papers/dreieck/teil0.tex index bcf2cf8..65eff7a 100644 --- a/buch/papers/dreieck/teil0.tex +++ b/buch/papers/dreieck/teil0.tex @@ -3,7 +3,48 @@ % % (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil % -\section{Testprinzip\label{dreieck:section:testprinzip}} -\rhead{Testprinzip} +\section{Problemstellung\label{dreieck:section:problemstellung}} +\rhead{Problemstellung} +Es ist bekannt, dass das Fehlerintegral +\[ +\frac{1}{\sqrt{2\pi}\sigma} \int_{-\infty}^x e^{-\frac{t^2}{2\sigma}}\,dt +\] +nicht in geschlossener Form dargestellt werden kann. +Mit der in Kapitel~\ref{buch:chapter:integral} skizzierten Theorie von +Liouville und dem Risch-Algorithmus kann dies strengt gezeigt werden. +Andererseits gibt es durchaus Integranden, die $e^{-t^2}$ enthalten, +für die eine Stammfunktion in geschlossener Form gefunden werden kann. +Zum Beispiel folgt aus der Ableitung +\[ +\frac{d}{dt} e^{-t^2} += +-2te^{-t^2} +\] +die Stammfunktion +\[ +\int te^{-t^2}\,dt += +-\frac12 e^{-t^2}. +\] +Leitet man $e^{-t^2}$ zweimal ab, erhält man +\[ +\frac{d^2}{dt^2} e^{-t^2} += +(4t^2-2) e^{-t^2} +\qquad\Rightarrow\qquad +\int (t^2-{\textstyle\frac12}) e^{-t^2}\,dt += +{\textstyle\frac14} +e^{-t^2}. +\] +Es gibt also eine viele weitere Polynome $P(t)$, für die der Integrand +$P(t)e^{-t^2}$ eine Stammfunktion in geschlossener Form hat. +Damit stellt sich jetzt das folgende allgemeine Problem. + +\begin{problem} +\label{dreieck:problem} +Für welche Polynome $P(t)$ hat der Integrand $P(t)e^{-t^2}$ +eine elementare Stammfunktion? +\end{problem} diff --git a/buch/papers/dreieck/teil1.tex b/buch/papers/dreieck/teil1.tex index 4abe2e1..45c1a23 100644 --- a/buch/papers/dreieck/teil1.tex +++ b/buch/papers/dreieck/teil1.tex @@ -3,9 +3,92 @@ % % (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil % -\section{Ordnungsstatistik und Beta-Funktion -\label{dreieck:section:ordnungsstatistik}} -\rhead{} +\section{Hermite-Polynome +\label{dreieck:section:hermite-polynome}} +\rhead{Hermite-Polyome} +In Abschnitt~\ref{dreieck:section:problemstellung} hat sich schon angedeutet, +dass die Polynome, die man durch Ableiten von $e^{-t^2}$ erhalten +kann, bezüglich des gestellten Problems besondere Eigenschaften +haben. +Zunächst halten wir fest, dass die Ableitung einer Funktion der Form +$P(t)e^{-t^2}$ mit einem Polynom $P(t)$ +\begin{equation} +\frac{d}{dt} P(t)e^{-t^2} += +P'(t)e^{-t^2} -2tP(t)e^{-t^2} += +(P'(t)-2tP(t)) e^{-t^2} +\label{dreieck:eqn:ableitung} +\end{equation} +ist. +Insbesondere hat die Ableitung wieder die Form $Q(t)e^{-t^2}$ +mit einem Polynome $Q(t)$, welches man auch als +\[ +Q(t) += +e^{t^2}\frac{d}{dt}P(t)e^{-t^2} +\] +erhalten kann. +Die Polynome, die man aus der Funktion $H_0(t)=e^{-t^2}$ durch +Ableiten erhalten kann, wurden bereits in +Abschnitt~\ref{buch:orthogonalitaet:section:rodrigues} +bis auf ein Vorzeichen hergeleitet, sie heissen die Hermite-Polynome +\index{Hermite-Polynome}% +und es gilt +\[ +H_n(t) += +(-1)^n +e^{t^2} \frac{d^n}{dt^n} e^{-t^2}. +\] +Das Vorzeichen dient dazu sicherzustellen, dass der Leitkoeffizient +immer $1$ ist. +Das Polynom $H_n(t)$ hat den Grad $n$. + +In Abschnitt wurde auch gezeigt, dass die Polynome $H_n(t)$ +bezüglich des Skalarproduktes +\[ +\langle f,g\rangle_{w} += +\int_{-\infty}^\infty f(t)g(t)e^{-t^2}\,dt, +\qquad +w(t)=e^{-t^2}, +\] +orthogonal sind. +Ausserdem folgt aus \eqref{dreieck:eqn:ableitung} +die Rekursionsbeziehung +\begin{equation} +H_{n}(t) += +2tH_{n-1}(t) +- +H_{n-1}'(t) +\label{dreieck:eqn:rekursion} +\end{equation} +für $n>0$. + +Im Hinblick auf die Problemstellung ist jetzt die Frage interessant, +ob die Integranden $H_n(t)e^{-t^2}$ eine Stammfunktion in geschlossener +Form haben. +Mit Hilfe der Rekursionsbeziehung~\eqref{dreieck:eqn:rekursion} +kann man für $n>0$ unmittelbar verifizieren, dass +\begin{align*} +\int H_n(t)e^{-t^2}\,dt +&= +\int \bigl( 2tH_{n-1}(t) - H'_{n-1}(t)\bigr)e^{-t^2}\,dt +\\ +&= +-\int \bigl( \exp'(-t^2) H_{n-1}(t) + H'_{n-1}(t)\bigr)e^{-t^2}\,dt +\\ +&= +-\int \bigl( e^{-t^2}H_{n-1}(t)\bigr)' \,dt += +-e^{-t^2}H_{n-1}(t) +\end{align*} +ist. +Für $n>0$ hat also $H_n(t)e^{-t^2}$ eine elementare Stammfunktion. +Die Hermite-Polynome sind also Lösungen für das +Problem~\ref{dreieck:problem}. diff --git a/buch/papers/dreieck/teil2.tex b/buch/papers/dreieck/teil2.tex index 83ea3cb..8e89f6a 100644 --- a/buch/papers/dreieck/teil2.tex +++ b/buch/papers/dreieck/teil2.tex @@ -3,7 +3,113 @@ % % (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil % -\section{Wahrscheinlichkeiten im Dreieckstest -\label{dreieck:section:wahrscheinlichkeiten}} -\rhead{Wahrscheinlichkeiten} +\section{Beliebige Polynome +\label{dreieck:section:beliebig}} +\rhead{Beliebige Polynome} +Im Abschnitt~\ref{dreieck:section:hermite-polynome} wurden die +Hermite-Polynome $H_n(t)$ mit $n>0$ als Lösungen des gestellten +Problems erkannt. +Eine Linearkombination von solchen Polynomen hat natürlich +ebenfalls eine elementare Stammfunktion. +Das Problem kann daher neu formuliert werden: + +\begin{problem} +\label{dreieck:problem2} +Welche Polynome $P(t)$ lassen sich aus den Hermite-Polynomen +$H_n(t)$ mit $n>0$ linear kombinieren? +\end{problem} + +Sei also +\[ +P(t) = p_0 + p_1t + \ldots + p_{n-1}t^{n-1} + p_nt^n +\] +ein beliebiges Polynom vom Grad $n$. +Eine elementare Stammfunktion von $P(t)e^{-t^2}$ existiert sicher, +wenn sich $P(t)$ aus den Funktionen $H_n(t)$ mit $n>0$ linear +kombinieren lässt. +Gesucht ist also zunächst eine Darstellung von $P(t)$ als Linearkombination +von Hermite-Polynomen. + +\begin{lemma} +Jedes Polynome $P(t)$ vom Grad $n$ lässt sich auf eindeutige Art und +Weise als Linearkombination +\begin{equation} +P(t) = a_0H_0(t) + a_1H_1(t) + \ldots + a_nH_n(t) += +\sum_{k=0}^n a_nH_n(t) +\label{dreieck:lemma} +\end{equation} +von Hermite-Polynomen schreiben. +\end{lemma} + +\begin{proof}[Beweis] +Zunächst halten wir fest, dass aus der +Rekursionsformel~\eqref{dreieck:eqn:rekursion} +folgt, dass der Leitkoeffizient bei jedem Rekursionsschnitt +mit $2$ multipliziert wird. +Der Leitkoeffizient von $H_n(t)$ ist also $2^n$. + +Wir führen den Beweis mit vollständiger Induktion. +Für $n=0$ ist $P(t)=p_0 = p_0 H_0(t)$ als Linearkombination von +Hermite-Polynomen darstellbar, dies ist die Induktionsverankerung. + +Wir nehmen jetzt im Sinne der Induktionsannahme an, +dass sich ein Polynom vom Grad $n-1$ als +Linearkombination der Polynome $H_0(t),\dots,H_{n-1}(t)$ schreiben +lässt und untersuchen ein Polynom $P(t)$ vom Grad $n$. +Da der Leitkoeffizient des Polynoms $H_n(t)$ ist $2^n$, ist zerlegen +wir +\[ +P(t) += +\underbrace{\biggl(P(t) - \frac{p_n}{2^n} H_n(t)\biggr)}_{\displaystyle = Q(t)} ++ +\frac{p_n}{2^n} H_n(t). +\] +Das Polynom $Q(t)$ hat Grad $n-1$, besitzt also nach Induktionsannahme +eine Darstellung +\[ +Q(t) = a_0H_0(t)+a_1H_1(t)+\ldots+a_{n-1}H_{n-1}(t) +\] +als Linearkombination der Polynome $H_0(t),\dots,H_{n-1}(t)$. +Somit ist +\[ +P(t) += a_0H_0(t)+a_1H_1(t)+\ldots+a_{n-1}H_{n-1}(t) + +\frac{p_n}{2^n} H_n(t) +\] +eine Darstellung von $P(t)$ als Linearkombination der Polynome +$H_0(t),\dots,H_n(t)$. +Damit ist der Induktionsschritt vollzogen und das Lemma für alle +$n$ bewiesen. +\end{proof} + +\begin{satz} +\label{dreieck:satz1} +Die Funktion $P(t)e^{-t^2}$ hat genau dann eine elementare Stammfunktion, +wenn in der Darstellung~\eqref{dreieck:lemma} +von $P(t)$ als Linearkombination von Hermite-Polynomen $a_0=0$ gilt. +\end{satz} + +\begin{proof}[Beweis] +Es ist +\begin{align*} +\int P(t)e^{-t^2}\,dt +&= +a_0\int e^{-t^2}\,dt ++ +\int +\sum_{k=1} a_kH_k(t)\,dt +\\ +&= +a_0 +\frac{\sqrt{\pi}}2 +\operatorname{erf}(t) ++ +\sum_{k=1} a_k\int H_k(t)\,dt. +\end{align*} +Da die Integrale in der Summe alle elementar darstellbar sind, +ist das Integral genau dann elementar, wenn $a_0=0$ ist. +\end{proof} + diff --git a/buch/papers/dreieck/teil3.tex b/buch/papers/dreieck/teil3.tex index e2dfd6b..c0c046a 100644 --- a/buch/papers/dreieck/teil3.tex +++ b/buch/papers/dreieck/teil3.tex @@ -3,8 +3,75 @@ % % (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil % -\section{Erweiterungen -\label{dreieck:section:erweiterungen}} -\rhead{Erweiterungen} +\section{Integralbedingung +\label{dreieck:section:integralbedingung}} +\rhead{Lösung} +Die Tatsache, dass die Hermite-Polynome orthogonal sind, erlaubt, das +Kriterium von Satz~\ref{dreieck:satz1} in einer besonders attraktiven +Integralform zu formulieren. + +Aus den Polynomen $H_n(t)$ lassen sich durch Normierung die +\index{orthogonale Polynome}% +\index{Polynome, orthogonale}% +orthonormierten Polynome +\[ +\tilde{H}_n(t) += +\frac{1}{\| H_n\|_w} H_n(t) +\qquad\text{mit}\quad +\|H_n\|_w^2 += +\int_{-\infty}^\infty H_n(t)e^{-t^2}\,dt +\] +bilden. +Da diese Polynome eine orthonormierte Basis des Vektorraums der Polynome +bilden, kann die gesuchte Zerlegung eines Polynoms $P(t)$ auch mit +Hilfe des Skalarproduktes gefunden werden: +\begin{align*} +P(t) +&= +\sum_{k=1}^n +\langle \tilde{H}_k, P\rangle_w +\tilde{H}_k(t) += +\sum_{k=1}^n +\biggl\langle \frac{H_k}{\|H_k\|_w}, P\biggr\rangle_w +\frac{H_k(t)}{\|H_k\|_w} += +\sum_{k=1}^n +\underbrace{ +\frac{ \langle H_k, P\rangle_w }{\|H_k\|_w^2} +}_{\displaystyle =a_k} +H_k(t). +\end{align*} +Die Darstellung von $P(t)$ als Linearkombination von Hermite-Polynomen +hat somit die Koeffizienten +\[ +a_k = \frac{\langle H_k,P\rangle_w}{\|H_k\|_w^2}. +\] +Aus dem Kriterium $a_0=0$ dafür, dass eine elementare Stammfunktion +von $P(t)e^{-t^2}$ existiert, wird daher die Bedingung, dass +$\langle H_0,P\rangle_w=0$ ist. +Da $H_0(t)=1$ ist, folgt als Bedingung +\[ +a_0 += +\langle H_0,P\rangle_w += +\int_{-\infty}^\infty P(t) e^{-t^2}\,dt += +0. +\] + +\begin{satz} +Ein Integrand der Form $P(t)e^{-t^2}$ mit einem Polynom $P(t)$ +hat genau dann eine elementare Stammfunktion, wenn +\[ +\int_{-\infty}^\infty P(t)e^{-t^2}\,dt = 0 +\] +ist. +\end{satz} + + diff --git a/vorlesungen/slides/hermite/hermiteentwicklung.tex b/vorlesungen/slides/hermite/hermiteentwicklung.tex index e1ced30..5f6e1c9 100644 --- a/vorlesungen/slides/hermite/hermiteentwicklung.tex +++ b/vorlesungen/slides/hermite/hermiteentwicklung.tex @@ -17,6 +17,7 @@ P(x) = p_0 + p_1x + p_2x^2 + \dots + p_nx^n \] +\uncover<2->{% als Linearkombination von Hermite-Polynome schreiben: \begin{align*} P(x) @@ -38,10 +39,11 @@ a_0\cdot 1 &\quad\;\;\vdots \\ &\quad + a_n(2^nx^n + \dots) -\end{align*} +\end{align*}} \end{block} \end{column} \begin{column}{0.48\textwidth} +\uncover<3->{% \begin{block}{Koeffizientenvergleich} führt auf ein Gleichungssystem \begin{center} @@ -58,11 +60,12 @@ a_0&a_1&a_2&a_3&a_4&\dots&\\ \hline \end{tabular} \end{center} -Dreiecksmatrix, Diagonalelement -$\ne 0$ -$\Rightarrow$ -$\exists$ eindeutige Lösung -\end{block} +\uncover<4->{% +Dreiecksmatrix}\uncover<5->{, Diagonalelement +$\ne 0$} +\uncover<6->{$\Rightarrow$ +$\exists$ eindeutige Lösung} +\end{block}} \end{column} \end{columns} \end{frame} diff --git a/vorlesungen/slides/hermite/loesung.tex b/vorlesungen/slides/hermite/loesung.tex index 7d4741f..68ee32e 100644 --- a/vorlesungen/slides/hermite/loesung.tex +++ b/vorlesungen/slides/hermite/loesung.tex @@ -20,36 +20,45 @@ P(t)e^{-\frac{t^2}2} \] in geschlossener Form angeben? \end{block} +\uncover<2->{% \begin{block}{``Hermite-Antwort''} \[ \int H_n(x)e^{-x^2}\,dx \] kann genau für $n>0$ in geschlossener Form angegeben werden. -\end{block} +\end{block}} \end{column} \begin{column}{0.48\textwidth} +\uncover<3->{% \begin{block}{Allgemein} \begin{align*} \int P(x)e^{-x^2}\,dx -&= -\int \sum_{k=0}^n a_kH_k(x)e^{-x^2}\,dx +&\uncover<4->{= +\int \sum_{k=0}^n a_kH_k(x)e^{-x^2}\,dx} \\ +\uncover<5->{ &= \sum_{k=0}^n a_k \int H_k(x)e^{-x^2}\,dx +} \\ +\uncover<6->{ &= a_0\operatorname{erf}(x) + C +} \\ +\uncover<6->{ &\hspace*{2mm} + \sum_{k=1}^n a_k\int H_k(x)e^{-x^2}\,dx +} \end{align*} -\end{block} +\end{block}} +\uncover<7->{% \begin{theorem} Das Integral von $P(x)e^{-x^2}$ ist genau dann elementar darstellbar, wenn $a_0=0$ -\end{theorem} +\end{theorem}} \end{column} \end{columns} \end{frame} diff --git a/vorlesungen/slides/hermite/normalhermite.tex b/vorlesungen/slides/hermite/normalhermite.tex index 16a314c..98721dc 100644 --- a/vorlesungen/slides/hermite/normalhermite.tex +++ b/vorlesungen/slides/hermite/normalhermite.tex @@ -19,6 +19,7 @@ H_n(x) \] \end{block} \vspace{-10pt} +\uncover<2->{% \begin{block}{Orthogonalität} $H_n(x)$ sind orthogonale Polynome bezüglich $w(x)=e^{-x^2}$, d.~h. \begin{align*} @@ -37,8 +38,9 @@ $H_n(x)$ sind orthogonale Polynome bezüglich $w(x)=e^{-x^2}$, d.~h. = \delta_{mn} \end{align*} -\end{block} +\end{block}} \vspace{-10pt} +\uncover<3->{% \begin{block}{Rekursion: Auf-/Absteigeoperatoren} Rekursionsformel: \[ @@ -46,33 +48,46 @@ H_n(x) = 2x\cdot H_{n-1}(x) - H_{n-1}'(x) \] -\end{block} +\end{block}} \end{column} \begin{column}{0.48\textwidth} +\uncover<4->{% \begin{block}{Stammfunktion} \begin{align*} -\int H_n(x) e^{-x^2}\,dx -&= -\int \bigl({\color{red}2x}H_{n-1}(x) +\uncover<4->{ +\int H_n(x) e^{-x^2}\,dx} +&\uncover<5->{= +\int \bigl({\color{red}2x}H_{n-1}(x)} \\ +\uncover<5->{ &\qquad -H_{n-1}'(x)\bigr) e^{-x^2}\,dx +} \\ +\uncover<6->{ {\color{gray}((e^{-x^2})'=-2x)} &= {\color{red}-}\int {\color{red}(e^{-x^2})'} H_{n-1}(x)\,dx +} \\ +\uncover<6->{ &\qquad - \int H_{n-1}'(x) e^{-x^2}\,dx +} \\ +\uncover<7->{ \text{\color{gray}(Produktregel)} &= \int (e^{-x^2}H_{n-1}(x))'\,dx +} \\ +\uncover<8->{ \text{\color{gray}(Ableitung)} &= e^{-x^2}H_{n-1}(x) +} \end{align*} +\uncover<9->{% ausser für $n=0$: \[ \int @@ -80,8 +95,8 @@ H_0(x)e^{-x^2}\,dx = \int e^{-x^2}\,dx -\] -\end{block} +\]} +\end{block}} \end{column} \end{columns} \end{frame} diff --git a/vorlesungen/slides/hermite/normalintegrale.tex b/vorlesungen/slides/hermite/normalintegrale.tex index 88abbe8..32333cd 100644 --- a/vorlesungen/slides/hermite/normalintegrale.tex +++ b/vorlesungen/slides/hermite/normalintegrale.tex @@ -20,12 +20,14 @@ P(t)e^{-t^2} \] in geschlossener Form angeben? \end{block} +\uncover<4->{% \begin{block}{Allgemeine Antwort} Satz von Liouville und Risch- Algorithmus können entscheiden, ob es eine elementare Stammfunktion gibt -\end{block} +\end{block}} \end{column} \begin{column}{0.48\textwidth} +\uncover<2->{% \begin{block}{Negativbeispiel} $P(t) = 1$, das Normalverteilungsintegral \[ @@ -34,7 +36,8 @@ F(x) \frac{1}{\sqrt{2\pi}} \int_{-\infty}^x e^{-t^2}\,dt \] ist nicht elementar darstellbar. -\end{block} +\end{block}} +\uncover<3->{% \begin{block}{Positivbeispiel} $P(t)=t$. Wegen \begin{align*} @@ -47,7 +50,7 @@ $P(t)=t$. Wegen -e^{-x^2}+C \end{align*} elementar darstellbar. -\end{block} +\end{block}} \end{column} \end{columns} \end{frame} diff --git a/vorlesungen/slides/hermite/skalarprodukt.tex b/vorlesungen/slides/hermite/skalarprodukt.tex index 32b933f..a51e9f6 100644 --- a/vorlesungen/slides/hermite/skalarprodukt.tex +++ b/vorlesungen/slides/hermite/skalarprodukt.tex @@ -18,6 +18,7 @@ Orthogonale $H_k$ normalisieren: \] mit Gewichtsfunktion $w(x)=e^{-x^2}$ \end{block} +\uncover<2->{% \begin{block}{``Hermite''-Analyse} \begin{align*} P(x) @@ -26,46 +27,55 @@ P(x) = \sum_{k=1}^\infty \tilde{a}_k \tilde{H}_k(x) \\ +\uncover<3->{ \tilde{a}_k &= \| H_k\|_w\, a_k +} \\ +\uncover<4->{ a_k &= \frac{1}{\|H_k\|} \langle \tilde{H}_k, P\rangle_w -= +}\uncover<5->{= \frac{1}{\|H_k\|^2} \langle H_k, P\rangle_w +} \end{align*} -\end{block} +\end{block}} \end{column} \begin{column}{0.48\textwidth} +\uncover<6->{% \begin{block}{Integrationsproblem} Bedingung: \begin{align*} a_0=0 +\uncover<7->{% \qquad\Leftrightarrow\qquad \langle H_0,P\rangle_w &= -0 +0} \\ +\uncover<8->{% \int_{-\infty}^\infty P(t) w(t) \,dt +}\uncover<9->{% = \int_{-\infty}^\infty P(t) e^{-t^2} \,dt &= -0 +0} \end{align*} -\end{block} +\end{block}} +\uncover<10->{% \begin{theorem} Das Integral von $P(t)e^{-t^2}$ ist in geschlossener Form darstellbar genau dann, wenn \[ \int_{-\infty}^\infty P(t)e^{-t^2}\,dt = 0 \] -\end{theorem} +\end{theorem}} \end{column} \end{columns} \end{frame} |