Cycle slip detection and repair for BeiDou-3 triple-frequency signals

When carrier phase observations are applied to high-precision positioning, how to handle the cycle slip is an unavoidable problem. For cycle slip correction, detection combination noise and the ionospheric delay are two crucial factors. Specifically, the drastic changes in the ionosphere and the increased noise of code observations will increase the failure probability of cycle slip detection. To reduce the influence of code observation noise and ionospheric bias, a novel cycle slip detection method for BDS-3 satellites is proposed. Considering that code measurement noise is closely related to the satellite elevation angle, an elevation-based model is built to evaluate the code measurement noise. Firstly, two modified code-phase combinations are selected optimally based on 1% missed detection rate and 99% success detection rate to minimize the effects of code measurement noise. However, the second modified code-phase combination is more affected by ionospheric delay bias, so ionospheric bias of current epoch needs to be corrected. To reduce the influence of ionospheric bias, two moving windows of time-differenced ionospheric delay are introduced to correct the ionospheric bias of the second code-phase combination. Experiments with BeiDou-3 data are implemented in three different scenarios. To verify the effectiveness of the algorithm in the environment of high code observations noise, Gaussian noise is added to the code observations in the first scenario, and the results demonstrate that the success rate of cycle slip detection and repair is still greater than 95% when the standard deviation of Gaussian noise is 0.8 m. The second scenario is carried out under low ionospheric activity, and results indicate that the proposed method significantly reduces the times of failed detection and repair. Moreover, in the third scenario, BeiDou-3 data with cycle slips of different types under high ionospheric activity are tested, and all cycle slips can be correctly detected and corrected.


Introduction
Global navigation satellite system (GNSS) plays an important role in military confrontations, 1 transportation, 2-4 and measurement. 5 The Chinese BeiDou Navigation Satellite System (BDS) includes regional system (BDS-2) and global system (BDS-3). BDS-2 has already been accomplished in 2012 and triple-frequency signals can be transmitted by all satellites. 6,7 BDS-3 can provide global services and it is expected to be completely built in 2020. The release of four signals including B1I, B3I, B1C, and B2a is beneficial to the ambiguity resolution, cycle slip detection, and highprecision positioning. 6 Millimeter precision is the advantage of carrier phase observations compared to code observations. However, the existence of cycle slips is an unavoidable problem in the application of carrier phase observations. Cycle slips often occur in poor environment, such as signal interruption and interference. The quality of satellite data preprocessing is directly related to the accuracy and reliability of navigation and positioning. If cycle slips are ignored or cannot be correctly repaired, it will bring deviations to the ultimate results.
For cycle slip detection, the basic idea is to obtain a sequence of detections that reflects the change of the cycle slip, and the position of the cycle slip can be determined from the detection sequence. For single-frequency observations, Beutler et al. 8 proposed a polynomial fitting method to detect cycle slip. Meanwhile Hofmann et al. 9 put forward a high-order inter-epoch phase differentiation method. However, they need several observations from previous epochs for detection. Zangeneh-Nejad et al. 10 improved the single-frequency GPS cycle slip correction method, which is based on the generalized likelihood ratio test. Momoh et al. 11 proposed to detect cycle slip through the adaptive difference sequences of observations, which can achieve single-frequency cycle slip detection and correction with the slips of one cycle.
For GPS dual-frequency data, Hatch-Melbourne-Wübbena (HMW) combinations [12][13][14] are extensively applied in cycle slip detection due to their excellent characteristics. Blewitt 15 developed the Turbo Edit algorithm and used undifferenced observations to detect cycle slip. Specifically, Turbo Edit algorithm consists of ionospheric combination and HMW combination. Considering the gross error of code observations may result in a wrong detection result, a method is proposed by Li et al., 16 which is relied on satellite orbit and smoothed code observations. Ionospheric total electron contents rate (TECR) method was proposed by Liu,17 which was based on smoothly changing ionosphere. However, the TECR method does not work well under high ionospheric activity. To solve this problem, Cai et al. 18 proposed improved TurboEdit method, which based on the second-order, time-difference phase ionospheric residual (STPIR) algorithm, and the proposed approach was able to determine the cycle slip of each signal under high ionospheric activities. Banville and Langley 19 developed the dual-frequency model, and least-squares adjustment was adopted to reduce the impact of the ionosphere in the case of ionospheric disturbance. Hu et al. 20 improved Turbo Edit algorithm and a polynomial function to fit the ionospheric delays in short times.
Up to now, many satellites have supported the transmission of triple-frequency signals, which benefits ambiguity resolution 21,22 and cycle slip detection. Triple-frequency signals can present more high-performance observation combinations with longer wavelengths, lower noise, and ionospheric errors. For GPS triple-frequency observations, Dai et al. 23 constructed geometry-free (GF) phase observations as detection sequence. De Lacy et al. 24 proposed a real-time detection method and five first-order timedifference GF linear combinations were selected to detect different types of cycle slips. Zhao et al. 25 and Gu and Zhu 26 extended three carrier ambiguity resolution algorithm from ambiguity resolution to cycle slip correction, and extra-wide-lane, wide-lane, and narrow-lane (NL) combinations are utilized to determine cycle slip.
For BDS-2 triple-frequency cycle slip detection, three GF combinations are used by Huang et al. 27 to determine the position of the cycle slip, and the size of the cycle slip will be determined by least-squares ambiguity decorrelation adjustment. 28 Xiao et al. 29 modified the geometrybased detection model and a new strategy was put forward to reduce false alarms. Zeng et al. 30 divided the special cycle slip groups into two types and used three GF linear combinations which were collected to detect and repair. Li and Melachroinos 31 proposed an enhanced repair algorithm, and the real-time implementation can be realized by Kalman filter. To decrease the influence of code observation noise, Zhao et al. 32 improved the HMW combination and determined different coefficients for satellites of different orbital types. Xiao et al. 33 applied three GF combinations and 20 was applied as the threshold of the elevation angle to ensure a high success rate. To eliminate the bias of ionosphere, prediction functions were used by Yao et al. 34 and Li et al. 35 to predict ionospheric variation. Pu et al. 36 used moving window filter function to predict ionospheric variation, and the cycle slip will be re-search when the predicted value exceeds the threshold. Xiao et al. 33 used previous DI to compensate the NL combination of current epoch when the second-order time-differenced ionospheric delay (DDI) became abnormal. The second-order time-difference algorithm was used by Liu et al. 37 to identify cycle slip during high ionospheric activity; however, the approach will increase noise of the detection combinations.
In general, there are many methods for cycle slip detection; code measurement noise and ionospheric bias are two crucial factors. Many methods assumed that code measurement noise is constant; however, code measurement noise will significantly increase when the satellite elevation angle becomes lower. 38 Based on the assumption that the ionosphere changes smoothly, the failure probability of cycle slip detection will increase. Therefore, a new method for cycle slip detection and correction is presented in this article. Firstly, in order to evaluate the code measurement noise, an elevation-based model is built. Secondly, HMW combination is used as the first detection combination due to its extremely high success rate. To minimize the combined observation noise, two optimally code-phase combinations are selected. Then, in order to decrease the impact of ionospheric bias, two moving windows of time-differenced ionospheric delay are introduced to correct the ionospheric bias.
The structure of this article is as follows. Triplefrequency linear combination model is introduced in the second section. The third section analyzes the error affecting the detection of cycle slip, and elevation-based model is built. In fourth and fifth sections, three cycle slip detection combinations and corresponding coefficients are determined. In the sixth section, two moving windows are proposed to correct the ionospheric bias of the third detection combination. In the seventh section, the approach is tested in three different scenarios. Finally, conclusions are summarized in the last section. Appendix 1 shows all symbols and their definitions used in this article.

Triple-frequency linear combination
Code and phase observations at frequency f i ði ¼ 1; 2; 3Þ can be expressed as 39 where P and F are code and phase observations, respectively. r denotes the geometric distance from satellite to receiver, and it is a frequency-independent item. c is the speed of light, dt r is the receiver clock error, and dt s is the satellite clock error. T denotes the error of the troposphere. N is the integer ambiguity and the unit is cycle. I is the ionospheric bias of frequency f 1 , and k 1i ¼ f 2 1 =f 2 i denotes ionospheric scale factor (ISF) for different frequencies. l is the wavelength; e P and e F denote the measurement noise of code and phase observations, respectively.
The triple-frequency carrier phase combination can be expressed as 21 The wavelength and integer ambiguity of the combination are

Error analysis
In the time-differenced model, the bias of troposphere, the clock error, and the hardware delay will be eliminated. However, the code measurement noise and the ionospheric bias are difficult to be eliminated by the time-differenced model, which are the main obstacle to the detection of cycle slip. To evaluate the code measurement noise and ionospheric bias, we analyzed the data of the six stations from International GNSS Service (IGS). The information of the six stations is listed in Table 1.

Satellite elevation-based measurement noise model of code observations
The measurement noise is correlated with observation environments, signal-noise ratio, and satellite elevation angle. Multipath effect is the main cause of the measurement noise. In this article, code multipath errors are used to evaluate code measurement noise. Multipath combination related to code range can be expressed as 40 where f is the frequency. i; jði; j ¼ 1; 2; 3Þ represent different frequency and i 6 ¼ j : B i mainly contains the linear combination of the phase ambiguities and the constant part of hardware delays. If no cycle slip occurs, B i can be considered a constant.
The standard deviation of multipath error corresponding to each elevation angle is selected to create multipathelevation sequence for modeling.
By analyzing and processing real data of C20-C36 at six IGS stations, we obtain the elevation-based measurement noise model of code observations where s p is standard deviation of code observations measurement noise, q is the satellite elevation angle, and c 0 , c 1 are coefficients of the model. 5 is set as the cutoff elevation angle, because the data are easily interrupted when the elevation angle is less than 5 . Table 2 lists the values of c 0 and c 1 for BDS-3 satellites. The unit of c 0 and c 1 is meter.
In this article, all of the carrier phase measurement noises s F 1 ; s F 2 , and s F 3 are assumed to be equal (s F 1 ¼ s F 2 ¼ s F 3 ¼ 0:003 m) for each satellite of BDS. B1C, B3I, and B2a are selected as three frequencies of cycle slip detection. f 1 ; f 2 , and f 3 represent the frequency of B1C, B3I, and B2a, respectively.

Time-differenced ionospheric delay
The ionospheric variation is a non-negligible factor for cycle slip detection. The ionospheric variation of two consecutive between-epoch on the B1C signal can be calculated as follows 25 where D denotes the first-order time-difference. DI is a steady change in high ionospheric activity. The second-order time-differenced ionospheric delay is where DD represents the second-order time-difference and t is the current epoch. De Lacy et al. 24 considered DI as normal distribution with mean zero, and standard deviation of DI is used to evaluate the magnitude of the ionospheric variation. Standard deviation of DI in 24 h for BDS-3 C20-C36 satellites is computed, and the result is shown in Figure 1. Figure 1 reflects the overall level of ionospheric variation within a day, and it shows that standard deviation of DI is less than 0.02 m for BDS-3 C20-C36 satellites.

First detection combination
The first detection combination is HMW combination. HMW combination is an optimal choice due to longer wavelength and higher reliability. In HMW combination, geometric distance, the bias of troposphere, and ionosphere will be eliminated. For HMW combination, the cycle slip can be accurately detected. HMW combination for BDS-3 satellites can be expressed as By differencing equation (11) between adjacent epochs, the magnitude of the cycle slip can be derived as   where DN ð0;1;À1Þ is the integer cycle slip of the first combination. When jDN ð0;1;À1Þ j > 0:5 cycle, we recognize that a cycle slip occurs. The value of cycle slip can be easily determined, and HMW combination can be repaired immediately. The combined signal DF ð0;1;À1Þ without cycle slip can be expressed as DF ð0;1;À1Þ ¼ DF ð0;1;À1Þ þ l ð0;1;À1Þ DN ð0;1;À1Þ ð13Þ Assuming DN ð0;À1;1Þ is normally distributed with zero mean and standard deviation at s DN ð0;1;À1Þ . 24 The success detection rate of HMW combination is P ð0;1;À1Þ ¼ P jDN ð0;1;À1Þ À DN ð0;1;À1Þ j < 0:5 ð14Þ From Figure 2, we can know that P ð0;1;À1Þ can achieve 100% when the elevation degree exceeds 30 . When the elevation degree is 5 , P ð0;1;À1Þ is 99.14%. It can be seen that HMW combination is still reliable at low elevation angles.

Second and third detection combinations
The traditional code-phase cycle slip detection combination is defined as DN ði;j;kÞ ¼ ðlDP 1 þ mDP 2 þ nDP 3 À DF ði;j;kÞ Þ=l ði;j;kÞ ð15Þ where l; m; n are the coefficients of the differential code observations. The disadvantage of this detection combination is that it will be seriously affected by the measurement noise of code observations. Two modified code-phase combinations are applied as the second and third detection combinations. Two virtual carrier phase observations are used. The cycle slip of the two combinations can be shown as follows DN ði;j;kÞ ¼ ðaDP þ bDP 3 þ cDF ð0;1;À1Þ À DF ði;j;kÞ Þ=l ði;j;kÞ DN ði;j;kÞ ¼ round½DN ði;j;kÞ ( ð16Þ where DP ¼ ðDP 1 þ DP 2 þ DP 3 Þ=3, a; b; c are the coefficients, and DF ð0;1;À1Þ is the phase combination whose cycle slip has been repaired in the first detection combination. By rounding DN ði;j;kÞ , an integer cycle slip DN ði;j;kÞ can be obtained.
From Figure 3, we can know that standard deviations of measurement noises DP, DP 3 , and DF ð0;1;À1Þ are smaller than those of DP 1 ,DP 2 , and DP 3 , especially at low elevation angles, which effectively reduces the influence of code observation noise on detection combinations. Since cycle slip of the first combination can always be successfully fixed, s DF ð0;1;À1Þ can be guaranteed at a normal level. The second and third detection combinations are GF combinations; hence, the coefficients a; b; c satisfy the following equation Determination of the second and third detection combinations In this section, the second and third detection combinations and the corresponding coefficients will be determined. The level of the time-differenced ionospheric delay and combination noise must be considered when choosing the optimal detection combinations. According to equations (4) and (16), the ionospheric delay bias of the second and third combinations can be derived as where b DN ði;j;kÞ is the ISF of the second and third detection combinations.
We can obtain the noise of the second and third detection combinations as The missed detection rate P M represents the probability that the detection combination does not detect the existence of a cycle slip when cycle slip occurs, it can be indicated as 32 We analyzed the relationship among d DN ði;j;kÞ , s DN ði;j;kÞ , and P S or P M , and the results can be seen in Figure 4. Figure 4(a) shows the relationship among d DN ði;j;kÞ , s DN ði;j;kÞ , and P S . The yellow surface represents equation (20); two red lines are boundaries between the surface and the plane P M ¼ 99%. The absolute value of slopes of the two red lines jkj ¼ jk 0 j ¼ 1:6753. If 99% is adopted as the acceptable degree of success rate, the following condition must be satisfied: jd DN ði;j;kÞ j 1:6753s DN ði;j;kÞ ð22Þ Figure 4(b) shows the relationship among d DN ði;j;kÞ , s DN ði;j;kÞ , andP M . It is obvious that the missed detection rate is only related to the noise of detection combination. If the noise of detection noise is too large, the cycle slip is likely to be overwhelmed and cannot be detected. The blue vertical line represents that the missed detection rate achieves 1% when s DN ði;j;kÞ ¼ 0:1581. To ensure P M 1%, the following condition must be satisfied: Selection of optimal detection combinations According to the above analysis, 99% success detection rate and 1% missed detection rate are adopted as the selection criterion of the second and third detection combinations. Hence, the second and third detection combinations must satisfy equations (17), (22), and (23). In addition to these three conditions, the three detection combinations must be linearly independent. To select the optimal virtual signals F ði;j;kÞ and the corresponding coefficients a; b; c, the values i; j; k of F ði;j;kÞ are tested within the range of [À20, 20], as shown in Figure 5. Here, DI is configured as 0.02 m. As shown in Figure 5, each group of i; j; k will be tested for satellite elevation angle q ¼ 5 . If i; j; k and a; b; c can satisfy the lowest elevation angle, other elevation angles can also be satisfied. The minimum of s DN ði;j;kÞ will be calculated. Only when the two conditions (22) and (23) are satisfied at the same time, the values of i; j; k and a; b; c can be output. Otherwise, the next group of i; j; k will be tested.
When the above process is completed, all possible combinations will be output. Notice that the three detection combinations must be linearly independent, so the signals (1, À4, 3) and (À4, 6, À1) are selected as the second and third detection combinations, respectively, and the specific results are shown in Table 3. Figure 6 shows P S and P M of the second and third detection combinations.
According to Figure 6, it is known that the success detection rate of the second detection combination (1, À4, 3) is close to 100%, and the missed detection rate is close to 0%. The second detection combination benefits from longer wavelength and smaller ISF, which are minimally affected by the code measurement noise and ionosphere. Due to shorter wavelength and larger ISF, the combined signal (À4, 6, À1) is more vulnerable to the ionospheric bias. When DI ¼ 0:03 m, P S of (À4, 6, À1) will drop below 94% for all possible elevation angles. Since the missed detection rate is only related to s DN ði;j;kÞ , it can maintain a low level regardless of the ionospheric bias. To eliminate the ionospheric bias of (À4, 6, À1) and improve the success detection rate, the ionospheric bias must be compensated. The detection value of the detection combination (À4, 6, À1) corrected for ionospheric bias can be expressed as     slips DN 1 ¼ DN 2 ¼ DN 3 occur on the three carrier phase observations, the value of DDI will become abnormal.

First moving window
The first moving window of DDI is used to check whether DDI has a significant increase. When there is no cycle slip, we assume that DDI is normally distributed with zero mean. W Fs DDI is the standard deviation of DDI in the moving window. Cycle slips including DN 1 ¼ DN 2 ¼ DN 3 may occur when DDI cur satisfies the following condition: where DDI cur is DDI of current epoch, and 4 Á W Fs DDI is adopted as the threshold. The number of moving windows is set to 20. If DDI cur satisfies equation (26) or DN ð0;1;À1Þ 6 ¼ 0 or DN ð1;À4;3Þ 6 ¼ 0, equation (24) cannot be updated by DI of current epoch, and the second moving window is adopted to evaluate DI of current epoch.

Second moving window
To evaluate DI of the current epoch, the second moving window of DI is adopted. The average value of DI in the second moving window is used to evaluate DI of current epoch. The mean of DI in the moving window can be expressed as where W S DI is the average value of DI in the moving window, and the number of moving window is n ¼ 10.
When jDDI cur j > 4 Á W Fs DDI or DN ð0;1;À1Þ 6 ¼ 0 or DN ð1;À4;3Þ 6 ¼ 0, W S DI is used as DI of current epoch to update equation (24). Figure 7 shows that the process of the third combination compensates for the ionospheric bias. When DN ð0;1;À1Þ 6 ¼ 0 or DN ð1;4;À3Þ 6 ¼ 0, we can be sure that there must be a cycle slip due to success detection rates of the first and second detection combinations which are close to 100%. However, the first and second detection combinations are insensitive to cycle slip combination DN 1 ¼ DN 2 ¼ DN 3 ; hence, the first moving window is introduced to detect identical cycle slips. When jDDI cur j > 4 Á W Fs DDI , cycle slips including identical cycle slips may occur, and equation (24) will be updated by W S DI . If all three conditions in Figure 7 are not satisfied, equation (24) will be updated by DI cur . Figure 8 shows DI and DDI of C29 from MIZU station. From Figure 8(a), we can see that the threshold DDI þ 4W Fs DDI and DDI À 4W Fs DDI of the first moving window is constantly changing with DDI, and a cycle slip may occur when DDI cur exceeds the threshold. Figure 8(b) shows the evaluation value of the second moving window for DI of the current epoch. It can be seen that DI and W S DI have the same trend. Except few epochs, the difference between DI and W S DI is small. From Figure 8(a) and (b), we can know that as the elevation angle  of satellite decreases, jDIj and jDDIj will increase, and the corresponding threshold of the first moving window will also increase.

Numerical tests and analysis
To verify the effectiveness of our proposed approach, we test it in three different scenarios.

Test with simulated noise of code observations
In a poor observation environment, the noise of code observations will increase significantly. Due to the existence of code observation in the detection combination, the noise of code observations is bound to affect the detection of cycle slip. To test the ability of the proposed method to detect cycle slip in the case of high code observation noise, the noise is added to the original code observations. Gaussian  noise with zero mean is added to the three code observations in the second experiment. We simulated two scenarios, in which the standard deviation of the added noise is 0.5 and 0.8 m, respectively. The data of BDS-3 satellites from SGOC, ULAB, and STUM stations on June 19, 2019 are used and the data are clean without cycle slip. The sampling interval is 30 s. At every 10 epochs, cycle slip groups from (0, 0, 0) to (2, 2, 2) will be added to the data in turn. The results of SGOC station are shown in Figure 9. Table 4 shows the total number of simulated cycle slips and the number of incorrect detection and repair.
In Table 4, scenario #1 is normal and no noise is added. Simulated noise is added to scenario #2 with standard deviation of 0.5 m. In scenario #3, simulated noise is added using noise with standard deviation of 0.8 m. Scenario #3 has higher noise level than scenario #2.
From Figure 9 and Table 4, we can know that the number of failures increases as the noise level increases, but the success rate is still greater than 95%. From Table 3, it is known that the code observations have a very small weight in the detection combination; hence, the proposed approach is less affected by code observations.

Test under normal circumstances
The second experiment was carried out under normal circumstances, the noise of code observations is at normal level, and the ionosphere changes slowly. We compare the proposed approach with the approach proposed by Zhao  et al. 32 They used three code-phase combinations to detect and correct cycle slip, and in order to eliminate the ionospheric bias in the third combination, they propose DÑ ði;j;kÞ ¼ DN ði;j;kÞ À k Á b DN ði;j;kÞ DI=l ði;j;kÞ ð28Þ where kð0 k 1Þ is scale factor to balance the corrected percentage of the ionospheric bias and the amplification of the measurement noise. When k ¼ 1, the ionospheric bias is totally corrected. When k ¼ 0, the ionospheric bias is not corrected. When 0 k 1, the ionospheric bias is partly corrected. The real BDS-3 data from SGOC station on June 19, 2019 are used to test the effectiveness of these two methods. The original observations have been repaired and have no cycle slips beforehand. Cycle slip groups from (0, 0, 0) to (2,2,2) are added to the data in turn, at every 10 epochs. The sampling interval is 30 s. The approach of Zhao et al. 32 is tested when k ¼ 0; 0:5; 1, respectively. The results of C20-C36 satellites can be seen in Figure 10, in which C31 satellite of SGOC station has no data.
In Figure 10, the graph consisting of two red triangles and a green square represents failure of detection and repair. The number in the right coordinate area represents the times of failed detection and correction. As shown in Figure 10, the times of failed detection and correction in Figure 9(a) have a significant reduction compared with   . Therefore, the method proposed by this article has stronger ability to eliminate the effect of ionospheric bias. From Figure 10, it can be known that all of the failed detection and repair of cycle slip occurs at low elevation angle. At low elevation angles, the ionospheric bias and code measurement noise may increase. If the ionospheric bias cannot be handled correctly, it will lead to erroneous results.

Test under high ionospheric activity
In theory, as long as the prediction of DI under high ionospheric activity can reach a required precision, the cycle slip will be detected. Many methods use STPIR algorithm to perform the second-order time-difference on detection combination, which will increase the noise of detection combination. To verify the validity of the approach, the data of BDS-3 during the period of active ionosphere will be used. The observations of BDS-3 are collected from WUH2 station on March 21, 2019, and the ionosphere is active for a period of time. The sampling interval is 30 s. Figure 11 shows DI and DDI of C28. In Figure 11(a), the ionospheric delay variation of C28 can reach 0.18 m. The prediction values W S D I for ionospheric delay variation are shown in Figure 11(b). Different types of cycle slips are added to original observations, and the results are shown in Figure 12 and Table 5. Small, particular, and large cycle slips are simulated in the third experiment. In Figure 12, b DN ðÀ4;6;À1Þ DI=l ðÀ4;6;À1Þ denotes the deviation of the cycle slip detection caused by the ionospheric delay variation, and the unit is cycle. From Figure 12, we can see that jb DN ðÀ4;6;À1Þ DI=l ðÀ4;6;À1Þ j is more than 0.5 cycles in many epochs. The bias caused by ionospheric delay can reach two cycles in 80th and 90th epochs. Although the ionosphere has a significant impact on the third detection combination, all simulated cycle slips can be correctly corrected.

Conclusions
A new method of cycle slip detection and repair for BDS-3 triple-frequency observations is proposed. Firstly, two modified code-phase combinations are used to minimize the effect of noise of code measurement. Then, 1% missed detection rate and 99% success detection rate are used as the criterion to optimally select the second and third detection combinations by combining BDS elevationbased model and HMW combination. Since the third detection combination is easily biased by ionosphere, two moving windows are used to correct the ionospheric bias of the third combination. The method has been tested in three different scenarios. In the first scenario, simulated Gaussian noise is added to the three code observations and the proposed method is still effective in the environment of high code observation noise. Compared with the method proposed by Zhao et al., 32 the proposed method has obvious advantages in reducing the times of failed detection and correction. From the first and second experiments, we can know that most of the detection failures occur at low elevation angles. We can increase the cutoff elevation angle of satellites to reduce the times of failed detection and repair. Further, BDS-3 data under high ionospheric activity is used to verify the validity of the algorithm, cycle slips of different types are added to the original carrier phase observations, and all cycle slips can be correctly detected and fixed. work was supported by the National Natural Science Foundation of China (no. 41574010).