Выложил первую версию программы - 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)