Расчёт текущих координат звёзд и прецессия.

[Николай15]

Расчитываю текущие координаты звёзд, имея их координаты в эпоху J2000.
Сейчас ограничился прецессией земной оси, собственные движения звёзд пока не трогаю.

Расчёт делаю по книге П.Даффета-Смитта "Практическая астрономия с калькулятором".
Пересчёт простой:

Ra_новое=Ra_J2000+(3.07420+1.33589*sin(Ra_J2000)*tg(Dec_J2000))*14.66
Dec_новое=Dec_J2000+(20.0383*cos(Ra_J2000))*14.66
где:
14.66 - прошло лет с 2000-го года.
остальное - константы для J2000.

Увы, не нашёл данных о точности этих формул.

Попробовал сравнить с текущими координатами из StarCalc - есть расхождения от нескольких секунд до минуты.
Причём, даже если учесть собственные движения звёзд, то они всё равно не покроют расхождение.
Исходные координаты J2000 в моём каталоге и в каталоге StarCalc почти одни и те же - расхождения на уровне десятых долей секунды.

(код проверил, вроде нигде не напутал...)

Собственно, 2 вопроса:

1) для проверки, нужны эталонные, заведомо правильные, текущие координаты c хорошей точностью. Конечно, лучше для звёзд с минимальными собственными движениями. Где взять? Или скинуть бы сюда, если есть.
(взять можно откуда угодно, но нужны именно надёжные)

2) может быть, эти формулы упрощены и чего-то не учитывают?

[xd]

Возьмите каталог Tycho-2, отфильтруйте записи с минимальным собственным движением - ищите среди них. "Астрономия на персональном компьютере" содержит формулы для прецессии и нутации.
Учтите, что на видимое положение накладывается ещё и световая аберрация.

[Николай15]

Я спрашивал о текущих координатах.
Tycho-2 был издан в 2000-ом году, 14 лет прошло, земная ось убежала...

"Астрономия на персональном компьютере" - Монтенбрук? Да, надо ещё у него посмотреть. Скорее всего нутация не учтена...

[xd]

"Астрономия на персональном компьютере" - Монтенбрук?

Да.

Я спрашивал о текущих координатах.
Tycho-2 был издан в 2000-ом году, 14 лет прошло, земная ось убежала...

Берёте координаты эпохи, редуцируете к текущей эпохе.

[andyp]

Возьмите каталог Tycho-2, отфильтруйте записи с минимальным собственным движением - ищите среди них. "Астрономия на персональном компьютере" содержит формулы для прецессии и нутации.
Учтите, что на видимое положение накладывается ещё и световая аберрация.

С нутацией там проблема.
Если нужны точные координаты,
посмотрите в моей программе LinOccult.

[Николай15]

Берёте координаты эпохи, редуцируете к текущей эпохе.

спасибо, за ответ конечно, но это не то - суть вопроса от Вас ускользнула...
Я сделал вычисления, теперь мне нужен эталон для проверки их правильности.
Нужен ГОТОВЫЙ эталон, а не мною вычисленный.

andyp, благодарю, сейчас поищу.

[Upsilon]

может быть, эти формулы упрощены и чего-то не учитывают?

Формулы использовались исторически, прецессию учитывают неполностью, для полной модели используют матрицу прецессии. На низких склонениях точность этих формул не хуже 1 секунды дуги. На высоких склонениях точность падает, потому что в формуле для прямого восхождения есть тангенс склонения, а во что превращается тангенс при приближении аргумента к 90 градусам - в школе учат.
В Старкалке учитывается и нутация, постоянная нутации 9.21", в этих пределах могут изменяться координаты вследствие нутации.
Постоянная аберрации 20.47", аберрацию тоже нужно учитывать.

[Николай15]

Да, спасибо, теперь понятнее.
Хотя, у меня ошибка намного больше девяти секунд..... странно, может быть в джаваскрпте тригонометрические функции огрублены....
Тангенс действительно глючит около Полярной.

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

[Upsilon]

Матрица прецессии:
http://www.astronet.ru/db/msg/1190817/node49.html
Матрица нутации (формулы 23,24 - ограниченное разложение достаточной точности):
http://law.rufox.ru/view/standarti/1672.htm
Формулы учета годичной аберрации низкой  (0.002") точности:
http://www.astronet.ru/db/msg/1190817/node41.html
Последовательность преобразований:
http://www.astronet.ru/db/msg/1190817/node46.html

[Николай15]

Upsilon, спасибо!
Бегло посмотрел - вроде нормально написано, но осталось найти время спокойно разобраться, переварить и закодировать...
время, как всегда время...