1
0

Ultimate version and info about converting times.

This commit is contained in:
Markus Birth 1999-09-10 13:50:43 +02:00
parent 33b97e81df
commit d73ad1e94e
31 changed files with 771 additions and 2 deletions

43
info/JD.html Normal file
View File

@ -0,0 +1,43 @@
<HTML>
<HEAD>
<TITLE>Julianisches Datum</TITLE>
</HEAD>
<BODY>
<SCRIPT LANGUAGE="JavaScript"><!--
function juld(frm) {
var h = Math.floor(frm.h.value);
var m = Math.floor(frm.m.value);
var s = Math.floor(frm.s.value);
var day = Math.floor(frm.d.value);
var month = Math.floor(frm.mo.value);
var year = Math.floor(frm.y.value);
var hdec = Math.floor(h) + (Math.floor(m) / 60) + (Math.floor(s) / 3600);
frm.hdec.value = hdec;
year += 4713; // JD begun 01 Jan 4713 B.C. - 12:00h
hdec -= 12;
if (hdec<0) { hdec+=24; day--; }
var modelta=new Array(0,0,31,59,90,120,151,181,212,243,273,304,334);
jdate = Math.floor(year/4) - Math.floor(year/100) + Math.floor(year/400);
if (month<3 && year%4==0 && year/100-4*year/400!=0) { jdate--; }
jdate += year*365+modelta[month]+day+hdec/24;
if (year>1600) { jdate -= 328; }
frm.jdate.value = jdate;
}
// -->
</SCRIPT>
<FORM NAME="jd">
<INPUT TYPE=text SIZE=2 NAME="h">:<INPUT TYPE=text SIZE=2 NAME="m">.<INPUT TYPE=text SIZE=6 NAME="s" VALUE="00.000"><BR>
<INPUT TYPE=text SIZE=40 NAME="hdec"><BR>
<INPUT TYPE=text SIZE=2 NAME="d">.<INPUT TYPE=text SIZE=2 NAME="mo">.<INPUT TYPE=text SIZE=4 NAME="y">
<P>
<INPUT TYPE=text SIZE=40 NAME="jdate">
<P>
<INPUT TYPE=button VALUE="Calculate" onClick="juld(this.form);">
</FORM>

242
info/Zeitsysteme.html Normal file
View File

@ -0,0 +1,242 @@
<html>
<head>
<title>Vertiefte Informationen Geodätische Astronomie, Zeitsysteme</title>
<link rel=stylesheet href=/scherer/css/astro.css type="text/css">
</head>
<body>
<h1>Geodätische Astronomie, Zeitsysteme</h1>
<dl>
<dt>Zeitsysteme in der Geodätischen Astronomie
<dd><a href="#UT">Universal Time</a> UT
<dd><a href="#GMST">Greenwich Sidereal Time</a> GMST, GAST
<dd><a href="#LMST">Local Sidereal Time</a> LMST, LAST
<dd><a href="#Rechnung">Berechnungen</a>
</dl>
<p class=bs>
In der Astronomie benötigt man verschiedene Zeitsysteme für die Darstellung
der sich zeitlich verändernden Verhältnisse zwischen den Himmelskörpern und
der Situation auf der Erde als Beobachtungsplattform. Im alltäglichen Leben
spielt dabei die <strong>Sonnenzeit</strong> eine zentrale Rolle, während zur
Beschreibung der Himmelsmechanik sich die <strong>Sternzeit</strong> weit
besser eignet.
<p class=bs>
Der wesentliche Unterschied beider Zeitsysteme ist der Bezug: 24h Sonnenzeit
ist der Zeitabschnitt zweier aufeinanderfolgender oberer Kulminationen der
Sonne, 24h Sternzeit ist der Zeitabschnitt zweier aufeinanderfolgender oberer
Kulminationen eines Sterns. Infolge der Translation der Erde um die Sonne
relativ zum System der Sterne variieren die Zeitsysteme jährlich um 1 Tag, also
ca 4m pro Tag. Dabei gilt:
<p>
<a name="(1)"><table width="100%" border=0>
<tr>
<td class=l>24h Sternzeit = 23h56m04s Sonnenzeit</td>
<td class=r>(1)</td>
</table></a>
<p class=bs>
Als Bezugsrichtung der Sternzeit wird anstatt eines Sterns der Frühlingspunkt
verwendet. Er befindet sich auf der Schnittgeraden von Himmelsäquator- und
Ekliptikebene.
<p class=bs>
Hinsichtlich der Zeiten kann man noch zwischen <strong>mittleren</strong> und
<strong>wahren</strong> Zeitangaben unterscheiden. Im Rahmen der mittleren
Zeitangaben wird die Bewegung der Referenz (Sonne/Frühlingspunkt) unter der
Annahme
einer ungestörten Ekliptik beschrieben (einfaches Präzessionsmodell), während
die wahren Zeitangaben aus der zusätzlichen Berücksichtigung des
Nutationsmodells resultieren. Beide Angaben unterscheiden sich:
<dl>
<dt>im Fall der Sternzeit
<dd>um die <em>Equation of Equinox</em> EqE
<p>
<dt>im Fall der Sonnenzeit
<dd>um die <em>Equation of Time</em> EqT (Zeitgleichung)
</dl>
<hr>
<a name="UT"><h2>UT</h2></a>
<p class=bs>
UT (Universal Time) ist eine definierte gleichförmige Zeit (Atomuhren), die auf
das Zeitsystem Sonnenzeit bezogen ist. UT stellt eine Annäherung an die
mittlere Sonnenzeit UT1 dar, die jedoch nicht gleichförmig ist. Infolgedessen
muß die gleichförmige Zeit der Atomuhren um ganze Schaltsekunden verändert
werden, falls |UT-UT1|>0.7s wird. Die aktuellen Unterschiede zwischen UT und UT1
kann man den ERP-Dateien (earth rotation parameters) des
<a href="http://hpiers.obspm.fr/">IERS</a> entnehmen.
<p class=bs>
Zeitangaben in UT beziehen sich auf den Nullmeridian von Greenwich. UT1 wird
auch GMT (Greenwich Mean Time) genannt. Die zugehörige wahre Zeitangabe GAT
(Greenwich Apparent Time) unterscheidet sich von GMT um die Zeitgleichung EqT je
nach Jahreszeit um mehrere Minuten. Analog definiert sind die lokalen
Sonnenzeiten LMT (local mean time) und LAT (local apparent time). Sie bestimmen
sich zu:
<p>
<a name="(2)"><table width="100%" border=0>
<tr>
<td class=l>{LMT,LAT} = {GMT,GAT} + <i>lambda</i>/15</td>
<td class=r>(2)</td>
</table></a>
<p class=bs>
mit der geografischen Länge <i>lambda</i> des Ortsmeridians in [°]. Der aktuelle
Wert der Zeitgleichung EqT kann z.B. dem <quote>Star Almanac for Land
Surveyors</quote> entnommen werden.
<p class=bs>
Bei der astro-geodätischen Beobachtung wird im folgenden davon ausgegangen, daß
UT durch den Einsatz entsprechender Zeitzeichenempfänger (z.B. Zeitzeichen DCF77
der <a href="http://www.ptb.de/deutsch/org/4/43/hp.htm">Physikalisch-Technischen
Bundesanstalt (PTB) in Braunschweig</a>) vorliegt.
<p>
<hr>
<a name="GMST"><h2>GMST / GAST</h2></a>
<p class=bs>
Die Greenwicher Sternzeitangaben GMST und GAST berechnen sich aus dem Winkel
(genannt <quote>Stundenwinkel</quote> zwischen der raumfesten Richtung zum
Frühlingspunkt und dem Nullmeridian von Greenwich, der an der Erdrotation
teilnimmt. Der Unterschied zwischen der mittleren und der wahren Zeitdarstellung
ergibt sich aus dem Grad der Modellierung der Erdrotation (ohne/mit
Nutationsmodell) und wird durch die EqE beschrieben.
<p class=bs>
Wie die Sonnenzeit <a href="#UT">UT</a> bildet der Nullmeridian von Greenwich
die Bezugsrichtung im äquatorialen Koordinatensystem. Zwischen UT und GMST
existiert eine einfache Beziehung, die auf den Bezugsepoche J.2000 (01-Jan-2000,
12.00UT) aufsetzt:
<p>
<a name="(3)"><table width="100%" border=0>
<tr>
<td class=l>GMST = 6h41m50.548s + UT + (8640184.812866 <i>T</i>
+ 0.093104 <i>T</i><sup>2</sup>
- 0.0000062 <i>T</i><sup>3</sup>)s</td>
<td class=r>(3)</td>
</table></a>
<p class=bs>
Dabei bezeichnet <i>T</i> den zeitlichen Abstand des aktuellen Zeitpunktes UT
von der Bezugsepoche in Julianischen Jahrhunderten zu je 36525 Tagen, wie er
nach <a href="#(4)">(4)</a> berechnet werden kann. Für den 01-Jul-1998
beispielsweise beträgt <i>T</i> = -1.500.
<p>
<a name="(4)"><table width="100%" border=0>
<tr>
<td class=l><i>T</i> = (<i>JD</i> - 2451545.0)/36525</td>
<td class=r>(4)</td>
</table></a>
<p class=bs>
Zur Berechnung eignet sich besonders die Verwendung des Julianischen
Datums <i>JD</i> <a href="#(5)">(5)</a> für den aktuellen Zeitpunkt. Für die
Ableitung von <i>T</i> ist es ausreichend, die UT-Zeit der Atomuhr einzusetzen.
Strenggenommen muß das Zeitmaß UT1 in die Formeln <a href="#(5)">(5)</a> Eingang
finden.
<p class=bs>
Die Eingabe der UT-Zeit wird in ihre Bestandteile zerlegt:
<ul>
<li><i>h</i> = UT-Zeit in dezimalen Stunden
<li><i>d</i> = Tag
<li><i>m</i> = Monat
<li><i>j</i> = Jahr
</ul>
<p class=bs>
Falls Umwandlungen für die Monate Januar (<i>m</i>=1) oder Februar (<i>m</i>=2)
berechnet werden sollen, müssen folgende Zuordnungen für die Monatszahl <i>m</i>
und Jahreszahl <i>j</i> geändert werden (Problematik: Schaltjahre):
<ul>
<li><i>m</i> := <i>m</i> + 12
<li><i>j</i> := <i>j</i> - 1
</ul>
<p class=bs>
Dann gilt für das Julianische Datum <i>JD</i>:
<p>
<a name="(5)"><table width="100%" border=0>
<tr>
<td class=l><i>A</i> = int(<i>j</i>/100<br>
<i>B</i> = 2 - <i>A</i> + int(<i>A</i>/4)<br>
<i>C</i> = <i>h</i>/24<br>
<i>JD</i> = int(365.25<i>j</i>) + int(30.6001(<i>m</i>+1) + <i>d</i>
+ 1720994.5 + <i>B</i> + <i>C</i></td>
<td class=r id=mehrfach>(5)</td>
</table></a>
<p>
<hr>
<a name="LMST"><h2>LMST / LAST</h2></a>
<p class=bs>
Die lokalen Sternzeiten LMST und LAST sind analog zu den Greenwicher Sternzeiten
definiert. Lediglich der Bezugsmeridian ist verschieden, im Fall der lokalen
Sternzeiten bildet der Ortsmeridian mit der geografischen Länge <i>lambda</i>
die Referenz für den Stundenwinkel.
<p class=bs>
Die Umrechnung zwischen lokalen und Greenwicher Sternzeiten erfolgt unter
sinngemäßer Anwendung von <a href="#(2)">(2)</a>.
<p>
<hr>
<a name="Rechnung"><h2>Berechnungen</h2></a>
<ul>
<li>Konversion der Zeitsysteme
<a href="http://habicht.bauv.unibw-muenchen.de/ex-scherer/ast1-UTST.html">UT
&lt;-&gt; Sternzeit</a>
<li>Berechnung
<a href="http://habicht.bauv.unibw-muenchen.de/ex-scherer/ast1-JD.html">
Julianisches Datum <i>JD</i></a>
</ul>
<hr>
<address>Institut für Geodäsie | Universität der Bundeswehr München</address>
Bernd Scherer <b>a91alvub@unibw-muenchen.de</b> [01-Jul-1998]

52
info/ast1-JD.html Normal file
View File

@ -0,0 +1,52 @@
<title>Umwandlung Datum in JD</title>
<h1>Umwandlung Datum in JD</h1>
<hr>
<form action="http://habicht.bauv.unibw-muenchen.de/ex-scherer/ast2-JD.html"
method=POST>
<dl>
<dt><b>Eingabewerte:</b>
<p>
<dd>
<b>UT</b> <input name="uth" size=2 maxlength=2 value="13"> h
<input name="utm" size=2 maxlength=2 value="24"> m
<input name="uts" size=6 maxlength=6 value="50.000"> s (<font size=-1>UT1, nicht UTC</font>)
<dd>
<b>Datum</b> <input name="dd" size=2 maxlength=2 value="18"> .
<input name="dm" size=2 maxlength=2 value="06"> .
<input name="dj" size=4 maxlength=4 value="1999">
<p>
<dd><b>Julianisches Datum</b> <i>JD</i> = gesucht
</dl>
<hr>
<input type="submit" value="Berechnen">
<input type="reset" value="Rücksetzen">
</form>
<hr>
<a href="http://habicht.bauv.unibw-muenchen.de/ex-scherer/scherer.html"
onmouseover="window.status='HomePage Scherer'"
onmouseout="window.status=''">
<img border=0
alt="HomePage Scherer"
src=http://habicht.bauv.unibw-muenchen.de/gif/Home-scherer.gif
align=middle></a>
Ende der Berechnungen, zurück zu meiner
<a href="http://habicht.bauv.unibw-muenchen.de/ex-scherer/scherer.html">
HomePage</a>
<hr>
Aufruf Nummer <b>
0002773
</b> seit dem <b>
10-May-1995
</b>.
<hr>
<address>Institut für Geodäsie | Universität der Bundeswehr München</address>
Bernd Scherer <b>a91alvub@unibw-muenchen.de</b> [20-Oct-1997]

66
info/ast1-UTST.html Normal file
View File

@ -0,0 +1,66 @@
<title>Umwandlung UT in Sternzeit</title>
<h1>Umwandlung UT in Sternzeit</h1>
<dl>
<dt><b>Sie geben ein</b>:
<dd><i>UT</i> (Universal Time)
(<font size=-1>UT1, nicht UTC: |UT1-UTC|&lt;0.7s</font>)
<dd>[evtl] gewünschtes Datum
<dd>[evtl] geographische Länge <i>L</i> des Ortsmeridians
<p>
<dt><b>Sie erhalten als Ergebnis</b>:
<dd>mittlere Sternzeit: <i>GMST</i> (Greenwich), <i>LMST</i> (lokal)
<dd>Equation of Equinoxes <i>EQE</i> zur Ermittlung der entsprechenden
scheinbaren Sternzeit
<dd>scheinbare Sternzeit: <i>GAST</i> (Greenwich), <i>LAST</i> (lokal)
<dd>den Zeiten entsprechende Stundenwinkel
</dl>
<hr>
<form action="http://habicht.bauv.unibw-muenchen.de/ex-scherer/ast2-UTST.html"
method=POST>
<dl>
<dt>
<b>Eingabewerte:</b> (Jahreszahlen gültig ab <u>1985</u>)
<p>
<dd>
<b>UT</b> <input name="uth" size=2 maxlength=2 value="13"> h
<input name="utm" size=2 maxlength=2 value="23"> m
<input name="uts" size=6 maxlength=6 value="15.000"> s (<font size=-1>UT1, nicht UTC</font>)
<dd>
<b>Datum</b> <input name="dd" size=2 maxlength=2 value="18"> .
<input name="dm" size=2 maxlength=2 value="06"> .
<input name="dj" size=4 maxlength=4 value="1999">
<dd>
<b>L</b> <input name="lg" size=3 maxlength=3 value="000"> °
<input name="lm" size=2 maxlength=2 value="00"> '
<input name="ls" size=7 maxlength=7 value="00.0000"> "
<input name="lx" size=1 maxlength=1 value="E">,
geographische Länge Ortsmeridian
</dl>
<hr>
<input type="submit" value="Berechnen">
<input type="reset" value="Rücksetzen">
</form>
<hr>
<a href="http://habicht.bauv.unibw-muenchen.de/ex-scherer/scherer.html"
onmouseover="window.status='HomePage Scherer'"
onmouseout="window.status=''">
<img border=0
alt="HomePage Scherer"
src=http://habicht.bauv.unibw-muenchen.de/gif/Home-scherer.gif
align=middle></a>
Ende der Berechnungen, zurück zu meiner
<a href="http://habicht.bauv.unibw-muenchen.de/ex-scherer/scherer.html">
HomePage</a>
<hr>
<address>Institut für Geodäsie | Universität der Bundeswehr München</address>
Bernd Scherer <b>a91alvub@unibw-muenchen.de</b> [20-Oct-1997]

35
info/ast2-JD.html Normal file
View File

@ -0,0 +1,35 @@
<title>Umwandlung Datum in JD</title>
<h1>Umwandlung Datum in JD</h1>
<hr>
<form action="http://habicht.bauv.unibw-muenchen.de/ex-scherer/ast2-JD.html"
method=POST>
<dl><dt><b>Eingabewerte:</b><p>
<dd><b>UT</b> <input name="uth" size=2 maxlength=2 value="13"> h
<input name="utm" size=2 maxlength=2 value="24"> m
<input name="uts" size=6 maxlength=6 value="50.000"> s (<font size=-1>UT1, nicht UTC</font>)
<dd><b>Datum</b> <input name="dt" size=2 maxlength=2 value="18"> . <input name="dm" size=2 maxlength=2 value="06"> . <input name="dj" size=4 maxlength=4 value="1999">, <tt>Freitag </tt>
<p>
<dd><b>Julianisches Datum</b> <i>JD</i> =<tt> 2451348.05891204</tt>
<img src=http://habicht.bauv.unibw-muenchen.de/gif/k-gruen.gif>
</dl><hr>
<input type="submit" value="Berechnen">
</form>
<hr>
<a href="http://habicht.bauv.unibw-muenchen.de/ex-scherer/scherer.html" onmouseover="window.status= 'HomePage Scherer'" onmouseout="window.status=' '">
<img border=0 alt="HomePage Scherer" src=http://habicht.bauv.unibw-muenchen.de/gif/Home-scherer.gif align=middle></a>
Ende der Berechnungen, zurück zu meiner
<a href="http://habicht.bauv.unibw-muenchen.de/ex-scherer/scherer.html">
HomePage</a>
<hr>
Aufruf Nummer <b>
0002774
</b> seit dem <b>
10-May-1995
</b>.
<hr>
<address>Institut für Geodäsie | Universität der Bundeswehr München</address>
Bernd Scherer <b>a91alvub@unibw-muenchen.de</b> [
18-Jun-1999
]

47
info/ast2-UTST.html Normal file
View File

@ -0,0 +1,47 @@
<title>Ergebnis der Zeitumwandlung</title>
<center><h1>Ergebnis der Zeitumwandlung UT &gt; Sternzeit</h1></center>
<b>Eingabewerte:</b><p>
<table width="70%" border=1>
<td align=center><i>UT</i></td><td align=left><b>13h23m15.000s</b> (<font size=-1>UT1, nicht UTC</font>)<br></td>
<tr><td align=center>Datum</td><td align=left><b>18-Jun-1999</b><br></td>
<tr><td align=center>Jul. Datum <i>JD</i></td><td align=left><b><font color="#788c14"> 2451348.05781250</font></b><br></td>
<tr><td align=center>Ortsmeridian</td><td align=left><b> 0°00' 0.0000"E</b></td></table>
<hr>
<b>Ausgabewerte:</b><p>
<table width="70%" border=1>
<td align=right>mittlere Sternzeit Greenwich:</td><td align=center><b>07h08m37.817s</b></td><td align=center><i>GMST</i><br></td>
<tr><td align=right>mittlere Sternzeit lokal:</td><td align=center><b>07h08m37.817s</b></td><td align=center><i>LMST</i><br></td>
<tr><td align=right>Equation of Equinoxes:</td><td align=center><b><font color="#ffff00">-0.739s</font></b></td><td align=center><a href="#EQE"><i>EQE</i></a><br></td>
<tr><td align=right>wahre Sternzeit Greenwich:</td><td align=center><b>07h08m37.078s</b></td><td align=center><i>GAST</i><br></td>
<tr><td align=right>wahre Sternzeit lokal:</td><td align=center><b>07h08m37.078s</b></td><td align=center><i>LAST</i></td></table>
<dl><dt><a name="EQE">Was bedeutet <i>EQE</i>?</a>
<dd>Übergang vom System der mittleren Sternzeiten zu dem der wahren (scheinbaren) Sternzeiten
</dl><hr>
<b>Entsprechende Stundenwinkel:</b><p>
<table width="35%" border=1>
<td align=center><font color="#ff0000"><i>GMST</i></font></td><td align=center><b>107°09'27.2525"</b><br></td>
<tr><td align=center><font color=#ff0000"><i>GAST</i></font></td><td align=center><b>107°09'16.1652"</b><br></td>
<tr><td align=center><font color="#ff0000"><i>LMST</i></font></td><td align=center><b>107°09'27.2525"</b><br></td>
<tr><td align=center><font color="ff0000"><i>LAST</i></font></td><td align=center><b>107°09'16.1652"</b></td></table>
<hr>
<a href="http://habicht.bauv.unibw-muenchen.de/ex-scherer/ast1-UTST.html"
onmouseover="window.status='erneute Berechnung'"
onmouseout="window.status=''">
<img border=0
alt="erneute Berechnung"
src=http://habicht.bauv.unibw-muenchen.de/gif/pfeil.gif></a>
erneute
<a href="http://habicht.bauv.unibw-muenchen.de/ex-scherer/ast1-UTST.html">
Berechnung</a> (zurück zum Haupttext)
<hr>
Aufruf Nummer <b>
0002123
</b> seit dem <b>
10-May-1995
</b>.
<hr>
<address>Institut für Geodäsie | Universität der Bundeswehr München</address>
Bernd Scherer <b>a91alvub@unibw-muenchen.de</b> [20-Oct-1997]

Binary file not shown.

After

Width:  |  Height:  |  Size: 283 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 870 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 113 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 323 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 739 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 686 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 423 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 335 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 324 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 249 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 303 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 126 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 143 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 200 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 623 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 113 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 113 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 551 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 299 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 341 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 286 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 684 B

Binary file not shown.

After

Width:  |  Height:  |  Size: 161 B

View File

@ -0,0 +1,280 @@
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="Microsoft Word 97">
<META NAME="GENERATOR" CONTENT="Mozilla/4.01 [de] (Win95; I) [Netscape]">
<TITLE>Universal- und Siderial Time</TITLE>
</HEAD>
<BODY TEXT="#000000" BGCOLOR="#DCF3DF" LINK="#7C333F" VLINK="#7C333F" ALINK="#FF0000">
<CENTER><FONT SIZE=+4>U</FONT><FONT SIZE=+3>niversal- und </FONT><FONT SIZE=+4>S</FONT><FONT SIZE=+3>iderial
</FONT><FONT SIZE=+4>T</FONT><FONT SIZE=+3>ime</FONT></CENTER>
<HR>
<P><B><FONT SIZE=+1>1.Aufgabe:</FONT></B>
<P>Zur Umrechnung der Universal Time in das Julianische Datum ist ein Programm
zu schreiben, dabei soll eine Berechnung f&uuml;r alle Tage nach dem 01.03.1600
m&ouml;glich sein.
<PRE>/*****************************************************************
**&nbsp; Titel:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Julianisches Datum&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; **
**&nbsp; Programmiersprache: Ansi-C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; **
**&nbsp; Autor:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Michael Gendt&nbsp; 171766&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; **
**&nbsp; Quellcode:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; julian.c&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; **
**&nbsp; Erstellungsdatum:&nbsp;&nbsp; 26.06.1998&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; **
**&nbsp; Letzte Aenderung:&nbsp;&nbsp; 26.06.1998&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; **
**&nbsp; Grund der Programmierung: Uebung 6 - Erdmessung I&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; **
*****************************************************************/
#include "stdio.h"
#include "stdlib.h"
#include "math.h"
typedef struct
{&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; int day,mon,year,h,min;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; double sec;
}&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; UT;
double UT_to_JD(UT time);
double UT_to_JD2(UT time);
void main()
{
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; UT time_UT={1,1,1600,12,00,00.0};
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; double time_JD;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time_JD=UT_to_JD(time_UT);
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; printf("\n\ntime_UT = %02i.%02i.%02i %02i.%02i.%06.3lf Uhr",time_UT.day,time_UT.mon
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ,time_UT.year,time_UT.h,time_UT.min,time_UT.sec);
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; printf("\ntime_JD = %13.3lf",time_JD);
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time_JD=UT_to_JD2(time_UT);
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; printf("\ntime_JD = %13.3lf",time_JD);
}
double UT_to_JD(UT time)
{
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; int mon[]={0,0,31,59,90,120,151,181,212,243,273,304,334},y_d,sj;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; double time_JD=2305447.0; /*** 00.01.1600 12.00.00 Uhr ***/
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; y_d=time.year-1600;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (y_d&lt;0) {printf("Berechnung nicht m”glich year&lt;1600");exit(0);}
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time_JD += time.day + mon[time.mon] + y_d*365.0;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time_JD += floor(y_d/4) - floor(y_d/100) + floor(y_d/400) + 1;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if ((int)fmod(y_d,4)==0) sj=1;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if ((int)fmod(y_d,100)==0) sj=0;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if ((int)fmod(y_d,400)==0) sj=1;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if ((time.mon&lt;3)&amp;&amp;(sj==1)) time_JD--;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time_JD += (time.h+time.min/60.+time.sec/3600.-12.)/24;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; return time_JD;
}</PRE>
Nutz man die Eigenschaft aus, da&szlig; man nach Teilung eines Integerwertes
auf jeden Fall ein Integer-Wert erh&auml;lt (Nachkommastellen werden abgehackt),
kann man folgende Funktion verwenden. Zudem wurde eine andere M&ouml;glichkeit
der Abfrage eines m&ouml;glichen Schaltjahres genutzt (schnellere Rechenzeit).
<PRE>double UT_to_JD2(UT time)
{
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; int mon[]={0,0,31,59,90,120,151,181,212,243,273,304,334},y_d;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; double time_JD=2305447.0; /*** 00.01.1600 12.00.00 Uhr ***/
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; y_d=time.year-1600;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (y_d&lt;0)&nbsp;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; {&nbsp; printf("Berechnung nicht moeglich year&lt;1600");
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; exit(0);
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; }
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time_JD += time.day + mon[time.mon]+ y_d*365.0;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time_JD += y_d/4 - y_d/100 + y_d/400 + 1;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if ((time.mon&lt;3)&amp;&amp;((int)fmod(y_d,4)==0)&amp;&amp;((y_d/100-4*y_d/400)!=0))
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time_JD--; /*** Abfrage Schaltjahr***/
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time_JD += (time.h+time.min/60.+time.sec/3600.-12.)/24;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; return time_JD;
}
Ergebnis: time_JD = 2445889,906</PRE>
<B><FONT SIZE=+1>2.Aufgabe:</FONT></B>
<P>F&uuml;r einen Punkt mit der L&auml;nge&nbsp;<IMG SRC="Image49.gif" HEIGHT=18 WIDTH=13>=13&deg;0537.000
und der Zeit 25.06.1998 23<SUP>h</SUP>06<SUP>min</SUP>35<SUP>s</SUP>.000
MESZ (Mitteleurop&auml;innsche Sommerzeit) soll die scheinbare (wahre)
Sternzeit (LAST) berechnet werden.
<P>F&uuml;r die Berechnung der Local Apperent Sidereal Time (LAST) l&auml;&szlig;t
sich folgende Formel angeben:
<UL>
<TABLE BORDER=0 COLS=1 WIDTH="165" BGCOLOR="#000000" >
<TR><TD>
<TABLE BORDER=0 COLS=1 WIDTH="163" BGCOLOR="#FFBB44" >
<TR><TD>
<IMG SRC="Image50.gif" HEIGHT=18 WIDTH=161>
</TD></TR></TABLE></TD></TR></TABLE>
</UL>
dabei sind:
<UL>
<TABLE BORDER=0 COLS=1 WIDTH="484" BGCOLOR="#000000" >
<TR><TD>
<TABLE BORDER=0 COLS=1 WIDTH="482" BGCOLOR="#FFBB44" >
<TR><TD>
<IMG SRC="Image51.gif" HEIGHT=21 WIDTH=480>
</TD></TR></TABLE></TD></TR></TABLE>
<BR>&nbsp;
<UL>Umrechnung der Universal Time in die Greenwich Mean Sidereal Time (GAST
= scheinbare Sternzeit von Greenwich) f&uuml;r die ganzen Tage
<P><IMG SRC="Image52.gif" HEIGHT=25 WIDTH=348></UL>
&nbsp;
<P>
<TABLE BORDER=0 COLS=1 WIDTH="200" BGCOLOR="#000000" >
<TR><TD><TABLE BORDER=0 COLS=1 WIDTH="198" BGCOLOR="#FFBB44" >
<TR><TD>
<IMG SRC="Image53.gif" HEIGHT=24 WIDTH=196>
</TD></TR></TABLE></TD></TR></TABLE>
</UL>
<DIR>
<DIR>Umrechnung der Tagesbruchteile von UT nach GAST, dabei wird ber&uuml;cksichtigt,
da&szlig; ein Sterntag rund 4<SUP>min</SUP> k&uuml;rzer als ein Sonnentag
ist. Dies beruht auf der Tatsache, da&szlig; sich die Erde pro Tag (360/365)&deg;
auf ihrer Bahn fortbewegt (Rotation um die Sonne).
<P>&nbsp;</DIR>
<TABLE BORDER=0 COLS=1 WIDTH="156" BGCOLOR="#000000" >
<TR><TD><TABLE BORDER=0 COLS=1 WIDTH="154" BGCOLOR="#FFBB44" >
<TR><TD>
<IMG SRC="Image54.gif" HEIGHT=22 WIDTH=152>
</TD></TR></TABLE></TD></TR></TABLE>
</DIR>
<DIR>
<DIR>Ber&uuml;cksichtigung der Nutation (Berechnung der Winkel siehe &Uuml;bung
4 - Nutation und Pr&auml;zession)
<BR>&nbsp;</DIR>
<P>
<TABLE BORDER=0 COLS=1 WIDTH="145" BGCOLOR="#000000" >
<TR><TD><TABLE BORDER=0 COLS=1 WIDTH="143" BGCOLOR="#FFBB44" >
<TR><TD>
<IMG SRC="Image55.gif" HEIGHT=21 WIDTH=141>
</TD></TR></TABLE></TD></TR></TABLE>
</DIR>
<DIR>
<DIR>Transformation von GAST in die lokale Zeit (LAST)
<BR>&nbsp;</DIR>
</DIR>
Zun&auml;chst m&uuml;ssen wir beachten, da&szlig; uns die Zeit nicht in
UT gegeben ist, sondern f&uuml;r die MESZ, es sollte deshalb die Transformation
<UL>UTC = MESZ - 2h</UL>
durchgef&uuml;hrt werden. In unserem Fall erhalten wir:
<UL>UTC = 25.06.1998 21<SUP>h</SUP>06<SUP>min</SUP>35<SUP>s</SUP>.000</UL>
F&uuml;r die Berechnung von A wird als erstes das julianische Datum von&nbsp;<IMG SRC="Image56.gif" HEIGHT=25 WIDTH=82 ALIGN=ABSCENTER>
<BR>ben&ouml;tigt.
<UL><IMG SRC="Image57.gif" HEIGHT=25 WIDTH=113 ALIGN=ABSBOTTOM>= 2450989,5</UL>
Damit ergibt sich f&uuml;r&nbsp;<IMG SRC="Image58.gif" HEIGHT=18 WIDTH=20 ALIGN=ABSCENTER>
<UL><IMG SRC="Image113.gif" HEIGHT=21 WIDTH=124></UL>
und f&uuml;r A folgt
<UL>A = 18<SUP>h</SUP>11<SUP>min</SUP>44<SUP>s</SUP>,0416.</UL>
F&uuml;r B erhalten wir
<UL>B = 1,00273790935079 <SUP>.</SUP> 21<SUP>h</SUP>06<SUP>min</SUP>35<SUP>s</SUP>.000
= 21<SUP>h</SUP>10<SUP>min</SUP>03<SUP>s</SUP>.0674</UL>
Daraus ergibt sich die Greenwich mean sidereal time
<UL>GMST = A+B = 15<SUP>h</SUP>21<SUP>min</SUP>47<SUP>s</SUP>.1090</UL>
Im Anschlu&szlig; k&ouml;nnen wir den Einflu&szlig; der Nutation berechnen
und auf die GMST anbringen. Die Nutationswinkel&nbsp;<IMG SRC="Image60.gif" HEIGHT=21 WIDTH=26 ALIGN=ABSCENTER>
und&nbsp;<IMG SRC="Image61.gif" HEIGHT=22 WIDTH=64 ALIGN=ABSCENTER>lassen
sich wie in &Uuml;bung 4 gezeigt in Abh&auml;ngigkeit von der Zeit und
mit Hilfe von Koeffizienten berechnen. F&uuml;r die Zeit t nehmen wir die
tats&auml;chliche Zeitdifferenz zur Basisepoche an, also nicht die Differenz
zu 0.00 Uhr.
<P>F&uuml;r unseren Fall erhalten wir:
<DIR><IMG SRC="Image62.gif" HEIGHT=24 WIDTH=330></DIR>
und damit
<UL><IMG SRC="Image114.gif" HEIGHT=45 WIDTH=364></UL>
Daraus folgt f&uuml;r EoE:
<UL>EoE = -00<SUP>h</SUP>00<SUP>min</SUP>00<SUP>s</SUP>.4393</UL>
Durch Addition mit der GMST erhalten wir die GAST (Greenwich Apparent Sidereal
Time)
<UL>
<TABLE BORDER=0 COLS=1 WIDTH="302" BGCOLOR="#000000" >
<TR><TD><TABLE BORDER=0 COLS=1 WIDTH="300" BGCOLOR="#FFDD44" >
<TR><TD>
GAST = GMST + EoE = 15<SUP>h</SUP>21<SUP>min</SUP>46<SUP>s</SUP>.6697
</TD></TR></TABLE></TD></TR></TABLE>
,</UL>
die um die L&auml;nge gedreht werden mu&szlig;, um LAST zu erhalten.
<UL>
LAST = GAST +&nbsp;<IMG SRC="Image64.gif" HEIGHT=18 WIDTH=13> = 15<SUP>h</SUP>21<SUP>min</SUP>46<SUP>s</SUP>.6697
+ 00<SUP>h</SUP>52<SUP>min</SUP>22<SUP>s</SUP>.4667<U></U>
<P>
<TABLE BORDER=0 COLS=1 WIDTH="202" BGCOLOR="#000000" >
<TR><TD><TABLE BORDER=0 COLS=1 WIDTH="200" BGCOLOR="#FFDD44" >
<TR><TD>
LAST = 16<SUP>h</SUP>14<SUP>min</SUP>09<SUP>s</SUP>.136
</TD></TR></TABLE></TD></TR></TABLE></UL>
<P><B><FONT SIZE=+1>3.Aufgabe:</FONT></B>
<P>F&uuml;r einen Punkt mit der L&auml;nge&nbsp;<IMG SRC="Image65.gif" HEIGHT=18 WIDTH=13>=13&deg;2050.000
und der scheinbaren Sternzeit 16.07.1998 12<SUP>h</SUP>04<SUP>min</SUP>37<SUP>s</SUP>.739
soll die MESZ berechnet werden. Die Berechnung ist die Umkehrung der Aufgabe
2 gezeigten Transformation, sie lautet:
<UL>
<TABLE BORDER=0 COLS=1 WIDTH="333" BGCOLOR="#000000" >
<TR><TD><TABLE BORDER=0 COLS=1 WIDTH="331" BGCOLOR="#FFBB44" >
<TR><TD>
<IMG SRC="Image66.gif" HEIGHT=21 WIDTH=329>
</TD></TR></TABLE></TD></TR></TABLE>
</UL>
Die julianischen Daten f&uuml;r die ben&ouml;tigten Zeiten sind:
<UL>UT(16.07.1998 12<SUP>h</SUP>04<SUP>min</SUP>37<SUP>s</SUP>.739) = 2451011,00321
<BR>UT(16.07.1998 00<SUP>h</SUP>00<SUP>min</SUP>00<SUP>s</SUP>.000) = 2451010,50000
<BR>UT(01.01.2000 00<SUP>h</SUP>00<SUP>min</SUP>00<SUP>s</SUP>.000) = 2451545,00000</UL>
In unserem Fall lauten die einzelnen Gr&ouml;&szlig;en:
<UL>LAST = 12<SUP>h</SUP>04<SUP>min</SUP>37<SUP>s</SUP>.739 (das Datum
mu&szlig; hierbei weggelassen werden)
<P><IMG SRC="Image64.gif" HEIGHT=18 WIDTH=13> = 13&deg;2050.000 = 00<SUP>h</SUP>53<SUP>min</SUP>23<SUP>s</SUP>.3333
<P>EoE = -00<SUP>h</SUP>00<SUP>min</SUP>00<SUP>s</SUP>.4351</UL>
<DIR>
<DIR>mit t = -0,014633812,<IMG SRC="Image67.gif" HEIGHT=21 WIDTH=141 ALIGN=ABSCENTER>
und&nbsp;<IMG SRC="Image68.gif" HEIGHT=22 WIDTH=158 ALIGN=ABSCENTER></DIR>
<P>A = -19<SUP>h</SUP>34<SUP>min</SUP>31<SUP>s</SUP>.7043
<DIR>
<P>mit&nbsp;<IMG SRC="Image69.gif" HEIGHT=21 WIDTH=125></DIR>
</DIR>
F&uuml;r Tagesbruchteil am 16.07.1998 erhalten wir folglich:
<UL><IMG SRC="Image70.gif" HEIGHT=21 WIDTH=414></UL>
Da wir nicht die Zeit UTC suchen, m&uuml;ssen wir jetzt noch 2h addieren
um die MESZ zu erhalten.
<UL>
<TABLE BORDER=0 COLS=1 WIDTH="180" BGCOLOR="#000000" >
<TR><TD><TABLE BORDER=0 COLS=1 WIDTH="178" BGCOLOR="#FFBB44" >
<TR><TD>
MESZ = 17<SUP>h</SUP>34<SUP>min</SUP>09<SUP>s</SUP>.678
</TD></TR></TABLE></TD></TR></TABLE>
</UL>
<P>(ohne Gew&auml;hr auf Richtigkeit)&nbsp;
<HR><B><FONT SIZE=-1>Stand:&nbsp;<SCRIPT language="JavaScript">
<!--
document.write(document.lastModified);
//-->
</SCRIPT>
&copy; <FONT FACE="cosmic,arial">(<A HREF="http://www.fpk.tu-berlin.de/~max/">Michael
Gendt</A>,Frank Oheim), Hinweise und Erg&auml;nzungen an <A HREF="mailto:max@fpk.tu-berlin.de">mg</A>
oder&nbsp;<A HREF="http://www.fpk.tu-berlin.de/~fsgeod/grafik/anre.htm"><IMG SRC="button2.gif" ALT="weiter" BORDER=0 ></A></FONT></FONT></B>
<BR>&nbsp;
</BODY>
</HTML>

View File

@ -48,10 +48,12 @@ $BIN = "$SETIDIR/setiathome";
@weekdays = ("Mon","Tue","Wed","Thu","Fri","Sat","Sun");
@months = ("Jan","Feb","Mar","Apr","May","Jun","Jul","Aug","Sep","Oct","Nov","Dec");
open (GPID, "pidof $BIN|");
#open (GPID, "pidof $BIN|");
open (GPID, "ps ax|grep setiathome|grep -v grep|cut -c 1-6|");
$pid = <GPID>;
close GPID;
chomp $pid;
$pid = int $pid;
if ($pid) {
open (PINF, "ps lw $pid|");
@ -60,7 +62,9 @@ if ($pid) {
close PINF;
chomp $pinfh, $pinfd;
$nice=substr $pinfd, (index $pinfh, " NI"), 3;
$nice=int $nice;
$memu=substr $pinfd, (index $pinfh, " SIZE"), 6;
$memu=int $memu;
$cmdl=substr $pinfd, (index $pinfh, "COMMAND"), (length $pinfd)-(index $pinfh, "COMMAND");
}
@ -172,7 +176,7 @@ sub StateInfo {
if ($memu<1000) {
$memut = "$memu kB";
} else {
$memut = substr ($memu,0,3) . "." . substr ($memu,3,3) . " kB";
$memut = substr ($memu,0,(length $memu)-3) . "." . substr ($memu,(length $memu)-3,3) . " kB";
}
$CS = "$CS (mem usage: $memut)";
}