임의의 순간 인공위성의 위치와 속도를 정밀하게 계산하기 위해서는 섭동력을 일으키는 우주공간의 환경을 정확하게 이해하고 분석하여 정량화함으로써 섭동력에 대한 수리적인 모형을 만들어야 한다. 이들 우주환경모델에 의해서 인공위성이 받는 총가속도는 2계 미분방정식으로 표현되며, 이 방정식을 두 번 적분함으로써 원하는 시각에서의 인공위성의 위치와 속도를 얻는 코웰방법을 사용하여 궤도예측 알고리즘을 완성하였다. 정지위성의 궤도에 미치는 주요한 섭동력으로는 지구의 비대칭 중력 포텐셜에 의한 섭동력, 태양과 달의 중력에 의한 섭동력, 태양의 복사압에 의한 섭동력들이 있는데, 그것들의 정밀성을 최대한 높이기 위해 spherical harmonic 계수들을 40
To calculate the position and velocity of the artificial satellite precisely, one has to build a mathematical model concerning the perturbations by understanding and analysing the space environment correctly and then quantifying. Due to these space environment model, the total acceleration of the artificial satellite can be expressed as the 2nd order differential equation and we build an orbit propagation algorithm by integrating twice this equation by using the Cowell's method which gives the position and velocity of the artificial satellite at any given time. Perturbations important for the orbits of geostationary spacecraft are the Earth's gravitational potential, the gravitational influences of the sun and moon, and the solar radiation pressure. For precise orbit propagation in Cowell' method, 40 x 40 spherical harmonic coefficients can be applied and the JPL DE403 ephemeris files were used to generate the range from earth to sun and moon and 8th order Runge-Kutta single step method with variable step-size control is used to integrate the the orbit propagation equations.
