auf-_und_untergangszeiten
Unterschiede
Hier werden die Unterschiede zwischen zwei Versionen der Seite angezeigt.
Beide Seiten, vorherige ÜberarbeitungVorherige ÜberarbeitungNächste Überarbeitung | Vorherige Überarbeitung | ||
auf-_und_untergangszeiten [2024/05/30 20:41] – [Numerisches Interpolationsverfahren] quern | auf-_und_untergangszeiten [2025/06/12 01:16] (aktuell) – [Alternative Version von J. Meeus] hcgreier | ||
---|---|---|---|
Zeile 12: | Zeile 12: | ||
| Sonstige: | | Sonstige: | ||
- | Hat man die Höhen gefunden, so berechnet man als nächstes den halben Tagbogen: | + | Hat man die Höhen gefunden, so berechnet man als nächstes den halben Tagbogen, der zwischen $0^h$ und $12^h$ liegen muß: |
$$\cos(15^h\cdot t) = \frac{\sin(h) - \sin(\beta_0) \cdot \sin(\delta)}{\cos(\beta_0) \cdot \cos(\delta)}\tag{1}$$ | $$\cos(15^h\cdot t) = \frac{\sin(h) - \sin(\beta_0) \cdot \sin(\delta)}{\cos(\beta_0) \cdot \cos(\delta)}\tag{1}$$ | ||
Zeile 61: | Zeile 61: | ||
---- | ---- | ||
- | Für die Sonne gilt als Auf/ | + | Für die Sonne gilt als Auf/ |
<WRAP center round tip 100%> | <WRAP center round tip 100%> | ||
- | Man kann die [[: | + | Man kann die [[: |
</ | </ | ||
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. |
</ | </ | ||
Zeile 136: | Zeile 136: | ||
\end{align}\) | | \end{align}\) | | ||
- | Dies ist ein Stundenwert und wird durch Mulitplikation | + | Dies ist ein Stundenwert und wird durch Multiplikation |
\(\begin{align} | \(\begin{align} | ||
Zeile 143: | Zeile 143: | ||
\end{align}\) | \end{align}\) | ||
- | Nun hat man sämliche | + | Nun hat man sämtliche |
\(\begin{align} | \(\begin{align} | ||
Zeile 205: | Zeile 205: | ||
\end{align}\) | \end{align}\) | ||
- | Diese Zeitpunkte | + | Diese Zeitpunkte |
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-/ | + | Bei Höhen von wenigen Grad bzw. wenn sich Objekte an ihrem Auf-/ |
</ | </ | ||
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. | + | im nachfolgenden Durchgang. |
$$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 | + | Hat man den neuen Wert $T_i$ berechnet, so muss man nun die geozentrisch äquatorialen |
<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# | + | Das iterative Verfahren wird in **Stunden** berechnet, nicht in Grad. Man sollte während der Iteration eine [[mathematische grundlagen# |
</ | </ | ||
Zeile 381: | Zeile 381: | ||
\end{align}\) | \end{align}\) | ||
- | Damit berechnet man nun sukzessive die neuen Mondkoordinaten sowie den neuen Stundenwinkel analog zu Schritt " | + | Damit berechnet man nun sukzessive die neuen Mondkoordinaten, sowie den neuen Stundenwinkel analog zu Schritt " |
\(\begin{align} | \(\begin{align} | ||
Zeile 469: | Zeile 469: | ||
< | < | ||
- | 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? | {{tablelayout? | ||
Zeile 512: | Zeile 511: | ||
==== Beispiel 3 ==== | ==== Beispiel 3 ==== | ||
- | {{: | + | {{: |
** | ** | ||
---- | ---- | ||
- | < | + | Man berechnet als erstes die Höhe $h$ im einem Intervall von 26 Stunden mit Hilfe von Gleichung (16): |
- | Beim Mond kann es vorkommen, dass er an einem bestimmten Datum nicht aufgeht | + | |
+ | {{tablelayout? | ||
+ | ^ $t(h)$ | ||
+ | | $\color{# | ||
+ | | $0$ | $-1.4008879269623948$ | ||
+ | | $1$ | $+7.6619905974511200$ | ||
+ | | $2$ | $+16.264024232490367$ | ||
+ | | $3$ | $+23.966091982009797$ | ||
+ | | $4$ | $+30.230720083427663$ | ||
+ | | $5$ | $+34.439308630383543$ | ||
+ | | $6$ | $+36.030017284660161$ | ||
+ | | $7$ | $+34.743031979369547$ | ||
+ | | $8$ | $+30.770848659482951$ | ||
+ | | $9$ | $+24.639260861006612$ | ||
+ | | $10$ | $+16.959030561803836$ | ||
+ | | $11$ | $+8.2681369546535599$ | ||
+ | | $12$ | $-0.99967795779209334$ | ||
+ | | $13$ | $-10.481806306675004$ | ||
+ | | $14$ | $-19.832722193158880$ | ||
+ | | $15$ | $-28.660280873493534$ | ||
+ | | $16$ | $-36.452091347905856$ | ||
+ | | $17$ | $-42.511934105372944$ | ||
+ | | $18$ | $-46.002569282772313$ | ||
+ | | $19$ | $-46.256329565749390$ | ||
+ | | $20$ | $-43.232253621070868$ | ||
+ | | $21$ | $-37.554964310685151$ | ||
+ | | $22$ | $-30.076882290628159$ | ||
+ | | $23$ | $-21.531024404194824$ | ||
+ | | $\color{# | ||
+ | | $\color{# | ||
+ | |||
+ | $t(h) = -1^h$ liegt im Vortag ($23$ Uhr), $t(h) = 24^h$ und $t(h) = 25^h$ liegen im nächsten Tag ($0$ Uhr bzw. $1$ Uhr) | ||
+ | |||
+ | Man kann bereits erkennen, zur welcher Uhrzeit der Mond aufgeht. | ||
+ | |||
+ | $a = -0.23042244468713324$ \\ | ||
+ | $b = +8.8324560797263807$ \\ | ||
+ | $c = +7.6619905974511200$ | ||
+ | |||
+ | bekommt | ||
+ | |||
+ | $y_- = -1.4008879269623948$ \\ | ||
+ | $y_0 = +7.6619905974511200$ \\ | ||
+ | $y_+ = +16.264024232490367$ \\ | ||
+ | |||
+ | Die Nullstellen ermittelt man mit der Mitternachtsformel (18): | ||
+ | |||
+ | $x_2 = -0.84869064901308278$, $x_1$ ist mit $+39.180269379197370$ zu groß und entfällt somit. | ||
+ | |||
+ | Jetzt muss festgestellt | ||
+ | |||
+ | Es gilt: $T_a = t_0 + x_2 = 1^h - 0\overset{h}{.}84869064901308278 = 0\overset{h}{.}151309351 = 00^h 09^m$, in Übereinstimmung mit dem Resultat aus Beispiel 2. | ||
</ | </ | ||
- | Man berechnet als erstes die Höhe h im einem Intervall | + | ==== Alternative Version |
- | {{tablelayout? | + | J. Meeus geht ebenfalls den Weg über die Interpolation, |
- | ^ t[h] ^ Höhe h[°] ^ | + | |
- | | $-$1 | $-$9.9174564728524164 | + | |
- | | 0 | $-$1.4008879269623948 | | + | |
- | | 1 | 7.6619905974511200 | | + | |
- | | 2 | 16.264024232490367 | | + | |
- | | 3 | 23.966091982009797 | | + | |
- | | 4 | 30.230720083427663 | | + | |
- | | 5 | 34.439308630383543 | | + | |
- | | 6 | 36.030017284660161 | | + | |
- | | 7 | 34.743031979369547 | | + | |
- | | 8 | 30.770848659482951 | | + | |
- | | 9 | 24.639260861006612 | | + | |
- | | 10 | 16.959030561803836 | | + | |
- | | 11 | 8.2681369546535599 | | + | |
- | | 12 | $-$0.99967795779209334 | | + | |
- | | 13 | $-$10.481806306675004 | | + | |
- | | 14 | $-$19.832722193158880 | | + | |
- | | 15 | $-$28.660280873493534 | | + | |
- | | 16 | $-$36.452091347905856 | | + | |
- | | 17 | $-$42.511934105372944 | | + | |
- | | 18 | $-$46.002569282772313 | | + | |
- | | 19 | $-$46.256329565749390 | | + | |
- | | 20 | $-$43.232253621070868 | | + | |
- | | 21 | $-$37.554964310685151 | | + | |
- | | 22 | $-$30.076882290628159 | | + | |
- | | 23 | $-$21.531024404194824 | | + | |
- | | 24 | $-$1.4008879269623948 | | + | |
- | | 25 | $-$3.9230286546594542 | | + | |
- | -1 liegt im Vortag (23h Uhr), 24h und 25 liegen im nächsten | + | \[\begin{align} |
+ | \alpha_1,& | ||
+ | \alpha_2,& \delta_2 \quad\textsf{am | ||
+ | \alpha_3,& | ||
+ | \end{align}\] | ||
- | Man kann bereits erkennen, zur welcher Uhrzeit | + | Mit Hilfe der maßgebenden Koordinaten $\alpha_2$, $\delta_2$ wird über Gleichung (1) der halbe Tagbogen |
- | a = $-$0.23042244468713324 | + | \[\begin{split} |
- | b = 8.8324560797263807 | + | m_0 &= K\\ |
- | c = 7.6619905974511200 | + | m_1 &= K - t\\ |
+ | m_2 &= K + t | ||
+ | \end{split}\tag{23}\] | ||
- | bekommt | + | $m_k$ steht für $m_0$ , $m_1$, $m_2$ und muß zwischen $0^h$ und $24^h$ liegen ([[: |
+ | $$\Theta_k = LMST + 1.00273790931 \cdot m_k\tag{24}$$ | ||
- | $y_-$ = $-$1.4008879269623948 | + | Nun wird die Interpolation der drei Mondpositionen gebraucht: |
- | $y_0$ = +7.6619905974511200 | + | $$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}$$ |
- | $y_+$ = +16.264024232490367 \\ | + | |
- | Die Nullstellen ermittelt man mit der Mitternachtsformel: | + | 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: |
- | $x_2$ = $-$0.84869064901308278, | + | Zähle |
+ | $\quad$ Falls |$\alpha_{i+1} | ||
+ | $\qquad$ Falls $\alpha_{i+1} > \alpha_i$ dann gilt: \\ | ||
+ | $\quad\qquad$ $\alpha_i = \alpha_i + 24^h$ \\ | ||
+ | $\qquad$ sonst gilt: \\ | ||
+ | $\quad\qquad$ $\alpha_{i+1} = \alpha_{i+1} + 24^h$ \\ | ||
+ | $\qquad$ Ende \\ | ||
+ | $\quad$ Ende \\ | ||
+ | Nächstes $i$ | ||
- | Jetzt muss festgestellt | + | Ab hier werden |
+ | $$\tau | ||
- | Es gilt: $T_a = t_0 + x_2$ = 1 $-$0.84869064901308278 | + | Durchgang: |
+ | $$\Delta m = - \tau\tag{27}$$ | ||
+ | |||
+ | Auf-/ | ||
+ | $$\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: | ||
+ | $$\sin(H) = \sin(\beta_0)\cdot\sin(\delta) | ||
+ | |||
+ | Der gesuchte Zeitpunkt $T$ ist dann: | ||
+ | $$T = m_k + \Delta m\tag{30}$$ | ||
+ | |||
+ | <WRAP center round box 100%> | ||
+ | |||
+ | ==== Beispiel 4 ==== | ||
+ | |||
+ | {{: | ||
+ | ** | ||
+ | ---- | ||
+ | |||
+ | 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 " | ||
+ | |||
+ | \(\begin{align} | ||
+ | 15.01.2001: | ||
+ | 16.01.2001: | ||
+ | 17.01.2001: | ||
+ | \end{align}\) | ||
+ | |||
+ | Der halbe Tagbogen und die Durchgangszeit | ||
+ | |||
+ | \(\begin{align} | ||
+ | t &= 5\overset{h}{.}7865382720225993\\ | ||
+ | m_0 &= 5\overset{h}{.}9320637476530429 | ||
+ | \end{align}\) | ||
+ | |||
+ | Daraus ergibt sich die erste Abschätzung der Aufgangszeit: | ||
+ | |||
+ | $m_1 = 0\overset{h}{.}14552547563044360 | ||
+ | |||
+ | Die interpolierten Koordinaten des Erdtrabanten sind: | ||
+ | |||
+ | \(\begin{align} | ||
+ | \alpha &= 13\overset{h}{.}300991425177983\\ | ||
+ | \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: | ||
+ | |||
+ | \(\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: | ||
+ | |||
+ | \(\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 | ||
</ | </ | ||
+ | |||
+ | ===== Sichtweite und Kimmtiefe ===== | ||
+ | |||
+ | Diesen beiden Themen ist eine eigene [[sichtweite_und_kimmtiefe|Seite]] gewidmet. |
auf-_und_untergangszeiten.1717094484.txt.gz · Zuletzt geändert: 2024/12/20 01:32 (Externe Bearbeitung)