Как построить орбиту двойной звезды и рассчитать элементы ее орбиты?

У меня есть набор дат, углов положения ( θ ) и угловых разносов ( р ) для визуальной двойной звезды. Например:

1994.3132  57.0    0.097
1996.1783  73.4    0.051
2006.1938 182.0    0.126
2009.2597 192.9    0.152
etc.

Каковы алгоритмы расчета элементов орбиты и построения орбиты двойной звезды?

Существуют ли доступные коды (или программы) для этих расчетов?

Будет ли sccomp.stackexchange.com лучшим местом для ответа на этот вопрос?
Ссылка была изменена на следующую: astro.gsu.edu/~gudehus/binary.html

Ответы (2)

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

Симулятор орбиты, см. здесь ; замените центральную звезду барицентром двойных.

Расчет орбиты по трем наблюдениям . Задача Кеплера в Википедии .

Если вы не найдете готового к использованию программного обеспечения, попробуйте решить его численно, моделируя орбиту двойных систем в соответствии с законом Кеплера, и варьируйте предположения о массе, расстоянии и эксцентриситете, пока они не будут соответствовать наблюдениям. Используйте методы оптимизации, например , алгоритмы подъема в гору или градиентные методы .

Или см., например: sourceforge.net/projects/phoebe
Требуется некоторое обобщение на незатменные двойные системы.

Хорошим программным обеспечением для расчета орбиты является «Binary» (Комбинированный пакет решений Binary Star) от команды CHARA:

http://www.astro.gsu.edu/~gudehus/binary.html

Ссылка мертва.