Емельянов Эдуард Владимирович (eddy_em) wrote,
Емельянов Эдуард Владимирович
eddy_em

Categories:

Сравнение астрометрических библиотек - 2

Как я уже писал, получается, что разные библиотеки для вычисления apparent place дают разные данные. С расхождением в минуты!
Один из авторов любезно поделился со мной статьей, где таких расхождений нет: там говорится о 5-6mas, а не минутах!!!

Берем все тот же comparelibNOVA-NOVAS-SOFA-ERFA. Вводим туда координаты БТА, сегодняшнюю поправку DUT1 и смотрим, что поменяется, если в первом случае атмосферное давление сделать равным 100кПа, а во втором — 78кПа (обычное давление на горе). Давление как на уровне моря:
                                                          d J2000             d NOVAS          d SOFA Atco13
                  ra deg           ,dec deg          ra''      ,dec''    ra''      ,dec''    ra''      ,dec''
J2000             345.986294056000,28.150389198100
libNOVA equ_prec2 346.228458943796,28.258597189533   871.794,  389.549    -0.009,    0.039   324.298, -266.516
NOVAS precession  346.228461429350,28.258586387865   871.803,  389.510                       324.306, -266.555
SOFA Atci13       346.221443558472,28.258637867687   846.538,  389.695   -25.264,    0.185   299.042, -266.370
SOFA Atco13       346.138376291446,28.332629488577   547.496,  656.065  -324.306,  266.555
SOFA Atco13 NO r  346.221423118462,28.258608573980   846.465,  389.590   -25.338,    0.080   298.969, -266.475
SOFA Atco13 NO eo 345.886031202046,28.332629488577  -360.946,  656.065 -1232.749,  266.555  -908.442,    0.000
ERFA Atco13       346.138376291446,28.332629488577   547.496,  656.065  -324.306,  266.555     0.000,    0.000

Давление как на БТА:
                                                          d J2000             d NOVAS          d SOFA Atco13
                  ra deg           ,dec deg          ra''      ,dec''    ra''      ,dec''    ra''      ,dec''
J2000             345.986294056000,28.150389198100
libNOVA equ_prec2 346.228458943945,28.258597189600   871.794,  389.549    -0.009,    0.039   258.673, -208.067
NOVAS precession  346.228461429500,28.258586387931   871.803,  389.510                       258.682, -208.106
SOFA Atci13       346.221443558383,28.258637867562   846.538,  389.695   -25.264,    0.185   233.418, -207.921
SOFA Atco13       346.156605302540,28.316393568600   613.120,  597.616  -258.682,  208.106
SOFA Atco13 NO r  346.221423120371,28.258608573613   846.465,  389.590   -25.338,    0.080   233.344, -208.026
SOFA Atco13 NO eo 345.904260212808,28.316393568600  -295.322,  597.616 -1167.124,  208.106  -908.442,    0.000
ERFA Atco13       346.156605302540,28.316393568600   613.120,  597.616  -258.682,  208.106     0.000,    0.000

Запустил однострочником, чтобы как можно меньше было расхождение по JD.
В принципе, судя по коду, сравнение делается некорректно. В libnova вызывается функция ln_get_equ_prec2, никакого отношения к apparent place не имеющая, да и julian date не вычисляется нормально. В NOVAS julian date вычисляется, исходя из поправок для TT/TAI, не учитывая DUT1, но все равно идет просто пересчет координат на новую эпоху, без apparent place.

В общем, сравнивать их надо корректно. Но, судя по тому, что NOVA не использует DUT1, ее результаты таки должны быть менее правильными. Хотя, понятное дело: даже для БТА ошибка вычисления координат в ±1'' — мелочь по сравнению с ошибками наведения и сопровождения.
Tags: c, астрофизика
Subscribe

  • DHT22/DHT11 на STM32F103

    Добил шайтана! Сначала ожидал, что нужно будет полноценным захватом ШИМ пользоваться, но т.к. в протоколе неинформативная часть имеет постоянную…

  • Свеженькие железячки

    Получил сегодня с али ожидаемые железки, в т.ч. для восстановления моей файлопомойки. Во-первых, это блок питания на 16.8В, который я брал для…

  • BMP180 на STM32F103

    Добавил еще один сниппет — работа с BMP180 (датчик температуры и давления). Опять в даташите формула для вычисления "компенсированных" значений…

promo eddy_em august 17, 2019 12:33 3
Buy for 10 tokens
Юра намедни напечатал корпус для хронометра. Для первого блина получилось неплохо: И еще немного фотографий:
  • Post a new comment

    Error

    Anonymous comments are disabled in this journal

    default userpic

    Your reply will be screened

  • 0 comments