From f5047d4d780e996a8b8f7738c1ac7c884a07f135 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Sun, 13 Mar 2022 23:26:58 +0100 Subject: new stuff about beta, test2 --- buch/aufgaben2.tex | 4 + buch/chapters/040-rekursion/Makefile.inc | 1 + buch/chapters/040-rekursion/beta.tex | 104 ++------- buch/chapters/040-rekursion/images/Makefile | 16 +- buch/chapters/040-rekursion/images/beta.pdf | Bin 0 -> 109772 bytes buch/chapters/040-rekursion/images/beta.tex | 236 +++++++++++++++++++++ buch/chapters/040-rekursion/images/betadist.m | 58 +++++ buch/chapters/040-rekursion/images/order.m | 119 +++++++++++ buch/chapters/040-rekursion/images/order.pdf | Bin 0 -> 32692 bytes buch/chapters/040-rekursion/images/order.tex | 125 +++++++++++ buch/chapters/070-orthogonalitaet/Makefile.inc | 1 + buch/chapters/070-orthogonalitaet/chapter.tex | 2 +- .../070-orthogonalitaet/uebungsaufgaben/701.tex | 137 ++++++++++++ buch/chapters/090-pde/Makefile.inc | 1 + buch/chapters/090-pde/chapter.tex | 12 +- buch/chapters/090-pde/kreis.tex | 2 +- buch/chapters/090-pde/uebungsaufgaben/901.tex | 82 +++++++ buch/papers/dreieck/images/beta.pdf | Bin 100791 -> 109717 bytes buch/papers/dreieck/images/beta.tex | 208 ++++++++++-------- buch/papers/dreieck/images/betadist.m | 24 ++- 20 files changed, 938 insertions(+), 194 deletions(-) create mode 100644 buch/chapters/040-rekursion/images/beta.pdf create mode 100644 buch/chapters/040-rekursion/images/beta.tex create mode 100644 buch/chapters/040-rekursion/images/betadist.m create mode 100644 buch/chapters/040-rekursion/images/order.m create mode 100644 buch/chapters/040-rekursion/images/order.pdf create mode 100644 buch/chapters/040-rekursion/images/order.tex create mode 100644 buch/chapters/070-orthogonalitaet/uebungsaufgaben/701.tex create mode 100644 buch/chapters/090-pde/uebungsaufgaben/901.tex (limited to 'buch') diff --git a/buch/aufgaben2.tex b/buch/aufgaben2.tex index bed14fb..f98562e 100644 --- a/buch/aufgaben2.tex +++ b/buch/aufgaben2.tex @@ -8,4 +8,8 @@ %\input chapters/40-eigenwerte/uebungsaufgaben/4004.tex %\item %\input chapters/40-eigenwerte/uebungsaufgaben/4005.tex +\item +\input{chapters/090-pde/uebungsaufgaben/901.tex} +\item +\input{chapters/070-orthogonalitaet/uebungsaufgaben/701.tex} diff --git a/buch/chapters/040-rekursion/Makefile.inc b/buch/chapters/040-rekursion/Makefile.inc index ed8fd51..a222b1c 100644 --- a/buch/chapters/040-rekursion/Makefile.inc +++ b/buch/chapters/040-rekursion/Makefile.inc @@ -9,6 +9,7 @@ CHAPTERFILES = $(CHAPTERFILES) \ chapters/040-rekursion/bohrmollerup.tex \ chapters/040-rekursion/integral.tex \ chapters/040-rekursion/beta.tex \ + chapters/040-rekursion/betaverteilung.tex \ chapters/040-rekursion/linear.tex \ chapters/040-rekursion/hypergeometrisch.tex \ chapters/040-rekursion/uebungsaufgaben/401.tex \ diff --git a/buch/chapters/040-rekursion/beta.tex b/buch/chapters/040-rekursion/beta.tex index ea847bc..ff59bad 100644 --- a/buch/chapters/040-rekursion/beta.tex +++ b/buch/chapters/040-rekursion/beta.tex @@ -3,11 +3,17 @@ % % (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule % -\subsection{Die Beta-Funktion -\label{buch:rekursion:gamma:subsection:beta}} +\section{Die Beta-Funktion +\label{buch:rekursion:gamma:section:beta}} Die Eulersche Integralformel für die Gamma-Funktion in -Definition~\ref{buch:rekursion:def:gamma} wurde bisher nicht -gerechtfertigt. +Definition~\ref{buch:rekursion:def:gamma} wurde in +Abschnitt~\ref{buch:subsection:integral-eindeutig} +mit dem Satz von Mollerup gerechtfertigt. +Man kann Sie aber auch als Grenzfall der Beta-Funktion verstehen, +die in diesem Abschnitt dargestellt wird. + + +\subsection{Beta-Integral} In diesem Abschnitt wird das Beta-Integral eingeführt, eine Funktion von zwei Variablen, welches eine Integral-Definition mit einer reichaltigen Menge von Rekursionsbeziehungen hat, die sich direkt auf @@ -233,6 +239,16 @@ B(x,y) = \frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)} berechnet werden. \end{satz} +% +% Info über die Beta-Verteilung +% +\input{chapters/040-rekursion/betaverteilung.tex} + +\subsection{Weitere Eigenschaften der Gamma-Funktion} +Die nahe Verwandtschaft der Gamma- mit der Beta-Funktion ermöglicht +nun, weitere Eigenschaften der Gamma-Funktion mit Hilfe der Beta-Funktion +herzuleiten. + \subsubsection{Nochmals der Wert von $\Gamma(\frac12)$?} Der Wert von $\Gamma(\frac12)=\sqrt{\pi}$ wurde bereits in \eqref{buch:rekursion:gamma:wert12} @@ -484,83 +500,3 @@ Setzt man $x=\frac12$ in die Verdoppelungsformel ein, erhält man in Übereinstimmung mit dem aus \eqref{buch:rekursion:gamma:gamma12} bereits bekannten Wert. -\subsubsection{Beta-Funktion und Binomialkoeffizienten} -Die Binomialkoeffizienten können mit Hilfe der Fakultät als -\begin{align*} -\binom{n}{k} -&= -\frac{n!}{(n-k)!\,k!} -\intertext{geschrieben werden. -Drückt man die Fakultäten durch die Gamma-Funktion aus, erhält man} -&= -\frac{\Gamma(n+1)}{\Gamma(n-k+1)\Gamma(k+1)}. -\intertext{Schreibt man $x=k-1$ und $y=n-k+1$, wird daraus -wegen $x+y=k+1+n-k+1=n+2=(n+1)+1$} -&= -\frac{\Gamma(x+y-1)}{\Gamma(x)\Gamma(y)}. -\intertext{Die Rekursionsformel für die Gamma-Funktion erlaubt, -den Zähler umzuwandeln in $\Gamma(x+y-1)=\Gamma(x+y)/(x+y-1)$, so dass -der Binomialkoeffizient schliesslich} -&= -\frac{\Gamma(x+y)}{(x+y-1)\Gamma(x)\Gamma(y)} -= -\frac{1}{(n-1)B(n-k+1,k+1)} -\label{buch:rekursion:gamma:binombeta} -\end{align*} -geschrieben werden kann. -Die Rekursionsbeziehung -\[ -\binom{n+1}{k} = \binom{n}{k-1} + \binom{n}{k} -\] -der Binomialkoeffizienten erzeugt das vertraute Pascal-Dreieck, -die Formel \eqref{buch:rekursion:gamma:binombeta} für die -Binomialkoeffizienten macht daraus -\[ -\frac{n-1}{B(n-k,k-1)} -= -\frac{n-2}{B(n-k,k-2)} -+ -\frac{n-2}{B(n-k-1,k-1)}, -\] -die für ganzzahlige Argumente gilt. -Wir wollen nachrechnen, dass dies für beliebige Argumente gilt. -\begin{align*} -\frac{(n-1)\Gamma(n-1)}{\Gamma(n-k)\Gamma(k-1)} -&= -\frac{(n-2)\Gamma(n-2)}{\Gamma(n-k)\Gamma(k-2)} -+ -\frac{(n-2)\Gamma(n-2)}{\Gamma(n-k-1)\Gamma(k-1)} -\\ -\frac{\Gamma(n)}{\Gamma(n-k)\Gamma(k-1)} -&= -\frac{\Gamma(n-1)}{\Gamma(n-k)\Gamma(k-2)} -+ -\frac{\Gamma(n-1)}{\Gamma(n-k-1)\Gamma(k-1)} -\intertext{Durch Zusammenfassen der Faktoren im Zähler mit Hilfe -der Rekursionsformel für die Gamma-Funktion und Multiplizieren -mit dem gemeinsamen Nenner -$\Gamma(n-k)\Gamma(k-1)=(n-k-1)\Gamma(n-k-1)(k-2)\Gamma(k-2)$ wird daraus} -\Gamma(n) -&= -(k-2) -\Gamma(n-1) -+ -(n-k-1) -\Gamma(n-1) -\intertext{Indem wir die Rekursionsformel für die Gamma-Funktion auf -die rechte Seite anwenden können wir erreichen, dass in allen Termen -ein Faktor -$\Gamma(n-1)$ auftritt:} -(n-1)\Gamma(n-1) -&= -(k-2)\Gamma(n-1) -+ -(n+k-1)\Gamma(n-1) -\\ -n-1 -&= -k-2 -+ -n-k-1 -\end{align*} - diff --git a/buch/chapters/040-rekursion/images/Makefile b/buch/chapters/040-rekursion/images/Makefile index 9608a94..86dfa1e 100644 --- a/buch/chapters/040-rekursion/images/Makefile +++ b/buch/chapters/040-rekursion/images/Makefile @@ -3,7 +3,7 @@ # # (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule # -all: gammaplot.pdf fibonacci.pdf +all: gammaplot.pdf fibonacci.pdf order.pdf beta.pdf gammaplot.pdf: gammaplot.tex gammapaths.tex pdflatex gammaplot.tex @@ -16,3 +16,17 @@ fibonaccigrid.tex: fibonacci.m fibonacci.pdf: fibonacci.tex fibonaccigrid.tex pdflatex fibonacci.tex + +order.pdf: order.tex orderpath.tex + pdflatex order.tex + +orderpath.tex: order.m + octave order.m + +beta.pdf: beta.tex betapaths.tex + pdflatex beta.tex + +betapaths.tex: betadist.m + octave betadist.m + + diff --git a/buch/chapters/040-rekursion/images/beta.pdf b/buch/chapters/040-rekursion/images/beta.pdf new file mode 100644 index 0000000..0e6567b Binary files /dev/null and b/buch/chapters/040-rekursion/images/beta.pdf differ diff --git a/buch/chapters/040-rekursion/images/beta.tex b/buch/chapters/040-rekursion/images/beta.tex new file mode 100644 index 0000000..1e1a1b3 --- /dev/null +++ b/buch/chapters/040-rekursion/images/beta.tex @@ -0,0 +1,236 @@ +% +% beta.tex -- display some symmetric beta distributions +% +% (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +% +\documentclass[tikz]{standalone} +\usepackage{amsmath} +\usepackage{times} +\usepackage{txfonts} +\usepackage{pgfplots} +\usepackage{csvsimple} +\usetikzlibrary{arrows,intersections,math,calc} +\input{betapaths.tex} +\begin{document} +\def\skala{12} +\definecolor{colorone}{rgb}{1.0,0.6,0.0} +\definecolor{colortwo}{rgb}{1.0,0.0,0.0} +\definecolor{colorthree}{rgb}{0.6,0.0,0.6} +\definecolor{colorfour}{rgb}{0.6,0.0,1.0} +\definecolor{colorfive}{rgb}{0.0,0.0,1.0} +\definecolor{colorsix}{rgb}{0.4,0.6,1.0} +\definecolor{colorseven}{rgb}{0.0,0.0,0.0} +\definecolor{coloreight}{rgb}{0.0,0.8,0.8} +\definecolor{colornine}{rgb}{0.0,0.8,0.2} +\definecolor{colorten}{rgb}{0.2,0.4,0.0} +\definecolor{coloreleven}{rgb}{0.6,1.0,0.0} +\definecolor{colortwelve}{rgb}{1.0,0.8,0.4} + +\def\achsen{ + \foreach \x in {0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9}{ + \draw ({\x*\dx},{-0.1/\skala}) -- ({\x*\dx},{0.1/\skala}); + \node at ({\x*\dx},{-0.1/\skala}) [below] {$\x$}; + } + \foreach \y in {1,2,3,4}{ + \draw ({-0.1/\skala},{\y*\dy}) -- ({0.1/\skala},{\y*\dy}); + \node at ({-0.1/\skala},{\y*\dy}) [left] {$\y$}; + } + \def\x{1} + \draw ({\x*\dx},{-0.1/\skala}) -- ({\x*\dx},{0.1/\skala}); + \node at ({\x*\dx},{-0.1/\skala}) [below] {$\x$}; + \def\x{0} + \node at ({\x*\dx},{-0.1/\skala}) [below] {$\x$}; + + \draw[->] ({-0.1/\skala},0) -- ({1*\dx+0.4/\skala},0) + coordinate[label={$x$}]; + \draw[->] (0,{-0.1/\skala}) -- (0,{\betamax*\dy+0.4/\skala},0) + coordinate[label={right:$\beta(a,b,x)$}]; +} + +\def\farbcoord#1#2{ + ({\dx*(0.63+((#1)/5)*0.27)},{\dx*(0.18+((#2)/5)*0.27)}) +} +\def\farbviereck{ + \foreach \x in {1,2,3,4}{ + \draw[color=gray!30] \farbcoord{\x}{0} -- \farbcoord{\x}{4}; + \draw[color=gray!30] \farbcoord{0}{\x} -- \farbcoord{4}{\x}; + } + \draw[->] \farbcoord{0}{0} -- \farbcoord{4.4}{0} + coordinate[label={$a$}]; + \draw[->] \farbcoord{0}{0} -- \farbcoord{0}{4.4} + coordinate[label={left: $b$}]; + \foreach \x in {1,2,3,4}{ + \node[color=gray] at \farbcoord{4}{\x} [right] {\tiny $b=\x$}; + %\fill[color=white,opacity=0.7] + % \farbcoord{(\x-0.1)}{3.3} + % rectangle + % \farbcoord{(\x+0.1)}{4}; + \node[color=gray] at \farbcoord{\x}{4} [right,rotate=90] + {\tiny $a=\x$}; + } +} +\def\farbpunkt#1#2#3{ + \fill[color=#3] \farbcoord{#1}{#2} circle[radius={0.1/\skala}]; +} + +\begin{tikzpicture}[>=latex,thick,scale=\skala] + +\def\dx{1.15} +\def\dy{0.1} +\def\opa{0.1} + +\def\betamax{4.9} + +\begin{scope} +\clip (0,0) rectangle ({1*\dx},{\betamax*\dy}); +\fill[color=colorone,opacity=\opa] (0,0) -- \betaaa -- (\dx,0) -- cycle; +\fill[color=colortwo,opacity=\opa] (0,0) -- \betabb -- (\dx,0) -- cycle; +\fill[color=colorthree,opacity=\opa] (0,0) -- \betacc -- (\dx,0) -- cycle; +\fill[color=colorfour,opacity=\opa] (0,0) -- \betadd -- (\dx,0) -- cycle; +\fill[color=colorfive,opacity=\opa] (0,0) -- \betaee -- (\dx,0) -- cycle; +\fill[color=colorsix,opacity=\opa] (0,0) -- \betaff -- (\dx,0) -- cycle; +\fill[color=colorseven,opacity=\opa] (0,0) -- \betagg -- (\dx,0) -- cycle; +\fill[color=coloreight,opacity=\opa] (0,0) -- \betahh -- (\dx,0) -- cycle; +\fill[color=colornine,opacity=\opa] (0,0) -- \betaii -- (\dx,0) -- cycle; +\fill[color=colorten,opacity=\opa] (0,0) -- \betajj -- (\dx,0) -- cycle; +\fill[color=coloreleven,opacity=\opa] (0,0) -- \betakk -- (\dx,0) -- cycle; +\fill[color=colortwelve,opacity=\opa] (0,0) -- \betall -- (\dx,0) -- cycle; + +\draw[color=colorone] \betaaa; +\draw[color=colortwo] \betabb; +\draw[color=colorthree] \betacc; +\draw[color=colorfour] \betadd; +\draw[color=colorfive] \betaee; +\draw[color=colorsix] \betaff; +\draw[color=colorseven] \betagg; +\draw[color=coloreight] \betahh; +\draw[color=colornine] \betaii; +\draw[color=colorten] \betajj; +\draw[color=coloreleven] \betakk; +\draw[color=colortwelve] \betall; + +\end{scope} + +\achsen + +\farbviereck + +\farbpunkt{\alphatwelve}{\betatwelve}{colortwelve} +\farbpunkt{\alphaeleven}{\betaeleven}{coloreleven} +\farbpunkt{\alphaten}{\betaten}{colorten} +\farbpunkt{\alphanine}{\betanine}{colornine} +\farbpunkt{\alphaeight}{\betaeight}{coloreight} +\farbpunkt{\alphaseven}{\betaseven}{colorseven} +\farbpunkt{\alphasix}{\betasix}{colorsix} +\farbpunkt{\alphafive}{\betafive}{colorfive} +\farbpunkt{\alphafour}{\betafour}{colorfour} +\farbpunkt{\alphathree}{\betathree}{colorthree} +\farbpunkt{\alphatwo}{\betatwo}{colortwo} +\farbpunkt{\alphaone}{\betaone}{colorone} + + +\def\betamax{4.9} + +\begin{scope}[yshift=-0.6cm] + +\begin{scope} +\clip (0,0) rectangle ({1*\dx},{\betamax*\dy}); +\fill[color=colorone,opacity=\opa] (0,0) -- \betaea -- (\dx,0) -- cycle; +\fill[color=colortwo,opacity=\opa] (0,0) -- \betaeb -- (\dx,0) -- cycle; +\fill[color=colorthree,opacity=\opa] (0,0) -- \betaec -- (\dx,0) -- cycle; +\fill[color=colorfour,opacity=\opa] (0,0) -- \betaed -- (\dx,0) -- cycle; +\fill[color=colorfive,opacity=\opa] (0,0) -- \betaee -- (\dx,0) -- cycle; +\fill[color=colorsix,opacity=\opa] (0,0) -- \betaef -- (\dx,0) -- cycle; +\fill[color=colorseven,opacity=\opa] (0,0) -- \betaeg -- (\dx,0) -- cycle; +\fill[color=coloreight,opacity=\opa] (0,0) -- \betaeh -- (\dx,0) -- cycle; +\fill[color=colornine,opacity=\opa] (0,0) -- \betaei -- (\dx,0) -- cycle; +\fill[color=colorten,opacity=\opa] (0,0) -- \betaej -- (\dx,0) -- cycle; +\fill[color=coloreleven,opacity=\opa] (0,0) -- \betaek -- (\dx,0) -- cycle; +\fill[color=colortwelve,opacity=\opa] (0,0) -- \betael -- (\dx,0) -- cycle; + +\draw[color=colorone] \betaea; +\draw[color=colortwo] \betaeb; +\draw[color=colorthree] \betaec; +\draw[color=colorfour] \betaed; +\draw[color=colorfive] \betaee; +\draw[color=colorsix] \betaef; +\draw[color=colorseven] \betaeg; +\draw[color=coloreight] \betaeh; +\draw[color=colornine] \betaei; +\draw[color=colorten] \betaej; +\draw[color=coloreleven] \betaek; +\draw[color=colortwelve] \betael; +\end{scope} + +\achsen + +\farbviereck + +\farbpunkt{\alphafive}{\betatwelve}{colortwelve} +\farbpunkt{\alphafive}{\betaeleven}{coloreleven} +\farbpunkt{\alphafive}{\betaten}{colorten} +\farbpunkt{\alphafive}{\betanine}{colornine} +\farbpunkt{\alphafive}{\betaeight}{coloreight} +\farbpunkt{\alphafive}{\betaseven}{colorseven} +\farbpunkt{\alphafive}{\betasix}{colorsix} +\farbpunkt{\alphafive}{\betafive}{colorfive} +\farbpunkt{\alphafive}{\betafour}{colorfour} +\farbpunkt{\alphafive}{\betathree}{colorthree} +\farbpunkt{\alphafive}{\betatwo}{colortwo} +\farbpunkt{\alphafive}{\betaone}{colorone} + +\end{scope} + +\begin{scope}[yshift=-1.2cm] + +\begin{scope} +\clip (0,0) rectangle ({1*\dx},{\betamax*\dy}); +\fill[color=colorone,opacity=\opa] (0,0) -- \betaal -- (\dx,0) -- cycle; +\fill[color=colortwo,opacity=\opa] (0,0) -- \betabl -- (\dx,0) -- cycle; +\fill[color=colorthree,opacity=\opa] (0,0) -- \betacl -- (\dx,0) -- cycle; +\fill[color=colorfour,opacity=\opa] (0,0) -- \betadl -- (\dx,0) -- cycle; +\fill[color=colorfive,opacity=\opa] (0,0) -- \betael -- (\dx,0) -- cycle; +\fill[color=colorsix,opacity=\opa] (0,0) -- \betafl -- (\dx,0) -- cycle; +\fill[color=colorseven,opacity=\opa] (0,0) -- \betagl -- (\dx,0) -- cycle; +\fill[color=coloreight,opacity=\opa] (0,0) -- \betahl -- (\dx,0) -- cycle; +\fill[color=colornine,opacity=\opa] (0,0) -- \betail -- (\dx,0) -- cycle; +\fill[color=colorten,opacity=\opa] (0,0) -- \betajl -- (\dx,0) -- cycle; +\fill[color=coloreleven,opacity=\opa] (0,0) -- \betakl -- (\dx,0) -- cycle; +\fill[color=colortwelve,opacity=\opa] (0,0) -- \betall -- (\dx,0) -- cycle; + +\draw[color=colorone] \betaal; +\draw[color=colortwo] \betabl; +\draw[color=colorthree] \betacl; +\draw[color=colorfour] \betadl; +\draw[color=colorfive] \betael; +\draw[color=colorsix] \betafl; +\draw[color=colorseven] \betagl; +\draw[color=coloreight] \betahl; +\draw[color=colornine] \betail; +\draw[color=colorten] \betajl; +\draw[color=coloreleven] \betakl; +\draw[color=colortwelve] \betall; +\end{scope} + +\achsen + +\farbviereck + +\farbpunkt{\alphatwelve}{\betatwelve}{colortwelve} +\farbpunkt{\alphaeleven}{\betatwelve}{coloreleven} +\farbpunkt{\alphaten}{\betatwelve}{colorten} +\farbpunkt{\alphanine}{\betatwelve}{colornine} +\farbpunkt{\alphaeight}{\betatwelve}{coloreight} +\farbpunkt{\alphaseven}{\betatwelve}{colorseven} +\farbpunkt{\alphasix}{\betatwelve}{colorsix} +\farbpunkt{\alphafive}{\betatwelve}{colorfive} +\farbpunkt{\alphafour}{\betatwelve}{colorfour} +\farbpunkt{\alphathree}{\betatwelve}{colorthree} +\farbpunkt{\alphatwo}{\betatwelve}{colortwo} +\farbpunkt{\alphaone}{\betatwelve}{colorone} + +\end{scope} + +\end{tikzpicture} +\end{document} + diff --git a/buch/chapters/040-rekursion/images/betadist.m b/buch/chapters/040-rekursion/images/betadist.m new file mode 100644 index 0000000..5b466a6 --- /dev/null +++ b/buch/chapters/040-rekursion/images/betadist.m @@ -0,0 +1,58 @@ +# +# betadist.m +# +# (c) 2022 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +# +global N; +N = 201; +global nmin; +global nmax; +nmin = -4; +nmax = 7; +n = nmax - nmin + 1 +A = 3; + +t = (nmin:nmax) / nmax; +alpha = 1 + A * t .* abs(t) +#alpha(1) = 0.01; + +#alpha = [ 1, 1.03, 1.05, 1.1, 1.25, 1.5, 2, 2.5, 3, 4, 5 ]; +beta = alpha; +names = [ "one"; "two"; "three"; "four"; "five"; "six"; "seven"; "eight"; + "nine"; "ten"; "eleven"; "twelve" ] + +function retval = Beta(a, b, x) + retval = x^(a-1) * (1-x)^(b-1) / beta(a, b); + if (retval > 100) + retval = 100 + end +end + +function plotbeta(fn, a, b, name) + global N; + fprintf(fn, "\\def\\beta%s{\n", strtrim(name)); + fprintf(fn, "\t({%.4f*\\dx},{%.4f*\\dy})", 0, Beta(a, b, 0)); + for x = (1:N-1)/(N-1) + X = (1-cos(pi * x))/2; + fprintf(fn, "\n\t--({%.4f*\\dx},{%.4f*\\dy})", + X, Beta(a, b, X)); + end + fprintf(fn, "\n}\n"); +end + +fn = fopen("betapaths.tex", "w"); + +for i = (1:n) + fprintf(fn, "\\def\\alpha%s{%f}\n", strtrim(names(i,:)), alpha(i)); + fprintf(fn, "\\def\\beta%s{%f}\n", strtrim(names(i,:)), beta(i)); +end + +for i = (1:n) + for j = (1:n) + printf("working on %d,%d:\n", i, j); + plotbeta(fn, alpha(i), beta(j), + char(['a' + i - 1, 'a' + j - 1])); + end +end + +fclose(fn); diff --git a/buch/chapters/040-rekursion/images/order.m b/buch/chapters/040-rekursion/images/order.m new file mode 100644 index 0000000..762f458 --- /dev/null +++ b/buch/chapters/040-rekursion/images/order.m @@ -0,0 +1,119 @@ +# +# order.m +# +# (c) 2022 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +# +global N; +N = 10; +global subdivisions; +subdivisions = 100; +global P; +P = 0.5 + +function retval = orderF(p, n, k) + retval = 0; + for i = (k:n) + retval = retval + nchoosek(n,i) * p^i * (1-p)^(n-i); + end +end + +function retval = orderd(p, n, k) + retval = 0; + for i = (k:n) + s = i * p^(i-1) * (1-p)^(n-i); + s = s - p^i * (n-i) * (1-p)^(n-i-1); + retval = retval + nchoosek(n,i) * s; + end +end + +function retval = orders(p, n, k) + retval = k * nchoosek(n, k) * p^(k-1) * (1-p)^(n-k); +end + +function orderpath(fn, k, name) + fprintf(fn, "\\def\\order%s{\n\t(0,0)", name); + global N; + global subdivisions; + for i = (0:subdivisions) + p = i/subdivisions; + fprintf(fn, "\n\t-- ({%.4f*\\dx},{%.4f*\\dy})", + p, orderF(p, N, k)); + end + fprintf(fn, "\n}\n"); +end + +function orderdpath(fn, k, name) + fprintf(fn, "\\def\\orderd%s{\n\t(0,0)", name); + global N; + global subdivisions; + for i = (1:subdivisions-1) + p = i/subdivisions; + fprintf(fn, "\n\t-- ({%.4f*\\dx},{%.4f*\\dy})", + p, orderd(p, N, k)); + end + fprintf(fn, "\n\t-- ({1*\\dx},0)"); + fprintf(fn, "\n}\n"); +end + +function orderspath(fn, k, name) + fprintf(fn, "\\def\\orders%s{\n\t(0,0)", name); + global N; + global subdivisions; + for i = (1:subdivisions-1) + p = i/subdivisions; + fprintf(fn, "\n\t-- ({%.4f*\\dx},{%.4f*\\dy})", + p, orders(p, N, k)); + end + fprintf(fn, "\n\t-- ({1*\\dx},0)"); + fprintf(fn, "\n}\n"); +end + +fn = fopen("orderpath.tex", "w"); + +orderpath(fn, 0, "zero"); +orderdpath(fn, 0, "zero"); +orderspath(fn, 0, "zero"); + +orderpath(fn, 1, "one"); +orderdpath(fn, 1, "one"); +orderspath(fn, 1, "one"); + +orderpath(fn, 2, "two"); +orderdpath(fn, 2, "two"); +orderspath(fn, 2, "two"); + +orderpath(fn, 3, "three"); +orderdpath(fn, 3, "three"); +orderspath(fn, 3, "three"); + +orderpath(fn, 4, "four"); +orderdpath(fn, 4, "four"); +orderspath(fn, 4, "four"); + +orderpath(fn, 5, "five"); +orderdpath(fn, 5, "five"); +orderspath(fn, 5, "five"); + +orderpath(fn, 6, "six"); +orderdpath(fn, 6, "six"); +orderspath(fn, 6, "six"); + +orderpath(fn, 7, "seven"); +orderdpath(fn, 7, "seven"); +orderspath(fn, 7, "seven"); + +orderpath(fn, 8, "eight"); +orderdpath(fn, 8, "eight"); +orderspath(fn, 8, "eight"); + +orderpath(fn, 9, "nine"); +orderdpath(fn, 9, "nine"); +orderspath(fn, 9, "nine"); + +orderpath(fn, 10, "ten"); +orderdpath(fn, 10, "ten"); +orderspath(fn, 10, "ten"); + +fclose(fn); + + diff --git a/buch/chapters/040-rekursion/images/order.pdf b/buch/chapters/040-rekursion/images/order.pdf new file mode 100644 index 0000000..cc175a9 Binary files /dev/null and b/buch/chapters/040-rekursion/images/order.pdf differ diff --git a/buch/chapters/040-rekursion/images/order.tex b/buch/chapters/040-rekursion/images/order.tex new file mode 100644 index 0000000..9a2511c --- /dev/null +++ b/buch/chapters/040-rekursion/images/order.tex @@ -0,0 +1,125 @@ +% +% order.tex -- Verteilungsfunktion für Ordnungsstatistik +% +% (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +% +\documentclass[tikz]{standalone} +\usepackage{amsmath} +\usepackage{times} +\usepackage{txfonts} +\usepackage{pgfplots} +\usepackage{csvsimple} +\usetikzlibrary{arrows,intersections,math} +\begin{document} +\def\skala{8} +\definecolor{darkgreen}{rgb}{0,0.6,0} + +\def\n{10} +\def\E#1#2{ + \draw[color=#2] + ({\dx*#1/(\n+1)},{-0.1/\skala}) -- ({\dx*#1/(\n+1)},{4.4*\dy}); + \node[color=#2] at ({\dx*#1/(\n+1)},{3.2*\dy}) + [rotate=90,above right] {$k=#1$}; +} +\def\var#1#2{ + \pgfmathparse{\dx*sqrt(#1*(\n-#1+1)/((\n+1)*(\n+1)*(\n+2)))} + \xdef\var{\pgfmathresult} + \fill[color=#2,opacity=0.5] + ({\dx*#1/(\n+1)-\var},0) rectangle ({\dx*#1/(\n+1)+\var},{4.4*\dy}); +} + +\input{orderpath.tex} +\begin{tikzpicture}[>=latex,thick,scale=\skala] + +\def\dx{1.6} +\def\dy{0.5} + +\def\pfad#1#2{ +\draw[color=#2,line width=1.4pt] ({-0.1/\skala},0) + -- + #1 + -- + ({1*\dx+0.1/\skala},0.5); +} + +\pfad{\orderzero}{darkgreen!20} +\pfad{\orderone}{darkgreen!20} +\pfad{\ordertwo}{darkgreen!20} +\pfad{\orderthree}{darkgreen!20} +\pfad{\orderfour}{darkgreen!20} +\pfad{\orderfive}{darkgreen!20} +\pfad{\ordersix}{darkgreen!20} +\pfad{\ordereight}{darkgreen!20} +\pfad{\ordernine}{darkgreen!20} +\pfad{\orderten}{darkgreen!20} +\pfad{\orderseven}{darkgreen} + +\draw[->] ({-0.1/\skala},0) -- ({1.03*\dx},0) coordinate[label={$x$}]; +\draw[->] (0,{-0.1/\skala}) -- (0,0.6) coordinate[label={right:$F(X)$}]; +\foreach \x in {0,0.2,0.4,0.6,0.8,1}{ + \draw ({\x*\dx},{-0.1/\skala}) -- ({\x*\dx},{0.1/\skala}); + \node at ({\x*\dx},{-0.1/\skala}) [below] {$\x$}; +} +\foreach \y in {0.5,1}{ + \draw ({-0.1/\skala},{\y*\dy}) -- ({0.1/\skala},{\y*\dy}); + \node at ({-0.1/\skala},{\y*\dy}) [left] {$\y$}; +} + +\node[color=darkgreen] at (0.65,{0.5*\dy}) [above,rotate=55] {$k=7$}; + +\begin{scope}[yshift=-0.7cm] +\def\dy{0.125} + +\foreach \k in {1,2,3,4,5,6,8,9,10}{ + \E{\k}{blue!30} +} +\def\k{7} +\var{\k}{orange!40} +\node[color=blue] at ({\dx*\k/(\n+1)},{4.3*\dy}) [above] {$E(X_{7:n})$}; + +\def\pfad#1#2{ + \draw[color=#2,line width=1.4pt] ({-0.1/\skala},0) + -- + #1 + -- + ({1*\dx+0.1/\skala},0.0); +} + +\begin{scope} +\clip ({-0.1/\skala},{-0.1/\skala}) + rectangle ({1*\dx+0.1/\skala},{0.56+0.1/\skala}); + +\pfad{\orderdzero}{red!20} +\pfad{\orderdone}{red!20} +\pfad{\orderdtwo}{red!20} +\pfad{\orderdthree}{red!20} +\pfad{\orderdfour}{red!20} +\pfad{\orderdfive}{red!20} +\pfad{\orderdsix}{red!20} +\pfad{\orderdeight}{red!20} +\pfad{\orderdnine}{red!20} +\pfad{\orderdten}{red!20} +\E{\k}{blue} +\pfad{\orderdseven}{red} + +\end{scope} + +\draw[->] ({-0.1/\skala},0) -- ({1.03*\dx},0) coordinate[label={$x$}]; +\draw[->] (0,{-0.1/\skala}) -- (0,0.6) coordinate[label={right:$\varphi(X)$}]; +\foreach \x in {0,0.2,0.4,0.6,0.8,1}{ + \draw ({\x*\dx},{-0.1/\skala}) -- ({\x*\dx},{0.1/\skala}); + \node at ({\x*\dx},{-0.1/\skala}) [below] {$\x$}; +} +\foreach \y in {1,2,3,4}{ + \draw ({-0.1/\skala},{\y*\dy}) -- ({0.1/\skala},{\y*\dy}); + \node at ({-0.1/\skala},{\y*\dy}) [left] {$\y$}; +} + +\node[color=red] at ({0.67*\dx},{2.7*\dy}) [above] {$k=7$}; + + +\end{scope} + +\end{tikzpicture} +\end{document} + diff --git a/buch/chapters/070-orthogonalitaet/Makefile.inc b/buch/chapters/070-orthogonalitaet/Makefile.inc index 48e5356..286ab2e 100644 --- a/buch/chapters/070-orthogonalitaet/Makefile.inc +++ b/buch/chapters/070-orthogonalitaet/Makefile.inc @@ -13,4 +13,5 @@ CHAPTERFILES = $(CHAPTERFILES) \ chapters/070-orthogonalitaet/jacobi.tex \ chapters/070-orthogonalitaet/sturm.tex \ chapters/070-orthogonalitaet/gaussquadratur.tex \ + chapters/070-orthogonalitaet/uebungsaufgaben/701.tex \ chapters/070-orthogonalitaet/chapter.tex diff --git a/buch/chapters/070-orthogonalitaet/chapter.tex b/buch/chapters/070-orthogonalitaet/chapter.tex index 5ebb795..4756844 100644 --- a/buch/chapters/070-orthogonalitaet/chapter.tex +++ b/buch/chapters/070-orthogonalitaet/chapter.tex @@ -25,7 +25,7 @@ \rhead{Übungsaufgaben} \aufgabetoplevel{chapters/070-orthogonalitaet/uebungsaufgaben} \begin{uebungsaufgaben} -%\uebungsaufgabe{0} +\uebungsaufgabe{701} %\uebungsaufgabe{1} \end{uebungsaufgaben} diff --git a/buch/chapters/070-orthogonalitaet/uebungsaufgaben/701.tex b/buch/chapters/070-orthogonalitaet/uebungsaufgaben/701.tex new file mode 100644 index 0000000..dad489f --- /dev/null +++ b/buch/chapters/070-orthogonalitaet/uebungsaufgaben/701.tex @@ -0,0 +1,137 @@ +Für Funktionen auf dem Interval $(-\frac{\pi}2,\frac{\pi}2)$ ist +\[ +\langle f,g\rangle += +\frac12\int_{-\frac{\pi}2}^{\frac{\pi}2} f(x)g(x)\cos x\,dx +\] +ein Skalarprodukt. +Bestimmen Sie bezüglich dieses Skalarproduktes orthogonale Polynome +bis zum Grad $2$. + +\begin{hinweis} +Verwenden Sie +\begin{align*} +\int_{-\frac{\pi}2}^{\frac{\pi}2} 1\cos x\,dx +&= +1, +& +\int_{-\frac{\pi}2}^{\frac{\pi}2} x^2\cos x\,dx +&= +\frac{\pi^2-8}{2}, +& +\int_{-\frac{\pi}2}^{\frac{\pi}2} x^4\cos x\,dx +&= +\frac{\pi^4-48\pi^2+384}{8}. +\end{align*} +\end{hinweis} + +\begin{loesung} +Wir müssen den Gram-Schmidt-Orthogonalisierungsprozess für die +Polynome $f_0(x)=1$, $f_1(x)=x$ und $f_2(x)=x^2$ durchführen. +Zunächst halten wir fest, dass +\[ +\langle f_0,f_0\rangle += +\frac12 +\int_{-\frac{\pi}2}^{\frac{\pi}2} \cos x\,dx += +1, +\] +das Polynom $g_0(x)=f_0(x)$ ist hat also Norm $1$. + +Ein dazu orthogonales Polynom ist +\( +f_1(x) - \langle g_0,f_1\rangle g_0(x), +\) +wir müssen also das Skalarprodukt +\[ +\langle g_0,f_1\rangle += +\frac{1}{2} +\int_{-\frac{\pi}2}^{\frac{\pi}2} +x\cos x\,dx +\] +bestimmen. +Es verschwindet, weil die Funktion $x\cos x$ ungerade ist. +Somit ist die Funktion $f_1(x)=x$ orthogonal zu $f_0(x)=1$, um sie auch zu +normieren berechnen wir das Integral +\[ +\| f_1\|^2 += +\frac12\int_{-\frac{\pi}2}^{\frac{\pi}2} x^2\cos x\,dx += +\frac{\pi^2-8}{4}, +\] +und +\[ +g_1(x) += +\frac{2}{\sqrt{\pi^2-8}} x. +\] + +Zur Berechnung von $g_2$ müssen wir die Skalarprodukte +\begin{align*} +\langle g_0,f_2\rangle +&= +\frac{1}{2} +\int_{-\frac{\pi}2}^{\frac{\pi}2} +x^2 +\cos x +\,dx += +\frac{\pi^2-8}{4} +\\ +\langle g_1,f_2\rangle +&= +\frac{1}{2} +\int_{-\frac{\pi}2}^{\frac{\pi}2} +\frac{2}{\sqrt{\pi^2-8}} +x +\cdot x^2 +\cos x +\,dx += +0 +\end{align*} +bestimmen. +Damit wird das dritte Polynom +\[ +f_2(x) +- g_0(x)\langle g_0,f_2\rangle +- g_1(x)\langle g_1,f_2\rangle += +x^2 - \frac{\pi^2-8}{4}, +\] +welches bereits orthogonal ist zu $g_0$ und $g_1$. +Wir können auch noch erreichen, obwohl das nicht verlangt war, +dass es normiert ist, indem wir die Norm berechnen: +\[ +\left\| x^2-\frac{\pi^2-8}{4} \right\|^2 += +\frac12 +\int_{-\frac{\pi}2}^{\frac{\pi}2} +\biggl(x^2-\frac{\pi^2-8}{4}\biggr)^2 +\cos x\,dx += +20-2\pi^2 +\] +woraus sich +\[ +g_2(x) += +\frac{1}{\sqrt{20-2\pi^2}} +\biggl( +x^2 - \frac{\pi^2-8}{4} +\biggr). +\] +Damit haben wir die ersten drei bezüglich des obigen Skalarproduktes +orthogonalen Polynome +\begin{align*} +g_0(x)&=1, +& +g_1(x)&=\frac{2x}{\sqrt{\pi^2-8}}, +& +g_2(x)&=\frac{1}{\sqrt{20-2\pi^2}}\biggl(x^2-\frac{\pi^2-8}{4}\biggr) +\end{align*} +gefunden. +\end{loesung} diff --git a/buch/chapters/090-pde/Makefile.inc b/buch/chapters/090-pde/Makefile.inc index a9ef74a..c64af06 100644 --- a/buch/chapters/090-pde/Makefile.inc +++ b/buch/chapters/090-pde/Makefile.inc @@ -10,4 +10,5 @@ CHAPTERFILES = $(CHAPTERFILES) \ chapters/090-pde/rechteck.tex \ chapters/090-pde/kreis.tex \ chapters/090-pde/kugel.tex \ + chapters/090-pde/uebungsaufgaben/901.tex \ chapters/090-pde/chapter.tex diff --git a/buch/chapters/090-pde/chapter.tex b/buch/chapters/090-pde/chapter.tex index db909ee..a393da5 100644 --- a/buch/chapters/090-pde/chapter.tex +++ b/buch/chapters/090-pde/chapter.tex @@ -21,11 +21,11 @@ deren Lösungen spezielle Funktionen sind. \input{chapters/090-pde/kreis.tex} \input{chapters/090-pde/kugel.tex} -%\section*{Übungsaufgaben} -%\rhead{Übungsaufgaben} -%\aufgabetoplevel{chapters/020-exponential/uebungsaufgaben} -%\begin{uebungsaufgaben} -%\uebungsaufgabe{0} +\section*{Übungsaufgaben} +\rhead{Übungsaufgaben} +\aufgabetoplevel{chapters/090-pde/uebungsaufgaben} +\begin{uebungsaufgaben} +\uebungsaufgabe{901} %\uebungsaufgabe{1} -%\end{uebungsaufgaben} +\end{uebungsaufgaben} diff --git a/buch/chapters/090-pde/kreis.tex b/buch/chapters/090-pde/kreis.tex index a24b6bb..b4ce8d7 100644 --- a/buch/chapters/090-pde/kreis.tex +++ b/buch/chapters/090-pde/kreis.tex @@ -120,7 +120,7 @@ für $\Phi(\varphi)$. Die Gleichung für $\Phi$ hat für $\mu\ne 0$ die Lösungen \begin{align*} \Phi(\varphi) &= \cos\mu\varphi -\text{und}\qquad +&&\text{und}& \Phi(\varphi) &= \sin\mu\varphi. \end{align*} Die Lösung muss aber auch stetig sein, d.~h.~es muss $\Phi(0)=\Phi(2\pi)$ diff --git a/buch/chapters/090-pde/uebungsaufgaben/901.tex b/buch/chapters/090-pde/uebungsaufgaben/901.tex new file mode 100644 index 0000000..67fa8e5 --- /dev/null +++ b/buch/chapters/090-pde/uebungsaufgaben/901.tex @@ -0,0 +1,82 @@ +Die Differentialgleichung +\begin{equation} +\frac{\partial u}{\partial t} = \kappa \frac{\partial^2 u}{\partial x^2} +\qquad +\text{im Gebiet} +\qquad +(t,x)\in \Omega=\mathbb{R}^+\times (0,l) +\label{505:waermeleitungsgleichung} +\end{equation} +beschreibt die Änderung der Temperatur eines Stabes der Länge $l$. +Die homogene Randbedingung +\begin{equation} +u(t,0)= +u(t,l)=0 +\label{505:homogene-randbedingung} +\end{equation} +besagt, dass der Stab an seinen Enden auf Temperatur $0$ gehalten. +Zur Lösung dieser Differentialgleichung muss auch die Temperatur +zur Zeit $t=0$ in Form einer Randbedingung +\[ +u(0,x) = T_0(x) +\] +gegeben sein. +Führen Sie Separation für die +Differentialgleichung~\eqref{505:waermeleitungsgleichung} +durch und bestimmen Sie die zulässigen Werte der Separationskonstanten. + +\begin{loesung} +Man verwendet den Ansatz $u(t,x)= T(t)\cdot X(x)$ und setzt diesen +in die Differentialgleichung ein, die dadurch zu +\[ +T'(t)X(x) = \kappa T(t) X''(x) +\] +wird. +Division durch $T(t)X(x)$ wird dies zu +\[ +\frac{T'(t)}{T(t)} += +\kappa +\frac{X''(x)}{X(x)}. +\] +Da die linke Seite nur von $t$ abhängt, die rechte aber nur von $x$, müssen +beide Seiten konstant sein. +Wir bezeichnen die Konstante mit $-\lambda^2$, so dass wir die beiden +gewöhnlichen Differentialgleichungen +\begin{align*} +\frac{1}{\kappa} +\frac{T'(t)}{T(t)}&=-\lambda^2 +& +\frac{X''(x)}{X(x)}&=-\lambda^2 +\\ +T'(t)&=-\lambda^2\kappa T(t) +& +X''(x) &= -\lambda^2 X(x) +\intertext{welche die Lösungen} +T(t)&=Ce^{-\lambda^2\kappa t} +& +X(x)&= A\cos\lambda x + B\sin\lambda x +\end{align*} +haben. +Die Lösung $X(x)$ muss aber auch die homogene Randbedingung +\eqref{505:homogene-randbedingung} erfüllen. +Setzt man $x=0$ und $x=l$ ein, folgt +\begin{align*} +0 = X(0)&=A\cos 0 + B\sin 0 = A +& +0 = X(l)&=B\sin \lambda l, +\end{align*} +woraus man schliessen kann, dass $\lambda l$ ein ganzzahliges +Vielfaches von $\pi$ ist, wir schreiben $\lambda l = k\pi$ oder +\[ +\lambda = \frac{k\pi}{l}. +\] +Damit sind die möglichen Werte $\lambda$ bestimmt und man kann jetzt +auch die möglichen Lösungen aufschreiben, sie sind +\[ +u(t,x) += +\sum_{k=1}^\infty b_k e^{-k^2\pi^2\kappa t/l^2}\sin\frac{k\pi x}{l}. +\qedhere +\] +\end{loesung} diff --git a/buch/papers/dreieck/images/beta.pdf b/buch/papers/dreieck/images/beta.pdf index c3ab4f6..cd5ed80 100644 Binary files a/buch/papers/dreieck/images/beta.pdf and b/buch/papers/dreieck/images/beta.pdf differ diff --git a/buch/papers/dreieck/images/beta.tex b/buch/papers/dreieck/images/beta.tex index 50509ee..f0ffdf0 100644 --- a/buch/papers/dreieck/images/beta.tex +++ b/buch/papers/dreieck/images/beta.tex @@ -23,7 +23,8 @@ \definecolor{coloreight}{rgb}{0.0,0.8,0.8} \definecolor{colornine}{rgb}{0.0,0.8,0.2} \definecolor{colorten}{rgb}{0.2,0.4,0.0} -\definecolor{coloreleven}{rgb}{1.0,0.8,0.4} +\definecolor{coloreleven}{rgb}{0.6,1.0,0.0} +\definecolor{colortwelve}{rgb}{1.0,0.8,0.4} \def\achsen{ \foreach \x in {0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9}{ @@ -47,24 +48,24 @@ } \def\farbcoord#1#2{ - ({\dx*(0.7+((#1-1)/4)*0.27)},{\dx*(0.15+((#2-1)/4)*0.27)}) + ({\dx*(0.63+((#1)/5)*0.27)},{\dx*(0.18+((#2)/5)*0.27)}) } \def\farbviereck{ - \foreach \x in {1,2,3,4,5}{ - \draw[color=gray!30] \farbcoord{\x}{1} -- \farbcoord{\x}{5}; - \draw[color=gray!30] \farbcoord{1}{\x} -- \farbcoord{5}{\x}; + \foreach \x in {1,2,3,4}{ + \draw[color=gray!30] \farbcoord{\x}{0} -- \farbcoord{\x}{4}; + \draw[color=gray!30] \farbcoord{0}{\x} -- \farbcoord{4}{\x}; } - \draw[->] \farbcoord{1}{1} -- \farbcoord{5.4}{1} + \draw[->] \farbcoord{0}{0} -- \farbcoord{4.4}{0} coordinate[label={$a$}]; - \draw[->] \farbcoord{1}{1} -- \farbcoord{1}{5.4} + \draw[->] \farbcoord{0}{0} -- \farbcoord{0}{4.4} coordinate[label={left: $b$}]; - \foreach \x in {1,2,3,4,5}{ - \node[color=gray] at \farbcoord{5}{\x} [right] {\tiny $b=\x$}; - \fill[color=white,opacity=0.7] - \farbcoord{(\x-0.1)}{4.3} - rectangle - \farbcoord{(\x+0.1)}{5}; - \node[color=gray] at \farbcoord{\x}{5} [left,rotate=90] + \foreach \x in {1,2,3,4}{ + \node[color=gray] at \farbcoord{4}{\x} [right] {\tiny $b=\x$}; + %\fill[color=white,opacity=0.7] + % \farbcoord{(\x-0.1)}{3.3} + % rectangle + % \farbcoord{(\x+0.1)}{4}; + \node[color=gray] at \farbcoord{\x}{4} [right,rotate=90] {\tiny $a=\x$}; } } @@ -74,23 +75,26 @@ \begin{tikzpicture}[>=latex,thick,scale=\skala] -\def\dx{1} +\def\dx{1.1} \def\dy{0.1} \def\opa{0.1} -\def\betamax{4.2} - -\fill[color=colorone,opacity=\opa] (0,0) -- \betaaa -- (\dx,0) -- cycle; -\fill[color=colortwo,opacity=\opa] (0,0) -- \betabb -- (\dx,0) -- cycle; -\fill[color=colorthree,opacity=\opa] (0,0) -- \betacc -- (\dx,0) -- cycle; -\fill[color=colorfour,opacity=\opa] (0,0) -- \betadd -- (\dx,0) -- cycle; -\fill[color=colorfive,opacity=\opa] (0,0) -- \betaee -- (\dx,0) -- cycle; -\fill[color=colorsix,opacity=\opa] (0,0) -- \betaff -- (\dx,0) -- cycle; -\fill[color=colorseven,opacity=\opa] (0,0) -- \betagg -- (\dx,0) -- cycle; -\fill[color=coloreight,opacity=\opa] (0,0) -- \betahh -- (\dx,0) -- cycle; -\fill[color=colornine,opacity=\opa] (0,0) -- \betaii -- (\dx,0) -- cycle; -\fill[color=colorten,opacity=\opa] (0,0) -- \betajj -- (\dx,0) -- cycle; +\def\betamax{4.9} + +\begin{scope} +\clip (0,0) rectangle ({1*\dx},{\betamax*\dy}); +\fill[color=colorone,opacity=\opa] (0,0) -- \betaaa -- (\dx,0) -- cycle; +\fill[color=colortwo,opacity=\opa] (0,0) -- \betabb -- (\dx,0) -- cycle; +\fill[color=colorthree,opacity=\opa] (0,0) -- \betacc -- (\dx,0) -- cycle; +\fill[color=colorfour,opacity=\opa] (0,0) -- \betadd -- (\dx,0) -- cycle; +\fill[color=colorfive,opacity=\opa] (0,0) -- \betaee -- (\dx,0) -- cycle; +\fill[color=colorsix,opacity=\opa] (0,0) -- \betaff -- (\dx,0) -- cycle; +\fill[color=colorseven,opacity=\opa] (0,0) -- \betagg -- (\dx,0) -- cycle; +\fill[color=coloreight,opacity=\opa] (0,0) -- \betahh -- (\dx,0) -- cycle; +\fill[color=colornine,opacity=\opa] (0,0) -- \betaii -- (\dx,0) -- cycle; +\fill[color=colorten,opacity=\opa] (0,0) -- \betajj -- (\dx,0) -- cycle; \fill[color=coloreleven,opacity=\opa] (0,0) -- \betakk -- (\dx,0) -- cycle; +\fill[color=colortwelve,opacity=\opa] (0,0) -- \betall -- (\dx,0) -- cycle; \draw[color=colorone] \betaaa; \draw[color=colortwo] \betabb; @@ -103,11 +107,15 @@ \draw[color=colornine] \betaii; \draw[color=colorten] \betajj; \draw[color=coloreleven] \betakk; +\draw[color=colortwelve] \betall; + +\end{scope} \achsen \farbviereck +\farbpunkt{\alphatwelve}{\betatwelve}{colortwelve} \farbpunkt{\alphaeleven}{\betaeleven}{coloreleven} \farbpunkt{\alphaten}{\betaten}{colorten} \farbpunkt{\alphanine}{\betanine}{colornine} @@ -124,88 +132,102 @@ \def\betamax{4.9} \begin{scope}[yshift=-0.6cm] -\fill[color=colorone,opacity=\opa] (0,0) -- \betaaa -- (\dx,0) -- cycle; -\fill[color=colortwo,opacity=\opa] (0,0) -- \betaab -- (\dx,0) -- cycle; -\fill[color=colorthree,opacity=\opa] (0,0) -- \betaac -- (\dx,0) -- cycle; -\fill[color=colorfour,opacity=\opa] (0,0) -- \betaad -- (\dx,0) -- cycle; -\fill[color=colorfive,opacity=\opa] (0,0) -- \betaae -- (\dx,0) -- cycle; -\fill[color=colorsix,opacity=\opa] (0,0) -- \betaaf -- (\dx,0) -- cycle; -\fill[color=colorseven,opacity=\opa] (0,0) -- \betaag -- (\dx,0) -- cycle; -\fill[color=coloreight,opacity=\opa] (0,0) -- \betaah -- (\dx,0) -- cycle; -\fill[color=colornine,opacity=\opa] (0,0) -- \betaai -- (\dx,0) -- cycle; -\fill[color=colorten,opacity=\opa] (0,0) -- \betaaj -- (\dx,0) -- cycle; -\fill[color=coloreleven,opacity=\opa] (0,0) -- \betaak -- (\dx,0) -- cycle; -\draw[color=colorone] \betaaa; -\draw[color=colortwo] \betaab; -\draw[color=colorthree] \betaac; -\draw[color=colorfour] \betaad; -\draw[color=colorfive] \betaae; -\draw[color=colorsix] \betaaf; -\draw[color=colorseven] \betaag; -\draw[color=coloreight] \betaah; -\draw[color=colornine] \betaai; -\draw[color=colorten] \betaaj; -\draw[color=coloreleven] \betaak; +\begin{scope} +\clip (0,0) rectangle ({1*\dx},{\betamax*\dy}); +\fill[color=colorone,opacity=\opa] (0,0) -- \betaea -- (\dx,0) -- cycle; +\fill[color=colortwo,opacity=\opa] (0,0) -- \betaeb -- (\dx,0) -- cycle; +\fill[color=colorthree,opacity=\opa] (0,0) -- \betaec -- (\dx,0) -- cycle; +\fill[color=colorfour,opacity=\opa] (0,0) -- \betaed -- (\dx,0) -- cycle; +\fill[color=colorfive,opacity=\opa] (0,0) -- \betaee -- (\dx,0) -- cycle; +\fill[color=colorsix,opacity=\opa] (0,0) -- \betaef -- (\dx,0) -- cycle; +\fill[color=colorseven,opacity=\opa] (0,0) -- \betaeg -- (\dx,0) -- cycle; +\fill[color=coloreight,opacity=\opa] (0,0) -- \betaeh -- (\dx,0) -- cycle; +\fill[color=colornine,opacity=\opa] (0,0) -- \betaei -- (\dx,0) -- cycle; +\fill[color=colorten,opacity=\opa] (0,0) -- \betaej -- (\dx,0) -- cycle; +\fill[color=coloreleven,opacity=\opa] (0,0) -- \betaek -- (\dx,0) -- cycle; +\fill[color=colortwelve,opacity=\opa] (0,0) -- \betael -- (\dx,0) -- cycle; + +\draw[color=colorone] \betaea; +\draw[color=colortwo] \betaeb; +\draw[color=colorthree] \betaec; +\draw[color=colorfour] \betaed; +\draw[color=colorfive] \betaee; +\draw[color=colorsix] \betaef; +\draw[color=colorseven] \betaeg; +\draw[color=coloreight] \betaeh; +\draw[color=colornine] \betaei; +\draw[color=colorten] \betaej; +\draw[color=coloreleven] \betaek; +\draw[color=colortwelve] \betael; +\end{scope} \achsen \farbviereck -\farbpunkt{\alphaone}{\betaeleven}{coloreleven} -\farbpunkt{\alphaone}{\betaten}{colorten} -\farbpunkt{\alphaone}{\betanine}{colornine} -\farbpunkt{\alphaone}{\betaeight}{coloreight} -\farbpunkt{\alphaone}{\betaseven}{colorseven} -\farbpunkt{\alphaone}{\betasix}{colorsix} -\farbpunkt{\alphaone}{\betafive}{colorfive} -\farbpunkt{\alphaone}{\betafour}{colorfour} -\farbpunkt{\alphaone}{\betathree}{colorthree} -\farbpunkt{\alphaone}{\betatwo}{colortwo} -\farbpunkt{\alphaone}{\betaone}{colorone} +\farbpunkt{\alphafive}{\betatwelve}{colortwelve} +\farbpunkt{\alphafive}{\betaeleven}{coloreleven} +\farbpunkt{\alphafive}{\betaten}{colorten} +\farbpunkt{\alphafive}{\betanine}{colornine} +\farbpunkt{\alphafive}{\betaeight}{coloreight} +\farbpunkt{\alphafive}{\betaseven}{colorseven} +\farbpunkt{\alphafive}{\betasix}{colorsix} +\farbpunkt{\alphafive}{\betafive}{colorfive} +\farbpunkt{\alphafive}{\betafour}{colorfour} +\farbpunkt{\alphafive}{\betathree}{colorthree} +\farbpunkt{\alphafive}{\betatwo}{colortwo} +\farbpunkt{\alphafive}{\betaone}{colorone} \end{scope} \begin{scope}[yshift=-1.2cm] -\fill[color=colorone,opacity=\opa] (0,0) -- \betaak -- (\dx,0) -- cycle; -\fill[color=colortwo,opacity=\opa] (0,0) -- \betabk -- (\dx,0) -- cycle; -\fill[color=colorthree,opacity=\opa] (0,0) -- \betack -- (\dx,0) -- cycle; -\fill[color=colorfour,opacity=\opa] (0,0) -- \betadk -- (\dx,0) -- cycle; -\fill[color=colorfive,opacity=\opa] (0,0) -- \betaek -- (\dx,0) -- cycle; -\fill[color=colorsix,opacity=\opa] (0,0) -- \betafk -- (\dx,0) -- cycle; -\fill[color=colorseven,opacity=\opa] (0,0) -- \betagk -- (\dx,0) -- cycle; -\fill[color=coloreight,opacity=\opa] (0,0) -- \betahk -- (\dx,0) -- cycle; -\fill[color=colornine,opacity=\opa] (0,0) -- \betaik -- (\dx,0) -- cycle; -\fill[color=colorten,opacity=\opa] (0,0) -- \betajk -- (\dx,0) -- cycle; -\fill[color=coloreleven,opacity=\opa] (0,0) -- \betakk -- (\dx,0) -- cycle; -\draw[color=colorone] \betaak; -\draw[color=colortwo] \betabk; -\draw[color=colorthree] \betack; -\draw[color=colorfour] \betadk; -\draw[color=colorfive] \betaek; -\draw[color=colorsix] \betafk; -\draw[color=colorseven] \betagk; -\draw[color=coloreight] \betahk; -\draw[color=colornine] \betaik; -\draw[color=colorten] \betajk; -\draw[color=coloreleven] \betakk; +\begin{scope} +\clip (0,0) rectangle ({1*\dx},{\betamax*\dy}); +\fill[color=colorone,opacity=\opa] (0,0) -- \betaal -- (\dx,0) -- cycle; +\fill[color=colortwo,opacity=\opa] (0,0) -- \betabl -- (\dx,0) -- cycle; +\fill[color=colorthree,opacity=\opa] (0,0) -- \betacl -- (\dx,0) -- cycle; +\fill[color=colorfour,opacity=\opa] (0,0) -- \betadl -- (\dx,0) -- cycle; +\fill[color=colorfive,opacity=\opa] (0,0) -- \betael -- (\dx,0) -- cycle; +\fill[color=colorsix,opacity=\opa] (0,0) -- \betafl -- (\dx,0) -- cycle; +\fill[color=colorseven,opacity=\opa] (0,0) -- \betagl -- (\dx,0) -- cycle; +\fill[color=coloreight,opacity=\opa] (0,0) -- \betahl -- (\dx,0) -- cycle; +\fill[color=colornine,opacity=\opa] (0,0) -- \betail -- (\dx,0) -- cycle; +\fill[color=colorten,opacity=\opa] (0,0) -- \betajl -- (\dx,0) -- cycle; +\fill[color=coloreleven,opacity=\opa] (0,0) -- \betakl -- (\dx,0) -- cycle; +\fill[color=colortwelve,opacity=\opa] (0,0) -- \betall -- (\dx,0) -- cycle; + +\draw[color=colorone] \betaal; +\draw[color=colortwo] \betabl; +\draw[color=colorthree] \betacl; +\draw[color=colorfour] \betadl; +\draw[color=colorfive] \betael; +\draw[color=colorsix] \betafl; +\draw[color=colorseven] \betagl; +\draw[color=coloreight] \betahl; +\draw[color=colornine] \betail; +\draw[color=colorten] \betajl; +\draw[color=coloreleven] \betakl; +\draw[color=colortwelve] \betall; +\end{scope} \achsen \farbviereck -\farbpunkt{\alphaeleven}{\betaeleven}{coloreleven} -\farbpunkt{\alphaten}{\betaeleven}{colorten} -\farbpunkt{\alphanine}{\betaeleven}{colornine} -\farbpunkt{\alphaeight}{\betaeleven}{coloreight} -\farbpunkt{\alphaseven}{\betaeleven}{colorseven} -\farbpunkt{\alphasix}{\betaeleven}{colorsix} -\farbpunkt{\alphafive}{\betaeleven}{colorfive} -\farbpunkt{\alphafour}{\betaeleven}{colorfour} -\farbpunkt{\alphathree}{\betaeleven}{colorthree} -\farbpunkt{\alphatwo}{\betaeleven}{colortwo} -\farbpunkt{\alphaone}{\betaeleven}{colorone} +\farbpunkt{\alphatwelve}{\betatwelve}{colortwelve} +\farbpunkt{\alphaeleven}{\betatwelve}{coloreleven} +\farbpunkt{\alphaten}{\betatwelve}{colorten} +\farbpunkt{\alphanine}{\betatwelve}{colornine} +\farbpunkt{\alphaeight}{\betatwelve}{coloreight} +\farbpunkt{\alphaseven}{\betatwelve}{colorseven} +\farbpunkt{\alphasix}{\betatwelve}{colorsix} +\farbpunkt{\alphafive}{\betatwelve}{colorfive} +\farbpunkt{\alphafour}{\betatwelve}{colorfour} +\farbpunkt{\alphathree}{\betatwelve}{colorthree} +\farbpunkt{\alphatwo}{\betatwelve}{colortwo} +\farbpunkt{\alphaone}{\betatwelve}{colorone} \end{scope} diff --git a/buch/papers/dreieck/images/betadist.m b/buch/papers/dreieck/images/betadist.m index 9ff78ed..5b466a6 100644 --- a/buch/papers/dreieck/images/betadist.m +++ b/buch/papers/dreieck/images/betadist.m @@ -5,24 +5,32 @@ # global N; N = 201; -global n; -n = 11; +global nmin; +global nmax; +nmin = -4; +nmax = 7; +n = nmax - nmin + 1 +A = 3; -t = (0:n-1) / (n-1) -alpha = 1 + 4 * t.^2 +t = (nmin:nmax) / nmax; +alpha = 1 + A * t .* abs(t) +#alpha(1) = 0.01; #alpha = [ 1, 1.03, 1.05, 1.1, 1.25, 1.5, 2, 2.5, 3, 4, 5 ]; beta = alpha; names = [ "one"; "two"; "three"; "four"; "five"; "six"; "seven"; "eight"; - "nine"; "ten"; "eleven" ] + "nine"; "ten"; "eleven"; "twelve" ] function retval = Beta(a, b, x) retval = x^(a-1) * (1-x)^(b-1) / beta(a, b); + if (retval > 100) + retval = 100 + end end function plotbeta(fn, a, b, name) global N; - fprintf(fn, "\\def\\beta%s{\n", name); + fprintf(fn, "\\def\\beta%s{\n", strtrim(name)); fprintf(fn, "\t({%.4f*\\dx},{%.4f*\\dy})", 0, Beta(a, b, 0)); for x = (1:N-1)/(N-1) X = (1-cos(pi * x))/2; @@ -35,8 +43,8 @@ end fn = fopen("betapaths.tex", "w"); for i = (1:n) - fprintf(fn, "\\def\\alpha%s{%f}\n", names(i,:), alpha(i)); - fprintf(fn, "\\def\\beta%s{%f}\n", names(i,:), beta(i)); + fprintf(fn, "\\def\\alpha%s{%f}\n", strtrim(names(i,:)), alpha(i)); + fprintf(fn, "\\def\\beta%s{%f}\n", strtrim(names(i,:)), beta(i)); end for i = (1:n) -- cgit v1.2.1 From 9f0ff73d26b3e096d848abfd20cd124433e2a4a7 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Sun, 13 Mar 2022 23:27:43 +0100 Subject: cleanup --- buch/aufgaben2.tex | 4 ---- 1 file changed, 4 deletions(-) (limited to 'buch') diff --git a/buch/aufgaben2.tex b/buch/aufgaben2.tex index f98562e..8073f26 100644 --- a/buch/aufgaben2.tex +++ b/buch/aufgaben2.tex @@ -4,10 +4,6 @@ % (c) 2022 Prof. Dr. Andreas Mueller, OST % -%\item -%\input chapters/40-eigenwerte/uebungsaufgaben/4004.tex -%\item -%\input chapters/40-eigenwerte/uebungsaufgaben/4005.tex \item \input{chapters/090-pde/uebungsaufgaben/901.tex} \item -- cgit v1.2.1 From 18e46179f2da76a3147d3f3b466206c6b5405859 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Mon, 14 Mar 2022 08:20:28 +0100 Subject: describe link between Jacobi-weights and Beta-distribution --- buch/chapters/070-orthogonalitaet/jacobi.tex | 22 ++ buch/chapters/070-orthogonalitaet/orthogonal.tex | 51 +++ buch/papers/dreieck/teil1.tex | 411 +---------------------- 3 files changed, 74 insertions(+), 410 deletions(-) (limited to 'buch') diff --git a/buch/chapters/070-orthogonalitaet/jacobi.tex b/buch/chapters/070-orthogonalitaet/jacobi.tex index 042d466..f776c03 100644 --- a/buch/chapters/070-orthogonalitaet/jacobi.tex +++ b/buch/chapters/070-orthogonalitaet/jacobi.tex @@ -189,6 +189,28 @@ rechten Rand haben. \label{buch:orthogonal:fig:jacobi-parameter}} \end{figure} +\subsection{Jacobi-Gewichtsfunktion und Beta-Verteilung +\label{buch:orthogonal:subsection:beta-verteilung}} +Die Jacobi-Gewichtsfunktion entsteht aus der Wahrscheinlichkeitsdichte +der Beta-Verteilung, die in +Abschnitt~\ref{buch:rekursion:subsection:beta-verteilung} +eingeführt wurde mit Hilfe der Variablen-Transformation $x = 2t-1$ +oder $t=(x+1)/2$. +Das Integral mit der Jacobi-Gewichtsfunktion $w^{(\alpha,\beta)}(x)$ +kann damit umgeformt werden in +\[ +\int_{-1}^1 +f(x)\,w^{(\alpha,\beta)}(x)\,dx += +\int_0^1 +f(2t-1) w^{(\alpha,\beta)}(2t-1)\,2\,dt += +\int_0^1 +f(2t-1) +(1-(2t-1))^\alpha (1+(2t-1))^\beta +\,2\,dt +\] + % % % diff --git a/buch/chapters/070-orthogonalitaet/orthogonal.tex b/buch/chapters/070-orthogonalitaet/orthogonal.tex index d06f46e..a84248a 100644 --- a/buch/chapters/070-orthogonalitaet/orthogonal.tex +++ b/buch/chapters/070-orthogonalitaet/orthogonal.tex @@ -737,6 +737,57 @@ rechten Rand haben. \label{buch:orthogonal:fig:jacobi-parameter}} \end{figure} +\subsubsection{Jacobi-Gewichtsfunktion und Beta-Verteilung +\label{buch:orthogonal:subsection:beta-verteilung}} +Die Jacobi-Gewichtsfunktion entsteht aus der Wahrscheinlichkeitsdichte +der Beta-Verteilung, die in +Abschnitt~\ref{buch:rekursion:subsection:beta-verteilung} +eingeführt wurde mit Hilfe der Variablen-Transformation $x = 2t-1$ +oder $t=(x+1)/2$. +Das Integral mit der Jacobi-Gewichtsfunktion $w^{(\alpha,\beta)}(x)$ +kann damit umgeformt werden in +\begin{align*} +\int_{-1}^1 +f(x)\,w^{(\alpha,\beta)}(x)\,dx +&= +\int_0^1 +f(2t-1) w^{(\alpha,\beta)}(2t-1)\,2\,dt +\\ +&= +\int_0^1 +f(2t-1) +(1-(2t-1))^\alpha (1+(2t-1))^\beta +\,2\,dt +\\ +&= +2^{\alpha+\beta+1} +\int_0^1 +f(2t-1) +\, +t^\beta +(1-t)^\alpha +\,dt +\\ +&= +2^{\alpha+\beta+1} +B(\alpha+1,\beta+1) +\int_0^1 +f(2t-1) +\, +\frac{ +t^\beta +(1-t)^\alpha +}{B(\alpha+1,\beta+1)} +\,dt. +\end{align*} +Auf der letzten Zeile steht ein Integral mit der Wahrscheinlichkeitsdichte +der Beta-Verteilung. +Orthogonale Funktionen bezüglich der Jacobischen Gewichtsfunktion +$w^{(\alpha,\beta)}$ werden mit der genannten Substitution also +zu orthogonalen Funktionen bezüglich der Beta-Verteilung mit +Parametern $\beta+1$ und $\alpha+1$. + + % % Tschebyscheff-Gewichtsfunktion % diff --git a/buch/papers/dreieck/teil1.tex b/buch/papers/dreieck/teil1.tex index 5e7090b..4abe2e1 100644 --- a/buch/papers/dreieck/teil1.tex +++ b/buch/papers/dreieck/teil1.tex @@ -5,416 +5,7 @@ % \section{Ordnungsstatistik und Beta-Funktion \label{dreieck:section:ordnungsstatistik}} -\rhead{Ordnungsstatistik und Beta-Funktion} -In diesem Abschnitt ist $X$ eine Zufallsvariable mit der Verteilungsfunktion -$F_X(x)$, und $X_i$, $1\le i\le n$ sei ein Stichprobe von unabhängigen -Zufallsvariablen, die wie $X$ verteilt sind. -Ziel ist, die Verteilungsfunktion und die Wahrscheinlichkeitsdichte -des grössten, zweitgrössten, $k$-t-grössten Wertes in der Stichprobe -zu finden. -Wir schreiben $[n]=\{1,\dots,n\}$ für die Menge der natürlichen -Zahlen von zwischen $1$ und $n$. +\rhead{} -\subsection{Verteilung von $\operatorname{max}(X_1,\dots,X_n)$ und -$\operatorname{min}(X_1,\dots,X_n)$ -\label{dreieck:subsection:minmax}} -Die Verteilungsfunktion von $\operatorname{max}(X_1,\dots,X_n)$ hat -den Wert -\begin{align*} -F_{\operatorname{max}(X_1,\dots,X_n)}(x) -&= -P(\operatorname{max}(X_1,\dots,X_n) \le x) -\\ -&= -P(X_1\le x\wedge \dots \wedge X_n\le x) -\\ -&= -P(X_1\le x) \cdot \ldots \cdot P(X_n\le x) -\\ -&= -P(X\le x)^n -= -F_X(x)^n. -\end{align*} -Für die Gleichverteilung ist -\[ -F_{\text{equi}}(x) -= -\begin{cases} -0&\qquad x< 0 -\\ -x&\qquad 0\le x\le 1 -\\ -1&\qquad 1 X_1\wedge \dots \wedge x > X_n) -\\ -&= -1- -(1-P(x\le X_1)) \cdot\ldots\cdot (1-P(x\le X_n)) -\\ -&= -1-(1-F_X(x))^n, -\end{align*} -Im Speziellen für im Intervall $[0,1]$ gleichverteilte $X_i$ ist die -Verteilungsfunktion des Minimums -\[ -F_{\operatorname{min}(X_1,\dots,X_n)}(x) -= -\begin{cases} -0 &\qquad x<0 \\ -1-(1-x)^n&\qquad 0\le x\le 1\\ -1 &\qquad 1 < x -\end{cases} -\] -mit Wahrscheinlichkeitsdichte -\[ -\varphi_{\operatorname{min}(X_1,\dots,X_n)} -= -\frac{d}{dx} -F_{\operatorname{min}(X_1,\dots,X_n)} -= -\begin{cases} -n(1-x)^{n-1}&\qquad 0\le x\le 1\\ -0 &\qquad \text{sonst} -\end{cases} -\] -und Erwartungswert -\begin{align*} -E(\operatorname{min}(X_1,\dots,X_n) -&= -\int_{-\infty}^\infty x\varphi_{\operatorname{min}(X_1,\dots,X_n)}(x)\,dx -= -\int_0^1 x\cdot n(1-x)^{n-1}\,dx -\\ -&= -\bigl[ -x(1-x)^n \bigr]_0^1 + \int_0^1 (1-x)^n\,dx -= -\biggl[ -- -\frac{1}{n+1} -(1-x)^{n+1} -\biggr]_0^1 -= -\frac{1}{n+1}. -\end{align*} -Es ergibt sich daraus als natürlich Verallgemeinerung die Frage nach -der Verteilung des zweitegrössten oder zweitkleinsten Wertes unter den -Werten $X_i$. - -\subsection{Der $k$-t-grösste Wert} -Sie wieder $X_i$ eine Stichprobe von $n$ unabhängigen wie $X$ verteilten -Zufallsvariablen. -Diese werden jetzt der Grösse nach sortiert, die sortierten Werte werden -mit -\[ -X_{1:n} \le X_{2:n} \le \dots \le X_{(n-1):n} \le X_{n:n} -\] -bezeichnet. -Die Grössen $X_{k:n}$ sind Zufallsvariablen, sie heissen die $k$-ten -Ordnungsstatistiken. -Die in Abschnitt~\ref{dreieck:subsection:minmax} behandelten Zufallsvariablen -$\operatorname{min}(X_1,\dots,X_n)$ -und -$\operatorname{max}(X_1,\dots,X_n)$ -sind die Fälle -\begin{align*} -X_{1:n} &= \operatorname{min}(X_1,\dots,X_n) \\ -X_{n:n} &= \operatorname{max}(X_1,\dots,X_n). -\end{align*} - -Um den Wert der Verteilungsfunktion von $X_{k:n}$ zu berechnen, müssen wir -die Wahrscheinlichkeit bestimmen, dass $k$ der $n$ Werte $X_i$ $x$ nicht -übersteigen. -Der $k$-te Wert $X_{k:n}$ übersteigt genau dann $x$ nicht, wenn -mindestens $k$ der Zufallswerte $X_i$ $x$ nicht übersteigen, also -\[ -P(X_{k:n} \le x) -= -P\left( -|\{i\in[n]\,|\, X_i\le x\}| \ge k -\right). -\] - -Das Ereignis $\{X_i\le x\}$ ist eine Bernoulli-Experiment, welches mit -Wahrscheinlichkeit $F_X(x)$ eintritt. -Die Anzahl der Zufallsvariablen $X_i$, die $x$ übertreffen, ist also -Binomialverteilt mit $p=F_X(x)$. -Damit haben wir gefunden, dass mit Wahrscheinlichkeit -\begin{equation} -F_{X_{k:n}}(x) -= -P(X_{k:n}\le x) -= -\sum_{i=k}^n \binom{n}{i}F_X(x)^i (1-F_X(x))^{n-i} -\label{dreieck:eqn:FXkn} -\end{equation} -mindestens $k$ der Zufallsvariablen den Wert $x$ überschreiten. - -\subsubsection{Wahrscheinlichkeitsdichte der Ordnungsstatistik} -Die Wahrscheinlichkeitsdichte der Ordnungsstatistik kann durch Ableitung -von \eqref{dreieck:eqn:FXkn} gefunden, werden, sie ist -\begin{align*} -\varphi_{X_{k:n}}(x) -&= -\frac{d}{dx} -F_{X_{k:n}}(x) -\\ -&= -\sum_{i=k}^n -\binom{n}{i} -\bigl( -iF_X(x)^{i-1}\varphi_X(x) (1-F_X(x))^{n-i} -- -F_X(x)^k -(n-i) -(1-F_X(x))^{n-i-1} -\varphi_X(x) -\bigr) -\\ -&= -\sum_{i=k}^n -\binom{n}{i} -\varphi_X(x) -F_X(x)^{i-1}(1-F_X(x))^{n-i-1} -\bigl( -iF_X(x)-(n-i)(1-F_X(x)) -\bigr) -\\ -&= -\varphi_X(x) -\biggl( -\sum_{i=k}^n i\binom{n}{i} F_X(x)^{i-1}(1-F_X(x))^{n-i} -- -\sum_{j=k}^n (n-j)\binom{n}{j} F_X(x)^{j}(1-F_X(x))^{n-j-1} -\biggr) -\\ -&= -\varphi_X(x) -\biggl( -\sum_{i=k}^n i\binom{n}{i} F_X(x)^{i-1}(1-F_X(x))^{n-i} -- -\sum_{i=k+1}^{n+1} (n-i+1)\binom{n}{i-1} F_X(x)^{i-1}(1-F_X(x))^{n-i} -\biggr) -\\ -&= -\varphi_X(x) -\biggl( -k\binom{n}{k}F_X(x)^{k-1}(1-F_X(x))^{n-k} -+ -\sum_{i=k+1}^{n+1} -\left( -i\binom{n}{i} -- -(n-i+1)\binom{n}{i-1} -\right) -F_X(x)^{i-1}(1-F_X(x))^{n-i} -\biggr) -\end{align*} -Mit den wohlbekannten Identitäten für die Binomialkoeffizienten -\begin{align*} -i\binom{n}{i} -- -(n-i+1)\binom{n}{i-1} -&= -n\binom{n-1}{i-1} -- -n -\binom{n-1}{i-1} -= -0 -\end{align*} -folgt jetzt -\begin{align*} -\varphi_{X_{k:n}}(x) -&= -\varphi_X(x)k\binom{n}{k} F_X(x)^{k-1}(1-F_X(x))^{n-k}(x). -\intertext{Im Speziellen für gleichverteilte Zufallsvariablen $X_i$ ist -} -\varphi_{X_{k:n}}(x) -&= -k\binom{n}{k} x^{k-1}(1-x)^{n-k}. -\end{align*} -Dies ist die Wahrscheinlichkeitsdichte einer Betaverteilung -\[ -\beta(k,n-k+1)(x) -= -\frac{1}{B(k,n-k+1)} -x^{k-1}(1-x)^{n-k}. -\] -Tatsächlich ist die Normierungskonstante -\begin{align} -\frac{1}{B(k,n-k+1)} -&= -\frac{\Gamma(n+1)}{\Gamma(k)\Gamma(n-k+1)} -= -\frac{n!}{(k-1)!(n-k)!}. -\label{dreieck:betaverteilung:normierung1} -\end{align} -Andererseits ist -\[ -k\binom{n}{k} -= -k\frac{n!}{k!(n-k)!} -= -\frac{n!}{(k-1)!(n-k)!}, -\] -in Übereinstimmung mit~\eqref{dreieck:betaverteilung:normierung1}. -Die Verteilungsfunktion und die Wahrscheinlichkeitsdichte der -Ordnungsstatistik sind in Abbildung~\ref{dreieck:fig:order} dargestellt. - -\begin{figure} -\centering -\includegraphics{papers/dreieck/images/order.pdf} -\caption{Verteilungsfunktion und Wahrscheinlichkeitsdichte der -Ordnungsstatistiken $X_{k:n}$ einer gleichverteilung Zuvallsvariable -mit $n=10$. -\label{dreieck:fig:order}} -\end{figure} - -\subsubsection{Erwartungswert} -Mit der Wahrscheinlichkeitsdichte kann man jetzt auch den Erwartungswerte -der $k$-ten Ordnungsstatistik bestimmen. -Die Rechnung ergibt: -\begin{align*} -E(X_{k:n}) -&= -\int_0^1 x\cdot k\binom{n}{k} x^{k-1}(1-x)^{n-k}\,dx -= -k -\binom{n}{k} -\int_0^1 -x^{k}(1-x)^{n-k}\,dx. -\intertext{Dies ist das Beta-Integral} -&= -k\binom{n}{k} -B(k+1,n-k+1) -\intertext{welches man durch Gamma-Funktionen bzw.~durch Fakultäten wie in} -&= -k\frac{n!}{k!(n-k)!} -\frac{\Gamma(k+1)\Gamma(n-k+1)}{n+2} -= -k\frac{n!}{k!(n-k)!} -\frac{k!(n-k)!}{(n+1)!} -= -\frac{k}{n+1} -\end{align*} -ausdrücken kann. -Die Erwartungswerte haben also regelmässige Abstände, sie sind in -Abbildung~\ref{dreieck:fig:order} als blaue vertikale Linien eingezeichnet. - -\subsubsection{Varianz} -Auch die Varianz lässt sich einfach berechnen, dazu muss zunächst -der Erwartungswert von $X_{k:n}^2$ bestimmt werden. -Er ist -\begin{align*} -E(X_{k:n}^2) -&= -\int_0^1 x^2\cdot k\binom{n}{k} x^{k-1}(1-x)^{n-k}\,dx -= -k -\binom{n}{k} -\int_0^1 -x^{k+1}(1-x)^{n-k}\,dx. -\intertext{Auch dies ist ein Beta-Integral, nämlich} -&= -k\binom{n}{k} -B(k+2,n-k+1) -= -k\frac{n!}{k!(n-k)!} -\frac{(k+1)!(n-k)!}{(n+2)!} -= -\frac{k(k+1)}{(n+1)(n+2)}. -\end{align*} -Die Varianz wird damit -\begin{align} -\operatorname{var}(X_{k:n}) -&= -E(X_{k:n}^2) - E(X_{k:n})^2 -\notag -\\ -& -= -\frac{k(k+1)}{(n+1)(n+2)}-\frac{k^2}{(n+1)^2} -= -\frac{k(k+1)(n+1)-k^2(n+2)}{(n+1)^2(n+2)} -= -\frac{k(n-k+1)}{(n+1)^2(n+2)}. -\label{dreieck:eqn:ordnungsstatistik:varianz} -\end{align} -In Abbildung~\ref{dreieck:fig:order} ist die Varianz der -Ordnungsstatistik $X_{k:n}$ für $k=7$ und $n=10$ als oranges -Rechteck dargestellt. - -\begin{figure} -\centering -\includegraphics[width=0.84\textwidth]{papers/dreieck/images/beta.pdf} -\caption{Wahrscheinlichkeitsdichte der Beta-Verteilung -$\beta(a,b,x)$ -für verschiedene Werte der Parameter $a$ und $b$. -Die Werte des Parameters für einen Graphen einer Beta-Verteilung -sind als Punkt im kleinen Quadrat rechts -im Graphen als Punkt mit der gleichen Farbe dargestellt. -\label{dreieck:fig:betaverteilungn}} -\end{figure} - -Die Formel~\eqref{dreieck:eqn:ordnungsstatistik:varianz} -besagt auch, dass die Varianz der proportional ist zu $k((n+1)-k)$. -Dieser Ausdruck ist am grössten für $k=(n+1)/2$, die Varianz ist -also grösser für die ``mittleren'' Ordnungstatistiken als für die -extremen $X_{1:n}=\operatorname{min}(X_1,\dots,X_n)$ und -$X_{n:n}=\operatorname{max}(X_1,\dots,X_n)$. -- cgit v1.2.1 From 8ecf2cd93564d76aee7dd81736d3cd5908b273cd Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Mon, 14 Mar 2022 12:41:22 +0100 Subject: move command configuration to common/Makefile.inc --- buch/Makefile | 28 ++++++++++++++-------------- buch/common/Makefile.inc | 8 +++++++- 2 files changed, 21 insertions(+), 15 deletions(-) (limited to 'buch') diff --git a/buch/Makefile b/buch/Makefile index e2ad4c0..00fcf42 100755 --- a/buch/Makefile +++ b/buch/Makefile @@ -18,15 +18,15 @@ ALLTEXFILES = $(TEXFILES) $(CHAPTERFILES) # Buchblock für Druckerei # buch.pdf: buch.tex $(TEXFILES) buch.ind $(BLXFILES) - pdflatex buch.tex - bibtex buch + $(pdflatex) buch.tex + $(bibtex) buch buch.idx: buch.tex $(TEXFILES) images - touch buch.ind - pdflatex buch.tex + $(touch) buch.ind + $(pdflatex) buch.tex buch.ind: buch.idx - makeindex buch.idx + $(makeindex) buch.idx # # Papers in einzelne PDF-Files separieren für digitales Feedback @@ -39,16 +39,16 @@ separate: buch.aux buch.pdf # SeminarSpezielleFunktionen.pdf: SeminarSpezielleFunktionen.tex $(TEXFILES) \ SeminarSpezielleFunktionen.ind $(BLXFILES) - pdflatex SeminarSpezielleFunktionen.tex - bibtex SeminarSpezielleFunktionen + $(pdflatex) SeminarSpezielleFunktionen.tex + $(bibtex) SeminarSpezielleFunktionen SeminarSpezielleFunktionen.idx: SeminarSpezielleFunktionen.tex $(TEXFILES) \ images - touch SeminarSpezielleFunktionen.ind - pdflatex SeminarSpezielleFunktionen.tex + $(touch) SeminarSpezielleFunktionen.ind + $(pdflatex) SeminarSpezielleFunktionen.tex SeminarSpezielleFunktionen.ind: SeminarSpezielleFunktionen.idx - makeindex SeminarSpezielleFunktionen + $(makeindex) SeminarSpezielleFunktionen # # This Makefile can also construct the short tests @@ -56,17 +56,17 @@ SeminarSpezielleFunktionen.ind: SeminarSpezielleFunktionen.idx tests: test1.pdf test2.pdf test3.pdf test1.pdf: common/test-common.tex common/test1.tex aufgaben1.tex - pdflatex common/test1.tex + $(pdflatex) common/test1.tex test2.pdf: common/test-common.tex common/test1.tex aufgaben2.tex - pdflatex common/test2.tex + $(pdflatex) common/test2.tex test3.pdf: common/test-common.tex common/test1.tex aufgaben3.tex - pdflatex common/test3.tex + $(pdflatex) common/test3.tex # # Errata # errata.pdf: errata.tex - pdflatex errata.tex + $(pdflatex) errata.tex diff --git a/buch/common/Makefile.inc b/buch/common/Makefile.inc index c8b0f6e..b9461e5 100755 --- a/buch/common/Makefile.inc +++ b/buch/common/Makefile.inc @@ -4,9 +4,15 @@ # (c) 2021 Prof Dr Andreas Mueller, OST Ostschweizer Fachhochschule # - SUBDIRECTORIES = chapters +# change the following variables to suit your environment + +pdflatex = pdflatex +bibtex = bibtex +makeindex = makeindex +touch = touch + .PHONY: images images: -- cgit v1.2.1 From bfaa02edc586b384691399c81aad54d92b5de986 Mon Sep 17 00:00:00 2001 From: Nao Pross Date: Tue, 15 Mar 2022 09:48:35 +0100 Subject: kugel: Add paper title and notes --- buch/papers/kugel/Makefile.inc | 8 ++---- buch/papers/kugel/main.tex | 57 ++++++++++++++++++++++-------------------- buch/papers/kugel/teil0.tex | 22 ---------------- buch/papers/kugel/teil1.tex | 55 ---------------------------------------- buch/papers/kugel/teil2.tex | 40 ----------------------------- buch/papers/kugel/teil3.tex | 40 ----------------------------- 6 files changed, 32 insertions(+), 190 deletions(-) delete mode 100644 buch/papers/kugel/teil0.tex delete mode 100644 buch/papers/kugel/teil1.tex delete mode 100644 buch/papers/kugel/teil2.tex delete mode 100644 buch/papers/kugel/teil3.tex (limited to 'buch') diff --git a/buch/papers/kugel/Makefile.inc b/buch/papers/kugel/Makefile.inc index d926229..50d6825 100644 --- a/buch/papers/kugel/Makefile.inc +++ b/buch/papers/kugel/Makefile.inc @@ -4,11 +4,7 @@ # (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule # dependencies-kugel = \ - papers/kugel/packages.tex \ + papers/kugel/packages.tex \ papers/kugel/main.tex \ - papers/kugel/references.bib \ - papers/kugel/teil0.tex \ - papers/kugel/teil1.tex \ - papers/kugel/teil2.tex \ - papers/kugel/teil3.tex + papers/kugel/references.bib diff --git a/buch/papers/kugel/main.tex b/buch/papers/kugel/main.tex index 0e632ec..06368af 100644 --- a/buch/papers/kugel/main.tex +++ b/buch/papers/kugel/main.tex @@ -1,36 +1,39 @@ % + % main.tex -- Paper zum Thema % % (c) 2020 Hochschule Rapperswil % -\chapter{Thema\label{chapter:kugel}} -\lhead{Thema} +\chapter{Recurrence Relations for Spherical Harmonics in Quantum Mechanics\label{chapter:kugel}} +\lhead{Recurrence Relations in Quantum Mechanics} \begin{refsection} -\chapterauthor{Hans Muster} - -Ein paar Hinweise für die korrekte Formatierung des Textes -\begin{itemize} -\item -Absätze werden gebildet, indem man eine Leerzeile einfügt. -Die Verwendung von \verb+\\+ ist nur in Tabellen und Arrays gestattet. -\item -Die explizite Platzierung von Bildern ist nicht erlaubt, entsprechende -Optionen werden gelöscht. -Verwenden Sie Labels und Verweise, um auf Bilder hinzuweisen. -\item -Beginnen Sie jeden Satz auf einer neuen Zeile. -Damit ermöglichen Sie dem Versionsverwaltungssysteme, Änderungen -in verschiedenen Sätzen von verschiedenen Autoren ohne Konflikt -anzuwenden. -\item -Bilden Sie auch für Formeln kurze Zeilen, einerseits der besseren -Übersicht wegen, aber auch um GIT die Arbeit zu erleichtern. -\end{itemize} - -\input{papers/kugel/teil0.tex} -\input{papers/kugel/teil1.tex} -\input{papers/kugel/teil2.tex} -\input{papers/kugel/teil3.tex} +\chapterauthor{Manuel Cattaneo, Naoki Pross} + +\begin{verbatim} + +Ideas and current research goals +-------------------------------- + +- Recurrence relations for spherical harmonics +- Associated Legendre polynomials +- Rodrigues' type formula aka Rodrigues' formula +- Applications: + * Quantization of angular momentum + * Gravitational field measurements (NASA ebb and flow, ESA goce) + * Literally anything that needs basis functions on the surface of a sphere + +Literature +---------- + +- Nichtkommutative Bildverarbeitung, T. Mendez, p57+ +- Linear Algebra Done Right, S. Axler, p212,221,231,237 +- Introduction to Quantum Mechanics, D. J. Griffith, p201+ +- Seminar Quantenmechanik, A. Müller, p101,106,114,121 +- Introduction to Partial Differential Equations, J. Oliver, p510+ +- Partial Differential Equations in Engineering Problems, K. Miller, p175,190 + +\end{verbatim} + \printbibliography[heading=subbibliography] \end{refsection} diff --git a/buch/papers/kugel/teil0.tex b/buch/papers/kugel/teil0.tex deleted file mode 100644 index f921a82..0000000 --- a/buch/papers/kugel/teil0.tex +++ /dev/null @@ -1,22 +0,0 @@ -% -% einleitung.tex -- Beispiel-File für die Einleitung -% -% (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil -% -\section{Teil 0\label{kugel:section:teil0}} -\rhead{Teil 0} -Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam -nonumy eirmod tempor invidunt ut labore et dolore magna aliquyam -erat, sed diam voluptua \cite{kugel:bibtex}. -At vero eos et accusam et justo duo dolores et ea rebum. -Stet clita kasd gubergren, no sea takimata sanctus est Lorem ipsum -dolor sit amet. - -Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam -nonumy eirmod tempor invidunt ut labore et dolore magna aliquyam -erat, sed diam voluptua. -At vero eos et accusam et justo duo dolores et ea rebum. Stet clita -kasd gubergren, no sea takimata sanctus est Lorem ipsum dolor sit -amet. - - diff --git a/buch/papers/kugel/teil1.tex b/buch/papers/kugel/teil1.tex deleted file mode 100644 index e56bb18..0000000 --- a/buch/papers/kugel/teil1.tex +++ /dev/null @@ -1,55 +0,0 @@ -% -% teil1.tex -- Beispiel-File für das Paper -% -% (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil -% -\section{Teil 1 -\label{kugel:section:teil1}} -\rhead{Problemstellung} -Sed ut perspiciatis unde omnis iste natus error sit voluptatem -accusantium doloremque laudantium, totam rem aperiam, eaque ipsa -quae ab illo inventore veritatis et quasi architecto beatae vitae -dicta sunt explicabo. -Nemo enim ipsam voluptatem quia voluptas sit aspernatur aut odit -aut fugit, sed quia consequuntur magni dolores eos qui ratione -voluptatem sequi nesciunt -\begin{equation} -\int_a^b x^2\, dx -= -\left[ \frac13 x^3 \right]_a^b -= -\frac{b^3-a^3}3. -\label{kugel:equation1} -\end{equation} -Neque porro quisquam est, qui dolorem ipsum quia dolor sit amet, -consectetur, adipisci velit, sed quia non numquam eius modi tempora -incidunt ut labore et dolore magnam aliquam quaerat voluptatem. - -Ut enim ad minima veniam, quis nostrum exercitationem ullam corporis -suscipit laboriosam, nisi ut aliquid ex ea commodi consequatur? -Quis autem vel eum iure reprehenderit qui in ea voluptate velit -esse quam nihil molestiae consequatur, vel illum qui dolorem eum -fugiat quo voluptas nulla pariatur? - -\subsection{De finibus bonorum et malorum -\label{kugel:subsection:finibus}} -At vero eos et accusamus et iusto odio dignissimos ducimus qui -blanditiis praesentium voluptatum deleniti atque corrupti quos -dolores et quas molestias excepturi sint occaecati cupiditate non -provident, similique sunt in culpa qui officia deserunt mollitia -animi, id est laborum et dolorum fuga \eqref{000tempmlate:equation1}. - -Et harum quidem rerum facilis est et expedita distinctio -\ref{kugel:section:loesung}. -Nam libero tempore, cum soluta nobis est eligendi optio cumque nihil -impedit quo minus id quod maxime placeat facere possimus, omnis -voluptas assumenda est, omnis dolor repellendus -\ref{kugel:section:folgerung}. -Temporibus autem quibusdam et aut officiis debitis aut rerum -necessitatibus saepe eveniet ut et voluptates repudiandae sint et -molestiae non recusandae. -Itaque earum rerum hic tenetur a sapiente delectus, ut aut reiciendis -voluptatibus maiores alias consequatur aut perferendis doloribus -asperiores repellat. - - diff --git a/buch/papers/kugel/teil2.tex b/buch/papers/kugel/teil2.tex deleted file mode 100644 index cb9e427..0000000 --- a/buch/papers/kugel/teil2.tex +++ /dev/null @@ -1,40 +0,0 @@ -% -% teil2.tex -- Beispiel-File für teil2 -% -% (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil -% -\section{Teil 2 -\label{kugel:section:teil2}} -\rhead{Teil 2} -Sed ut perspiciatis unde omnis iste natus error sit voluptatem -accusantium doloremque laudantium, totam rem aperiam, eaque ipsa -quae ab illo inventore veritatis et quasi architecto beatae vitae -dicta sunt explicabo. Nemo enim ipsam voluptatem quia voluptas sit -aspernatur aut odit aut fugit, sed quia consequuntur magni dolores -eos qui ratione voluptatem sequi nesciunt. Neque porro quisquam -est, qui dolorem ipsum quia dolor sit amet, consectetur, adipisci -velit, sed quia non numquam eius modi tempora incidunt ut labore -et dolore magnam aliquam quaerat voluptatem. Ut enim ad minima -veniam, quis nostrum exercitationem ullam corporis suscipit laboriosam, -nisi ut aliquid ex ea commodi consequatur? Quis autem vel eum iure -reprehenderit qui in ea voluptate velit esse quam nihil molestiae -consequatur, vel illum qui dolorem eum fugiat quo voluptas nulla -pariatur? - -\subsection{De finibus bonorum et malorum -\label{kugel:subsection:bonorum}} -At vero eos et accusamus et iusto odio dignissimos ducimus qui -blanditiis praesentium voluptatum deleniti atque corrupti quos -dolores et quas molestias excepturi sint occaecati cupiditate non -provident, similique sunt in culpa qui officia deserunt mollitia -animi, id est laborum et dolorum fuga. Et harum quidem rerum facilis -est et expedita distinctio. Nam libero tempore, cum soluta nobis -est eligendi optio cumque nihil impedit quo minus id quod maxime -placeat facere possimus, omnis voluptas assumenda est, omnis dolor -repellendus. Temporibus autem quibusdam et aut officiis debitis aut -rerum necessitatibus saepe eveniet ut et voluptates repudiandae -sint et molestiae non recusandae. Itaque earum rerum hic tenetur a -sapiente delectus, ut aut reiciendis voluptatibus maiores alias -consequatur aut perferendis doloribus asperiores repellat. - - diff --git a/buch/papers/kugel/teil3.tex b/buch/papers/kugel/teil3.tex deleted file mode 100644 index 734fff9..0000000 --- a/buch/papers/kugel/teil3.tex +++ /dev/null @@ -1,40 +0,0 @@ -% -% teil3.tex -- Beispiel-File für Teil 3 -% -% (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil -% -\section{Teil 3 -\label{kugel:section:teil3}} -\rhead{Teil 3} -Sed ut perspiciatis unde omnis iste natus error sit voluptatem -accusantium doloremque laudantium, totam rem aperiam, eaque ipsa -quae ab illo inventore veritatis et quasi architecto beatae vitae -dicta sunt explicabo. Nemo enim ipsam voluptatem quia voluptas sit -aspernatur aut odit aut fugit, sed quia consequuntur magni dolores -eos qui ratione voluptatem sequi nesciunt. Neque porro quisquam -est, qui dolorem ipsum quia dolor sit amet, consectetur, adipisci -velit, sed quia non numquam eius modi tempora incidunt ut labore -et dolore magnam aliquam quaerat voluptatem. Ut enim ad minima -veniam, quis nostrum exercitationem ullam corporis suscipit laboriosam, -nisi ut aliquid ex ea commodi consequatur? Quis autem vel eum iure -reprehenderit qui in ea voluptate velit esse quam nihil molestiae -consequatur, vel illum qui dolorem eum fugiat quo voluptas nulla -pariatur? - -\subsection{De finibus bonorum et malorum -\label{kugel:subsection:malorum}} -At vero eos et accusamus et iusto odio dignissimos ducimus qui -blanditiis praesentium voluptatum deleniti atque corrupti quos -dolores et quas molestias excepturi sint occaecati cupiditate non -provident, similique sunt in culpa qui officia deserunt mollitia -animi, id est laborum et dolorum fuga. Et harum quidem rerum facilis -est et expedita distinctio. Nam libero tempore, cum soluta nobis -est eligendi optio cumque nihil impedit quo minus id quod maxime -placeat facere possimus, omnis voluptas assumenda est, omnis dolor -repellendus. Temporibus autem quibusdam et aut officiis debitis aut -rerum necessitatibus saepe eveniet ut et voluptates repudiandae -sint et molestiae non recusandae. Itaque earum rerum hic tenetur a -sapiente delectus, ut aut reiciendis voluptatibus maiores alias -consequatur aut perferendis doloribus asperiores repellat. - - -- cgit v1.2.1 From 9f542be91a61e6eb1fcb3c29e8a908eb9e5b8949 Mon Sep 17 00:00:00 2001 From: Marc Benz Date: Tue, 15 Mar 2022 14:00:47 +0100 Subject: initial commit from Marc Benz --- buch/papers/transfer/main.tex | 23 ++--------------------- 1 file changed, 2 insertions(+), 21 deletions(-) (limited to 'buch') diff --git a/buch/papers/transfer/main.tex b/buch/papers/transfer/main.tex index 2aae635..ed16998 100644 --- a/buch/papers/transfer/main.tex +++ b/buch/papers/transfer/main.tex @@ -3,29 +3,10 @@ % % (c) 2020 Hochschule Rapperswil % -\chapter{Thema\label{chapter:transfer}} +\chapter{Transferfunktionen\label{chapter:transfer}} \lhead{Thema} \begin{refsection} -\chapterauthor{Hans Muster} - -Ein paar Hinweise für die korrekte Formatierung des Textes -\begin{itemize} -\item -Absätze werden gebildet, indem man eine Leerzeile einfügt. -Die Verwendung von \verb+\\+ ist nur in Tabellen und Arrays gestattet. -\item -Die explizite Platzierung von Bildern ist nicht erlaubt, entsprechende -Optionen werden gelöscht. -Verwenden Sie Labels und Verweise, um auf Bilder hinzuweisen. -\item -Beginnen Sie jeden Satz auf einer neuen Zeile. -Damit ermöglichen Sie dem Versionsverwaltungssysteme, Änderungen -in verschiedenen Sätzen von verschiedenen Autoren ohne Konflikt -anzuwenden. -\item -Bilden Sie auch für Formeln kurze Zeilen, einerseits der besseren -Übersicht wegen, aber auch um GIT die Arbeit zu erleichtern. -\end{itemize} +\chapterauthor{Marc Benz} \input{papers/transfer/teil0.tex} \input{papers/transfer/teil1.tex} -- cgit v1.2.1 From 3381bfe1e6ebcd66249cd4c6f49bdd820643a5be Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Tue, 15 Mar 2022 17:35:00 +0100 Subject: add some stuff about separation --- buch/chapters/090-pde/kreis.tex | 2 +- buch/chapters/090-pde/kugel.tex | 146 ++++++++++++++++++++++++++++++++++++++ buch/papers/000template/main.tex | 5 +- buch/papers/000template/teil0.tex | 3 + buch/papers/000template/teil1.tex | 3 + buch/papers/000template/teil2.tex | 3 + buch/papers/000template/teil3.tex | 3 + 7 files changed, 163 insertions(+), 2 deletions(-) (limited to 'buch') diff --git a/buch/chapters/090-pde/kreis.tex b/buch/chapters/090-pde/kreis.tex index b4ce8d7..c60fd44 100644 --- a/buch/chapters/090-pde/kreis.tex +++ b/buch/chapters/090-pde/kreis.tex @@ -32,7 +32,7 @@ Der Laplace-Operator hat in Polarkoordinaten die Form \frac1r \frac{\partial}{\partial r} + -\frac{1}{r 2} +\frac{1}{r^2} \frac{\partial^2}{\partial\varphi^2}. \label{buch:pde:kreis:laplace} \end{equation} diff --git a/buch/chapters/090-pde/kugel.tex b/buch/chapters/090-pde/kugel.tex index 0e3524f..c081029 100644 --- a/buch/chapters/090-pde/kugel.tex +++ b/buch/chapters/090-pde/kugel.tex @@ -5,4 +5,150 @@ % \section{Kugelfunktionen \label{buch:pde:section:kugel}} +Kugelsymmetrische Probleme können oft vorteilhaft in Kugelkoordinaten +beschrieben werden. +Die Separationsmethode kann auf partielle Differentialgleichungen +mit dem Laplace-Operator angewendet werden. +Die daraus resultierenden gewöhnlichen Differentialgleichungen führen +einerseits auf die Laguerre-Differentialgleichung für den radialen +Anteil sowie auf Kugelfunktionen für die Koordinaten der +geographischen Länge und Breite. + +\subsection{Kugelkoordinaten} +Wir verwenden Kugelkoordinaten $(r,\vartheta,\varphi)$, wobei $r$ +der Radius ist, $\vartheta$ die geographische Breite gemessen vom +Nordpol der Kugel und $\varphi$ die geographische Breite. +Der Definitionsbereich für Kugelkoordinaten ist +\[ +\Omega += +\{(r,\vartheta,\varphi) +\;|\; +r\ge 0\wedge +0\le \vartheta\le \pi\wedge +0\le \varphi< 2\pi +\}. +\] +Die Entfernung eines Punktes von der $z$-Achse ist $r\sin\vartheta$. +Daraus lassen sich die karteischen Koordinaten eines Punktes mit Hilfe +von +\[ +\begin{pmatrix}x\\y\\z\end{pmatrix} += +\begin{pmatrix} +r\cos\vartheta\\ +r\sin\vartheta\cos\varphi\\ +r\sin\vartheta\sin\varphi +\end{pmatrix}. +\] +Man beachte, dass die Punkte auf der $z$-Achse keine eindeutigen +Kugelkoordinaten haben. +Sie sind charakterisiert durch $r\sin\vartheta=0$, was $\cos\vartheta=\pm1$ +impliziert. +Entsprechend führen alle Werte von $\varphi$ auf den gleichen Punkt +$(0,0,\pm r)$. + +\subsection{Der Laplace-Operator in Kugelkoordinaten} +Der Laplace-Operator in Kugelkoordinaten lautet +\begin{align} +\Delta +&= +\frac{1}{r^2} \frac{\partial}{\partial r}r^2\frac{\partial}{\partial r} ++ +\frac{1}{r^2\sin\vartheta}\frac{\partial}{\partial\vartheta} +\sin\vartheta\frac{\partial}{\partial\vartheta} ++ +\frac{1}{r^2\sin^2\vartheta}\frac{\partial^2}{\partial\varphi^2}. +\label{buch:pde:kugel:laplace1} +\intertext{Dies kann auch geschrieben werden als} +&= +\frac{\partial^2}{\partial r^2} ++ +\frac{2}{r}\frac{\partial}{\partial r} ++ +\frac{1}{r^2\sin\vartheta}\frac{\partial}{\partial\vartheta} +\sin\vartheta\frac{\partial}{\partial\vartheta} ++ +\frac{1}{r^2\sin^2\vartheta}\frac{\partial^2}{\partial\varphi^2} +\label{buch:pde:kugel:laplace2} +\intertext{oder} +&= +\frac{1}{r} +\frac{\partial^2}{\partial r^2} r ++ +\frac{1}{r^2\sin\vartheta}\frac{\partial}{\partial\vartheta} +\sin\vartheta\frac{\partial}{\partial\vartheta} ++ +\frac{1}{r^2\sin^2\vartheta}\frac{\partial^2}{\partial\varphi^2}. +\label{buch:pde:kugel:laplace3} +\end{align} +Dabei ist zu berücksichtigen, dass mit der Notation gemeint ist, +dass ein Ableitungsoperator auf alles wirkt, was rechts im gleichen +Term steht. +Der Operator +\[ +\frac{1}{r} +\frac{\partial^2}{\partial r^2}r +\quad\text{wirkt daher als}\quad +\frac{1}{r} +\frac{\partial^2}{\partial r^2}rf += +\frac{1}{r} +\frac{\partial}{\partial r}\biggl(f + r\frac{\partial f}{\partial r}\biggr) += +\frac{1}{r} +\frac{\partial f}{\partial r} ++ +\frac{1}{r} +\frac{\partial f}{\partial r} ++ +\frac{\partial^2f}{\partial r^2}. += +\frac{2}{r}\frac{\partial f}{\partial r} ++ +\frac{\partial^2f}{\partial r^2}, +\] +was die Äquivalenz der beiden Formen +\eqref{buch:pde:kugel:laplace2} +und +\eqref{buch:pde:kugel:laplace3} +rechtfertigt. +Auch die Äquivalenz mit +\eqref{buch:pde:kugel:laplace1} +kann auf ähnliche Weise verstanden werden. + +Die Herleitung dieser Formel ist ziemlich aufwendig und soll hier +nicht dargestellt werden. +Es sei aber darauf hingewiesen, dass sich für $\vartheta=\frac{\pi}2$ +wegen $\sin\vartheta=\sin\frac{\pi}2=1$ +der eingeschränkte Operator +\[ +\Delta += +\frac{1}{r^2}\frac{\partial }{\partial r} r^2\frac{\partial}{\partial r} ++ +\frac{1}{r^2}\frac{\partial^2}{\partial\varphi^2} +\] +ergibt. +Wendet man wie oben die Produktregel auf den ersten Term an, entsteht die +Form +\[ +\frac{\partial^2}{\partial r^2} ++ +\frac{2}{r} +\frac{\partial}{\partial r} ++ +\frac{1}{r^2}\frac{\partial^2}{\partial\varphi^2} +\] +die {\em nicht} übereinstimmt mit dem Laplace-Operator in +Polarkoordinaten~\eqref{buch:pde:kreis:laplace}. +Der Unterschied rührt daher, dass der Laplace-Operator die Krümmung +der Koordinatenlinien berücksichtigt, in diesem Fall der Meridiane. + + +\subsection{Separation} + + + + diff --git a/buch/papers/000template/main.tex b/buch/papers/000template/main.tex index 87a5685..91b6d6e 100644 --- a/buch/papers/000template/main.tex +++ b/buch/papers/000template/main.tex @@ -1,7 +1,10 @@ % % main.tex -- Paper zum Thema <000template> % -% (c) 2020 Hochschule Rapperswil +% (c) 2020 Autor, OST Ostschweizer Fachhochschule +% +% !TEX root = ../../buch.tex +% !TEX encoding = UTF-8 % \chapter{Thema\label{chapter:000template}} \lhead{Thema} diff --git a/buch/papers/000template/teil0.tex b/buch/papers/000template/teil0.tex index 7b9f088..65d7ae1 100644 --- a/buch/papers/000template/teil0.tex +++ b/buch/papers/000template/teil0.tex @@ -3,6 +3,9 @@ % % (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil % +% !TEX root = ../../buch.tex +% !TEX encoding = UTF-8 +% \section{Teil 0\label{000template:section:teil0}} \rhead{Teil 0} Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam diff --git a/buch/papers/000template/teil1.tex b/buch/papers/000template/teil1.tex index 00d3058..0f8dfae 100644 --- a/buch/papers/000template/teil1.tex +++ b/buch/papers/000template/teil1.tex @@ -3,6 +3,9 @@ % % (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil % +% !TEX root = ../../buch.tex +% !TEX encoding = UTF-8 +% \section{Teil 1 \label{000template:section:teil1}} \rhead{Problemstellung} diff --git a/buch/papers/000template/teil2.tex b/buch/papers/000template/teil2.tex index 471adae..496557f 100644 --- a/buch/papers/000template/teil2.tex +++ b/buch/papers/000template/teil2.tex @@ -3,6 +3,9 @@ % % (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil % +% !TEX root = ../../buch.tex +% !TEX encoding = UTF-8 +% \section{Teil 2 \label{000template:section:teil2}} \rhead{Teil 2} diff --git a/buch/papers/000template/teil3.tex b/buch/papers/000template/teil3.tex index 4697813..ef2aa75 100644 --- a/buch/papers/000template/teil3.tex +++ b/buch/papers/000template/teil3.tex @@ -3,6 +3,9 @@ % % (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil % +% !TEX root = ../../buch.tex +% !TEX encoding = UTF-8 +% \section{Teil 3 \label{000template:section:teil3}} \rhead{Teil 3} -- cgit v1.2.1 From 61202dbd7a7762ceeee673cf27da26e47d72b966 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Wed, 16 Mar 2022 09:42:00 +0100 Subject: Kugelfunktionen --- buch/chapters/090-pde/kugel.tex | 237 +++++++++++++++++++++++++++++++++++++++- 1 file changed, 236 insertions(+), 1 deletion(-) (limited to 'buch') diff --git a/buch/chapters/090-pde/kugel.tex b/buch/chapters/090-pde/kugel.tex index c081029..d466e26 100644 --- a/buch/chapters/090-pde/kugel.tex +++ b/buch/chapters/090-pde/kugel.tex @@ -145,8 +145,243 @@ Polarkoordinaten~\eqref{buch:pde:kreis:laplace}. Der Unterschied rührt daher, dass der Laplace-Operator die Krümmung der Koordinatenlinien berücksichtigt, in diesem Fall der Meridiane. - \subsection{Separation} +In Abschnitt~\ref{buch:pde:subsection:eigenwertproblem} +wurde bereits gzeigt, wie die Wellengleichung +\[ +\frac{1}{c^2} +\frac{\partial^2 U}{\partial t^2} +-\Delta U += +0 +\] +durch Separation der Zeit auf ein Eigenwertproblem für eine +Funktion $u$ reduziert werden kann, die nur von den Ortskoordinaten +abhängt. +Es geht also nur noch darum, dass Eigenwertproblem +\[ +\Delta u = -\lambda^2 u +\] +mit geeigneten Randbedingungen zu lösen. +Dazu gehören einerseits eventuelle Gebietsränder, die im Moment +nicht interessieren. +Andererseits muss sichergestellt sein, dass die Lösungsfunktionen +stetig und differentierbar sind an den Orten, wo das Koordinatensystem +singulär ist. +So müssen $u(r,\vartheta,\varphi)$ $2\pi$-periodisch in $\varphi$ sein. +% XXX Ableitungen + +\subsubsection{Separation des radialen Anteils} +Für das Eigenwertproblem verwenden wir den Ansatz +\[ +u(r,\vartheta,\varphi) += +R(r) \Theta(\vartheta) \Phi(\varphi), +\] +den wir in die Differentialgleichung einsetzen. +So erhalten wir +\[ +\biggl(\frac{1}{r^2}R''(r)+\frac{2}{r}R'(r) \biggr) +\Theta(\vartheta)\Phi(\varphi) ++ +R(r) +\frac{1}{r^2\sin\vartheta} +\frac{\partial}{\partial\vartheta}(\sin\vartheta \Theta'(\vartheta)) +\Phi(\varphi) ++ +R(r)\Theta(\vartheta) +\frac{1}{r^2\sin\vartheta} \Phi''(\varphi) += +-\lambda^2 R(r)\Theta(\vartheta)\Phi(\varphi). +\] +Die Gleichung lässt sich nach Multiplikation mit $r^2$ und +Division durch $u$ separieren in +\begin{equation} +\frac{R''(r)+2rR'(r)+\lambda^2r^2}{R(r)} ++ +\frac{1}{\Theta(\vartheta) \sin\vartheta} +\frac{\partial}{\partial\vartheta}\sin\vartheta\Theta'(\vartheta) ++ +\frac{1}{\sin^2\vartheta}\frac{\Phi''(\varphi)}{\Phi(\varphi)} += +0 +\label{buch:pde:kugel:separiert2} +\end{equation} +Der erste Term hängt nur von $r$ ab, die anderen nur von $\vartheta$ und +$\varphi$, daher muss der erste Term konstant sein. +Damit ergbit sich für den Radialanteil die gewöhnliche Differentialgleichung +\[ +R''(r) + 2rR'(r) +\lambda^2 r^2 = \mu^2 R(r), +\] +die zum Beispiel mit der Potenzreihenmethode gelöst werden kann. +Sie kann aber durch eine geeignete Substition nochmals auf die +Laguerre-Differentialgleichung reduziert werden, wie in +Kapitel~\ref{chapter:laguerre} dargelegt wird. + +\subsubsection{Kugelflächenanteil} +Für die Separation der verbleibenden winkelabhängigen Teile muss die +Gleichung +\[ +\frac{1}{\Theta(\vartheta) \sin\vartheta} +\frac{\partial}{\partial\vartheta}\sin\vartheta\Theta'(\vartheta) ++ +\frac{1}{\sin^2\vartheta}\frac{\Phi''(\varphi)}{\Phi(\varphi)} += +-\mu^2 +\] +mit $\sin^2\vartheta$ multipliziert werden, was auf +\[ +\frac{\sin\vartheta}{\Theta(\vartheta)} +\frac{\partial}{\partial\vartheta}\sin\vartheta\Theta'(\vartheta) ++ +\frac{\Phi''(\varphi)}{\Phi(\varphi)} += +-\mu^2\sin^2\vartheta +\quad\Rightarrow\quad +\frac{\sin\vartheta}{\Theta(\vartheta)} +\frac{\partial}{\partial\vartheta}\sin\vartheta\Theta'(\vartheta) ++ +\mu^2\sin^2\vartheta += +- +\frac{\Phi''(\varphi)}{\Phi(\varphi)} +\] +führt. +Die linke Seite der letzten Gleichung hängt nur von $\vartheta$ +ab, die rechte nur von $\varphi$, beide Seiten müssen daher +konstant sein, wir bezeichnen diese Konstante mit $\alpha^2$. +So ergibt sich die Differentialgleichung +\[ +\alpha^2 += +-\frac{\Phi''(\varphi)}{\Phi(\varphi)} +\] +für die Abhängigkeit von $\varphi$, mit der allgemeinen Lösung +\[ +\Phi(\varphi) += +A\cos\alpha \varphi ++ +B\sin\alpha \varphi. +\] +Die Randbedingungen verlangen, dass $\Phi(\varphi)$ eine $2\pi$-periodische +Funktion ist, was genau dann möglich ist, wenn $\alpha=m$ ganzzahlig ist. +Damit ergibt sich für die $\vartheta$-Abhängigkeit die Differentialgleichung +\begin{equation} +\frac{\sin\vartheta}{\Theta(\vartheta)} +\frac{\partial}{\partial\vartheta}\sin\vartheta\Theta'(\vartheta) ++ +\mu^2\sin^2\vartheta += +m^2. +\label{buch:pde:kugel:eqn:thetaanteil} +\end{equation} + +\subsubsection{Abhängigkeit von $\vartheta$} +Die Differentialgleichung~\eqref{buch:pde:kugel:eqn:thetaanteil} +ist etwas unhandlich, daher verwenden wir die Substitution $z=\cos\vartheta$, +um die trigonometrischen Funktionen los zu werden. +Wegen +\[ +\frac{dz}{d\vartheta} = -\sin\vartheta =-\sqrt{1-z^2} +\] +können die Ableitungen nach $\vartheta$ auch durch Ableitungen nach $z$ +ausgedrückt werden. +Wir schreiben dazu $Z(z)=\Theta(\vartheta)$ und berechnen +\[ +\Theta'(\vartheta) += +\frac{d\Theta}{d\vartheta} += +\frac{dZ}{dz}\frac{dz}{d\vartheta} += +- +\sqrt{1-z^2} +Z'(z). +\] +Dies bedeutet auch, dass +\[ +\sin\vartheta\frac{d}{d\vartheta} += +- +(1-z^2)\frac{d}{dz}, +\] +damit lässt sich die Differentialgleichung für $\Theta(\vartheta)$ umschreiben +in eine Differentialgleichung für $Z(z)$, nämlich +\[ +(1-z^2)\frac{d}{dz}(1-z^2)\frac{d}{dz} Z(z) ++ +\mu^2 +(1-z^2) +Z(z) += +m^2 +Z(z). +\] +Indem man die Ableitung im ersten Term mit Hilfe der Produktregel +ausführt, kann man die Gleichung +\[ +(1-z^2)\biggl( +-2zZ'(z) + (1-z^2)Z''(z) +\biggr) ++ +\mu^2(1-z^2)Z(z) += +-m^2 Z(z) +\] +bekommen. +Division durch $1-z^2$ ergibt die +{\em Legendre-Differentialgleichung} +\begin{equation} +(1-z^2)Z''(z) +-2zZ'(z) ++ +\biggl( +\mu^2 - \frac{m^2}{1-z^2} +\biggr) +Z(z) += +0. +\label{buch:pde:kugel:eqn:legendre-dgl} +\end{equation} +Eine Diskussion der Lösungen dieser Differentialgleichung erfolgt im +Kapitel~\ref{chapter:kugel}. + +\subsection{Kugelfunktionen} +Die Legendre-Differentialgleichung~\eqref{buch:pde:kugel:eqn:legendre-dgl} +hat Lösungen für Werte von $\mu$ derart, dass $\mu^2=l(l+1)$ für natürliche +Zahlen $l$. +Die Lösungen sind sogar Polynome, die wir mit $P_l^{(m)}(z)$ +bezeichnen, dabei ist $m$ eine ganze Zahl mit $-l\le m\le l$. +Die Funktionen $P_l^{(m)}(\cos\vartheta)e^{im\varphi}$ +sind daher alle Lösungen des von $\vartheta$ und $\varphi$ +abhängigen Teils der Lösungen des Eigenwertproblems. +Mit einer geeigneten Normierung kann man zudem eine Familie von +bezüglich des Skalarproduktes +\[ +\langle f,g\rangle_{S^2} += +\int_{-\pi}^{\pi} +\int_{0}^{\pi} +\overline{f(\vartheta,\varphi)} +g(\vartheta,\varphi) +\sin\vartheta +\,d\vartheta +\,d\varphi +\] +orthonormiete Funktionen auf der Kugeloberfläche erhalten, die +man normalerweise als +\[ +Y_{lm}(\vartheta,\varphi) += +\frac{1}{\sqrt{2\pi}} +\sqrt{ +\frac{2l+1}{2}\cdot +\frac{(l-m)!}{(l+m)!} +} +P_{l}^{(m)}(\cos\vartheta)e^{im\varphi} +\] +bezeichnet. -- cgit v1.2.1 From fd8e4b595f3f2c4245e7ba201a727585e34cfc82 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Wed, 16 Mar 2022 14:36:37 +0100 Subject: add missing file --- buch/chapters/040-rekursion/betaverteilung.tex | 487 +++++++++++++++++++++++++ buch/chapters/090-pde/gleichung.tex | 1 + buch/chapters/090-pde/kreis.tex | 1 + buch/chapters/090-pde/kugel.tex | 1 + buch/chapters/090-pde/rechteck.tex | 1 + buch/chapters/090-pde/separation.tex | 1 + 6 files changed, 492 insertions(+) create mode 100644 buch/chapters/040-rekursion/betaverteilung.tex (limited to 'buch') diff --git a/buch/chapters/040-rekursion/betaverteilung.tex b/buch/chapters/040-rekursion/betaverteilung.tex new file mode 100644 index 0000000..979d04c --- /dev/null +++ b/buch/chapters/040-rekursion/betaverteilung.tex @@ -0,0 +1,487 @@ +% +% teil1.tex -- Beispiel-File für das Paper +% +% (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil +% +\subsection{Ordnungsstatistik und Beta-Funktion +\label{buch:rekursion:ordnung:section:ordnungsstatistik}} +\rhead{Ordnungsstatistik und Beta-Funktion} +In diesem Abschnitt ist $X$ eine Zufallsvariable mit der Verteilungsfunktion +$F_X(x)$, und $X_i$, $1\le i\le n$ sei ein Stichprobe von unabhängigen +Zufallsvariablen, die wie $X$ verteilt sind. +Ziel ist, die Verteilungsfunktion und die Wahrscheinlichkeitsdichte +des grössten, zweitgrössten, $k$-t-grössten Wertes in der Stichprobe +zu finden. +Wir schreiben $[n]=\{1,\dots,n\}$ für die Menge der natürlichen +Zahlen von zwischen $1$ und $n$. + +\subsubsection{Verteilung von $\operatorname{max}(X_1,\dots,X_n)$ und +$\operatorname{min}(X_1,\dots,X_n)$ +\label{buch:rekursion:ordnung:subsection:minmax}} +Die Verteilungsfunktion von $\operatorname{max}(X_1,\dots,X_n)$ hat +den Wert +\begin{align*} +F_{\operatorname{max}(X_1,\dots,X_n)}(x) +&= +P(\operatorname{max}(X_1,\dots,X_n) \le x) +\\ +&= +P(X_1\le x\wedge \dots \wedge X_n\le x) +\\ +&= +P(X_1\le x) \cdot \ldots \cdot P(X_n\le x) +\\ +&= +P(X\le x)^n += +F_X(x)^n. +\end{align*} +Für die Gleichverteilung ist +\[ +F_{\text{equi}}(x) += +\begin{cases} +0&\qquad x< 0 +\\ +x&\qquad 0\le x\le 1 +\\ +1&\qquad 1 X_1\wedge \dots \wedge x > X_n) +\\ +&= +1- +(1-P(x\le X_1)) \cdot\ldots\cdot (1-P(x\le X_n)) +\\ +&= +1-(1-F_X(x))^n, +\end{align*} +Im Speziellen für im Intervall $[0,1]$ gleichverteilte $X_i$ ist die +Verteilungsfunktion des Minimums +\[ +F_{\operatorname{min}(X_1,\dots,X_n)}(x) += +\begin{cases} +0 &\qquad x<0 \\ +1-(1-x)^n&\qquad 0\le x\le 1\\ +1 &\qquad 1 < x +\end{cases} +\] +mit Wahrscheinlichkeitsdichte +\[ +\varphi_{\operatorname{min}(X_1,\dots,X_n)} += +\frac{d}{dx} +F_{\operatorname{min}(X_1,\dots,X_n)} += +\begin{cases} +n(1-x)^{n-1}&\qquad 0\le x\le 1\\ +0 &\qquad \text{sonst} +\end{cases} +\] +und Erwartungswert +\begin{align*} +E(\operatorname{min}(X_1,\dots,X_n) +&= +\int_{-\infty}^\infty x\varphi_{\operatorname{min}(X_1,\dots,X_n)}(x)\,dx += +\int_0^1 x\cdot n(1-x)^{n-1}\,dx +\\ +&= +\bigl[ -x(1-x)^n \bigr]_0^1 + \int_0^1 (1-x)^n\,dx += +\biggl[ +- +\frac{1}{n+1} +(1-x)^{n+1} +\biggr]_0^1 += +\frac{1}{n+1}. +\end{align*} +Es ergibt sich daraus als natürlich Verallgemeinerung die Frage nach +der Verteilung des zweitegrössten oder zweitkleinsten Wertes unter den +Werten $X_i$. + +\subsubsection{Der $k$-t-grösste Wert} +Sie wieder $X_i$ eine Stichprobe von $n$ unabhängigen wie $X$ verteilten +Zufallsvariablen. +Diese werden jetzt der Grösse nach sortiert, die sortierten Werte werden +mit +\[ +X_{1:n} \le X_{2:n} \le \dots \le X_{(n-1):n} \le X_{n:n} +\] +bezeichnet. +Die Grössen $X_{k:n}$ sind Zufallsvariablen, sie heissen die $k$-ten +Ordnungsstatistiken. +Die in Abschnitt~\ref{buch:rekursion:ordnung:subsection:minmax} behandelten Zufallsvariablen +$\operatorname{min}(X_1,\dots,X_n)$ +und +$\operatorname{max}(X_1,\dots,X_n)$ +sind die Fälle +\begin{align*} +X_{1:n} &= \operatorname{min}(X_1,\dots,X_n) \\ +X_{n:n} &= \operatorname{max}(X_1,\dots,X_n). +\end{align*} + +Um den Wert der Verteilungsfunktion von $X_{k:n}$ zu berechnen, müssen wir +die Wahrscheinlichkeit bestimmen, dass $k$ der $n$ Werte $X_i$ $x$ nicht +übersteigen. +Der $k$-te Wert $X_{k:n}$ übersteigt genau dann $x$ nicht, wenn +mindestens $k$ der Zufallswerte $X_i$ $x$ nicht übersteigen, also +\[ +P(X_{k:n} \le x) += +P\left( +|\{i\in[n]\,|\, X_i\le x\}| \ge k +\right). +\] + +Das Ereignis $\{X_i\le x\}$ ist eine Bernoulli-Experiment, welches mit +Wahrscheinlichkeit $F_X(x)$ eintritt. +Die Anzahl der Zufallsvariablen $X_i$, die $x$ übertreffen, ist also +Binomialverteilt mit $p=F_X(x)$. +Damit haben wir gefunden, dass mit Wahrscheinlichkeit +\begin{equation} +F_{X_{k:n}}(x) += +P(X_{k:n}\le x) += +\sum_{i=k}^n \binom{n}{i}F_X(x)^i (1-F_X(x))^{n-i} +\label{buch:rekursion:ordnung:eqn:FXkn} +\end{equation} +mindestens $k$ der Zufallsvariablen den Wert $x$ überschreiten. + +\subsubsection{Wahrscheinlichkeitsdichte der Ordnungsstatistik} +Die Wahrscheinlichkeitsdichte der Ordnungsstatistik kann durch Ableitung +von \eqref{buch:rekursion:ordnung:eqn:FXkn} gefunden, werden, sie ist +\begin{align*} +\varphi_{X_{k:n}}(x) +&= +\frac{d}{dx} +F_{X_{k:n}}(x) +\\ +&= +\sum_{i=k}^n +\binom{n}{i} +\bigl( +iF_X(x)^{i-1}\varphi_X(x) (1-F_X(x))^{n-i} +- +F_X(x)^k +(n-i) +(1-F_X(x))^{n-i-1} +\varphi_X(x) +\bigr) +\\ +&= +\sum_{i=k}^n +\binom{n}{i} +\varphi_X(x) +F_X(x)^{i-1}(1-F_X(x))^{n-i-1} +\bigl( +iF_X(x)-(n-i)(1-F_X(x)) +\bigr) +\\ +&= +\varphi_X(x) +\biggl( +\sum_{i=k}^n i\binom{n}{i} F_X(x)^{i-1}(1-F_X(x))^{n-i} +- +\sum_{j=k}^n (n-j)\binom{n}{j} F_X(x)^{j}(1-F_X(x))^{n-j-1} +\biggr) +\\ +&= +\varphi_X(x) +\biggl( +\sum_{i=k}^n i\binom{n}{i} F_X(x)^{i-1}(1-F_X(x))^{n-i} +- +\sum_{i=k+1}^{n+1} (n-i+1)\binom{n}{i-1} F_X(x)^{i-1}(1-F_X(x))^{n-i} +\biggr) +\\ +&= +\varphi_X(x) +\biggl( +k\binom{n}{k}F_X(x)^{k-1}(1-F_X(x))^{n-k} ++ +\sum_{i=k+1}^{n+1} +\left( +i\binom{n}{i} +- +(n-i+1)\binom{n}{i-1} +\right) +F_X(x)^{i-1}(1-F_X(x))^{n-i} +\biggr) +\end{align*} +Mit den wohlbekannten Identitäten für die Binomialkoeffizienten +\begin{align*} +i\binom{n}{i} +- +(n-i+1)\binom{n}{i-1} +&= +n\binom{n-1}{i-1} +- +n +\binom{n-1}{i-1} += +0 +\end{align*} +folgt jetzt +\begin{align*} +\varphi_{X_{k:n}}(x) +&= +\varphi_X(x)k\binom{n}{k} F_X(x)^{k-1}(1-F_X(x))^{n-k}(x). +\intertext{Im Speziellen für gleichverteilte Zufallsvariablen $X_i$ ist +} +\varphi_{X_{k:n}}(x) +&= +k\binom{n}{k} x^{k-1}(1-x)^{n-k}. +\end{align*} +Dies ist die Wahrscheinlichkeitsdichte einer Betaverteilung +\[ +\beta(k,n-k+1)(x) += +\frac{1}{B(k,n-k+1)} +x^{k-1}(1-x)^{n-k}. +\] +Tatsächlich ist die Normierungskonstante +\begin{align} +\frac{1}{B(k,n-k+1)} +&= +\frac{\Gamma(n+1)}{\Gamma(k)\Gamma(n-k+1)} += +\frac{n!}{(k-1)!(n-k)!}. +\label{buch:rekursion:ordnung:betaverteilung:normierung1} +\end{align} +Andererseits ist +\[ +k\binom{n}{k} += +k\frac{n!}{k!(n-k)!} += +\frac{n!}{(k-1)!(n-k)!}, +\] +in Übereinstimmung mit~\eqref{buch:rekursion:ordnung:betaverteilung:normierung1}. +Die Verteilungsfunktion und die Wahrscheinlichkeitsdichte der +Ordnungsstatistik sind in Abbildung~\ref{buch:rekursion:ordnung:fig:order} dargestellt. + +\begin{figure} +\centering +\includegraphics{chapters/040-rekursion/images/order.pdf} +\caption{Verteilungsfunktion und Wahrscheinlichkeitsdichte der +Ordnungsstatistiken $X_{k:n}$ einer gleichverteilung Zuvallsvariable +mit $n=10$. +\label{buch:rekursion:ordnung:fig:order}} +\end{figure} + +% +% Die Beta-Funktion +% +\subsection{Die Beta-Verteilung +\label{buch:rekursion:subsection:beta-verteilung}} +Die Wahrscheinlichkeitsdichte, die im +Abschnitt~\ref{buch:rekursion:ordnung:section:ordnungsstatistik} +gefunden worden ist, ist nicht nur für ganzzahlige Exponenten +definiert. + +\begin{figure} +\centering +\includegraphics[width=0.92\textwidth]{chapters/040-rekursion/images/beta.pdf} +\caption{Wahrscheinlichkeitsdichte der Beta-Verteilung +$\beta(a,b,x)$ +für verschiedene Werte der Parameter $a$ und $b$. +Die Werte des Parameters für einen Graphen einer Beta-Verteilung +sind im kleinen Quadrat rechts im Graphen +als Punkt mit der gleichen Farbe dargestellt. +\label{buch:rekursion:ordnung:fig:betaverteilungn}} +\end{figure} + +\begin{definition} +Die Beta-Verteilung ist die Verteilung mit der Wahrscheinlichkeitsdichte +\[ +\beta_{a,b}(x) += +\begin{cases} +\displaystyle +\frac{1}{B(a,b)} +x^{a-1}(1-x)^{b-1}&\qquad 0\le x \le 1\\ +0&\qquad\text{sonst.} +\end{cases} +\] +\end{definition} + +Die Beta-Funktion ist also die Normierungskonstante der Beta-Verteilung. +Die wichtigsten Kennzahlen der Beta-Verteilung wie Erwartungswert und +Varianz lassen sich alle ebenfalls als Werte der Beta-Funktion ausdrücken. + +\subsubsection{Erwartungswert} +Mit der Wahrscheinlichkeitsdichte kann man jetzt auch den Erwartungswerte +der $k$-ten Ordnungsstatistik bestimmen. +Die Rechnung ergibt: +\begin{align*} +E(X_{k:n}) +&= +\int_0^1 x\cdot k\binom{n}{k} x^{k-1}(1-x)^{n-k}\,dx += +k +\binom{n}{k} +\int_0^1 +x^{k}(1-x)^{n-k}\,dx. +\intertext{Dies ist das Beta-Integral} +&= +k\binom{n}{k} +B(k+1,n-k+1) +\intertext{welches man durch Gamma-Funktionen bzw.~durch Fakultäten wie in} +&= +k\frac{n!}{k!(n-k)!} +\frac{\Gamma(k+1)\Gamma(n-k+1)}{n+2} += +k\frac{n!}{k!(n-k)!} +\frac{k!(n-k)!}{(n+1)!} += +\frac{k}{n+1} +\end{align*} +ausdrücken kann. +Die Erwartungswerte haben also regelmässige Abstände, sie sind in +Abbildung~\ref{buch:rekursion:ordnung:fig:order} als blaue vertikale Linien eingezeichnet. + +Für die Beta-Verteilung lässt sich die Rechnung noch allgemeiner +durchführen. +Der Erwartungswert einer $\beta_{a,b}$-verteilten Zufallsvariablen $X$ +ist +\begin{align*} +E(X) +&= +\int_0^1 x \beta_{a,b}(x)\,dx += +\frac{1}{B(a,b)} +\int_0^1 x\cdot x^{a-1}(1-x)^{b-1}\,dx += +\frac{B(a+1,b)}{B(a,b)} += +\frac{a}{a+b}. +\end{align*} +Durch Einsetzen von $a=k+1$ und $b=n-k+1$ lassen sich die für die +Ordnungsstatistik berechneten Werte wiederfinden. + +\subsubsection{Varianz} +Auch die Varianz lässt sich einfach berechnen, dazu muss zunächst +der Erwartungswert von $X_{k:n}^2$ bestimmt werden. +Er ist +\begin{align*} +E(X_{k:n}^2) +&= +\int_0^1 x^2\cdot k\binom{n}{k} x^{k-1}(1-x)^{n-k}\,dx += +k +\binom{n}{k} +\int_0^1 +x^{k+1}(1-x)^{n-k}\,dx. +\intertext{Auch dies ist ein Beta-Integral, nämlich} +&= +k\binom{n}{k} +B(k+2,n-k+1) += +k\frac{n!}{k!(n-k)!} +\frac{(k+1)!(n-k)!}{(n+2)!} += +\frac{k(k+1)}{(n+1)(n+2)}. +\end{align*} +Die Varianz wird damit +\begin{align} +\operatorname{var}(X_{k:n}) +&= +E(X_{k:n}^2) - E(X_{k:n})^2 +\notag +\\ +& += +\frac{k(k+1)}{(n+1)(n+2)}-\frac{k^2}{(n+1)^2} += +\frac{k(k+1)(n+1)-k^2(n+2)}{(n+1)^2(n+2)} += +\frac{k(n-k+1)}{(n+1)^2(n+2)}. +\label{buch:rekursion:ordnung:eqn:ordnungsstatistik:varianz} +\end{align} +In Abbildung~\ref{buch:rekursion:ordnung:fig:order} ist die Varianz der +Ordnungsstatistik $X_{k:n}$ für $k=7$ und $n=10$ als oranges +Rechteck dargestellt. + +Auch die Varianz kann ganz allgemein für die Beta-Verteilung +bestimmt werden. +Dazu berechnen wir zunächst +\begin{align*} +E(X^2) +&= +\frac{1}{B(a,b)} +\int_0^1 +x^2\cdot x^{a-1}(1-y)^{b-1}\,dx += +\frac{B(a+2,b)}{B(a,b)}. +\end{align*} +Daraus folgt dann +\[ +\operatorname{var}(X) += +E(X^2)-E(X)^2 += +\frac{B(a+2,b)B(a,b)-B(a+1,b)^2}{B(a,b)^2}. +\] + +Die Formel~\eqref{buch:rekursion:ordnung:eqn:ordnungsstatistik:varianz} +besagt auch, dass die Varianz der proportional ist zu $k((n+1)-k)$. +Dieser Ausdruck ist am grössten für $k=(n+1)/2$, die Varianz ist +also grösser für die ``mittleren'' Ordnungstatistiken als für die +extremen $X_{1:n}=\operatorname{min}(X_1,\dots,X_n)$ und +$X_{n:n}=\operatorname{max}(X_1,\dots,X_n)$. + diff --git a/buch/chapters/090-pde/gleichung.tex b/buch/chapters/090-pde/gleichung.tex index 7f65f06..583895d 100644 --- a/buch/chapters/090-pde/gleichung.tex +++ b/buch/chapters/090-pde/gleichung.tex @@ -5,6 +5,7 @@ % \section{Gleichungen und Randbedingungen \label{buch:pde:section:gleichungen-und-randbedingungen}} +\rhead{Gebiete, Gleichungen und Randbedingungen} \subsection{Gebiete, Differentialoperatoren, Randbedingungen} diff --git a/buch/chapters/090-pde/kreis.tex b/buch/chapters/090-pde/kreis.tex index c60fd44..a8cab3e 100644 --- a/buch/chapters/090-pde/kreis.tex +++ b/buch/chapters/090-pde/kreis.tex @@ -5,6 +5,7 @@ % \section{Kreisförmige Membran \label{buch:pde:section:kreis}} +\rhead{Kreisförmige Membran} In diesem Abschnitt soll die Differentialgleichung einer kreisförmigen Membran mit Hilfe der Separationsmethode gelöst werden. Dabei werden die Bessel-Funktionen als Lösungsfunktionen diff --git a/buch/chapters/090-pde/kugel.tex b/buch/chapters/090-pde/kugel.tex index d466e26..ee56316 100644 --- a/buch/chapters/090-pde/kugel.tex +++ b/buch/chapters/090-pde/kugel.tex @@ -5,6 +5,7 @@ % \section{Kugelfunktionen \label{buch:pde:section:kugel}} +\rhead{Kugelfunktionen} Kugelsymmetrische Probleme können oft vorteilhaft in Kugelkoordinaten beschrieben werden. Die Separationsmethode kann auf partielle Differentialgleichungen diff --git a/buch/chapters/090-pde/rechteck.tex b/buch/chapters/090-pde/rechteck.tex index 72e2806..b7dfe11 100644 --- a/buch/chapters/090-pde/rechteck.tex +++ b/buch/chapters/090-pde/rechteck.tex @@ -5,6 +5,7 @@ % \section{Rechteckige Membran \label{buch:pde:section:rechteck}} +\rhead{Rechteckige Membran} Als Beispiel für die Lösung des in Abschnitt~\ref{buch:pde:subsection:eigenwertproblem} aus der Wellengleichung abgeleiteten Eigenwertproblems diff --git a/buch/chapters/090-pde/separation.tex b/buch/chapters/090-pde/separation.tex index 6faceaa..e5e144a 100644 --- a/buch/chapters/090-pde/separation.tex +++ b/buch/chapters/090-pde/separation.tex @@ -5,6 +5,7 @@ % \section{Separationsmethode \label{buch:pde:section:separation}} +\rhead{Separationsmethode} Die Existenz der Lösung einer gewöhnlichen Differentialgleichung ist unter einigermassen milden Bedingungen in der Nähe der Anfangsbedingung garantiert. -- cgit v1.2.1 From e245386805c7bb6d7515154f2bd50aba3a7713a8 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Fri, 18 Mar 2022 08:45:05 +0100 Subject: typos --- buch/chapters/070-orthogonalitaet/gaussquadratur.tex | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'buch') diff --git a/buch/chapters/070-orthogonalitaet/gaussquadratur.tex b/buch/chapters/070-orthogonalitaet/gaussquadratur.tex index 55f9700..b7a5643 100644 --- a/buch/chapters/070-orthogonalitaet/gaussquadratur.tex +++ b/buch/chapters/070-orthogonalitaet/gaussquadratur.tex @@ -135,7 +135,7 @@ p(x)&=x^2\colon& \frac23 &= A_0x_0^2 + A_1x_1^2\\ p(x)&=x^3\colon& 0 &= A_0x_0^3 + A_1x_1^3. \end{aligned} \] -Dividiert man die zweite und vierte Gleichung in der Form +Dividiert man die zweite und dritte Gleichung in der Form \[ \left. \begin{aligned} @@ -155,7 +155,7 @@ x_1=-x_0. \] Indem wir dies in die zweite Gleichung einsetzen, finden wir \[ -0 = A_0x_0 + A_1x_1 = A_0x_1 -A_1x_0 = (A_0-A_1)x_0 +0 = A_0x_0 + A_1x_1 = A_0x_0 -A_1x_0 = (A_0-A_1)x_0 \quad\Rightarrow\quad A_0=A_1. \] -- cgit v1.2.1 From 8ea715ed37adfbaf6369c7fbd889380857b574e0 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Fri, 18 Mar 2022 08:56:45 +0100 Subject: fix Gauss quadrature example --- buch/chapters/070-orthogonalitaet/gaussquadratur.tex | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) (limited to 'buch') diff --git a/buch/chapters/070-orthogonalitaet/gaussquadratur.tex b/buch/chapters/070-orthogonalitaet/gaussquadratur.tex index b7a5643..acfdb1a 100644 --- a/buch/chapters/070-orthogonalitaet/gaussquadratur.tex +++ b/buch/chapters/070-orthogonalitaet/gaussquadratur.tex @@ -135,12 +135,12 @@ p(x)&=x^2\colon& \frac23 &= A_0x_0^2 + A_1x_1^2\\ p(x)&=x^3\colon& 0 &= A_0x_0^3 + A_1x_1^3. \end{aligned} \] -Dividiert man die zweite und dritte Gleichung in der Form +Dividiert man die vierte durch die zweite Gleichung in der Form \[ \left. \begin{aligned} -A_0x_0 &= -A_1x_1\\ -A_0x_0^2 &= -A_1x_1^2 +A_0x_0^3 &= -A_1x_1^3 &\qquad&\text{(vierte Gleichung)}\\ +A_0x_0 &= -A_1x_1 &\qquad&\text{(zweite Gleichung)} \end{aligned} \quad \right\} -- cgit v1.2.1 From 9b8e3f91862b5acd367cf7bbf9d52cb17563390d Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Sat, 26 Mar 2022 22:08:42 +0100 Subject: add new lecture notes --- buch/chapters/110-elliptisch/images/jacobiplots.pdf | Bin 57192 -> 56975 bytes buch/chapters/110-elliptisch/images/jacobiplots.tex | 2 +- 2 files changed, 1 insertion(+), 1 deletion(-) (limited to 'buch') diff --git a/buch/chapters/110-elliptisch/images/jacobiplots.pdf b/buch/chapters/110-elliptisch/images/jacobiplots.pdf index d11bde8..88cf119 100644 Binary files a/buch/chapters/110-elliptisch/images/jacobiplots.pdf and b/buch/chapters/110-elliptisch/images/jacobiplots.pdf differ diff --git a/buch/chapters/110-elliptisch/images/jacobiplots.tex b/buch/chapters/110-elliptisch/images/jacobiplots.tex index 4fc572e..fec04fc 100644 --- a/buch/chapters/110-elliptisch/images/jacobiplots.tex +++ b/buch/chapters/110-elliptisch/images/jacobiplots.tex @@ -31,7 +31,7 @@ \fill[color=gray!50] (-0.2,1.65) rectangle (7.0,2.3); \draw[line width=0.5pt] (-0.2,-6) rectangle (7.0,2.3); \begin{scope}[scale=0.5] -\node at (6.5,{\dy+2}) {$m = #1$}; +\node at (6.5,{\dy+2}) {$k^2 = #1$}; \end{scope} } \def\jacobiplot#1#2#3#4{ -- cgit v1.2.1 From 4a49ccec57384ba582c1c132a33942c938bc1b43 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Mon, 28 Mar 2022 13:26:54 +0200 Subject: new stuff about bessel --- buch/chapters/075-fourier/2d.tex | 19 + buch/chapters/075-fourier/Makefile.inc | 2 + buch/chapters/075-fourier/bessel.tex | 620 +++++++++++++++++++++++++++++++++ buch/chapters/075-fourier/chapter.tex | 3 +- 4 files changed, 643 insertions(+), 1 deletion(-) create mode 100644 buch/chapters/075-fourier/2d.tex create mode 100644 buch/chapters/075-fourier/bessel.tex (limited to 'buch') diff --git a/buch/chapters/075-fourier/2d.tex b/buch/chapters/075-fourier/2d.tex new file mode 100644 index 0000000..cc019c7 --- /dev/null +++ b/buch/chapters/075-fourier/2d.tex @@ -0,0 +1,19 @@ +% +% 2d.tex +% +% (c) 2022 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +% +\section{Zweidimensionale Fourier-Transformation +\label{buch:fourier:section:2d}} +\rhead{Zweidimensionale Fourier-Transformation} + +\subsection{Fourier-Transformation und partielle Differentialgleichungen} + +\subsection{Fourier-Transformation in kartesischen Koordinaten} + +\subsection{Basisfunktionen in Polarkoordinaten} + + + + + diff --git a/buch/chapters/075-fourier/Makefile.inc b/buch/chapters/075-fourier/Makefile.inc index ee9641c..c153dc4 100644 --- a/buch/chapters/075-fourier/Makefile.inc +++ b/buch/chapters/075-fourier/Makefile.inc @@ -5,4 +5,6 @@ # CHAPTERFILES = $(CHAPTERFILES) \ + chapters/075-fourier/bessel.tex \ + chapters/075-fourier/2d.tex \ chapters/075-fourier/chapter.tex diff --git a/buch/chapters/075-fourier/bessel.tex b/buch/chapters/075-fourier/bessel.tex new file mode 100644 index 0000000..7e978f7 --- /dev/null +++ b/buch/chapters/075-fourier/bessel.tex @@ -0,0 +1,620 @@ +% +% bessel.tex +% +% (c) 2022 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +% +\section{Fourier-Transformation und Bessel-Funktionen +\label{buch:fourier:section:fourier-und-bessel}} +\rhead{Fourier-Transformation und Bessel-Funktionen} + +Sei $f\colon \mathbb{R}^2\to\mathbb{C}$ eine auf $\mathbb{R}$ definierte +Funktion. +Die Fourier-Transformation von $f$ ist das Integral +\begin{equation} +(\mathscr{F}f)(u,v) += +F(u,v) += +\frac{1}{2\pi} +\int_{-\infty}^\infty +\int_{-\infty}^\infty +f(x,y) e^{i(xu+yv)} +\,dx\,dy. +\label{buch:fourier:eqn:2dfourier} +\end{equation} +Die Funktionen $e_{u,v}\colon (x,y)\mapsto e^{i(xu+yv)}$ +sind die Eigenfunktionen des Laplace-Operators in kartesischen Koordinaten, +sie erfüllen +\[ +\Delta e_{u,v} = (u^2+v^2) \Delta e_{u,v}. +\] +Die Fourier-Integrale sind die Skalarprodukte +\[ +(\mathscr{F}f)(u,v) += +\langle +e_{u,v}, +f +\rangle, +\] +wobei das Skalarprodukt durch +\[ +\langle f,g\rangle += +\int_{-\infty}^\infty +\int_{-\infty}^\infty +\overline{f(x)} g(x) +\,dx\,dy +\] +definiert ist. + +Jede Funktion in der Ebene kann auch in Polarkoordinaten ausgedrückt werden. +Die kartesischen Koordinaten können mittels +\begin{align*} +x&=r\cos\varphi +y&=r\sin\varphi +\end{align*} +durch die Polarkoordinaten $(r,\varphi)$ ausgedrückt werden. +Wir schreiben +\[ +\tilde{f}(r,\varphi) += +f(r\cos\varphi,r\sin\varphi) +\] +für die Funktion $f$ ausgedrückt in Polarkoordinaten. + +In Polarkoordinaten wird das Skalarprodukt +\[ +\langle f,g\rangle += +\int_0^\infty \int_{0}^{2\pi} e^{in\varphi} +\overline{ +\tilde{f}(r,\varphi) +} +\tilde{g}(r,\varphi) +r\,dr\,d\varphi. +\] +Auch die Fouriertransformation kann jetzt durch Berechnung eines +doppelten Integrals in Polarkoordinaten ermittelt werden. +Ziel dieses Abschnitts ist zu zeigen, dass auch diese Berechnung auf +Bessel-Funktionen führt. +Im Gegenzug werden sich neue Eigenschaften und Darstellungen derselben +ergeben. + + +\subsection{Berechnung der Fourier-Transformation in Polarkoordinaten} +Die Fourier-Transformation $(\mathscr{F}f)(u,v)$ ist eine Funktion +$\mathbb{R}^2\to\mathbb{C}$, die vom Wellenvektor $(u,v)$ abhängt. +Auch dieser Vektor kann in Polarkoordinaten ausgedrückt werden. +Für die Polarkoordinaten in der Wellenvektor-Ebene soll die Bezeichnung +$(R,\vartheta)$ verwendet werden, was auf die Transformationsgleichungen +\begin{align*} +u&=R\cos\vartheta\\ +v&=R\sin\vartheta +\end{align*} +führt. +Im Exponenten der Exponentialfunktion +des Fourier-Integrals~\eqref{buch:fourier:eqn:2dfourier} +steht der Ausdruck +\[ +xu+yv += +r\cos\varphi\cdot R\cos\vartheta ++ +r\sin\varphi\cdot R\sin\vartheta += +rR\cos(\varphi-\vartheta). +\] +Mit diesen Bezeichnungen wird das +Fourier-Integral~\eqref{buch:fourier:eqn:2dfourier} +zu +\begin{align} +\tilde{F}(R,\vartheta) +&= +\frac{1}{2\pi} +\int_{0}^{\infty} +\int_{0}^{2\pi} +f(r\cos\varphi,r\sin\varphi) +e^{irR\cos(\varphi-\vartheta)} +\,d\varphi\,r\, dr +\notag +\\ +&= +\frac{1}{2\pi} +\int_{0}^{\infty} +\int_{0}^{2\pi} +\tilde{f}(r,\varphi) +e^{irR\cos(\varphi-\vartheta)} +\,d\varphi\,r\, dr. +\label{buch:fourier:eqn:fouriertrafopolar} +\end{align} +Die partielle Funktion $\varphi\mapsto \tilde{f}(r,\varphi)$ +ist eine $2\pi$-periodische Funktion, sie lässt sich also als +komplexe Fourier-Reihe +\begin{equation} +\tilde{f}(r,\varphi) += +\sum_{n\in\mathbb{Z}} \hat{f}_n(r) e^{in\varphi} +\label{buch:fourier:eqn:fourierkoef} +\end{equation} +schreiben, die Funktionen $\hat{f}_n(r)$ sind die komplexen +Fourier-Koeffizienten. +Setzt man \eqref{buch:fourier:eqn:fourierkoef} in die Fourier-Transformation +\eqref{buch:fourier:eqn:fouriertrafopolar} ein, erhält man +\begin{align*} +\tilde{F}(R,\vartheta) +&= +\sum_{n\in\mathbb{Z}} +\int_0^\infty +\hat{f}_n(r) +\frac{1}{2\pi} +\int_0^{2\pi} +e^{in\varphi+irR\cos(\varphi-\vartheta)} +\,d\varphi +\, +r\,dr. +\end{align*} +Der Exponent im inneren Integral kann als +\[ +in\varphi+irR\cos(\varphi-\vartheta) += +i(n(\varphi-\vartheta)+rR\cos(\varphi-\vartheta)) ++ +in\vartheta, +\] +oder im Integral als +\[ +\tilde{F}(R,\vartheta) += +\sum_{n\in\mathbb{Z}} +\int_0^\infty +\hat{f}_n(r) +\frac{1}{2\pi} +\int_0^{2\pi} +e^{in(\varphi-\vartheta)+irR\cos(\varphi-\vartheta)} +e^{in\vartheta} +\,d\varphi +\, +r\,dr +\] +geschrieben werden. +Der zweite Exonentialfaktor hängt nicht von $\varphi$ ab und kann daher +aus dem Integral herausgezogen werden. +Der erste Exponentialfaktor hängt nur von $\varphi-\vartheta$ ab. +Da die Exponentialfunktion $2\pi$-periodisch ist, hat die Verschiebung +um $\vartheta$ keinen Einfluss auf den Wert des Integrals. +Die Fourier-Transformation ist daher auch +\[ +\tilde{F}(R,\vartheta) += +\sum_{n\in\mathbb{Z}} +\int_0^\infty +\hat{f}_n(r) +e^{in\vartheta} +\underbrace{ +\frac{1}{2\pi} +\int_0^{2\pi} +e^{in\varphi+irR\cos\varphi} +\,d\varphi +}_{\displaystyle =:F_n(rR)} +\, +r\,dr. +\] +Die Beziehung zu den Besselfunktionen können wir daraus herstellen, +indem wir zunächst $\xi = rR$ abkürzen und dann das innere Integral +\begin{equation} +F_n(\xi) += +\frac{1}{2\pi} +\int_{0}^{2\pi} +e^{in\varphi+i\xi\cos\varphi} +\,d\varphi += +\frac{1}{2\pi} +\int_{0}^{2\pi} +e^{in\varphi}e^{i\xi\cos\varphi} +\,d\varphi +\label{buch:fourier:eqn:Fncosphi} +\end{equation} +auswerten. +Exponentialfunktion als Potenzreihe entwickeln: +\[ +F_n(\xi) += +\frac{1}{2\pi} +\int_0^{2\pi} +e^{in\varphi} +\sum_{k=0}^\infty +\frac{ +i^k\xi^k \cos^k\varphi +}{k!} +\,d\varphi += +\sum_{k=0}^\infty +\frac{i^k\xi^k}{k!} +\underbrace{ +\frac{1}{2\pi} +\int_0^{2\pi} +e^{in\varphi} +\cos^k\varphi +\,d\varphi}_{\displaystyle =c_{n,k}}. +\] +Das Integral auf der rechten Seite ist im Wesentlichen ein +Fourier-Koeffizient der Funktion $\varphi\mapsto \cos^k\varphi$. + +\subsubsection{Berechnung der Fourier-Koeffizienten von $\cos^k\varphi$} +Indem man die Kosinus-Funktion als die Linearkombination +\[ +\cos\varphi += +\frac{e^{i\varphi}+e^{-i\varphi}}2 +\] +von Exponentialfunktionen ausdrückt, kann man auch die $k$-te Potenz +mit Hilfe des binomischen Satzes als +\[ +\cos^k\varphi += +\sum_{m=0}^k +\frac{1}{2^k} +\binom{k}{m} +e^{im\varphi}e^{i(m-k)\varphi} += +\sum_{m=0}^k +\frac{1}{2^k} +\binom{k}{m} +e^{i(2m-k)\varphi} +\] +ausdrücken. +Der Fourier-Koeffizient von $\cos^k\varphi$ ist daher das Integral +\begin{align*} +c_{n,k} +&= +\frac{1}{2\pi} +\int_0^{2\pi} +e^{in\varphi}\cos^k\varphi\,d\varphi +\\ +&= +\frac{1}{2^k} +\sum_{m=0}^k +\binom{k}{m} +\frac{1}{2\pi} +\int_0^{2\pi} +e^{in\varphi}e^{i(2m-k)\varphi} +\,d\varphi +\\ +&= +\frac{1}{2^k} +\sum_{m=0}^k +\binom{k}{m} +\frac{1}{2\pi} +\int_0^{2\pi} +e^{i(2m-k+n)\varphi} +\,d\varphi. +\end{align*} +Für $2m-k+n=0$ ist das Integral ein Integral der Funktion $1$ über +ein Intervall der Länge $2\pi$, zusammen mit dem Faktor $1/2\pi$ hat +es daher den Wert $1$. +Für $2m-k+n\ne 0$ ist das Integral +\[ +\frac{1}{2\pi} +\int_0^{2\pi} +e^{i(2m-k+n)\varphi} +\,d\varphi += +\frac{1}{i} +\biggl[ +\frac{e^{i(2m-k+n)\varphi}}{2m-k+n} +\biggr]_0^{2\pi} += +0 +\] +weil die Exponentialfunktion $2\pi$-periodisch ist. +Nur für $k=2m+n$ ergibt sich ein nicht verschwindender +Fourier-Koeffizient. +Eine Summe über $k\in\mathbb{N}$ kann daher auch als Summe über +$m\in\mathbb{N}$ interpretiert werden, in der $k$ durch die Formel +$k=2m+n$ gegeben wird. +Mit dieser Konvention wird +\[ +c_{n,k} += +c_{n,2m+n} +%= +%\frac{1}{2\pi} +%\int_0^{2\pi} +%e^{-i(2m+n)\varphi} +%\cos^{2m+n}\varphi +%\,d\varphi += +\frac{1}{2^{2m+n}} +\binom{2m+n}{m} +\] +schreiben lässt. + +\subsubsection{Berechnung von $F_n(\xi)$} +Die Reihe für $F_n(\xi)$ lässt sich weiter vereinfachen. +Wir verwenden wieder die Tatsache, dass sich nur für $n=-2m-k$ +ein Beitrag ergibt. +Dies bedeutet, dass $k=2m+n$ sein muss, die Summe kann damit als +Summe über $m$ statt über $k$ geschrieben werden. +Somit ist +\begin{align*} +F_n(\xi) +&= +\sum_{k=0}^\infty +\frac{i^k\xi^k}{k!} +c_{n,k} += +\sum_{m=0}^\infty +\frac{i^{2m+n}\xi^{2m+n}}{(2m+n)!} +c_{n,2m+n} +\\ +&= +\sum_{m=0}^\infty +\frac{1}{2^{2m+n}} +\binom{2m+n}{m} +\frac{i^{2m+n}\xi^{2m+n}}{(2m+n)!} +\\ +&= +i^n +\sum_{m=0}^\infty +\frac{(-1)^m}{(2m+n)!} +\frac{(2m+n)!}{m!\,(2m+n-m)!} +\biggl(\frac{\xi}{2}\biggr)^{2m+n} +\\ +&= +i^n +\sum_{m=0}^\infty +\frac{(-1)^m} +{m!\,\Gamma(m+n+1)} +\biggl(\frac{\xi}{2}\biggr)^{2m+n} += +i^n J_n(\xi). +\end{align*} +Die Funktionen $F_n(\xi)$ sind daher bis auf einen Phasenfaktor der +Wert $J_n(\xi)$ einer Bessel-Funktion. + +\subsubsection{Berechnung der Fourier-Transformation mit Bessel-Funktionen} +Mit allen oben zusammengestellten Notationen kann die Fourier-Transformation +jetzt in Polarkoordinaten als +\[ +\tilde{F}(R,\vartheta) += +\sum_{n\in\mathbb{Z}} +e^{in\vartheta} +\int_0^\infty +\hat{f}_n(r) +i^n +J_n(rR) +r\,dr +\] +geschrieben werden. +Dies hat tatsächlich die Form eines Skalarproduktes der Funktion +$\tilde{f}(r,\varphi)$ mit einer Funktion der Form +\[ +\tilde{e}_{n,R}(r,\varphi) += +e^{in\varphi} +J_n(rR). +\] +Letzeres sind die in Abschnitt~\ref{buch:fourier:section:2d} +versprochenen Basisfunktionen. + +\subsubsection{Fourier-Reihe von $e^{i\xi\cos\varphi}$} +Die Funktionen $F_n(\xi)$ sind wegen +\[ +F_n(\xi) += +\frac{1}{2\pi} +\int_0^{2\pi} +e^{in\varphi} +e^{i\xi\cos\varphi} +\,d\varphi, +\] +daraus kann man die Fourier-Reihe von $e^{i\xi\cos\varphi}$ +berechnen, dies wird im folgenden Satz durchgeführt. + + +\begin{satz} +\label{buch:fourier:satz:expinphi} +Die komplexe Fourier-Reihe der Funktion +$\varphi\mapsto \exp(i\xi\cos\varphi)$ +ist +\begin{align} +e^{i\xi\cos\varphi} +&= +J_0(\xi) ++ +2\sum_{n=1}^\infty i^n J_n(\xi) \cos n\varphi. +\label{buch:fourier:eqn:expinphicomplex}. +\intertext{Real- und Imaginärteil davon sind die Fourier-Reihen} +\cos(\xi\cos\varphi) +&= +J_0(\xi) + 2\sum_{m=1}^\infty (-1)^m J_{2m}(\xi) \cos2m\varphi +\label{buch:fourier:eqn:expinphireal} +\\ +\sin(\xi\cos\varphi) +&= +2\sum_{m=0}^\infty (-1)^m J_{2m+1}(\xi) \cos(2m+1)\varphi. +\label{buch:fourier:eqn:expinphiimaginary} +\end{align} +\end{satz} + +\begin{proof}[Beweis] +Die Fourier-Koeffizienten $F_n(\xi)$ der Funktion $e^{i\xi\cos\varphi}$ +führen auf die Fourier-Reihe +\begin{align*} +e^{i\xi\cos\varphi} +&= +\sum_{n\in\mathbb{Z}} F_n(\xi) e^{in\varphi} += +\sum_{n\in\mathbb{Z}} i^n J_n(\xi) e^{in\varphi}. +\end{align*} +Terme mit $\pm n$ können wegen +\[ +\left. +\begin{aligned} +J_{-n}(\xi) &= (-1)^n J_n(\xi) +\\ +i^{-n}&=(-1)^n i^n +\end{aligned} +\quad +\right\} +\qquad\Rightarrow\qquad +i^{-n}J_{-n}(\xi) = i^n J_n(\xi) +\] +zusammengefasst werden, auf diese Weise erhält man +\begin{align*} +e^{i\xi\cos\varphi} +&= +J_0(\xi) ++ +\sum_{n=1}^\infty i^n J_n(\xi) (e^{in\varphi}+e^{-in\varphi}) += +2\sum_{n=1}^\infty i^n J_n(\xi) \cos n\varphi. +\end{align*} +Dies beweist +\eqref{buch:fourier:eqn:expinphicomplex}. + +Indem man Real- und Imaginärteil trennt, kann man daraus auch +die Fourier-Reihen von $\cos(\xi\cos\varphi)$ und +$\sin(\xi\cos\varphi)$ gewinnen, sie sind +\begin{align*} +\exp(\xi\cos\varphi) +&= +J_0(\xi) + 2\sum_{n=1}^\infty i^{n} J_{n}(\xi) \cos n\varphi +\\ +&= +J_0(\xi) ++ +2\sum_{m=1}^\infty i^{2m}J_{2m}(\xi)\cos 2m\varphi ++ +2\sum_{m=0}^\infty i^{2m+1}J_{2m+1}(\xi)\cos(2m+1)\varphi +\\ +&= +J_0(\xi) ++ +2\sum_{m=1}^\infty (-1)^{m}J_{2m}(\xi)\cos 2m\varphi ++ +2i\sum_{m=0}^\infty (-1)^{m}J_{2m+1}(\xi)\cos(2m+1)\varphi +\\ +\cos(\xi\cos\varphi) +&= +J_0(\xi) ++ +2\sum_{m=1}^\infty (-1)^{m}J_{2m}(\xi)\cos 2m\varphi +\\ +\sin(\xi\cos\varphi) +&= +2\sum_{m=0}^\infty (-1)^m J_{2m+1}(\xi) \cos(2m+1)\varphi. +\end{align*} +Damit sind auch die Formeln +\eqref{buch:fourier:eqn:expinphireal} +und +\eqref{buch:fourier:eqn:expinphiimaginary} +für die reellen Fourier-Reihen bewiesen. +\end{proof} + +% +% Integraldarstellung der Bessel-Funktion +% +\subsection{Integraldarstellung der Bessel-Funktion} +Aus \eqref{buch:fourier:eqn:Fncosphi} kann jetzt die Integraldarstelltung +der Bessel-Funktionen gewonnen werden. +Dazu substituiert man $\varphi$ durch $\tau$ mit +$\varphi = \frac{\pi}2-\tau$ +oder +$\tau=\frac{\pi}2-\varphi$ +und $d\tau = -d\varphi$ +im Integral und berechnet +\begin{align*} +J_n(\xi) +&= +(-i)^n +\frac{1}{2\pi} +\int_0^{2\pi} +e^{in\varphi+i\xi \cos\varphi} +\,d\varphi +\\ +&= +- +(-i)^n +\frac{1}{2\pi} +\int_{\frac{\pi}2}^{-\frac{3\pi}2} +e^{in(\frac{\pi}2-\tau) + i\xi\cos(\frac{\pi}2-\tau)} +\,d\tau +\\ +&= +(-i)^n +\frac{1}{2\pi} +\int^{\frac{\pi}2}_{-\frac{3\pi}2} +i^n +e^{-in\tau + i\xi\sin\tau)} +\,d\tau. +\intertext{Da der Integrand $2\pi$-periodisch ist, kann das +Integrationsintervall auf $[-\pi,\pi]$ verschoben werden, was} +&= +\frac{1}{2\pi} +\int_{-\pi}^{\pi} +e^{-in\tau + i\xi\sin\tau)} +\,d\tau. +\intertext{ergibt. +Das Integral kann in zwei Integrale} +&= +\frac{1}{2\pi} +\int_0^\pi +e^{-in\tau + i\xi\sin\tau} +\,d\tau ++ +\frac{1}{2\pi} +\int_0^\pi +e^{in\tau - i\xi\sin\tau} +\,d\tau +\intertext{aufgeteilt werden, +} +&= +\frac{1}{\pi} +\int_0^\pi +\frac{ +e^{-in\tau + i\xi\sin\tau} ++ +e^{in\tau - i\xi\sin\tau} +}{2} +\,d\tau +\\ +&= +\frac{1}{\pi} +\int_0^\pi +\frac{ +e^{i(-n\tau + \xi\sin\tau)} ++ +e^{-i(-n\tau + \xi\sin\tau)} +}{2} +\,d\tau +\\ +&= +\frac{1}{\pi} +\int_0^\pi +\cos(n\tau - \xi\sin\tau) +\,d\tau. +\end{align*} +Damit haben wir den folgenden Satz bewiesen: + +\begin{satz}[Integraldarstelltung der Bessel-Funktionen] +\label{buch:fourier:satz:bessel-integraldarstellung} +Die Bessel-Funktionen $J_n$ mit ganzzahliger Ordnung $n$ haben +die Integraldarstellung +\begin{equation} +J_n(\xi) += +\frac{1}{\pi} +\int_0^\pi +\cos(n\tau - \xi\sin\tau) +\,d\tau. +\label{buch:fourier:eqn:bessel-integraldarstellung} +\end{equation} +\end{satz} + + + + diff --git a/buch/chapters/075-fourier/chapter.tex b/buch/chapters/075-fourier/chapter.tex index 341d8df..681a1c0 100644 --- a/buch/chapters/075-fourier/chapter.tex +++ b/buch/chapters/075-fourier/chapter.tex @@ -13,7 +13,8 @@ führen zu neuen speziellen Funktionen. In diesem Kapitel soll als Beispiel die Fourier-Transformation der Bessel-Funktionen untersucht werden. -%\input{chapters/075-fourier/bessel.tex} +\input{chapters/075-fourier/2d.tex} +\input{chapters/075-fourier/bessel.tex} %\section{TODO} %\begin{itemize} -- cgit v1.2.1 From 45241da09f3dc8b8d700d505edea2d38a26a517c Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Fri, 1 Apr 2022 17:45:58 +0200 Subject: add Vorlesung Komplexe Funktionen --- buch/chapters/110-elliptisch/Makefile.inc | 1 + buch/chapters/110-elliptisch/chapter.tex | 12 ++++++------ 2 files changed, 7 insertions(+), 6 deletions(-) (limited to 'buch') diff --git a/buch/chapters/110-elliptisch/Makefile.inc b/buch/chapters/110-elliptisch/Makefile.inc index 0ca1392..538db68 100644 --- a/buch/chapters/110-elliptisch/Makefile.inc +++ b/buch/chapters/110-elliptisch/Makefile.inc @@ -8,4 +8,5 @@ CHAPTERFILES = $(CHAPTERFILES) \ chapters/110-elliptisch/ellintegral.tex \ chapters/110-elliptisch/jacobi.tex \ chapters/110-elliptisch/lemniskate.tex \ + chapters/110-elliptisch/uebungsaufgaben/001.tex \ chapters/110-geometrie/chapter.tex diff --git a/buch/chapters/110-elliptisch/chapter.tex b/buch/chapters/110-elliptisch/chapter.tex index a03ce24..e09fa53 100644 --- a/buch/chapters/110-elliptisch/chapter.tex +++ b/buch/chapters/110-elliptisch/chapter.tex @@ -20,11 +20,11 @@ aufgebaute Integrale in dieser Familie zu finden. \input{chapters/110-elliptisch/jacobi.tex} \input{chapters/110-elliptisch/lemniskate.tex} -%\section*{Übungsaufgaben} -%\rhead{Übungsaufgaben} -%\aufgabetoplevel{chapters/020-exponential/uebungsaufgaben} -%\begin{uebungsaufgaben} +\section*{Übungsaufgaben} +\rhead{Übungsaufgaben} +\aufgabetoplevel{chapters/110-elliptisch/uebungsaufgaben} +\begin{uebungsaufgaben} %\uebungsaufgabe{0} -%\uebungsaufgabe{1} -%\end{uebungsaufgaben} +\uebungsaufgabe{1} +\end{uebungsaufgaben} -- cgit v1.2.1 From f4ce26a24fbb50621ca52316209bbffd25a60794 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Sat, 2 Apr 2022 08:38:00 +0200 Subject: add missing exercise --- buch/chapters/110-elliptisch/uebungsaufgaben/1.tex | 312 +++++++++++++++++++++ .../110-elliptisch/uebungsaufgaben/Makefile | 8 + .../uebungsaufgaben/anharmonisch.pdf | Bin 0 -> 19279 bytes .../uebungsaufgaben/anharmonisch.tex | 62 ++++ 4 files changed, 382 insertions(+) create mode 100644 buch/chapters/110-elliptisch/uebungsaufgaben/1.tex create mode 100644 buch/chapters/110-elliptisch/uebungsaufgaben/Makefile create mode 100644 buch/chapters/110-elliptisch/uebungsaufgaben/anharmonisch.pdf create mode 100644 buch/chapters/110-elliptisch/uebungsaufgaben/anharmonisch.tex (limited to 'buch') diff --git a/buch/chapters/110-elliptisch/uebungsaufgaben/1.tex b/buch/chapters/110-elliptisch/uebungsaufgaben/1.tex new file mode 100644 index 0000000..8e4b39f --- /dev/null +++ b/buch/chapters/110-elliptisch/uebungsaufgaben/1.tex @@ -0,0 +1,312 @@ +In einem anharmonische Oszillator oszilliert eine Masse $m$ unter dem +Einfluss einer Kraft, die nach dem Gesetz +\[ +F(x) = -\kappa x + \delta x^3 +\] +von der Auslenkung aus der Ruhelage abhängt. +Nehmen Sie im Folgenden an, dass $\delta >0$ ist, +dass also die rücktreibende Kraft $F(x)$ kleiner ist als bei einem +harmonischen Oszillator. +Ziel der folgenden Teilaufgaben ist, die Lösung $x(t)$ schrittweise +dadurch zu bestimmen, dass die Bewegungsgleichung in die Differentialgleichung +der Jacobischen elliptischen Funktion $\operatorname{sn}(u,k)$ umgeformt +wird. +\begin{teilaufgaben} +\item +Berechnen Sie die Auslenkung $x_0$, bei der die rücktreibende Kraft +verschwindet. +Eine beschränkte Schwingung kann diese Amplitude nicht überschreiten. +\item +Berechnen Sie die potentielle Energie in Abhängigkeit von der +Auslenkung. +\item +\label{buch:1101:basic-dgl} +Formulieren Sie den Energieerhaltungssatz für die Gesamtenergie $E$ +dieses Oszillators. +Leiten Sie daraus eine nichtlineare Differentialgleichung erster Ordnung +for den anharmonischen Oszillator ab, die sie in der Form +$\frac12m\dot{x}^2 = f(x)$ schreiben. +\item +Die Amplitude der Schwingung ist derjenige $x$-Wert, für den die +Geschwindigkeit verschwindet. +Leiten Sie die Amplitude aus der Differentialgleichung von +\ref{buch:1101:basic-dgl} ab. +Sie erhalten zwei Werte $x_{\pm}$, wobei der kleinere $x_-$ +die Amplitude einer beschränkten Schwingung beschreibt, +während die $x_+$ die minimale Ausgangsamplitude einer gegen +$\infty$ divergenten Lösung ist. +\item +Rechnen Sie nach, dass +\[ +\frac{x_+^2+x_-^2}{2} += +x_0^2 +\qquad\text{und}\qquad +x_-^2x_+^2 += +\frac{4E}{\delta}. +\] +\item +Faktorisieren Sie die Funktion $f(x)$ in der Differentialgleichung +von Teilaufgabe c) mit Hilfe der in Teilaufgabe d) bestimmten +Nullstellen $x_{\pm}^2$. +\item +Dividieren Sie die Differentialgleichung durch $x_-^2$, schreiben +Sie $X=x/x_-$ und bringen Sie die Differentialgleichung in die +Form +\begin{equation} +A \dot{X}^2 += +(1-X^2) +(1-k^2X^2), +\label{buch:1101:eqn:dgl3} +\end{equation} +wobei $k^2=x_-^2/x_+^2$ und $A$ geeignet gewählt werden müssen. +\item +\label{buch:1101:teilaufgabe:dgl3} +Verwenden Sie $t(\tau) = \alpha\tau$ +und +$Y(\tau)=X(t(\tau))$ um eine Differentialgleichung für die Funktion +$Y(\tau)$ zu gewinnen, die die Form der Differentialgleichung +von $\operatorname{sn}(u,k)$ hat, für die also $A=0$ in +\eqref{buch:1101:eqn:dgl3} ist. +\item +Verwenden Sie die Lösung $\operatorname{sn}(u,k)$ der in +\ref{buch:1101:teilaufgabe:dgl3} erhaltenen Differentialgleichung, +um die Lösung $x(t)$ der ursprünglichen Gleichung aufzuschreiben. +\end{teilaufgaben} + +\begin{loesung} +\begin{figure} +\centering +\includegraphics{chapters/110-elliptisch/uebungsaufgaben/anharmonisch.pdf} +\caption{Rechte Seite der Differentialgleichung +\eqref{buch:1101:eqn:dglf}. +Eine beschränkte Lösung bewegt sich im Bereich $xx_+$ die Kraft abstossend ist und zu einer +divergenten Lösung führt. +\label{buch:1101:fig:potential} +} +\end{figure} +\begin{teilaufgaben} +\item +Wegen +\[ +F(x) += +-\kappa x\biggl(1-\frac{\delta}{\kappa}x^2\biggr) += +-Ix +\biggl(1-\sqrt{\frac{\delta}{\kappa}}x\biggr) +\biggl(1+\sqrt{\frac{\delta}{\kappa}}x\biggr) +\] +folgt, dass die rücktreibende Kraft bei der Auslenkung $\pm x_0$ mit +\[ +x_0^2 += +\frac{\kappa}{\delta} +\qquad\text{oder}\qquad +x_0 = \sqrt{\frac{\kappa}{\delta}} +\] +verschwindet. +\item +Die potentielle Energie ist die Arbeit, die gegen die rücktreibende Kraft +geleistet wird, um die Auslenkung $x$ zu erreichen. +Sie entsteht durch Integrieren der Kraft über +das Auslenkungsinterval, also +\[ +E_{\text{pot}} += +- +\int_0^x F(\xi) \,d\xi += +\int_0^x \kappa\xi-\delta\xi^3\,d\xi += +\biggl[ +\kappa\frac{\xi^2}{2} +- +\delta +\frac{\xi^4}{4} +\biggr]_0^x += +\kappa\frac{x^2}{2} +- +\delta\frac{x^4}{4}. +\] +\item +Die kinetische Energie ist gegeben durch +\[ +E_{\text{kin}} += +\frac12m\dot{x}^2. +\] +Die Gesamtenergie ist damit +\[ +E += +\frac12m\dot{x}^2 ++ +\kappa +\frac{x^2}{2} +- +\delta +\frac{x^4}{4}. +\] +Die verlangte Umformung ergibt +\begin{align} +\frac12m\dot{x}^2 +&= +E +- +\kappa\frac{x^2}{2} ++ +\delta\frac{x^4}{4} +\label{buch:1101:eqn:dglf} +\end{align} +als Differentialgleichung für $x$. +Die Ableitung $\dot{x}$ hat positives Vorzeichen wenn die Kraft +abstossend ist und negatives Vorzeichen dort, wo die Kraft anziehend ist. +% +\item +Die Amplitude der Schwingung ist derjenige $x$-Wert, für den +die Geschwindigkeit verschwindet, also eine Lösung der Gleichung +\[ +0 += +\frac{2E}{m} -\frac{\kappa}{m}x^2 + \frac{\delta}{2m}x^4. +\] +Der gemeinsame Nenner $m$ spielt offenbar keine Rolle. +Die Gleichung hat die zwei Lösungen +\[ +x_{\pm}^2 += +\frac{\kappa \pm \sqrt{\kappa^2-4E\delta}}{\delta} += +\frac{\kappa}{\delta} +\pm +\sqrt{ +\biggl(\frac{\kappa}{\delta}\biggr)^2 +- +\frac{4E}{\delta} +}. +\] +Die Situation ist in Abbildung~\ref{buch:1101:fig:potential} +Für $x>x_+$ ist die Kraft abstossend, die Lösung divergiert. +Die Lösung mit dem negativen Zeichen $x_-$ bleibt dagegen beschränkt, +dies ist die Lösung, die wir suchen. + +\item +Die beiden Formeln ergeben sich aus den Regeln von Vieta für die +Lösungen einer quadratischen Gleichungg der Form $x^4+px^2+q$. +Die Nullstellen haben den Mittelwert $-p/2$ und das Produkt $q$. + +\item +Die rechte Seite der Differentialgleichung lässt sich mit Hilfe +der beiden Nullstellen $x_{\pm}^2$ faktorisieren und bekommt die Form +\[ +\frac12m\dot{x}^2 += +\frac{\delta}{4}(x_+^2-x^2)(x_-^2-x^2). +\] + +\item +Indem die ganze Gleichung durch $x_-^2$ dividiert wird, entsteht +\[ +\frac12m +\biggl(\frac{\dot{x}}{x_-}\biggr)^2 += +\frac{\delta}{4} +(x_+^2-x^2) +\biggl(1-\frac{x^2}{x_-^2}\biggr). +\] +Schreiben wir $X=x/x_-$ wird daraus +\[ +\frac1{2}m\dot{X}^2 += +\frac{\delta}{4} +\biggl(x_+^2-x_-^2 X^2\biggr) +(1-X^2). +\] +Durch Ausklammern von $x_+^2$ im ersten Faktor wir daraus +\[ +\frac1{2}m\dot{X}^2 += +\frac{\delta}{4} +x_+^2 +\biggl(1-\frac{x_-^2}{x_+^2} X^2\biggr) +(1-X^2). +\] +Mit der Schreibweise $k^2 = x_-^2/x_+^2$ wird die Differentialgleichung +zu +\begin{equation} +\frac{2m}{\delta x_+^2} \dot{X}^2 += +(1-X^2)(1-k^2X^2), +\label{buch:1101:eqn:dgl2} +\end{equation} +was der Differentialgleichung für die Jacobische elliptische Funktion +$\operatorname{sn}(u,k)$ bereits sehr ähnlich sieht. +\item +Bis auf den Faktor vor $\dot{X}^2$ ist +\eqref{buch:1101:eqn:dgl2} +die Differentialgleichung +von +$\operatorname{sn}(u,k)$. +Um den Faktor zum Verschwinden zu bringen, schreiben wir +$t(\tau) = \alpha\tau$. +Die Ableitung von $Y(\tau)=X(t(\tau))$ nach $\tau$ ist +\[ +\frac{dY}{d\tau} += +\dot{X}(t(\tau))\frac{dt}{d\tau} += +\alpha +\dot{X}(t(\tau)) +\qquad\Rightarrow\qquad +\frac{1}{\alpha^2}\frac{dY}{d\tau} += +\dot{X}(t(\tau)). +\] +Die Differentialgleichung für $Y(\tau)$ ist +\[ +\frac{2mk^2}{\delta x_+^2\alpha^2} +\frac{dY}{d\tau} += +(1-Y^2)(1-k^2Y^2). +\] +Der Koeffizient vor der Ableitung wird $1$, wenn man +\[ +\alpha^2 += +\frac{2mk^2}{\delta x_+^2} +\] +wählt. +Diese Differentialgleichug hat die Lösung +\[ +Y(\tau) = \operatorname{sn}(\tau,k). +\] +\item +Indem man die gefunden Grössen einsetzt kann man jetzt die Lösung +der Differentialgleichung in geschlossener Form als +\begin{align*} +x(t) +&= +x_- X(t) += +x_- \operatorname{sn}\biggl( +t\sqrt{\frac{\delta x_+^2}{2mk^2} } +,k +\biggr) +\end{align*} +Das Produkt $\delta x_+^2$ kann auch als +\[ +\delta x_+^2 += +\kappa+\sqrt{\kappa -4\delta E} +\] +geschrieben werden. +\qedhere +\end{teilaufgaben} +\end{loesung} + + diff --git a/buch/chapters/110-elliptisch/uebungsaufgaben/Makefile b/buch/chapters/110-elliptisch/uebungsaufgaben/Makefile new file mode 100644 index 0000000..0ca5234 --- /dev/null +++ b/buch/chapters/110-elliptisch/uebungsaufgaben/Makefile @@ -0,0 +1,8 @@ +# +# Makefile +# +# (c) 2022 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +# + +anharmonisch.pdf: anharmonisch.tex + pdflatex anharmonisch.tex diff --git a/buch/chapters/110-elliptisch/uebungsaufgaben/anharmonisch.pdf b/buch/chapters/110-elliptisch/uebungsaufgaben/anharmonisch.pdf new file mode 100644 index 0000000..4b00f4d Binary files /dev/null and b/buch/chapters/110-elliptisch/uebungsaufgaben/anharmonisch.pdf differ diff --git a/buch/chapters/110-elliptisch/uebungsaufgaben/anharmonisch.tex b/buch/chapters/110-elliptisch/uebungsaufgaben/anharmonisch.tex new file mode 100644 index 0000000..a00c393 --- /dev/null +++ b/buch/chapters/110-elliptisch/uebungsaufgaben/anharmonisch.tex @@ -0,0 +1,62 @@ +% +% anharmonisch.tex -- Potential einer anharmonischen Schwingung +% +% (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +% +\documentclass[tikz]{standalone} +\usepackage{amsmath} +\usepackage{times} +\usepackage{txfonts} +\usepackage{pgfplots} +\usepackage{csvsimple} +\usetikzlibrary{arrows,intersections,math} +\begin{document} +\def\skala{1} +\definecolor{darkgreen}{rgb}{0,0.6,0} +\begin{tikzpicture}[>=latex,thick,scale=\skala] + +\def\E{3} +\def\K{0.2} +\def\D{0.0025} + +\pgfmathparse{sqrt(\K/\D)} +\xdef\xnull{\pgfmathresult} + +\pgfmathparse{sqrt((\K+sqrt(\K*\K-4*\E*\D))/\D)} +\xdef\xplus{\pgfmathresult} +\pgfmathparse{sqrt((\K-sqrt(\K*\K-4*\E*\D))/\D)} +\xdef\xminus{\pgfmathresult} + +\def\xmax{13} + +\fill[color=darkgreen!20] (0,-1.5) rectangle (\xminus,4.7); +\node[color=darkgreen] at ({0.5*\xminus},4.7) [below] {anziehende Kraft\strut}; + +\fill[color=orange!20] (\xplus,-1.5) rectangle (\xmax,4.7); +\node[color=orange] at ({0.5*(\xplus+\xmax)},4.7) [below] {abstossende\strut}; +\node[color=orange] at ({0.5*(\xplus+\xmax)},4.3) [below] {Kraft\strut}; + +\node[color=gray] at (\xnull,4.7) [below] {verbotener Bereich\strut}; + +\draw (-0.1,\E) -- (0.1,\E); +\node at (-0.1,\E) [left] {$E$}; + +\draw[color=red,line width=1pt] + plot[domain=0:13,samples=100] + ({\x},{\E-(0.5*\K-0.25*\D*\x*\x)*\x*\x}); + +\draw[->] (-0.1,0) -- ({\xmax+0.3},0) coordinate[label={$x$}]; +\draw[->] (0,-1.5) -- (0,5) coordinate[label={right:$f(x)$}]; + +\fill[color=blue] (\xminus,0) circle[radius=0.08]; +\node[color=blue] at (\xminus,0) [below left] {$x_-\mathstrut$}; + +\fill[color=blue] (\xplus,0) circle[radius=0.08]; +\node[color=blue] at (\xplus,0) [below right] {$x_+\mathstrut$}; + +\fill[color=blue] (\xnull,0) circle[radius=0.08]; +\node[color=blue] at (\xnull,0) [below] {$x_0\mathstrut$}; + +\end{tikzpicture} +\end{document} + -- cgit v1.2.1