The paper presents a new step-by-step procedure for constructing a navigation satellite motion model. At the first stage of the procedure, the parameters of the radiation pressure model are estimated using the maximum likelihood method. The statistic estimator based on the continuous-discrete adaptive unscented Kalman filter is proposed for the solar radiation model parameters estimation. Step-by-step scheme of filtering algorithm used for the software development are given. At the second stage, the parameters of the unaccounted perturbations model are estimated based on the results of residual differences measurements. The obtained results lead to significant improvement of prediction quality of the satellite trajectory.