Odszkodowanie za dziurÄ™ w drodze
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³
:o 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 ;)
zanotowane.pldoc.pisz.plpdf.pisz.plmizuyashi.htw.pl