From 0c784258872a10b196382edf8fd3f05be97980c9 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Fri, 18 Feb 2022 22:27:29 +0100 Subject: more elliptic function stuff --- buch/chapters/110-elliptisch/images/Makefile | 7 +- .../chapters/110-elliptisch/images/jacobiplots.pdf | Bin 57192 -> 57192 bytes buch/chapters/110-elliptisch/images/sncnlimit.pdf | Bin 0 -> 24327 bytes buch/chapters/110-elliptisch/images/sncnlimit.tex | 79 +++++ buch/chapters/110-elliptisch/jacobi.tex | 336 ++++++++++++++++++++- 5 files changed, 412 insertions(+), 10 deletions(-) create mode 100644 buch/chapters/110-elliptisch/images/sncnlimit.pdf create mode 100644 buch/chapters/110-elliptisch/images/sncnlimit.tex (limited to 'buch/chapters/110-elliptisch') diff --git a/buch/chapters/110-elliptisch/images/Makefile b/buch/chapters/110-elliptisch/images/Makefile index d1e0afe..68322b6 100644 --- a/buch/chapters/110-elliptisch/images/Makefile +++ b/buch/chapters/110-elliptisch/images/Makefile @@ -4,7 +4,8 @@ # (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule # all: lemniskate.pdf ellipsenumfang.pdf unvollstaendig.pdf rechteck.pdf \ - ellipse.pdf pendel.pdf jacobiplots.pdf jacobidef.pdf jacobi12.pdf + ellipse.pdf pendel.pdf jacobiplots.pdf jacobidef.pdf jacobi12.pdf \ + sncnlimit.pdf lemniskate.pdf: lemniskate.tex pdflatex lemniskate.tex @@ -66,3 +67,7 @@ jacobidef.pdf: jacobidef.tex jacobi12.pdf: jacobi12.tex pdflatex jacobi12.tex + +sncnlimit.pdf: sncnlimit.tex + pdflatex sncnlimit.tex + diff --git a/buch/chapters/110-elliptisch/images/jacobiplots.pdf b/buch/chapters/110-elliptisch/images/jacobiplots.pdf index 8ebd501..d11bde8 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/sncnlimit.pdf b/buch/chapters/110-elliptisch/images/sncnlimit.pdf new file mode 100644 index 0000000..c721427 Binary files /dev/null and b/buch/chapters/110-elliptisch/images/sncnlimit.pdf differ diff --git a/buch/chapters/110-elliptisch/images/sncnlimit.tex b/buch/chapters/110-elliptisch/images/sncnlimit.tex new file mode 100644 index 0000000..2b284fb --- /dev/null +++ b/buch/chapters/110-elliptisch/images/sncnlimit.tex @@ -0,0 +1,79 @@ +% +% sncnlimit.tex -- template for standalon tikz images +% +% (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +% +\documentclass[tikz]{standalone} +\usepackage{amsmath} +\usepackage{times} +\usepackage{txfonts} +\usepackage{pgfplots} +\usepackage{csvsimple} +\definecolor{darkgreen}{rgb}{0,0.6,0} +\usetikzlibrary{arrows,intersections,math,calc} +\begin{document} +\def\skala{2} +\begin{tikzpicture}[>=latex,thick,scale=\skala] + +\def\pitick#1#2{ + \draw (#1,{-0.1/\skala}) -- (#1,{0.1/\skala}); + \node at (#1,{-0.1/\skala}) [below] {$\mathstrut #2$}; +} + +\def\achsen{ + \draw[->] ({-0.1/\skala},0) -- ({6.4+0.3/\skala},0) + coordinate[label=$u$]; + \pitick{3.1415}{\pi} + \pitick{2*3.1415}{2\pi} + \draw ({-0.1/\skala},1) -- ({0.1/\skala},1); + \node at ({-0.1/\skala},1) [left] {$1$}; + \node at ({-0.1/\skala},0) [left] {$0$}; +} + +\begin{scope} +\achsen +\draw[->] (0,{-1-0.1/\skala}) -- (0,{1+0.3/\skala}) + coordinate[label={right:$y$}]; +\draw[color=red!50,line width=1.4pt] + plot[domain=0:365,samples=100] ({3.1415*\x/180},{sin(\x)}); +\draw[color=red,line width=1.4pt] + plot[domain=0:6.4,samples=100] + ({\x},{(exp(\x)-exp(-\x))/(exp(\x)+exp(-\x))}); +\node[color=red] at ({1.5*3.1415},1) + [above] {$\operatorname{sn}(u,1)=\tanh u$}; +\node[color=red] at ({1.5*3.1415},{-1+0.14}) + [above] {$\operatorname{sn}(u,0)=\sin u$}; +\end{scope} + +\begin{scope}[yshift=-2.4cm] +\achsen +\draw[->] (0,{-1-0.1/\skala}) -- (0,{1+0.3/\skala}) + coordinate[label={right:$y$}]; + +\draw[color=blue!50,line width=1.4pt] + plot[domain=0:365,samples=100] ({3.1415*\x/180},{cos(\x)}); +\draw[color=blue,line width=1.4pt] + plot[domain=0:6.4,samples=100] ({\x},{2/(exp(\x)+exp(-\x))}); +\node[color=blue] at (3.1415,{-1+0.15}) + [above] {$\operatorname{cn}(u,0)=\cos u$}; +\node[color=blue] at (3.1415,{2/(exp(3.1415)+exp(-3.1415))+0.05}) + [above] {$\operatorname{sech}u$}; +\end{scope} + +\begin{scope}[yshift=-4.8cm] +\achsen +\draw[->] (0,{-0-0.1/\skala}) -- (0,{1+0.3/\skala}) + coordinate[label={right:$y$}]; +\node[color=darkgreen] at (3.1415,1) [above] {$\operatorname{dn}(u,0) = 1$}; +\draw[color=darkgreen,line width=1.4pt] + plot[domain=0:6.4,samples=100] ({\x},{2/(exp(\x)+exp(-\x))}); +\draw[color=darkgreen!50,line width=1.4pt] + (0,1) -- (6.4,1); +\node[color=darkgreen] at (3.1415,{2/(exp(3.1415)+exp(-3.1415))+0.05}) + [above] {$\operatorname{dn}(u,1)=\operatorname{sech}u$}; +\node at (0,{-0.1/\skala}) [below] {$\mathstrut 0$}; +\end{scope} + +\end{tikzpicture} +\end{document} + diff --git a/buch/chapters/110-elliptisch/jacobi.tex b/buch/chapters/110-elliptisch/jacobi.tex index e224490..f1e0987 100644 --- a/buch/chapters/110-elliptisch/jacobi.tex +++ b/buch/chapters/110-elliptisch/jacobi.tex @@ -489,6 +489,97 @@ Damit haben wir die grundlegenden Ableitungsregeln \end{align*} der elliptischen Funktionen nach Jacobi. +% +% Der Grenzfall $k=1$ +% +\subsubsection{Der Grenzwert $k\to1$} +\begin{figure} +\centering +\includegraphics{chapters/110-elliptisch/images/sncnlimit.pdf} +\caption{Grenzfälle der Jacobischen elliptischen Funktionen +für die Werte $0$ und $1$ des Parameters $k$. +\label{buch:elliptisch:fig:sncnlimit}} +\end{figure} +Für $k=1$ ist $k^{\prime2}=1-k^2=$ und es folgt aus den +Relationen~\eqref{buch:elliptisch:eqn:jacobi-relationen} +\[ +\operatorname{cn}^2(u,k) +- +k^2 +\operatorname{dn}^2(u,k) += +k^{\prime2} += +0 +\qquad\Rightarrow\qquad +\operatorname{cn}^2(u,1) += +\operatorname{dn}^2(u,1), +\] +die beiden Funktionen +$\operatorname{cn}(u,k)$ +und +$\operatorname{dn}(u,k)$ +fallen also zusammen. +Die Ableitungsregeln werden dadurch vereinfacht: +\begin{align*} +\operatorname{sn}'(u,1) +&= +\operatorname{cn}(u,1) +\operatorname{dn}(u,1) += +\operatorname{cn}^2(u,1) += +1-\operatorname{sn}^2(u,1) +&&\Rightarrow& y'&=1-y^2 +\\ +\operatorname{cn}'(u,1) +&= +- +\operatorname{sn}(u,1) +\operatorname{dn}(u,1) += +- +\operatorname{sn}(u,1)\operatorname{cn}(u,1) +&&\Rightarrow& +\frac{z'}{z}&=(\log z)' = -y +\end{align*} +Die erste Differentialgleichung für $y$ lässt sich separieren, man findet +die Lösung +\[ +\frac{y'}{1-y^2} += +1 +\quad\Rightarrow\quad +\int \frac{dy}{1-y^2} = \int \,du +\quad\Rightarrow\quad +\operatorname{artanh}(y) = u +\quad\Rightarrow\quad +\operatorname{sn}(u,1)=\tanh u. +\] +Damit kann man jetzt auch $z$ berechnen: +\begin{align*} +(\log \operatorname{cn}(u,1))' +&= +\tanh u +&&\Rightarrow& +\log\operatorname{cn}(u,1) +&= +-\int\tanh u\,du += +-\log\cosh u +\\ +& +&&\Rightarrow& +\operatorname{cn}(u,1) +&= +\frac{1}{\cosh u} += +\operatorname{sech}u. +\end{align*} +Die Grenzfunktionen sind in Abbildung~\ref{buch:elliptisch:fig:sncnlimit} +dargestellt. + % % Das Argument u % @@ -605,6 +696,7 @@ des Strahlensatzes geometrisch interpretieren lassen. \label{buch:elliptisch:fig:jacobi12}} \end{figure} \begin{table} +\centering \renewcommand{\arraystretch}{2.5} \begin{tabular}{|>{$\displaystyle}c<{$}|>{$\displaystyle}c<{$}>{$\displaystyle}c<{$}>{$\displaystyle}c<{$}>{$\displaystyle}c<{$}|} \hline @@ -616,34 +708,36 @@ des Strahlensatzes geometrisch interpretieren lassen. \\[5pt] \hline 1& -\operatorname{nn}(u,k)=\frac{1}{1} & +&%\operatorname{nn}(u,k)=\frac{1}{1} & \operatorname{ns}(u,k)=\frac{1}{\operatorname{sn}(u,k)} & \operatorname{nc}(u,k)=\frac{1}{\operatorname{cn}(u,k)} & \operatorname{nd}(u,k)=\frac{1}{\operatorname{dn}(u,k)} \\ \operatorname{sn}(u,k) & \operatorname{sn}(u,k)=\frac{\operatorname{sn}(u,k)}{1}& -\operatorname{ss}(u,k)=\frac{\operatorname{sn}(u,k)}{\operatorname{sn}(u,k)}& +&%\operatorname{ss}(u,k)=\frac{\operatorname{sn}(u,k)}{\operatorname{sn}(u,k)}& \operatorname{sc}(u,k)=\frac{\operatorname{sn}(u,k)}{\operatorname{cn}(u,k)}& \operatorname{sd}(u,k)=\frac{\operatorname{sn}(u,k)}{\operatorname{dn}(u,k)} \\ \operatorname{cn}(u,k) & \operatorname{cn}(u,k)=\frac{\operatorname{cn}(u,k)}{1} & \operatorname{cs}(u,k)=\frac{\operatorname{cn}(u,k)}{\operatorname{sn}(u,k)}& -\operatorname{cc}(u,k)=\frac{\operatorname{cn}(u,k)}{\operatorname{cn}(u,k)}& +&%\operatorname{cc}(u,k)=\frac{\operatorname{cn}(u,k)}{\operatorname{cn}(u,k)}& \operatorname{cd}(u,k)=\frac{\operatorname{cn}(u,k)}{\operatorname{dn}(u,k)} \\ \operatorname{dn}(u,k) & \operatorname{dn}(u,k)=\frac{\operatorname{dn}(u,k)}{1} & \operatorname{ds}(u,k)=\frac{\operatorname{dn}(u,k)}{\operatorname{sn}(u,k)}& \operatorname{dc}(u,k)=\frac{\operatorname{dn}(u,k)}{\operatorname{cn}(u,k)}& -\operatorname{dd}(u,k)=\frac{\operatorname{dn}(u,k)}{\operatorname{dn}(u,k)} +%\operatorname{dd}(u,k)=\frac{\operatorname{dn}(u,k)}{\operatorname{dn}(u,k)} \\[5pt] \hline \end{tabular} \caption{Zusammenstellung der abgeleiteten Jacobischen elliptischen -Funktionen als Quotienten der grundlegenden Jacobischen elliptischen -Funktionen. +Funktionen in hinteren drei Spalten als Quotienten der grundlegenden +Jacobischen elliptischen Funktionen. +Die erste Spalte zum Nenner $1$ enthält die grundlegenden +Jacobischen elliptischen Funktionen. \label{buch:elliptisch:table:abgeleitetjacobi}} \end{table} \subsubsection{Die abgeleiteten elliptischen Funktionen} @@ -675,8 +769,232 @@ Länge auf dem Strahl mit Polarwinkel $\vartheta$. Die Relationen~\ref{buch:elliptisch:eqn:jacobi-relationen} ermöglichen, jede Funktion $\operatorname{zn}(u,k)$ durch jede andere auszudrücken. +Die schiere Anzahl solcher Beziehungen macht es unmöglich, sie +übersichtlich in einer Tabelle zusammenzustellen, daher soll hier +nur an einem Beispiel das Vorgehen gezeigt werden: + +\begin{beispiel} +Die Funktion $\operatorname{sc}(u,k)$ soll durch $\operatorname{cd}(u,k)$ +ausgedrückt werden. +Zunächst ist +\[ +\operatorname{sc}(u,k) += +\frac{\operatorname{sn}(u,k)}{\operatorname{cn}(u,k)} +\] +nach Definition. +Im Resultat sollen nur noch $\operatorname{cn}(u,k)$ und +$\operatorname{dn}(u,k)$ vorkommen. +Daher eliminieren wir zunächst die Funktion $\operatorname{sn}(u,k)$ +mit Hilfe von \eqref{buch:elliptisch:eqn:jacobi-relationen} und erhalten +\begin{equation} +\operatorname{sc}(u,k) += +\frac{\sqrt{1-\operatorname{cn}^2(u,k)}}{\operatorname{cn}(u,k)}. +\label{buch:elliptisch:eqn:allgausdruecken} +\end{equation} +Nun genügt es, die Funktion $\operatorname{cn}(u,k)$ durch +$\operatorname{cd}(u,k)$ auszudrücken. +Aus der Definition und der +dritten Relation in \eqref{buch:elliptisch:eqn:jacobi-relationen} +erhält man +\begin{align*} +\operatorname{cd}^2(u,k) +&= +\frac{\operatorname{cn}^2(u,k)}{\operatorname{dn}^2(u,k)} += +\frac{\operatorname{cn}^2(u,k)}{k^{\prime2}+k^2\operatorname{cn}^2(u,k)} +\\ +\Rightarrow +\qquad +k^{\prime 2} +\operatorname{cd}^2(u,k) ++ +k^2\operatorname{cd}^2(u,k)\operatorname{cn}^2(u,k) +&= +\operatorname{cn}^2(u,k) +\\ +\operatorname{cn}^2(u,k) +- +k^2\operatorname{cd}^2(u,k)\operatorname{cn}^2(u,k) +&= +k^{\prime 2} +\operatorname{cd}^2(u,k) +\\ +\operatorname{cn}^2(u,k) +&= +\frac{ +k^{\prime 2} +\operatorname{cd}^2(u,k) +}{ +1 - k^2\operatorname{cd}^2(u,k) +} +\end{align*} +Für den Zähler brauchen wir $1-\operatorname{cn}^2(u,k)$, also +\[ +1-\operatorname{cn}^2(u,k) += +\frac{ +1 +- +k^2\operatorname{cd}^2(u,k) +- +k^{\prime 2} +\operatorname{cd}^2(u,k) +}{ +1 +- +k^2\operatorname{cd}^2(u,k) +} += +\frac{1-\operatorname{cd}^2(u,k)}{1-k^2\operatorname{cd}^2(u,k)} +\] +Einsetzen in~\eqref{buch:elliptisch:eqn:allgausdruecken} gibt +\begin{align*} +\operatorname{sc}(u,k) +&= +\frac{ +\sqrt{1-\operatorname{cd}^2(u,k)} +}{\sqrt{1-k^2\operatorname{cd}^2(u,k)}} +\cdot +\frac{ +\sqrt{1 - k^2\operatorname{cd}^2(u,k)} +}{ +k' +\operatorname{cd}(u,k) +} += +\frac{ +\sqrt{1-\operatorname{cd}^2(u,k)} +}{ +k' +\operatorname{cd}(u,k) +}. +\qedhere +\end{align*} +\end{beispiel} \subsubsection{Ableitung der abgeleiteten elliptischen Funktionen} +Aus den Ableitungen der grundlegenden Jacobischen elliptischen Funktionen +können mit der Quotientenregel nun auch beliebige Ableitungen der +abgeleiteten Jacobischen elliptischen Funktionen gefunden werden. +Als Beispiel berechnen wir die Ableitung von $\operatorname{sc}(u,k)$. +Sie ist +\begin{align*} +\frac{d}{du} +\operatorname{sc}(u,k) +&= +\frac{d}{du} +\frac{\operatorname{sn}(u,k)}{\operatorname{cn}(u,k)} += +\frac{ +\operatorname{sn}'(u,k)\operatorname{cn}(u,k) +- +\operatorname{sn}(u,k)\operatorname{cn}'(u,k)}{ +\operatorname{cn}^2(u,k) +} +\\ +&= +\frac{ +\operatorname{cn}^2(u,k)\operatorname{dn}(u,k) ++ +\operatorname{sn}^2(u,k)\operatorname{dn}(u,k) +}{ +\operatorname{cn}^2(u,k) +} += +\frac{( +\operatorname{sn}^2(u,k) ++ +\operatorname{cn}^2(u,k) +)\operatorname{dn}(u,k)}{ +\operatorname{cn}^2(u,k) +} +\\ +&= +\frac{1}{\operatorname{cn}(u,k)} +\cdot +\frac{\operatorname{dn}(u,k)}{\operatorname{cn}(u,k)} += +\operatorname{nc}(u,k) +\operatorname{dc}(u,k). +\end{align*} +Man beachte, dass das Quadrat der Nennerfunktion im Resultat +der Quotientenregel zur Folge hat, dass die +beiden Funktionen im Resultat beide den gleichen Nenner haben wie +die Funktion, die abgeleitet wird. + +Mit etwas Fleiss kann man nach diesem Muster alle Ableitungen +\begin{equation} +%\small +\begin{aligned} +\operatorname{sn}'(u,k) +&= +\phantom{-} +\operatorname{cn}(u,k)\,\operatorname{dn}(u,k) +&&\qquad& +\operatorname{ns}'(u,k) +&= +- +\operatorname{cs}(u,k)\,\operatorname{ds}(u,k) +\\ +\operatorname{cn}'(u,k) +&= +- +\operatorname{sn}(u,k)\,\operatorname{dn}(u,k) +&&& +\operatorname{nc}'(u,k) +&= +\phantom{-} +\operatorname{sc}(u,k)\,\operatorname{dc}(u,k) +\\ +\operatorname{dn}'(u,k) +&= +-k^2 +\operatorname{sn}(u,k)\,\operatorname{cn}(u,k) +&&& +\operatorname{nd}'(u,k) +&= +\phantom{-} +k^2 +\operatorname{sd}(u,k)\,\operatorname{cd}(u,k) +\\ +\operatorname{sc}'(u,k) +&= +\phantom{-} +\operatorname{dc}(u,k)\,\operatorname{nc}(u,k) +&&& +\operatorname{cs}'(u,k) +&= +- +\operatorname{ds}(u,k)\,\operatorname{ns}(u,k) +\\ +\operatorname{cd}'(u,k) +&= +-k^{\prime2} +\operatorname{sd}(u,k)\,\operatorname{nd}(u,k) +&&& +\operatorname{dc}'(u,k) +&= +\phantom{-} +k^{\prime2} +\operatorname{dc}(u,k)\,\operatorname{nc}(u,k) +\\ +\operatorname{ds}'(d,k) +&= +- +\operatorname{cs}(u,k)\,\operatorname{ns}(u,k) +&&& +\operatorname{sd}'(d,k) +&= +\phantom{-} +\operatorname{cd}(u,k)\,\operatorname{nd}(u,k) +\end{aligned} +\label{buch:elliptisch:eqn:alleableitungen} +\end{equation} +finden. +Man beachte, dass in jeder Identität alle Funktionen den gleichen +zweiten Buchstaben haben. \subsubsection{TODO} XXX algebraische Beziehungen \\ @@ -836,9 +1154,9 @@ Je nach Vorzeichen sind also eine andere elliptische Funktion als Lösung zu verwenden. % -% Jacobi elliptische Funktionen und elliptische Integrale +% Jacobischen elliptische Funktionen und elliptische Integrale % -\subsubsection{Jacobi elliptische Funktionen als elliptische Integrale} +\subsubsection{Jacobische elliptische Funktionen als elliptische Integrale} Die in Tabelle~\ref{buch:elliptisch:tabelle:loesungsfunktionen} zusammengestellten Differentialgleichungen ermöglichen nun, den Zusammenhang zwischen den Funktionen @@ -875,7 +1193,7 @@ ist daher \[ y(u) = F^{-1}(u+C). \] -Die Jacobi elliptischen Funktionen sind daher inverse Funktionen +Die Jacobischen elliptischen Funktionen sind daher inverse Funktionen der unvollständigen elliptischen Integrale. \subsubsection{Differentialgleichung zweiter Ordnung} -- cgit v1.2.1