Astronet > Forumy > Astronomiya i Internet > Re[11]: Smeshenie perigeliya Merkuriya i drugih planet

S. Yu. Yudin · 18.11.2007 9:08

Ya uzhe proboval otkryt' etu temu na etom forume, no togda u menya voznikli problemy s kodirovkoi i tekst okazalsya ne chitaemym. Seichas poprobuyu drugoi variant kodirovki, t.k. ...

V. V. Chazov · 17.01.2008 16:04

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.

S. Yu. Yudin · 18.01.2008 9:50

< tbody>

Citata:

Nablyudeniya prohozhdenii Merkuriya i Venery po disku Soln ca, planomerno nachatye vo F rancii uchenymi-enciklopedi stami 18 veka, ispol'zovany tol'ko dlya uluchsheniya ravnomernosti shkal y vremeni.
Nablyudeniya, s kotorymi rabotali metr Urban Leverr'e i ser Saimon N'yukom, voobshe ne prinyaty vo vnimanie.
Druz'ya, etot fakt, veroyatno, bolee vsego ogorchi t Vas. Pri slovah o smeshenii perigeliya ssylayutsya na dostizheniya 19 veka, no o nablyudeniyah te h let - ni gu- gu.

Ya ne tol'ko ogorchilsya, no i nich ego ne ponyal. A razve stat'yu ChISLENNYE I ANALITIChESKIE TEORII DVIZhENIYa BOL'ShIH PLANET SOLNEChN O' SISTEMY, kotoruyu ya skachal s Vashego saita http://vadimchazov.narod.ru/text_pdf/particl1.pdf ne Vy pi sali. A to tam napisano, chto efemeridy ochen' dazhe sovpadayut s analiticheskimi efemeridami N'yukoma na intervale vremeni ot 1850 do 1895 godov. Na vsyakii sluchai privedu citatu iz stat' i.

//Efemerida DE200/LE200 vklyuchaet v sebya vse noveishie nauchnye razrabotki. Chislen noe integri rovanie uravnenii dvizheniya vypolneno v inercial'noi sisteme otscheta s nachalom v baricentre So lnechnoi sistemy. V kachestve argumenta integrirovaniya ispol'zovano ravnomernoe baricen tricheskoe dinamicheskoe vremya pryamoi potomok efe meridnogo vremeni. Uchteny relyativistskie effekty. Sistema gravitiruyushih tel dopolnena pyat' yu massivnymi asteroidami. Za osnovnuyu ploskost' vybrana ploskost' ekvatora, fiksirovannogo na epohu J200 0.0. Ispol'zovany radiolokacionnye nablyudeniya planet zemnoi gruppy, lazernye nablyudeniya Luny i dannye o parametrah dvizheniya kosmicheskih apparatov pri sblizhenii ih s bol'shimi plane tami. Rezul'taty analiticheskoi teorii dvizheniya chetyreh vnutrennih planet, sozdannoi N' yukomom, p rinyaty vo vnimanie samym neposredstvennym obrazom: novye chislennye efemeridy sostavl eny tak, chto nailuchshim obrazo m sovpadayut s analiticheskimi efemeridami N'yukoma na intervale vremeni ot 1850 do 1895 god ov.//

S nailuchshimi pozhelaniyami Sergei Yudin.

V. V. Chazov · 18.01.2008 15:34

Dorogie druz'ya, vse pravil'no.
Nablyudeniya 19 veka ne ispol'zovalis', a nachal'nye usloviya na intervale vremeni 1850-1895 gody podobrali tak, chtoby rezul'taty chislennogo integrirovaniya na etom intervale predstavlyali v srednem koordinaty vnutrennih planet, vychislennye po tablicam sera Saimona N'yukoma.
Bolee togo, no, pozhaluista, ne delaite nikakih orgvyvodov, sovremennye uchenye Parizha, gordyashiesya svoei sobstvennoi teoriei dvizheniya planet, voobshe ne prinimali vo vnimanie nikakie nablyudeniya. Oni ran'she vseh poluchali po dogovorennosti rezul'taty Laboratorii reaktivnogo dvizheniya SShA i approksimirovali ih ryadami Fur'e.
Eto ochen' nauchnye metody. Vazhno otchitat'sya tak, chtoby parlamenty i nacional'nye sobraniya poverili i vydelili assignovaniya na "novye razrabotki".
S poklonom, Vadim.

V. V. Chazov · 18.01.2008 15:48

Dear Friends, eto prostoe dopolnenie k otvetu pod nomerom 20.
Predstav'te sebe, srabotala privychka prostogo sovetskogo inzhenera: postavil zadachu, dolgo i nudno rasskazyval, kak ee sleduet reshat', a rezul'tatov ne privel.
A rezul'taty to est', vot oni:
geliocentricheskie eklipticheskie koordinaty, dolgota i shirota v gradusah, minutah, sekundah i rasstoyanie v astronomicheskih edinicah
2008 03 02 18 00 00.000 date in UTC
162 33 29.5 -00 00 35.5 00.99003 Luna
162 25 30.3 +00 00 00.8 00.99132 Zemlya
232 17 46.1 -00 29 19.7 00.45553 00.190 Merkurii
124 41 02.5 +01 47 17.2 01.64113 00.198 Mars
281 53 08.7 -01 26 57.3 00.72756 -3.749 Venera
276 40 32.7 +00 05 10.1 05.22526 -2.028 Yupiter
155 18 01.1 +01 39 17.2 09.28244 00.456 Saturn
348 28 03.0 -00 46 11.3 20.09633 05.945 Uran
321 43 53.3 -00 18 23.9 30.04105 07.973 Neptun
269 03 18.8 +06 22 30.1 31.42581 13.993 Pluton
Yours ever, Vadim.

S. Yu. Yudin · 18.01.2008 22:12

< tbody>

Citata:

Dorogie druz'ya, vse pravil'no.< br />Nablyudeniya 19 veka ne ispol'zovalis', a nachal'nye usloviya na intervale vremeni 1850-1895 go dy podobrali tak, chtoby rezul'taty chislennogo integrirovani ya na etom intervale predstavlyali v srednem koordinaty vnutrennih planet, vychislennye po tablicam sera Saimona N'yukoma.

Vadim, a nel'zya li poproshe. Ya ved' n e astronom, a Vy otlichno znaete zachem mne vse eto nado. Tak vot otvet'te poproshe. Mogu li ya ispol'zovat' dannye JPL, kak dannye passivnogo naturnogo eksperimenta (nablyudat el'nye dannye) dlya opredeleniya po rezul'tatam aktivnogo mnogofaktornogo vychislitel'nogo eks perimenta skorosti Sol nca i skorosti rasprostraneniya gravitacii ili dannye JPL eto te zhe samye dannye, chto daet i moya programma, no tol'ko nemnogo potochnee, po tomu, chto inte grirovanie provodilos' na otn ositel'no malen'

S. Yu. Yudin · 18.01.2008 22:14

kih intervalah vremeni, kogda nachal'nye usloviya dlya integrirovaniya postoyanno podpravlyalis' (ne izvestno kak), t.e. dannye JPL eto ne approksimaciya eksperimental'nyh nablyudatel'nyh dannyh, a gibridnye dannye s neizvestnym nazvaniem.

S nailuchshimi pozhelaniyami Sergei Yudin.

V. V. Chazov · 19.01.2008 7:58

Dorogie druz'ya, uvazhaemyi Sergei.
S kazhdym chasom ubezhdayus' v svoei dvoistvennosti. Voprosy stanovyatsya slozhnee, a otvechat' vse trudnee.
Dannye JPL - eto approksimaciya eksperimental'nyh nablyudatel'nyh dannyh.
Sergei, v Vashei terminologii - eto dannye passivnogo naturnogo eksperimenta (nablyudatel'nye dannye). Konechno zhe, takie zamechatel'nye i legkodostupnye materialy mozhno ispol'zovat' i dlya Vashei zadachi.
Ispol'zuyu chislennye efemeridy Laboratorii reaktivnogo dvizheniya SShA dlya obrabotki i planirovaniya nablyudenii, dlya prognoza dvizheniya malyh tel Solnechnoi sistemy. V SSSR pervye magnitnye lenty s efemeridami poyavilis' v god Olimpiady. Eto sobytie bylo ochen' zametnym.
S poklonom, Vadim.

S. Yu. Yudin · 20.01.2008 11:06

< tbody>

Vadim, bol'shoe spasibo z a ochen' podrobnye rekomendacii (kstati i v al goritme, kotoryi Vy mne prislali ochen' mnogo kommentari ev, a to mne na Membrane tozhe prislali al goritm i tozhe na Turbo Passal 7, no absolyut no bez kommentariev nichego ponyat' nel'zya), no u menya est' eshe o din vopros. Dopustim, ugol naklo na ekliptiki k ekvatoru v stan dartnuyu epohu J2000.0 mne izvest en 23 gr 26 min 21,448 sek, no k ak konkretno perevesti ekvatorial'nye koor dinaty v eklipticheskie ya somnevayus'. Nap rimer, nado li uchityvat' nutaciyu v standartn uyu epohu J2000.0 ravnuyu 9,2025 sek i zachem v fai lah JPL dany dannye po nutacii na konkretnuyu datu . Napishite pozhaluista konkre tnye formuly kak perevesti ekvatorial'nye koordinaty v eklipticheskie, a to, dazh e esli ya i razberus' po uche bnikam v etih tonkostyah, ya vse ravno budu somnevat'sya v svo ih formulah.

P.S. Svoi variant algoritma rascheta koordinat p o dannym JPL, t.e. adaptirovannyi dlya progra mmy Solsys ya napisal, no koordinaty poka poluchayut sya hot' i pohozhimi, no ne te, chto nado. Budu iskat' oshibku .

S nailuchshimi pozhelaniyami Sergei Yudin.

﷯ ﷯﷯﷯﷯ ﷯﷯ ﷯﷯﷯﷯﷯﷯﷯﷯﷯﷯Astronet﷯ | ﷯Nauchnaya set'﷯ | ﷯GAISh MGU﷯ | ﷯Poisk po MGU﷯ | ﷯O proekte﷯ | ﷯Avtoram﷯ ﷯ ﷯﷯Kommentarii, voprosy? Pishite: ﷯﷯info@astronet.ru﷯﷯ ili ﷯syuda﷯ ﷯﷯﷯﷯﷯﷯﷯﷯﷯﷯﷯ ﷯﷯﷯﷯﷯﷯﷯﷯﷯ ﷯﷯﷯﷯ ﷯﷯﷯﷯﷯﷯}(xt(_c)),
s pomosh'yu koefficientov approksimacii vychislyaem po otdel'nosti t i koordinaty Merkuriya na zadannyi moment.
N my poluchim polozhenie Merkuriya otnositel'no baricentra S lnechnoi sistemy otnositel'no nepodvizhnogo ekvatora.
D ya vychisleniya polozheniya Merkuriya otnositel'no Solnca n do vychislit' koordinaty Solnca otnositel'no baricentra. V chisleniya vypolnyayutsya po tomu zhe algoritmu, tol'ko nomera yacheek v massive drugie.
I kazhdoi koordinaty Merkuriya nado vychest' sootvetstvuyushuyu k ordinatu Solnca.
D ya perehoda iz ploskosti ekvatora v ploskost' ekliptiki s edut vypolnit' povorot protiv chasovoi strelki na ugol \varepsilon - u ol naklona ekliptiki k ekvatoru.
S poklonom, Vash Vadim. < table> & bsp;

Citata:

Dlya perehoda iz ploskosti e kvatora v ploskost' ekliptiki sledut vy polnit' povorot protiv chasovoi strelki na ugol \varepsilon - ugo l naklona ekliptiki k ekvatoru.
S poklono m, Vash Vadim.

S. Yu. Yudin · 20.01.2008 15:42

Voobshe to sami formuly u menya est', no ya ne uveren, chto eto okonchatel'noe preobrazovanie. Tak, esli ugol naklona ekliptiki k ekvatoru v standartnuyu epohu J2000.0 Betta= 23 gr 26 min 21,448 sek, to koordinaty v ploskosti ekliptiki dolzhny byt'
Xecl=X
Yecl= Y * cos(Betta) + Z* sin(Betta)
Zecl= Z * cos(Betta) - Y* sin(Betta)
Ya vse pravil'no zapisal? Nichego utochnyat' bol'she ne nado?

S nailuchshimi pozhelaniyami Sergei Yudin.

V. V. Chazov · 20.01.2008 16:03

Dobryi den', Sergei.
Znachenie ugla naklona ekliptiki k ekvatoru v standartnuyu epohu \varepsilon sleduet perevesti v radiany
\varepsilon = (3.141592653589793/180)*(23+(26+21.448/60)/60)
i vychislit' sinus \sin\varepsilon
kosinus \cos\varepsilon etoi velichiny.
Efemerida Laboratorii reaktivnogo dvizheniya SShA pozvolyaet vychislit' pryamougol'nye koordinaty
x, y, z,
v sisteme standarnogo ekvatora na lyubuyu datu.
Pryamougol'nye koordinaty
X, Y, Z
otnositel'no ploskosti fiksirovannoi ekliptiki vychislyayutsya po formulam:
X = x;
Y = +y*\cos\varepsilon + z*\sin\varepsilon;
Z = -y*\sin\varepsilon + z*\cos\varepsilon;
to est' s pomosh'yu povorota protiv chasovoi strelki vokrug osi OX na ugol \varepsilon.
Ploskost' ekliptiki sohranyaet neizmennoe polozhenie otnositel'no udalennyh (nepodvizhnyh) radioistochnikov.
Poklon vsem druz'yam, do vstrechi, Vash Vadim.