Rambler's Top100Astronet    
  по текстам   по форуму  внутри темы 		 

args[0]=message
args[1]=DB::DB::Message=HASH(0x4d35e90)
Re[2]: Смещение перигелия Меркурия и других планет
5.03.2008 8:55 | С. Ю. Юдин

Наконец-то теоретическое обоснование моей методики обработки экспериментальных данных (наблюдательных и полученных при вычислительном эксперименте на модели) более менее устаканилась (хотя и остались некоторые пробелы) и я решил выложить некоторые полученные результаты для обсуждения. Ранее я уже приводил данные по смещениям перигелиев планет в таблице 4, а сейчас я уберу из нее свои данные, которые получены на модели, и данные НАСА, которые приводил со ссылкой на Хайдарова, т.к. последние отличаются от найденных мною данных НАСА, а конкретно от данных JPL, которая является подразделением НАСА. Но в то же время я и дополню эту табличку данными, полученными мною после обработки первичных данных JPL (из фонда DE405), которые признаются практически всеми именно данными наблюдений. А по причине, о которой я скажу позже, я уберу данные по Юпитеру и Сатурну.

Таблица 4а. Обработанные разными авторами экспериментальные данные наблюдений для определения вековых смещений перигелиев 4-х внутренних планет. В скобках указан источник, откуда взяты уже обработанные данные (у меня данные JPL это исходные данные).

______________Меркурий__Венера__Земля_Марс
JPL (JPL)__________577,7____9,7__1163,8__1599,9
Юдин (JPL)________572,2___37,8__1157,2__1599,5
Ньюком (Роузвер)___575,1___42,5__1162,9__1602,7
Данком (Брумберг)__575,2___34,3__1153,8___---
Брауэр, Вурком _____589,1___---_____---___1800,0
(Субботин)

Если не принимать во внимание данные Брауэра и Вуркома, которые я нашел у Субботина, то можно сказать, что, не смотря на разные исходные данные и примененные для этого методики обработки данных наблюдений, результаты у всех авторов получаются практически одни и те же (кроме данных JPL по Венере, но и здесь не все потеряно). Если учесть, что, приведенные значения вековых смещений, JPL рекомендует использовать для расчета параметров орбит в интервале с 1800 по 2050 год, то, наверное, и данные наблюдений для их определения они использовали за этот период, а у всех остальных авторов (насколько я об этом знаю) использованы данные наблюдений примерно с 17-го века. А в этом случае вполне может получить для 2000 года и значение приведенное JPL (у меня, например, в таблице 4 вообще приведено расчетное значение для 2000 года минус 5). Таким образом, можно сказать что, не смотря на некоторые существующие пока теоретические пробелы, моя методика прошла предварительную проверку практикой и, следовательно, ее можно использовать уже сейчас в практической деятельности.

Но прежде чем я перейду к практическому использованию моей методики для сложных расчетов по определению скорости распространения гравитации и абсолютной скорости Солнечной системы я бы хотел не только привести, полученные мною сейчас, после обработки экспериментальных данных наблюдений, предоставленных JPL, значения вековых смещений всех остальных (кроме полуоси эллипса) параметров внутренних планет, но и остановится на некоторых тонкостях определения вековых смещений планет с использованием моей методики. Первым делом я обработал по своей методике данные наблюдений JPL за период с 1.01.1601 года по 1.01.2001 года (использовано 2-а разных фонда DE200 и DE405), а также данные полученные за этот же период при проведение вычислительных экспериментов на разных моделях систем. Вообще то обе модели были чисто Ньютоновскими, но отличались массами планет и начальными параметрами орбит планет, которые задавались или по данным JPL или по данным Ньюкома.

При этом я привожу не только полученные мною значения вековых смещений параметров орбит (AlfaP перигелий, AlfaU узел восхождения, Betta наклон орбиты и Eks эксцентриситет), но привожу и доверительный интервал для этих значений при доверительной вероятности (надежности) 95%. Полученные смещения углов у меня даны в угловых секндах, а безразмерные значения вековых смещений эксцентриситетов увеличены мною в 3600 раз. Шаг численного решения дифференциальных уравнений и обработки данных наблюдений JPL для всех планет был постоянным 3600 с, но при определение угла в перигелии и угла узла восхождения он уменьшался так, чтобы интервал погрешности замера угла в перигелии не превышал одной угловой секунды (например, для Меркурия основной шаг уменьшался в 1000 раз). А при статистической обработке данных я всю выборку полученных данных по параметрам орбит разбивал на 5 групп, вычислял в них средние значения параметра, а затем эти пять точек аппроксимировал линейной зависимостью

Alfa=k0+k1*dT (2a)

а взяв первую производную я получал значения вековых смещений параметров орбит - k1 (если у нас время dT в формулах 2a будет измеряться столетиями). Т.к. многим все эти, полученные мною значения параметров орбит будут не интересны, я прикрепляю их в виде текстового файла таблица11, который можно скачать отсюда http://modsys.narod.ru/Arhiv/tabl11.zip , а здесь приведу только полученные мною данные по 4-м аномалиям, которые обнаружил Ньюком, и о которых мы сейчас и поговорим. Но прежде я сделаю некоторые общие выводы по данным, приведенным в таблице11. Во-первых, радуют довольно маленькие доверительные интервалы получившиеся при обработке данных, а во-вторых, не смотря на столь хорошую точность, значения вековых смещений полученных по фондам DE200 и DE405 почти все совпадают до последней цифры. Таким образом, я, во-первых, в следующей версии программы фонд DE200 уберу, т.к. не вижу, чтобы данные фонда DE405 были, как пишут, гораздо точнее данных фонда DE200 (а для меня более важно, что они не подогнаны под нужный результат), а во-вторых, и данные, полученные на двух разных моделях систем, тоже не очень сильно отличаются друг от друга и, по этому, при предварительном анализе я буду пользоваться только данными, полученными на модели с параметрами JPL, что является оправданным еще и по тому, что и экспериментальные данные у нас имеются пока только предоставленные JPL. А теперь давайте вернемся к аномалиям Ньюкома, где рядом с данными Ньюкома я приведу и полученные мною данные при обработке данных фонда DE405 и полученные на модели с параметрами JPL.

Таблица 12. Аномальные смещения параметров орбит планет.
Параметр_____Наблюдения_______Теория__________Разность
Eks1*dAlfaP1__118,24+/-0,40*___109,76+/-0,16*___8,48+/-0,43*
dAlfaP1______575,06+/-1,95*___533,82+/-0,78*___41,24+/-2,09*
dAlfaP1______572,17+/-0,30____529,38+/-0,31____42,79+/-0,31

Eks4*dAlfaP4__149,55+/-0,35*___148,80+/-0,04*___0,75+/-0,35*
dAlfaP4_____1602,69+/-3,75*__1594,65+/-0,43*___8,04+/-3,75*
dAlfaP4_____1599,53+/-6,25___1601,67+/-6,29____2,14+/-6,29

sin(Betta2)*dAlfaU2_-105,40+/-0,12*__-106,00+/-0,12*___0,60+/-0,17*
dAlfaU2_________-1780,68+/-2,03*__-1790,81+/-2,03*__10,14+/-2,87*
dAlfaU2__________-998,75+/-1,25___-998,96+/-1,44_____0,21+/-1,44

dEks1_______3,36+/-0,50*______4,24+/-0,01*_____0,88+/-0,50*
dEks1_____0,07375+/-0,0010__0,07402+/-0,0011__0,00027+/-0,0011

* - здесь Ньюком приводит средние ошибки, но для определения аномальности смещения он использует вероятные ошибки, которые получаются умножением средних ошибок на 0,67454. Вот только я нигде не нашел для какой доверительной вероятности (надежности) он приводит эти вероятные ошибки, а без указания доверительной вероятности все эти вероятные ошибки не несут никакой информации.

В приведенных выше данных, первая строка (или две первых) это данные Ньюкома, а последняя строка это мои данные. Обозначения вековых смещений параметров орбит соответствуют принятым в моей программе Solsys5 и в вышеприведенной таблице11 (AlfaP перигелий, AlfaU узел восхождения, Betta наклон орбиты и Eks эксцентриситет), а последняя цифра соответствует порядковому номеру планеты. Я уже писал о том, что Ньюком не всегда пользовался общепринятыми формулами, но приведенные в этой таблице его данные ставят меня просто в тупик. Что касается смещений перигелия, то здесь все понятно, и он дает, как это обычно делают при определение смещений перигелия с использованием теории возмущений, произведение эксцентриситета орбиты на смещение перигелия Eks*dAlfaP, что я элементарно перевожу в само смещение dAlfaP. А вот перевод произведения синуса угла наклона орбиты на смещение узла восхождения sin(Betta2)*dAlfaU2 для Венеры, у меня дает цифры очень отличающиеся от цифр, которые я получаю непосредственно замеряя угол восхождения в моей программе dAlfaU2 и я буду признателен всем, кто поможет мне с этим казусом разобраться. Ну и совсем я не понимаю почему у Ньюкома смещение эксцентриситета дано в угловых секундах. Правда, если он его определял как изменение параллакса орбиты, то тогда вроде все нормально и можно использовать и эти его значения для выводов об аномальности этого смещения, но если нет, то опять возникает вопрос.

Хорошо, пусть будем считать, что все эти казусы нами успешно разрешены и доверительная вероятность у Ньюкома никак не меньше 95%. Все равно ничего хорошего в этом случае у нас не получается и из 4-х аномалий Ньюкома я вижу по своим данным только одну аномалию для перигелия Меркурия. А аномального смещения узла восхождения у Венеры, перигелия Марса и эксцентриситета у Меркурия я не вижу, т.к., при заданной мною доверительной вероятности, полученные доверительные интервалы для этих параметров перекрывают и без того микроскопическую разницу в полученных мною наблюдаемых и расчетных значениях. Хотя, вот у Земли впору уже говорить об аномальном смещение перигелия, т.к. доверительные интервалы наблюдаемых и расчетных значений с большим трудом перекрывают друг друга. Априори будем считать, что у Ньюкома с полученными данными все более-менее в порядке и тогда остается только подумать, что мы можем сделать при обработке данных полученных мною на математической модели и из фонда DE405 (если там у JPL тоже все в порядке). И через пару дней я постараюсь что ни-будь придумать в этом направление и доложу свои соображения, а то получается, что, после приведенных мною данных, ОТО вообще может спать спокойно, т.к. ни смещение узлов Венеры, ни смещение эксцентриситета Меркурия ей теперь объяснять не надо. Да и небольшая проблема с перигелием Марса у нее тоже отпала.

С наилучшими пожеланиями Сергей Юдин.




Форумы >> Астрономия и Интернет
Список  /  Дерево
Заголовки  /  Аннотации  /  Текст
Астронет | Научная сеть | ГАИШ МГУ | Поиск по МГУ | О проекте | Авторам

Комментарии, вопросы? Пишите: info@astronet.ru или сюда

Rambler's Top100 Яндекс цитирования