Тематические рабочие группы > Исследования астероидов и спутников планет
Уточнение орбит астероидов, что покрывают звезды
ctac:
Ну, что сравним ;) :
Semimajor axis = 3.15715683238201E+00 AU Semimajor axis = 3.15712194
Eccentricity = 0.215590453070217 Eccentricity = 0.21558403
Inclination = 1.3496186753636 deg Inclination = 1.349615
Long. of node = 40.0585571287539 deg Long. of node = 40.058537
Arg. of pericenter = 356.1177640071788 deg Arg. of pericenter =356.116158
Mean anomaly = 38.7063966102194 deg М. аномалия=43.132685
Epoch of elements : MJD 53674.81800000 TDT (Oct 31, 2005, 19.632000 h) эпоха=20051126
Как видно, очень сильна разница в средней аномалии.
Ошибка расчетов, или норма?
andyp:
Нельзя так сравнивать. Эпохи разные.
У меня Oct 31, 2005, 19.632000
У тебя 2005 11 26
Алексей Лосюк:
--- Цитата: ctac от 20 октября 2005, 22:58:05 ---
Интересно! Надо попробовать скачать astrob.dat и сравнить ;)
--- Конец цитаты ---
Лучше сравнивать не с astrob.dat (там данные очень приблизительные), а с уточнениями Престона, выполненными на основе "исторической" астрометрии. Вот они:
--- Цитировать --- * orbital information:
orbit source: Preston
date of fit: 2005 Sep 27
source of used astrometry: FASTT MPC (via Astdys)
number of used observations: 587
number of rejected observations: 26
time covered by the observations: 1931 Oct 17 - 2005 Jan 18
1 sigma uncertainty ellipse (major, minor, PA): 0.058, 0.041, 100°
remarks:
* orbital elements for (1253) Frisia :
Mean anomaly = 38.70651095 deg
Arg. of pericenter = 356.11797722 deg
Long. of node = 40.05830015 deg
Inclination = 1.34963804 deg
Eccentricity = 0.21559083
Semimajor axis = 3.15715610 AU
Perihelion dist = 2.47650221 AU
Mag: H = 11.96, G = 0.15
Epoch of elements : MJD 53674.81800000 TDT
(2005 Oct 31.818)
--- Конец цитаты ---
Но нагляднее для наших целей сравнить положение полосы покрытия, рассчитанное по данным Андрея Плеханова и Стива Престона. Я постараюсь сделать это сегодня вечером.
Андрей! А не могли бы Вы произвести ещё и такой эксперимент: рассчитать элементы орбиты этого астероида без использования "кадаровской" астрометрии - только на основании "исторический"? Это позволило бы сравнить Ваши результаты с результатами Престона, выполненными на одном и том же массиве наблюдательных данных.
И ещё вопрос. Для нас очень важно знать не только положение полосы покрытия, но и эллипс ошибок положения астероида (от него зависит вероятность покрытия в той или иной точке). В Ваших результатах параметров этого эллипса нет. Не знаете ли Вы, откуда их берёт Стив?
--- Цитата: ctac от 20 октября 2005, 22:58:05 ---
А отвергнуты 2 дня, как лишние?
--- Конец цитаты ---
На эту тему я пытался строить предположения, когда мы обнаружили, что Престон при предвычислении покрытия 5 октября Эпштейном не использовал "кадаровской" астрометрии. Может быть, программа OrbFit автоматически отвергает те наблюдения, результаты которых сильно отличаются от рассчитанных на основании предшествующих наблюдений?
Алексей Лосюк
andyp:
Результаты без Кадаровских измерений
Semimajor axis = 3.15715680741299E+00 AU
Eccentricity = 0.215590433904910
Inclination = 1.3496181841768 deg
Long. of node = 40.0585783661317 deg
Arg. of pericenter = 356.1177554771783 deg
Mean anomaly = 38.7063895487084 deg
Epoch of elements : MJD 53674.81800000 TDT (Oct 31, 2005, 19.632000 h)
--- Цитата: ctac от 20 октября 2005, 22:58:05 --- Интересно! Надо попробовать скачать astrob.dat и сравнить ;)
А отвергнуты 2 дня, как лишние?
--- Конец цитаты ---
Нет.
Слишком большая ошибка.
Вообще пока надо расматривать эти результаты не очень критично,
поскольку уточнения - это шаманство, которым я еще не обладаю.
Кроме того, AstOrb довольно запутанная программа с множеством
настроек и иногда неадекватным поведением. Документации тоже
практически нет, на вопросы разработчики не отвечают.
Поэтому, чтобы понять что там к чему надо смотреть код.
Фортрановский...
Андрей
Алексей Лосюк:
--- Цитата: andyp от 21 октября 2005, 13:47:08 ---
Вообще пока надо рассматривать эти результаты не очень критично,
поскольку уточнения - это шаманство, которым я еще не обладаю.
--- Конец цитаты ---
Вечером в пятницу, 21 октября, форум не работал, поэтому только сегодня мне удалось переписать элементы орбиты, рассчитанные Андреем.
Сравнение его элементов с элементами Престона показывают, что уточнение элементов орбиты - дело действительно субъективное (а, следовательно, и к результатам Престона следует относиться критически). :(
Вот цифры. Полоса, рассчитанная по элементам Андрея, полученным без учёта наблюдений на "Ка-Даре" (то есть основанная на тех же самых исходных данных, которыми пользовался Престон), проходит на 102 км южнее, чем у Престона. Полоса, вычисленная по элементам Андрея, полученным с учётом наблюдений на "Ка-Даре", проходит на 94 км южнее, чем у Престона.
Кроме того, по элементам орбиты Андрея покрытие получается на 15 секунд раньше, чем по престоновским элементам.
Предварительный практический вывод для наблюдателей покрытий, вероятно, должен быть таким: погрешности Престона следует увеличивать. Для покрытия астероидом (1253) Frisia звезды TYC 1374-00552-1 31 октября 2005 г. "сигма" положения полосы у Престона равнялась 110 км, а "сигма" по времени - 10 секундам. Умножим их хотя бы на 1.5...
Алексей
Навигация
Перейти к полной версии