По этой скорости я не находил не одной ссыл ки, где бы данные отличались максимум на пару процентов. А вот по скоростям галактики и местной группы галактик очень разнятся в разных литературных источниках, но нигде не указ ывается, что ско

Скачал файлы с данными JPL (конкретно DE200 и DE405), но никак не могу разобраться как их использовать для вычисления координат планет и Солнца в своей программе, т.к. использовать уже готовые программы я не могу, хотя и скачал одну на С+. То что все данные записаны блоками с данными на 32 дня и каждая запись составляет 826 или 1652 числа по 8 байт я понял. Первые два числа это начало и конец периода в 32 дня в Юлианском исчисление, а далее должны идти коэффициенты для полиномов Чебышева, но не ясно только для координат или и для скоростей планет тоже, т.к. в разных местах написано по разному. Ясно и то, что для разных планет эти коэффициенты расположены в строго определенном месте. Меркурий с 3-ей позиции по 146, Венера с 147 по 182 и т.д. Для DE200 у меня есть расшифровка не только этих позиций, но и дано число коэффициентов аппроксимации, частный интервал аппроксимации в днях (для Меркурия он равен 8 дней) и количество аппроксимируемых данных (для всех планет указано 3, т.е. получается, что коэффициенты даны только для координат), а как это будет для DE405 найти не могу. А если даже и найду, то не знаю что за коэффициенты я считаю и куда их потом подставлять. Если кто-то в курсе, помогите.

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

Вадим, большое спасибо. Сильно меня выручили с полиномами Чебышева, но перейти к вычислениям я пока не могу, т.к. остались не выясненными некоторые моменты. Например, Вы для вычисления и координат и скоростей записали в формулах одни и те же коэффициенты 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, которые будут равны времени на концах частного интервала. Или я не прав и их надо брать в начале и конце общего интервала, а коэффициенты использовать из частного интервала.

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