Астрономический форум ФПА "Ка-Дар"
Тематические рабочие группы => Исследования астероидов и спутников планет => Тема начата: Starcat от 03 декабря 2005, 19:51:47
-
Астероид Huruhata.
Вероятное покрытие - 26 января 2006 года 19:48.9 UT.
Звезда - TYC 2351-00042-1u, 11.m3 находится в Персее.
Продолжительность 6.5s
Ссылка на файл .rwo http://hamilton.dm.unipi.it/astdys/mpcobs/numbered/3/3939.rwo
Последние астрометрические наблюдения: 03.11.2005 Code: 704
Наблюдения на Терсколе 30.11.2005 в аттаче.
-
Наблюдения (03939) 6 декабря.
-
В Find_orb получены следующие элементы:
Orbital elements:
(3939)
Perihelion 2003 Dec 1.021839 TT
Epoch 2005 Dec 7.0 TT = JDT 2453711.5 Find_Orb
M 132.16738 (2000.0) P Q
n 0.17933691 Peri. 254.30156 0.24309537 0.94343150
a 3.1142699 Node 32.59862 -0.64928228 0.33096002
e 0.0979161 Incl. 24.74153 -0.72065050 0.02006153
P 5.50 H 11.6 G 0.15 q 2.8093326
From 204 observations 1982 Sept. 19-2005 Dec. 6; RMS error 0.389 arcseconds
.rwo файл в формате mpc и с последними наблюдения прилагается
-
По элементам орбиты астероида, взятым из файла astorb.dat от 2 декабря 2005 г., получается, что полоса покрытия пройдёт через центр Европы (карта прилагается).
Кстати, и через Бакырлытепе (Турция), где установлен полутораметровый российско-турецкий телескоп. С сотрудниками обсерватории контактирует Денис Денисенко. Возможно, он мог бы обратиться к ним с просьбой пронаблюдать покрытие...
Алексей Лосюк
-
Кстати, и через Бакырлытепе (Турция), где установлен полутораметровый российско-турецкий телескоп. С сотрудниками обсерватории контактирует Денис Денисенко. Возможно, он мог бы обратиться к ним с просьбой пронаблюдать покрытие...
Да, давно, я что-то Дениса не видел :( Просто теперь работаю в другом городе по средам и не могу ходить на собрания Московского Астрономического Клуба - ведь именно там мы чаще всего и пересекались. Надо бы письмице ему чиркнуть :)
-
Вот результаты вычисления элементов орбиты:
Semimajor axis = 3.11426988335481E+00 AU
Eccentricity = 0.097916548457838
Inclination = 24.7415480631589 deg
Long. of node = 32.5986098618435 deg
Arg. of pericenter = 254.3016678368601 deg
Mean anomaly = 132.1672897593710 deg
Epoch of elements : MJD 53711.00000000 TDT (Dec 7, 2005, 0.000000 h)
0.099" 0.090" 52.6
-
Вот результаты вычисления элементов орбиты:
<...>
Epoch of elements : MJD 53711.00000000 TDT (Dec 7, 2005, 0.000000 h)
0.099" 0.090" 52.6
Андрей!
У меня просьба, касающаяся как данного покрытия, так и любых других.
Не могли бы Вы рассчитывать параметры эллипса ошибок чуть иначе:
- во-первых, вычислять их именно на момент покрытия (чтобы мы могли оценить погрешность положения полосы именно рассматриваемого покрытия) - в данном случае на 2006 Jan 25, 19 h 53.8 m, то есть MJD=53760.82903;
- во-вторых, значения полуосей выводить с 4 цифрами после запятой (нужна "запасная" цифра при вычислениях). Динамический сайт астероидов даёт их c шестью значащими цифрами, например: 8.254091E-02 arcsec, 0.105676 arcsec.
Попутно замечу, что Максим ошибся в одной цифре, когда давал название этой теме: правильно будет "Уточнение орбиты астероида (3939) Huruhata (25.01.2006)".
-
Я ставлю ту же дату, которая использовалась в FindOrb для
сравнения результатов.
Больше количество цифр, увы, OrbFit не выводит.
Вот уточнения для момента покрытия:
Semimajor axis = 3.11413634253370E+00 AU
Eccentricity = 0.097958494963160
Inclination = 24.7416299334460 deg
Long. of node = 32.5987625215677 deg
Arg. of pericenter = 254.2884485838171 deg
Mean anomaly = 141.1157588755332 deg
Epoch of elements : MJD 53760.82903000 TDT (Jan 25, 2006, 19.896720 h)
0.082" 0.078"-152.2
-
Поправил название темы ;)
-
Я ставлю ту же дату, которая использовалась в FindOrb для сравнения результатов.
Понятно!
Но непосредственное сравнение значений элементов орбиты получается не очень наглядным. Изменение положения полосы покрытия, наверное, интереснее.
Вот уточнения для момента покрытия:
<…>
Спасибо! Буду теперь вычислять подробные обстоятельства покрытия - как обычно, в трёх вариантах: с данными из astorb.dat, от Find_orb и от OrbFit.
Больше количество цифр, увы, OrbFit не выводит.
Вот это сюрприз! Причём неприятный... :(
Мне-то казалось (по косвенным признакам), что в Динамический сайт астероидов "встроена" программа OrbFit. Выходит, что это не так: если и встроена, то версия, отличная от той, которой пользуется Андрей.
Кстати, если уж речь зашла об OrbFit. Андрей, давно хотел задать Вам вопрос: не попадался ли Вам дистрибутив "OrbFit2_3_win.zip" от 05.02.2004 г.? Похоже, это OrbFit для DOS и Windows. Мне этот архив дал Денис Денисенко. Можно ли его к чему-нибудь с пользой применить? :)
-
Этот дистрибутив можно взять сайта 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 даже отправил в названную рассылку красочные карты предстоящего покрытия (которые я здесь привожу), на которых отмечено положение потенциальных наблюдателей этого явления.
-
Я думаю, что в данном случае ошибся Стив.
Вот что он мне написал:
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.
Поэтому никто это покрытие и не наблюдал.
Андрей