Next Article in Journal
A Rigorous and Integrated On-Water Monitoring System for Performance and Technique Improvement in Rowing
Previous Article in Journal
Smart Contract Broker: Improving Smart Contract Reusability in a Blockchain Environment
Previous Article in Special Issue
Image Mapping Accuracy Evaluation Using UAV with Standalone, Differential (RTK), and PPP GNSS Positioning Techniques in an Abandoned Mine Site
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Extraction and Characterization of Pseudorange Multipath Based on BDS-3 Multi-Frequency Observations

1
Key Laboratory of Aviation-Aerospace-Ground Cooperative Monitoring and Early Warning of Coal Mining-Induced Disasters of Anhui Higher Education Institutes, Anhui University of Science and Technology, Huainan 232001, China
2
School of Environment and Surveying Engineering, Suzhou University, Suzhou 234000, China
3
School of Earth and Environment, Anhui University of Science and Technology, Huainan 232001, China
4
School of Geomatics, Anhui University of Science and Technology, Huainan 232001, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(13), 6151; https://doi.org/10.3390/s23136151
Submission received: 12 June 2023 / Revised: 27 June 2023 / Accepted: 3 July 2023 / Published: 4 July 2023
(This article belongs to the Special Issue Methods and Applications of Multi-GNSS PNT and Remote Sensing)

Abstract

:
Global Navigation Satellite System (GNSS) observations are subject to various errors during their propagation process. A reasonable correction of these errors can improve the positioning, navigation, and timing (PNT) service capability. The impact of multipaths on pseudorange observations can reach a decimeters or even meters level. However, their mechanism is complex and there is currently no universally accepted high-precision correction model. The correlation between the pseudorange multipaths (MP) of BDS-2 satellites and satellite elevation has been confirmed, while there have been fewer analyses of the MP characteristics for different frequencies of BDS-3 satellites. The broadcasting of multi-frequency observations in BDS-3 should theoretically make the extracted MP more accurate compared to traditional methods. Based on this, in this contribution, a multi-frequency MP extraction algorithm based on the least squares principle is proposed, which can simultaneously eliminate the influence of higher-order ionospheric delay. The analytical expression for only eliminating first-order ionospheric delay is successfully derived. Subsequently, the characteristics of the MPs extracted from different frequency combinations and the impact of combination noise on the extraction accuracy are discussed. The influence of second-order ionospheric delay on the MPs for each frequency under different combination noises, as well as the periodic behavior exhibited in long-term observations of the BDS-3 medium earth orbit (MEO) and inclined geosynchronous orbit (IGSO) satellites, are also analyzed. Finally, the correlations between the MPs of each frequency of BDS satellite and elevation are quantitatively analyzed based on observations from 35 stations. Overall, this work has positive implications for the study of the MP characteristics of BDS-3 and subsequent modeling efforts.

1. Introduction

GNSS has been widely used in PNT services due to its advantages of high accuracy, high frequency, and all-weather capabilities [1,2]. As the main member of navigation systems, the Beidou Satellite Navigation System (BDS) was independently designed and developed by China and has undergone a “three-step” development strategy. On 31 July 2020, it was officially opened for global service. Compared with other navigation systems, the BDS satellite constellation is composed of three types of orbiting satellites: geostationary orbit (GEO), IGSO, and MEO. Both the IGSO and MEO satellites of BDS-3 can transmit observations of five frequency signals [3,4]. Despite the increased complexity in processing the observations of multi-orbit and multi-frequency signals, this ability significantly improves the overall service performance and accuracy of the system [5,6].
Pseudorange and carrier phase observations are the basic measurements used for positioning in most navigation systems. These measurements are subject to various errors during signal propagation, such as atmospheric delay, ionospheric delay, and multipath effects [7,8,9,10]. Reasonably eliminating or reducing the effects of these errors can effectively improve the positioning accuracy [11,12,13,14]. Previous studies have shown that a reasonable combination of multiple frequency observations can reduce the impact of these errors, such as ionosphere-free combinations, ionosphere-reduced combinations, and geometric-free combinations [15,16,17,18,19,20]. With the increasing number of transmission signals, there are more favorable combinations available.
Due to their complex mechanism and the difficulty in modeling them, it is difficult to accurately correct multipath delay errors during the positioning process. Currently, this issue has been researched by many scholars [21,22,23,24,25]. However, some studies in recent years have shown that BDS is not only affected by the reflected signals from the surrounding environment on the ground, but also contains certain systematic biases in the pseudorange observations on its satellite end. In 2012, Hauschild et al. first discovered the systematic errors of the MPs in BDS-2 MEO satellite signals [26]. Subsequently, relevant scholars analyzed the MPs of three BDS-2 types of satellites and explored model correction algorithms [27,28,29,30]. The results showed that correcting the MP can effectively improve the positioning accuracy. With the gradual launch and deployment of BDS-3 satellites, some scholars have also investigated the signal qualities of various frequencies of BDS-3 and found that the MPs of BDS-3 satellites are significantly smaller than those of BDS-2 [31,32,33]. However, it is worth noting that, when extracting the MP of each BDS-3 satellite frequency, the traditional dual-frequency pseudorange and carrier phase observations combination method is still used and the multi-frequency signals provided by BDS-3 are not fully utilized. In theory, introducing additional carrier phase observations can better extract the true signal of an MP [20,34]. Therefore, how to reasonably use multi-frequency observations to extract the MPs for each frequency should be studied.
The traditional MP extraction method can only obtain a unique solution by eliminating the geometric distance term independent of the frequency and the first-order ionospheric delay term, as it only uses dual-frequency observations. However, when triple-, quad-, or five-frequency observations are used, redundant observations can be obtained, which leads to the optimal solution with the minimum combination noise. Moreover, the introduction of triple-, quad-, and five-frequency observations can make it possible to simultaneously eliminate higher-order ionospheric delay terms. When the first- and second-order ionospheric delay terms are simultaneously eliminated, a unique solution can be obtained using triple-frequency observations, while quad- and five-frequency observations can still achieve optimal solutions. The use of multi-frequency observations can also provide support for analyzing the influence of higher-order ionospheric delay terms on MP extraction.
The method of extracting MPs using BDS-3 multi-frequency observations is mainly discussed in this work, and based on this, the characteristics of the MPs on BDS-3 MEO and IGSO satellites are analyzed. The main structure of the article is arranged as follows. The formula for calculating the optimal combination coefficient of an MP based on multi-frequency observations is derived in Section 2, where a combination coefficient calculation method based on the least squares algorithm is proposed. In Section 3, the characteristics of MPs at various frequencies are studied considering three aspects: the impact of different frequency combinations on MP extraction, the influence of second-order ionospheric delay on MP extraction, and the correlation between MPs and elevation. Finally, some conclusions and the next research plans are given in Section 4.

2. Methods

Carrier phase and pseudorange observations are the basic observations broadcasted by each navigation system and their linearized observation equations can be expressed as [35]:
{ Φ i = ρ r s + c d t r c d t s + T r k i I 1 , 1 ν i I 1 , 2 + λ i N i + ε Φ i P i = ρ r s + c d t r c d t s + T r + k i I 1 , 1 + 2 ν i I 1 , 2 + MP i + ε P i
where the subscript i (i = 1, 2,…, n) represents different frequencies; Φ and P represent the carrier phase and pseudorange observations; ρ r s denotes the geometric distance between the satellite and the receiver; d t r and d t s represent the receiver and satellite clock errors, c is the speed of light; Tr represents the tropospheric delay; k i = f 1 2 / f i 2 and ν i = f 1 3 / f i 3 represent the first- and second-order ionospheric delay amplification factors of frequency i, respectively, I1,1 and I1,2 denote the first- and second-order ionospheric delay of the f1 frequency, λ i and N i represent the wavelength and ambiguity of frequency i, MPi is the pseudorange multipath, and ε Φ i and ε P i denote other errors of the carrier phase and pseudorange observations.
Usually, MPi can be extracted from dual-frequency observations and its expression is [26,31,36]:
MP i = P i f i 2 + f j 2 f i 2 f j 2 Φ i + 2 f j 2 f i 2 f j 2 Φ j B i , j + ξ
where i and j denote different frequencies, B i , j is the combination of the ambiguity and hardware delay, and ξ represents other errors. When the carrier phase observation does not include cycle slip, B i , j can be treated as a constant, that is, after removing B i , j from the sequence extracted from Equation (2), an MP sequence with certain noise can be obtained [37].
The magnitude of the noise can have a certain impact on the modeling and correction of an MP. As shown in Equation (2), the extraction of MPs at different frequencies only uses the carrier phase observations of two frequencies. However, with the increase in the frequencies broadcasted by various navigation systems, introducing additional frequency carrier phase observations should be able to obtain MP sequences with less noise. For BDS-3 MEO and IGSO satellites, five signal frequencies can now be broadcasted. If five-frequency carrier phase observations are used to extract an MP at a certain frequency, it can be represented as:
MP j = η 0 P j + η j Φ j + i = 1 , i k 5 η i Φ i
where η 0 and η j represent the combination coefficients of the pseudorange and carrier phase observations at frequency j, respectively, while η i denotes the combination coefficients of the carrier phase observations at other frequencies.
From Equation (1), it follows that, for the MP1 extraction using Equation (3), when only eliminating the influence of I1,1, η 0 ~ η 5 should satisfy:
{ η 0 = 1 η 1 + η 2 + η 3 + η 4 + η 5 = 1 η 0 k 1 η 1 k 1 η 2 k 2 + η 3 k 3 η 4 k 4 η 5 k 5 = 0 η 0 2 + η 1 2 + η 2 2 + η 3 2 + η 4 2 + η 5 2 = min
When the effects of I1,1 and I1,2 are considered simultaneously, Equation (4) can be combined with Equation (5).
2 η 0 ν 1 η 1 ν 1 η 2 ν 2 η 3 ν 3 η 4 ν 4 η 5 ν 5 = 0
According to Equations (4) and (5), η 1 ~ η 5 can be obtained by the Lagrange multiplier method. In fact, the solution of Equation (4) should introduce additional parameters, which increases the computational workload. In this work, starting from the observation Equations of the pseudorange and carrier phase, a method for calculating the multi-frequency MP combination coefficients based on the least squares principle is proposed, directly according to the characteristics of the combined parameters.
Equation (4) shows that the coefficients of the MP combination need to satisfy the condition that the combined geometric distance term and ionospheric delay term are both equal to zero. Therefore, when considering only the geometric distance term, I1,1, and MP, the observation equations for the pseudorange and carrier phase can be simplified as follows:
[ P 1 Φ 1 Φ 2 Φ 3 Φ 4 Φ 5 ] Z = [ 1 k 1 1 1 k 1 0 1 k 2 0 1 k 3 0 1 k 4 0 1 k 5 0 ] H [ ρ ˜ I 1 , 1 MP 1 ] X + [ ε P 1 ε Φ 1 ε Φ 2 ε Φ 3 ε Φ 4 ε Φ 5 ] ε
where ρ ˜ denotes the sum of the geometric distances of the all frequency-independent terms from the satellite to the receiver.
Assuming that the covariance matrix of the observation is Q, the least squares solution of the estimated parameters X can be obtained from Equation (6), which is expressed as:
X = ( H T Q 1 H ) 1 H T Q 1 Z
Equation (7) can also be seen as a combination of the pseudorange observations and carrier phase observations at different frequencies to extract ρ ˜ , I1,1, and MP1, where the combination coefficients are:
x = ( H T Q 1 H ) 1 H T Q 1
In order to satisfy the condition of Equation (4), the transformation matrix R = [ R 1 R 2 R 3 ] T can be defined, where: R 1 = R 2 = [ 0 0 0 ] , R 3 = [ 0 0 1 ] . At this point, the combined ρ ˜ and I1,1 can be kept at 0, while the coefficient of MP1 remains at 1, that is:
x ^ = R x = R ( H T Q 1 H ) 1 H T Q 1
According to Equation (9), the combination coefficient of MP1 can be expressed as:
η = [ η 0 η 1 η 2 η 3 η 4 η 5 ] = R 3 ( H T Q 1 H ) 1 H T Q 1
The optimal coefficient of the MP combination for a certain frequency can be obtained through Equation (10). Assuming that the carrier phase observations of n frequencies are used, the MP combination coefficient can be expressed as:
{ η 0 = 1 η 1 = ( i = 2 n k i 2 n + 1 ) / C n η e = ( i = 2 n ( k i 2 k i k e + k i ) ( n + 1 ) k e + 2 ) / C n ( e = 2 , 3 , , n )
C n = i = 2 n 1 j = i + 1 n ( k i k j ) 2 + k = 2 n ( k k 1 ) 2
Equations (11) and (12) provide analytical expressions for the optimal coefficients of multi-frequency MP combinations. When extracting MPs using different frequency observations, it is only necessary to directly bring in the ki corresponding to different frequencies. However, it is important to note that ki is not fixed and constant. Although the same carrier phase observations are used, when the frequency of extracting the MP is different, ki also needs to be updated again.
If the influence of I1,2 needs to be eliminated when extracting an MP, Equation (6) can be rephrased as:
[ P 1 Φ 1 Φ 2 Φ 3 Φ 4 Φ 5 ] Z = [ 1 k 1 2 ν 1 1 1 k 1 ν 1 0 1 k 2 ν 2 0 1 k 3 ν 3 0 1 k 4 ν 4 0 1 k 5 ν 5 0 ] H [ ρ ˜ I 1 , 1 I 1 , 2 MP 1 ] X + [ ε P 1 ε Φ 1 ε Φ 2 ε Φ 3 ε Φ 4 ε Φ 5 ] ε
By updating the transformation matrix R( R 4 = [ 0 0 0 1 ] ) accordingly, the MP combination coefficients eliminating I1,2 can be quickly obtained. Taking triple-frequency observations as an example, the combination coefficient when simultaneously eliminating I1,1 and I1,2 can be expressed as:
{ η 0 = 1 η 1 = ( 2 ( k 2 k 3 ) + ( ν 2 ν 3 ) ( k 2 ν 3 k 3 ν 2 ) ) / D 3 η 2 = ( 3 k 3 + 2 ν 3 + 1 ) / D 3 η 3 = ( 3 k 2 2 ν 2 1 ) / D 3
D 3 = ( k 2 k 3 ) + ( ν 2 ν 3 ) + ( k 2 ν 3 k 3 ν 2 )
Taking the B1I frequency as an example, Table 1 shows the combination coefficients and carrier phase observation noise amplification factors ( Ω = η 1 2 + η 2 2 + + η n 2 ) for extracting MPs using different frequency combinations. The DF, TF, QF, and FF represent dual-frequency, triple-frequency, quad-frequency, and five-frequency, respectively. Figure 1 shows the combination noise of the B1C/B1I/B3I/B2a frequency when extracting MPs through different frequency combinations.
From Table 1 and Figure 1, it can be seen that the Ω varies significantly under different frequency combinations. For example, the Ω of B1I/B1C can be 25–35 times that of other dual-frequency combinations and about 45 times that of quad- and five-frequency combinations. Related studies have pointed out that, when the Ω is large, it will amplify the observation error implied by the carrier phase observations [20,34]. Therefore, combinations of Ω similar to B1I/B1C should be careful not to be used for extracting MPs. The commonly used B1I/B3I combination has a combined noise of 6.245. From Figure 1, it can be seen that the Ω is greater than most combinations. In theory, when the Ω is smaller, the extracted MP sequence contains less noise. Therefore, combinations such as triple-frequency or quad-frequency should be able to extract MP sequences with smaller error fluctuations.
Additionally, when both the I1,1 and I1,2 are eliminated simultaneously, the Ω significantly increases and the optimal combination is different from when only I1,1 is considered. For example, the combination of B1I/B1C/B2a for B1I is relatively optimal for the triple-frequency combination. However, when both I1,1 and I1,2 are simultaneously considered, this combination actually has the largest combination noise. Through the analysis of the coefficients of the four frequency combinations, it is found that the contribution of the B3I frequency to the combination is relatively small when considering only I1,1, while it has the largest contribution when eliminating both I1,1 and I1,2 simultaneously. Another point to note is that, when extracting MPs with the same frequency of different frequencies, the combination coefficients are different. Although the expression of the coefficients is the same, due to the change in the first frequency, the corresponding k2 and k3 will also change. For instance, when using the combination of B1C/B1I/B2a to extract the MP of each frequency, the corresponding combination noises are Ω B 1 C = 3.593 , Ω B 1 I = 3.621 , and Ω B 2 a = 4.820 , respectively.

3. Experiments and Analyses

In this section, the correctness of the theories proposed in this work and the characteristics of the MPs for each frequency of BDS-3 will be discussed through experiments. The observations of 35 stations (18 MGEX stations and 17 iGMAS stations) over 7 days (DOY 1–7, 2021) were selected, and all stations could receive observations of the B1I, B1C, B3I, and B2a frequencies of BDS-3. The sampling interval for the observations at each station was 30 s. Figure 2 shows the distribution of these stations. Additionally, a set of dynamic vehicle observations (CAR1 station) was also selected to further validate the theories proposed in this work. The data acquisition scheme can be found in reference [38]. Standard deviation (STD), mean, and range were selected as statistics to measure the accuracy.

3.1. Comparison of Characteristics of MP Extracted with Different Frequency Combinations

In order to compare the differences in the extracted MP sequences when different frequency combinations were used, taking the B1I and B1C frequencies as examples, four modes with varying combination noise from dual-frequency to quad-frequency were selected to extract the MPs for each frequency. The specific combination frequencies and coefficients are shown in Table 2. Although there were differences in the combined noise of the four modes, each mode listed was the optimal combination coefficient calculated for the corresponding frequency. This section mainly verifies the difference between the different frequency combinations on the obtained MP sequences, without considering the difference between the optimal combination and non-optimal combination among each frequency. In fact, the experimental results can still reflect the difference between them, as the combination noise of the non-optimal combinations also increased correspondingly, and the comparison results should be similar to the results.
Taking the observation of the KUN1 station as an example, Figure 3 presents the MP sequences extracted from the C28 and C40 satellites using different frequencies. In each sub-figure, the upper portion represents the extracted MP sequences for the B1I and B1C frequencies, while the lower portion shows the first-order difference between each mode and the QF mode. Table 3 provides the accuracy statistics for the MP sequences extracted using different combination modes, as well as the improvement relative to DF1.
By comparing the four modes, it can be observed that the MP sequences extracted using the triple-frequency and quad-frequency combinations exhibited a similar trend to those extracted using the traditional dual-frequency method, but with smaller fluctuations. This indicates that the coefficient calculated using Equation (11) was correct and optimal. The results also show that if the Ω was small, the fluctuation of the extracted MP sequence was smaller. By comparing this with the QF mode, it can be observed that, when the difference of the Ω was relatively small, the first-order difference between the two sequences generally fluctuated within 3 cm. However, the difference between the QF mode and the DF1 mode still fluctuated at the meter level. This indicates that, when the Ω was small, the extracted sequences were essentially consistent with each other. Additionally, Table 3 also shows that the improvement in the accuracy, relative to the DF1 mode, was consistent across the different combination modes. Therefore, when certain frequency observations were missing, other combination modes could be adaptively selected to extract the MP. Figure 4 presents histograms of the frequency distributions for the MP sequences extracted using different modes. It can be seen that the sequences from the DF2, TF, and QF modes are relatively more concentrated around zero. However, due to the favorable observation environment at the KUN1 station, the representation in the histograms of the frequency distribution is not significant.
To further validate the theories proposed in this article and compare the characteristics of the MP sequences extracted from the different modes, the same frequency combinations and coefficients shown in Table 2 were used to process the dynamic observations. Similar to Figure 3, Figure 5 presents the results and difference sequences for satellites C23 and C38, while Figure 6 displays histograms of the frequency distribution for the MP sequences extracted using the different modes. Table 4 provides the accuracy statistics for the two satellites.
It can be seen that, due to the significant influence of the surrounding environment during the collection of the dynamic observations, the selected multi-frequency combinations in this work could effectively reduce the fluctuation level and magnitude of the MP sequences. Figure 5 provides a more intuitive representation of the numerical variation range of the MPs extracted using different combinations compared to the DF1 mode. Although the improvement in the accuracy achieved using the DF2, TF, and QF modes for the dynamic observations was not significantly different from that of the DF1 mode, this accuracy improvement was more pronounced compared to the KUN1 station. The B1I and B1C frequencies of the two satellites at the KUN1 station showed accuracy improvements ranging from 38–42% to 13–15%, respectively, while the CAR1 station exhibited even higher improvements, reaching over 30–39% for the B1I frequency and 43–50% for the B1C frequency. The comparison between the two stations highlights the increased importance of using optimal multi-frequency combinations in dynamic experiments.

3.2. Analysis of the Influence of Second-Order Ionospheric Delay

In Section 3.1, the influence of the four different frequency combination modes on the MP extraction was compared using the B1I and B1C frequencies as examples. However, the coefficients used were the optimal solutions only when I1,1 was eliminated, and the influence of the simultaneous elimination of I1,2 on the extracted MPs was not considered. From the derivation process of the optimal coefficients in Section 2, it can be found that, when only I1,1 was considered, the dual-frequency combination could only have a unique solution, while the triple-frequency, quad-frequency, and five-frequency combinations could obtain the optimal solution with the minimum combination noise due to the presence of redundant observations. Similarly, when considering both I1,1 and I1,2, Equations (4) and (5) need at least three frequency observations to achieve a unique solution, and the quad-frequency and five-frequency combinations could obtain the optimal solution. Based on the findings in Section 3.1, when the combination noise was similar, the extracted MP sequences were also very similar. In this section, a further analysis is conducted to examine the impact of ionospheric delay on these MP sequences. Taking the triple-frequency and quad-frequency combinations as examples, the influence of second-order ionospheric delay on the MP sequences extracted using different frequency combinations is analyzed. Table 5 provides the selected combination modes and their corresponding coefficients for the B1C, B3I, and B2a frequencies. The combination modes for the B1I frequency can be found in Table 1, where B1I/B1C/B3I, B1I/B1C/B2a, B1I/B3I/B2a, and B1I/B1C/B3I/B2a are represented as B1ITF1, B1ITF2, B1ITF3, and B1IQF1, respectively.
Theoretically, regardless of whether the I1,2 was eliminated or not, the coefficients of the pseudorange observations involved in the combination were all 1, which meant that the truth values contained in the extracted sequences were consistent. However, due to the different combination coefficients of the two methods, the errors contained in the extracted sequence also differed. In order to analyze the impact of eliminating I1,2 on the extraction of the MP sequences, the optimal coefficients of the same mode in Table 5 were used to extract the sequences for eliminating I1,1 and I1,1/I1,2, and a difference was made between the two sets of sequences. The difference between the two sets of sequences could be considered as the influence of I1,2 under that mode. Taking the continuous 7-day observations from the WUH2 station as an example, Figure 7 and Figure 8, respectively, show the difference sequences for C28 (MEO satellite) and C38 (IGSO satellite) at various frequencies. B1CTF1 in the figures represents the first-order difference sequence for the B1CTF1 mode, and the meanings of the other sequences are similar. Table 6 provides statistical information on the accuracy of all the BDS-3 MEO and IGSO satellites at the experimental stations for the 7 days under each mode.
Based on Figure 7 and Figure 8, as well as Table 6, it can be observed that:
(1) There was a clear correlation between the magnitude of the fluctuations in the difference sequences and the Δ Ω . The Δ Ω for the B1ITF1, B1ITF2, B1CTF1, B1CTF2, B3ITF1, and B2aTF1 modes were all above 140. The STD and range values for each sequence were also greater than 0.24 m and 1.5 m, respectively, which were 3–5 times larger than those for the other modes. Among them, the B2aTF1 mode had the largest Δ Ω , which was 245.626. Its STD and range were 0.401 m and 2.478 m, respectively, which were also significantly higher than those of the other modes. With the exception of the MP extraction mode for the B2a frequency, the fluctuation range of the differenced sequences for other modes was less than 0.5 m and the STD was basically around 0.05–0.07 m.
(2) According to reference [32,33], when only I1,1 is considered, the anti-multipath performance of the B2a frequency is better than that of other BDS-3 frequencies. However, the experiments showed that, when eliminating I1,2 simultaneously, the MPs of the B2a frequency extracted using different modes exhibited larger fluctuations. Table 6 also indicates that, except for modes such as B2aTF1 with significant noise, the STDs for the other three modes were approximately 2–3 times larger than those for the other frequency modes. From the statistical results, it can be seen that the averages of the STDs and ranges of the B2aTF2, B2aTF3, and B2aQF1 modes were 0.113 m and 0.726 m, respectively, which were about 1.8, 1.7, and 2.2 times and 1.6, 1.7, and 1.9 times greater than those of the B1I, B1C, and B3I frequency combinations. The anti-multipath performance of the B3I frequency was also relatively superior, and the combination with a smaller Δ Ω corresponded to a lower STD and range compared to those of the other frequencies. This characteristic was contrary to that of the B2a frequency, and further research is needed in the future.
(3) By comparing the differenced MP sequences of the C28 and C38 satellites at various frequencies, it could be found that, regardless of the magnitude of the fluctuations in the differenced sequences, the differenced sequences of the C38 satellite at each frequency exhibited a clear periodicity. However, the periodicity of the C28 satellite was not prominent. Compared with other IGSO and MEO satellites, they all showed the same periodicity pattern. Additionally, when comparing the four frequencies, it could be observed that their periodic patterns were consistent.
In this section, the analysis of the MP difference sequences under the different modes was also conducted using the dynamic observations from Section 3.1. Figure 9 and Figure 10 provide the MP difference sequences extracted with the consideration of I1,2 for various frequencies of the C23 and C38 satellites at the CAR1 station. Table 7 presents the accuracy statistical information of all the MEO and IGSO satellites at the CAR1 station using different modes.
According to the experimental results, it can be concluded that the fluctuation magnitude of the differenced MP sequences of the BDS-3 IGSO satellites at the CAR1 station was smaller than that of the MEO satellites. The STD and range for each frequency were approximately half of those for the MEO satellites. In addition, similar to the WUH2 station, the fluctuation amplitude of each sequence was positively correlated with the Δ Ω . The fluctuation magnitude of the differenced sequences for the B2a frequency in various modes remained higher than that of the B1I, B1C, and B3I frequencies and the STD and range were also around 1.5–2.0 times larger. The accuracy of the B3I frequency was also better than that of the other frequencies. Due to the observation time being only 2 h, all the sequences fluctuated around the zero value, and the mean values were in the millimeter or sub-millimeter level. The IGSO satellites did not exhibit periodicity. Compared with the WUH2 station, although the observation at the CAR1 station was more affected by the environment, overall, the magnitude of the second-order ionospheric delay remained consistent. For modes with a high Δ Ω , the STD and range of the CAR1 station were slightly higher than those of the WUH2 station. However, when the Δ Ω was small, the STD and range were actually better than that of the WUH2 station.

3.3. Analysis of the Correlation between MP and Elevation of BDS-3 Satellite

Based on the previous experiments, it is known that the impact of the BDS-3 satellite’s MP on the pseudorange observations can reach the order of decimeters or even meters, and it needs to be corrected during precise positioning. Currently, there have been numerous studies discussing MP correction methods for BDS-2 satellites, but relatively fewer studies have been conducted for BDS-3 satellites [30,39]. Referring to the relevant research on BDS-2, this section focuses on analyzing the relationship between the MPs and satellite elevations for BDS-3 MEO and IGSO satellites. It should be noted that this analysis only examined their correlation and did not investigate the modeling algorithms, as that would require other processing methods, which will be addressed in future research. For the analysis of this correlation, the MP sequences extracted using the QF mode in Section 3.1 were utilized. Furthermore, in order to avoid MPs being applied to specific satellites, uniform processing was applied to a specific type of satellite at each station to reflect its statistical characteristics. As of February 2023, the distribution of BDS-3 MEO and IGSO satellites can be found in Table 8. The satellites marked with bold numbers indicate their current status as on-orbit testing and do not participate in statistics.
Taking IGSO satellites as an example, assuming that a station observes a total of n IGSO satellites during a certain period and that the effective dataset of the fi frequency for each satellite is X f i n = { X 1 , X 2 , , X n 1 , X n } , the dataset composed of the satellite elevation (ele) and MP for all satellites and all epochs is Q B D S 2 I G S O = { e l e ( X f i n ) M P ( X f i n ) } . Due to the large amount of data, the elevation is divided into different ranges with steps of 0.1°, namely: d0 = 0°, d1 = 0.1°, d2 = 0.2°, dk = (k/10)°, …, d900 = 90°. The dataset of all the epochs in Q B D S 2 I G S O is selected with the elevation between dk−1 and dk, and a new dataset Q B D S 2 I G S O ( d k ) = { e l e ( d k ) M P ( d k ) } is formed. The average of M P ( d k ) is taken as the MP corresponding to the fi frequency when the height angle is dk’ = dk − 0.05°.
By using this method, a dataset consisting of elevation and MP can be obtained for each type of satellite and frequency, with both the elevation and MP containing 900 elements. Considering that different stations receive satellites at different elevations and that the maximum and minimum elevation for receiving a particular satellite may vary, there may be cases in the dataset where the MP is empty for certain elevation angles. In such cases, the maximum elevation in Q ( d k ) can be adaptively set based on the specific satellite reception conditions at each station. Figure 11 and Figure 12 provide examples of the statistical results of the MPs and satellite elevations for MAYG and OWMG stations, respectively.
The analysis of the statistical results for all the stations in Figure 2 reveals that the MPs of the MEO and IGSO satellites of BDS-3 exhibited similar systematic biases related to elevation, similar to BDS-2 satellites. Both of them showed more significant changes in MEO satellite systematic biases with elevation. The correlation strength between the MP and elevation varied significantly among different types of satellites and frequencies. The anti-multipath performance of BDS-3 at various frequencies was better than that of BDS-2. The MP of BDS-3 fluctuated around −0.5 m-0.5 m, while BDS-2 mainly fluctuated around −1.0 m-1.0 m. Comparing the different frequencies, it can be seen that the MP fluctuation amplitude of the B2a and B3I frequencies was smaller, while that of the B1C and B1I frequencies was larger. Additionally, the magnitude of the MP fluctuations varied across different satellite orbits for the same frequency. Overall, for BDS-3 satellites, the MP magnitude of the MEO and IGSO satellites was roughly equivalent, while for BDS-2 satellites, the magnitude of the MEO satellites was greater than that of the IGSO.
In order to quantify the correlation between the MP and elevation for each frequency, the correlation coefficient was calculated for different frequencies at each station. The calculation equation is [40]:
r X Y = i = 1 n ( X X ¯ ) ( Y Y ¯ ) i = 1 n ( X X ¯ ) 2 i = 1 n ( Y Y ¯ ) 2
where X and Y are the two vectors for the correlation to be determined, X ¯ and Y ¯ are the average values corresponding to each vector, and n is the length of the vector. The value of rxy is between −1 and 1, with rxy > 0 indicating a positive correlation and rxy < 0 indicating a negative correlation.
Generally, the absolute of r indicates the strength of the correlation. Typically, | r | 0.8 indicates a very strong correlation, 0.6 | r | <   0.8 indicates a strong correlation, 0.4 | r | <   0.6 indicates a moderate correlation, 0.2 | r | <   0.4 indicates a weak correlation, and | r | <   0.2 is uncorrelated [41]. Table 9 shows the average correlation coefficient of each experimental station. From Table 9, it can be observed that, for GEO satellites, the absolute values of the correlation coefficients for each frequency ranged from 0.06 to 0.13, indicating a clear lack of correlation. For IGSO and MEO satellites, the absolute values of the correlation coefficients for each frequency ranged from 0.42 to 0.56 and 0.63 to 0.80, respectively, showing moderate and strong correlations. Additionally, the correlation between the MEO and IGSO satellites of BDS-2 was stronger at the B1I and B3I frequencies than that of BDS-3. Among the four frequencies of BDS-3, the B2a frequency exhibited the strongest correlation.

4. Summary and Conclusions

The broadcasting of the multi-frequency observations in various navigation systems presents new opportunities for high-precision PNT services. However, the pseudorange and carrier phase observations of each frequency are subject to various interferences during signal propagation. As one of the errors affecting the pseudorange observations, MPs have also been studied by many scholars. On the basis of existing research, the MP characteristics of BDS-3 satellites at various frequencies were studied based on multi-frequency observations in this work. Through the analysis of static and dynamic observations, the main conclusions can be drawn as follows:
(1) Starting from the pseudorange and carrier phase observation equations, a direct multi-frequency MP combination coefficient calculation method based on the least squares principle was proposed. This method was simple in its calculation and could effectively utilize the observation information of each frequency. On this basis, a detailed formula for calculating the MP combination coefficients for eliminating I1,1 and I1,2 was derived, and a unified analytical expression was summarized when only eliminating I1,1. Compared with the traditional dual-frequency MP combination extraction results, the trend of the MP sequence extracted by the TF and QF combination was basically consistent with it, and the fluctuation amplitude was smaller, which confirmed the correctness of the proposed method.
(2) The accuracy of the MP sequences extracted from the four modes with different Ω was compared by using static and dynamic observations. The results showed that the magnitude of Ω directly affected the fluctuation amplitude of the extracted MP sequences. When there was a significant difference in Ω , the difference in the extracted MP sequences could reach the decimeter level. When the difference of Ω was small, the difference in the extracted MP sequences was generally within 3 cm. Compared to the DF1 mode, the modes with a smaller Ω showed an average improvement of over 25% in the accuracy of the extracted MP sequences, with an even higher accuracy improvement for dynamic observation. Additionally, the selected combinations were the optimal combination coefficients for each frequency combination, which also indirectly reflected the comparison of the accuracy between the optimal and non-optimal combinations in the MP extraction.
(3) Taking the triple-frequency and quad-frequency combinations as examples, the influence of I1,2 on the extraction of the MPs for the different frequencies and modes was analyzed by taking the first-order difference with the MP sequence obtained by only eliminating I1,1. The fluctuation amplitude of the differenced MP sequence also showed a clear correlation with Δ Ω , with the STD and range differing by 3–5 times. The B2a and B3I frequencies exhibited better anti-multipath performances. However, under the same level of combination noise, the STD and range of the differenced MP sequence of the B2a frequency were both more than 1.5 times larger than those of the other three frequencies, while the STD and range of the B3I frequency were the smallest.
(4) The impact of I1,2 on the IGSO and MEO satellites of BDS-3 was different. In different modes, the first-order difference sequences of the four frequencies of the IGSO satellites always exhibited a significant periodicity, while the periodicity of the MEO satellites was not significant. In addition, for the static observation, the magnitude of the differenced MP sequences between the two types of satellites was basically consistent, but for the dynamically collected observation, the STD and range of the IGSO satellites were about half those of the MEO satellites.
(5) Based on the observations from all the experimental stations, the correlation between the MPs of the IGSO and MEO satellites and the satellite elevation was analyzed. The MP of the MEO satellites showed a strong correlation with the elevation, while the correlation for the IGSO satellites was of moderate strength. The correlation between the MEO and IGSO satellites of BDS-2 at the B1I and B3I frequencies was stronger than that of BDS-3. Among the four frequencies of BDS-3, the B2a frequency exhibited the strongest correlation. Overall, the anti-multipath performance of BDS-3 frequencies was superior to that of BDS-2.
This work focused on the method of extracting MPs using multi-frequency observations and analyzed the characteristics of MPs under different scenarios. However, in practical positioning processes, it is necessary to consider how to correct such errors to improve the positioning accuracy. This involves addressing issues such as the extraction of the true information from MPs and the construction of high-precision correction models. These aspects are important considerations for multi-frequency and multi-system precise positioning and will be one of the future research directions.

Author Contributions

Conceptualization, Z.G. and X.Y.; methodology, Z.G. and C.H.; validation, Z.G. and C.H.; formal analysis, C.J. and H.T.; data curation, Z.G. and C.J.; writing—original draft preparation, Z.G. and M.Z.; writing—review and editing, Z.G., X.Y. and S.X.; supervision, X.Y.; project administration, Z.G. and X.Y.; funding acquisition, X.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Key Laboratory of Aviation-aerospace-ground Cooperative Monitoring and Early Warning of Coal Mining-induced Disasters of Anhui Higher Education Institutes (Anhui University of Science and Technology), grant number KLAHEI202202; National Natural Science Foundation of China, grant number 41474026; Key Research and Development Program of Anhui Province, grant number 202104a07020014; Major science and technology projects of Anhui Province, grant number 202103a05020026; Anhui Natural Science Foundation, grant number 2008085MD114; Key scientific research project of Suzhou University, grant number 2021yzd03, 2021yzd09; Coal Industry Engineering Research Center of Mining Area Environmental Additionally, Disaster Cooperative Monitoring(Anhui University of Science and Technology), grant number KSXTJC202203; First Class Undergraduate Major of Suzhou University, grant number szxy2021ylzy05; Surveying and Mapping Engineering Professional Teaching Team of Anhui Province, grant number 2020jxtd285.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Thanks to the reviewers and their suggestions for this work. The authors appreciate the International GNSS Monitoring and Assessment Service (iGMAS) and Multi-GNSS Experiment (MGEX) for the provision of relevant data. Thanks to the teacher of China University of Mining and Technology for providing the dynamic vehicle observation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yang, Y.; Mao, Y.; Sun, B. Basic performance and future developments of BeiDou global navigation satellite system. Satell. Navig. 2020, 1, 1. [Google Scholar] [CrossRef] [Green Version]
  2. Ge, Y.; Cao, X.; Lyu, D.; He, Z.; Ye, F.; Xiao, G.; Shen, F. An investigation of PPP time transfer via BDS-3 PPP-B2b service. GPS Solut. 2023, 27, 61. [Google Scholar] [CrossRef]
  3. Guo, Z.; Yu, X.; Hu, C.; Yu, Z.; Jiang, C. Correction Model of BDS-3 Satellite Pseudorange Multipath Delays and Its Impact on Single-Frequency Precise Point Positioning. Math. Probl. Eng. 2021, 2021, 9189541. [Google Scholar] [CrossRef]
  4. Cai, H.; Meng, Y.; Geng, C.; Gao, W.; Zhang, T.; Li, G.; Shao, B.; Xin, J.; Lu, H.; Mao, Y.; et al. BDS3 performance assessment: PNT, SBAS, PPP, SMC and SAR. Acta Geod. Cartogr. Sin. 2021, 50, 427–435. [Google Scholar]
  5. Ge, Y.; Wang, Q.; Wang, Y.; Lyu, D.; Cao, X.; Shen, F.; Meng, X. A new receiver clock model to enhance BDS-3 real-time PPP time transfer with the PPP-B2b service. Satell. Navig. 2023, 4, 8. [Google Scholar] [CrossRef]
  6. Chai, Y.; Ou, J.; Yuan, Y. Overview of Data Processing for BDS/GNSS Precise Positioning. Front. Data Comput. 2022, 4, 25–33. [Google Scholar]
  7. Zhang, P.; Ding, W.; Qu, X.; Yuan, Y. Simulation Analysis of LEO Constellation Augmented GNSS (LeGNSS) Zenith Troposphere Delay and Gradients Estimation. IEEE Trans. Geosci. Remote Sens. 2023, 61, 5800312. [Google Scholar] [CrossRef]
  8. Zhang, X.; Hu, J.; Ren, X. New progress of PPP/PPP-RTK and positioning performance comparison of BDS/GNSS PPP. Acta Geod. Cartogr. Sin. 2020, 49, 1084–1100. [Google Scholar]
  9. Bahadur, B. An improved weighting strategy for tropospheric delay estimation with real-time single-frequency precise positioning. Earth Sci. Inform. 2022, 15, 1267–1284. [Google Scholar] [CrossRef]
  10. Li, J.; Zhu, H.; Xu, A.; Xu, Z. Estimating ionospheric power spectral density for long-range RTK positioning using uncombined observations. GPS Solut. 2023, 27, 109. [Google Scholar] [CrossRef]
  11. Wu, Z.; Wang, Q.; Hu, C.; Yu, Z.; Wu, W. Modeling and assessment of five-frequency BDS precise point positioning. Satell. Navig. 2022, 3, 14. [Google Scholar] [CrossRef]
  12. Jin, S.; Su, K. PPP models and performances from single- to quad frequency BDS observations. Satell. Navig. 2020, 1, 13. [Google Scholar] [CrossRef]
  13. Zhang, Z. Code and phase multipath mitigation by using the observation-domain parameterization and its application in five-frequency GNSS ambiguity resolution. GPS Solut. 2021, 25, 144. [Google Scholar] [CrossRef]
  14. Yang, H.; Yang, Y.; Zhang, Z.; Sun, B.; Qin, W. Evaluation of the Effect of Higher-Order Ionospheric Delay on GPS Precise Point Positioning Time Transfer. Remote Sens. 2020, 12, 2129. [Google Scholar] [CrossRef]
  15. Zhao, Q.; Wang, G.; Liu, Z.; Hu, Z.; Dai, Z.; Liu, J. Analysis of BeiDou Satellite Measurements with Code Multipath and Geometry-Free Ionosphere-Free Combinations. Sensors 2016, 16, 123. [Google Scholar] [CrossRef] [Green Version]
  16. Li, J.; Yang, Y.; He, H.; Guo, H. Benefits of BDS-3 B1C/B1I/B2a triple-frequency signals on precise positioning and ambiguity resolution. GPS Solut. 2020, 24, 100. [Google Scholar] [CrossRef]
  17. Yuan, H.; Zhang, Z.; He, X.; Xu, T.; Zang, N. Real-time Cycle Slip Detection and Repair Method for BDS-3 Five-frequency Data. IEEE Access 2021, 9, 51189–51201. [Google Scholar] [CrossRef]
  18. Wang, H.; Pan, S.; Gao, W.; Ye, F.; Ma, C.; Tao, J.; Wang, Y. Real-time cycle slip detection and repair for BDS-3 triple-frequency and quad-frequency B1C/B1I/B3I/B2a signals. Acta Geodyn. Geomater. 2021, 3, 363–377. [Google Scholar] [CrossRef]
  19. Chen, D.; Ye, S.; Xia, J.; Liu, Y.; Xia, P. A geometry-free and ionosphere-free multipath mitigation method for BDS three-frequency ambiguity resolution. J. Geod. 2016, 90, 703–714. [Google Scholar] [CrossRef]
  20. Guo, Z.; Yu, X.; Hu, C.; Jiang, C.; Zhu, M. Research on Linear Combination Models of BDS Multi-Frequency Observations and Their Characteristics. Sustainability 2022, 14, 8644. [Google Scholar] [CrossRef]
  21. Chen, X.; Morton, Y.; Yu, W.; Truong, T. GPS L1CA/BDS B1I Multipath Channel Measurements and Modeling for Dynamic Land Vehicle in Shanghai Dense Urban Area. IEEE Trans. Veh. Technol. 2020, 69, 14247–14263. [Google Scholar] [CrossRef]
  22. Liu, C.; Tao, Y.; Xin, H.; Zhao, X.; Zhou, T. A Single-Difference Multipath Hemispherical Map for Multipath Mitigation in BDS-2/BDS-3 Short Baseline Positioning. Remote Sens. 2021, 13, 304. [Google Scholar] [CrossRef]
  23. Lu, R.; Chen, W.; Zhang, C.; Li, L.; Peng, Y.; Zheng, Z. Characteristics of the BDS-3 multipath effect and mitigation methods using precise point positioning. GPS Solut. 2022, 26, 41. [Google Scholar] [CrossRef]
  24. Chen, C.; Chang, G.; Zheng, N.; Xu, T. GNSS Multipath Error Modeling and Mitigation by Using Sparsity-Promoting Regularization. IEEE Access 2019, 7, 24096–24108. [Google Scholar] [CrossRef]
  25. Jin, T.; Zhang, P.; Yuan, H.; Cong, L.; Lu, X.; Qu, J. Two-Step Unbiased Estimation on Small-Delay Multipath of GNSS Signals Based on S-Curve Bias. IEEE Trans. Instrum. Meas. 2023, 72, 6500721. [Google Scholar] [CrossRef]
  26. Hauschild, A.; Montenbruck, O.; Sleewaegen, J.; Huisman, L.; Teunissen, P. Characterization of compass M-1 signals. GPS Solut. 2012, 16, 117–126. [Google Scholar] [CrossRef]
  27. Wanninger, L.; Beer, S. BeiDou satellite-induced code pseudorange variations: Diagnosis and therapy. GPS Solut. 2015, 19, 639–648. [Google Scholar] [CrossRef] [Green Version]
  28. Wang, X.; Liu, W.; Li, B.; Sun, G. A modified sidereal filtering for mitigating the code measurement multipath effects of BeiDou GEO satellites. Surv. Rev. 2016, 49, 428–436. [Google Scholar] [CrossRef]
  29. Wang, J.; He, X. Correction Model of BeiDou code systematic multipath errors and its impacts on signal-frequency PPP. Acta Geod. Cartogr. Sin. 2017, 46, 841–847. [Google Scholar]
  30. Rao, P.; Wang, Q. An improved code bias correction strategy for BDS. Sci. Surv. Mapp. 2020, 45, 1–10. [Google Scholar]
  31. Zhang, X.; Wu, M.; Liu, W.; Li, X.; Yu, S.; Lu, C.; Wickert, J. Initial assessment of the COMPASS/BeiDou-3: New-generation navigation signal. J. Geod. 2017, 91, 1224–1240. [Google Scholar] [CrossRef]
  32. He, Y. Quality Analysis of Satellite Signal for BDS-3 Simplest System. Geomat. Inf. Sci. Wuhan Univ. 2020, 45, 394–402. [Google Scholar]
  33. Wang, X.; Xu, X.; Zhao, Q. Signal quality analysis and orbit accuracy verification of BDS-3. Geomat. Inf. Sci. Wuhan Univ. 2023, 48, 611–619. [Google Scholar]
  34. Jin, T.; Hu, B.; Sun, Y.; Huang, Z.; Wang, Q.; Wu, Q. Optimal Solution to Multi-Frequency BDS Code-Multipath Combination Measurement. J. Navig. 2019, 72, 1297–1314. [Google Scholar] [CrossRef]
  35. Zhang, B.; Liu, T.; Xu, L.; Gao, R. Model algorithm and applications of uncombined precise point positioning for GNSS. J. Navig. Position. 2023, 11, 1–12. [Google Scholar]
  36. Estey, L.; Meertens, C. TEQC: The multi-purpose toolkit for GPS/GLONASS data. GPS Solut. 1999, 3, 42–49. [Google Scholar] [CrossRef]
  37. Chang, G.; Chen, C.; Yang, Y.; Xu, T. Tikhonov Regularization Based Modeling and Sidereal Filtering Mitigation of GNSS Multipath Errors. Remote Sens. 2018, 10, 1801. [Google Scholar] [CrossRef] [Green Version]
  38. Hu, C.; Wang, Q.; Wu, Z.; Guo, Z. A mixed multi-frequency precise point positioning strategy based on the combination of BDS-3 and GNSS multi-frequency observations. Meas. Sci. Technol. 2022, 34, 025008. [Google Scholar] [CrossRef]
  39. Li, X.; Zhang, X. Comparison of satellite-induced code bias between BDS-2 and BDS-3 satellite. J. Geod. Geodyn. 2018, 38, 191–194. [Google Scholar]
  40. Hu, C.; Wang, Z.; Rao, P.; Cheng, T. One-step correction strategy for BDS-2/BDS-3 satellite observation code bias and multipath delay. Acta Geod. Geophys. 2021, 56, 29–59. [Google Scholar] [CrossRef]
  41. Hu, D.; Yang, G.; Suo, X.; Liu, T.; Shao, M.; Cao, C. Spatiotemporal variation of water reserves and its correlation with precipitation in Qinghai province based on GRACE. J. Navig. Position. 2022, 10, 181–186. [Google Scholar]
Figure 1. Ω of different frequency combinations for MP extracting.
Figure 1. Ω of different frequency combinations for MP extracting.
Sensors 23 06151 g001
Figure 2. Distribution of experimental stations (Pentagram: MGEX station, triangle: IGMAS station).
Figure 2. Distribution of experimental stations (Pentagram: MGEX station, triangle: IGMAS station).
Sensors 23 06151 g002
Figure 3. Comparison of MP sequences extracted by different modes of KUN1 station.
Figure 3. Comparison of MP sequences extracted by different modes of KUN1 station.
Sensors 23 06151 g003
Figure 4. Frequency distribution of MP sequences extracted by different modes of KUN1 stations.
Figure 4. Frequency distribution of MP sequences extracted by different modes of KUN1 stations.
Sensors 23 06151 g004
Figure 5. Comparison of MP sequences extracted by CAR1 station with different modes.
Figure 5. Comparison of MP sequences extracted by CAR1 station with different modes.
Sensors 23 06151 g005
Figure 6. Frequency distribution of MP sequences extracted by different modes of CAR1 stations.
Figure 6. Frequency distribution of MP sequences extracted by different modes of CAR1 stations.
Sensors 23 06151 g006
Figure 7. The differenced MP sequences of C28 satellite at WUH2 station under each mode.
Figure 7. The differenced MP sequences of C28 satellite at WUH2 station under each mode.
Sensors 23 06151 g007
Figure 8. The differenced MP sequences of C38 satellite at WUH2 station under each mode.
Figure 8. The differenced MP sequences of C38 satellite at WUH2 station under each mode.
Sensors 23 06151 g008aSensors 23 06151 g008b
Figure 9. The differenced MP sequences of C23 satellite at CAR1 station under each mode.
Figure 9. The differenced MP sequences of C23 satellite at CAR1 station under each mode.
Sensors 23 06151 g009
Figure 10. The differenced MP sequences of C38 satellite at CAR1 station under each mode.
Figure 10. The differenced MP sequences of C38 satellite at CAR1 station under each mode.
Sensors 23 06151 g010
Figure 11. Relationship between satellite elevation and MP of MAYG station.
Figure 11. Relationship between satellite elevation and MP of MAYG station.
Sensors 23 06151 g011
Figure 12. Relationship between satellite elevation and MP of OWMG station.
Figure 12. Relationship between satellite elevation and MP of OWMG station.
Sensors 23 06151 g012aSensors 23 06151 g012b
Table 1. MP combination coefficients of B1I frequency.
Table 1. MP combination coefficients of B1I frequency.
Frequency NumberP Φ Ω
B1IB1IB1CB3IB2bB2a
Eliminate I1,1
DF1.000109.502−110.502///155.568
1.000−4.887/3.887//6.245
1.000−3.974//2.974/4.964
1.000−3.629///2.6294.481
TF1.000−2.292−2.5073.799//5.096
1.000−1.741−1.844//2.5853.621
1.000−3.754/0.386/2.3684.455
QF1.000−1.871−1.9640.782/2.0533.491
FF1.000−1.870−1.9480.3561.0401.4223.244
Simultaneously eliminate I1,1 and I1,2
TF1.000105.310−106.4520.142//149.741
1.000106.580−107.647//0.068151.483
1.000−7.784/12.834/−6.05016.183
QF1.000−2.763−4.72612.270/−5.78114.627
FF1.000−2.755−4.68411.5461.583−6.68914.494
Table 2. MP combination of B1C and B1I frequencies selected in the experiment.
Table 2. MP combination of B1C and B1I frequencies selected in the experiment.
Combination ModeCombination FrequenciesCombination Coefficients Ω
B1I
DF1(B1I, B1C)(1, 109.502, −110.502)155.568
DF2(B1I, B3I)(1, −4.887, 3.887)6.245
TF(B1I, B1C, B3I)(1, −2.292, −2.507, 3.799)5.096
QF(B1I, B1C, B3I, B2a)(1, −1.871, −1.964, 0.782, 2.053)3.491
B1C
DF1(B1C, B1I)(1, −109.502, 108.502)154.154
DF2(B1C, B3I)(1, −4.687, 3.687)5.964
TF(B1C, B1I, B3I)(1, −2.489, −2.276, 3.765)5.054
QF(B1C, B1I, B3I, B2a)(1, −1.951, −1.858, 0.774, 2.035)3.464
Table 3. Accuracy statistics of MP sequences extracted by different modes of KUN1 stations.
Table 3. Accuracy statistics of MP sequences extracted by different modes of KUN1 stations.
PRNDF1DF2TFQF
B1I
STD/mC280.2700.161/40.3%0.160/40.3%0.157/41.8%
C420.3770.231/38.7%0.228/39.5%0.223/40.8%
mean/mC28−7.744 × 10−168.808 × 10−161.182 × 10−15−1.744 × 10−16
C42−1.145 × 10−14−1.001 × 10−15−1.695 × 10−152.270 × 10−15
B1C
STD/mC280.3960.342/13.6%0.341/13.9%0.341/13.9%
C420.4630.393/15.1%0.395/14.7%0.393/15.1%
mean/mC283.577 × 10−16−1.773 × 10−161.750 × 10−151.833 × 10−16
C424.006 × 10−15−3.303 × 10−162.608 × 10−16−1.716 × 10−15
Table 4. Accuracy statistics of MP sequences extracted by different modes of CAR1 stations.
Table 4. Accuracy statistics of MP sequences extracted by different modes of CAR1 stations.
PRNDF1DF2TFQF
B1I
STD/mC230.4310.270/37.4%0.272/36.9%0.263/39.0%
C380.2270.159/30.0%0.158/30.4%0.158/30.4%
mean/mC23−2.80 × 10−15−5.78 × 10−16−1.54 × 10−161.56 × 10−16
C382.67 × 10−146.48 × 10−164.70 × 10−188.41 × 10−16
B1C
STD/mC230.4250.228/46.4%0.220/48.2%0.213/49.9%
C380.2150.121/43.7%0.121/43.7%0.120/44.0%
mean/mC232.43 × 10−151.04 × 10−15−6.68 × 10−16−3.39 × 10−16
C338.28 × 10−151.30 × 10−17−1.22 × 10−16−6.51 × 10−16
Table 5. Combination modes of simultaneous eliminate I1,1 and I1,2 for B1C, B3I, and B2a frequencies.
Table 5. Combination modes of simultaneous eliminate I1,1 and I1,2 for B1C, B3I, and B2a frequencies.
ModeCombination
Frequencies
Combination Coefficients Ω I 1 , 1 Ω I 1 , 2 Δ Ω
I1,1I1,1 + I1,2
B1CTF1B1C, B1I, B3I1, −2.489, −2.276, 3.7651, −113.491, 112.631, −0.1405.054159.894154.84
B1CTF2B1C, B1I, B2a1, −1.832, −1.729, 2.5611, −112.314, 111.381, −0.0673.593158.178154.585
B1CTF3B1C, B3I, B2a1, −3.644, 0.389, 2.2551, −7.474, 12.499, −6.0254.30315.76011.457
B1CQF1B1C, B1I, B3I, B2a1, −1.951, −1.858, 0.774, 2.0351, −4.845, −2.793, 12.812, −6.1743.46415.28211.818
B3ITF1B3I, B1C, B1I1, 4.781, −3.023, −2.7581, 10.181, 150.478, −161.6596.293221.090214.797
B3ITF2B3I, B1C, B2a1, 0.550, −4.379, 2.8291, −7.960, −1.688, 8.6485.24211.8746.632
B3ITF3B3I, B1I, B2a1,0.538,−4.479,2.9411, −7.759, −1.793, 8.5525.38511.6856.3
B3IQF1B3I, B1C, B1I, B2a1,1.012,−2.346,−2.231,2.5651, −7.786, −0.230, −1.549, 8.5654.25211.6807.428
B2aTF1B2a, B1C, B1I1, 3.573, −2.355, −2.2181, 7.734, 172.586, −181.3204.820250.446245.626
B2aTF2B2a, B1I, B3I1, 3.215, −4.827, 0.6121, 17.543, 2.033, −20.5765.83227.11521.283
B2aTF3B2a, B1C, B3I1, 3.094, −4.718, 0.6241, 17.434, 1.913, −20.3475.67726.86321.186
B2aQF1B2a, B1C, B1I, B3I1, 2.810, −2.528, −2.404, 1.2121, 17.394, 2.614, −0.744, −20.2644.61826.84322.225
Table 6. Accuracy statistics of MP difference sequences in different modes of experimental stations (unit: m).
Table 6. Accuracy statistics of MP difference sequences in different modes of experimental stations (unit: m).
ModeMeanSTDRangeModeMeanSTDRange
B1ITF1−0.0060.2471.599B1CTF1−0.0130.2681.689
B1ITF2−0.0060.2521.633B1CTF2−0.0120.2551.605
B1ITF3−0.0010.0630.451B1CTF30.0020.0670.441
B1IQF10.0000.0590.441B1CQF10.0010.0650.429
B3ITF1−0.0080.3662.275B2aTF10.0030.4012.478
B3ITF2−0.0010.0520.389B2aTF2−0.0030.1130.724
B3ITF3−0.0050.0490.367B2aTF3−0.0030.1120.713
B3IQF1−0.0050.0510.382B2aQF1−0.0040.1140.740
Table 7. Accuracy statistics of MP sequences in different modes of CAR1 station (unit: m).
Table 7. Accuracy statistics of MP sequences in different modes of CAR1 station (unit: m).
ModeMeanSTDRangeModeMeanSTDRange
B1ITF1−0.0020.2751.665B1CTF10.0000.2881.712
B1ITF2−0.0030.2801.683B1CTF20.0000.2741.626
B1ITF3−0.0030.0600.393B1CTF3−0.0010.0610.396
B1IQF10.0000.0560.359B1CQF10.0000.0600.389
B3ITF10.0010.3812.274B2aTF10.0030.4152.466
B3ITF20.0000.0440.294B2aTF2−0.0010.1000.610
B3ITF3−0.0010.0410.284B2aTF30.0010.1040.633
B3IQF10.0000.0450.307B2aQF10.0010.1050.635
Table 8. Satellite distribution in each orbit of BDS-3.
Table 8. Satellite distribution in each orbit of BDS-3.
BDS-2 (B1I, B3I,B2I)BDS-3 (IGSO/MEO: B1C,B1I,B3I,B2a)
IGSOC06, C07, C08, C09,
C10, C13, C16
C31, C38, C39, C40, C56
MEOC11, C12, C14C19, C20, C21, C22, C23, C24, C25, C26, C27, C28, C29, C30, C32, C33, C34, C35, C36, C37, C41, C42, C43, C44, C45, C46, C57, C58
Table 9. The correlation coefficient between MP and satellite elevation.
Table 9. The correlation coefficient between MP and satellite elevation.
BDS-2BDS-3
B1IB3IB1CB1IB3IB2a
GEO−0.06−0.13−0.09−0.10
IGSO−0.56−0.540.42−0.42−0.440.51
MEO−0.79−0.800.63−0.700.630.78
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, Z.; Yu, X.; Hu, C.; Jiang, C.; Tan, H.; Zhu, M.; Xie, S. The Extraction and Characterization of Pseudorange Multipath Based on BDS-3 Multi-Frequency Observations. Sensors 2023, 23, 6151. https://doi.org/10.3390/s23136151

AMA Style

Guo Z, Yu X, Hu C, Jiang C, Tan H, Zhu M, Xie S. The Extraction and Characterization of Pseudorange Multipath Based on BDS-3 Multi-Frequency Observations. Sensors. 2023; 23(13):6151. https://doi.org/10.3390/s23136151

Chicago/Turabian Style

Guo, Zhongchen, Xuexiang Yu, Chao Hu, Chuang Jiang, Hao Tan, Mingfei Zhu, and Shicheng Xie. 2023. "The Extraction and Characterization of Pseudorange Multipath Based on BDS-3 Multi-Frequency Observations" Sensors 23, no. 13: 6151. https://doi.org/10.3390/s23136151

APA Style

Guo, Z., Yu, X., Hu, C., Jiang, C., Tan, H., Zhu, M., & Xie, S. (2023). The Extraction and Characterization of Pseudorange Multipath Based on BDS-3 Multi-Frequency Observations. Sensors, 23(13), 6151. https://doi.org/10.3390/s23136151

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