Sonnenaufgang

[Daniel K. (danielk91)]

Hallo User,

hoffentlich kann mir hier jemand weiter helfen. Ich versuche den
Zeitpunkt des Sonnenaufgangs und des Sonnenuntergangs zu ermitteln.
Leider berechnet mir mein ATMega 32 nicht den erwuenschten Wert. Hab es
manuell nachgerechnet und komme auf in etwa eine Stunde Unterschied zu
dem Ergebnis des Mikrocontrollers.

Sieht jemand auf Anhieb den Fehler. Wuerde mich freuen wenn mir hier
jemand weiterhelfen kann! Alle Variablen sind als float deklariert.

Gruss
Daniel



void sunset(void)
{
  unsigned int dd;

  b=3.141592654*(-26.4)/180;
  hh=-0.1145;

  dekli=0.4095*sin(0.016906*(330-80.086));

  timedifferenz=12*acos((sin(hh)-sin(b)*sin(dekli))/(cos(b)*cos(dekli)))/3 
.141592654;

  realtime=12-timedifferenz;

  sunrisetime=realtime+(15-28.08)*4/60;

[Klaus R. (klaus2)]

...gib doch einfach mal zwischenwerte per uart aus?!?!

mit soviel math.h knallt es bestimmt irgendwo.

Klaus.

[Benjamin S. (recycler)]

Ich würde auch auf die math.h setzen. Berechne lieber die meisten Werte 
im Vorraus und nur das nötigste in der Routine. Dann gibts weniger 
Fehler.

Mit welchem Rechner rechnest du. Hast du nicht vergessen auf RAD 
umzustellen, und nicht mit DEG zu rechnen?

[Marvin M. (Gast)]

Da sind einige Integer-Konstanten drin - schreib die mal als Float, also 
z.B. anstelle 180 eine 180.0
Könnte sein, dass der Compiler das nicht richtig casted.

[Daniel K. (danielk91)]

@ Benjamin  Wo stell ich den Rechner auf Rad um. Ich dachte math.h 
rechnet schon in Rad. Weiterhin kann ich nicht viel vorher berechnen, da 
die Tage (hier sind es 330Tage) sich ja aendern und dann alles wieder 
neu berechnet werden muss.

Die Umschreibung in float hat leider nichts bewirkt.

Danke schon mal!!!!

Gruss
Daniel

[Daniel C. (cagara)]

timedifferenz kannst du zusammenfassen in

 ArcCos( Sec(b) Sec(dekli) Sin(hh) - Tan(b) Tan(dekli)) *  7.63944


(falls nicht bekannt, Sec ist der Sekans = 1/cos)

[Johann L. (gjlayde)]

Du weisst ja, daß hier nur mit float gerechnet wird.

Evtl. findest du ne Formulierung, die besser konditioniert ist.

http://de.wikipedia.org/wiki/Gleitkommazahl#Ausl.C3.B6schung
http://de.wikipedia.org/wiki/Fehlerschranke

[Johann L. (gjlayde)]

Die Deklination ist offenbar konstant --> Ausrechnen und Konstante ins 
Programm eintragen. Rumrechnen vergrößert nur die Ungenauigkeit. 
Gleiches gilt für ihre sin- und cos-Werte und 12/pi

d.h. anstatt

$$\Delta t=\frac{12}{\pi}\operatorname{acos}\left(\frac{\sin(hh)-\sin(b) \cdot\sin(\delta)}{\cos(b)\cdot\cos(\delta)}\right)$$

Berechnest du

$$\Delta t=K_1\cdot\operatorname{acos}\left(\frac{\sin(hh)-K_2\cdot\sin b }{K_3\cdot\cos b }\right)$$

...und das gilt auch für b

Was an der Formel ist denn überhaupt variabel? Nur hh ? In dem Fall 
wäre zu berechnen

$$\Delta t=K_1\cdot\operatorname{acos}\left(\frac{\sin(hh)-K_4 }{K_5}\right)$$

[Benjamin S. (recycler)]

Die Rechnung mit RAD führt die math.h richtig aus.
Ich meine ob du in deinem Taschenrechner oder in Windowsrechner oder 
auch unter Linux mit kcalc usw. den richtigen Modus für die Verifikation 
eingestellt hast.
Oder liefert nur das Programm fehlerhafte Werte.

[Daniel K. (danielk91)]

Danke an alle. Ich hab das Problem wieder mal mit eurer Hilfe geloesst. 
Ich hab teilweise eure Ratschlaege umgesetzt und die Genauigkeit der 
Werte erhoeht und nun ist die Rechnung Praezise.

Danke und schoenen Tag
Daniel

[Ansgar S. (ansi)]

Moin Daniel

Daniel K. wrote:
> Danke an alle. Ich hab das Problem wieder mal mit eurer Hilfe geloesst.
> Ich hab teilweise eure Ratschlaege umgesetzt und die Genauigkeit der
> Werte erhoeht und nun ist die Rechnung Praezise.

Ich bin auch gerade auf der Suche nach einer entsprechenden 
Berechnungsroutine. Kannst du deine Routine eventuell zur Verfuegung 
stellen oder mir zumindest privat zuschicken?
Tausend Dank.

Gruss Ansgar

[900ss (900ss)]

Ich habe diese Routine mal getestet und sie funktionierte:
http://www.mikrocontroller.net/attachment/30393/sun.c

Gruß 900ss