Эклиптические координаты Солнца (а может и Луны)

[filiasfox]

Добрый день
Прошу поделиться алгоритмами или указать где есть достаточно точные формулы для определения эклиптических координат Солнца. Просмотрел форум. Увидел что ссылки идут в основном на "Астрономию с персональным компьютером" или на "Астрономию с калькулятором". К сожалению погрешность после составления софта составила около двадцати минут (например при определении весеннего равноденствия), что несколько великовато.

[xd]

К сожалению погрешность после составления софта составила около двадцати минут (например при определении весеннего равноденствия), что несколько великовато.

Каким образом это проверялось?

[filiasfox]

Каким образом это проверялось?

Например на основании вот этих данных:
https://ru.wikipedia.org/wiki/%D0%A0%D0%B0%D0%B2%D0%BD%D0%BE%D0%B4%D0%B5%D0%BD%D1%81%D1%82%D0%B2%D0%B8%D0%B5

[filiasfox]

Каким образом это проверялось?

Сверкой результата.  Сверялось и поэтапно. Ошибка вылезает где-то в расчете эклиптических координат например в sun200 (астрономия с персональным компьютером, вариант для pascal)

[xd]

Дайте больше информации: что вводилось, что было эталонными данными, как проверялось, что учитывалось, а что - нет.
20 минут - это время или угол? Какие использовались координаты: геоцентрические или топоцентрические?

[filiasfox]

Дайте больше информации: что вводилось, что было эталонными данными, как проверялось, что учитывалось, а что - нет.
20 минут - это время или угол? Какие использовались координаты: геоцентрические или топоцентрические?

Дата: 20.03.2017.
Рассчитывалась модифицированная юлианская дата. На тоже время рассчитывались эклиптические координаты солнца согласно sun200 (собственно я писал, что рассчитывал эклиптические координаты). Учитывалось всё, что учитывалось в этом куске кода. Как то возмущения от венеры, марса, юпитера, сатурна. возмущения связанные с луной. 20 минут - это время. на выходе получалось время 10.06 вместо 10.28

[xd]

Окей, давайте по очереди. Берём дату 20.03.2017, какое время? Считаем MJD. Какая получилась дата?
Считаем положение Солнца. В какой системе координат? Барицентрическая, гелиоцентрическая, геоцентрическая, топоцентрическая?

[filiasfox]

Окей, давайте по очереди. Берём дату 20.03.2017, какое время? Считаем MJD. Какая получилась дата?
Считаем положение Солнца. В какой системе координат? Барицентрическая, гелиоцентрическая, геоцентрическая, топоцентрическая?

время 0 ч 0 м.
модифицированная юлианская дата = 57832
эклиптическая долгота - 358.7310854265




время 12 ч 0 м.
модифицированная юлианская дата = 57832.5
эклиптическая долгота - 359.21505689308

[xd]

А самое интересное?

[filiasfox]

А самое интересное?

А самое интересное, что при времени 12 часов - долгота должна быть явно в переделах 0-1 градус. что по расчетам согласно проге не выходит.

[xd]

Ладно, прямой вопрос. Вы световую задержку не забыли?

[filiasfox]

Ладно, прямой вопрос. Вы световую задержку не забыли?

А она по вашему составляет 20 минут?

[xd]

Так собственно я и пытаюсь понять, где Вы столько ошибок наловили.

[filiasfox]

Так собственно я и пытаюсь понять, где Вы столько ошибок наловили.

вы сам код из той книги (или с диска) видели? Может тогда стоит по нему пойти и разобраться  - где и что считается не так?
в личке вам отписал то что там получается согласно использованию формул без всех выкладок в коде. От кода пришлось отходить, поскольку (1) отсутствуют контрольные точки для проверки кода (2) перенос кода привел к ухудшению результатов расчетов

[xd]

Исходники видел в сишной версии, но под рукой их нет. Если приложите к сообщению, можем посмотреть вместе.

[filiasfox]

Исходники видел в сишной версии, но под рукой их нет. Если приложите к сообщению, можем посмотреть вместе.

Исходники прилагаются
параметр Т рассчитывается
  T:= (mjd-51544.5)/36525.0;
mjd=57832

При указании времени 10 ч 18 м 38.25 с (57832.429609375)  эклиптическая долгота получается практически равна нулю, что слегка навевает подозрения что что-то тут не так.

[xd]

Как я давно паскаля в глаза не видел...

[xd]

Окей, запишем для начала вот это:

Код: [Выделить]

(*-----------------------------------------------------------------------*)
(* SUN200: ecliptic coordinates L,B,R (in deg and AU) of the             *)
(*         Sun referred to the mean equinox of date                      *)
(*         (T: time in Julian centuries since J2000)                     *)
(*         (   = (JED-2451545.0)/36525             )                     *)
(*-----------------------------------------------------------------------*)

Как данная функция используется?

[filiasfox]

Окей, запишем для начала вот это:

Код: [Выделить]

(*-----------------------------------------------------------------------*)
(* SUN200: ecliptic coordinates L,B,R (in deg and AU) of the             *)
(*         Sun referred to the mean equinox of date                      *)
(*         (T: time in Julian centuries since J2000)                     *)
(*         (   = (JED-2451545.0)/36525             )                     *)
(*-----------------------------------------------------------------------*)

Как данная функция используется?

На вход отправляем T. На выходе получаем l,b,r.
Кстати сишный код повторяет паскалевский полностью (я не про внешний вид, а про точность и пр.)?

[xd]

Кстати сишный код повторяет паскалевский полностью (я не про внешний вид, а про точность и пр.)?

А что, не должен?

Собственно, в том издании, что у меня под руками (4-е) после описания данного алгоритма идут следующие главы:
6.4. Видимые и астрометрические координаты
 6.4.1 Аберрация и световая задержка
 6.4.2 Нутация

И дальше описание программы Planpos, где помимо указанных действий выполняется ещё много всего для редукций эпох.