Simulation and modeling of Pseudo ranges measurements was applied widely in GPS system analysis and error's mitigation, and the simulated result of residues from some foregone modeling methods was not satisfied with the needs of special application. An approach of non-stationary time series analysis was given to reduce the simulation residues. The pseudo ranges data gathered from GPS were analyzed, and ARIMA model was modeled, and the model parameter was estimated using invert function. At last simulating with Matlab, and the results show using the statistical models that reproduce the pseudo range measurements with residues less than 3 meters, the satisfying results are obtained.