 |
Астронет: В. В. Нестеров/ГАИШ Основные алгоритмы спутниковой геодинамики http://variable-stars.ru/db/msg/1178726/node15.html |
1.11. Без эфемериды - в никуда
Стандарт МАС1976 рекомендует положения Солнца, Луны и больших планет вычислять на основе эфемериды DE200/ LE200.
Эфемеридные данные DE200/LE200, распространяемые на магнитных носителях, состоят из записей длиной 6608 байт. Каждая запись содержит 826 чисел двойной точности в двоичном формате по 8 байт. В одной записи упакована информация о положениях и скоростях небесных объектов на интервале времени, равном 32 суткам. Первое число - юлианская дата начала интервала прогнозирования, второе число - конечная юлианская дата данного интервала, образованная прибавлением 32 суток к начальной дате. Положения и скорости планет и Солнца даны относительно барицентра Солнечной системы. Вместо координат Земли упакованы коэффициенты для вычисления положения центра масс системы Земля-Луна. Положение Луны дано относительно центра Земли. Весь массив из 826 чисел расшифровывается с помощью пяти массивов целых чисел:
| объект | 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 |
где
- начальный номер в массиве,
- конечный номер в массиве,
- число коэффициентов аппроксимации,
- интервал частной аппроксимации в днях внутри общего
интервала,
- количество аппроксимируемых переменных.
Рассмотрим, для примера, алгоритм вычисления положения Венеры относительно центра Земли.
Пусть 
- дата в юлианских днях в шкале динамического барицентрического времени TDB,
выбранная для вычислений, 
, 
- начальная и конечная даты эфемеридных данных. Прежде
всего определим порядковый номер того тридцатидвухсуточного интервала, в который попадает
выбранный нами момент
:
целому значению выражения 
и центра
масс Земля-Луна 
относительно барицентра Солнечной системы и положение Луны
относительно центра Земли. Вычисляем вектор положения Земли
относительно
барицентра
равен отношению массы Луны к массе Земли, а затем этот вектор
вычитаем из соответствующего вектора для планеты Венера. Если эти же расчеты проводятся с учетом
скорости распространения света, то положение Венеры относительно барицентра вычисляется в момент
времени, предшествующий данному и определяемый двумя-тремя итерациями.
<< 1.10. Начальные условия или ... | Оглавление | 1.12. Цели ясны, нужны указания >>