Тематические рабочие группы > Исследования астероидов и спутников планет

(3939) Huruhata

<< < (3/3)

andyp:
Этот дистрибутив можно взять сайта AstDys, только он без танцев с бубном
работать не будет. Я уже в этом убедился.

Алексей Лосюк:
Стив Престон 22 января 2006 г. рассчитал элементы орбиты Huruhata с использованием результатов наблюдений Максима Андреева, выполненных им в ноябре и декабре 2005 г. Вот что получилось у Престона (http://www.asteroidoccultation.com/2006_01/0125_3939_7373.htm ):


 * orbital elements for (3939) Huruhata :
      Mean anomaly       = 141.04669536 deg
      Arg. of pericenter = 254.28849790 deg
      Long. of node      =  32.59872853 deg
      Inclination        =  24.74160342 deg
      Eccentricity       =   0.09795809
      Semimajor axis     =   3.11413760 AU
      Perihelion dist    =   2.80908262 AU
      Mag:             H = 11.04,  G = 0.15
      Epoch of elements  : MJD   53760.44400000 TDT
                           (2006 Jan 25.444)

Я рассчитал положение полосы покрытия по четырём наборам элементов орбиты:

      1 - из файла astorb.dat от 02 декабря 2005 г. (эллипс ошибок - с сайта http://hamilton.dm.unipi.it/ )
      2 - рассчитаны Максимом Андреевым по программе Find_orb (http://www.kadar.ru/forum/index.php/topic,204.msg1303.html#msg1303 )
      3 - рассчитаны Андреем Плехановым по программе OrbFit (http://www.kadar.ru/forum/index.php/topic,204.msg1484.html#msg1484 )
      4 - рассчитаны Стивом Престоном

Давайте сравним промежуточные и окончательные результаты.

Эллипс ошибок астероида - большая полуось ошибок, малая полуось ошибок, позиционный угол большой полуоси (mas, mas, градусы):

      1) 92.21782, 86.30528, -144.5459116
      2) 389, 389, --
      3) 82, 78, -152.2
      4) 77, 76, 175=-185

Положение полос относительно полосы, рассчитанной по элементам орбиты из astorb (км):

      1) 0
      2) 619 км к югу
      3) 137 км к югу
      4) 270 км к югу

Погрешность положения полосы (1 "сигма", км):

      1) 204
      2) 852
      3) 183
      4) 172

Смещение моментов наступления середины покрытия относительно момента, рассчитанного по элементам орбиты из astorb (с):

      1) 0
      2) +9
      3) -11
      4) -14

Погрешность по времени (3 "сигмы", с):

      1) 121
      2) 532
      3) 109
      4) 106

Вероятность покрытия в середине полосы (%):

      1) 6.1
      2) 1.5
      3) 6.8
      4) 7.2

Какие можно сделать выводы? К сожалению, никто не производил контрольного вычисления элементов орбиты без учёта наблюдений в Терсколе. Поэтому невозможно сказать, насколько возросла вероятность покрытия в середине полосы в результате новых наблюдений.

Что касается сравнения расчётов, то программа Find_Orb опять продемонстрировала, что не может конкурировать с OrbFit по точности. Или, как вариант, Find_Orb даёт объективную информацию о погрешности элементов орбиты, а OrbFit её бессовестно занижает. :)

Результаты, полученные с помощью OrbFit Андреем Плехановым и Стивом Престоном, опять (второй раз из трёх проведённых попыток сравнения - смотрите уточнение обстоятельств покрытия 2006.01.06 и 2006.01.08 астероидами Pushkin и Zhangguoxi) заметно различаются. Расхождение в положении полосы на 133 км (при ширине полосы 31 км) - это многовато. :(

Ну и, наконец, о прикладном значении проведённой работы: собственно о подготовке к наблюдениям этого покрытия.

После того как Стив Престон выложил на свой сайт результаты предвычислений покрытия, в среде европейских наблюдателей покрытий обнаружился явный интерес к нему (ведь полоса проходит через центр Европы) - сужу об этом по сообщениям в рассылку [PLANOCCULT]. Европейские наблюдатели сразу отметили, что покрытие будет происходить в то время, когда они будут наблюдать покрытие звезды 7.9 m Сатурном. И начали прикидывать, как бы переложить телескопы на две-три минуты для наблюдения покрытия астероидом, чтобы затем продолжить наблюдения Сатурна.

Oliver Kloes даже отправил в названную рассылку красочные карты предстоящего покрытия (которые я здесь привожу), на которых отмечено положение потенциальных наблюдателей этого явления.

andyp:
Я думаю, что в данном случае ошибся Стив.
Вот что он мне написал:

You found a bug in my code :)  Thanks...

Obscode 809 is La Silla. This is one of the sources that I treat as "high
accuracy".  But my code should only reduce the rms for specific instances
and it was reducing the rms for ALL La Silla observations... Opps!

I have fixed this now so it shouldn't happen in future orbit fits.

Поэтому никто это покрытие и не наблюдал.

Андрей

Навигация

[0] Главная страница сообщений

[*] Предыдущая страница

Перейти к полной версии