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

args[0]=message
args[1]=DB::DB::Message=HASH(0x2de4950)
Re[18]: Смещение перигелия Меркурия и других планет
17.01.2008 14:19 | С. Ю. Юдин

Вадим, большое спасибо. Сильно меня выручили с полиномами Чебышева, но перейти к вычислениям я пока не могу, т.к. остались не выясненными некоторые моменты. Например, Вы для вычисления и координат и скоростей записали в формулах одни и те же коэффициенты a_0, a_1, a_2 a_m, а в файле DE405 содержится в два раза больше записей (чисел) в каждом блоке данных, чем в DE200, т.е. 1652 против 826. И, если количество аппроксимируемых переменных и число коэффициентов аппроксимации осталось то же самое, то не понятно, что аппроксимируют лишние 826 записей в файле DE405. Я думал, что это коэффициенты для скоростей, а у Вас выходит, что скорости можно определить по тем же самым коэффициентам, что и координаты. Какова вообще структура данных в DE405. Как я понял в DE200, например, для Меркурия с третьей позиции идут коэффициенты a_0_X, a_0_Y, a_0_Z a_m_X, a_m_Y, a_m_Z, а количество коэффициентов и другие параметры определяются по схеме, которую я прочитал в статье, расположенной на Вашем сайте, а конкретно // Весь массив из 826 чисел расшифровывается с помощью пяти массивов целых чисел:
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
где i - начальный номер в массиве, j - конечный номер в массиве, k - число коэффициентов аппроксимации, l интервал частной аппроксимации в днях внутри общего интервала, m - количество аппроксимируемых переменных.// Желательно бы взглянуть на подобную схему и для файла DE405. И потом, в связи с тем, что для Меркурия внутри общего интервала времени для данных существуют еще и частные интервалы, то получается, что x(t_c)=[t_c-0.5*(t_b+t_a)]/[0.5*(t_b-t_a)] надо рассчитывать по значениям t_a и t_b, которые будут равны времени на концах частного интервала. Или я не прав и их надо брать в начале и конце общего интервала, а коэффициенты использовать из частного интервала.

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




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

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

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