Rambler's Top100Astronet    
  po tekstam   po forumu  vnutri temy
 

args[0]=message
args[1]=DB::DB::Message=HASH(0x315b380)
Re[28]: Smeshenie perigeliya Merkuriya i drugih planet
4.05.2008 13:01 | S. Yu. Yudin

Citata:

No dazhe esli i ne v nih, to cherez neskol'ko nedel' ya vse ravno vyskazhu VSE, chto ya dumayu ob etih efemiridah.

A nachnu ya svoi rasskaz, pozhalui, s konca. Na 99,9% eti efemeridy yavlyayutsya polnost'yu raschetnymi dannymi, poluchennymi na matematicheskoi modeli Solnechnoi sistemy, gde uskoreniya planet opredelyayutsya po teoreticheskoi formule http://ser.t- k.ru/Ris/....gif (zerkalo http://modsys.narod.ru/Ris/....gif) s uchetom trebovanii OTO, izlozhennoi v rabote http://iau-comm4.jpl.nasa.gov/XSChap8.pdf (E. Mails Stendish, Dzheims G. Vil'yams, Orbital'nye efemeridy Solnca, Luny i planet). T.e. dannye etih efemerid ne tol'ko ne yavlyayutsya obrabotannymi metodom naimen'shih kvadratov dannymi nablyudenii, kak ya i podozreval v samom nachale, no dazhe ne yavlyayutsya dannymi nablyudenii, sglazhennymi na matematicheskoi modeli, kak ya ostorozhno predpolozhil pozzhe. A teper' o tom, kak ya prishel k takomu vyvodu.

Posle togo, kak ya dopolnil programmu Solsys5 eshe i vozmozhnost'yu opredelyat' vekovye smesheniya bol'shoi poluosi ellipsa (srednego radiusa) i perioda obrasheniya, a takzhe nemnogo utochnil raschet nachal'nyh dannyh, ya reshil obnovit' tablicu 11 i dopolnit' ee i etimi smesheniyami (obnovlennuyu tablicu 11b mozhno skachat' otsyuda http://modsys.narod.ru/Arhiv/tabl11b.zip ). Ispol'zuya matematicheskii imitator v vide nabora approksimacii, kotoryi u astronomov nazyvaetsya teoriei N'yukoma (sm. u Duboshina str. 487 ili u Chebotareva str. 322), mozhno vychislit' parametry orbit na lyubuyu datu (v razumnyh predelah) s uchetom vekovyh smeshenii parametrov orbit, pri uslovie, chto planety dvizhutsya po ellipsam. U JPL tozhe est' takie approksimacii http://ssd.jpl.nasa.gov/txt/p_elem_t1.txt , no oni ne nazyvayutsya teoriei JPL, i parametry orbit na zadannuyu datu vychislyayutsya v etih approksimaciyah tol'ko po lineinym zavisimostyam
Alfa=k0+k1*dT
gde k1 eto vekovoe smeshenie parametra, a dT eto vremya v Yulianskih stoletiyah proshedshee s daty na kotoruyu fiksirovana epoha. U N'yukoma approksimacii bolee slozhnye. Est' chleny i s dT^2 i s dT^3, no oni dayut neznachitel'nye popravki, a osnovnoe smeshenie opredelyaetsya koefficientom k1 pri dT. I, t.k. JPL krugom pishet, chto ee efemeridy maksimal'no priblizheny k dannym N'yukoma, to ya oprometchivo ne stal sravnivat' koefficienty k1 v approksimaciyah JPL i N'yukoma, predpologaya, chto dannye JPL, esli i budut neznachitel'no otlichat'sya ot dannyh N'yukoma, to tol'ko v storonu bol'shei dostovernosti, t.e. budut prosto utochnyat' dannye N'yukoma. I dazhe, kogda ya v tablice 12 privel dannye N'yukoma po poluchennym im 4-m anomaliyam i dannye, poluchennye mnoyu pri obrabotke dannyh efemeridy DE405, i u menya poluchilos' bol'shoe rashozhdenie (v dva raza) po smesheniyu uzlov Venery, to ya ne zapodozril nichego ne ladnogo v dannyh JPL, a obvinil v etom N'yukoma. Vot chto ya togda napisal

//Ya uzhe pisal o tom, chto N'yukom ne vsegda pol'zovalsya obsheprinyatymi formulami, no privedennye v etoi tablice ego dannye stavyat menya prosto v tupik. Chto kasaetsya smeshenii perigeliya, to zdes' vse ponyatno, i on daet, kak eto obychno delayut pri opredelenie smeshenii perigeliya s ispol'zovaniem teorii vozmushenii, proizvedenie ekscentrisiteta orbity na smeshenie perigeliya Eks*dAlfaP, chto ya elementarno perevozhu v samo smeshenie dAlfaP. A vot perevod proizvedeniya sinusa ugla naklona orbity na smeshenie uzla voshozhdeniya sin(Betta2)*dAlfaU2 dlya Venery, u menya daet cifry ochen' otlichayushiesya ot cifr, kotorye ya poluchayu neposredstvenno zameryaya ugol voshozhdeniya v moei programme dAlfaU2 i ya budu priznatelen vsem, kto pomozhet mne s etim kazusom razobrat'sya. Nu i sovsem ya ne ponimayu pochemu u N'yukoma smeshenie ekscentrisiteta dano v uglovyh sekundah. Pravda, esli on ego opredelyal kak izmenenie parallaksa orbity, to togda vrode vse normal'no i mozhno ispol'zovat' i eti ego znacheniya dlya vyvodov ob anomal'nosti etogo smesheniya, no esli net, to opyat' voznikaet vopros.//

No vot, kogda ya stal opredelyat' koefficienty k1 po vekovym smesheniyam srednego radiusa i perioda obrasheniya planet, to u menya nikak ne poluchalis' znacheniya k1 pri doveritel'nyh intervalah men'she chem samo znachenie k1. Hotya ya proboval approksimirovat' dannye, poluchennye na imitatorah DE200 i DE405, ne tol'ko lineinoi, no i kvadratichnoi zavisimost'yu. I eto ne smotrya na to, chto na klassicheskoi modeli s uchetom skorosti gravitacii u menya poluchalis' otlichnye dannye i po smesheniyam srednego radiusa i po smesheniyam perioda obrasheniya pri ih approksimacii lineinoi zavisimost'yu. I tol'ko, kogda ya stal approksimirovat' parametry, poluchennye na imitatorah DE200 i DE405, konstantoi, t.e. prinyal vekovye smesheniya srednego radiusa ravnye nulyu, to ya poluchil samye malen'kie doveritel'nye intervaly. No eto rezko protivorechilo znacheniyam smeshenii srednego radiusa, rekomenduemym JPL, kak nablyudatel'nye dannye dlya rascheta parametrov orbit na zadannuyu datu. I eto ne smotrya na to, chto ya obrabatyval dannye efemerid DE200 i DE405, t.e. te zhe samye dannye JPL, kotorye oni nazyvayut nablyudatel'nymi dannymi. Vot tut u menya vozniklo podozrenie, chto znacheniya vekovyh smeshenii srednih radiusov v JPL poluchili po dannym, kotorye ochen' otlichayutsya ot dannyh approksimirovannyh v efemeridah DE200 i DE405, t.e. u nih imeyutsya i kakie to drugie nablyudatel'nye dannye, kotorye oni skryvayut.

Togda ya i reshil vernut'sya k opredeleniyu smeshenii, privedennyh v kombinirovannyh formulah N'yukomom v ego knige, posvyashennoi analizu primenennyh metodov i prinyatyh gipotez pri postroenie ego teorii, ot kotoryh ya otkazalsya posle togo, kak poluchil u N'yukoma smesheniya uzlov Venery ochen' otlichayusheesya ot dannyh JPL, a privedennye u nego smesheniya ekscentrisitetov v uglovyh sekundah postavili menya togda prosto v tupik. Na etot raz ya reshil snachala vychislit' smesheniya ne iz kombinirovannyh znachenii smeshenii privedennyh im dlya 4-h vnutrennih planet v knige, posvyashennoi analizu, a opredelil ih neposredstvenno iz approksimacii v ego teorii, prinyav smesheniya ravnymi koefficientu k1. Pri etom, t.k. u nego v teorii parametry orbit vychislyayutsya s uchetom vekovoi precessii, ya pri opredelenie smeshenii perigeliya i uzla voshozhdeniya vychital iz privedennyh im znachenii precessiyu, t.e. 5024,7 uglovyh sekund. A, chto kasaetsya smesheniya ekscentrisiteta dannogo v uglovyh sekundah, i kotoroe postavilo menya ranee v tupik, to ya dogadalsya, chto smeshenie ekscentrisiteta v bezrazmernyh edinicah N'yukom prinyal kak smeshenie v radianah, a zatem uzhe radiany, kak meru uglov, perevel v uglovye sekundy. Poluchennye, kak iz ego teorii, tak i iz tablicy dlya 4-h planet v knige posvyashennoi analizu, znacheniya vseh smeshenii pochti sovpali, no, t.k. oni ochen' otlichalis' ot teh, chto privodit JPL, ya reshil ih pereproverit' eshe i dannymi, privedennymi de Sitter (kstati bol'shim populyarizatorom OTO). Ego dannye po smesheniyam perigeliev i uzlov planet, kotorye ya nashel u Rouzvera, tozhe prakticheski sovpali s dannymi N'yukoma i, sledovatel'no, dannye N'yukoma sootvetstvuyut tem znacheniyam, kotorye byli v rasporyazhenie astronomov v konce 19-go i nachale 20-go veka. Za dostovernost' etih dannyh govorit i to, chto vse astronomy pochti ves' 20-yi vek ispol'zovali v svoih raschetah imenno teoriyu N'yukoma.

No togda, poluchaetsya, chto dannye JPL, ochen' otlichayushiesya ot dannyh N'yukoma, eto ne dannye nablyudenii, a kakie to drugie dannye, t.k. rashozhdeniya v sotni uglovyh sekund eto uzhe ne pogreshnost'. Kak te, tak i drugie dannye ya pomestil v utochnennuyu tablicu 11b, gde horosho vidno, chto po smesheniyam perigeliev i ekscentrisitetov rashozhdeniya imeyutsya, no bukval'no v neskol'ko procentov. A vot po smesheniyam uzlov i uglov naklona imeyutsya rashozhdeniya ne prosto v razy, no v nekotoryh sluchayah znacheniya imeyut dazhe protivopolozhnye znaki. Posle etogo ya dlya naglyadnosti isklyuchil iz tablicy 11b dannye, poluchennye na imitatore DE200 i na modeli s dannymi N'yukoma, a takzhe ubral doveritel'nye intervaly i oformil vse eto v vide tablicy 11bb.

Tablica 11bb. Vekovye smesheniya parametrov orbit rekomenduemye dlya rascheta parametrov orbit JPL i N'yukomom, t.e. smesheniya, poluchennye pri obrabotke dannyh nablyudenii, i smesheniya, poluchennye na programme Solsys5 s ispol'zovaniem imitatora JPL (efemerida DE405) i s ispol'zovaniem klassicheskoi matematicheskoi modeli Solnechnoi sistemy s nachal'nymi dannymi, rasschitannymi po parametram JPL za period s 1601 po 2001 gody (dlya Zemli s 1601 po 1951, t.k. pri priblizhenie k ekliptike epohi J2000 ugol naklona orbity stanovitsya ochen' malen'kim i voznikayut bol'shie pogreshnosti pri opredelenie perigeliya i uzla voshozhdeniya).

Parametr__Nablyudatel'nye_znacheniya___Poluchennye_na_programme_Solsys
___________JPL __________N'yukom**______imitator_______model'
___________JPL______teoriya__analiz_____DE4 05_____/R^2_____/R^n*
dAlfaP1___+577,7_____+575,1___+575,1_____+572,2_____+529,2___+573, 0
dAlfaU1___-451,2______-758,0___-753,7______-450,0_____-450,0___- 449,9
dBetta1____-21,4______+6,70____+7,14______-21,4______-21,4_____- 21,4
dEks1_____+19,1______+20,5____+16,3______+20,5_____+20,5_____+20,5

dAlfaP2___ _+9,7______+44,3_____+42,5______+41,2_____+32,0____+47,7
dAlfaU2___-999,7_____- 1785,2___-1780,6_____-998,5_____-998,2___-998,5
dBetta2____- 2,84______+3,62____+3,87______-2,50______-2,50_____-2,50
dEks2_____-41,1______- 47,7_____-45,9_______-48,3______-49,1____- 49,2

dAlfaP3___+1163,8____+1164,3__+1162,9___+1155,9____+1135,4__+1146,0
dAlfaU3__ ___---________---________---______-836,4_____-546,6____-546,3
dBetta3____-46,6_______--- ______-47,1______-47,2______-47,2_____-47,2
dEks3_____-43,9______-41,8_____-41,5______- 42,0______-42,7_____- 42,7

dAlfaP4___+1599,9____+1602,0__+1602,7____+1600,2____+1599,2__+1604,8
dAlfaU4_ __-1053,3_____-2249,1__-2248,8_____-1050,3_____-1049,5__-1049,5
dBetta4____-29,3_______- 2,43____-2,26_______-28,9______-28,9_____- 28,9
dEks4_____+78,8______+92,1____+92,1______+91,4______+95,8_____+95,8

Smesheniya uglov perigeliya dAlfaP, uzla voshozhdeniya dAlfaU i naklona dBetta dany v uglovyh sekundah. Smesheniya ekscentrisiteta dEks dany v bezrazmernyh edinicah i uvelicheny v 10^6 raz. Vekovye smesheniya parametrov (dAlfa) na programme Solsys5 polucheny pri approksimacii parametrov lineinoi zavisimost'yu Alfa=k0+k1*dT (otsyuda dAlfa=k1).

* - pokazatel' stepeni n vo vtoroi modeli v formule N'yutona F=G*m*M/R^n bralsya takim, kakim on byl u N'yukoma, t.e. n=2,0000001612 (u Holla bylo n=2,0000001574).

** - v pervoi kolonke dany smesheniya, kotorye poluchayutsya iz approksimacii N'yukoma (teoriya N'yukoma), dlya 1900 goda i pri vychitanie iz lineinyh chlenov uglov perigeliya i voshodyashego uzla vekovoi precessii 5024,7 uglovyh sekund (Spravochnoe rukovodstvo po nebesnoi mehanike i astrodinamike, pod redakciei G.N. Duboshina, izdanie vtoroe, dopolnennoe i pererabotannoe - M.: Nauka, 1976. 864 s. //sm. str. 487-494//), a vo vtoroi kolonke dany smesheniya, kotorye vzyaty iz knigi N'yukoma posvyashennoi analizu primenennyh metodov i prinyatyh gipotez pri postroenie ego teorii (Rouzver N.T. Perigelii Merkuriya. Ot Lever'e do Einshteina: Per. s angl. M.: Mir, 1985. 246 s. //sm. tabl. 3.1//). Znachenie vekovoi precessii ya prinyal 5024,7 uglovyh sekund, chtoby poluchennoe iz teorii nablyudaemoe znachenie po smesheniyu perigeliya Merkuriya, sovpalo so znacheniem privedennym v knige analiza, t.e. chtoby poluchilos' 575,06 uglovyh sekund.

Vot tut to, glyadya na etu tablicu, ya i ponyal, chto efemeridy DE200 i DE405 voobshe ne imeyut nikakogo otnosheniya k nablyudatel'nym dannym, a yavlyayutsya polnost'yu raschetnymi dannymi, poluchennymi na matematicheskoi modeli, v kotoroi uchteny trebovaniya OTO. I, t.k. OTO pozvolyaet tol'ko nemnogo smestit' perigelii v klassicheskoi modeli, po etomu kak raz po perigeliyu bolee-menee i sovpadayut nablyudatel'nye dannye, privedennye N'yukomom, i dannye, poluchennye na modeli JPL. A chto kasaetsya znacheniya dAlfaP2, to zdes' bystree vsego v JPL prosto zavralis', napisav, chto rekomenduemoe znachenie vekovogo smesheniya budet +9,7 uglovyh sekund, t.e. primerno stol'ko skol'ko trebuet OTO, t.k. u menya po dannym DE405 poluchaetsya znachenie, sovpadayushee s dannymi N'yukoma (spravedlivosti radi, nado zametit', chto opredelit' tochno dAlfaP2 dovol'no taki slozhno). Chto zhe kasaetsya vekovyh smeshenii ekscentrisitetov, to zdes' hot' i poluchayutsya rashozhdeniya mezhdu nablyudatel'nymi dannymi i raschetnymi, no nikakih opredelennyh vyvodov o modelyah po etim smesheniyam sdelat' nel'zya. A vot uchest' vekovye smesheniya uzlov voshozhdeniya i naklonov orbit OTOshnaya model' JPL ne mozhet tochno i po etomu dannye, kotorye polucheny na etoi modeli i zalozhennye zatem v efemeridy DE200 i DE405, sovershenno i ne otrazhayut vekovye smesheniya etih parametrov.

No, esli delo obstoit tak, kak ya narisoval, to togda voznikaet vopros a zachem JPL potrebovalas' slozhneishaya matematicheskaya model', postroennaya s uchetom trebovanii OTO, esli pochti takie zhe raschetnye dannye mozhno bylo poluchit' i na klassicheskoi modeli. Dlya etogo nado bylo vsego-navsego ispol'zovat' zakon Holla, t.e. prosto zamenit' v formule N'yutona pokazatel' stepeni v znamenatele na znachenie nemnogo otlichnoe ot 2, kak eto sdelal N'yukom dlya sostavleniya svoih tablic dvizheniya planet i kotorye sluzhili dlya vychisleniya efemerid v astronomicheskih ezhegodnikah s 1901 po 1959 god. N'yukom prinyal znachenie pokazatelya stepeni 2,0000001612, chto dolzhno bylo po ego dannym dat' dopolnitel'noe smeshenie perigeliya dlya 4-h vnutrennih planet, sootvetstvenno, 43,37 16,98 10,45 5,55 uglovyh sekund. Vzglyanite na poslednyuyu kolonku v tablice 11bb, gde dany smesheniya parametrov, poluchennye mnoyu na klassicheskoi modeli s dannymi JPL s modernizirovannym po Hollu zakonom N'yutona i u menya popravki poluchilis' 43,8 15,7 10,6 - 5,6 uglovyh sekund. Pri etom vidno, chto po perigeliyam smesheniya poluchayutsya primerno takie zhe, kak i pri modernizacii zakona N'yutona preobrazovaniyami OTO, t.e. po imitatoru DE405. Pravda pri etom, takzhe, kak i s ispol'zovaniem OTO, nam ne udaetsya poluchit' nuzhnye smesheniya uzlov i naklona orbity, no formula Holla, v otlichie ot OTO, i ne pretenduet ni na kakuyu fundamental'nost', a prosto pozvolyaet poluchit' nuzhnyi chastnyi rezul'tat, t.e. dopolnitel'noe smeshenie po perigeliyam planet.

Takim obrazom, poluchaetsya, chto vsya eta slozhneishaya matematicheskaya abstrakciya, t.e. OTO, byla ispol'zovana JPL tol'ko dlya togo, chtoby pridat' naukoobraznost' nebol'shomu uvelicheniyu vekovyh smeshenii perigeliev planet, po sravneniyu s rezul'tatami, poluchennymi na klassicheskoi modeli Solnechnoi sistemy i, sledovatel'no, nikakoi nauchnoi i prakticheskoi cennosti ne predstavlyaet. T.e. poluchaetsya, chto JPL prosto naveshala lapshi na ushi kongressmenam SShA o tom kakuyu bol'shuyu nauchnuyu i prakticheskuyu cennost' predstavlyayut ih issledovaniya, chtoby poluchit' finansirovanie. Vprochem, kak ya ponyal iz ob'yasnenii Vadima Chazova, tozhe samoe seichas delayut i francuzy, razlagaya dannye, poluchennye v JPL v nikomu ne nuzhnye ryady s tysyachami chlenov, kolichestvo kotoryh dolzhno pridat' naukoobraznosti uzhe ih issledovaniyam, chtoby poluchit' finansirovanie uzhe ot svoego pravitel'stva. No eto, kak govorit'sya ih lichnye problemy, no ved' JPL teper' obmanyvaet vse nauchnye organizacii mira (a mozhet byt' oni sami hotyat obmanyvat'sya), zastavlyaya ih pol'zovat'sya ih efemeridami, kotorye ona ob'yavila obrabotannymi dannymi nablyudenii, kotorye k tomu zhe predskazyvayut s vysokoi tochnost'yu polozheniya planet v budushem, t.e. predstavlyayut i prakticheskuyu cennost'. No, kak vidno iz nizheprivedennoi tablichki, efemeridy JPL, poluchennye na OTOshnoi modeli, ne verno otrazhayut dazhe smesheniya perigeliev planet, t.e. togo parametra, kotoryi yavlyaetsya odnim iz glavnyh eksperimental'nyh podtverzhdenii vernosti OTO kak teorii, i kotoryi oni prosto obyazany otrazhat' verno.

Tablica 3b. Anomal'nye ostatki nablyudatel'nyh dannyh po smesheniyu perigeliev planet, kotorye ne ob'yasnyayutsya teoriei N'yutona, i ostatok, rasschitannyi analiticheski po formule OTO i chislennymi metodami na OTOshnoi modeli JPL, t.e. po efemeridam DE405. Ostatki ot nablyudatel'nyh dannyh, opredelennye po raschetnym dannym samih avtorov, dany v chislitele, a po moim raschetnym dannym v znamenatele, gde u menya dlya dannyh JPL raschetnoe znachenie polucheno s sistemoi mass JPL, a dlya ostal'nyh dannyh s sistemoi mass N'yukoma.

____________________________Merkurii___Venera___Zemlya___Mars
ostatok Lever'e (Chebotarev)____38,3/---_____---______---______8.0*
ostatok N'yukoma (Rouzver) _____41,2/45,6__-7,3/9,9__6,0/25,9__8,0/0,3
ostatok Dankoma (Brumberg) ____43,1/45,7___8,1/1,6__5,0/16,7___---
ostatok JPL (Yudin)___________---/48,5___---/- 22,3__---/28,4__---/0,7
OTO Einshteina (Subbotin) ______43,0_______8,6______3,8______1,4
OTO DE405 (Yudin) _____________43,0_______9,2______20,5_____1,0
* - Subbotin (str.60) pishet, chto u Lever'e poluchilos' by znachenie blizkoe k znacheniyu N'yukoma esli by Lever'e ispol'zoval bolee tochnye znacheniya mass planet.

(sm. prodolzhenie)




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