Телескопы покупают здесь


A A A A Автор Тема: Эклиптические координаты Солнца (а может и Луны)  (Прочитано 2212 раз)

0 Пользователей и 1 Гость просматривают эту тему.

Оффлайн filiasfoxАвтор темы

  • Новичок
  • *
  • Сообщений: 37
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от filiasfox
Добрый день
Прошу поделиться алгоритмами или указать где есть достаточно точные формулы для определения эклиптических координат Солнца. Просмотрел форум. Увидел что ссылки идут в основном на "Астрономию с персональным компьютером" или на "Астрономию с калькулятором". К сожалению погрешность после составления софта составила около двадцати минут (например при определении весеннего равноденствия), что несколько великовато.

Оффлайн xd

  • *****
  • Сообщений: 17 977
  • Благодарностей: 378
    • Skype - deimos.belastro.net
  • Награды Открытие комет, астероидов, сверхновых звезд, научно значимые исследования.
    • Сообщения от xd
    • Белорусская любительская астрономическая сеть
К сожалению погрешность после составления софта составила около двадцати минут (например при определении весеннего равноденствия), что несколько великовато.
Каким образом это проверялось?
У природы нет плохой погоды, у неё просто на нас аллергия.

Учение без размышления бесполезно, но и размышление без учения опасно /Конфуций/
Слово есть поступок. /Л. Толстой/

Оффлайн filiasfoxАвтор темы

  • Новичок
  • *
  • Сообщений: 37
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от 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Автор темы

  • Новичок
  • *
  • Сообщений: 37
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от filiasfox
Каким образом это проверялось?

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

Оффлайн xd

  • *****
  • Сообщений: 17 977
  • Благодарностей: 378
    • Skype - deimos.belastro.net
  • Награды Открытие комет, астероидов, сверхновых звезд, научно значимые исследования.
    • Сообщения от xd
    • Белорусская любительская астрономическая сеть
Дайте больше информации: что вводилось, что было эталонными данными, как проверялось, что учитывалось, а что - нет.
20 минут - это время или угол? Какие использовались координаты: геоцентрические или топоцентрические?
У природы нет плохой погоды, у неё просто на нас аллергия.

Учение без размышления бесполезно, но и размышление без учения опасно /Конфуций/
Слово есть поступок. /Л. Толстой/

Оффлайн filiasfoxАвтор темы

  • Новичок
  • *
  • Сообщений: 37
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от filiasfox
Дайте больше информации: что вводилось, что было эталонными данными, как проверялось, что учитывалось, а что - нет.
20 минут - это время или угол? Какие использовались координаты: геоцентрические или топоцентрические?

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

Оффлайн xd

  • *****
  • Сообщений: 17 977
  • Благодарностей: 378
    • Skype - deimos.belastro.net
  • Награды Открытие комет, астероидов, сверхновых звезд, научно значимые исследования.
    • Сообщения от xd
    • Белорусская любительская астрономическая сеть
Окей, давайте по очереди. Берём дату 20.03.2017, какое время? Считаем MJD. Какая получилась дата?
Считаем положение Солнца. В какой системе координат? Барицентрическая, гелиоцентрическая, геоцентрическая, топоцентрическая?
У природы нет плохой погоды, у неё просто на нас аллергия.

Учение без размышления бесполезно, но и размышление без учения опасно /Конфуций/
Слово есть поступок. /Л. Толстой/

Оффлайн filiasfoxАвтор темы

  • Новичок
  • *
  • Сообщений: 37
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от filiasfox
Окей, давайте по очереди. Берём дату 20.03.2017, какое время? Считаем MJD. Какая получилась дата?
Считаем положение Солнца. В какой системе координат? Барицентрическая, гелиоцентрическая, геоцентрическая, топоцентрическая?

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




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

Оффлайн xd

  • *****
  • Сообщений: 17 977
  • Благодарностей: 378
    • Skype - deimos.belastro.net
  • Награды Открытие комет, астероидов, сверхновых звезд, научно значимые исследования.
    • Сообщения от xd
    • Белорусская любительская астрономическая сеть
А самое интересное?
У природы нет плохой погоды, у неё просто на нас аллергия.

Учение без размышления бесполезно, но и размышление без учения опасно /Конфуций/
Слово есть поступок. /Л. Толстой/

Оффлайн filiasfoxАвтор темы

  • Новичок
  • *
  • Сообщений: 37
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от filiasfox
А самое интересное?

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

Оффлайн xd

  • *****
  • Сообщений: 17 977
  • Благодарностей: 378
    • Skype - deimos.belastro.net
  • Награды Открытие комет, астероидов, сверхновых звезд, научно значимые исследования.
    • Сообщения от xd
    • Белорусская любительская астрономическая сеть
Ладно, прямой вопрос. Вы световую задержку не забыли?
У природы нет плохой погоды, у неё просто на нас аллергия.

Учение без размышления бесполезно, но и размышление без учения опасно /Конфуций/
Слово есть поступок. /Л. Толстой/

Оффлайн filiasfoxАвтор темы

  • Новичок
  • *
  • Сообщений: 37
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от filiasfox
Ладно, прямой вопрос. Вы световую задержку не забыли?

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

Оффлайн xd

  • *****
  • Сообщений: 17 977
  • Благодарностей: 378
    • Skype - deimos.belastro.net
  • Награды Открытие комет, астероидов, сверхновых звезд, научно значимые исследования.
    • Сообщения от xd
    • Белорусская любительская астрономическая сеть
Так собственно я и пытаюсь понять, где Вы столько ошибок наловили.
У природы нет плохой погоды, у неё просто на нас аллергия.

Учение без размышления бесполезно, но и размышление без учения опасно /Конфуций/
Слово есть поступок. /Л. Толстой/

Оффлайн filiasfoxАвтор темы

  • Новичок
  • *
  • Сообщений: 37
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от filiasfox
Так собственно я и пытаюсь понять, где Вы столько ошибок наловили.

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

Оффлайн xd

  • *****
  • Сообщений: 17 977
  • Благодарностей: 378
    • Skype - deimos.belastro.net
  • Награды Открытие комет, астероидов, сверхновых звезд, научно значимые исследования.
    • Сообщения от xd
    • Белорусская любительская астрономическая сеть
Исходники видел в сишной версии, но под рукой их нет. Если приложите к сообщению, можем посмотреть вместе.
У природы нет плохой погоды, у неё просто на нас аллергия.

Учение без размышления бесполезно, но и размышление без учения опасно /Конфуций/
Слово есть поступок. /Л. Толстой/

Оффлайн filiasfoxАвтор темы

  • Новичок
  • *
  • Сообщений: 37
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от filiasfox
Исходники видел в сишной версии, но под рукой их нет. Если приложите к сообщению, можем посмотреть вместе.

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

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



Оффлайн xd

  • *****
  • Сообщений: 17 977
  • Благодарностей: 378
    • Skype - deimos.belastro.net
  • Награды Открытие комет, астероидов, сверхновых звезд, научно значимые исследования.
    • Сообщения от xd
    • Белорусская любительская астрономическая сеть
Как я давно паскаля в глаза не видел...
У природы нет плохой погоды, у неё просто на нас аллергия.

Учение без размышления бесполезно, но и размышление без учения опасно /Конфуций/
Слово есть поступок. /Л. Толстой/

Оффлайн xd

  • *****
  • Сообщений: 17 977
  • Благодарностей: 378
    • Skype - deimos.belastro.net
  • Награды Открытие комет, астероидов, сверхновых звезд, научно значимые исследования.
    • Сообщения от 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Автор темы

  • Новичок
  • *
  • Сообщений: 37
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от 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

  • *****
  • Сообщений: 17 977
  • Благодарностей: 378
    • Skype - deimos.belastro.net
  • Награды Открытие комет, астероидов, сверхновых звезд, научно значимые исследования.
    • Сообщения от xd
    • Белорусская любительская астрономическая сеть
Кстати сишный код повторяет паскалевский полностью (я не про внешний вид, а про точность и пр.)?
А что, не должен?

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

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

Учение без размышления бесполезно, но и размышление без учения опасно /Конфуций/
Слово есть поступок. /Л. Толстой/