diff options
Diffstat (limited to 'buch')
22 files changed, 1053 insertions, 30 deletions
diff --git a/buch/chapters/110-elliptisch/images/Makefile b/buch/chapters/110-elliptisch/images/Makefile index a7c9e74..c8f98cb 100644 --- a/buch/chapters/110-elliptisch/images/Makefile +++ b/buch/chapters/110-elliptisch/images/Makefile @@ -5,7 +5,7 @@ # all: lemniskate.pdf ellipsenumfang.pdf unvollstaendig.pdf rechteck.pdf \ ellipse.pdf pendel.pdf jacobiplots.pdf jacobidef.pdf jacobi12.pdf \ - sncnlimit.pdf slcl.pdf + sncnlimit.pdf slcl.pdf torusschnitt.pdf kegelpara.pdf lemniskate.pdf: lemniskate.tex pdflatex lemniskate.tex @@ -78,3 +78,27 @@ slcldata.tex: slcl ./slcl --outfile=slcldata.tex --a=0 --b=13.4 --steps=200 slcl.pdf: slcl.tex slcldata.tex pdflatex slcl.tex + +KEGELSIZE = -W256 -H256 +KEGELSIZE = -W128 -H128 +KEGELSIZE = -W1080 -H1080 +kegelpara.png: kegelpara.pov + povray +A0.1 $(KEGELSIZE) -Okegelpara.png kegelpara.pov + +kegelpara.jpg: kegelpara.png Makefile + convert -extract 1080x1040+0+0 kegelpara.png \ + -density 300 -units PixelsPerInch kegelpara.jpg + +kegelpara.pdf: kegelpara.tex kegelpara.jpg + pdflatex kegelpara.tex + +torusschnitt.png: torusschnitt.pov + povray +A0.1 -W1920 -H1080 -Otorusschnitt.png torusschnitt.pov + +torusschnitt.jpg: torusschnitt.png Makefile + convert -extract 1560x1080+180+0 torusschnitt.png \ + -density 300 -units PixelsPerInch torusschnitt.jpg + +torusschnitt.pdf: torusschnitt.tex torusschnitt.jpg + pdflatex torusschnitt.tex + diff --git a/buch/chapters/110-elliptisch/images/jacobiplots.pdf b/buch/chapters/110-elliptisch/images/jacobiplots.pdf Binary files differindex f0e6e78..c11affc 100644 --- a/buch/chapters/110-elliptisch/images/jacobiplots.pdf +++ b/buch/chapters/110-elliptisch/images/jacobiplots.pdf diff --git a/buch/chapters/110-elliptisch/images/kegelpara.pdf b/buch/chapters/110-elliptisch/images/kegelpara.pdf Binary files differnew file mode 100644 index 0000000..2f76593 --- /dev/null +++ b/buch/chapters/110-elliptisch/images/kegelpara.pdf diff --git a/buch/chapters/110-elliptisch/images/kegelpara.pov b/buch/chapters/110-elliptisch/images/kegelpara.pov new file mode 100644 index 0000000..13b66cc --- /dev/null +++ b/buch/chapters/110-elliptisch/images/kegelpara.pov @@ -0,0 +1,329 @@ +// +// kegelpara.pov +// +// (c) 2022 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +// +#version 3.7; +#include "colors.inc" + +#declare O = <0,0,0>; + +global_settings { + assumed_gamma 1 +} + +#declare imagescale = 0.08; + +camera { + location <28, 20, -40> + look_at <0, 0.1, 0> + right x * imagescale + up y * imagescale +} + +light_source { + <30, 10, -40> color White + area_light <1,0,0> <0,0,1>, 10, 10 + adaptive 1 + jitter +} + +sky_sphere { + pigment { + color rgb<1,1,1> + } +} + + +// +// draw an arrow from <from> to <to> with thickness <arrowthickness> with +// color <c> +// +#macro arrow(from, to, arrowthickness, c) +#declare arrowdirection = vnormalize(to - from); +#declare arrowlength = vlength(to - from); +union { + sphere { + from, 1.1 * arrowthickness + } + cylinder { + from, + from + (arrowlength - 5 * arrowthickness) * arrowdirection, + arrowthickness + } + cone { + from + (arrowlength - 5 * arrowthickness) * arrowdirection, + 2 * arrowthickness, + to, + 0 + } + pigment { + color c + } + finish { + specular 0.9 + metallic + } +} +#end + +arrow(<-2.6,0,0>,<2.5,0,0>,0.02,White) +arrow(<0,-2,0>,<0,2.3,0>,0.02,White) +arrow(<0,0,-3.2>,<0,0,3.7>,0.02,White) + +#declare epsilon = 0.0001; +#declare l = 1.5; + +#macro Kegel(farbe) +union { + difference { + cone { O, 0, <l, 0, 0>, l } + cone { O + <epsilon, 0,0>, 0, <l+epsilon, 0, 0>, l } + } + difference { + cone { O, 0, <-l, 0, 0>, l } + cone { O + <-epsilon, 0, 0>, 0, <-l-epsilon, 0, 0>, l } + } + pigment { + color farbe + } + finish { + specular 0.9 + metallic + } +} +#end + +#macro Kegelpunkt(xx, phi) + < xx, xx * sin(phi), xx * cos(phi) > +#end + +#macro Kegelgitter(farbe, r) +union { + #declare s = 0; + #declare smax = 2 * pi; + #declare sstep = pi / 6; + #while (s < smax - sstep/2) + cylinder { Kegelpunkt(l, s), Kegelpunkt(-l, s), r } + #declare s = s + sstep; + #end + #declare phimax = 2 * pi; + #declare phisteps = 100; + #declare phistep = phimax / phisteps; + #declare xxstep = 0.5; + #declare xxmax = 2; + #declare xx = xxstep; + #while (xx < xxmax - xxstep/2) + #declare phi = 0; + #while (phi < phimax - phistep/2) + cylinder { + Kegelpunkt(xx, phi), + Kegelpunkt(xx, phi + phistep), + r + } + sphere { Kegelpunkt(xx, phi), r } + cylinder { + Kegelpunkt(-xx, phi), + Kegelpunkt(-xx, phi + phistep), + r + } + sphere { Kegelpunkt(-xx, phi), r } + #declare phi = phi + phistep; + #end + #declare xx = xx + xxstep; + #end + pigment { + color farbe + } + finish { + specular 0.9 + metallic + } +} +#end + +#macro F(w, r) + <r * cos(w), r * r/sqrt(2), r * sin(w) > +#end + +#macro Paraboloid(farbe) +mesh { + #declare phi = 0; + #declare phimax = 2 * pi; + #declare phisteps = 100; + #declare phistep = pi / phisteps; + #declare rsteps = 100; + #declare rmax = 1.5; + #declare rstep = rmax / rsteps; + #while (phi < phimax - phistep/2) + #declare r = rstep; + #declare h = r * r / sqrt(2); + triangle { + O, F(phi, r), F(phi + phistep, r) + } + #while (r < rmax - rstep/2) + // ring + triangle { + F(phi, r), + F(phi + phistep, r), + F(phi + phistep, r + rstep) + } + triangle { + F(phi, r), + F(phi + phistep, r + rstep), + F(phi, r + rstep) + } + #declare r = r + rstep; + #end + #declare phi = phi + phistep; + #end + pigment { + color farbe + } + finish { + specular 0.9 + metallic + } +} +#end + +#macro Paraboloidgitter(farbe, gr) +union { + #declare phi = 0; + #declare phimax = 2 * pi; + #declare phistep = pi / 6; + + #declare rmax = 1.5; + #declare rsteps = 100; + #declare rstep = rmax / rsteps; + + #while (phi < phimax - phistep/2) + #declare r = rstep; + #while (r < rmax - rstep/2) + cylinder { F(phi, r), F(phi, r + rstep), gr } + sphere { F(phi, r), gr } + #declare r = r + rstep; + #end + #declare phi = phi + phistep; + #end + + #declare rstep = 0.2; + #declare r = rstep; + + #declare phisteps = 100; + #declare phistep = phimax / phisteps; + #while (r < rmax) + #declare phi = 0; + #while (phi < phimax - phistep/2) + cylinder { F(phi, r), F(phi + phistep, r), gr } + sphere { F(phi, r), gr } + #declare phi = phi + phistep; + #end + #declare r = r + rstep; + #end + pigment { + color farbe + } + finish { + specular 0.9 + metallic + } +} +#end + +#declare a = sqrt(2); +#macro G(phi,sg) + < a*sg*sqrt(cos(2*phi))*cos(phi), a*cos(2*phi), a*sqrt(cos(2*phi))*sin(phi)> +#end + +#macro Lemniskate3D(s, farbe) +union { + #declare phi = -pi / 4; + #declare phimax = pi / 4; + #declare phisteps = 100; + #declare phistep = phimax / phisteps; + #while (phi < phimax - phistep/2) + sphere { G(phi,1), s } + cylinder { G(phi,1), G(phi+phistep,1), s } + sphere { G(phi,-1), s } + cylinder { G(phi,-1), G(phi+phistep,-1), s } + #declare phi = phi + phistep; + #end + pigment { + color farbe + } + finish { + specular 0.9 + metallic + } +} +#end + +#declare a = sqrt(2); +#macro G2(phi,sg) + a * sqrt(cos(2*phi)) * < sg * cos(phi), 0, sin(phi)> +#end + +#macro Lemniskate(s, farbe) +union { + #declare phi = -pi / 4; + #declare phimax = pi / 4; + #declare phisteps = 100; + #declare phistep = phimax / phisteps; + #while (phi < phimax - phistep/2) + sphere { G2(phi,1), s } + cylinder { G2(phi,1), G2(phi+phistep,1), s } + sphere { G2(phi,-1), s } + cylinder { G2(phi,-1), G2(phi+phistep,-1), s } + #declare phi = phi + phistep; + #end + pigment { + color farbe + } + finish { + specular 0.9 + metallic + } +} +#end + +#macro Projektion(s, farbe) +union { + #declare phistep = pi / 16; + #declare phi = -pi / 4 + phistep; + #declare phimax = pi / 4; + #while (phi < phimax - phistep/2) + cylinder { G(phi, 1), G2(phi, 1), s } + cylinder { G(phi, -1), G2(phi, -1), s } + #declare phi = phi + phistep; + #end + pigment { + color farbe + } + finish { + specular 0.9 + metallic + } +} +#end + +#declare kegelfarbe = rgbf<0.2,0.6,0.2,0.2>; +#declare kegelgitterfarbe = rgb<0.2,0.8,0.2>; +#declare paraboloidfarbe = rgbf<0.2,0.6,1.0,0.2>; +#declare paraboloidgitterfarbe = rgb<0.4,1,1>; + +//intersection { +// union { + Paraboloid(paraboloidfarbe) + Paraboloidgitter(paraboloidgitterfarbe, 0.004) + + Kegel(kegelfarbe) + Kegelgitter(kegelgitterfarbe, 0.004) +// } +// plane { <0, 0, -1>, 0.6 } +//} + + +Lemniskate3D(0.02, rgb<0.8,0.0,0.8>) +Lemniskate(0.02, Red) +Projektion(0.01, Yellow) diff --git a/buch/chapters/110-elliptisch/images/kegelpara.tex b/buch/chapters/110-elliptisch/images/kegelpara.tex new file mode 100644 index 0000000..8fcefbf --- /dev/null +++ b/buch/chapters/110-elliptisch/images/kegelpara.tex @@ -0,0 +1,41 @@ +% +% kegelpara.tex +% +% (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +% +\documentclass[tikz]{standalone} +\usepackage{times} +\usepackage{amsmath} +\usepackage{txfonts} +\usepackage[utf8]{inputenc} +\usepackage{graphics} +\usetikzlibrary{arrows,intersections,math} +\usepackage{ifthen} +\begin{document} + +\newboolean{showgrid} +\setboolean{showgrid}{false} +\def\breite{4} +\def\hoehe{4} + +\begin{tikzpicture}[>=latex,thick] + +% Povray Bild +\node at (0,0) {\includegraphics[width=8cm]{kegelpara.jpg}}; + +% Gitter +\ifthenelse{\boolean{showgrid}}{ +\draw[step=0.1,line width=0.1pt] (-\breite,-\hoehe) grid (\breite, \hoehe); +\draw[step=0.5,line width=0.4pt] (-\breite,-\hoehe) grid (\breite, \hoehe); +\draw (-\breite,-\hoehe) grid (\breite, \hoehe); +\fill (0,0) circle[radius=0.05]; +}{} + +\node at (4.1,-1.4) {$X$}; +\node at (0.2,3.8) {$Z$}; +\node at (4.0,1.8) {$Y$}; + +\end{tikzpicture} + +\end{document} + diff --git a/buch/chapters/110-elliptisch/images/torusschnitt.pdf b/buch/chapters/110-elliptisch/images/torusschnitt.pdf Binary files differnew file mode 100644 index 0000000..11bd353 --- /dev/null +++ b/buch/chapters/110-elliptisch/images/torusschnitt.pdf diff --git a/buch/chapters/110-elliptisch/images/torusschnitt.pov b/buch/chapters/110-elliptisch/images/torusschnitt.pov new file mode 100644 index 0000000..43d50c6 --- /dev/null +++ b/buch/chapters/110-elliptisch/images/torusschnitt.pov @@ -0,0 +1,265 @@ +// +// kegelpara.pov +// +// (c) 2022 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +// +#version 3.7; +#include "colors.inc" + +#declare O = <0,0,0>; + +global_settings { + assumed_gamma 1 +} + +#declare imagescale = 0.060; + +camera { + location <28, 20, -40> + look_at <0, 0.55, 0> + right (16/9) * x * imagescale + up y * imagescale +} + +light_source { + <30, 10, -40> color White + area_light <1,0,0> <0,0,1>, 10, 10 + adaptive 1 + jitter +} + +sky_sphere { + pigment { + color rgb<1,1,1> + } +} + + +// +// draw an arrow from <from> to <to> with thickness <arrowthickness> with +// color <c> +// +#macro arrow(from, to, arrowthickness, c) +#declare arrowdirection = vnormalize(to - from); +#declare arrowlength = vlength(to - from); +union { + sphere { + from, 1.1 * arrowthickness + } + cylinder { + from, + from + (arrowlength - 5 * arrowthickness) * arrowdirection, + arrowthickness + } + cone { + from + (arrowlength - 5 * arrowthickness) * arrowdirection, + 2 * arrowthickness, + to, + 0 + } + pigment { + color c + } + finish { + specular 0.9 + metallic + } +} +#end + +arrow(<-2,0,0>,<2,0,0>,0.02,White) +arrow(<0,-1.1,0>,<0,2.2,0>,0.02,White) +arrow(<0,0,-1.6>,<0,0,2.4>,0.02,White) + +#declare epsilon = 0.001; +#declare l = 1.5; + + +#declare a = sqrt(2); +#macro G2(phi,sg) + a * sqrt(cos(2*phi)) * < sg * cos(phi), 0, sin(phi)> +#end + +#macro Lemniskate(s, farbe) +union { + #declare phi = -pi / 4; + #declare phimax = pi / 4; + #declare phisteps = 100; + #declare phistep = phimax / phisteps; + #while (phi < phimax - phistep/2) + sphere { G2(phi,1), s } + cylinder { G2(phi,1), G2(phi+phistep,1), s } + sphere { G2(phi,-1), s } + cylinder { G2(phi,-1), G2(phi+phistep,-1), s } + #declare phi = phi + phistep; + #end + pigment { + color farbe + } + finish { + specular 0.9 + metallic + } +} +#end + +#macro Projektion(s, farbe) +union { + #declare phistep = pi / 16; + #declare phi = -pi / 4 + phistep; + #declare phimax = pi / 4; + #while (phi < phimax - phistep/2) + cylinder { G(phi, 1), G2(phi, 1), s } + cylinder { G(phi, -1), G2(phi, -1), s } + #declare phi = phi + phistep; + #end + pigment { + color farbe + } + finish { + specular 0.9 + metallic + } +} +#end + +#macro Ebene(l, b, farbe) +mesh { + triangle { <-l, 0, -b>, < l, 0, -b>, < l, 0, b> } + triangle { <-l, 0, -b>, < l, 0, b>, <-l, 0, b> } + pigment { + color farbe + } + finish { + specular 0.9 + metallic + } +} +#end + +#macro Ebenengitter(l, b, s, r, farbe) +union { + #declare lmax = floor(l / s); + #declare ll = -lmax; + #while (ll <= lmax) + cylinder { <ll * s, 0, -b>, <ll * s, 0, b>, r } + #declare ll = ll + 1; + #end + #declare bmax = floor(b / s); + #declare bb = -bmax; + #while (bb <= bmax) + cylinder { <-l, 0, bb * s>, <l, 0, bb * s>, r } + #declare bb = bb + 1; + #end + pigment { + color farbe + } + finish { + specular 0.9 + metallic + } +} +#end + +#declare b = 0.5; +#macro T(phi, theta) + b * < (2 + cos(theta)) * cos(phi), (2 + cos(theta)) * sin(phi) + 1, sin(theta) > +#end + +#macro breitenkreis(theta, r) + #declare phi = 0; + #declare phimax = 2 * pi; + #declare phisteps = 200; + #declare phistep = phimax / phisteps; + #while (phi < phimax - phistep/2) + cylinder { T(phi, theta), T(phi + phistep, theta), r } + sphere { T(phi, theta), r } + #declare phi = phi + phistep; + #end +#end + +#macro laengenkreis(phi, r) + #declare theta = 0; + #declare thetamax = 2 * pi; + #declare thetasteps = 200; + #declare thetastep = thetamax / thetasteps; + #while (theta < thetamax - thetastep/2) + cylinder { T(phi, theta), T(phi, theta + thetastep), r } + sphere { T(phi, theta), r } + #declare theta = theta + thetastep; + #end +#end + +#macro Torusgitter(farbe, r) +union { + #declare phi = 0; + #declare phimax = 2 * pi; + #declare phistep = pi / 6; + #while (phi < phimax - phistep/2) + laengenkreis(phi, r) + #declare phi = phi + phistep; + #end + #declare thetamax = pi; + #declare thetastep = pi / 6; + #declare theta = thetastep; + #while (theta < thetamax - thetastep/2) + breitenkreis(theta, r) + breitenkreis(thetamax + theta, r) + #declare theta = theta + thetastep; + #end + breitenkreis(0, 1.5 * r) + breitenkreis(pi, 1.5 * r) + pigment { + color farbe + } + finish { + specular 0.9 + metallic + } +} +#end + +#macro Torus(farbe) +mesh { + #declare phi = 0; + #declare phimax = 2 * pi; + #declare phisteps = 200; + #declare phistep = phimax/phisteps; + #while (phi < phimax - phistep/2) + #declare theta = 0; + #declare thetamax = 2 * pi; + #declare thetasteps = 200; + #declare thetastep = thetamax / thetasteps; + #while (theta < thetamax - thetastep/2) + triangle { + T(phi, theta), + T(phi + phistep, theta), + T(phi + phistep, theta + thetastep) + } + triangle { + T(phi, theta), + T(phi + phistep, theta + thetastep), + T(phi, theta + thetastep) + } + #declare theta = theta + thetastep; + #end + #declare phi = phi + phistep; + #end + pigment { + color farbe + } + finish { + specular 0.9 + metallic + } +} +#end + +#declare torusfarbe = rgbt<0.2,0.6,0.2,0.2>; +#declare ebenenfarbe = rgbt<0.2,0.6,1.0,0.2>; + +Lemniskate(0.02, Red) +Ebene(1.8, 1.4, ebenenfarbe) +Ebenengitter(1.8, 1.4, 0.5, 0.005, rgb<0.4,1,1>) +Torus(torusfarbe) +Torusgitter(Yellow, 0.005) diff --git a/buch/chapters/110-elliptisch/images/torusschnitt.tex b/buch/chapters/110-elliptisch/images/torusschnitt.tex new file mode 100644 index 0000000..3053ac5 --- /dev/null +++ b/buch/chapters/110-elliptisch/images/torusschnitt.tex @@ -0,0 +1,41 @@ +% +% torusschnitt.tex +% +% (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +% +\documentclass[tikz]{standalone} +\usepackage{times} +\usepackage{amsmath} +\usepackage{txfonts} +\usepackage[utf8]{inputenc} +\usepackage{graphics} +\usetikzlibrary{arrows,intersections,math} +\usepackage{ifthen} +\begin{document} + +\newboolean{showgrid} +\setboolean{showgrid}{false} +\def\breite{6} +\def\hoehe{4} + +\begin{tikzpicture}[>=latex,thick] + +% Povray Bild +\node at (0,0) {\includegraphics[width=11.4cm]{torusschnitt.jpg}}; + +% Gitter +\ifthenelse{\boolean{showgrid}}{ +\draw[step=0.1,line width=0.1pt] (-\breite,-\hoehe) grid (\breite, \hoehe); +\draw[step=0.5,line width=0.4pt] (-\breite,-\hoehe) grid (\breite, \hoehe); +\draw (-\breite,-\hoehe) grid (\breite, \hoehe); +\fill (0,0) circle[radius=0.05]; +}{} + +\node at (4.4,-2.4) {$X$}; +\node at (3.5,0.6) {$Y$}; +\node at (0.3,3.8) {$Z$}; + +\end{tikzpicture} + +\end{document} + diff --git a/buch/chapters/110-elliptisch/lemniskate.tex b/buch/chapters/110-elliptisch/lemniskate.tex index f750a82..fceaadf 100644 --- a/buch/chapters/110-elliptisch/lemniskate.tex +++ b/buch/chapters/110-elliptisch/lemniskate.tex @@ -12,6 +12,9 @@ veröffentlich hat. In diesem Abschnitt soll die Verbindung zu den Jacobischen elliptischen Funktionen hergestellt werden. +% +% Lemniskate +% \subsection{Lemniskate \label{buch:gemotrie:subsection:lemniskate}} \begin{figure} @@ -71,6 +74,165 @@ Sie gilt für Winkel $\varphi\in[-\frac{\pi}4,\frac{\pi}4]$ für das rechte Blatt und $\varphi\in[\frac{3\pi}4,\frac{5\pi}4]$ für das linke Blatt der Lemniskate. +% +% Schnitt eines Kegels mit einem Paraboloid +% +\subsubsection{Schnitt eines Kegels mit einem Paraboloid} +\begin{figure} +\center +\includegraphics{chapters/110-elliptisch/images/kegelpara.pdf} +\caption{Leminiskate (rot) als Projektion (gelb) der Schnittkurve (pink) +eines geraden +Kreiskegels (grün) mit einem Rotationsparaboloid (hellblau). +\label{buch:elliptisch:lemniskate:kegelpara}} +\end{figure}% +Schreibt man in der Gleichung~\eqref{buch:elliptisch:eqn:lemniskate} +für die Klammer auf der rechten Seite $Z^2 = X^2 - Y^2$, dann wird die +Lemniskate die Projektion in die $X$-$Y$-Ebene der Schnittmenge der Flächen, +die durch die Gleichungen +\begin{equation} +X^2-Y^2 = Z^2 +\qquad\text{und}\qquad +(X^2+Y^2) = R^2 = \sqrt{2}aZ +\label{buch:elliptisch:eqn:kegelparabolschnitt} +\end{equation} +beschrieben wird. +Die linke Gleichung in +\eqref{buch:elliptisch:eqn:kegelparabolschnitt} +beschreibt einen geraden Kreiskegel, die rechte ist ein Rotationsparaboloid. +Die Schnittkurve ist in Abbildung~\ref{buch:elliptisch:lemniskate:kegelpara} +dargestellt. + +\subsubsection{Schnitt eines Torus mit einer Ebene} +\begin{figure} +\centering +\includegraphics{chapters/110-elliptisch/images/torusschnitt.pdf} +\caption{Die Schnittkurve (rot) eines Torus (grün) +mit einer zur Torusachse parallelen Ebene (blau), +die den inneren Äquator des Torus berührt, ist eine Lemniskate. +\label{buch:elliptisch:lemniskate:torusschnitt}} +\end{figure} +Schneidet man einen Torus mit einer Ebene, die zur Achse des Torus +parallel ist und den inneren Äquator des Torus berührt, entsteht +ebenfalls eine Lemniskate. +Die Situation ist in Abbildung~\ref{buch:elliptisch:lemniskate:torusschnitt} +dargestellt. + +Der Torus kann mit den Radien $2$ und $1$ mit der $y$-Achse als Torusachse +kann mit der Parametrisierung +\[ +(s,t) +\mapsto +\begin{pmatrix} +(2+\cos s) \cos t \\ +\sin s \\ +(2+\cos s) \sin t + 1 +\end{pmatrix} +\] +beschrieben werden. +Die Gleichung $z=1$ beschreibt eine +achsparallele Ebene, die den inneren Äquator berührt. +Die Schnittkurve erfüllt daher +\[ +(2+\cos s)\sin t + 1 = 0, +\] +was wir auch als $2 +\cos s = -1/\sin t$ schreiben können. +Wir müssen nachprüfen dass die Koordinaten +$X=(2+\cos s)\cos t$ und $Y=\sin s$ die Gleichung einer Lemniskate +erfüllen. + +Zunächst können wir in der $X$-Koordinate den Klammerausdruck durch +\begin{equation} +X += +(2+\cos s) \cos t += +-\frac{1}{\sin t}\cos t += +-\frac{\cos t}{\sin t} +\qquad\Rightarrow\qquad +X^2 += +\frac{\cos^2t}{\sin^2 t} += +\frac{1-\sin^2t}{\sin^2 t} +\label{buch:elliptisch:lemniskate:Xsin} +\end{equation} +ersetzen. +Auch die $Y$-Koordinaten können wir durch $t$ ausdrücken, +nämlich +\begin{equation} +Y^2=\sin^2 s = 1-\cos^2 s += +1- +\biggl( +\frac{1}{\sin t} +-2 +\biggr)^2 += +\frac{-3\sin^2 t+4\sin t-1}{\sin^2 t}. +\label{buch:elliptisch:lemniskate:Ysin} +\end{equation} +Die Gleichungen +\eqref{buch:elliptisch:lemniskate:Xsin} +und +\eqref{buch:elliptisch:lemniskate:Ysin} +zeigen, dass man $X^2$ und $Y^2$ sogar einzig durch $\sin t$ +parametrisieren kann. +Um die Ausdrücke etwas zu vereinfachen, schreiben wir $S=\sin t$ +und erhalten zusammenfassend +\begin{equation} +\begin{aligned} +X^2 +&= +\frac{1-S^2}{S^2} +\\ +Y^2 +&= +\frac{-3S^2+4S-1}{S^2}. +\end{aligned} +\end{equation} +Daraus kann man jetzt die Summen und Differenzen der Quadrate +berechnen, sie sind +\begin{equation} +\begin{aligned} +X^2+Y^2 +&= +\frac{-4S^2+4S}{S^2} += +\frac{4S(1-S)}{S^2} += +\frac{4(1-S)}{S} += +4\frac{1-S}{S} +\\ +X^2-Y^2 +&= +\frac{2-4S+2S^2}{S^2} += +\frac{2(1-S)^2}{S^2} += +2\biggl(\frac{1-S}{S}\biggr)^2. +\end{aligned} +\end{equation} +Die Berechnung des Quadrates von $X^2+Y^2$ ergibt +die Gleichung +\[ +(X^2+Y^2) += +16 +\biggl(\frac{1-S}{S}\biggr)^2 += +8 \cdot 2 +\biggl(\frac{1-S}{S}\biggr)^2 += +2\cdot 2^2\cdot (X-Y)^2. +\] +Sie ist eine Lemniskaten-Gleichung für $a=2$. + +% +% Bogenlänge der Lemniskate +% \subsection{Bogenlänge} Die Funktionen \begin{equation} diff --git a/buch/papers/nav/beispiel.txt b/buch/papers/nav/beispiel.txt index 12c309a..b8716fc 100644 --- a/buch/papers/nav/beispiel.txt +++ b/buch/papers/nav/beispiel.txt @@ -5,31 +5,15 @@ Sternzeit: 7h 54m 26.593s 7.90738694h Deneb RA 20h 42m 12.14s 20.703372h -DEC 45g 21' 40.3" 45.361194 +DEC 45 21' 40.3" 45.361194 H 50g 15' 17.1" 50.254750 Azi 59g 36' 02.0" 59.600555 -Altair - -RA 19h 51' 53.39" 19.864831h -DEC 8g 55' 42.3 8.928416 - -H 45g 27' 48.1" 45.463361 -Azi 117g 16' 14.1" 117.270583 - -Arktur - -RA 14h 16' 42.14" 14.278372 -DEC 19g 03' 47.6 19.063222 - -H 47g 25' 38.8" 47.427444 -Azi 259g 09' 38.4" 259.160666 - Spica RA 13h 26m 23.44s 13.439844h -DEC -11g 16' 46.8" -11.279666 +DEC -11g 16' 46.8" 11.279666 H 18g 27' 30.0" 18.458333 Azi 240g 23' 52.5" 240.397916 diff --git a/buch/papers/nav/bilder/beispiele1.pdf b/buch/papers/nav/bilder/beispiele1.pdf Binary files differnew file mode 100644 index 0000000..d0fe3dc --- /dev/null +++ b/buch/papers/nav/bilder/beispiele1.pdf diff --git a/buch/papers/nav/bilder/beispiele2.pdf b/buch/papers/nav/bilder/beispiele2.pdf Binary files differnew file mode 100644 index 0000000..8579ee5 --- /dev/null +++ b/buch/papers/nav/bilder/beispiele2.pdf diff --git a/buch/papers/nav/bilder/position1.pdf b/buch/papers/nav/bilder/position1.pdf Binary files differnew file mode 100644 index 0000000..ba7755f --- /dev/null +++ b/buch/papers/nav/bilder/position1.pdf diff --git a/buch/papers/nav/bilder/position2.pdf b/buch/papers/nav/bilder/position2.pdf Binary files differnew file mode 100644 index 0000000..3333dd4 --- /dev/null +++ b/buch/papers/nav/bilder/position2.pdf diff --git a/buch/papers/nav/bilder/position3.pdf b/buch/papers/nav/bilder/position3.pdf Binary files differnew file mode 100644 index 0000000..fae0b85 --- /dev/null +++ b/buch/papers/nav/bilder/position3.pdf diff --git a/buch/papers/nav/bilder/position4.pdf b/buch/papers/nav/bilder/position4.pdf Binary files differnew file mode 100644 index 0000000..ac80c46 --- /dev/null +++ b/buch/papers/nav/bilder/position4.pdf diff --git a/buch/papers/nav/bilder/position5.pdf b/buch/papers/nav/bilder/position5.pdf Binary files differnew file mode 100644 index 0000000..afe120e --- /dev/null +++ b/buch/papers/nav/bilder/position5.pdf diff --git a/buch/papers/nav/bsp.tex b/buch/papers/nav/bsp.tex new file mode 100644 index 0000000..d544588 --- /dev/null +++ b/buch/papers/nav/bsp.tex @@ -0,0 +1,181 @@ +\section{Beispielrechnung} + +\subsection{Einführung} +In diesem Abschnitt wird die Theorie vom Abschnitt 21.6 in einem Praxisbeispiel angewendet. +Wir haben die Deklination, Rektaszension, Höhe der beiden Planeten Deneb und Arktur und die Sternzeit von Greenwich als Ausgangslage. +Die Deklinationen und Rektaszensionen sind von einem vergangenen Datum und die Höhe der Gestirne und die Sternzeit wurden von unserem Dozenten digital in einer Stadt in Japan mit den Koordinaten 35.716672 N, 140.233336 E bestimmt. +Wir werden rechnerisch beweisen, dass wir mit diesen Ergebnissen genau auf diese Koordinaten kommen. +\subsection{Vorgehen} + +\begin{center} + \begin{tabular}{l l l} + 1. & Koordinaten der Bildpunkte der Gestirne bestimmen \\ + 2. & Dreiecke aufzeichnen und richtig beschriften\\ + 3. & Dreieck ABC bestimmmen\\ + 4. & Dreieck BPC bestimmen \\ + 5. & Dreieck ABP bestimmen \\ + 6. & Geographische Breite bestimmen \\ + 7. & Geographische Länge bestimmen \\ + \end{tabular} +\end{center} + +\subsection{Ausgangslage} +\begin{wrapfigure}{R}{5.6cm} + \includegraphics{papers/nav/bilder/position1.pdf} + \caption{Ausgangslage} +\end{wrapfigure} +Die Rektaszension und die Sternzeit sind in der Regeln in Stunden angegeben. +Für die Umrechnung in Grad kann folgender Zusammenhang verwendet werden: +\[ Stunden \cdot 15 = Grad\]. +Dies wurde hier bereits gemacht. +\begin{center} + \begin{tabular}{l l l} + Sternzeit $s$ & $118.610804^\circ$ \\ + Deneb&\\ + & Rektaszension $RA_{Deneb}$& $310.55058^\circ$ \\ + & Deklination $DEC_{Deneb}$& $45.361194^\circ$ \\ + & Höhe $h_c$ & $50.256027^\circ$ \\ + Arktur &\\ + & Rektaszension $RA_{Arktur}$& $214.17558^\circ$ \\ + & Deklination $DEC_{Arktur}$& $19.063222^\circ$ \\ + & Höhe $h_b$ & $47.427444^\circ$ \\ + \end{tabular} +\end{center} +\subsection{Koordinaten der Bildpunkte} +Als erstes benötigen wir die Koordinaten der Bildpunkte von Arktur und Deneb. +$\delta$ ist die Breite, $\lambda$ die Länge. +\begin{align} +\delta_{Deneb}&=DEC_{Deneb} = \underline{\underline{45.361194^\circ}} \nonumber \\ +\lambda_{Deneb}&=RA_{Deneb} - s = 310.55058^\circ -118.610804^\circ =\underline{\underline{191.939776^\circ}} \nonumber \\ +\delta_{Arktur}&=DEC_{Arktur} = \underline{\underline{19.063222^\circ}} \nonumber \\ +\lambda_{Arktur}&=RA_{Arktur} - s = 214.17558^\circ -118.610804^\circ = \underline{\underline{5.5647759^\circ}} \nonumber +\end{align} + + +\subsection{Dreiecke definieren} +\begin{figure} + \begin{center} + \includegraphics[width=6cm]{papers/nav/bilder/beispiele1.pdf} + \includegraphics[width=6cm]{papers/nav/bilder/beispiele2.pdf} + \caption{Arktur-Deneb; Spica-Altiar} +\end{center} +\end{figure} +Das Festlegen der Dreiecke ist essenziell für die korrekten Berechnungen. +Ein Problem, welches in der Theorie nicht berücksichtigt wurde ist, dass der Punkt $P$ nicht zwingend unterhalb der Seite $a$ sein muss. +Wenn man das nicht berücksichtigt, erhält man falsche oder keine Ergebnisse. +In der Realität weiss man jedoch ungefähr auf welchem Breitengrad man ist, so kann man relativ einfach entscheiden, ob der eigene Standort über $a$ ist oder nicht. +Beim unserem genutzten Paar Arktur-Deneb ist dies kein Problem, da der Punkt unterhalb der Seite $a$ liegt. +Würde man aber das Paar Altair-Spica nehmen, liegt $P$ über $a$ (vgl. Abbildung 21.11) und man müsste trigonometrisch anders vorgehen. + +\subsection{Dreieck $ABC$} +\begin{wrapfigure}{R}{5.6cm} + \includegraphics{papers/nav/bilder/position2.pdf} + \caption{Dreieck ABC} +\end{wrapfigure} +Nun berechnen wir alle Seitenlängen $a$, $b$, $c$ und die Innnenwinkel $\alpha$, $\beta$ und $\gamma$ +Wir können $b$ und $c$ mit den geltenten Zusammenhängen des nautischen Dreiecks wie folgt bestimmen: +\begin{align} + b=90^\circ-DEC_{Deneb} = 90^\circ - 45.361194^\circ = \underline{\underline{44.638806^\circ}}\nonumber \\ + c=90^\circ-DEC_{Arktur} = 90^\circ - 19.063222^\circ = \underline{\underline{70.936778^\circ}} \nonumber +\end{align} +Um $a$ zu bestimmen, benötigen wir zuerst den Winkel \[\alpha= RA_{Deneb} - RA_{Arktur} = 310.55058^\circ -214.17558^\circ = \underline{\underline{96.375^\circ}}.\] +Danach nutzen wir den sphärischen Winkelkosinussatz, um $a$ zu berechnen: +\begin{align} + a &= \cos^{-1}(\cos(b) \cdot \cos(c) + \sin(b) \cdot \sin(c) \cdot \cos(\alpha)) \nonumber \\ + &= \cos^{-1}(\cos(44.638806) \cdot \cos(70.936778) + \sin(44.638806) \cdot \sin(70.936778) \cdot \cos(96.375)) \nonumber \\ + &= \underline{\underline{80.8707801^\circ}} \nonumber +\end{align} +Für $\beta$ und $\gamma$ nutzen wir den sphärischen Seitenkosinussatz: +\begin{align} + \beta &= \cos^{-1} \bigg[\frac{\cos(b)-\cos(a) \cdot \cos(c)}{\sin(a) \cdot \sin(c)}\bigg] \nonumber \\ + &= \cos^{-1} \bigg[\frac{\cos(44.638806)-\cos(80.8707801) \cdot \cos(70.936778)}{\sin(80.8707801) \cdot \sin(70.936778)}\bigg] \nonumber \\ + &= \underline{\underline{45.0115314^\circ}} \nonumber +\end{align} + + \begin{align} + \gamma &= \cos^{-1} \bigg[\frac{\cos(c)-\cos(b) \cdot \cos(a)}{\sin(a) \cdot \sin(b)}\bigg] \nonumber \\ + &= \cos^{-1} \bigg[\frac{\cos(70.936778)-\cos(44.638806) \cdot \cos(80.8707801)}{\sin(80.8707801) \cdot \sin(44.638806)}\bigg] \nonumber \\ + &=\underline{\underline{72.0573328^\circ}} \nonumber +\end{align} +\subsection{Dreieck $BPC$} +\begin{wrapfigure}{R}{5.6cm} + \includegraphics{papers/nav/bilder/position3.pdf} + \caption{Dreieck BPC} +\end{wrapfigure} +Als nächstes berechnen wir die Seiten $h_b$, $h_c$ und die Innenwinkel $\beta_1$ und $\gamma_1$. +\begin{align} + h_b&=90^\circ - h_b \nonumber \\ + &= 90^\circ - 47.42744^\circ \nonumber \\ + &= \underline{\underline{42.572556^\circ}} \nonumber +\end{align} +\begin{align} + h_c &= 90^\circ - h_c \nonumber \\ + &= 90^\circ - 50.256027^\circ \nonumber \\ + &= \underline{\underline{39.743973^\circ}} \nonumber +\end{align} +\begin{align} + \beta_1 &= \cos^{-1} \bigg[\frac{\cos(h_c)-\cos(a) \cdot \cos(h_b)}{\sin(a) \cdot \sin(h_b)}\bigg] \nonumber \\ + &= \cos^{-1} \bigg[\frac{\cos(39.743973)-\cos(80.8707801) \cdot \cos(42.572556)}{\sin(80.8707801) \cdot \sin(42.572556)}\bigg] \nonumber \\ + &=\underline{\underline{12.5211127^\circ}} \nonumber +\end{align} +\begin{align} + \gamma_1 &= \cos^{-1} \bigg[\frac{\cos(h_b)-\cos(a) \cdot \cos(h_c)}{\sin(a) \cdot \sin(h_c)}\bigg] \nonumber \\ + &= \cos^{-1} \bigg[\frac{\cos(42.572556)-\cos(80.8707801) \cdot \cos(39.743973)}{\sin(80.8707801) \cdot \sin(39.743973)}\bigg] \nonumber \\ + &=\underline{\underline{13.2618475^\circ}} \nonumber +\end{align} + +\subsection{Dreieck $ABP$} +\begin{wrapfigure}{R}{5.6cm} + \includegraphics{papers/nav/bilder/position4.pdf} + \caption{Dreieck ABP} +\end{wrapfigure} +Als erster müssen wir den Winkel $\beta_2$ berechnen: +\begin{align} + \beta_2 &= \beta + \beta_1 = 45.011513^\circ + 12.5211127^\circ \nonumber \\ + &=\underline{\underline{44.6687451^\circ}} \nonumber +\end{align} +Danach können wir mithilfe von $\beta_2$, $c$ und $h_b$ die Seite $l$ berechnen: +\begin{align} + l &= \cos^{-1}(\cos(c) \cdot \cos(h_b) + \sin(c) \cdot \sin(h_b) \cdot \cos(\beta_2)) \nonumber \\ + &= \cos^{-1}(\cos(70.936778) \cdot \cos(42.572556) + \sin(70.936778) \cdot \sin(42.572556) \cdot \cos(57.5326442)) \nonumber \\ + &= \underline{\underline{54.2833404^\circ}} \nonumber +\end{align} +Damit wir gleich den Längengrad berechnen können, benötigen wir noch den Winkel $\omega$: +\begin{align} + \omega &= \cos^{-1} \bigg[\frac{\cos(h_b)-\cos(c) \cdot \cos(l)}{\sin(c) \cdot \sin(l)}\bigg] \nonumber \\ + &=\cos^{-1} \bigg[\frac{\cos(42.572556)-\cos(70.936778) \cdot \cos(54.2833404)}{\sin(70.936778) \cdot \sin(54.2833404)}\bigg] \nonumber \\ + &= \underline{\underline{44.6687451^\circ}} \nonumber +\end{align} + +\subsection{Längengrad und Breitengrad bestimmen} + +\begin{align} + \delta &= 90^\circ - l \nonumber \\ + &= 90^\circ - 54.2833404 \nonumber \\ + &= \underline{\underline{35.7166596^\circ}} \nonumber +\end{align} +\begin{align} + \lambda &= \lambda_{Arktur} + \omega \nonumber \\ + &= 95.5647759^\circ + 44.6687451^\circ \nonumber \\ + &= \underline{\underline{140.233521^\circ}} \nonumber +\end{align} +Wie wir sehen, stimmen die berechneten Koordinaten mit den Koordinaten des Punktes, an welchem gemessen wurde überein. + +\subsection{Fazit} +Die theoretische Anleitung im Abschnitt 21.6 scheint grundsätzlich zu funktionieren. +Allerdings gab es zwei interessante Probleme. + +Einerseits das Problem, ob der Punkt P sich oberhalb oder unterhalb von $a$ befindet. +Da wir eigentlich wussten, wo der gesuchte Punkt P ist, konnten wir das Dreieck anhand der Koordinaten der Bildpunkte richtig aufstellen. +In der Praxis muss man aber schon wissen, auf welchem Breitengrad man ungefähr ist. +Dies weis man in der Regeln aber, da die eigene Breite die Höhe des Polarsterns ist. +Diese Höhe wird mit dem Sextant gemessen. + +Andererseits ist da noch ein Problem mit dem Sinussatz. +Beim Sinussatz gibt es immer zwei Lösungen, weil \[ \sin(\pi-a)=\sin(a).\] +Da kann es sein (und war in unserem Fall auch so), dass man das falsche Ergebnis erwischt. +Durch diese Erkenntnis haben wir nur Kosinussätze verwendet und dies ebenfalls im Abschnitt 21.6 abgeändert, da es für den Leser auch relevant sein kann, wenn er es Probieren möchte. + + + + diff --git a/buch/papers/nav/images/position/test.tex b/buch/papers/nav/images/position/test.tex index 8f4b341..3247ed1 100644 --- a/buch/papers/nav/images/position/test.tex +++ b/buch/papers/nav/images/position/test.tex @@ -17,7 +17,7 @@ \usepackage{wrapfig} \begin{document} -\begin{wrapfigure}{R}{5.2cm} +\begin{wrapfigure}{R}{5.6cm} \includegraphics{position1-small.pdf} \end{wrapfigure} Lorem ipsum dolor sit amet, consectetuer adipiscing elit. diff --git a/buch/papers/nav/main.tex b/buch/papers/nav/main.tex index 4c52547..37bc83a 100644 --- a/buch/papers/nav/main.tex +++ b/buch/papers/nav/main.tex @@ -15,6 +15,7 @@ \input{papers/nav/sincos.tex} \input{papers/nav/trigo.tex} \input{papers/nav/nautischesdreieck.tex} +\input{papers/nav/bsp.tex} \printbibliography[heading=subbibliography] diff --git a/buch/papers/nav/nautischesdreieck.tex b/buch/papers/nav/nautischesdreieck.tex index d8a14af..44153bd 100644 --- a/buch/papers/nav/nautischesdreieck.tex +++ b/buch/papers/nav/nautischesdreieck.tex @@ -97,7 +97,6 @@ Mithilfe dieser Dreiecken können wir die einfachen Sätze der sphärischen Trig \end{center} \end{figure} - \subsubsection{Dreieck $ABC$} \begin{center} @@ -140,12 +139,9 @@ können wir nun die dritte Seitenlänge bestimmen. Es ist darauf zu achten, dass hier natürlich die Seitenlängen in Bogenmass sind und dementsprechend der Kosinus und Sinus verwendet wird. Jetzt fehlen noch die beiden anderen Innenwinkel $\beta$ und\ $\gamma$. -Diese bestimmen wir mithilfe des Sinussatzes \[\frac{\sin (a)}{\sin (\alpha)} =\frac{\sin (b)}{\sin (\beta)} = \frac{\sin (c)}{\sin (\gamma)}.\] -Hier muss man aufpassen, dass man Seite von Winkel unterscheiden kann. -Im Zähler sind die Seiten, im Nenner die Winkel. -Somit ist \[\beta =\sin^{-1} [\sin(b) \cdot \frac{\sin(\alpha)}{\sin(a)}].\] +Diese bestimmen wir mithilfe des Kosinussatzes: \[\beta=\cos^{-1} \bigg[\frac{\cos(b)-\cos(a) \cdot \cos(c)}{\sin(a) \cdot \sin(c)}\bigg]\] und \[\gamma = \cos^{-1} \bigg[\frac{\cos(c)-\cos(b) \cdot \cos(a)}{\sin(a) \cdot \sin(b)}.\bigg]\] -Schlussendlich haben wir die Seiten $a,b\ und \ c$, die Ecken A,B und C und die Winkel $\alpha$, $\beta$ und $\gamma$ bestimmt und somit das ganze Kugeldreieck $ABC$ berechnet. +Schlussendlich haben wir die Seiten $a$ $b$ und $c$, die Ecken A,B und C und die Winkel $\alpha$, $\beta$ und $\gamma$ bestimmt und somit das ganze Kugeldreieck $ABC$ berechnet. \subsubsection{Dreieck $BPC$} Wir bilden nun ein zweites Dreieck, welches die Ecken $B$ und $C$ des ersten Dreiecks besitzt. @@ -167,8 +163,7 @@ und \delta =\cos^{-1} [\cos(c) \cdot \cos(pb) + \sin(c) \cdot \sin(pb) \cdot \cos(\kappa)]. \] -Für die Geographische Länge $\lambda$ des eigenen Standortes muss man den Winkel $\omega$, welcher sich im Dreieck $ACP$ in der Ecke bei $A$ befindet. -Mithilfe des Sinussatzes \[\frac{\sin (a)}{\sin (\alpha)} =\frac{\sin (b)}{\sin (\beta)} = \frac{\sin (c)}{\sin (\gamma)}\] können wir das bestimmen. -Somit ist \[ \omega=\sin^{-1}[\sin(pc) \cdot \frac{\sin(\gamma)}{\sin(l)}] \]und unsere gesuchte geographische Länge schlussendlich +Für die Geographische Länge $\lambda$ des eigenen Standortes nutzt man den Winkel $\omega$, welcher sich im Dreieck $ACP$ in der Ecke bei $A$ befindet. +Mithilfe des Kosinussatzes können wir \[\omega = \cos^{-1} \bigg[\frac{\cos(pb)-\cos(c) \cdot \cos(l)}{\sin(c) \cdot \sin(l)}.\bigg]\] berechnen und schlussentlich dann \[\lambda=\lambda_1 - \omega\] wobei $\lambda_1$ die Länge des Bildpunktes $X$ von $C$ ist. diff --git a/buch/papers/nav/packages.tex b/buch/papers/nav/packages.tex index f2e6132..bedaccd 100644 --- a/buch/papers/nav/packages.tex +++ b/buch/papers/nav/packages.tex @@ -9,4 +9,4 @@ %\usepackage{packagename} \usepackage{amsmath} -\usepackage{cancel}
\ No newline at end of file +\usepackage{cancel} |