Вычисление координат спутника
Для эллиптической орбиты используется следующий алгоритм вычислений координат спутника на момент времени (эпоху) t:
1. Вычисление средней аномалии М по формулам
(18)
где wc – угловая скорость обращения спутника на орбите, Т – период обращения, τ – время прохождения через перигей, fMЗ – геоцентрическая гравитационная постоянная Земли.
2. Вычисление эксцентрической аномалии Е последовательными итерациями по формуле
(19)
3. Определение радиус-вектора R спутника по формуле:
(20)
4. Вычисление истинной аномалии υ и аргумента широты u по формулам
(21)
(22)
5. Определение прямоугольных геоцентрических координат спутника по формуле
(23)
Алгоритм и формулы упрощаются для круговых орбит, эксцентриситет которых e равен 0.