Next Article in Journal
Experimental Study of Multispectral Characteristics of an Unmanned Aerial Vehicle at Different Observation Angles
Previous Article in Journal
Amyloid Beta Detection by Faradaic Electrochemical Impedance Spectroscopy Using Interdigitated Microelectrodes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Real-Time Cycle Slip Detection and Repair Method under High Ionospheric Activity for a Triple-Frequency GPS/BDS Receiver

1
School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China
2
Collaborative Innovation Center for Geospatial Technology, Wuhan 430079, China
3
Key Laboratory of Geophysical Geodesy, National Administration of Surveying, Mapping and Geoinformation, Wuhan 430079, China
4
Faculty of Information Engineering, China University of Geosciences, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Submission received: 28 November 2017 / Revised: 21 January 2018 / Accepted: 29 January 2018 / Published: 1 February 2018
(This article belongs to the Section Remote Sensors)

Abstract

:
Cycle slip detection and repair is a prerequisite for high-precision global navigation satellite system (GNSS)-based positioning. With the modernization and development of GNSS systems, more satellites are available to transmit triple-frequency signals, which allows the introduction of additional linear combinations and provides new opportunities for cycle slip detection and repair. In this paper, we present a new real-time cycle slip detection and repair method under high ionospheric activity for undifferenced Global Positioning System (GPS)/BeiDou Navigation Satellite System (BDS) triple-frequency observations collected with a single receiver. First, three optimal linearly independent geometry-free pseudorange minus phase combinations are selected to correctly and uniquely determine the cycle slips on the original triple-frequency carrier phase observations. Then, a second-order time-difference algorithm is employed for the pseudorange minus phase combinations to mitigate the impact of between-epoch ionospheric residuals on cycle slip detection, which is especially beneficial under high ionospheric activity. The performance of the approach is verified with static GPS/BDS triple-frequency observations that are collected with a 30 s sampling interval under active ionospheric conditions, and observations are manually inserted with simulated cycle slips. The results show that the method can correctly detect and repair cycle slips at a resolution as small as 1 cycle. Moreover, kinematic data collected from car-driven and airborne experiments are also processed to verify the performance of the method. The experimental results also demonstrate that the method is effective in processing kinematic data.

1. Introduction

The availability of high-quality carrier phase observations is a prerequisite for high-precision global navigation satellite system (GNSS)-based positioning. However, under conditions that cause signal interruption, a low signal-to-noise ratio or high receiver dynamics, GNSS receivers may lose track of navigation signals, and cycle slips occur, which result in discontinuities in the original carrier phase observations. Without careful consideration of such cycle slips, the precision and reliability of position and velocity estimates are seriously affected.
Over the past three decades, a number of cycle slip detection and repair methods have been developed and can be generally divided into three types according to the number of available signals. Polynomial fitting [1] and high-order between-epoch phase differentiation [2] can be implemented based on only single-frequency observations. However, they cannot detect small cycle slips due to the presence of clock error. Based on multi-constellation signals, Qian et al. [3] further proposed a cycle slip detection and repair method for a single-frequency Global Positioning System (GPS)/BeiDou Navigation Satellite System (BDS) receiver. But it is not applicable for single-system receiver because both GPS and BDS signals are required. Outlier detection methods, such as statistical hypothesis test (e.g., χ 2 test) [4,5], robust estimation [6], the Bayesian method [7,8], and generalized likelihood ratio (CLR) test [9], can be utilized to detect and correct cycle slips as well. There are also many cycle slip detection and repair methods based on dual-frequency observations. Methods to detect and repair cycle slips using dual-frequency double-differenced or triple-differenced observations have been presented by Gao and Li [10], Bisnath and Langley [11], and Kim and Langley [12,13]. However, these methods require at least two receivers to realize the double-differenced or triple-differenced operation. To meet the requirement of high-precision positioning, such as Precise Point Positioning (PPP), it is more convenient and effective to implement cycle slip detection and repair based on undifferenced observations. The TurboEdit algorithm developed by Blewitt [14] is recognized as the first approach for undifferenced cycle slip detection and repair that uses dual-frequency observations. Both Melbourne–Wübbena (MW) linear combination [15,16,17] and ionospheric residual combination are used in TurboEdit. Although the TurboEdit algorithm is widely used, it has two disadvantages: First, the ionospheric residual combination is sensitive to variations in ionospheric delay residuals, which degrade the efficiency of the method. Furthermore, the wavelength of wide-lane (WL) combination is still relatively small, which also degrades the performance of the method when there is substantial pseudorange noise. de Lacy et al. [18] introduced the Bayesian criterion and polynomial fitting to detect cycle slips for dual-frequency observations with high sampling interval (e.g., 1 s or 5 s). At low sampling rate or under high ionospheric activity, their method is not suitable due to the large residuals produced by the ionospheric delay terms. Based on the TurboEdit algorithm, Liu [19] and Cai et al. [20] specially study cycle slip detection and repair under high ionospheric activity. However, their methods may not work well in situations with significant measurement noise. Based on a time-differenced model, the dual-frequency cycle slip detection and repair method has also been investigated by Zhang and Li [21], and Banville and Langley [22]. The external information such as precise orbit and clock corrections are required in these methods. There are also some cycle slip detection and repair methods that integrate the GNSS and other sensors. Aided by inertial navigation system (INS), cycle slip detection can be implemented with a single-frequency receiver [23,24] and dual-frequency receiver [25]. However, these methods are not flexible in application scenario due to the cost and complexity of GPS/INS.
With the modernization and development of GNSS systems, triple-frequency signals are now available, which allows the introduction of additional linear combinations with longer wavelengths, weaker ionospheric delays, and smaller noise [26,27]. The availability of triple-frequency observations provides new opportunities for cycle slip detection and repair. Dai et al. [28] proposed an instantaneous triple-frequency cycle slip detection and repair method that applies two geometry-free phase combinations with the least-squares ambiguity decorrelation adjustment (LAMBDA) [29] algorithm to search for cycle slip candidates. Since only two geometry-free phase combinations are used in the method, some undetectable cycle slips remain. Based on triple-frequency observations combination, Huang et al. [30], Li et al. [31], and de Lacy et al. [32] developed real-time cycle slip detection and repair methods. These methods assume that ionospheric delay varies smoothly over time. The performance of these methods might be reduced under high ionospheric activity. Zhao et al. [33], Yao et al. [34], and Li et al. [35] further proposed to compensate for the ionospheric bias by using ionospheric prediction. However, several epochs of continuous phase data that lack cycle slips are required to satisfy the criteria for the phase connection, e.g., 5 or 10 historic epochs are required in method presented by Li et al. [35]. Huang et al. [36] used three geometry-free and ionosphere-free (GIF) pseudorange minus phase combinations to eliminate the impact of ionospheric delay. Although GIF combinations are not affected by ionospheric delay, the cycle slip detection success rate of GIF combinations is decreased when there is significant pseudorange noise. On the other hand, Zhang [37] extended the geometry-based approach by integrating time-differenced pseudorange and carrier phase observations to estimate the integer number of triple-frequency cycle slips.
In this paper, we present a new real-time cycle slip detection and repair method under high ionospheric activity for undifferenced GPS/BDS triple-frequency observations collected with a single receiver. In this new approach, we jointly use three linearly independent geometry-free pseudorange minus phase combinations to efficiently detect and correctly repair the cycle slips on original triple-frequency carrier phase observations. Furthermore, a second-order time-difference algorithm is employed to mitigate the impact of between-epoch ionospheric variations, which improves the rounding success rate of estimated float cycle slips under high ionospheric activity. As observations of only three consecutive epochs (i.e., the current epoch and two previous epochs) are needed, cycle slips occurred at the current epoch can be efficiently detected and repaired in real time. Finally, we verify the performance of the approach by conducting static experiments during a magnetic storm and kinematic (car-driven and airborne) experiments. It is demonstrated that the presented method can detect and correctly repair cycle slips that are as small as 1 cycle in real time under high ionospheric activity. Compared with the existing double-frequency cycle slip detection and repair methods under high ionospheric activity [19,20], the impact of pseudorange noise will be reduced as the wavelength of triple-frequency combinations used in our proposed method is longer than that of dual-frequency WL combination. Compared with the triple-frequency cycle slip detection and repair methods utilizing ionospheric prediction [33,34,35], less consecutive epochs are required in our proposed method. Additionally, no external information is required in the proposed method with respect to the geometry-based time-differential positioning method, which requires precise orbit and clock correction [21,22].
The paper is organized as follows. In Section 2, we present the cycle slip detection and repair method based on geometry-free triple-frequency pseduorange minus phase combinations and second-order time-difference algorithm. In Section 3, we select optimal cycle slip detection and repair combinations for GPS and BDS, respectively. In Section 4, the proposed method is verified with real receiver data and simulated cycle slip, and the conclusions follow in Section 5.

2. Method of Triple-Frequency Cycle Slip Detection and Repair

The undifferenced pseudorange and carrier phase observation equations can be expressed as [38]:
P i = ρ + q i I i + ε P i
λ i φ i = ρ q i I 1 + λ i N i + λ i ε φ i
where P i and φ i are pseudorange observations (in meters) and carrier phase observations (in cycles), respectively, on frequency f i . ρ is the geometric distance from satellite to receiver, and includes the satellite and receiver clock errors, the tropospheric error, and the receiver and satellite hardware delays (in meters). I 1 is the first-order ionospheric delay at frequency f 1 (in meters) and it is multiplied by q i = f 1 2 / f i 2 to obtain the corresponding delay at frequency f i . λ i is the wavelength (in meters), N i is the integer ambiguity, and ε P i and ε φ i are pseudorange noise (in meters) and carrier phase observation noise (in cycles), respectively.
The geometry-based, linearly combined triple-frequency carrier phase observation equation is defined as [39]:
λ i j k φ i j k = i f 1 λ 1 φ 1 + j f 2 λ 2 φ 2 + k f 3 λ 3 φ 3 i f 1 + j f 2 + k f 3 = ρ β i j k I 1 + λ i j k N i j k + λ i j k ε φ i j k
where the combination coefficients i , j , k are integers. The linearly combined wavelength, integer ambiguity, and the ionospheric scale factor (ISF) are [39]:
λ i j k = c i f 1 + j f 2 + k f 3
N i j k = i N 1 + j N 2 + k N 3
β i j k = f 1 2 ( i / f 1 + j / f 2 + k / f 3 ) i f 1 + j f 2 + k f 3
Similarly, the geometry-based, linearly combined triple-frequency pseudorange observation equation is defined as:
P l m n = l P 1 + m P 2 + n P 3 = ρ + β l m n I 1 + ε P l m n
where the combination coefficients l , m , n are real numbers with a sum of 1 (i.e., l + m + n = 1 ). β l m n = l + m f 1 2 / f 2 2 + n f 1 2 / f 3 2 is the ISF of P l m n .
According to Equations (3) and (4), the geometry-free pseduorange minus phase combination can be expressed as:
φ i j k P l m n λ i j k = K i j k , l m n I 1 + N i j k + ε φ i j k ε P l m n / λ i j k
where K i j k , l m n = ( β i j k + β l m n ) / λ i j k . The geometry-free pseduorange minus phase combination can cancel out the geometric distance between satellite and receiver, the tropospheric delay, and satellite and receiver clock errors and hardware delays. Consequently, it is extremely suitable for cycle slip detection and repair for undifferenced observations under kinematic conditions.
The single-differenced observation equation for geometry-free pseudorange minus phase combinations between two consecutive epochs is defined as:
φ ˙ i j k P ˙ l m n λ i j k = K i j k , l m n I ˙ 1 + N ˙ i j k + ε ˙ φ i j k ε ˙ P l m n / λ i j k
where the single dot denotes the first-order time-difference algorithm. The inter-frequency biases for both pseduorange and phase of both receiver and satellite are very stable such that they can be deemed as constant over a short time span. Moreover, the inter-observation-type bias (between pseudorange and phase) is stable as well. Consequently, the inter-frequency biases and inter-observation-type bias can be omitted through first-order time difference algorithm [35]. According to Equation (6), the first-order ionospheric delay variation ( I ˙ 1 ) remains in the single-differenced observation of the geometry-free pseudorange minus phase combination. It can be too large to ignore under conditions of high ionospheric activity. The impact of the ionospheric variation is further discussed in Section 3. To further reduce the impact of the residual ionospheric error, the second-order time-difference algorithm can be employed.
The second-order, time-differenced, geometry-free, pseduorange minus phase combination observation equation is defined as:
φ ¨ i j k P ¨ l m n λ i j k = K i j k , l m n I ¨ 1 + N ¨ i j k + ε ¨ φ i j k ε ¨ P l m n / λ i j k
where the double dot denotes the second-order time-difference algorithm. According to Equation (7), for three consecutive epochs t 0 , t 1 , and t 2 , ignoring small residual ionospheric error and noise term, the estimated combined cycle slip on the triple-frequency phase combination φ i j k at epoch t 2 , can be defined as:
N ˙ ^ i j k = φ ¨ i j k P ¨ l m n λ i j k = [ φ i j k ( t 2 ) φ i j k ( t 1 ) ( φ i j k ( t 1 ) φ i j k ( t 0 ) ) ] [ P l m n ( t 2 ) P l m n ( t 1 ) ( P l m n ( t 1 ) P l m n ( t 0 ) ) ] / λ i j k
Assuming that the noise terms of the carrier phase observations on each frequency are independent in time and are identical in variance (i.e., σ φ 1 = σ φ 2 = σ φ 3 = σ φ ), and the same is true for the pseudorange noise terms (i.e., σ P 1 = σ P 2 = σ P 3 = σ P ), the standard deviations (STDs) of the estimated combined cycle slip N ˙ ^ i j k is calculated as:
σ N ˙ ^ i j k = 2 ( l 2 + m 2 + n 2 ) σ P 2 / λ i j k 2 + ( i 2 + j 2 + k 2 ) σ φ 2
Assuming that N ˙ ^ i j k is normally distributed, the cycle slip can be detected when the following condition is satisfied:
| φ ˙ ^ i j k P ˙ ^ l m n λ i j k | > s σ N ˙ ^ i j k
where s σ N ˙ ^ i j k is the critical value for the cycle slip detection criterion, and the scalar s denotes the multiple of the standard deviation. s = 3 and s = 4 denote a 99.7% confidence level and a 99.9% confidence level, respectively.
To detect and repair cycle slips on the three original carrier phase observations, i.e., N ˙ 1 , N ˙ 2 , N ˙ 3 , three sets of linearly independent triple-frequency phase combinations are required.
Assuming the coefficients for three linear combinations are ( i 1 , j 1 , k 1 ) , ( i 2 , j 2 , k 2 ) , and ( i 3 , j 3 , k 3 ) , the corresponding combined cycle slips are N ˙ i 1 j 1 k 1 , N ˙ i 2 j 2 k 2 , and N ˙ i 3 j 3 k 3 , respectively. In such a case, the relationship between the combined cycle slip and the original cycle slips can be defined as:
L = | N ˙ i 1 j 1 k 1 N ˙ i 2 j 2 k 2 N ˙ i 3 j 3 k 3 | = | i 1 j 1 k 1 i 2 j 2 k 2 i 3 j 3 k 3 | | N ˙ 1 N ˙ 2 N ˙ 3 | = A x
According to Equation (11), in order to ensure that the cycle slips on the three original carrier phase observations can be recovered, the A-matrix must be reversible. In such a case, the cycle slips on the three original carrier phase observations can be obtained using x = A 1 L . As the combined cycle slip N ˙ i j k is obtained by rounding the estimated combined cycle slip N ˙ ^ i j k , the elements of the L-matrix are all integers. To ensure that the cycle slips on the three original carrier phase observations that were computed with equation x = A 1 L are also integers, the elements of the inverse of the A-matrix must be integers as well. This means that the A-matrix must satisfy the condition that the matrix elements are all integers and the determinant is equal to plus or minus 1 (i.e., det ( A ) = ± 1 ).
Generally, cycle slip validation is the step that follows the cycle slip determination. Similar to Equation (10), the cycle slip validation formula is:
| φ ¨ i j k r e p a i r P ¨ l m n λ i j k | < s σ N ˙ ^ i j k
where superscript “repair” indicates that the original carrier phase observations are already corrected by subtracting the integer cycle slips. If Equation (12) is not satisfied, the determined cycle slip is considered to be false.

3. Selection of Optimal Geometry-Free Pseudorange Minus Phase Combinations

According to Equation (8), the integral combined cycle slip N ˙ i j k can be determined by rounding the estimated float combined cycle slip N ˙ ^ i j k . Consequently, a rounding success rate is introduced to evaluate the performance of the cycle slip detection.
Assuming that N ˙ ^ i j k is normally distributed,
N ˙ ^ i j k ( N ˙ i j k , σ N ˙ ^ i j k 2 )
where N ˙ i j k is the true value of the combined cycle slip. Assuming that N ˙ i j k is the nearest integer close to N ˙ ^ i j k (i.e., N ˙ i j k = int ( N ˙ ^ i j k ) ), the probability of N ˙ i j k taking an integer i is [40]:
P ( N ˙ i j k = i ) = ( i N ˙ i j k ) 0.5 ( i N ˙ i j k ) + 0.5 1 2 π σ N ˙ ^ i j k exp ( 1 2 z 2 / σ N ˙ ^ i j k 2 ) d z
Based on this, the probability of obtaining the correct combined cycle slip N ˙ i j k by rounding the estimated float combined cycle slip N ˙ ^ i j k can be computed using [40]:
P ( N ˙ i j k = N ˙ i j k ) = P ( | N ˙ ^ i j k N ˙ i j k | 0.5 ) = 2 Φ ( 0.5 σ N ˙ ^ i j k 1 ) 1
where Φ ( x ) = x 1 2 π exp ( 1 2 z 2 ) d z . It is shown in Equation (15) that the probability of rounding to the correct combined cycle slip N ˙ i j k increases as σ N ˙ ^ i j k (STD of N ˙ ^ i j k ) decreases. The relationship between the rounding success rate and σ N ˙ ^ i j k is shown in Figure 1, and the rounding success rates when σ N ˙ ^ i j k = 0.15 , 0.2 , 0.25 (cycles) are listed in Table 1.
According to Equation (6), the residual ionospheric variation still remains in first-order time-difference pseudorange minus phase combination. Taking ionospheric variation into account, the bias-affected cycle slip rounding success rate can be computed by [41]:
P b i a s = Φ ( 1 + 2 b i a s i o n 2 σ N ˙ ^ i j k ) + Φ ( 1 2 b i a s i o n 2 σ N ˙ ^ i j k ) 1
where b i a s i o n = K i j k , l m n I ˙ 1 and denotes ionospheric bias in cycles. When the between-epoch ionospheric variation I ˙ 1 is larger than 0.02 m, the bias-affected cycle slip rounding success rate will be reduced to about 97% whereas it can be larger than 99% when ionospheric bias is so small that it can be ignored. The ionospheric total electron content (TEC) change rate can be over 0.03 TECU/s during ionosphere disturbances in the low latitude region [42,43]. In such cases, the impact of residual ionospheric variation in the first-order time-difference pseudorange and phase combination cannot be ignored. However, the second-order time-difference algorithm can further reduce the impact of ionospheric variation, resulting in improved cycle slip rounding success rate.
According to Equations (7) and (10), the cycle slip detection criterion for the proposed second-order, time-difference, pseudorange minus phase combination method can be expressed as:
| φ ¨ i j k P ¨ l m n / λ i j k + f ( o ) | > s σ N ˙ ^ i j k
where f ( o ) = ε ¨ i j k + ε ¨ l m n / λ i j k K i j k , l m n I ¨ 1 .
To reduce the impact of f ( o ) and improve the rounding success rate, we propose the following three criteria for choosing the optimal pseudorange minus phase combinations:
① It should have a relatively larger wavelength ( λ i j k ) to reduce the impact of pseudorange noise.
② It should have a relatively smaller coefficient K i j k , l m n of I ¨ 1 to reduce the impact of the second-order, between-epoch, ionospheric delay variations.
③ It should have relatively smaller σ N ˙ ^ i j k to ensure a higher rounding success rate for combined cycle slip.
For criterion ②, K i j k , l m n can be expressed as:
K i j k , l m n = ( β i j k + β l m n ) / λ i j k = ( f 1 / c + f 1 β l m n / c ) [ i + j f 1 2 + f 2 2 β l m n f 1 f 2 + f 1 f 2 β l m n + k f 1 2 + f 3 2 β l m n f 1 f 3 + f 1 f 3 β l m n ]
The K i j k , l m n of the geometry-free pseudorange minus phase combinations for BDS and GPS are 12.475 ( i + 0.990 j + 0.987 k ) and 13.033 ( i + 0.982 j + 0.986 k ) , respectively. Generally, K i j k , l m n is small in cases where | i + j + k | 2 .
For criterion ③, according to Equation (9), pseudorange noise is the largest source of error in estimated combined cycle slips. To minimize the noise associated with estimated combined cycle slips (i.e., σ N ˙ ^ i j k = min ), the combination coefficients l , m , n should be l = m = n = 1 / 3 , where l 2 + m 2 + n 2 = min can be interpreted as the minimum distance from the origin point ( 0 , 0 , 0 ) to plane l + m + n = 1 .
The noise of pseudorange and carrier phase is about 1% of the code width and carrier phase wavelength respectively. This means, e.g., the noise for the GPS C/A code and P code are approximately 3 m and 0.3 m, respectively, and the noise for the carrier phase is about 0.003 m. When smoothing the code with the carrier phase, the C/A code noise might be reduced to approximately 0.6 m. In addition, the pseudorange noise also depends on the signal strength, which varies with the elevation angle. Finally, the pseudorange noise is generally larger in kinematic conditions as well. Take this into account, assuming that the carrier phase noise is σ φ = 0.003 m, the pseudorange noise is σ P = ( 0.3 , 0.6 , 3 ) m, the combined wavelength is larger than 4 m, and | i + j + k | 2 , the optimal geometry-free pseudorange minus phase linear combinations can be searched by setting i ,   j ,   and   k within the range of [ 10 , 10 ] . Table 2 and Table 3 list the first seven optimal pseudorange minus phase combinations for BDS and GPS, respectively. Columns 1 through 3 are coefficients of pseudorange minus phase combinations. Columns 4 and 5 are combined wavelength and ionospheric scale factors, respectively. Columns 6 through 8 are STDs of estimated combined cycle slip where σ P = 0.3 , 0.6 ,   and   3 m, respectively.
Table 2 and Table 3 show that the listed pseudorange minus phase linear combinations are not affected as much by the pseudorange noise due to their larger wavelength. The STDs of estimated combined cycle slip are generally smaller than 0.2 cycles in the case of 0.3 m and 0.6 m pseudorange noise, which results in at least a 98.8% rounding success rate for cycle slip detection. The STDs can be less than 0.15 cycles in certain cases, which can result in at least a 99.9% rounding success rate. However, the STDs for all the pseudorange minus phase linear combinations are larger than 0.2 cycles in the case of 3 m pseudorange noise. For combinations with wavelengths that are less than 10 m, the impact of pseudorange noise can be significant, and the STDs can be larger than 0.3 cycles. In such cases, the rounding success rate can be relatively small.
It should be noted that the impact of ionospheric delay can be effectively reduced by using a second-order time difference of pseudorange minus phase combinations, but at the same time, the observation noise is significantly amplified. The impact of observation noise can be reduced by choosing pseudorange minus phase combinations with larger wavelengths. In the case of low pseudorange noise, combinations with small STDs can be derived. Contrary to this, the combined noise can be relatively large in the case of large pseudorange noise, which makes it unsuitable for cycle slip detection.
Table 4 and Table 5 list the optimal pseudorange minus phase combinations for BDS and GPS cycle slip detection and repair, respectively. Columns 1 through 3 list the coefficients of the first combination, columns 4 through 6 list the coefficients of the second combination, and columns 7 through 9 list the coefficients of the third combination. The final three columns are the minimum probability for cycle slip detection when σ P = 0.3 , 0.6 ,   and   3 m, respectively. Table 4 and Table 5 show that the cycle slip detection success rates can reach up to 99.9% and 99.2% for BDS and GPS if the pseudorange noise is 0.3 m and 0.6 m, respectively. However, the highest cycle slip detection success rate is only approximately 95% if the pseudorange noise is 3 m.

4. Data Tests and Analysis

4.1. Static Test

According to [44], a moderate geomagnetic storm occurred in the middle of March 2013. Figure 2 shows the geomagnetic Kp index from 17 March 2013. The average Kp index is 5.25, which also indicates that the level of ionospheric activity was high. The Earth’s magnetic field was disturbed dramatically in the near-equatorial, mid-latitude, and auroral regions. Triple-frequency BDS/GPS observations collected from the multi-GNSS experiment (MGEX) station JFNG (30.52° N, 114.49° E)—which is in the area where the Total Election Content (TEC) changed significantly on 17 March 2013—were used to verify the performance of the proposed method. The observations were collected with a 30 s sampling interval using a Trimble NetR9 receiver, which was connected to a TRM59800 antenna. BDS C03, C09, C12 and GPS G01, G24, G25 were tested, and simulated cycle slips were previously added to their corresponding observations.
We first analyzed ionospheric delay variations. The first-order, between-epoch ionospheric delay variations can be estimated using I ˙ = f 3 2 ( λ 1 φ ˙ 1 λ 3 φ ˙ 3 ) / ( f 1 2 f 3 2 ) , and the second-order, between-epoch time-difference ionospheric delay variations can be calculated using I ¨ i = I ˙ i I ˙ i 1 .
Figure 3 displays the between-epoch, first-order and second-order ionospheric delay variations of BDS C03, C09, C12 and GPS G01, G24, G25. The first-order ionospheric delay variations for all the involved satellites, except for C03, were fairly significant and cannot be ignored due to the impact of the magnetic storm. This is reasonable considering that the signal propagation path of C03 in the ionosphere varies much more slowly than the other satellites due to its slower movement in the Geostationary Earth Orbit (GEO). However, when I ˙ is further differenced in time, the trend is removed and the second-order ionospheric delay variations (i.e., I ¨ ) for all the involved satellites generally vary within ±2 cm. This indicates that the impact of ionospheric delay can be significantly reduced using the between-epoch, second-order, time-difference algorithm. Consequently, the rounding success rate for the estimated float combined cycle slips using the second-order time-difference algorithm can be significantly improved compared to the first-order time-difference algorithm.
According to Table 4 and Table 5, the optimal geometry-free pseudorange minus phase combinations were [−4, 1, 4], [−3, 6, −2], [4, −2, −3] for BDS and [−6, 1, 7], [3, 0, −4], [4, −8, 3] for GPS, respectively. To evaluate the effectiveness of the proposed method, two types of evenly distributed cycle slips were simulated.
① The combined cycle slips satisfy the condition of | i N ˙ 1 + j N ˙ 2 + k N ˙ 3 | 1 to verify the performance in identifying small combined cycle slips.
② The cycle slips for each of the original carrier phase observations were defined as 1 cycle to verify the ability of identifying small cycle slips.
First, we used three independent, geometry-free, pseudorange minus phase combinations to detect and repair the cycle slips on the original carrier phase observations using Equation (6) and ignored the between-epoch ionospheric delay variations. According to Equation (6), when I ˙ is sufficiently large, misjudgment may occur even if no cycle slips occur on the original carrier phase observations (i.e., N ˙ i j k = 0 ( N ˙ 1 = N ˙ 2 = N ˙ 3 = 0 )). Consequently, the performance will be reduced under conditions of high ionospheric activity. The cycle slip detection results for raw observations are shown in Figure 4. It is apparent that misjudgments occur in all cases due to the presence of relatively large between-epoch ionospheric delay variations.
Second, the proposed method was tested with data that were inserted with simulated cycle slips. Cycle slip detection and repair results are shown in Figure 5 and Table 6. It is apparent that the proposed method can correctly detect and repair simulated cycle slips that are as small as 1 cycle. Moreover, there were no insensitive cycle slip pairs that occurred in response to the proposed method. For example, although BDS cycle slips pair (5, 4, 4) cannot be identified by geometry-free pseudorange minus phase combinations [−4, 1, 4] and [4, −2, −3], it was obviously identified by [−3, 6, 2]. Consequently, it can still be uniquely determined by the three independent cycle slip detection combinations. Finally, the simulated cycle slips can be correctly identified even under conditions of high levels of ionospheric activity. Using C09 for example, as shown in Figure 5, the between-epoch ionospheric delay variations can exceed 0.1 m from GPST 9:00 to 10:00. However, after applying the second-order time difference, the ionospheric variations are close to zero, and the simulated cycles-slips at GPST 09:59:30 can be correctly detected and repaired (as shown in Figure 5 and Table 6).

4.2. Kinematic Test

In this section, BDS observations from the car-driven experiment and GPS observations from the airborne experiment are used to further evaluate the real-time cycle slip detection and repair performance of the proposed method under kinematic conditions. The ionosphere is quiet for both the car-driven experiment and airborne experiment, and the simulated cycle slips are previously inserted to observations. To reduce the impact of relatively large pseudorange noise in the kinematic environments, geometry-free pseudorange minus phase combinations with relatively large wavelengths were chosen. Consequently, according to Table 2 through Table 5, the chosen BDS cycle slip detection combinations were [−4, 1, 4], [−3, 6, −2], and [4, −2, −3], which had combined wavelengths of 8.14 m, 13.321 m, and 12.211 m, respectively. The chosen GPS cycle slip detection combinations were [−6, 1, 7], [3, 0, −4], and [4, −8, 3], which had combined wavelengths of 29.305 m, 14.653 m, and 29.305 m, respectively.

4.2.1. Kinematic Test for BDS

The BDS kinematic data were collected from the car-driven experiment conducted in Wuhan on 27 October 2016. The sampling interval was 0.5 s. The type of GNSS receiver and antenna that were used in the experiment were a ComNav M300 receiver and a NovAtel 704-WB antenna. The observational environments and GNSS antennas in car-driven experiment are shown in Figure 6. The trajectory of the car is shown in Figure 7. The speed of car was approximately 50 km/h. Observations from C01, C06, and C14 were analyzed to verify the cycle slip detection and repair results of the proposed method and are shown in Figure 8, Figure 9 and Figure 10 and in Table 7.

4.2.2. Kinematic Test for GPS

The GPS kinematic data were collected from the airborne experiment conducted at Xi’an in China on 24 April 2015. The sampling interval was 0.2 s. The type of GNSS receiver and antenna that were used in the experiment were a NovAtel ProPak6 receiver and an Antcom aviation GNSS antenna (743GNSSA-XT-1). The airplane and antenna used in the airborne experiment are shown in Figure 11. The trajectory of the airplane is shown in Figure 12. The speed of airplane was approximately 300 km/h. Observations from G03, G09, and G26 were analyzed to verify the cycle slip detection and repair results of the proposed method and are shown in Figure 13, Figure 14 and Figure 15 and in Table 8.

4.2.3. Analysis of Results

It is apparent in Figure 8, Figure 9 and Figure 10 and Figure 13, Figure 14 and Figure 15 that the proposed method can also be applied to process kinematic observations to effectively identify and repair cycle slips that are as small as 1 cycle.
The pseudorange noise observed during the kinematic experiments is generally larger than that observed during the static experiments. According to Equation (7), the impact of pseudorange noise decreases when divided by a larger combined wavelength. Optimal combinations with larger wavelengths can be constructed by using triple-frequency observations compared to dual-frequency observations. As such, the availability of triple-frequency observations plays an important role in cycle slip detection and repair, especially for preprocessing observations collected from kinematic experiments.

4.3. Summary of Static and Kinematic Test Results

In summary, the proposed method can detect and repair as small as 1 cycle slip for both static and kinematic observations. Moreover, the cycle slip detection possibility is improved under high ionospheric activity, and there are no insensitive cycle slips existing in our method.

5. Conclusions

In this paper, we present a new real-time cycle slip detection and repair method for undifferenced GPS/BDS triple-frequency observations collected by a single receiver under conditions of high ionospheric activity. In the approach, the cycle slips on the original triple-frequency carrier phase observations were correctly and uniquely determined by applying the second-order time difference algorithm to three linearly independent geometry-free pseudorange minus phase combinations. The proposed approach has multiple advantages. First, no insensitive cycle slips exist in our method. Second, the impact of ionospheric delay variations can be significantly reduced by using the second-order time-difference algorithm even if the between-epoch ionospheric delay varies rapidly, resulting in an improved cycle slip rounding success rate. Finally, triple-frequency carrier phase and pseudorange observations of only three consecutive epochs are required for data processing. Therefore, the method enables the implementation of high-efficiency cycle slip detection and repair in real time.
The performance of the method was verified by both static and kinematic tests. It was shown that the method can detect and correctly repair cycle slips as small as 1 cycle under conditions of high ionospheric activity. Additionally, the method is effective in processing kinematic data. However, the performance of the method may be influenced by the pseudorange noise level. The cycle slip combinations can be determined with a rounding success rate of more than 99% in the case of the 0.3 m and 0.6 m pseudorange noise, whereas the success rate is only approximately 95% in the case of the 3 m pseudorange noise. These results indicate that the cycle slip detection performance of the presented method is reduced when pseudorange noise is particularly large (e.g., 3 m or larger).

Acknowledgments

The authors gratefully acknowledge the support of the National Key R&D Program of China (No. 2016YFB0501803), the National Natural Science Foundation of China (No. 91638203, 41774031), the Natural Science Foundation for Distinguished Young Scholars of Hubei Province of China (No. 2015CFA039), and LIESMARS Special Research Funding. Constructive comments and suggestions from the anonymous reviewers are acknowledged as well.

Author Contributions

W.L. and X.J. provided the initial idea and conception; M.W. analyzed the feasibility of the proposed method in theory; J.H. and Y.W. analyzed the quality of collected data; X.J. implemented the software and carried out the experimental validation; W.L., M.W., X.J., J.H. and Y.W. wrote the paper. All authors reviewed the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Beutler, G.; Davidson, D.A.; Langley, R.B.; Santerre, R.; Vanicek, P.; Wells, D.E. Some Theoretical and Practical Aspects of Geodetic Positioning Using Carrier Phase Difference Observations of GPS Satellites; Surveying Engineering, University of New Brunswick: Fredericton, NB, Canada, 1985. [Google Scholar]
  2. Hofmann, B.; Lichtenegger, H.; Collins, J. GPS Theory and Practice; Springer Wien: New York, NY, USA, 2001. [Google Scholar]
  3. Qian, C.; Liu, H.; Zhang, M.; Shu, B.; Xu, L.; Zhang, R. A geometry-based cycle slip detection and repair method with time-differenced carrier phase (TDCP) for a single frequency global position system (GPS)+ BeiDou navigation satellite system (BDS) receiver. Sensors 2016, 16, 2064. [Google Scholar] [CrossRef] [PubMed]
  4. Fujita, S.; Saito, S.; Yoshihara, T. Cycle slip detection and correction methods with time-differenced model for single frequency GNSS applications. Trans. Inst. Syst. Control Inf. Eng. 2013, 26, 8–15. [Google Scholar] [CrossRef]
  5. Kirkko-Jaakkola, M.; Traugott, J.; Odijk, D.; Collin, J.; Sachs, G.; Holzapfel, F. A raim approach to GNSS outlier and cycle slip detection using L1 carrier phase time-differences. In Proceedings of the 2009 IEEE Workshop on Signal Processing Systems (SiPS 2009), Tampere, Finland, 7–9 October 2009; pp. 273–278. [Google Scholar]
  6. Song, W.; Yao, Y. Pre-process strategy in complex movement using single-frequency GPS data. Geomat. Inf. Sci. Wuhan Univ. 2009, 34, 1005–1008. [Google Scholar]
  7. Zhang, Q.; Gui, Q.; Li, J.; Gong, Y.; Han, S. Bayesian methods for cycle slips detection based on autoregressive model. In China Satellite Navigation Conference (CSNC) 2012 Proceedings; Springer: Berlin/Heidelberg, Germany, 2012; pp. 317–335. [Google Scholar]
  8. Zhang, Q.; Gui, Q. Bayesian methods for outliers detection in GNSS time series. J. Geodesy 2013, 87, 609–627. [Google Scholar]
  9. Zangeneh-Nejad, F.; Amiri-Simkooei, A.; Sharifi, M.; Asgari, J. Cycle slip detection and repair of undifferenced single-frequency GPS carrier phase observations. GPS Solut. 2017, 21, 1593–1603. [Google Scholar] [CrossRef]
  10. Gao, Y.; Li, Z. Cycle slip detection and ambiguity resolution algorithms for dual-frequency GPS data processing. Mar. Geodesy 1999, 22, 169–181. [Google Scholar]
  11. Bisnath, S.B.; Langley, R.B. Efficient, automated cycle-slip correction of dual-frequency kinematic GPS data. In Proceedings of the 13th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GPS 2000), Salt Lake City, UT, USA, 19–22 September 2000; pp. 145–154. [Google Scholar]
  12. Kim, D.; Langley, R.B. Instantaneous real-time cycle-slip correction of dual frequency GPS data. In Proceedings of the International Symposium on Kinematic Systems in Geodesy, Geomatics and Navigation, Banff, AB, Canada, 5–8 June 2001; pp. 255–264. [Google Scholar]
  13. Kim, D.; Langley, R.B. Instantaneous real-time cycle-slip correction for quality control of GPS carrier-phase measurements. Navigation 2002, 49, 205–222. [Google Scholar] [CrossRef]
  14. Blewitt, G. An automatic editing algorithm for GPS data. Geophys. Res. Lett. 1990, 17, 199–202. [Google Scholar] [CrossRef]
  15. Hatch, R. The synergism of GPS code and carrier measurements. In Proceedings of the International Geodetic Symposium on Satellite Doppler Positioning, Las Cruces, NM, USA, 8–12 February 1982; pp. 1213–1231. [Google Scholar]
  16. Melbourne, W.G. The case for ranging in GPS-based geodetic systems. In Proceedings of the First International Symposium on Precise Positioning with the Global Positioning System, Rockville, MD, USA, 15–19 April 1985. [Google Scholar]
  17. Wübbena, G. Software developments for geodetic positioning with GPS using TI-4100 code and carrier measurements. In Proceedings of the First International Symposium on Precise Positioning with the Global Positioning System, Rockville, MD, USA, 15–19 April 1985; Volume 19. [Google Scholar]
  18. De Lacy, M.C.; Reguzzoni, M.; Sansò, F.; Venuti, G. The bayesian detection of discontinuities in a polynomial regression and its application to the cycle-slip problem. J. Geodesy 2008, 82, 527–542. [Google Scholar] [CrossRef]
  19. Liu, Z. A new automated cycle slip detection and repair method for a single dual-frequency GPS receiver. J. Geodesy 2010, 85, 171–183. [Google Scholar] [CrossRef]
  20. Cai, C.; Liu, Z.; Xia, P.; Dai, W. Cycle slip detection and repair for undifferenced GPS observations under high ionospheric activity. GPS Solut. 2012, 17, 247–260. [Google Scholar] [CrossRef]
  21. Zhang, X.; Li, X. Instantaneous re-initialization in real-time kinematic PPP with cycle slip fixing. GPS Solut. 2012, 16, 315–327. [Google Scholar]
  22. Banville, S.; Langley, R.B. Mitigating the impact of ionospheric cycle slips in GNSS observations. J. Geodesy 2013, 87, 179–193. [Google Scholar] [CrossRef]
  23. Lee, H.K.; Wang, J.; Rizos, C. Effective cycle slip detection and identification for high precision GPS/INS integrated systems. J. Navig. 2003, 56, 475–486. [Google Scholar] [CrossRef]
  24. Kim, Y.; Song, J.; Kee, C.; Park, B. GPS cycle slip detection considering satellite geometry based on TDCP/INS integrated navigation. Sensors 2015, 15, 25336–25365. [Google Scholar] [CrossRef] [PubMed]
  25. Du, S.; Gao, Y. Inertial aided cycle slip detection and identification for integrated PPP GPS and INS. Sensors 2012, 12, 14344–14362. [Google Scholar] [CrossRef] [PubMed]
  26. Cocard, M.; Bourgon, S.; Kamali, O.; Collins, P. A systematic investigation of optimal carrier-phase combinations for modernized triple-frequency GPS. J. Geodesy 2008, 82, 555–564. [Google Scholar] [CrossRef]
  27. Richert, T.; El-Sheimy, N. Optimal linear combinations of triple frequency carrier phase data from future global navigation satellite systems. GPS Solut. 2006, 11, 11–19. [Google Scholar] [CrossRef]
  28. Dai, Z.; Knedlik, S.; Loffeld, O. Real-Time Cycle-Slip Detection and Determination for Multiple Frequency GNSS. In Proceedings of the 2008 5th Workshop on Positioning, Navigation and Communication (WPNC 2008), Hannover, Germany, 27 March 2008; pp. 37–43. [Google Scholar]
  29. Teunissen, P.J. The least-squares ambiguity decorrelation adjustment: A method for fast GPS integer ambiguity estimation. J. Geod. 1995, 70, 65–82. [Google Scholar] [CrossRef]
  30. Huang, L.; Lu, Z.; Zhai, G.; Ouyang, Y.; Huang, M.; Lu, X.; Wu, T.; Li, K. A new triple-frequency cycle slip detecting algorithm validated with BDS data. GPS Solut. 2015, 20, 761–769. [Google Scholar] [CrossRef]
  31. Li, J.; Yang, Y.; Xu, J. Real-time cycle-slip detection and repair based on code-phase combinations for GNSS triple-frequency un-differenced observations. Acta Geod. Cartogr. Sin. 2011, 40, 717–722. [Google Scholar]
  32. De Lacy, M.C.; Reguzzoni, M.; Sansò, F. Real-time cycle slip detection in triple-frequency GNSS. GPS Solut. 2011, 16, 353–362. [Google Scholar] [CrossRef]
  33. Zhao, Q.; Sun, B.; Dai, Z.; Hu, Z.; Shi, C.; Liu, J. Real-time detection and repair of cycle slips in triple-frequency GNSS measurements. GPS Solut. 2014, 19, 381–391. [Google Scholar] [CrossRef]
  34. Yao, Y.; Gao, J.; Wang, J.; Hu, H.; Li, Z. Real-time cycle-slip detection and repair for BeiDou triple-frequency undifferenced observations. Surv. Rev. 2016, 48, 367–375. [Google Scholar] [CrossRef]
  35. Li, B.; Qin, Y.; Li, Z.; Lou, L. Undifferenced cycle slip estimation of triple-frequency BeiDou signals with ionosphere prediction. Mar. Geodesy 2016, 39, 348–365. [Google Scholar] [CrossRef]
  36. Huang, L.; Zhai, G.; Ouyang, Y. Ionospheric cycle slip processing in triple-frequency GNSS. Acta Geod. Cartogr. Sin. 2015, 44, 717–725. [Google Scholar]
  37. Zhang, X.; Li, P. Benefits of the third frequency signal on cycle slip correction. GPS Solut. 2016, 20, 451–460. [Google Scholar] [CrossRef]
  38. Eueler, H.J.; Goad, C.C. On optimal filtering of GPS dual frequency observations without using orbit information. Bull. Geodesique 1991, 65, 130–143. [Google Scholar] [CrossRef]
  39. Feng, Y. GNSS three carrier ambiguity resolution using ionosphere-reduced virtual signals. J. Geodesy 2008, 82, 847–862. [Google Scholar] [CrossRef]
  40. Teunissen, P.J. Success probability of integer GPS ambiguity rounding and bootstrapping. J. Geodesy 1998, 72, 606–612. [Google Scholar] [CrossRef]
  41. Teunissen, P.J. Integer estimation in the presence of biases. J. Geodesy 2001, 75, 399–407. [Google Scholar] [CrossRef]
  42. Basu, S.; Groves, K.M.; Quinn, J.M.; Doherty, P. A comparison of TEC fluctuations and scintillations at Ascension Island. J. Atmos. Sol.-Terr. Phys. 1999, 61, 1219–1226. [Google Scholar] [CrossRef]
  43. Liu, Z.; Chen, W. Study of the ionospheric tec rate in Hong Kong region and its GPS/GNSS application. In Proceedings of the International Technical Meeting on GNSS Global Navigation Satellite System-Innovation and Application, Beijing, China, 8–9 August 2009; pp. 129–137. [Google Scholar]
  44. Jin, R.; Jin, S.; Tao, X. Ionospheric anomalies during the March 2013 geomagnetic storm from BeiDou Navigation Satellite System (BDS) observations. In China Satellite Navigation Conference (CSNC) 2014 Proceedings: Volume I; Springer: Berlin/Heidelberg, Germany, 2014; pp. 97–104. [Google Scholar]
Figure 1. Relationship between rounding success rate and σ N ˙ ^ i j k .
Figure 1. Relationship between rounding success rate and σ N ˙ ^ i j k .
Sensors 18 00427 g001
Figure 2. Geomagnetic Kp index from 17 March 2013.
Figure 2. Geomagnetic Kp index from 17 March 2013.
Sensors 18 00427 g002
Figure 3. Between-epoch first-order and second-order ionospheric delay variations for C03, C09, C12, G01, G24, and G25.
Figure 3. Between-epoch first-order and second-order ionospheric delay variations for C03, C09, C12, G01, G24, and G25.
Sensors 18 00427 g003
Figure 4. Cycle slip detection by first-order between-epoch difference pseudorange minus phase combination for C03, C09, C12, G01, G24, and G25.
Figure 4. Cycle slip detection by first-order between-epoch difference pseudorange minus phase combination for C03, C09, C12, G01, G24, and G25.
Sensors 18 00427 g004
Figure 5. Cycle slip detection by second-order between-epoch difference pseudorange minus phase combination for C03, C09, C12, G01, G24, and G25.
Figure 5. Cycle slip detection by second-order between-epoch difference pseudorange minus phase combination for C03, C09, C12, G01, G24, and G25.
Sensors 18 00427 g005
Figure 6. Observational environments and GNSS antennas in the car-driven experiment.
Figure 6. Observational environments and GNSS antennas in the car-driven experiment.
Sensors 18 00427 g006
Figure 7. Trajectory of the car.
Figure 7. Trajectory of the car.
Sensors 18 00427 g007
Figure 8. Cycle slip detection for C01 with sampling interval of 0.5 s.
Figure 8. Cycle slip detection for C01 with sampling interval of 0.5 s.
Sensors 18 00427 g008
Figure 9. Cycle slip detection for C06 with sampling interval of 0.5 s.
Figure 9. Cycle slip detection for C06 with sampling interval of 0.5 s.
Sensors 18 00427 g009
Figure 10. Cycle slip detection for C14 with sampling interval of 0.5 s.
Figure 10. Cycle slip detection for C14 with sampling interval of 0.5 s.
Sensors 18 00427 g010
Figure 11. (a) Airplane used in the airborne experiment; (b) GNSS antenna used in the airborne experiment.
Figure 11. (a) Airplane used in the airborne experiment; (b) GNSS antenna used in the airborne experiment.
Sensors 18 00427 g011
Figure 12. Trajectory of the airplane.
Figure 12. Trajectory of the airplane.
Sensors 18 00427 g012
Figure 13. Cycle slip detection for G03 with sampling interval of 0.2 s.
Figure 13. Cycle slip detection for G03 with sampling interval of 0.2 s.
Sensors 18 00427 g013
Figure 14. Cycle slip detection for G09 with sampling interval of 0.2 s.
Figure 14. Cycle slip detection for G09 with sampling interval of 0.2 s.
Sensors 18 00427 g014
Figure 15. Cycle slip detection for G26 with sampling interval of 0.2 s.
Figure 15. Cycle slip detection for G26 with sampling interval of 0.2 s.
Sensors 18 00427 g015
Table 1. Rounding success rates for σ N ˙ ^ i j k = 0.15 , 0.2 , 0.25 .
Table 1. Rounding success rates for σ N ˙ ^ i j k = 0.15 , 0.2 , 0.25 .
σ N ˙ ^ i j k (Cycles)Rounding Success Rate
0.1599.9%
0.298.8%
0.2595.5%
Table 2. BeiDou Navigation Satellite System (BDS) geometry-free pseudorange minus phase linear combinations.
Table 2. BeiDou Navigation Satellite System (BDS) geometry-free pseudorange minus phase linear combinations.
i j k λ ( m ) K STD (Cycles)
σ P = 0.3   m σ P = 0.6   m σ P = 3   m
−4148.1411.710.1230.1430.441
−36−213.32112.0710.1420.1490.295
−1−5620.932−0.3620.1580.1610.228
0−114.884−0.040.0760.1450.71
4−2−312.211−11.750.1110.1220.303
53−929.305−11.3880.2150.2160.245
7−8−1146.526−23.8210.2140.2140.215
STD, standard deviation.
Table 3. Global Positioning System (GPS) geometry-free pseudorange minus phase linear combinations.
Table 3. Global Positioning System (GPS) geometry-free pseudorange minus phase linear combinations.
i j k λ ( m ) K STD (Cycles)
σ P = 0.3   m σ P = 0.6   m σ P = 3   m
−61729.30524.5250.1860.1870.22
−3139.76812.2420.0940.1120.365
−18−729.305−0.5130.2140.2150.244
01−15.861−0.0410.0660.1220.592
30−414.653−12.2830.1030.1110.257
4−8329.305−11.770.1890.190.223
7−8−19.768−24.0520.2160.2250.414
Table 4. Pseudorange minus phase combinations for BDS cycle slip detection and repair.
Table 4. Pseudorange minus phase combinations for BDS cycle slip detection and repair.
i 1 j 1 k 1 i 2 j 2 k 2 i 3 j 3 k 3 P σ P = 0.3   m P σ P = 0.6   m P σ P = 3   m
−414−36−20−1199.96%99.92%51.88%
−414−36−24−2−399.96%99.92%74.33%
−36−20−114−2−399.96%99.92%51.88%
−1−5653−97−8−198.01%97.95%95.88%
Table 5. Pseudorange minus phase combinations for GPS cycle slip detection and repair.
Table 5. Pseudorange minus phase combinations for GPS cycle slip detection and repair.
i 1 j 1 k 1 i 2 j 2 k 2 i 3 j 3 k 3 P σ P = 0.3   m P σ P = 0.6   m P σ P = 3   m
−61730−44−8399.18%99.14%94.86%
−31330−44−8399.18%99.14%82.91%
−617−3134−8399.18%99.14%82.91%
−617−18−74−8398.06%98.00%95.95%
Table 6. Cycle slip detection and repair results.
Table 6. Cycle slip detection and repair results.
PRNTime (HH:MM:SS)Simulated Cycle Slip N ˙ ^ i 1 j 1 k 1 N ˙ ^ i 2 j 2 k 2 N ˙ ^ i 3 j 3 k 3 Computed Integer Cycle Slip
C0304:09:30(1,1,1)1.020.96−1.02(1,1,1)
12:29:30(5,4,4)−0.021.02−0.01(5,4,4)
20:49:30(22,17,18)1.020.01−0.06(22,17,18)
C0901:39:30(1,1,1)1.041.01−1.03(1,1,1)
05:49:30(5,4,4)−0.030.940.01(5,4,4)
09:59:30(22,17,18)1.050.09−0.01(22,17,18)
C1213:19:30(1,1,1)1.181.08−1.10(1,1,1)
15:49:30(5,4,4)0.041.07−0.06(5,4,4)
18:19:30(22,17,18)0.80−0.040.08(22,17,18)
G0114:09:30(1,1,1)1.77−0.91−0.91(1,1,1)
15:49:30(5,4,4)0.08−0.200.98(5,4,4)
16:39:30(22,17,18)−1.011.03−1.08(22,17,18)
G2401:39:30(1,1,1)2.09−1.06−1.04(1,1,1)
03:44:30(4,3,3)−0.020.010.97(4,3,3)
05:49:30(23,18,17)−1.011.03−1.08(23,18,17)
G2506:39:30(1,1,1)1.96−1.03−1.13(1,1,1)
07:54:30(4,3,3)0.07−0.080.95(4,3,3)
09:09:30(23,18,17)−1.141.08−0.91(23,18,17)
PRN, pseudorandom noise code.
Table 7. Cycle slip detection and repair results.
Table 7. Cycle slip detection and repair results.
PRNTime (HH:MM:SS)Simulated Cycle Slip N ˙ ^ i 1 j 1 k 1 N ˙ ^ i 2 j 2 k 2 N ˙ ^ i 3 j 3 k 3 Computed Integer Cycle Slip
C0102:53:40(1,1,1)0.970.98−0.97(1,1,1)
03:00:20(5,4,4)−0.041.140.02(5,4,4)
03:07:00(22,17,18)1.050.02−0.03(22,17,18)
C0602:53:40(1,1,1)1.011.08−1.03(1,1,1)
03:00:20(5,4,4)−0.051.030.04(5,4,4)
03:07:00(22,17,18)1.130.04−0.17(22,17,18)
C1402:53:40(1,1,1)1.031.00−1.03(1,1,1)
03:00:20(5,4,4)−0.011.000.00(5,4,4)
03:07:00(22,17,18)1.11−0.04−0.08(22,17,18)
Table 8. Cycle slip detection and repair results.
Table 8. Cycle slip detection and repair results.
PRNTime (HH:MM:SS)Simulated Cycle Slip N ˙ ^ i 1 j 1 k 1 N ˙ ^ i 2 j 2 k 2 N ˙ ^ i 3 j 3 k 3 Computed Integer Cycle Slip
G0305:13:51(1,1,1)1.81−0.89−1.08(1,1,1)
05:16:31(4,3,3)−0.060.030.99(4,3,3)
05:19:11(23,18,17)−0.850.91−0.92(23,18,17)
G0906:00:31(1,1,1)2.01−1.01−1.00(1,1,1)
06:03:11(4,3,3)−0.210.120.90(4,3,3)
06:05:51(23,18,17)−1.171.10−1.06(23,18,17)
G2605:13:51(1,1,1)1.80−0.88−1.09(1,1,1)
05:16:31(4,3,3)0.010.000.99(4,3,3)
05:19:11(23,18,17)−1.031.02−1.02(23,18,17)

Share and Cite

MDPI and ACS Style

Liu, W.; Jin, X.; Wu, M.; Hu, J.; Wu, Y. A New Real-Time Cycle Slip Detection and Repair Method under High Ionospheric Activity for a Triple-Frequency GPS/BDS Receiver. Sensors 2018, 18, 427. https://doi.org/10.3390/s18020427

AMA Style

Liu W, Jin X, Wu M, Hu J, Wu Y. A New Real-Time Cycle Slip Detection and Repair Method under High Ionospheric Activity for a Triple-Frequency GPS/BDS Receiver. Sensors. 2018; 18(2):427. https://doi.org/10.3390/s18020427

Chicago/Turabian Style

Liu, Wanke, Xueyuan Jin, Mingkui Wu, Jie Hu, and Yun Wu. 2018. "A New Real-Time Cycle Slip Detection and Repair Method under High Ionospheric Activity for a Triple-Frequency GPS/BDS Receiver" Sensors 18, no. 2: 427. https://doi.org/10.3390/s18020427

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop