ГОСТ Р 52865— 2007
р</»=а<’»)(1-е2),п =0.1.2
....
где v =-<!>:
/= /ср+ л/;
и = т<р+*т-
За начальное приближение принимают а|0>= |ц(I2zf
I^dJm od 86400’
Wk= Г /где W— целая часть Wk\
f* = t,-tk + 86400 (N0- Na);
Xk=X+(Ji-C03)(ГгрИ^+ДПУ2);
Приближение заканчивается при выполнении условия |а<л-,)_а’п)| < 10-3 км.
Для этого обычно достаточно двух — трех итераций;
2) рассчитывают момент прохождения восходящего узла орбиты tkна витке, к которому относится момент f.
и долгота восходящего узла на этом витке равна А ■
cos/
(1-е2)2 ’
л =2тт/Г йр, Q =Ak+ S.S =S0+to3( ^ -10800);
где— коэффициент при второй зональной гармонике разложения геопотенциала Земли в ряд по сферичес
ким функциям, равный 1082.63 - 10*6;
afl — экваториальный радиус Земли, равный 6378,136 км;
SQ — истинное звездное время на гринвичскую полночь даты N0 , к которой относится время t’,
со3 — угловая скорость вращения Земли, равная 0,7292115 • 10~* с-1;
ц — константа гравитационного поля Земли, равная 398600,44 км3/ с2;
3) вычисляют константы интегрирования на момент fJJr
Sa’1"’ / a =2J{ae/ a)2(1-3/2 sin2!)(/cos 1) +J(aQ/ a)2sin i (1/2 ft sin 1- 1/2/ cos 1 +cos 2X+ 7/2/cos 31 * 7/2 ft sin 31):
8ft"”>= J{a01a f (1 - 3/2 sin2/) (/л т + sin 31 + 3/2 / sin 2L - 3/2 ft cos 21) - 1/43(3,, / a)2sin2/ (sin 1 - 7/3 sin 3/. +
+ 5 / sin 21 - 17/2 / sin 4L + 17/2 ft cos 4L + /) cos 21) + J{ae/ a)2 cos2/(/л x - 1/2/sin 2L);
8/<m>= J(ae/ a)2(1 - 3<2sin2!} ( -ftn x+ cos 1 + 3/2 / cos 2L +3/2 ft sin 21)- 1/4J(a0/ a)2sin2/ (-cos 1 - 7/3 cos 31 -
- 5/
j
sin 21 - 17/2 / cos 41 - 17/2 /»sin 41 + /cos 21) + J\ac/ a)2 cos2/(-ftn т + 1/2 ft sin 21);
8£2<m| = J(aeI a)2cos / (m + 7/2 / sin 1 - 7/2 ft cos 1 - 1/2sin 21 - 7/6 sin 31 + 7/6 ft cos 31);
5/,m| = 1/2 J(a0/ a)2sin / cos / (-1 cos 1 + ft sin 1 +cos 21 + 7/3 / cos 31 + 7/3 ft sin 31):
61(m| =2J(a„ / a)2(1 —3/2 sin /) (лг + 7/4 / sin 1 - 7/4 ft cos 1) + 3J(ac/ a)2sin / (-7/24 ft cos 1 -
- 7/24 / sin 1 - 49/72 ft cos 1 + 49/72 / sin 31 + 1/4 sin 21) + J(ac/ a)2 cos / (л г + 7/2/ sin 1 -
- 5/2 ft cos 1 - 1/2 sin 21 - 7/6 /sin 31 + 7/6 ft cos 31).
где 1 =M +со;
Af =E - esin £;
18