14 / 23
Д.Ю. Виноградов, Е.А. Давыдов
14
Инженерный журнал: наука и инновации
# 6·2017
( 1)
( )
( )
1
3
j
j
j
a
a
T
a
+
µ
= + ∆
π
и осуществляется переход к п. 3. Номер итерации
1
j j
= +
.
8.
По формулам (14) определяются значения элементов
0 0
,
e
ω
.
9.
Формируются результаты первого этапа решения задачи:
н
др др 0 0 0 0
0
, , , ,
0,
.
=
Ω = ω
T T a e i
По формулам (21)
проводится расчет
элементов фазового вектора и выполняется настройка модели дви-
жения КА.
10.
Осуществляется уточнение элемента
0
a
. Для этого расчеты
согласно пп. 3–9 продолжаются до тех пор, пока не будет выполнено
условие (22). При этом значение
др( )
j
T
рассчитывается путем численно-
го прогнозирования движения КА в принятой модели движения.
11.
Проводится уточнение наклонения плоскости ССО. Для этого
при заданных значениях элементов
н
др др 0 0 0 0
, , , ,
T T a e i
=
ω
итераци-
онно решается уравнение (20) относительно
i
с использованием ме-
тода Ньютона первого порядка. На каждой итерации вычисляется
уточненное значение наклонения по формуле
2
( )
др
0( 1)
0( )
0( )
,
2 365, 2422sin
j
j
j
j
f
p T
i
i
i
+
µ
= +
πε ⋅
где
( )
j
f
рассчитывается согласно выражению (20). Итерационный
процесс завершается при выполнении условия
( 1)
,
j
f
f
+
≤ ε
где
6
10 .
f
−
ε =
12.
Осуществляется осреднение элементов
1 2
,
e e
(
0 0
,
e
ω
) с
использованием принятой модели движения КА по принципу, изло-
женному в работе [7], а именно проводится численное интегрирова-
ние движения КА на интервале времени, равном продолжительности
долгопериодического возмущения элементов
1 2
,
e e
. В таком случае
начальные средние значения элементов
0
1
e
и
0
2
e
будут определяться как
max min
0 1
1
1
;
2
e
e
e
+ =
max min
0 2
2
2
2
e
e
e
+ =
.
13.
Формируются результаты второго этапа решения задачи в виде
массива параметров программной орбиты:
{
н
др др 0 0 0 0
, , , ,
0,
T T a e i
=
Ω =
}
0
м.в сут c пр
.
, ,
,
,
ω Ω ∆ ∆
L L N
При этом значение величины
Ω
находит-