Soooooooo zufrieden!!!!!!!!!!!!!!!!!
Soooooooo zufrieden!!!!!!!!!!!!!!!!!
Hallo,
ich habe nach http://lexikon.astronomie.info/zeitgleichung/ ein BASCOM-Programm geschrieben, das allerdings noch nicht ausführlich getestet ist. Als Zeitgeber soll eine DCF-Uhr dienen.
'Sonnentest 13.4.12
'lexikon.astronomie.info/zeitgleichung/
'für Radebeul: 13°37' Ost, 51°06' Nord
Code:$crystal = 1000000 $regfile = "m48def.dat" $hwstack = 32 ' default use 32 for the hardware stack $swstack = 10 ' default use 10 for the SW stack $framesize = 40 Dim Woz_moz As Single , Term1 As Single , Term2 As Single , Term3 As Single Dim Deklination As Single , Zeitdiff As Single , Sa As Single , Su As Single Dim Samin As Single , Sumin As Single Dim Tag As Integer , Diff2 As Integer Dim Sahour As Byte , Samin1 As Byte , Suhour As Byte , Sumin1 As Byte Dim Sastrmin As String * 3 , Sustrmin As String * 3 Anfang: Input "Tagzahl " , Tag Term1 = 0.0337 * Tag 'Berechnung Differenz Wahre + Mittlere Term1 = Term1 + 0.465 'Ortszeit Term1 = Sin(term1) Term1 = -0.171 * Term1 Term2 = 0.01787 * Tag Term2 = Term2 - 0.168 Term2 = Sin(term2) Term2 = 0.1299 * Term2 Woz_moz = Term1 - Term2 'Zeitgleichung /in Stunden Diff2 = Woz_moz * 60 'auf Minuten gerechnet Term1 = Tag - 80.086 'Berechnung Deklination Term1 = Term1 * 0.016906 Term1 = Sin(term1) Deklination = Term1 * 0.4095 Term1 = Sin(deklination) 'Berechnung Zeitdifferenz Term1 = 0.7771 * Term1 'B=0,89 für Radebeul Term1 = -0.0145 - Term1 Term2 = Cos(deklination) Term2 = 0.6294 * Term2 Zeitdiff = Term1 / Term2 Zeitdiff = Acos(zeitdiff) Zeitdiff = 3.8197 * Zeitdiff Sa = 12 - Zeitdiff 'Sonnenaufgang Sa = Sa - Woz_moz Sa = Sa - 0.9067 Sa = Sa + 1 'MEZ, bei MESZ +2 Samin = Frac(sa) Sahour = Sa - Samin Samin1 = Samin * 60 Sastrmin = Str(samin1) Sastrmin = Format(sastrmin , "00") Su = 12 + Zeitdiff Su = Su - Woz_moz Su = Su - 0.9067 Su = Su + 1 Sumin = Frac(su) Suhour = Su - Sumin Sumin1 = Sumin * 60 Sustrmin = Str(sumin1) Sustrmin = Format(sustrmin , "00") Print "WOZ-MOZ: " ; Woz_moz Print "min: " ; Diff2 Print "Deklination: " ; Deklination Print "Zeitdifferenz: " ; Zeitdiff Print "Sonnenaufgang: " ; Sahour ; ":" ; Sastrmin Print "Sonnenuntergang: " ; Suhour ; ":" ; Sustrmin Print Goto Anfang End
Geändert von radbruch (20.08.2012 um 15:54 Uhr) Grund: Code-Tag eingefügt
Hi
Schau mal hier :http://suncalc.net/#/51.1116,13.6494...12.08.20/16:42
Da kannst du deine zahlen überprüfen
achso eine Sache Term1 = Term1 - 1 dann ist es Perfekt
Gruß
Andreas
Geändert von U-1110 (20.08.2012 um 16:03 Uhr)
Hallo, da ich oft term1 nehme, finde ich nicht die Stelle, die Du meinst .
Im Moment komme ich nicht zum Testen, ist aber vorgemerkt.
Dolivo
Hi dolivo
Term1 = 0.0337 * Tag 'Berechnung Differenz Wahre + Mittlere
Term1 = Term1 + 0.465 'Ortszeit
Term1 = Sin(term1)
Term1 = -0.171 * Term1
Term1 = Term1 - 1 'dann stimmt die Uhrzeit(Sommerzeit)
Habs grad mal laufen lassen ,um eine Minute daneben als im Internet
Echt Super!!!!!
Hi Dolivo,
Dein Code sieht ganz chick aus! Nun würde ich diesen gern bei mir in meine Haussteuerung einsetzen. Leider nur kann ich nicht einen Bezugspunkt im Code finden, bei dem ich meine Standortkoordinaten einfügen kann. Könntest Du mir bitte mal kurz unter die Arme greifen und helfen .
Danke und Grüße MAT
Hallo,
ich habe die im lexikon.astronomie aufgeführten Formeln einfach nur "aufgetröselt", damit sie in Bascom verarbeitet werden können. Wenn Du Dir das Zahlenbeispiel von Berlin ansiehst, wird auch manches klarer. Ich hänge aber noch einmal mein Programm an, in dem ich ein paar Erläuterungen mehr eingefügt habe.
Hi,
dank für Deine Antwort!
habe mal versucht den Code im Simulator laufen zu lassen, natürlich mit meinen Koordinaten.
Was ich bei Deinem Code nicht so recht verstehe ist:
'Sonnenaufgang
Sa = 12 - Zeitdiff
Sa = Sa - Woz_moz
Sa = Sa - 0.9067 '13°37'/15°=0,9067
...wie kommst Du auf 0,9067?
Und dann ist da noch... welche Zeit wird genau ausgerechnet? Wenn ich die Ergebnisse mit http://www.suncalc.net/#/51.7236,12....12.10.01/18:37 vergleiche, gibt es mehrere Arten der Sonnen auf/untergänge...
hmmmmm.... werde wohl noch weiter testen.
Hallo, mat-sche,
ja, da habe ich falsch gerechnet. Richtig ist 0,9078.
Ich habe meine Werte bisher nur mit den Zahlen verglichen, die im Kalender stehen und die Zeitdifferenz zwischen meinem Ort und dem angegebenen (hier Kassel, ermittelt mit der Längengraddifferenz) abgezogen (da östlich). Das stimmt recht gut überein innerhalb des genannten Fehlerbereichs.
Nun noch zum Beitrag von U-1110 v. 28.8.: Ich berücksichtige die Sommer- bzw. Normalzeit beim Sonnenauf- bzw. -untergang. An der Stelle interessiert mich ja nur das Zwischenergebnis.
dolivo
Hi Dolivo,
nun endlich habe ich Deinen Code in meine Haussteuerung integriert und nun läuft die Berechnung dank dcf77 um die 1min genau.
Danke für die Codezeilen!
Grüße MAT
Lesezeichen