wschody i zachody Księzyca

[sybic]

Witam,

 

Ponieważ Piotr nalega co bym mu jednak trochę teorii wyłożył

to postanowiłem przedstawić to tu, bo może ktoś zechciał

by coś uzupełnić/poprawić

 

zaczynamy od rocznika z niego odczytujemy 3 podstawowe

parametry Księżyca

RA - rektascencja

Dec - deklinacja

HP - paralaksa horyzontalna

dodatkowo odczytujemy też Czas Gwiazdowy

wszystko to dla momentu wschodu lub zachodu ...

oczywiście nie wiemy jeszcze kiedy to ma nastąpić, wszak to

mamy przecież obliczyć, więc odczytujemy to dla godz. 12 UT

 

 

kat godzinny wschodu lub zachodu wyrażony jest wzorem :

 

cos(L) = (-sin(F)*sin(Dec)-sin((Hz)) /(cos(F))*cos(Dec));

 

gdzie :

F - szerokosc geograficzna

Hz - wysokość Księzyca przy której obserwator

zobaczy jego górną krawędź (Def)

Hz = HP - 34' - D/2 (D - średnica tarczy)

można tez spotkać inny wzór :

-9.890199094634534e-3 + (1 - 0.2725076)*rad(HP/60)

- 3.874e-7 (wynik w radianach)

 

rzecz jasna wyrażenie po prawie stronie musi

być z przedziału -1 .. 1

jeśli tak nie jest to :

gdy >=1 to Księżyc nie wschodzi

gdy <=-1 to Księżyc nie zachodzi

 

kąt godzinny wschodu = 360 - L

kąt godzinny zachodu = L

 

 

teraz liczymy kąt godzinny kulminancji górnej

Lg = długość -15*strefa + Czas Gwiazdowy - 180 - RA

 

strefa : 1 - czas zimowy 2 - czas letni 0 - UT

 

 

i teraz godziny :

 

Gw = (Lw - Lg)/15

Gz = (Lz - Lg)/15

 

oczywiście są to momenty przybliżone,

teraz Gw trzeba przeliczyć na UT i dla tego momentu odczytać

dane z rocznika i jeszcze raz przeliczyć Gw

- tak ze 3 razy

 

jako uzupełnienie można tez obliczyć azymut

na to jest wiele wzorów, ale najlepszy jest ten :

 

Az:=2*arcsin( sqrt( (cos(F + Hz)-sin(F) ) / (2*cos(F)*cos(Hz)) ));

 

Rafał

[gazda]

O... O... O... O... O raju ale to wszystko fajnie wyglada ;D - jak dla mnie to gdybym nad tym pomyslal z godzinke to bym skapowal