 |
Астронет: В. В. Нестеров/ГАИШ Стандарт основных вычислений астрономии http://www.astronet.ru/db/msg/1178709/node12.html |
1.8. Запинка с преобразованием
Все вычисления должны быть стандартными и использовать выражения для прямоугольных координат.
Сферические координаты объектов 
, 
, 
,
,
, 
можно преобразовать к
декартовым координатам 
, 
, 
, 
, 
, 
по формулам
Здесь , - не обязательно экваториальные координаты, это могут быть и
эклиптические долгота и широта.
Расстояние до звeзды может быть вычислено при помощи параллакса или приравнено единице, если
параллакс не известен. Второй вариант используется для экваториальных координат Опорного каталога
Тихо. В этом случае шестимерный вектор принимает вид , , 
,
,
, 0.
Обратное преобразование
Когда объект находится в полюсе системы, его координата не определена,
при составлении программы ей можно дать любое значение.
Еще один входной параметр - это время, может быть выражено в самых разнообразных формах. Наиболее удобная форма - сквозная нумерация суток, предложена Скалигером. Счет дат идeт от 1 января 4713 года до нашей эры. Юлианские сутки начинаются в полдень.
В спутниковой геодезии употребляется модифицированный юлианский день, MJD=JD-2400000.5, начало суток смещено в полночь.
Запишем алгоритм, который по заданной дате григорианского
календаря позволяет вычислять юлианскую дату. Для этого введем
обозначения 
, 
, 
для дня, месяца и года. Будем
использовать функцию
для обозначения целой
части числа .
Если 
то 
и
, иначе 
, 
.
Вспомогательные величины
- целое число
столетий и
- учет того, что последний
год только одного столетия из четырех окажется високосным.
Юлианская дата равна
Если 
то 
иначе
.
Если 
то
, иначе
.
<< 1.7. Затраты огромны | Оглавление | 1.9. Оттуда и сюда >>