Друк

[Гість]

Ласкаво просимо, Гість. Будь ласка, увійдіть або зареєструйтеся.
18 Квітня 2024, 09:21:38

1 година 1 день 1 тиждень 1 місяць Назавжди

Увійти

[Новини]

•    Новини
•      Астрономія
•      Астрозаходи
•    Життя клуба
•      Новини
•      Проекти
•      Заходи
•      Звіти про спостереження
•      Документи
•    Статті
•      Обладнання та телескопи
•      Спостереження
•      Астрофотографія
•      Астрософт
•      Астрономія в Україні
•      Астрономи шуткують
•      Інше
•     
•      Файловий архів

Наш канал

Наші сайти

Астро-сайти

Програми

Наш канал

Наші сайти

Астро-сайти

Програми

[uallxx]

Кто знает что такое "Besselian elements" и как по ним рассчитать параметры тени на поверхности Земли?

Например программа Occult 4.0.7.0 выдает следующие результаты:

2.7080  0.01472  0.33104 19.59958 -1.02606 -0.00405  0.00180

С такими пояснениями:

Besselian elements Field Definitions Content
1 UT at closest approach (decimal hours)
2 X coords of shadow at closest approach
3 Y coords of shadow at closest approach
4 hourly change in X
5 hourly change in Y
6 second-order hourly change in X
7 second-order hourly change in Y

А на сайте NASA можно узнать что

            Polynomial Besselian Elements for:   2009 Jul 22   03:00:00.0 TDT  (=t0)

              n        x          y         d          l1         l2          μ

              0   0.239961  -0.003276   20.26424   0.530427  -0.015633  223.38821
              1   0.5563955 -0.1774579  -0.007874  0.0000063  0.0000063  15.001004
              2  -0.0000576 -0.0001344  -0.000005 -0.0000128 -0.0000127
              3  -0.0000094  0.0000032

                             Tan ƒ1 = 0.0046014    Tan ƒ2 = 0.0045784

            At time t1 (decimal hours), each Besselian element is evaluated by:

               a = a0 + a1*t + a2*t^2 + a3*t^3    (or a = Σ [an*t^n]; n = 0 to 3)

                  where:    a = x, y, d, l1, l2, or μ
                            t = t1 - t0  (decimal hours) and t0 =  3.000 TDT


            The Besselian elements were derived from a least-squares fit to elements
            calculated at five uniformly spaced times over a six hour period centered at t0.
            The Besselian elements are valid over the period  0.00 ≤ t0 ≤  6.00 TDT.

            Note that all times are expressed in Terrestrial Dynamical Time (TDT).

Таким образом Occult выдал начальное время и коэффициенты квадратичного полинома для расчета X и Y. Но как их перевести в широту и долготу?

P.S.:
Вот такие результаты ожидаются на выходе:
(взято из kml, генерируемого программой Occult)

         -16.0000,  16.7458, 0
         -14.0000,  18.7030, 0
         -12.0000,  20.6610, 0
         -10.0000,  22.6122, 0
          -8.0000,  24.5494, 0
          -6.0000,  26.4652, 0
          -4.0000,  28.3527, 0
          -2.0000,  30.2053, 0
           0.0000,  32.0170, 0
           2.0000,  33.7825, 0
           4.0000,  35.4970, 0
           6.0000,  37.1567, 0
           8.0000,  38.7583, 0
          10.0000,  40.2995, 0
          12.0000,  41.7786, 0
          14.0000,  43.1944, 0
          16.0000,  44.5466, 0
          18.0000,  45.8353, 0
          20.0000,  47.0612, 0
          22.0000,  48.2250, 0
          24.0000,  49.3282, 0
          26.0000,  50.3723, 0
          28.0000,  51.3591, 0
          30.0000,  52.2904, 0
          32.0000,  53.1682, 0
          34.0000,  53.9944, 0
          36.0000,  54.7712, 0
          38.0000,  55.5005, 0
          40.0000,  56.1844, 0
          42.0000,  56.8247, 0
          44.0000,  57.4232, 0
          46.0000,  57.9818, 0
          48.0000,  58.5022, 0
          50.0000,  58.9858, 0
          52.0000,  59.4343, 0...

[IgorZ]

Похоже что это координаты в бесселевой системе, точнее аппроксимация полиномом движения центра тени в данной системе координат. Применяется при рассчете покрытий/затмений. Ось Z всегда параллельна направлению Луна - Солнце (для затмений), ось X лежит в экваториальной плоскости, Y дополняет систему до правой.
На калькуляторе это посчитать вряд ли получится.

[IgorZ]

Нашел описание в постоянной части АК, стр. 111. Еще у Монтенбрука была глава про затмения и покрытия, можете и там посмотреть.
А зачем Вам самим считать? Occult вроде бы и сам неплохо считает рисует.

[logrus]

Желающие разобраться с бесселевскими элементами могут почитать книжечку Text-Book on Spherical Astronomy:
http://books.google.com/books?id=I0k4AAAAIAAJ&pg=PA373&lpg=PA373

Содержание:
http://books.google.com/books?id=I0k4AAAAIAAJ&pg=PP10&lpg=PP1

Если отключить JavaScripts, то сканы можно сохранять.


А кому облом читать по страничке, могут найти эту книгу в формате DjVu на рапидшаре (11 мегабайт). Ссылку не даю, так как, кажется, это запрещено правилами форума.

[uallxx]

Нашел описание в постоянной части АК, стр. 111. Еще у Монтенбрука была глава про затмения и покрытия, можете и там посмотреть.
А зачем Вам самим считать? Occult вроде бы и сам неплохо считает рисует.

Спасибо за предоставленную информацию. Occult и вправду все неплохо делает, но иногда у него проскакивают цифры типа -336.9590 вместо 20.7757. И это только очевидные ошибки...

[IgorZ]

Пожалуйста Я знаю, что у него с точностью не очень, но чтобы так... Может бесселевы координаты тоже не лучше считаются?

[uallxx]

Это покажет дальнейшая работа (если мне поддастся алгоритм расчета тени).

[uallxx]

Кто сталкивался с эфемеридами из Лаборатории Реактивного Движения (JPL)? ftp://ssd.jpl.nasa.gov/pub/eph/export/ascii/
Читаю их инструкцию превращения свинца в золото convert the ASCII blocks into a binary ephemeris, а чуда не выходит.

Сказано сделать следующее:

An example for running asc2eph on a PC is the following:

   C:\> copy header.405+ascp2000.405  infile.405

   C:\> asc2eph < infile.405

Если первое проходит без затруднений, то при вызове "asc2eph < infile.405" читаем:

1) Name of the header file    = 2) Start of the julian date   = 3) End of the julian date     = 4) Name of binary output file =

Изменение команды на "asc2eph  < infile.405 2305424.50 2525008.50 out.405" ничего не меняет.