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 perigel