Performance Evaluation of Precise Point Positioning for BeiDou-3 B1c/B2a Signals in the Global Range

With the construction and development of the BeiDou navigation satellite system (BDS), the precise point positioning (PPP) performance of the BDS is worthy of research. In this study, observational data from 17 stations around the world across 20 days are used to comprehensively evaluate the PPP performance of BDS B1c/B2a signals. For greater understanding, the results are also compared with the Global Positioning System (GPS) and BDS PPP performance of different signals and system combinations. The evaluation found root mean square (RMS) values of the static PPP in the north (N), east (E), and upward (U) components, based on the B1c/B2a frequency of BDS-3, to be 6.9 mm, 4.7 mm, and 26.6 mm, respectively. Similar to the static positioning, the RMS values of kinematic PPP in the three directions of N, E, and U are 2.6 cm, 6.0 cm, and 8.5 cm, respectively. Besides this, the static PPP of BDS-3 (B1cB2a) and BDS-2 + BDS-3 (B1IB3I) have obvious system bias. Compared with static PPP, kinematic PPP is more sensitive to the number of satellites, and the coordinate accuracy in three dimensions can be increased by 27% with the combination of GPS (L1L2) and BDS. Compared with BDS-2+BDS-3 (B1IB3I), the convergence time of BDS-3 (B1CB2a) performs better in both static and kinematic modes. The antenna model does not show a significant difference in terms of the effect of the convergence speed, though the number of satellites observed has a certain influence on the convergence time.


Introduction
The BeiDou Navigation Satellite System (BDS) is a global satellite navigation system independently developed and operated by China. On 31 July 2020, The BDS-3 global satellite navigation system was officially announced to provide a global service [1,2]. In addition to the traditional B1I and B3I signals, BDS-3 also transmits several new signals with advanced signal structures, namely, B1C, B2a, B2b, and B2a + b. In terms of signal design, B1C signals are compatible with GPS L1C and Galileo E1 signals, while B2a signals are compatible with GPS L5C and Galileo E5a signals. The new B2b signal is used for smooth transition and compatibility with BDS-2, which can achieve better performance and is competitive with other GNSS signals [3]. Compared with BDS-2, the BDS-3 constellation is composed of three Geostationary Earth Orbit (GEO) satellites, three Inclined Geosynchronous Satellite Orbit (IGSO) satellites, and twenty-four Medium Earth Orbit (MEO) satellites, and these satellites of BDS-3 are manufactured by the China Academy of Space Technology (CAST) and Shanghai Engineering Center for Microsatellites (SECM) [4]. regarding the mathematical relationship between observations, observation errors, and parameters [36][37][38]. At present, the widely used PPP models include ionosphere-free PPP and ionosphere-weighted PPP [29,[39][40][41][42]. Each model has its own advantages. In this paper, a combination of ionospheric-free PPP models is used.
The basic observation equations of the BDS ionosphere-free carrier phase and pseudorange can be expressed as P s,T r,IF = ρ s,T + c · dt r − c · dt s,T + T s,T r + B T r,P IF − B s,T P IF + Re s,T r + Er s,T r + Pco s,T r,IF +Pcv s,T r,IF + Tid s,T r + ε s,T P IF where the subscript r and superscript s refer to the receiver and satellite, respectively; the subscript IF denotes ionosphere-free; the superscript T denotes the satellite system; P and L are the pseudo-range observations (m) and the phase observation (m), respectively; ρ is the geometric distance from the satellite to receiver (m); c is the speed of light in a vacuum (s); dt r and dt s are the receiver clock offset and satellite clock offset (m/s), respectively; T is the tropospheric delay on the signal propagation path (m); λ is the wavelength (m); Φ is the phase observations in cycles; N is the ambiguity of the carrier phase; B r and B s are the pseudo-range hardware delays at the receiver and satellite sides (m), respectively; W is the phase wind-up error in cycles; Re is the relativistic effect caused by satellite track eccentricity (m); Er is the error caused by the Earth's rotation; Pco and Pcv are the phase center offset (PCO) and phase center variation (PCV) of the receiver and satellite; Tid is the earth tide error; and ε are the unmodeled errors (multipath, noise). ρ = (x r − x s ) 2 + (y r − y s ) 2 + (z r − z s ) 2 (3) where [x r , y r , z r ] and [x s , y s , z s ] are the coordinates of the receiver and satellite respectively. Then, linearize ρ at its approximate coordinate [x 0 , y 0 , z 0 ], which can be written as where ρ 0 is the geometric distance from the satellite to the receiver's approximate coordinates.
Observations are used to establish the PPP model, but the accuracy of the tropospheric model does not meet the requirements of PPP. Hence, the residual part of the zenith tropospheric wet delay is normally estimated together with station coordinates and carrier phase ambiguity. However, other error corrections, such as PCO, PCV, phase wind-up, relativistic effect, etc., are corrected by corresponding methods. The unknown parameters are the receiver coordinates and clock errors, zenith tropospheric wet delay, and phase ambiguity. Therefore, the observation equation of PPP is formed as: where M wet is the tropospheric wet delay mapping function; dZTD ω is the zenith tropospheric wet delay; B is the ambiguity parameter. Formula (5) is abbreviated as y = Gx.

Kalman Filtering
Kalman filtering is the most commonly used algorithm in GNSS positioning. Kalman filtering is based on the state estimation of the previous epoch combined with observations of the current epoch to recursive new state parameters with the state transition matrix. It is mainly divided into two steps: state prediction and parameter update. The prediction model can be expressed as where k is the epoch; x is the estimated parameters; Φ is the state transition matrix from the previous epoch to the current epoch; Q is the process noise matrix; P is the covariance matrix of the parameters; * is the estimated value; * − is a priori value. In Kalman filtering, the setting of kinematic and static modes is different. Differences are mainly manifested in the setting of state transition matrix Φ and process noise matrix Q [43,44].
For static PPP, the coordinate parameters are unchanged, the clock error parameter is a random noise model, and the tropospheric parameters are generally a random walk model. When there is no cycle slip, the ambiguity parameter remains unchanged, so the included coordinates and receiver clock error and the state transition matrix Φ and process noise matrix Q of the tropospheric wet delay and ambiguity parameters can be written as: where σ δt = 1 ms and σ trop = 20 mm/ √ hour; when the phase cycle slips, the coefficient in the state transition matrix is 0.
For kinematic PPP, if its moving speed is unknown and the coordinate parameters are random noise models, its state transition matrix Φ and process noise matrix Q can be defined as: where σ dx , σ dy , σ dz is the coordinate noise, and the kinematic PPP is set to 100 m. The parameter update step can be expressed as where y is the observation minus correction; G is the observation design matrix; K(k) the Kalman gain, which can be expressed as

Data source and Configuration Analysis
To verify the performance of BDS PPP, observation data of 17 global stations from 1 December 2020 to 20 December 2020 were used. BDS B1CB2a, B1IB3I, and GPS signals can be received by the test stations. Taking a day at one station as a sample, there are in total 340 data samples. The observation data interval is 30 s. Since BDS-3 was built to serve the world, to ensure the objectivity of the evaluation, the selection of the stations around the world should be as uniform as possible. The detailed distribution of the stations is shown in Figure 1.   The precise ephemeris and orbit provided by the IGS analysis centers include final, rapid, ultra-rapid, and real-time (RTS) products. The accuracy of the products provided by each analysis center is slightly different [45]. The precise orbit and clock products used in this research are the GBM products provided by GFZ (Geo Forschungs Zentrum, Potsdam, Germany). The coordinate reference of the stations is the SNX file provided by IGS. It should be noted that the PCO and PCV corrections of some antennas of the receiver only provide corrections for GPS and GLONASS [27]. Therefore, the PCO and PCV corrections of GPS are used to correct for BDS.
The GNSS multi-function precision positioning software Net_Diff is used for PPP in this research (https://github.com/YizeZhang/Net_Diff) (accessed on 20 August 2021). The specific settings for PPP are shown in Table 1.

Experimental Process
In this study, the abnormal value of a single sample has been excluded, which refers to a single sample with a small number of satellites participating in the PPP calculation. The excluded samples for static PPP and kinematic PPP are shown in Table 2, where G, C12, and C13, respectively, represent the positioning modes of GPS(L1L2), BDS (B1cB2a), and BDS (B1IB3I). The number (DOY) + station name represents the day of year and the station. For example, 337PIE1 represents the sample of PIE1 on the 337th day in 2020, and "all" represents all days from 336-355.
For the same observation data, regardless of whether static PPP or kinematic PPP, the number of satellites involved in positioning is the same. This article counts the average number of satellites participating in BDS (B1CB2a), BDS (B1IB31), and GPS (L1L2) PPP for 17 stations, as shown in Figure 2. The PPP availability of this station can be better reflected by the average number of satellites. The stations with a small average number of satellites shown in Figure 2 can basically correspond to the stations excluded in Table 2. For the overall performance of static PPP, due to the limited number of satellites, there are large errors, which have a great impact on positioning statistics and lack statistical objectivity. Figure 3 shows the static PPP error at station PIE1 for BDS-2+3 (B1IB3I), and the number of satellites involved in positioning on DOY 343 in 2020. In Figure 3, there are less than five satellites participating in static PPP for a long time, and the position information cannot be calculated in static PPP. For the overall performance of static PPP, due to the limited number of satellites, there are large errors, which have a great impact on positioning statistics and lack statistical objectivity. Figure 3 shows the static PPP error at station PIE1 for BDS-2+3 (B1IB3I), and the number of satellites involved in positioning on DOY 343 in 2020. In Figure 3, there are less than five satellites participating in static PPP for a long time, and the position information cannot be calculated in static PPP.  For kinematic PPP, the number of satellites has a very large impact on positioning. When the number of satellites is less than five, abnormal jumps occur, as shown in Figure  4. For kinematic PPP, the number of satellites has a very large impact on positioning. When the number of satellites is less than five, abnormal jumps occur, as shown in Figure 4.
Therefore, in the subsequent data analysis, the above cases are excluded.  For kinematic PPP, the number of satellites has a very large impact on positioning. When the number of satellites is less than five, abnormal jumps occur, as shown in Figure  4. Therefore, in the subsequent data analysis, the above cases are excluded.

Static PPP
The precise orbit and clock offset provided by GBM are used to perform static PPP of BDS (B1CB2a). Figure 5 demonstrates a static PPP error series for 20 days at 17 stations utilizing GBM's precise orbit/clock corrections. All results for the 20 days at 17 stations are plotted in the same figure. The premise is to eliminate the abnormal values in Table 2, where S represents the results for static PPP. It can be seen that all PPP solutions can reach convergence but with significantly different convergence times. In the meantime, the positioning errors in the N and E directions are obviously smaller than those of the U component.

Static PPP
The precise orbit and clock offset provided by GBM are used to perform static PPP of BDS (B1CB2a). Figure 5 demonstrates a static PPP error series for 20 days at 17 stations utilizing GBM's precise orbit/clock corrections. All results for the 20 days at 17 stations are plotted in the same figure. The premise is to eliminate the abnormal values in Table 2, where S represents the results for static PPP. It can be seen that all PPP solutions can reach convergence but with significantly different convergence times. In the meantime, the positioning errors in the N and E directions are obviously smaller than those of the U component. When performing static PPP, the last epoch of a one-day solution is taken as the re sult of the static positioning of that day (mark with location point; no explanation is given for the subsequent static PPP position error statistics, and the same processing method i used). Figure 6 shows the distribution of the north (N), east (E), and upward (U) posi tioning errors of all stations on all days (BDS(B1cB2a)). It can be seen that the accuracie of BDS (B1cB2a) static PPP in the three directions of N, E, and U are 6.9 mm, 4.7 mm, and When performing static PPP, the last epoch of a one-day solution is taken as the result of the static positioning of that day (mark with location point; no explanation is given for the subsequent static PPP position error statistics, and the same processing method is used). Figure 6 shows the distribution of the north (N), east (E), and upward (U) positioning errors of all stations on all days (BDS(B1cB2a)). It can be seen that the accuracies of BDS (B1cB2a) static PPP in the three directions of N, E, and U are 6.9 mm, 4.7 mm, and 26.6 mm, respectively. Compared with the horizontal error, the distribution of vertical error is more scattered.  It can be seen from Table 3 and Figure 7 that the upward error of BDS-3 (B1cB2a) is larger than that of GPS and BDS-2+3 (B1IB3I). Compared with BDS (B1IB3I) and BDS-3 (B1cB2a), the horizontal error of GPS is 16.8% and 26.5% better, and the vertical component is 46.3% and 69.5% better, respectively. On the one hand, the precise orbit determination of BDS satellite PCO on-orbit estimation, precise attitude model, solar radiation pressure model, and other aspects of precise orbit determination strategies are not perfect yet, and using the correction values of GPS to correct the PCV and PCO of BDS will lead to significantly lower accuracy than that of GPS, especially in the vertical component. On the other hand, the number of satellites with B1C and B2a signals involved in static PPP is less than that of GPS and BDS-2+3 (B1IB3I). BDS-3 (B1cB2a) is significantly different from BDS-2+3 (B1IB3I) in terms of mean error (M-RMS). The mean error of BDS-2+3 (B1IB3I) is within 2mm and there is no large system error, while BDS-2+3 (B1IB3I) reached -23.2mm in the vertical direction, indicating a difference between PCO and PCV at different frequencies of BDS. Compared with the BDS-only system, the positioning accuracy of combined GPS and BDS is significantly improved. However, due to the influence of the BDS PCO correction issue, the performance of the vertical error is not ideal, and even slightly lower than the accuracy of GPS-only static PPP. Compared with BDS-3 (B1cB2a), GPS+BDS-3 (B1cB2a) increased by 27% and 52% in the horizonal and vertical directions, respectively.  Table 3.

Kinematic PPP
The kinematic solutions of all stations are shown in Figure 8, where K represents the results for kinematic PPP. Unlike the static solutions in Figure 5, the Y-axis in Figure 8 is expanded to ±2 m for a clearer view. Similar to static solutions, the kinematic positioning  It can be seen from Table 3 and Figure 7 that the upward error of BDS-3 (B1cB2a) is larger than that of GPS and BDS-2+3 (B1IB3I). Compared with BDS (B1IB3I) and BDS-3 (B1cB2a), the horizontal error of GPS is 16.8% and 26.5% better, and the vertical component is 46.3% and 69.5% better, respectively. On the one hand, the precise orbit determination of BDS satellite PCO on-orbit estimation, precise attitude model, solar radiation pressure model, and other aspects of precise orbit determination strategies are not perfect yet, and using the correction values of GPS to correct the PCV and PCO of BDS will lead to significantly lower accuracy than that of GPS, especially in the vertical component. On the other hand, the number of satellites with B1C and B2a signals involved in static PPP is less than that of GPS and BDS-2+3 (B1IB3I). BDS-3 (B1cB2a) is significantly different from BDS-2+3 (B1IB3I) in terms of mean error (M-RMS). The mean error of BDS-2+3 (B1IB3I) is within 2mm and there is no large system error, while BDS-2+3 (B1IB3I) reached -23.2mm in the vertical direction, indicating a difference between PCO and PCV at different frequencies of BDS. Compared with the BDS-only system, the positioning accuracy of combined GPS and BDS is significantly improved. However, due to the influence of the BDS PCO correction issue, the performance of the vertical error is not ideal, and even slightly lower than the accuracy of GPS-only static PPP. Compared with BDS-3 (B1cB2a), GPS+BDS-3 (B1cB2a) increased by 27% and 52% in the horizonal and vertical directions, respectively.

Kinematic PPP
The kinematic solutions of all stations are shown in Figure 8, where K represents the results for kinematic PPP. Unlike the static solutions in Figure 5, the Y-axis in Figure 8 is expanded to ±2 m for a clearer view. Similar to static solutions, the kinematic positioning errors in the vertical component are larger than those of the horizontal component. However, the kinematic positioning errors are obviously larger when compared with the static solutions.
Similar to the evaluation method of static PPP, the positioning error within one hour after positioning is taken as the result of convergence, and the RMS of the PPP error after convergence is counted. The kinematic positioning results of BDS-3 (B1cB2a) are shown in Figure 9.
The kinematic solutions of all stations are shown in Figure 8, where K represents the results for kinematic PPP. Unlike the static solutions in Figure 5, the Y-axis in Figure 8 is expanded to ±2 m for a clearer view. Similar to static solutions, the kinematic positioning errors in the vertical component are larger than those of the horizontal component. However, the kinematic positioning errors are obviously larger when compared with the static solutions.   It can be seen from Figure 9 that the RMS values of the BDS-3 (B1cB2a) kinematic PPP results in the three directions of N, E, and U are 2.6, 6.0, and 8.5cm, respectively. When the horizontal error is converged, the positioning accuracy reaches 6.5cm, which proves that kinematic PPP can reach centimeter-level positioning accuracy.
To better compare the results of BDS-3 (B1cB2a) kinematic PPP, the RMS of kinematic PPP are calculated for BDS-3 (B1cB2a), GPS, BDS-2+3 (B1IB3I), GPS + BDS-3 (B1cB2a), and GPS + BDS-2+3 (B1IB3I) in the N, E, U, and horizontal directions, and in 3D, as shown in Figure 10. It can be seen from Figure 9 that the RMS values of the BDS-3 (B1cB2a) kinematic PPP results in the three directions of N, E, and U are 2.6, 6.0, and 8.5 cm, respectively. When the horizontal error is converged, the positioning accuracy reaches 6.5 cm, which proves that kinematic PPP can reach centimeter-level positioning accuracy.
PPP results in the three directions of N, E, and U are 2.6, 6.0, and 8.5cm, respectively. When the horizontal error is converged, the positioning accuracy reaches 6.5cm, which proves that kinematic PPP can reach centimeter-level positioning accuracy.
To better compare the results of BDS-3 (B1cB2a) kinematic PPP, the RMS of kinematic PPP are calculated for BDS-3 (B1cB2a), GPS, BDS-2+3 (B1IB3I), GPS + BDS-3 (B1cB2a), and GPS + BDS-2+3 (B1IB3I) in the N, E, U, and horizontal directions, and in 3D, as shown in Figure 10.   Figure 10 shows that compared with BDS-2+3 (B1IB3I), kinematic PPP of BDS-3 (B1cB2a) reduces the positioning errors of N, E, U, and 3-D error by 29%, 36%, 12%, and 24%, respectively. It can be seen that the observation value of B1cB2a is more accurate. The kinematic PPP accuracy of BDS-3 (B1cB2a) is slightly lower than that of GPS. The reason is that the satellite orbits and clock products provided by GBM do not include all BDS satellites, resulting in a relatively small number of satellites that can be used in positioning. On the other hand, the accuracy of the BDS orbit and clock is worse than that of GPS. Different from the impact of static PPP, the number of satellites has a larger impact on kinematic PPP. For the dual-system combination, due to the increase in the number of satellites, the 3-D error is increased by 59% and 70%, respectively, relative to BDS-3 (B1cB2a) and BDS-2+3 (B1IB3I). The effect of multi-GNSS on improving the positioning accuracy is more obvious in kinematic positioning. Considering that the receiver coordinates are estimated as white noise, the receiver coordinate parameters cannot be constrained by the explicit kinematics model as in the static processing mode. Therefore, increasing redundant observations can increase the reliability of kinematic solutions [9].

Convergency
Compared with RTK, PPP requires a longer convergence time to achieve centimeterlevel positioning accuracy. In this section, factors of kinematic and static PPP convergence time are discussed from the previous hourly convergence performance and receiver antenna type.
vergence time are discussed from the previous hourly convergence performance and receiver antenna type.
(a) Static PPP convergence comparison (b) Kinematic PPP convergence comparison  Figure 11 shows that the convergence efficiency of the static PPP of a single system is significantly higher than that of a kinematic one. In addition, the differences between different systems are also obvious. Among them, whether it is static or kinematic PPP, the convergence effect of BDS (B1IB3I) is worse those that of GPS and BDS (B1CB2a), mainly due to the number of satellites observed by individual stations is small, which leads to poor overall convergence, as shown in Figures 2 and 12. When the number of observation satellites is the same, because the observation noise of BDS (B1IB3I) is relatively large, the error of BDS (B1IB3I) PPP is greater than that of BDS (B1CB2a) PPP. Compared with the single system, the dual system has significantly improved the con-  Figure 11 shows that the convergence efficiency of the static PPP of a single system is significantly higher than that of a kinematic one. In addition, the differences between different systems are also obvious. Among them, whether it is static or kinematic PPP, the convergence effect of BDS (B1IB3I) is worse those that of GPS and BDS (B1CB2a), mainly due to the number of satellites observed by individual stations is small, which leads to poor overall convergence, as shown in Figures 2 and 12. When the number of observation satellites is the same, because the observation noise of BDS (B1IB3I) is relatively large, the error of BDS (B1IB3I) PPP is greater than that of BDS (B1CB2a) PPP. Compared with the single system, the dual system has significantly improved the convergence efficiency.  To further analyze the convergence of BDS-3 (B1cB2a), the receiver types were classified and counted, and the positioning errors of each station at different epochs were observed. Table 4 shows the receiver model and antenna model at each station.  To further analyze the convergence of BDS-3 (B1cB2a), the receiver types were classified and counted, and the positioning errors of each station at different epochs were observed. Table 4 shows the receiver model and antenna model at each station. The 20-day data are used for static PPP, and the positioning error of each station during the first hour is analyzed according to the receiver type, as shown in Figure 13.  Figure 13 shows that there is no obvious difference between the receiver type and the convergence time. The main reason is the technological advancement of various manufacturers, and the antenna of the high-end receiver has little effect on BDS (B1CB2a) PPP. The PPP convergence speed is closely related to the number of satellites. Figure 14 shows the percentage with less than five satellites participating in static PPP at the PIE1 station in a day (20 day in total). Compared with other stations, the static PPP of PIE1 station has a poorer convergence performance, as shown in Figure 13. The percentage of PIE1 stations with less than five satellites participating in the 20-day PPP solutions is more than 40%, and the static PPP does not converge to less than 0.1 m within 60 min.  Figure 13 shows that there is no obvious difference between the receiver type and the convergence time. The main reason is the technological advancement of various manufacturers, and the antenna of the high-end receiver has little effect on BDS (B1CB2a) PPP. The PPP convergence speed is closely related to the number of satellites. Figure 14 shows the percentage with less than five satellites participating in static PPP at the PIE1 station in a day (20 day in total). Compared with other stations, the static PPP of PIE1 station has a poorer convergence performance, as shown in Figure 13. The percentage of PIE1 stations with less than five satellites participating in the 20-day PPP solutions is more than 40%, and the static PPP does not converge to less than 0.1 m within 60 min.
PPP. The PPP convergence speed is closely related to the number of satellites. Figure 14 shows the percentage with less than five satellites participating in static PPP at the PIE1 station in a day (20 day in total). Compared with other stations, the static PPP of PIE1 station has a poorer convergence performance, as shown in Figure 13. The percentage of PIE1 stations with less than five satellites participating in the 20-day PPP solutions is more than 40%, and the static PPP does not converge to less than 0.1 m within 60 min.

Conclusions
In this study, through the GNSS multifunctional precision positioning software Net_Diff, 20-day data for 17 stations in global were assessed to analyze the BDS-3 (B1cB2a) static and kinematic PPP performance. The BDS-2+3 (B1IB3I), GPS, GPS+BDS-3 (B1cB2a), and GPS + BDS-2+3 (B1IB3I) positioning error, static and kinematic convergence, and convergence efficiency between different receiver models were compared. The conclusions are summarized as follows.
(1) In terms of PPP error, the results show that the RMS accuracies of BDS (B1cB2a) static PPP are 6.9 mm, 5.2 mm, and 27.2 mm, and of kinematic PPP are 2.6 cm, 6.0 cm, and 8.5 cm in the North, East and Upward directions, respectively. Whether it is static PPP or kinematic PPP, the positioning accuracy of BDS (B1CB2a) and BDS (B1IB3I) PPP are worse than that of GPS. The main reason is that the accuracy of BeiDou's precise orbit and clock difference provided by GBM is worse than that of GPS, and the BeiDou PCO correction model uses GPS correction parameters. In static PPP, BDS (B1cB2a) and BDS (B1IB3I) have obvious inter-system bias, which is mainly reflected in the upward. Kinematic PPP has a greater impact on the number of satellites. (2) In terms of Convergency, for static PPP or kinematic PPP, BDS + GPS has better convergence than GPS only or BDS only. The main reason is the increase in the number of satellites involved in positioning. For BDS (BICB2a) PPP, there is no obvious difference in convergence between different antenna models (high-end receiver), though the number of satellites observed has a great impact on convergence. Because the observation noise of BDS (B1IB3I) is larger than that of BDS (B1CB2a), when the number of satellites is the same, the convergence of BDS (B1IB3I) is worse than that of BDS (B1CB2a). When the number of satellites participating in static PPP is less than five satellites and there are too many epochs of the total epochs in a day, the static PPP convergence does not reach 0.1 m within 60 min.
To sum up, BDS (B1CB2a) PPP can meet the PPP needs of BeiDou users in the global range, but compared with GPS, the elevation direction BDS (B1CB2a) PPP performs poorly. With the advancement of algorithms, these gaps will gradually narrow. BDS (B1CB2a) is equivalent to GPS in terms of convergence.