Abstract:A spline interpolation can obtain higherlevel smoothness using lowerlevel polynomial. In GPS position calculation, if we use suitable method to supply boundary conditions, the spline interpolation function which has better convergence and stability can be proposed. The quintic spline interpolation algorithm is introduced and standardized satellite orbit based on GPS broadcast ephemeris is calculated. The experimental results show this method has higher accuracy of approximation and better smoothness, can avoid some shortcomings compared with other interpolation methods.