Robotrontechnik-Forum

Registrieren || Einloggen || Hilfe/FAQ || Suche || Mitglieder || Home || Statistik || Kalender || Admins

Willkommen Gast!

Robotrontechnik-Forum » Technische Diskussionen » Mond- Sonnenberechnungen » Themenansicht

Autor

Thread - Seiten: -1-

000
10.10.2020, 21:22 Uhr
Lippi

Moin Leute...

Ich suche (BASIC) ein Stück Programm, welches z.B. aktuelle Mondphase oder Sonnenauf- untergangszeit berechnet. Ohne Grafik oder Sternhimmel....

Für Sonnenauf- und Untergangszeiten hatte ich mal ausreichend genaue und einfache Formeln früher im Netz gefunden. Das hatte ich zusammen mit einigen anderen Unterprogrammen dann in Arexx für den Amiga umgeschrieben und benutzt mit InfoChannel zusammen zur aktuellen Einblendung dieser Werte im Kabelfernsehen.
Nun habe ich fast alle Amigas verkauft und leider von dem Script nirgends mehr eine Kopie ! Sonst könnte ich das vielleicht in Basic umschreiben.

Aber vielleicht hat einer von Euch ein Stück Programm dafür ?
Ich glaube, für den KC gabs das mal - man gibt das Datum ein und geografische Koordinaten und man bekam die Werte berechnet.

Vielen Dank schonmal

--
MfG Mario

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

001
10.10.2020, 21:26 Uhr
schlaub_01

Hallo Lippi,

hab mal kurz hier geschaut: https://www.iee.et.tu-dresden.de/~kc-club/09/RUBRIK04.HTM
Such mal nach Mond.sss - könnte es das sein?

Viele Grüße,
Sven

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

002
10.10.2020, 21:31 Uhr
Lippi

Das hatte ich gefunden. Das kann die Halb- Voll, Neumondtage berechnen und auch Finsternisse. Das ist leider zu "aufwändig".
Oder ich muss mir die Mühe machen, die entsprechenden Programmstücke dort zu extrahieren .... Und Sonnenaufgang ist leider nicht dabei :-(

Nachtrag: Ist leider auch im Programm nichts dokumentiert :-(
--
MfG Mario

Dieser Beitrag wurde am 10.10.2020 um 21:32 Uhr von Lippi editiert.

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

003
10.10.2020, 21:36 Uhr
schlaub_01

okay, schade. Hab dann noch das hier gefunden: https://www.iee.et.tu-dresden.de/~kc-club/02/KCN98-03/KCN98-03-03.HTML

Grüße,
Sven.

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

004
10.10.2020, 21:42 Uhr
Lippi

Oh, liest sich gut, leider aber nicht BASIC.

Im Archiv habe ich auch schon nach allem, was mit Astro, Sterne usw. passen könnte gesucht...

OK. Ist nicht lebenswichtig - weiss aber nicht, ob ich heute schlafen kann, wenn ich nix finde

--
MfG Mario

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

005
10.10.2020, 22:39 Uhr
Besserwisser

Avatar von Besserwisser

Lies mal das hier:

https://webcache.googleusercontent.com/search?q=cache:OsH-Vkx7H4cJ:https://lexikon.astronomie.info/zeitgleichung/+&cd=2&hl=de&ct=clnk&gl=de

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

006
10.10.2020, 22:59 Uhr
Lippi

Danke, das ist ein Teil, den ich gesucht habe.
Auf die Idee, im Cache zu suchen, bin ich leider nicht gekommen :-( Die ganze Hompage ist leider seit einigen Tagen endgültig abgeschaltet....
--
MfG Mario

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

007
11.10.2020, 02:41 Uhr
Besserwisser

Avatar von Besserwisser

Zitat:

Lippi schrieb
Auf die Idee, im Cache zu suchen, bin ich leider nicht gekommen.

Das Internet vergisst nichts.
Es weiß aber auch nicht alles.

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

008
11.10.2020, 08:55 Uhr
Lippi

Verdammt, Mathe ist schon lange her

Was mache ich denn nun, wenn der KC keine arccos - Funktion hat ? Das finde ich zwar selsam ATAN hat er ja auch, aber halt nicht zu ändern. Arkus ist ja nix anderes, als die Umkehr der Funktion cos, sin usw. Doch wie berechnen in BASIC ?

Sorry für die blöde Frage.....
--
MfG Mario

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

009
11.10.2020, 10:11 Uhr
Besserwisser

Avatar von Besserwisser

https://en.wikipedia.org/w/index.php?title=Inverse_trigonometric_functions§ion=7#Relationships_among_the_inverse_trigonometric_functions

Dieser Beitrag wurde am 11.10.2020 um 10:14 Uhr von Besserwisser editiert.

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

010
11.10.2020, 11:44 Uhr
MS05

Zitat:

Lippi schrieb
Was mache ich denn nun, wenn der KC keine arccos - Funktion hat ? Das finde ich zwar selsam ATAN hat er ja auch, aber halt nicht zu ändern. Arkus ist ja nix anderes, als die Umkehr der Funktion cos, sin usw. Doch wie berechnen in BASIC?

Das ist umrechenbar

arccos x = arctan(wurzel(1-x²)/x)

Dieser Beitrag wurde am 11.10.2020 um 11:49 Uhr von MS05 editiert.

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

011
11.10.2020, 12:29 Uhr
wolle1945

Zitat:

Lippi schrieb
Verdammt, Mathe ist schon lange her

schau mal hier:
https://www.matheretter.de/wiki/sinus-kosinus-berechnen
dann unter "Sinus konkret berechnen (Taylorreihe)"
Es gibt auch noch andere Seiten, wo Formeln dargestellt werden,
die man nun mit Basic lösen kann.
Ich hatte mitte der 80-ziger mal eine Computer-AG mit Schülern
geleitet. Wir hatten sin, cos, tan, ... mit kleinen Programmen in
Basic erstellt. Muß mal auf dem Boden nachschauen, ob noch
etwas zu finden ist.
--
mfG wolle1945

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

012
11.10.2020, 15:10 Uhr
Lippi

Leute, Ihr seid Spitze

Das "arccos x = arctan(wurzel(1-x²)/x)" wars

Habe die anderen Beiträge auch gelesen und versucht zu verstehen - die kurze Formel brachte aber den Sieg.
Beschämend für mich, dass ich sowas nicht weiß. Ich hatte in Mathe ne glatte Eins - leider ist das schon ne ganze Weile her.... und Sachen, die man nicht benötigt, geraten halt irgendwann in Vergessenheit.
Leider hab ich auch durch Gurgel und Co im www nichts darüber gefunden. Erklärungen und Onlinerechner haufenweise, aber halt nicht die Berechnungsgrundlage.

Jetzt bin ich soweit, dass der KC irgendwas berechnet, was plausibel scheint. Fehlt nur noch die Kleinigkeit der Anpassung der Ortszeit an die MEZ / MESZ.

Ganz lieben Dank an alle !

Dann aber gleich die nächste - noch blödere - Frage: Wie zum Geier nutzt man die Exponentialfunktion in BASIC ??

Laut Handbuch EXP(x) - dabei soll x der Exponent sein, sicherlich. Ja bitte und die Zahl, die ich damit rechnen lassen will ?? Ich stehe grade auf dem Schlauch. Print EXP(2) gibt 7.38906 ??

Was mache ich hier wieder fürn Denkfehler ?

(ok. ich habe für x hoch 2 dann einfach x * x gerechnet - klar, das geht :-)
--
MfG Mario

Dieser Beitrag wurde am 11.10.2020 um 15:14 Uhr von Lippi editiert.

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

013
11.10.2020, 15:49 Uhr
volkerp
Default Group and Edit

exp(x) = e^x, Exponentialfunktion

Normale Potenzfunktion ist "^": PRINT 2^3 --> 8

Winkelfunktionen s. Heblik, "Wissenspeicher BASIC", letzte Umschlagseite

arcussinus: ATN(X/SQR(-X*X+1))
arcuscosinus: -ATN(X/SQR(-X*X+1))+1.5708
--
VolkerP

http://hc-ddr.hucki.net
(Z9001, Z1013, LC-80, ...)

Dieser Beitrag wurde am 11.10.2020 um 16:01 Uhr von volkerp editiert.

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

014
11.10.2020, 15:51 Uhr
MS05

Gerne. Ich habe das nach den Jahren auch nicht mehr alles auf der Zunge liegen. Die beste Hilfe ist mir hier meist "Arithmetische Algorithmen der Mikrorechentechnik" von lorke/Lampe/Wengel. Das stand immer in Griffweite.

Nuja... da ist das Basic etwas "speziell".
Gemeint ist hier die Basis des natürlichen Logarithmus, die eulersche Zahl e = 2,71828182....
EXP(X) ist hier e^x

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

015
11.10.2020, 16:19 Uhr
Lippi

Danke Euch.achje, da habe ich das mit EXP verwechselt.
Also, das klappt jetzt.
Nur noch bissel formatieren in h und min......

Macht Spass mit dem alten Infokanalprogramm -)

Wenn das mal jemanden als ganzes interessieren sollte, würde ichs auf meiner HP zum Download anbieten. Sven und Mario kennen es ja schon.....
--
MfG Mario

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

016
11.10.2020, 18:10 Uhr
Lippi

Mist, zu früh gefreut

Vielleicht mache ich noch einen Denkfehler:

Die Formel für die Zeitdifferenz zwischen Meridiandurchgang der Sonne und der (in etwa) Horizonthöhe lautet: (h ist die Höhe Horizont minus 50 Bodenminuten)

Zeitdifferenz = 12 * arccos((sin(h) - sin(Breitengrad) * sin(Deklination_der_Sonne))) / (cos((Breitengrad) * cos(Deklination_der_Sonne)))) / Pi

Um arccos zu ersetzen habe ich nun geschrieben:

80 x = ((sin(h) - sin(Breitengrad) * sin(Deklination_der_Sonne))) / (cos(Breitengrad) * cos(Deklination_der_Sonne)))
90 ZD = 12 * atn(sqr(1-x*x) / x) / Pi

Was mache ich nun wieder falsch ??

Wenn ich das Programm in einer Schleife durchlaufen lasse, scheinen die Werte plausibel.
Das fängt mit Tag 1 an, was der 01.01. ist. Aufgangs- und Untergangszeiten verschieben sich, bis wir zur Frühjahrs Tag- und Nachtgleiche kommen. Dann spinnen die Zeiten total.
Ich habe mir ZD mit ausdrucken lassen und dieser Wert läuft hoch bis 5,98... und springt am nächsten Tag auf fast den gleichen Wert negativ !

Bis etwa zur Herbst-Tag-und Nachtgleiche, dann springts wieder auf positiv.

Da ich die gleiche Berechnung (glaube zumindest) damals in dem AREXX Programm genutzt hatte, nehme ich an, ich mache einen Fehler in der Ersatzrechnung für den arccos - Befehl ?

Kennt Ihr das ? Wenn ein was total eigendlich Unwichtiges nicht funktioniert ?? - Furchtbar

--
MfG Mario

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

017
11.10.2020, 18:14 Uhr
volkerp
Default Group and Edit

Schau mal in meinen Beitrag: du hast die Formel f. Arcsin, nicht für arccos genommen.
--
VolkerP

http://hc-ddr.hucki.net
(Z9001, Z1013, LC-80, ...)

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

018
11.10.2020, 18:25 Uhr
Lippi

Ach, verdammich !!!

Danke Dir.

Ich glaube, ich gehe erstmal ne kurze Pause machen. Mir raucht die Birne ....

--
MfG Mario

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

019
05.12.2020, 11:46 Uhr
Lippi

Moin Leute.

Ich brauche jetzt mal wieder Eure Hilfe....

Und zwar zur Berechnung der Mondphase und der Berechnung der beleuchteten Fläche.

Bis zur verhältnismässigen Berechnung des Mondalters bin ich fertig.
Ich gehe von einem bekannten Vollmondzeitpunkt aus, berechne die Tagesdifferenz und bekomme das Mondalter vom letzten Vollmond in %, was mit wenigen Stunden Differenz stimmt und damit ausreichend genau ist.

Das bedeutet ich habe MA(lter) von 0 bis 100%
0% = Vollmond
25% = abnehmender Halbmond
50% = Neumond
75% = zunehmender Halbmond
100% = Vollmond
nach kurzer Umrechnung ergeben sich die Werte für % in beleuchteter Fläche:
0% = 100%
25% = 50%
50% = 0%
75% = 50%
100% = 100%
jeweils links Mondalter, rechts Mondphase in bel. %

Nun stimmt das leider nur genau bei Voll-, Halb- und Neumond, da der Mond ja eine (sichtbare) Halbkugel ist und die prozentuale Änderung der sichtbaren beleuchteten Fläche in der Nähe von Neu- und Vollmond "langsamer" verläuft, als im Bereich Halbmond.
Oder anders ausgedrückt ist die tägliche Änderung im Bereich Halbmond schneller, als bei Neu- und Vollmond, da der Terminator (Grenze Hell / Dunkel) auf dem Mond ja da zunehmend von der Erde in einem spitzeren Winkel zu sehen ist.

Jetzt habe ich schon tagelang mit sin und cos - Funktionen herumgedoktort, komme aber zu keinem Ergebnis.
Mit cos bekomme ich das hin, dass z.B nach Vollmond die tägliche Änderung des Mondalters z.B. 99%, 94%, 89% usw. ist, die Mondphase korrekt "langsamer weniger wird.
So einfach geht es aber nicht, da bei Halbmond der cos 75 anstatt 50% ergibt.
Im Neumondbereich zeigt es dann sogar zu viel an, bis es bei genau Neumond auch bei 0% ist.

Das geht bestimmt mit cas oder sin recht easy zu berechnen - ich komme nur nicht drauf !!

Im Netz finde ich nur Anleitungen zur Mondalterberechnung, das ist recht einfach und funktioniert bereits ausreichend genau.

Vielleicht steckt von Euch einer (noch) mehr in der Mathematik drin, als ich - die 1 in Mathe ist leider über 45 Jahre her :-)

Danke schonmal...

--
MfG Mario

Dieser Beitrag wurde am 05.12.2020 um 11:47 Uhr von Lippi editiert.

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

020
06.12.2020, 08:08 Uhr
Bert

Zitat:

Lippi schrieb
So einfach geht es aber nicht, da bei Halbmond der cos 75 anstatt 50% ergibt.
Im Neumondbereich zeigt es dann sogar zu viel an, bis es bei genau Neumond auch bei 0% ist.

Hallo Mario,

ich weiß nicht, ob ich Dich da richtig verstehe:
Der Cosinus hat ja als Wertebereich -1 bis +1. Das muß ja sowieso an den 'Wunschbereich Mondphase" von 0 bis 1 bzw. 0% bis 100% angepasst werden.

Entweder man halbiert den Cosinus und addiert 1/2, oder man quadriert.
Beim Quadrieren verdoppelt sich gewissermaßen die Frequenz, was mit dem Faktor 1/2 berücksichtigt wird.

Beide Ergebnisse sind identisch und haben bei einem Mondalter von 0,25 und 0,50 als Ergebnis 50%:

Ist es das, was Du wolltest?

Viele Grüße,
Bert

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

021
06.12.2020, 15:05 Uhr
Lippi

Danke Bert, ich probiere das mit einem Testlauf dann .
--
MfG Mario

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

022
06.12.2020, 19:06 Uhr
Lippi

Funktioniert - nochmals danke !
--
MfG Mario

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

023
15.01.2021, 07:46 Uhr
Lippi

So Leute, das nächste Mathematikproblem

Ich habe jetzt in ein altes Windows-Präsentationsprogramm (SCALA IC) mit einigen Tricks (da die mathematischen Funktionen hier nicht so umfangreich, wie im KC Basic sind und leider die Rechengenaueigkeit auch nicht) die Berechnung für die Mondphase hinbekommen - im KC hatte das ja dann auch perfekt geklappt.

Nun wollte ich Sonnenauf- und Untergang auch noch gern dem Windows-Prpgramm beibringen.
Leider gibts da schon wieder eine Hürde !!

Das Ding hat kein ATN .... im Netzt habe ich leider nichts gefunden.
Wie könnte man die Arctangensfunktion ersetzen, wenn die Software den Befehl nicht hat ?

Sin, Cos, Tan, Quadratwurzel hat er, dann hörts aber schon fast auf.

Vielleicht habt Ihr eine Idee

Nachtrag: Jetzt habe ich es herausgefunden :-)
--
MfG Mario

Dieser Beitrag wurde am 15.01.2021 um 10:23 Uhr von Lippi editiert.

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

024
15.01.2021, 10:36 Uhr
volkerp
Default Group and Edit

s. https://mathepedia.de/Arkustangens_und_Arkuskotangens.html

atn(x) ~ 3*x/(3+x*x) für x < 1

Die Arkustangenswerte über 1 lassen sich aus den Werten zwischen 0 und 1 ableiten:

atn(1/x)=PI/2-atn(x)
--
VolkerP

http://hc-ddr.hucki.net
(Z9001, Z1013, LC-80, ...)

Dieser Beitrag wurde am 15.01.2021 um 10:37 Uhr von volkerp editiert.

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

025
16.01.2021, 09:13 Uhr
Lippi

Guten Morgen

Ich habe folgende Näherung benutzt:

ATN = x - x * x * x / 3 + x * x * x * x * x / 5

Da die komische Software für Windows (ich hatte jahrelang SCALA Infochannel für Amiga im Einsatz: Das konnte "besser" rechnen und hatte eine AREXX - Schnittstelle - die Version für Windows ist dagegen eine Krücke!!) auch ein Minus vor dem Ausdruck nicht akzeptiert, muss danach noch das hier folgen:

ATN = ATN - ATN - ATN
--
MfG Mario

Profil || Private Nachricht || Suche

Zitatantwort || Editieren || Löschen

Seiten: -1- [ Technische Diskussionen ]