Through mobile communication data transmission, the precise Ephemeris and clock offset of the IGS website are obtained in real time, and the SBAS differential enhanced data is obtained through the receiver antenna; then The dual-frequency ionospheric non-differential observation value combination is used to eliminate the influence of the ionosphere, and use the extended Kalman filter (EKF) for ambiguity and coordinate smoothing; then use the mean clustering(FCM) algorithm to classify the positioning results, get the center point of each cluster, and add the center points of each cluster to find average on, the starting point of PPP convergence is obtained, and then the Hatch algorithm is used to quickly converge the PPF positioning result and the positioning result of cm is obtained. Compared with the previous PPP positioning method, this method shortens the convergence time by more than 90%, saves the time in the calibration process, and meets the engineering requirements of PPP rapid and high-precision positioning.