Avtor |
Soobshenie |
S. Yu. Yudin
|
Smeshenie perigeliya Merkuriya i drugih planet
|
18.11.2007 10: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. ...
|
|
Naverh |
|
|
S. Yu. Yudin
|
Re[13]: Smeshenie perigeliya Merkuriya i drugih planet
|
21.12.2007 21:47 |
|
<
tbody>Citata: |
Togda skorost' dvizheniya Solnca okolo 20 km/s, a eklipticheskaya sh
irota apeksa bol'she 50 gradusov. | <
p> Po etoi skorosti ya ne nahodil ne odnoi ssyl
ki, gde by dannye otlichalis' maksimum na paru procentov. A vot po skorostyam galaktiki
i mestnoi gruppy galaktik ochen' raznyatsya v raznyh literaturnyh istochnikah, no nigde ne ukaz
yvaetsya,
chto sko
|
|
Naverh |
|
|
S. Yu. Yudin
|
Re[14]: Smeshenie perigeliya Merkuriya i drugih planet
|
16.01.2008 16:05 |
|
Skachal faily s dannymi JPL (konkretno DE200 i DE405), no nikak ne mogu razobrat'sya kak ih ispol'zovat' dlya vychisleniya koordinat planet i Solnca v svoei programme, t.k.
ispol'zovat' uzhe gotovye programmy ya ne mogu, hotya i skachal odnu na S+. To chto vse dannye zapisany blokami s dannymi na 32 dnya i kazhdaya zapis' sostavlyaet 826 ili 1652 chisla
po 8 bait ya ponyal. Pervye dva chisla eto nachalo i konec perioda v 32 dnya v Yulianskom ischislenie, a dalee dolzhny idti koefficienty dlya polinomov Chebysheva, no ne yasno tol'ko
dlya koordinat ili i dlya skorostei planet tozhe, t.k. v raznyh mestah napisano po raznomu. Yasno i to, chto dlya raznyh planet eti koefficienty raspolozheny v strogo opredelennom
meste. Merkurii s 3-ei pozicii po 146, Venera s 147 po 182 i t.d. Dlya DE200 u menya est' rasshifrovka ne tol'ko etih pozicii, no i dano chislo koefficientov approksimacii, chastnyi
interval approksimacii v dnyah (dlya Merkuriya on raven 8 dnei) i kolichestvo approksimiruemyh dannyh (dlya vseh planet ukazano 3, t.e. poluchaetsya, chto koefficienty dany tol'ko
dlya koordinat), a kak eto budet dlya DE405 naiti ne mogu. A esli dazhe i naidu, to ne znayu chto za koefficienty ya schitayu i kuda ih potom podstavlyat'. Esli kto-to v kurse, pomogite.
S nailuchshimi pozhelaniyami Sergei Yudin.
|
|
Naverh |
|
|
V. V. Chazov
|
Re[15]: Smeshenie perigeliya Merkuriya i drugih planet
|
16.01.2008 22:31 |
|
Dorogie druz'ya. Uvazhaemyi Sergei, prostite velikodushno, ezheli nadoel Vam.
V efemeridnyh dannyh Laboratorii reaktivnogo dvizheniya soderzhatsya koefficienty polinomov
Chebysheva dlya vychisleniya rasstoyanii v kilometrah.
Differencirovanie kazhdogo polinoma po vremeni pozvolyaet vychislit' skorosti,
povtornoe differencirovanie daet znacheniya uskorenii nebesnyh tel.
Koordinaty vseh planet i Solnca (krome pary Zemlya-Luna) poluchayutsya otnositel'no
baricentra Solnechnoi sistemy. Dlya pary Zemlya-Luna dany koefficienty polinomov dlya
centra mass Zemlya-Luna. Otdel'no, pod nomerom 11, privodyatsya koefficienty dlya opredeleniya
polozheniya Luny otnositel'no Zemli.
Uvazhaemyi Sergei, mogu pereslat' po pochte podrobnye teksty i algoritmy. Moi adres prost:
vadimchazov@yandex.ru .
S poklonom, Vash Vadim.
|
|
Naverh |
|
|
S. Yu. Yudin
|
Re[16]: Smeshenie perigeliya Merkuriya i drugih planet
|
17.01.2008 0:43 |
|
<
tbody>Citata: |
Uvazhaemyi Sergei, prostite velikodushno, ezheli nadoel Vam. <
/td><
/tbody> A vot i ne ugodali (nichut' ne nadoeli), t.k. ya uzhe dva dnya dumayu o Vas
, a Vy tut kak tut - legki na pomine (dolgo zhit' budete). Kon
echno zhe vysylaite bystree i teksty i algoritmy, a to ya vtoroi den' borozdyu prostory I
nterneta i ne mogu naiti nichego putnego. Kstati kakoi fail dannyh posovetuete
DE200 ili DE405. Mne nuzhny dannye tol'ko po koordinatam planet i Solnca (skorosti ne
nuzhny) i zhelatel'no tol'ko nadezhnye dannye i imenno nablyudenii (pust' nemnogo podpravlen
nye JPL, no imenno nablyudenii, a ne prognoz na budushee). Ya tak dumayu nadezhnymi mozhno schitat'
dannye gde-to s 1850 goda (a mozhet ran'
she?). P.S. Soobshenie dlya vernosti dubliruyu pis'mom (moi adres ser@t-k.ru) S
nailuchshimi pozhelaniyami Sergei Yudi
n.
|
|
|
Naverh |
|
|
V. V. Chazov
|
Re[17]: Smeshenie perigeliya Merkuriya i drugih planet
|
17.01.2008 12:16 |
|
Dorogie druz'ya.
Shag za shagom reshim zadachu:
Dany
chislennye efemeridy Laboratorii reaktivnogo dvizheniya SShA.
Vychislit'
vektor polozheniya i vektor skorosti
planety Merkurii otnositel'no Solnca
2 marta 2008 goda v 21 chas moskovskogo dekretnogo zimnego vremeni.
Reshenie
Polinomy Chebysheva T_n(x).
Opredeleny dlya znachenii x v intervale -1<=x<=+1.
T_0(x)=1, nachal'noe znachenie,
T_1(x)=x, nachal'noe znachenie,
a eto uzhe rekurrentnaya formula dlya lyubogo
celogo polozhitel'nogo znacheniya n:
T_{n+2}(x)=2*x*T_{n+1}(x)-T_{n}(x).
Proizvodnye ot polinomov Chebysheva po peremennoi x
opredeleny v tom zhe intervale -1<=x<=+1.
T'_0(x)=0, nachal'noe znachenie,
T'_1(x)=1, nachal'noe znachenie,
a eto uzhe rekurrentnaya formula dlya lyubogo
celogo polozhitel'nogo znacheniya n:
T'_{n+2}(x)=2*x*T'_{n+1}(x)+2*T_{n+1}(x)-T'_{n}(x).
Esli izvesten moment vremeni t_c
lezhashii v intervale t_a<=t_c<=t_b,
to perehod k peremennoi x
i intervalu -1<=x<=+1 vypolnyaetsya po formule:
x(t_c)=[t_c-0.5*(t_b+t_a)]/[0.5*(t_b-t_a)].
V efemeridah dlya kazhdoi iz koordinat
dany chislennye znacheniya koefficientov
a_0, a_1, a_2, ... , a_m,
znachenie polozheniya po koordinate opredeleno formuloi:
r_x(t_c)=a_0*T_0(x(t_c))+a_1*T_1(x(t_c))
+a_2*T_2(x(t_c))+...a_m*T_m(x(t_c)),
znachenie skorosti po koordinate opredeleno formuloi:
v_x(t_c)=[a_0*T'_0(x(t_c))+a_1*T'_1(x(t_c))
+a_2*T'_2(x(t_c))+...a_m*T'_m(x(t_c))]*[(2*t_c)/(t_b-t_a)].
Edinica izmereniya polozhenii - kilometry.
Edinica izmereniya skorostei - kilometry za sutki.
Do vstrechi, Vash Vadim.
|
|
Naverh |
|
|
V. V. Chazov
|
Re[17]: Smeshenie perigeliya Merkuriya i drugih planet
|
17.01.2008 12:46 |
|
Dorogie druz'ya, prodolzhim.
Shag za shagom reshaem zadachu:
Dany
chislennye efemeridy Laboratorii reaktivnogo dvizheniya SShA.
Vychislit'
vektor polozheniya i vektor skorosti
planety Merkurii otnositel'no Solnca
2 marta 2008 goda v 21 chas moskovskogo dekretnogo zimnego vremeni.
Reshenie
Perevedem moskovskoe vremya vo Vsemirnoe koordinirovannoe vremya
UTC=02.03.2008 18 hour 00 minute 00 second (minus 3 chasa zimoi). .
Perevedem moment vremeni iz shkaly UTC v shkalu zemnogo vremeni TT
TT=UTC+65.184 second=02.03.2008 18 hour 00 minute 65.184 second.
Vychislim yulianskuyu datu, sootvetstvuyushuyu nashemu momentu vremeni:
TT(jd)=2454528.25075444 (sutki i doli sutok v shkale zemnogo vremeni).
S pomosh'yu nebol'shoi popravki (menee 0.002 sekundy vremeni)
pereidem iz shkaly zemnogo vremeni TT v shkalu
ravnomernogo vremeni T_{eph} chislennoi
efemeridy Laboratorii reaktivnogo dvizheniya SShA:
d=(2454528.25075444-2451545.0)/36525;
g=(3.141592654/180)*(357.528+35999.050*d); (v radianah)
T_{eph}=TT(jd)+(0.001658*Sin(g+0.0167*Sin(g)))/86400.e0;
V marte uglovaya velichina g ravna priblizitel'no 0 gradusov,
to est' i popravka na ravnomernuyu shkalu - nulevaya.<\br>
Imenno na moment T_{eph}=2454528.25075444,
sootvetstvuyushii nashei ishodnoi date,
sleduet vychislyat' polozheniya i skorosti po dannym efemeridy.
Do vstrechi, Vash Vadim.
|
|
Naverh |
|
|
S. Yu. Yudin
|
Re[18]: Smeshenie perigeliya Merkuriya i drugih planet
|
17.01.2008 15:19 |
|
Vadim, bol'shoe spasibo. Sil'no menya vyruchili s polinomami Chebysheva, no pereiti k vychisleniyam ya poka ne mogu, t.k. ostalis' ne vyyasnennymi nekotorye momenty. Naprimer,
Vy dlya vychisleniya i koordinat i skorostei zapisali v formulah odni i te zhe koefficienty a_0, a_1, a_2 a_m, a v faile DE405 soderzhitsya v dva raza bol'she zapisei (chisel)
v kazhdom bloke dannyh, chem v DE200, t.e. 1652 protiv 826. I, esli kolichestvo approksimiruemyh peremennyh i chislo koefficientov approksimacii ostalos' to zhe samoe, to ne ponyatno,
chto approksimiruyut lishnie 826 zapisei v faile DE405. Ya dumal, chto eto koefficienty dlya skorostei, a u Vas vyhodit, chto skorosti mozhno opredelit' po tem zhe samym koefficientam,
chto i koordinaty. Kakova voobshe struktura dannyh v DE405. Kak ya ponyal v DE200, naprimer, dlya Merkuriya s tret'ei pozicii idut koefficienty a_0_X, a_0_Y, a_0_Z a_m_X, a_m_Y,
a_m_Z, a kolichestvo koefficientov i drugie parametry opredelyayutsya po sheme, kotoruyu ya prochital v stat'e, raspolozhennoi na Vashem saite, a konkretno // Ves' massiv iz 826 chisel
rasshifrovyvaetsya s pomosh'yu pyati massivov celyh chisel: N object i j k l m 1 Mercury 3 146 12 8 3 2 Venus 147 182 12 32 3 3 Earth-Moon 183 272 15 16
3 4 Mars 273 302 10 32 3 5 Jupiter 303 329 9 32 3 6 Saturn 330 353 8 32 3 7 Uranus 354 377 8 32 3 8 Neptune 378 395 6 32 3 9 Pluto 396 413 6
32 3 10 Moon 414 701 12 4 3 11 Sun 702 746 15 32 3 12 Nutation 747 826 10 8 2 gde i - nachal'nyi nomer v massive, j - konechnyi nomer v massive, k - chislo
koefficientov approksimacii, l interval chastnoi approksimacii v dnyah vnutri obshego intervala, m - kolichestvo approksimiruemyh peremennyh.// Zhelatel'no by vzglyanut' na podobnuyu
shemu i dlya faila DE405. I potom, v svyazi s tem, chto dlya Merkuriya vnutri obshego intervala vremeni dlya dannyh sushestvuyut eshe i chastnye intervaly, to poluchaetsya, chto x(t_c)=[t_c-0.5*(t_b+t_a)]/[0.5*(t_b-t_a)]
nado rasschityvat' po znacheniyam t_a i t_b, kotorye budut ravny vremeni na koncah chastnogo intervala. Ili ya ne prav i ih nado brat' v nachale i konce obshego intervala, a koefficienty
ispol'zovat' iz chastnogo intervala. S nailuchshimi pozhelaniyami Sergei Yudin.
|
|
Naverh |
|
|
V. V. Chazov
|
Re[19]: Smeshenie perigeliya Merkuriya i drugih planet
|
17.01.2008 15:35 |
|
Prostaya rabota nad oshibkami:
Sleduet chitat'
znachenie skorosti po koordinate opredeleno formuloi:
v_x(t_c)=[a_0*T'_0(x(t_c))+a_1*T'_1(x(t_c)) +a_2*T'_2(x(t_c))+...a_m*T'_m(x(t_c))]*[2/(t_b-
t_a)].
2/(t_b-t_a) - eto masshtabnyi mnozhitel'.
Prostite, trudno s formulami na Forume. Vadim.
|
|
Naverh |
|
|
V. V. Chazov
|
Re[19]: Smeshenie perigeliya Merkuriya i drugih planet
|
17.01.2008 16:03 |
|
Dorogie druz'ya, na nashem Forume voznikayut prekrasnye voprosy, na kotorye ne mogu otvetit'
odnoznachno. Naprimer
Kstati, kakoi fail dannyh posovetuete DE200 ili DE405?
Mezhdunarodnyi astronomicheskii soyuz rekomenduet ispol'zovat' model' DE405/LE405.
Pri sozdanii etoi modeli byli uchteny rezul'taty kosmicheskih missii i nazemnye nablyudeniya
planet posle 1973 goda (do 1997 goda). V DE405/LE405 nemnogo ispravlena sistema otscheta,
to est' osnovnaya ploskost', v kotoroi vychislyayutsya polozheniya planet, sovpadaet s ploskost'yu
standartnogo ekvatora, zadavaemogo polozheniyami dalekih kvazarov. No kto iz nas eti kvazary
videl?
Ispol'zuyu v prakticheskih celyah DE200 i ne vstrechayu nepriyatnostei.
Eshe odin zamechatel'nyi vopros o nablyudeniyah, ispol'zovannyh pri postroenii
modelei.
V modelyah serii DE ispol'zovany nablyudeniya planet-gigantov na meridiannom kruge
observatorii Vashington, vypolnennye s 1904 goda. Uchteno neskol'ko serii nablyudenii
vnutrennih planet na razlichnyh meridiannyh krugah i astrolyabiyah, poluchennyh takzhe v 20
veke.
No osnovnoi ves poluchili radiolokacionnye izmereniya rasstoyanii kak do samih planet,
tak i do apparatov kosmicheskih missii. V sluchae Luny osnovnymi yavlyayutsya rezul'taty
lazernoi lokacii, provodimye s 1971 goda.
Nablyudeniya prohozhdenii Merkuriya i Venery po disku Solnca, planomerno nachatye vo Francii
uchenymi-enciklopedistami 18 veka, ispol'zovany tol'ko dlya uluchsheniya ravnomernosti shkaly
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 ogorchit Vas. Pri slovah o
smeshenii perigeliya ssylayutsya na dostizheniya 19 veka, no o nablyudeniyah teh let - ni gu-
gu.
Vprochem, i metr Urban Leverr'e i ser Saimon N'yukom prinimali k vychisleniyam tol'ko
sobstvennye nablyudeniya.
Uchenye vseh vremen i narodov - ves'ma original'ny, vsegda sluzhat po kontraktu, obyazany
otchitat'sya, kto pered intendantom finansov, kto pered kongressom, kto pered fondom. V
svyazi s etim kazhdomu prihoditsya vystavlyat' vpered sobstvennye dostizheniya - nablyudeniya,
teorii, modeli. O predshestvennikah luchshe ne napominat', vse ravno grantov s nimi ne
razdelish'.
Vadim.
|
|
Naverh |
|
|
V. V. Chazov
|
Re[19]: Smeshenie perigeliya Merkuriya i drugih planet
|
17.01.2008 16:26 |
|
Dear Friends.
Odin massiv, rabotayushii na intervale vremeni
32 dnya, v modeli DE405/LE405 soderzhit 1018 chisel.
Rasshifrovat' dannye mozhno s pomosh'yu tablicy:
N object i j k l m
1 Mercury 3 170 14 8 3
2 Venus 171 230 10 16 3
3 Earth-Moon 231 308 13 16 3
4 Mars 309 341 11 32 3
5 Jupiter 342 365 8 32 3
6 Saturn 366 386 7 32 3
7 Uranus 387 404 6 32 3
8 Neptune 405 422 6 32 3
9 Pluto 423 440 6 32 3
10 Moon 441 752 13 4 3
11 Sun 753 718 11 16 3
12 Nutation 819 898 10 8 2
13 Libration 899 1018 10 8 3
Chislo l - kolichestvo sutok korotkogo intervala approksimacii vnutri vsego intervala
dlinoi 32 dnya.
Dlya Merkuriya, naprimer, interval 32 dnya razdelen na chetyre podintervala po 8 dnei
kazhdyi.
Dlya Luny, naprimer, interval 32 dnya razdelen na vosem' podintervalov po 4 dnya kazhdyi.
Dlya Solnca, naprimer, interval 32 dnya razdelen na dva podintervala po 16 dnei kazhdyi.
Chislo m- kolichestvo peremennyh dlya approksimacii. Dlya vseh ob'ektov eto chislo ravno
3 - ibo v dekartovom mire tri koordinaty. Dlya nutacii approksimirovany izmeneniya dvuh
velichin, nutacii v dolgote i nutacii v naklone. Punkt 13 - eto approksimaciya chislovyh
znachenii treh uglov Eilera vrasheniya Luny, dannye neobhodimy pri obrabotke lazernyh
nablyudenii ugolkovyh otrazhatelei na nashem sputnike.
Do vstrechi. Vash Vadim.
|
|
Naverh |
|
|