Тематические рабочие группы > Вопросы астрометрии
OrbSearch - приложение для расчета и уточнения элементов орбит астероидов
OMI:
Кто там дела обстоят с данным проектом?
Программу ещё не сделали?
А то у меня есть интерес к данному вопросу, в части алгоритмов.
andyp:
Основной каркас программы создан.
Реализован модуль вычисления предварительной
орбиты по трем точкам.
Результат совпадает с OrbFit.
andyp:
Выложил первую версию программы - OrbSearch 1.0.0 Beta.
Исходники можно найти на моей страничке:
http://andyplekhanov.narod.ru
и на
http://sourceforge.net/projects/orbsearch
Начальная орбита вычисляется с помощью метода Гаусса по трем наблюдениям,
а затем уточняется с помощью нахождения минимума функции отклонения Simplex методом.
В отличии от известных программ упор делается на статистическую обработку и вычисление
неопределености обриты путем вычисления нескольких вариантов.
В программе используются библиотеки из LinOccult и GNU scientific library.
Вот несколько сравнений результатов с OrbFit:
Астероид 49:
OrbSearch (best fitted):
Epoch (y m d) 2008/05/25.25 (54611.25)
Semi-major axis (a) 3.09061843 AU
Eccentricity (e) 0.23131746
Inclination (i) 3.18124886 deg
Long. of ascending node 286.13984314 deg
Argument of perihelion 110.03371310 deg
Mean anomaly (M) 240.29977619 deg
OrbFit:
Semimajor axis = 3.09061933737767E+00 AU
Eccentricity = 0.231307001433158
Inclination = 3.1814838945372 deg
Long. of node = 286.1301224096144 deg
Arg. of pericenter = 110.0438989194008 deg
Mean anomaly = 240.3045776164618 deg
Epoch of elements : MJD 54611.28400000 TDT (May 25, 2008, 6.816000
h)
Астероид 165:
OrbSearch (best fitted):
Epoch (y m d) 2008/05/20.00 (54606.00)
Semi-major axis (a) 3.12572643 AU
Eccentricity (e) 0.08338424
Inclination (i) 11.23946034 deg
Long. of ascending node 302.65132928 deg
Argument of perihelion 350.83658219 deg
Mean anomaly (M) 308.39276761 deg
OrbFit:
Semimajor axis = 3.12571555818036E+00 AU
Eccentricity = 0.083388649509875
Inclination = 11.2393034099170 deg
Long. of node = 302.6507359851024 deg
Arg. of pericenter = 350.8433943595194 deg
Mean anomaly = 308.4035998068666 deg
Epoch of elements : MJD 54606.09000000 TDT (May 20, 2008, 2.160000
h)
Lupus:
OrbSearch портирована и под Windows. Подробности ТУТ. Тестирование приветствуется. :)
Lupus:
Для расчёта орбит конкретных астероидов нужны файлы с данными их наблюдений. Но...
--- Цитировать ---The .obs files, in the format specified by the Minor Planets Center, are currently available from the MPC Extended Computer Service at http://cfa-www.harvard.edu/iau/services/WebCSAccess.html. However, note that at present this is a service reserved to paying customers.
--- Конец цитаты ---
Буржуины проклятые... :'( Закрывают доступ к научной информации, а потом деньгу сшибают... Это что ж -- чтобы получить эти данные, надо теперь региться к ним за деньги?!
Навигация
Перейти к полной версии