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


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

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

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

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

Planpos в данный момент пока не интересует, а вот куда выпадают несколько минут - вполне интересно

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

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

Ну в общем пару ошибок у себя нашел. Но Как бы всё равно нестыковка великовата слишком.
2017 3 20 Весеннее равноденствие 10 ч 18 м
2017 6 21 Летнее солнцестояние 4 ч 12 м
2017 9 22 Осеннее равноденствие 19 ч 50 м
2017 12 21 Зимнее солнцестояние 16 ч 16 м

Около десяти минут. явно больше, чем световая задержка

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

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

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

Оффлайн xd

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

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

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

  • Новичок
  • *
  • Сообщений: 37
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от filiasfox
Если с солнцестояниями и равноденствиями не так всё плохо,
В видимые экваториальные координаты эпохи объекта перевели? Уравнение времени учли? Радиус Солнца учли? Атмосферную рефракцию учли?

Брался модуль sunset. Там все учитывалось. В том числе и перевод в экваториальные координаты. Единственное отличие от того,  что там использовалось - для расчета эклиптических координат ипсользовался присланный вам код, а не minisun и не формулы, представленные в книге. При использовании формул из книги - точность была на уровне 5-10 минут или и того выше (но для восхода и захода, а не равноденствий и солнцестояний). Но честно говоря не хотелось бы использовать разные расчеты координат  для разных случаев, как делается в той книге. Это несколько абсурдно. И лля каждого из случаев они дают разную точность.

Я читаю, что вы пишите. Но после вашего предложения протестировать или проверить кусок кода - с вашей стороны тоже действий не было видно.

Готовый ответ получить конечно проще всего. Но мне хочется понять что и где делается и что именно делаю не так. В частности есть подозрения что перевод из эклиптических в экваториальные координаты происходит не совсем корректно. Хотя ошибка может пребывать и в попытке определить время восхода и захода для гражданских сумерек. Вот завтра и займусь проверкой перевода одних координат в другие. Хотя когда нет контрольных точек для расчёта чего-либо... ковыряться в коде - слегка стрёмно
« Последнее редактирование: 03 Мая 2017 [21:52:16] от filiasfox »

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

  • Новичок
  • *
  • Сообщений: 37
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от filiasfox
Возвращаемся к вопросу темы
Для расчётов координат используется кусок кода из Монтенбрука (все процедуры и функции перенесены в прогу из исходников. один к одному):
  DT:= (8.32/1440.0)/36525;
  SUN200(T-DT,L,B,R);          (* geocentric ecliptic coordinates      *)
  CART(R,B,L,X,Y,Z);           (* cartesian ecliptic coordinates       *)
  ECLEQU(T,X,Y,Z);             (* equatorial ecliptic coordinates      *)
  NUTEQU(T,X,Y,Z);             (* correction for nutation              *)
  POLAR(X,Y,Z,R,DEC,RA);       (* spherical coordinates;               *)

При попытке с помощью него найти время таких событий как равноденствия, солнцестояния  - получаем некрасивую картинку, где весеннее равноденствие не совпадает (если сверять с этими данными https://ru.wikipedia.org/wiki/%D0%A1%D0%BE%D0%BB%D0%BD%D1%86%D0%B5%D1%81%D1%82%D0%BE%D1%8F%D0%BD%D0%B8%D0%B5) на 2-3 минуты (по времени). Что касается хотя бы летнего солнцестояния, то ошибка идёт от двух минут до 12-13 минут.

Для расчёта эпохи используется модифицированная юлианская дата с собственно определение
  t:=(mjd_start-51544.5)/36525;

Где следует искать ошибку? Контрольных точек для проверки расчётов в исходниках и книге как-то не обнаружено

Оффлайн AN14440

  • ****
  • Сообщений: 291
  • Благодарностей: 16
    • Сообщения от AN14440
15 лет назад: Как рассчитываются
Там первой на вопрос дилетанта откликнулась Karen, у которой сегодня кстати день рождения, а потом долго и упорно всё разжёвывал, не теряя терпения, Анатолий Волчков, почётный член форума, его к сожалению с нами уже нет.

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

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

Спасибо. Посмотрел ответы Анатолия Волчкова. Там также делается сноска на Монтенбрука (в простейшем варианте расчёта координат солнца mini_sun). Но в описании минисан есть один момент. там прописана формула для определения T=(JD-2451545)/36525 ). Если ставить именно эту формулу, то расхождение будет несколько часов при определении тех же равноденствий.

При этом у Монтенбрука вызов данной процедуры осуществляется с параметром T=(JD-51544.5)/36525 ). Что дает расхождение в несколько минут. Вот это-то и смущает. То есть сырцы брались один к одному, но значения....

Кстати та ветка была в основном про время восхода и захода. Это в данный момент мне не совсем актуально.

Оффлайн xd

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

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

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

  • Новичок
  • *
  • Сообщений: 37
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от filiasfox
Это одна и та же формула, только во второй вместо JD должно быть MJD.
MJD = JD - 2400000.5
Формулы полностью эквивалентны.
JD считаются от полудня, MJD - от полуночи.

Это знаю. Используется mjd.
Соответственно используется формула для расчёта T=(mjd-51544.5)/36525;
Нестыковки от 2 до десяти минут в разные стороны при всех расчётах





Оффлайн xd

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

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

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

  • Новичок
  • *
  • Сообщений: 37
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от filiasfox
Тогда расскажите, что учитывается в процессе расчёта.
1. Какая система отсчёта используется в модели?
2. Какая шкала времени используется?
3. Какого порядка возмущения от планет учитывается? По каким координатам?
4. На какой промежуток времени выполняются расчёт и как ведут себя отклонения на краю интервала?
5. Расчёт геоцентрический или топоцентрический?
ну и так далее

Сырцы для расчёта экваториальных координат - в SUNEQU.PAS у Монтенбрука. Вы писали что у вас возможно есть исходник на Си (или c++). Там различий от версии на pascal почти нет.


(*-----------------------------------------------------------------------*)
(* SUNEQU: apparent equatorial coordinates of the Sun                    *)
(*         (right ascension RA, declination DEC in deg, R in AU)         *)
(*         T in Julian centuries since J2000 (T=(JD-2451545.0)/36525)    *)
(*-----------------------------------------------------------------------*)
PROCEDURE SUNEQU(T:REAL;VAR RA,DEC,R:REAL);
  VAR DT,L,B,X,Y,Z: REAL;
  BEGIN
    DT := (8.32/1440.0)/36525.0; (* light-time correction of 8.32 min    *)
    SUN200(T-DT,L,B,R);          (* geocentric ecliptic coordinates      *)
    CART(R,B,L,X,Y,Z);           (* cartesian ecliptic coordinates       *)
    ECLEQU(T,X,Y,Z);             (* equatorial ecliptic coordinates      *)
    NUTEQU(T,X,Y,Z);             (* correction for nutation              *)
    POLAR(X,Y,Z,R,DEC,RA);       (* spherical coordinates;               *)
  END;                           (* true equinox of date                 *)


Промежуток времени: 2017-2024
Контрольные значения для равноденствий и солнцестояний были взяты из таблицы:
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

Расхождения - по весенним равноденствиям около двух минус
По летним солнцестояниям - от 2 до 15 минут
по осенним равноденствиям около двух минус
По зимнему солнцестоянию - разница доходила до 10-15 минут.



Оффлайн xd

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

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

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

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

Я сказал что авторы там дураки? Если память не изменяет, то нет. Что оттуда использовал - также сказал. Вы можете помочь пройтись по контрольным точкам или только впустую поговорить хотите, упрекая собеседника в том, что он не говорил и не делал?

Могу сырец полностью сбросить. Но вроде в личку вы не отвечаете


Оффлайн xd

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

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

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

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

Первое. Но желательно с указанием на последовательности формул и пр. для вычислений.
Что касается второго, то за основу была взята работа Монтенбрука "Астрономия с персональным компьютером". Результаты вам уже сообщил.

В той ветке, на которую мне перед этим дали ссылку, была ссылка как раз на кусок из исходников Монтенбрука "Астрономия с персональным компьютером". В частности модуль mini_sun.
Теоретически sunequ выдает координаты точнее, поскольку там учитываются погрешности из-за разных планет и пр..

Собственно вопрос и был в том - где искать ошибку. Часть неких контрольных точек для расчётов я привел. Но на данный момент помощи в поиске пока не вижу

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

  • Новичок
  • *
  • Сообщений: 37
  • Благодарностей: 0
  • Мне нравится этот форум!
    • Сообщения от filiasfox
Ну вот в частности то, что получилось
Данные по Вики
2017    10:28    04:24    20:02    16:28
2018    16:15    10:07    01:54    22:23

Рассчитано по Монтенбруку
2017   10.30   04.22   20.04   16.39
2018   16.17   10.12   01.55   22.18

Вот из-за подобных расхождений и пытаюсь разобраться что и где не так.

Оффлайн xd

  • *****
  • Сообщений: 17 977
  • Благодарностей: 378
    • Skype - deimos.belastro.net
  • Награды Открытие комет, астероидов, сверхновых звезд, научно значимые исследования.
    • Сообщения от xd
    • Белорусская любительская астрономическая сеть
Прежде, чем браться за алгоритмы и ошибки, стоит ознакомиться с матчастью. Для этого хорошо подходит "Астрономический календарь. Постоянная часть."
В той ветке, на которую мне перед этим дали ссылку, была ссылка как раз на кусок из исходников Монтенбрука "Астрономия с персональным компьютером". В частности модуль mini_sun.
Теоретически sunequ выдает координаты точнее, поскольку там учитываются погрешности из-за разных планет и пр..
Вот собственно я и просил показать, где это учитывается. Очевидно, в функции SUN200, но хотелось бы это проверить.

Я в исходнике вижу приведение к эпохе объекта и учёт нутации.
У природы нет плохой погоды, у неё просто на нас аллергия.

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

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

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




Прежде, чем браться за алгоритмы и ошибки, стоит ознакомиться с матчастью. Для этого хорошо подходит
Вот собственно я и просил показать, где это учитывается. Очевидно, в функции SUN200, но хотелось бы это проверить.
Я в исходнике вижу приведение к эпохе объекта и учёт нутации.

Оффлайн xd

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

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