aboutsummaryrefslogtreecommitdiffstats
path: root/buch/chapters/110-elliptisch
diff options
context:
space:
mode:
Diffstat (limited to 'buch/chapters/110-elliptisch')
-rw-r--r--buch/chapters/110-elliptisch/images/Makefile7
-rw-r--r--buch/chapters/110-elliptisch/images/jacobiplots.pdfbin57192 -> 57192 bytes
-rw-r--r--buch/chapters/110-elliptisch/images/sncnlimit.pdfbin0 -> 24327 bytes
-rw-r--r--buch/chapters/110-elliptisch/images/sncnlimit.tex79
-rw-r--r--buch/chapters/110-elliptisch/jacobi.tex336
5 files changed, 412 insertions, 10 deletions
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
--- a/buch/chapters/110-elliptisch/images/jacobiplots.pdf
+++ b/buch/chapters/110-elliptisch/images/jacobiplots.pdf
Binary files 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
--- /dev/null
+++ b/buch/chapters/110-elliptisch/images/sncnlimit.pdf
Binary files 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
@@ -490,6 +490,97 @@ Damit haben wir die grundlegenden Ableitungsregeln
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
%
\subsubsection{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}