EphemPedia

Anhänger der rechnenden Astronomie

Benutzer-Werkzeuge

Webseiten-Werkzeuge

Zuletzt angesehen: • auf-_und_untergangszeiten
auf-_und_untergangszeiten

Unterschiede

Hier werden die Unterschiede zwischen zwei Versionen der Seite angezeigt.

Link zu der Vergleichsansicht

Beide Seiten, vorherige ÜberarbeitungVorherige Überarbeitung
Nächste Überarbeitung		Vorherige Überarbeitung
auf-_und_untergangszeiten [2024/06/01 22:53] – [Alternative Version von J. Meeus] quernauf-_und_untergangszeiten [2025/06/12 01:16] (aktuell) – [Alternative Version von J. Meeus] hcgreier
Zeile 61: Zeile 61:
 ---- ----
  
-Für die Sonne gilt als Auf/Untergangshöhe $h = -50' = -0\overset{\circ}{.}833333$. Man benötigt weiters die geozentrischen äquatorialen Koordinaten $\alpha_{\odot}, \delta_{\odot}$ für die Sonne um $\textrm{00:00}\;UT$ des gegebenen Tages und die mittlere Sternzeit $GMST$ in Greenwich für denselben Zeitpunkt, sowie die mittlere Schiefe der Ekliptik $\varepsilon_0$.+Für die Sonne gilt als Auf/Untergangshöhe $h = -50' = -0\overset{\circ}{.}833333$. Man benötigt weiter die geozentrischen äquatorialen Koordinaten $\alpha_{\odot}, \delta_{\odot}$ für die Sonne um $\textrm{00:00}\;UT$ des gegebenen Tages und die mittlere Sternzeit $GMST$ in Greenwich für denselben Zeitpunkt, sowie die mittlere Schiefe der Ekliptik $\varepsilon_0$.
  
 <WRAP center round tip 100%> <WRAP center round tip 100%>
-Man kann die [[:astronomische_begriffe#topozentrisch|topozentrischen Koordinaten]] der Sonne verwenden. Die Sonnenparallaxe ist aber relativ klein mit etwa $8\overset{''}{.}749$ und liegt weit unterhalb der Genauigkeit dieses Algorithmus. Topozentrische Koordinaten lohnen sich nur beim Erdmond.+Man kann die [[:astronomische_begriffe#topozentrisch|topozentrischen Koordinaten]] der Sonne verwenden. Die Sonnenparallaxe ist aber relativ klein mit etwa $8\overset{''}{.}749$ und liegt weit unterhalb der Genauigkeit dieses Algorithmus. Die topozentrischen Koordinaten lohnen sich nur beim Erdmond.
 </WRAP> </WRAP>
  
Zeile 106: Zeile 106:
  
 <WRAP center round tip 100%> <WRAP center round tip 100%>
-Da man sich nicht allzu weit von der Epoche $J2000$ entfernt befindet kann man die Terme für $T^2$ und $T^3$ hier auch vernachlässigen.+Da man sich nicht allzu weit von der Epoche $J2000$ entfernt befindet, kann man die Terme für $T^2$ und $T^3$ hier auch vernachlässigen.
 </WRAP> </WRAP>
  
Zeile 136: Zeile 136:
 \end{align}\) | \end{align}\) |
  
-Dies ist ein Stundenwert und wird durch Mulitplikation mit $15$ in Grad umgewandelt.+Dies ist ein Stundenwert und wird durch Multiplikation mit $15$ in Grad umgewandelt.
  
 \(\begin{align} \(\begin{align}
Zeile 143: Zeile 143:
 \end{align}\) \end{align}\)
  
-Nun hat man sämliche Werte für $\alpha_{\odot}, \delta_{\odot}, GMST$ sowie die geografische Länge $\lambda_0 = -10^{\circ}$ (Ost negativ!) in Grad gegeben. Die Kulminationszeit $K$ ist dann+Nun hat man sämtliche Werte für $\alpha_{\odot}, \delta_{\odot}, GMST$, sowie die geografische Länge $\lambda_0 = -10^{\circ}$ (Ost negativ!) in Grad gegeben. Die Kulminationszeit $K$ ist dann
  
 \(\begin{align} \(\begin{align}
Zeile 205: Zeile 205:
 \end{align}\) \end{align}\)
  
-Diese Zeitpunkte sind natürlich in Weltzeit $UT$. Die Umrechnung in die gewünschte [[:zeiteingabe#zeitzone|Zeitzone]] kann noch angebracht werden.+Diese Zeitpunkte gelten natürlich für die lokale [[:zeiteingabe#zeitzone|Zeit]] (hier: MOZ). Die Umrechnung in die Sommerzeit kann noch angebracht werden. Sie gilt dann allerdings für den lokalen Meridian.
  
 Die Astronomiesoftware GUIDE gibt für den gegebenen Tag die folgenden Daten an: Die Astronomiesoftware GUIDE gibt für den gegebenen Tag die folgenden Daten an:
Zeile 216: Zeile 216:
  
 <WRAP center round important 100%> <WRAP center round important 100%>
-Bei Höhen von wenigen Grad bzw. wenn sich Objekte an ihrem Auf-/Untergangspunkt befinden, sind die Ergebnisse der Formeln für die [[:koordinatenreduktion#refraktion|Refraktion]] mit Vorsicht zu beurteilen. In der Nähe des Horizonts spielen unvorhersehbare Störungen der Atmosphäre eine große Rolle. Nach Untersuchungen von [[:literaturhinweise#paper_schaelill|B. Schaefer und W. Liller]] (1990) schwankt die Refraktion am Horizont normalerweise um $0\overset{\circ}{.}3 = 18'$ um den Mittelwert, teilweise auch deutlich mehr. **Es macht daher keinen Sinn, die Auf- oder Untergangszeiten eines Himmelskörpers genauer als auf ganze Minuten anzugeben**.+Bei Höhen von wenigen Grad bzw. wenn sich Objekte an ihrem Auf-/Untergangspunkt befinden, sind die Ergebnisse der Formeln für die [[:koordinatenreduktion#refraktion|Refraktion]] mit Vorsicht zu beurteilen. In der Nähe des Horizonts spielen unvorhersehbare Störungen der Atmosphäre eine große Rolle. Nach Untersuchungen von [[:literaturhinweise#paper_schaelill|B. Schaefer und W. Liller]] (1990) schwankt die Refraktion am Horizont normalerweise um $0\overset{\circ}{.}3 = 18'$ um den Mittelwert, teilweise auch deutlich mehr. **Es macht daher keinen Sinn, die Auf- oder Untergangszeiten eines Himmelskörpers genauer als auf //ganze Minuten// anzugeben**.
 </WRAP> </WRAP>
  
Zeile 235: Zeile 235:
 $$\Theta_{i+1} = \Theta_i + 1.00273790931\cdot (T_{i+1} − T_i)\tag{8}$$ $$\Theta_{i+1} = \Theta_i + 1.00273790931\cdot (T_{i+1} − T_i)\tag{8}$$
  
-im nachfolgenden Durchgang. Die Höhe lautet:+im nachfolgenden Durchgang. Es wird die "//geozentrische//" Höhe verwendet, sie lautet für den Mond:
  
 $$h_0 = + 0^{\circ} 08' = 0\overset{\circ}{.}1333\dots\tag{9}$$ $$h_0 = + 0^{\circ} 08' = 0\overset{\circ}{.}1333\dots\tag{9}$$
Zeile 259: Zeile 259:
 $$T_{i+1} = T_i + \frac{t_i - \tau_i}{\mu_i}\tag{14}$$ $$T_{i+1} = T_i + \frac{t_i - \tau_i}{\mu_i}\tag{14}$$
  
-Hat man den neuen Wert $T_i$ berechnet, so muss man nun die geozentrisch äquatorialen Koordianten des Mondes $\alpha_i$ und $\delta_i$ erneut berechnen und dann mit den neuen Werten zur Sternzeit $\Theta_i$ springen. Der Durchlauf startet solange erneut, bis ein genügend geringer Differenzwert zwischen $T_{i+1}$ und $T_i$ gefunden wird.+Hat man den neuen Wert $T_i$ berechnet, so muss man nun die geozentrisch äquatorialen Koordinaten des Mondes $\alpha_i$ und $\delta_i$ erneut berechnen und dann mit den neuen Werten zur Sternzeit $\Theta_i$ springen. Der Durchlauf startet solange erneut, bis ein genügend geringer Differenzwert zwischen $T_{i+1}$ und $T_i$ gefunden wird.
  
 <WRAP center round info> <WRAP center round info>
-Das iterative Verfahren wird in Stunden berechnet, nicht in Grad. Man sollte während der Iteration eine [[mathematische grundlagen#reduktionsfunktion|Reduktion]] auf ein Intervall zwischen $0^h$ und $24^h$ vermeiden, denn das stört die Iteration erheblich. +Das iterative Verfahren wird in **Stunden** berechnet, nicht in Grad. Man sollte während der Iteration eine [[mathematische grundlagen#reduktionsfunktion|Reduktion]] auf ein Intervall zwischen $0^h$ und $24^h$ vermeiden, denn das stört die Iteration erheblich. 
 </WRAP> </WRAP>
  
Zeile 381: Zeile 381:
 \end{align}\) \end{align}\)
  
-Damit berechnet man nun sukzessive die neuen Mondkoordinaten sowie den neuen Stundenwinkel analog zu Schritt "0" und erhält+Damit berechnet man nun sukzessive die neuen Mondkoordinaten, sowie den neuen Stundenwinkel analog zu Schritt "0" und erhält
  
 \(\begin{align} \(\begin{align}
Zeile 469: Zeile 469:
 <imgcaption image2|Prinzip der Nullstellensuche mit quadratischer Interpolation>{{ :methodik_auf_unter_montenbruck.png?800 |}}</imgcaption> <imgcaption image2|Prinzip der Nullstellensuche mit quadratischer Interpolation>{{ :methodik_auf_unter_montenbruck.png?800 |}}</imgcaption>
  
-Man beginnt für beide Kurven mit den ersten Zeitpunkt $t = 1^h$ bis $t = 24^h$, berechnet den Stundenwinkel $\tau(t)$ als Startwert und ermittelt daraufhin die drei in der Graphik gezeigten Stützpunkte für +Man beginnt für beide Kurven mit den ersten Zeitpunkt $t = 1^h$ bis $t = 24^h$, berechnet den Stundenwinkel $\tau(t)$ als Startwert und ermittelt daraufhin die drei in der Graphik gezeigten Stützpunkte für den gesuchten Zeitpunkt $x$:
-den gesuchten Zeitpunkt $x$:+
  
 {{tablelayout?rowsHeaderSource=1&colwidth="220px,170px"&float=center}} {{tablelayout?rowsHeaderSource=1&colwidth="220px,170px"&float=center}}
Zeile 516: Zeile 515:
 ---- ----
  
-Man berechnet als erstes die Höhe h im einem Intervall von 26 Stunden mit Hilfe von Gleichung (16):+Man berechnet als erstes die Höhe $h$ im einem Intervall von 26 Stunden mit Hilfe von Gleichung (16):
  
 {{tablelayout?rowsHeaderSource=1&colwidth="100px,300px"&float=center}} {{tablelayout?rowsHeaderSource=1&colwidth="100px,300px"&float=center}}
Zeile 573: Zeile 572:
 ==== Alternative Version von J. Meeus ==== ==== Alternative Version von J. Meeus ====
  
-J. Meeus geht ebenfalls den Weg über die Interpolation, aber das Ziel sind die äquatorialen Koordinaten des Himmelsobjekts, die über einem kurzen Zeitraum von drei Tagen berechnet werden müssen. Die Vorgehensweise ist in Analogie zum Abschnitt [[:auf-_und_untergangszeiten#sonne_sterne_und_planeten]]. Es folgt nun die Berechnung der [[:koordinatentransformation#geozentrische_koordinaten|geozentrisch - äquatorialen Koordinaten]] des entsprechenden Himmelskörpers zu folgenden drei Zeitpunkten um 0$^h$ dynamischer Zeit:+J. Meeus geht ebenfalls den Weg über die Interpolation, aber das Ziel sind die äquatorialen Koordinaten des Himmelsobjekts, die über einem kurzen Zeitraum von drei Tagen berechnet werden müssen. Die Vorgehensweise ist in Analogie zum Abschnitt [[:auf-_und_untergangszeiten#sonne_sterne_und_planeten|Sonne, Sterne und Planeten]]. Es folgt nun die Berechnung der [[:koordinatentransformation#geozentrische_koordinaten|geozentrisch - äquatorialen Koordinaten]] des entsprechenden Himmelskörpers zu folgenden drei Zeitpunkten um 0$^h$ dynamischer Zeit:
  
-$\alpha_1$, $\delta_1$ am Tag d $-$ 1$^d$ um 0$^h$ \\ +\[\begin{align} 
-$\alpha_2$, $\delta_2$ am Tag d um 0$^h$ \\ +\alpha_1,& \delta_1 \quad\textsf{am Tag}\quad d - 1^d \quad\textsf{um}\quad 0^h \\ 
-$\alpha_3$, $\delta_3$ am Tag d + 1$^d$ um 0$^h$+\alpha_2,& \delta_2 \quad\textsf{am Tag}\quad d \quad\textsf{um}\quad 0^h \\ 
 +\alpha_3,& \delta_3 \quad\textsf{am Tag}\quad d + 1^d \quad\textsf{um}\quad 0^h 
 +\end{align}\]
  
 Mit Hilfe der maßgebenden Koordinaten $\alpha_2$, $\delta_2$ wird über Gleichung (1) der halbe Tagbogen t berechnet. Es gilt mit $\mu_p$ = 1 und den Gleichungen (2), (4) und (5): Mit Hilfe der maßgebenden Koordinaten $\alpha_2$, $\delta_2$ wird über Gleichung (1) der halbe Tagbogen t berechnet. Es gilt mit $\mu_p$ = 1 und den Gleichungen (2), (4) und (5):
 +
 \[\begin{split} \[\begin{split}
 m_0 &= K\\ m_0 &= K\\
Zeile 590: Zeile 592:
  
 Nun wird die Interpolation der drei Mondpositionen gebraucht: Nun wird die Interpolation der drei Mondpositionen gebraucht:
-$$y = y_2 + \frac{n}{2}\left(y_3 - y_1 + n \cdot (y_1 + y_3 - 2\cdot y_2)\right) \quad\text{mit}\quad n = \frac{m_k}{24^h}\tag{25}$$+$$y = y_2 + \frac{n}{2}\cdot\left[ y_3 - y_1 + n \cdot (y_1 + y_3 - 2\cdot y_2) \right] \quad\text{mit}\quad n = \frac{m_k}{24^h}\tag{25}$$
  
 Für $y_i$ ist der jeweilige Wert von $\alpha_i$ oder $\delta_i$ einzusetzen. Man erhält dann die interpolierten Koordinaten y = $\alpha$ oder y = $\delta$. Um einen Sprung in der Darstellung von $24^h$ auf $0^h$ zu vermeiden, muss folgende Relation beachtet werden: Für $y_i$ ist der jeweilige Wert von $\alpha_i$ oder $\delta_i$ einzusetzen. Man erhält dann die interpolierten Koordinaten y = $\alpha$ oder y = $\delta$. Um einen Sprung in der Darstellung von $24^h$ auf $0^h$ zu vermeiden, muss folgende Relation beachtet werden:
  
-Zähle i = 1 bis 2 \\ +Zähle $i = 1$ bis $2$ \\ 
-Falls |$\alpha_{i+1} - \alpha_i$| > 12$^h$ dann gilt: \\ +$\quad$ Falls |$\alpha_{i+1} - \alpha_i$| > $12^h$ dann gilt: \\ 
-Falls $\alpha_{i+1}$ > $\alpha_i$ dann gilt: \\ +$\qquad$ Falls $\alpha_{i+1} > \alpha_i$ dann gilt: \\ 
-$\alpha_i$ = $\alpha_i$ + 24$^h$ \\ +$\quad\qquad$ $\alpha_i = \alpha_i + 24^h$ \\ 
-sonst gilt: \\ +$\qquad$ sonst gilt: \\ 
-$\alpha_{i+1}$ = $\alpha_{i+1}$ + 24$^h$ \\ +$\quad\qquad$ $\alpha_{i+1} = \alpha_{i+1} + 24^h$ \\ 
-Ende \\ +$\qquad$ Ende \\ 
-Ende \\ +$\quad$ Ende \\ 
-Nächstes i \\+Nächstes $i$
  
 Ab hier werden nur noch die interpolierten Koordinaten des Mondes gebraucht. Der Stundenwinkel $\tau$ stammt aus Gleichung (10): Ab hier werden nur noch die interpolierten Koordinaten des Mondes gebraucht. Der Stundenwinkel $\tau$ stammt aus Gleichung (10):
Zeile 613: Zeile 615:
 $$\Delta m = \frac{H-h_0}{15^h\cdot\cos(\delta)\cdot\cos(\beta_0)\cdot\sin(15^h\cdot\tau)}\tag{28}$$ $$\Delta m = \frac{H-h_0}{15^h\cdot\cos(\delta)\cdot\cos(\beta_0)\cdot\sin(15^h\cdot\tau)}\tag{28}$$
  
-Die Höhe $h_0$ stammt aus Tabelle 1. Die Höhe H bekommt man mit:+Die Höhe $h_0$ stammt aus Tabelle 1. Die Höhe $H$ bekommt man mit:
 $$\sin(H) = \sin(\beta_0)\cdot\sin(\delta) + \cos(\beta_0)\cdot\cos(\delta) \cdot \cos(15^h\cdot\tau)\tag{29}$$ $$\sin(H) = \sin(\beta_0)\cdot\sin(\delta) + \cos(\beta_0)\cdot\cos(\delta) \cdot \cos(15^h\cdot\tau)\tag{29}$$
  
-Der gesuchte Zeitpunkt T ist dann:+Der gesuchte Zeitpunkt $T$ ist dann:
 $$T = m_k + \Delta m\tag{30}$$ $$T = m_k + \Delta m\tag{30}$$
  
 <WRAP center round box 100%> <WRAP center round box 100%>
 +
 ==== Beispiel 4 ==== ==== Beispiel 4 ====
  
Zeile 626: Zeile 629:
 ---- ----
  
-Als erstes werden die geozentrisch - äquatorialen koordinaten des Mondes für die drei Zeitpunkte am 15.1., am 16.1 und am 17.1. bestimmt: +Als erstes werden die geozentrisch - äquatorialen Koordinaten des Mondes für die drei Zeitpunkte am 15.1., am 16.1 und am 17.1. bestimmt, z.B. mit dem "schnellen" Algorithmus laut [[:erdmond#mondposition|dieser Methode]]: 
-\[\begin{array}{cc} + 
-15.01.2001:\quad \alpha_1 = 12\overset{h}{.}455588640675865 & \delta_1 = +2\overset{\circ}{.}7180513083004127 \\      +\(\begin{align} 
-16.01.2001:\quad \alpha_2 = 13\overset{h}{.}295955922809693 & \delta_2 = -2\overset{\circ}{.}5095094272111229 \\     +15.01.2001:\quad \alpha_1 &= 12\overset{h}{.}455588640675865 & \delta_1 = +2\overset{\circ}{.}7180513083004127 \\      
-17.01.2001:\quad \alpha_3 = 14\overset{h}{.}116613997835508 & \delta_3 = -7\overset{\circ}{.}4619042104156623 +16.01.2001:\quad \alpha_2 &= 13\overset{h}{.}295955922809693 & \delta_2 = -2\overset{\circ}{.}5095094272111229 \\     
-\end{array}\]+17.01.2001:\quad \alpha_3 &= 14\overset{h}{.}116613997835508 & \delta_3 = -7\overset{\circ}{.}4619042104156623 
 +\end{align}\)
  
 Der halbe Tagbogen und die Durchgangszeit $m_0$ wurde ermittelt zu: Der halbe Tagbogen und die Durchgangszeit $m_0$ wurde ermittelt zu:
-$$t = 5\overset{h}{.}7865382720225993$$ + 
-$$m_0 = 5\overset{h}{.}9320637476530429$$+\(\begin{align} 
 +t &= 5\overset{h}{.}7865382720225993\\ 
 +m_0 &= 5\overset{h}{.}9320637476530429 
 +\end{align}\)
    
 Daraus ergibt sich die erste Abschätzung der Aufgangszeit: Daraus ergibt sich die erste Abschätzung der Aufgangszeit:
-$$m_1 = 0\overset{h}{.}14552547563044360 = 00^h 08^m 44^s$$+ 
 +$m_1 = 0\overset{h}{.}14552547563044360 = 00^h 08^m 44^s$
  
 Die interpolierten Koordinaten des Erdtrabanten sind: Die interpolierten Koordinaten des Erdtrabanten sind:
-$$\alpha = 13\overset{h}{.}300991425177983\qquad\delta = -2\overset{\circ}{.}5403677618275222$$ 
  
-Der Stundenwinkel $\tau$ und die Höhe H sind: +\(\begin{align} 
-$$\tau = -5\overset{h}{.}7911753388363190$$ +\alpha &= 13\overset{h}{.}300991425177983\\ 
-$$H = 6\overset{\circ}{.}5070941042247865 \cdot 10^{-2} \approx 4'$$+\delta &= -2\overset{\circ}{.}5403677618275222 
 +\end{align}\) 
 + 
 +Der Stundenwinkel $\tau$ und die Höhe $H$ sind: 
 + 
 +\(\begin{align} 
 +\tau &= -5\overset{h}{.}7911753388363190\\ 
 +H &= 6\overset{\circ}{.}5070941042247865 \cdot 10^{-2} \approx 4' 
 +\end{align}\)
  
 Die Korrektur $\Delta m$ lautet: Die Korrektur $\Delta m$ lautet:
-$$\Delta m = 7\overset{h}{.}0973970892369157 \cdot 10^{-3} = 00^h 00^m 25,6^s$$+ 
 +\(\begin{align} 
 +\Delta m &= 7\overset{h}{.}0973970892369157 \cdot 10^{-3} \\ 
 +&= 00^h 00^m 25\overset{s}{.}6 
 +\end{align}\)
  
 Die Korrektur ist damit sehr klein. Daraus resultiert, daß die Meeus-Variante sehr effektiv ist. Und damit hat man die Aufgangszeit des Erdmondes am 16.1.2001: Die Korrektur ist damit sehr klein. Daraus resultiert, daß die Meeus-Variante sehr effektiv ist. Und damit hat man die Aufgangszeit des Erdmondes am 16.1.2001:
-$$T = m_1 + \Delta m$$ + 
-$$T = 0\overset{h}{.}14552547563044360 + 7\overset{h}{.}0973970892369157 \cdot 10^{-3} = 0\overset{h}{.}152622872719 = 00^h 09^m 10^s$$+\(\begin{align} 
 +T &= m_1 + \Delta m\\ 
 +T &= 0\overset{h}{.}14552547563044360 + 7\overset{h}{.}0973970892369157 \cdot 10^{-3}\\ 
 +&= 0\overset{h}{.}152622872719 = 00^h 09^m 10^s 
 +\end{align}\)
  
 in Übereinstimmung mit den Beispielen 2 und 3. in Übereinstimmung mit den Beispielen 2 und 3.
 </WRAP> </WRAP>
 +
 +===== Sichtweite und Kimmtiefe =====
 +
 +Diesen beiden Themen ist eine eigene [[sichtweite_und_kimmtiefe|Seite]] gewidmet.
auf-_und_untergangszeiten.1717275194.txt.gz · Zuletzt geändert: 2024/12/20 01:32 (Externe Bearbeitung)
Falls nicht anders bezeichnet, ist der Inhalt dieses Wikis unter der folgenden Lizenz veröffentlicht: CC Attribution-Share Alike 4.0 International
CC Attribution-Share Alike 4.0 International Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki