Precise Orbit Determination of Global Navigation Satellite System of ...
Precise Orbit Determination of Global Navigation Satellite System of ...
Precise Orbit Determination of Global Navigation Satellite System of ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Chapter 8 Geostationary <strong>Orbit</strong> <strong>Determination</strong> And Prediction During <strong>Satellite</strong> Maneuvers<br />
R GM J �3<br />
2 a 3�112 � 35 J4<br />
a 5�332 3 4 ��<br />
s = � 3 ( ) � − sin i�<br />
− 5 ( ) � − sin i+ sin i�<br />
�2<br />
a r �3<br />
2 � 8 a r �35<br />
7 8 �<br />
�<br />
(8-6)<br />
�<br />
and the periodic variation perturbation function by<br />
R GM J � a �<br />
� J a �<br />
�<br />
p = � � i v+<br />
� − � i− v+ − i v+ � i<br />
� a r �<br />
� a r �<br />
�<br />
3 2 3 1 2 3 4 15 2 3<br />
5 2<br />
3 ( ) sin cos 2(<br />
ω) 4 ( ) ( sin )sin( ω) sin sin 3(<br />
ω)<br />
sin<br />
2 2<br />
8 2<br />
8<br />
35 J4<br />
a 5�2 3 1 2 1 4<br />
��<br />
− 5 ( ) �sin<br />
i( − sin i)cos 2(<br />
v+ ω) + sin icos 4(<br />
v+<br />
ω ) �<br />
8 a r � 7 2<br />
8<br />
�<br />
�<br />
(8-7)<br />
�<br />
Inserting Eq.(8-6) and Eq.(8-7) into Eq.(8-1), the variations <strong>of</strong> Ω and i due to perturbations Eq.(8-6) and Eq.(8-<br />
7) can be obtained as<br />
dΩ<br />
=<br />
dt 2<br />
na<br />
1 ∂R<br />
2<br />
1−<br />
e sin i ∂i<br />
=<br />
n<br />
2<br />
1−<br />
e<br />
�3<br />
J2<br />
a 3 35 J4<br />
a 5�632 ��<br />
� 2 ( ) ( −cos i)<br />
− 4 ( ) � cos i+ sin icos i�<br />
�2<br />
a r<br />
8 a r �7<br />
2 �<br />
�<br />
�<br />
(8-8)<br />
di<br />
=<br />
dt 2<br />
na<br />
cosi<br />
∂R<br />
2<br />
1 − e sini<br />
∂ω<br />
=<br />
ncosi �3<br />
J2<br />
a 3 J3<br />
a 4�15 2 3<br />
15 2<br />
�<br />
� ( ) { − sinisin 2(<br />
v+<br />
) } − ( ) �(<br />
sin i− )cos( v+ ) − sin icos 3(<br />
v+<br />
) �<br />
2 2<br />
ω 3<br />
ω ω<br />
1 − e �2<br />
a r<br />
a r � 8 2<br />
8<br />
�<br />
35 J4<br />
a 5�312 1 3<br />
��<br />
+ 4 ( ) �2sin<br />
i( − sin i)sin 2(<br />
v+ ω) + sin isin 4(<br />
v+<br />
ω ) �<br />
8 a r � 7 2<br />
2<br />
�<br />
�<br />
(8-9)<br />
�<br />
From two-body theory <strong>of</strong> satellite movement,<br />
2<br />
a( 1−<br />
e )<br />
�<br />
r = = a( 1−ecos<br />
E)<br />
�<br />
1+<br />
ecos v<br />
�<br />
2<br />
1−<br />
e sin E<br />
�<br />
sin v =<br />
�<br />
1−<br />
ecos E �<br />
�<br />
M = E− esin E = n( t−T0) �<br />
dM = ndt<br />
�<br />
�<br />
dM = ( 1−ecos<br />
E) dE �<br />
2 �<br />
a( 1−<br />
e )sinvdv<br />
sin EdE =<br />
�<br />
2<br />
( 1+<br />
ecos v)<br />
�<br />
�<br />
where<br />
n unperturbed two-body mean motion<br />
E eccentric anomaly<br />
T0 time <strong>of</strong> perifocal passage<br />
then the following equation can be obtained<br />
dM r<br />
= ( )<br />
dv a − e<br />
2 1<br />
2<br />
1<br />
97<br />
(8-10)<br />
(8-11)<br />
In order to evaluate Ω and i variations under the influence <strong>of</strong> nonspherical earth gravitation, Eq.(8-8) and Eq.(8-<br />
9) are integrated over a revolution <strong>of</strong> geostationary satellite movement. The results are<br />
�<br />
te<br />
te<br />
�<br />
1 J<br />
a J<br />
a<br />
∆Ω s =− �3<br />
2 1 3 35 4 �6<br />
3 2 � 1 5<br />
cos i ndt − � i+ i i�<br />
ndt�<br />
− e<br />
� a � ( ) cos sin cos<br />
r a �<br />
� � ( )<br />
2 2<br />
4<br />
1 2 2π<br />
8 7 2<br />
2π<br />
r �<br />
�<br />
ts<br />
ts<br />
�<br />
=−<br />
�<br />
2<br />
2<br />
1 3 J<br />
�<br />
� 2 1 a 3 35 4 �6<br />
3 2 � 1 5<br />
− � + �<br />
�<br />
2 2 i<br />
− � �<br />
4<br />
�<br />
� �<br />
1 e 2 a 2 r 8 7 2 2 �<br />
�<br />
0<br />
0 �<br />
dM<br />
π<br />
π<br />
J<br />
a<br />
cos ( ) cosi sin icos i ( ) dM<br />
π a<br />
π r<br />
(8-12)