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

args[0]=message
args[1]=DB::DB::Message=HASH(0x31355c0)
Re[20]: Смещение перигелия Меркурия и других планет
17.01.2008 17:04 | В. В. Чазов

Дорогие друзья, продолжим.
На момент T_{eph}=2454528.25075444 надо вычислить положение Меркурия относительно Солнца.
Среди текстовых файлов эфемериды DE200/LE200 находим файл ascp2000.200. В этом файле есть массив из 826 чисел, необходимых для нашего момента времени.
Действительно, этот массив расположен под номером 94, начинается с юлианской даты 2454512.5, работает до юлианской даты 2454544.5.
Для Меркурия разбиваем интервал на четыре части, по 8 дней:
1. 2454512.5 - 2454520.5,
2. 2454520.5 - 2454528.5,
3. 2454528.5 - 2454536.5,
4. 2454536.5 - 2454544.5,
Наш момент t_c=2454528.25075444 попадает во второй интервал, в таком случае
(t_a=2454520.5) < t_c < t_b=2454528.5.
Вычисляем номер ячейки в массиве из 826 чисел, с которой начинаются коэффициенты аппроксимации для второго интервала.
Начальный номер i=3.
Число коэффициентов для каждой переменной k=12.
Число переменных m=3.
Общее число ячеек для одного подинтервала равно (k*m)=36. Второй интервал начинается с элемента массива номер 3+36=39.
Первые двенадцать чисел - коэффициенты
a_0, a_1, a_2, ... , a_11 для координаты x.
Элементы с номера 51 по номер 62 - для координаты y.
Элементы с номера 63 по номер 74 - для координаты z.
Преобразуем момент времени t_c к интервалу изменения [-1,+1]:
x(t_c)=[t_c-0.5*(t_b+t_a)]/[0.5*(t_b-t_a)],
по рекуррентному алгоритму вычисляем числовые значения полиномов Чебышева от T_0(x(t_c)) до T_{11}(xt(_c)),
с помощью коэффициентов аппроксимации вычисляем по отдельности три координаты Меркурия на заданный момент.
Но мы получим положение Меркурия относительно барицентра Солнечной системы относительно неподвижного экватора.
Для вычисления положения Меркурия относительно Солнца надо вычислить координаты Солнца относительно барицентра. Вычисления выполняются по тому же алгоритму, только номера ячеек в массиве другие.
Из каждой координаты Меркурия надо вычесть соответствующую координату Солнца.
Для перехода из плоскости экватора в плоскость эклиптики следут выполнить поворот против часовой стрелки на угол \varepsilon - угол наклона эклиптики к экватору.
С поклоном, Ваш Вадим.


Форумы >> Астрономия и Интернет
Список  /  Дерево
Заголовки  /  Аннотации  /  Текст