Rambler's Top100Astronet    
  po tekstam   po forumu  vnutri temy
 

args[0]=message
args[1]=DB::DB::Message=HASH(0x30786d0)
Re[20]: Smeshenie perigeliya Merkuriya i drugih planet
17.01.2008 17:04 | V. V. Chazov

Dorogie druz'ya, prodolzhim.
Na moment T_{eph}=2454528.25075444 nado vychislit' polozhenie Merkuriya otnositel'no Solnca.
Sredi tekstovyh failov efemeridy DE200/LE200 nahodim fail ascp2000.200. V etom faile est' massiv iz 826 chisel, neobhodimyh dlya nashego momenta vremeni.
Deistvitel'no, etot massiv raspolozhen pod nomerom 94, nachinaetsya s yulianskoi daty 2454512.5, rabotaet do yulianskoi daty 2454544.5.
Dlya Merkuriya razbivaem interval na chetyre chasti, po 8 dnei:
1. 2454512.5 - 2454520.5,
2. 2454520.5 - 2454528.5,
3. 2454528.5 - 2454536.5,
4. 2454536.5 - 2454544.5,
Nash moment t_c=2454528.25075444 popadaet vo vtoroi interval, v takom sluchae
(t_a=2454520.5) < t_c < t_b=2454528.5.
Vychislyaem nomer yacheiki v massive iz 826 chisel, s kotoroi nachinayutsya koefficienty approksimacii dlya vtorogo intervala.
Nachal'nyi nomer i=3.
Chislo koefficientov dlya kazhdoi peremennoi k=12.
Chislo peremennyh m=3.
Obshee chislo yacheek dlya odnogo podintervala ravno (k*m)=36. Vtoroi interval nachinaetsya s elementa massiva nomer 3+36=39.
Pervye dvenadcat' chisel - koefficienty
a_0, a_1, a_2, ... , a_11 dlya koordinaty x.
Elementy s nomera 51 po nomer 62 - dlya koordinaty y.
Elementy s nomera 63 po nomer 74 - dlya koordinaty z.
Preobrazuem moment vremeni t_c k intervalu izmeneniya [-1,+1]:
x(t_c)=[t_c-0.5*(t_b+t_a)]/[0.5*(t_b-t_a)],
po rekurrentnomu algoritmu vychislyaem chislovye znacheniya polinomov Chebysheva ot T_0(x(t_c)) do T_{11}(xt(_c)),
s pomosh'yu koefficientov approksimacii vychislyaem po otdel'nosti tri koordinaty Merkuriya na zadannyi moment.
No my poluchim polozhenie Merkuriya otnositel'no baricentra Solnechnoi sistemy otnositel'no nepodvizhnogo ekvatora.
Dlya vychisleniya polozheniya Merkuriya otnositel'no Solnca nado vychislit' koordinaty Solnca otnositel'no baricentra. Vychisleniya vypolnyayutsya po tomu zhe algoritmu, tol'ko nomera yacheek v massive drugie.
Iz kazhdoi koordinaty Merkuriya nado vychest' sootvetstvuyushuyu koordinatu Solnca.
Dlya perehoda iz ploskosti ekvatora v ploskost' ekliptiki sledut vypolnit' povorot protiv chasovoi strelki na ugol \varepsilon - ugol naklona ekliptiki k ekvatoru.
S poklonom, Vash Vadim.


Forumy >> Astronomiya i Internet
Spisok  /  Derevo
Zagolovki  /  Annotacii  /  Tekst

Astronet | Nauchnaya set' | GAISh MGU | Poisk po MGU | O proekte | Avtoram

Kommentarii, voprosy? Pishite: info@astronet.ru ili syuda

Rambler's Top100 Yandeks citirovaniya