Posmotrel kod i ne ponyal, kuda v nem delos' nazvanie samogo risunka, t.k. vmesto nego tol'ko pi-pi-pi. Seichas poprobuyu eshe raz.

Ne ponyal. Ved' v rezhime prosmotra kartinka byla vidna.

Ne, Vadim, davai s relyativizmom pogodim, t.k. ya nikak ne mogu ponyat', chto za dannye aproksimirovany polinomami Chebysheva v efemeridah DE405, po tomu, chto ya ne ponimayu kak eto dannye radiolokacionnyh izmerenii mogut byt' raschetnymi dannymi, kotorye korrektiruyut dannye nablyudenii. Ved' dannye radiolokacionnyh izmerenii eto i est' dannye nablyudenii, po etomu ya ne poimu kak dannye nablyudenii mogut otkorrektirovat' dannye nablyudenii. Eti dannye radiolokacionnyh izmerenii mogut otkorrektirovat' tol'ko dannye vychislitel'nogo eksperimenta, poluchennye na modeli Solnechnoi sistemy, t.e. po uravneniyu, kotoroe my nikak ne mozhem zdes' pokazat'. I takim obrazom poluchaetsya, chto v efemeridah DE405 predstavleny dannye, poluchennye pri vychislitel'nom eksperimente na modeli Solnechnoi sistemy, prichem ne tol'ko dannye posle 2000 goda, chto samo soboi razumeetsya, no i s 1600 goda po 2000 god. Ne ponyatno tol'ko po kakim dannym sozdateli efemerid DE405 korrektirovali dannye vychislitel'nogo eksperimenta s 1600 po 1970 gody, t.k. posle 1970 goda yasno, chto oni ispol'zovali dannye radiolokacionnyh izmerenii teh observatorii, o kotoryh Vy napisali. Ili Vy hotite skazat', chto raschetnymi dannymi otkorrektirovany tol'ko dannye radiolokacionnyh izmerenii, a dannye nablyudenii proshlyh vekov ne korrektirovalis' raschetnymi dannymi, t.e. dannymi poluchennymi na matematicheskoi modeli.

S nailuchshimi pozhelaniyami Sergei Yudin.

Nu, esli ya vse pravil'no rasshifroval, to teper' mne voobshe nichego ne ponyatno. Kak eto eksperimental'nye dannye (prichem takie tochnye) mogut podgonyat'sya pod teoriyu (dazhe esli ona odobrena samim prezidentom). Obychno v nauke teoriya podgonyaetsya pod eksperimental'nye dannye. I potom, mne sovershenno ne ponyatno kakova rol' OTOshnoi modeli Solnechnoi sistemy JPL v sozdanie efemerid DE405. Esli eti efemeridy do 2000 goda yavlyayutsya approksimaciei eksperimental'nyh dannyh, to eto mozhno bylo sdelat' i bez modeli Solnechnoi sistemy, approksimirovav eksperimental'nye dannye na intervalah ot 8 do 32 dnei obychnoi kvadratichnoi zavisimost'yu s primeneniem metoda naimen'shih kvadratov (togda i polinomy Chebysheva ne nuzhny). No, esli polinomami Chebysheva approksimirovany dannye vychislitel'nogo eksperimenta provedennogo na modeli Solnechnoi sistemy, to dannye nablyudenii pri etom mogut byt' ispol'zovany tol'ko dlya kontrolya za pravil'nost'yu raboty modeli i pri etom nikak nel'zya govorit', chto my approksimirovali imenno dannye nablyudenii. A vot rol' modeli dlya sozdaniya efemerid posle 2000 goda mne vpolne ponyatna. I zdes' v efemeridah approksimirovany konechno zhe dannye vychislitel'nogo eksperimenta provedennogo na modeli Solnechnoi sistemy, t.k. nikakih nablyudatel'nyh dannyh iz budushego u nas byt' ne mozhet. Vot i voznikaet vopros chto zhe tam v efemeridah DE405 my imeem.

S nailuchshimi pozhelaniyami Sergei Yudin.

Poisk eksperimental'nyh dannyh nablyudenii ya konechno zhe prodolzhu, a seichas ya nachal provodit' predvaritel'nye eksperimenty po opredeleniyu optimal'nyh skorostei Solnechnoi sistemy (po osyam koordinat) i skorosti rasprostraneniya gravitacii, po kriteriyu, kogda dannye po vekovym smesheniyam parametrov orbit, poluchennye na N'yutonovskoi matematicheskoi modeli Solnechnoi sistemy, gde uchtena skorost' rasprostraneniya gravitacii, budut minimal'no otlichat'sya ot dannyh nablyudenii. Parametry orbit reshil poka vychislyat' tol'ko dlya Merkuriya, t.k. eti eksperimenty yavlyayutsya poka bol'she poiskovymi. Pri etom, dlya opredeleniya optimal'nyh parametrov sistemy ya budu ispol'zovat' mnogofaktornoe planirovanie. I seichas uzhe vypolnil pervuyu seriyu eksperimentov dlya chetyrehfaktornogo plana Boksa, gde byli sleduyushie znacheniya nulevyh urovnei i intervalov var'irovaniya parametrov VXsys=0+/-100, VYsys=0+/-100, VZsys=0+/- 100 i Vgr=11*Vsv+/-10*Vsv, gde skorosti sistemy v km/s, a skorost' gravitacii v dolyah skorosti sveta Vsv. Voobshe to pri poiskovyh opytah ne ispol'zuyut mnogofaktornoe planirovanie, a ispol'zuyut gradientnye metody, no mne hotelos' pobystree proverit' chast' programmy, otvechayushuyu za mnogofaktornoe planirovanie, i po etomu ya srazu nachal s nego. A seichas, chtoby bylo bolee ponyatno o chem budet idti rech', ya korotko izlozhu osnovy mnogofaktornogo planirovaniya.

Mnogofaktornoe planirovanie otlichaetsya ot odnofaktornogo, kogda v kazhdom novom eksperimente izmenyaetsya znachenie tol'ko odnogo parametra (faktora), tem, chto v kazhdom novom eksperimente izmenyayutsya po opredelennomu planu vse parametry. A posle obrabotki eksperimental'nyh dannyh my poluchaem uravnenie regressii (chashe vsego kvadratichnuyu approksimaciyu), kotoroe otrazhaet vliyanie vseh parametrov na kriterii optimizacii (celevuyu funkciyu). Esli u nas odnofaktornyi eksperiment to my poluchim takoe uravnenie

Y=k0+k1*X1+k2*X1^2 (1)
gde Y kriterii optimizacii, kotoryi nado minimizirovat' ili maksimizirovat', X1 optimiziruemyi parametr, a k0,k1 i k2 koefficienty, kotorye my poluchaem metodom naimen'shih kvadratov pri statisticheskoi obrabotke dannyh. A vot, esli u nas budet 4-e parametra, to my poluchim uravnenie regressii sleduyushego vida

Y = k0 + k1*X1 + k2*X2 + k3*X3 + k4*X4 +
+ k5*X1*X2 + k6*X1*X3 + k7*X1*X4 + k8*X2*X3 + k9*X2*X4 + k10*X3*X4 +
+ k11*X1^2 + k12*X2^2 + k13*X3^2 + k14*X4^2 (2)

Principial'no eto uravnenie otlichaetsya ot uravneniya, poluchennogo pri odnofaktornom planirovanie, nalichiem smeshannyh koefficientov k5 k10, kotorye otrazhayut smeshannoe vliyanie parametrov na kriterii optimizacii. Zadavaya znacheniya parametrov v kazhdom eksperimente po opredelennomu planu, my poluchaem minimal'noe kolichestvo eksperimentov neobhodimoe dlya polucheniya zavisimosti (2). Plany byvayut raznye i, naprimer, v programme Plank ya ispol'zoval chetyre raznyh plana -ortogonal'nyi, rotatabel'nyi, uniformrotatabel'nyi i pochti D- optimal'nyi Boksa. Otlichayutsya oni raznym kolichestvom urovnei var'irovaniya faktorov, velichinoi zvezdnogo plecha i kolichestvom nulevyh tochek, a v konechnom itoge eto vliyaet na to kak budet razmazana informaciya po gipersfere. Mne, naprimer, bol'she nravyatsya pochti D- optimal'nye plany Boksa i imenno takoi plan my i budem ispol'zovat'.

Dlya polucheniya uravneniya (2) nam neobhodimo po planu Boksa vypolnit' 24 eksperimenta. U etogo plana faktory var'iruyutsya tol'ko na 3-h urovnyah (vklyuchaya nulevoi) i zvezdnoe plecho budet ravno 1. Obychno v planah 0 oznachaet, chto parametr beretsya na nulevom urovne, +1 na nulevom urovne plyus interval var'irovaniya i -1 - na nulevom urovne minus interval var'irovaniya, a t.k. u nas zvezdnoe plecho ravno 1 i urovnei vsego tri, mozhno pisat' prosto 0 + - i po etomu plan Boksa dlya chetyreh faktorov budet vyglyadet' sleduyushim obrazom

1 + + + +
2 + + + -
3 + + - +
4 + + - -
5 + - + +
6 + - + -
7 + - - +
8 + - - -
Sleduyushie 8 opytov povtoryayut pervye 8. Tol'ko pervyi parametr prinimaet krugom znachenie -1.
17 + 0 0 0
18 0 0 0
19 0 + 0 0
20 0 0 0
21 0 0 + 0
22 0 0 0
23 0 0 0 +
24 0 0 0

Posle obrabotki dannyh vychislitel'nyh eksperimentov my poluchim uravnenie (2), kotoroe (esli u nas poluchilas' poverhnost', ne imeyushaya sedlovyh tochek) differenciruem po X1, X2, X3, X4 i poluchaem 4-e algebraicheskih uravneniya. Reshaya teper' sistemu etih uravnenii, nahodim optimal'nye znacheniya parametrov X1, X2, X3, X4 pri kotoryh nasha celevaya funkciya (kriterii optimizacii) prinimaet minimal'noe ili maksimal'noe znachenie. V nashem sluchae nam nuzhno minimal'noe znachenie, t.k. optimizaciyu my vedem po minimal'noi raznosti poluchennogo pri vychislitel'nom eksperimente znacheniya vekovogo smesheniya parametra i znacheniya parametra poluchennogo pri obrabotke dannyh naturnyh passivnyh eksperimentov (dannyh nablyudenii). Konkretno u menya v programme kriterii optimizacii vychislyaetsya po formule

Yu0(U) = SUMi,j ( kVesa(I, J) * Abs((YRas(I, J, U) - YNab(I, J)) / YNab(I, J)) / 100) (3)

Gde: Yu0(U) raznica mezhdu raschetnym YRas(I, J, U) i nablyudaemym YNab(I, J) znacheniem vekovogo smesheniya v U-om eksperimente (U=1...24) dlya I oi planety i J-go parametra
kVesa(I, J) vesovoi koefficient J-ogo parametra dlya I oi planety (ot 1 do 100).

Programma Solsys5 pozvolyaet proizvol'no vybirat' planety i parametry, po kotorym nado proizvodit' optimizaciyu, a vesovye koefficienty zadayutsya vruchnuyu samim issledovatelem v faile SolsysPlan.dat, gde zadayutsya takzhe i nulevye urovni parametrov i ih intervaly var'irovaniya (zadayutsya programmno). Voobshe to poluchenie edinyh kriteriev optimizacii pri reshenie mnogokriterial'nyh zadach vopros dazhe bolee slozhnyi, chem problema vekovyh smeshenii parametrov orbit, i ya obychno terpet' ne mogu kogda ispol'zuyut edinye kriterii optimizacii, poluchennye po principu avtomotovelofotoprimusradioruzh'e ili s primeneniem vesovyh koefficientov. Pravda podobnymi fokusami, gde procvetaet mahrovyi sub'ektivizm, chashe vsego pol'zuyutsya issledovateli social'no-ekonomicheskih ili tehniko- ekonomicheskih sistem. I mne dazhe prishlos' korennym obrazom modernizirovat' funkcional'no- stoimostnoi analiz, chtoby mozhno bylo ob'ektivnym putem poluchit' edinyi kriterii optimizacii dlya social'no-ekonomicheskih ili tehniko-ekonomicheskih sistem. No v nashem sluchae, t.k. vliyanie vseh parametrov na edinyi kriterii optimizacii odnoi prirody i nosit kolichestvennyi harakter, ya schitayu vpolne priemlemym ispol'zovanie vesovyh koefficientov dlya polucheniya edinogo kriteriya optimizacii, hotya nizhe i privozhu dannye optimizacii po kazhdomu parametru otdel'no.

Seichas ya vypolnil vse 24 eksperimenta plana, no kak okazalos' dannye po nekotorym eksperimentam, kogda skorost' gravitacii byla na nizhnem urovne, t.e. ravna odnoi skorosti sveta, poluchilis' ne prigodnymi dlya polucheniya za zadannyi promezhutok vremeni (400 let) dannyh po vekovym smesheniyam parametrov orbit. Na privedennom nizhe risunke dany grafiki izmeneniya parametrov orbit dlya 2-go eksperimenta plana, t.e. s parametrami +,+,+,-.

http://ser.t-k.ru/Ris/Plan1.gif (zerkalo http://modsys.narod.ru/Ris/Plan1.gif)

Kak vidim, dannye po perigeliyu i ekscentrisitetu hot' i ne ochen' horoshego kachestva, no prigodny dlya obrabotki, a vot dannye po uglu voshozhdeniya i uglu naklonu mozhno ispol'zovat' tol'ko dlya nebol'shogo promezhutka vremeni po tomu, chto ugol naklona orbity Merkuriya za 150 let umen'shilsya ot 7 do 3 gradusov, a potom nachal opyat' uvelichivat'sya. A ugol voshodyashego uzla za eto vremya umen'shilsya s 48 gradusov do nulya, a potom programma stala vydavat' ne minus 1 gradus, a 359 gradusov. Estestvenno pri takih masshtabnyh izmeneniyah parametrov orbit ya ne mogu approksimirovat' ih ni lineinoi, ni dazhe kvadratichnoi zavisimost'yu, chtoby naiti vekovye smesheniya. I v sleduyushei versii programmy ya postarayus' chto nibud' pridumat', chtoby izbavit'sya hotya by ot skachka v ugle voshozhdeniya, naprimer, budu ego opredelyat' v radianah. A dlya togo, chtoby vypolnit' kachestvennye eksperimenty na etoi versii programmy nado budet ili umen'shit' intervaly var'irovaniya skorostei sistemy po osyam ili uvelichit' nulevoi uroven' skorosti gravitacii. No, ne smotrya na ne ochen' horoshee kachestvo eksperimental'nyh dannyh, dlya illyustracii teoreticheskogo materiala, ya ih vse taki obrabotal i nashel koefficienty dlya 2-h uravnenii (2) po smesheniyu perigeliya Merkuriya i po smesheniyu ego ekscentrisiteta, hotya, kak sleduet iz uravneniya (3) mozhno bylo by i srazu poluchit' odno uravnenie po edinomu kriteriyu optimizacii. I nizhe Vy vidite kak po uravneniyam (2) kazhdyi iz parametrov vliyaet na kriterii optimizacii, kogda ostal'nye parametry zafiksirovany na nulevyh urovnyah.

http://ser.t- k.ru/Ris/Box1.gif (zerkalo http://modsys.narod.ru/Ris/Box1.gif)

V principe nichego interesnogo na etih grafikah net, i my vidim to, o chem ya uzhe skazal vyshe. T.e. pri teh skorostyah rasprostraneniya gravitacii, kotorye ya ispol'zoval v vychislitel'nyh eksperimentah, vse optimal'nye skorosti sistemy po osyam ... Z zhmutsya k nulevomu urovnyu, a optimal'noe znachenie skorosti rasprostraneniya gravitacii nahoditsya yavno za predelami var'irovaniya (v storonu uvelicheniya). Hotya, pri bol'shom zhelanii, i po uzhe poluchennym uravneniyam (2) mozhno chisto matematicheski naiti optimal'nye skorosti sistemy i optimal'nuyu skorost' rasprostraneniya gravitacii. Vot tol'ko v nashem sluchae etogo delat' nikak nel'zya (dazhe esli ochen' hochetsya), t.k. optimal'noe znachenie skorosti rasprostraneniya gravitacii nahoditsya yavno za predelami var'irovaniya parametrov. A ya lichno, k tomu zhe, starayus' eshe poluchit' vsegda optimum ne prosto v predelah intervala var'irovaniya parametrov, no i pri minimal'nyh intervalah var'irovaniya parametrov i chtoby optimum poluchilsya vblizi nulevogo urovnya i po kazhdomu parametru otdel'no (ne ochen' to ya vse taki doveryayu logike matematicheskoi celesoobraznosti i kombinirovannym kriteriyam optimizacii). No pri pravil'nom podbore vesovyh koefficientov kombinirovannyi kriterii optimizacii mozhet okazat'sya ochen' dazhe polezen.

V dopolnenie k sdelannym vyvodam mogu takzhe skazat', chto, kak izvestno iz literaturnyh istochnikov, absolyutnaya skorost' Solnechnoi sistemy, po dannym reliktovogo izlucheniya, sostavlyaet VXsys= -358, VYsys= 57, VZsys= -73 km/s, t.e. intervaly var'irovaniya skorosti sistemy ya zadal v plane vpolne priemlemye i, sledovatel'no, esli moya model' pravil'no otrazhaet ob'ektivnuyu real'nost', prichina stol' sil'nyh izmenenii parametrov orbit zaklyuchaetsya v ochen' malen'koi skorosti gravitacii. Po etomu ya seichas vypolnyu novyi plan, gde skorost' gravitacii zadam Vgr=50*Vsv+/-30*Vsv.

S nailuchshimi pozhelaniyami Sergei Yudin.

Poproboval na saite Instituta nebesnoi mehaniki i vychisleniya efemerid Francii naiti efemeridy otlichayushiesya ot efemerid JPL i dazhe na ih FTP servere ftp://ftp.imcce.fr/ v papke ftp://ftp.imcce.fr/pub/ephem/sun vrode by nashel kakie to efemeridy, no chto eto za efemeridy i kak imi pol'zovat'sya ya s moim znaniem angliiskogo yazyka i astronomii ne ponyal. Budu blagodaren, esli kto ni bud' pomozhet mne v etom razobrat'sya.

S nailuchshimi pozhelaniyami Sergei Yudin.