Киевский клуб любителей астрономии "Астрополис"

astromagazin.net
* *
Добро пожаловать, Гость. Пожалуйста, войдите или зарегистрируйтесь.
15 Декабря 2017, 21:29:31


Автор Тема: Вопрос знатокам  (Прочитано 905 раз)

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

uallxx

  • Оффлайн Оффлайн
  • Сообщений: 82
  • Благодарностей: 0
Вопрос знатокам
« : 04 Ноября 2009, 10:53:48 »
. - .

Кто знает что такое "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

  • Оффлайн Оффлайн
  • Сообщений: 336
  • Благодарностей: 0
Re: Вопрос знатокам
« Ответ #1 : 04 Ноября 2009, 11:21:10 »
. - .

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

IgorZ

  • Оффлайн Оффлайн
  • Сообщений: 336
  • Благодарностей: 0
Re: Вопрос знатокам
« Ответ #2 : 04 Ноября 2009, 11:35:41 »
. - .

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

logrus

  • Оффлайн Оффлайн
  • Сообщений: 679
  • Благодарностей: 6
  • Максим
Re: Вопрос знатокам
« Ответ #3 : 04 Ноября 2009, 13:41:02 »
. - .

Желающие разобраться с бесселевскими элементами могут почитать книжечку 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

  • Оффлайн Оффлайн
  • Сообщений: 82
  • Благодарностей: 0
Re: Вопрос знатокам
« Ответ #4 : 04 Ноября 2009, 18:32:31 »
. - .

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

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

IgorZ

  • Оффлайн Оффлайн
  • Сообщений: 336
  • Благодарностей: 0
Re: Вопрос знатокам
« Ответ #5 : 04 Ноября 2009, 19:04:03 »
. - .

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

uallxx

  • Оффлайн Оффлайн
  • Сообщений: 82
  • Благодарностей: 0
Re: Вопрос знатокам
« Ответ #6 : 04 Ноября 2009, 19:19:19 »
. - .

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

uallxx

  • Оффлайн Оффлайн
  • Сообщений: 82
  • Благодарностей: 0
Re: Вопрос знатокам
« Ответ #7 : 08 Ноября 2009, 18:39:40 »
. - .

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

Сказано сделать следующее:
Цитировать
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" ничего не меняет.
Записан