Next Article in Journal
Agrobot Lala—An Autonomous Robotic System for Real-Time, In-Field Soil Sampling, and Analysis of Nitrates
Previous Article in Journal
Benchmarking Object Detection Deep Learning Models in Embedded Devices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Improving Satellite Positioning Precision Based on Multi-Frequency Navigation Signals

1
National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100101, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(11), 4210; https://doi.org/10.3390/s22114210
Submission received: 19 April 2022 / Revised: 26 May 2022 / Accepted: 30 May 2022 / Published: 31 May 2022
(This article belongs to the Section Remote Sensors)

Abstract

:
In satellite positioning systems, optimizing navigation satellite constellation and reducing the observation residuals are usually adopted to improve positioning precision and accuracy of the receiver. This paper presents a method to improve positioning precision by using multi-frequency navigation signals. The observation data of CAPS and GPS system are used to simulate the experiment. When the number of downlink frequencies is different, the root mean square of positioning error, improvement percentage, and standard deviation are calculated, respectively. When the number of descending frequencies is k, the root mean square of positioning error in three-dimensional space is 1/ k of that in single frequency. The theoretical derivation and experiment show that the precision of satellite positioning can be effectively improved by using multi-frequency navigation signals. The research work can provide theoretical support and data reference for the future research of satellite positioning.

1. Introduction

In satellite positioning systems, positioning precision of the receiver is mainly affected by spatial layout of navigation satellites constellation and measurement errors [1]. Measurement errors generally include satellite clock error, ephemeris error, and errors caused by ionospheric delay, troposphere delay, multipath effect, and receiver thermal noise. The methods such as wide-area difference technology, local-area difference technology, navigation message enhancement, high-precision orbit/clock correction, atmospheric delay correction, and measurement noise suppression are generally adopted to eliminate or weaken the pseudo-range observation error and improve the positioning precision of the receiver. In papers [2,3,4], they introduced the methods of wide-area difference and local-area difference and client performance evaluation. In papers [5,6,7], the observation residual caused by atmospheric delay was improved by improving modeling accuracy. In paper [8], a particle filter tracking algorithm was proposed to suppress multi-path errors under strong and weak signals. The spatial layout of satellite constellation has an important influence on the system performance and positioning precision, and dilution of precision (DOP) is the main parameter to measure the geometric relationship between satellite and receiver. In paper [9,10], the relationship between satellite positioning error and optimal constellation layout was analyzed. By adding satellites and improving the constellation structure, the DOP value and positioning precision can be improved. In the paper [11], a satellite selection method based on NSGA-II genetic algorithm is proposed, which takes into account the Geometric DOP (GDOP) value and the number of selected satellites. In paper [12], a GNSS multi-constellation selection method and an improved GPS real-time satellite selection algorithm were proposed. In paper [13], a constellation selection scheme for the regional satellite positioning system was proposed. In paper [14], it analyzed the influence of multi-constellation on the performance of the navigation system. In paper [15], it demonstrated that the Position DOP (PDOP) value of the GNSS system and the positioning precision can be improved by using LEO satellites. In paper [16], it introduced a GNSS analysis software based on MATLAB software for evaluating positioning results, such as positioning error, number of satellites, DOP value, etc. In paper [17], the weighted GDOP (WDOP) was proposed as an indicator for selecting stars and evaluating positioning precision.
Chinese Area Positioning System (CAPS) is a satellite positioning system based on communication satellites. Satellite navigation signals are broadcast from ground stations and forwarded by communication satellites to reach ground receivers to achieve positioning and navigation [18]. The principle of implementation of the system is shown in Figure 1. The communication satellites in the system include geostationary orbit (GEO) satellites and small inclined geostationary orbit (SIGSO) satellites. During the CAPS system test from 2019 to 2021, the ASIASAT-3S, CHINASAT-10, CHINASAT-12, and APSTAR-6C communication satellites were selected as navigation satellites. Each satellite has more than 20 communication transponders used to transmit CAPS navigation signals. The developed CAPS receiver has 48 signal channels and can simultaneously receive 48 downlink navigation signals. When the constellation layout is not good, three communication satellites can be used to achieve positioning with the aid of elevation [19,20]; when the constellation layout is good, the conventional least squares method can be used to achieve positioning. Combining the forwarding characteristics of the CAPS system and the resource advantages of communication satellites, Li et al. proposed the virtual clock method to effectively eliminate or greatly weaken the ephemeris error of the CAPS system [21]. Ai et al. proposed a method called physical augmentation factor of precision (PAFP) to equivalently improve the spatial layout [22].
The communication signal covers a wide range and has a wide signal frequency band. With the development of communication and positioning technology and as the variety of user needs continues to increase, the combination of communication and navigation has also become a hot research topic. This paper combines the characteristics of the CAPS system and the existing technology to carry out research on using multi-frequency navigation signals to improve positioning precision. Without increasing the number of satellites, by increasing the number of satellite downlink navigation signals, the pseudo-range equations of multi-frequency navigation signals are established to improve the positioning precision. The simulation experiments were carried out using the CAPS measured data and the observation data of the International GNSS Service (IGS) station in Fangshan, Beijing. According to the experimental results, the PDOP value and positioning precision are statistically analyzed to evaluate their improvement. This research provides theoretical support and data reference for the research on the positioning system combining communication satellites and navigation and has certain research significance and application value.

2. Methodology and Equations

Under the assumption that the measurement errors are independent and have equal variance, B. Parkinson et al. proposed that the optimal layout of six satellites is that four satellites are equally spaced on the horizontal plane, and the other two are on the zenith; at this time, the minimum value of the geometric dilution of precision (GDOP) is 2 , as shown in Figure 2 [23]. It is speculated from this that when the optimal layout of N satellites is known, and two satellites are placed in the same position where there was one satellite before, the optimal layout of 2N satellites can be obtained and, at this time, the GDOP is 1 / 2 times the previous.
The CAPS system adopts the virtual clock technology with the pseudo-range differential effect proposed by Li [21], which can eliminate the ephemeris error, clock error, and radial error caused by ionospheric delay and tropospheric delay in measurement errors. Only random errors are left, caused by multi-path effects or receiver thermal noise, etc., which makes the measurement errors of multi-frequency signals independent. It is technically difficult to arrange two satellites at the same position, but it is easy to achieve multi-frequency navigation signals from the same satellite, as shown in Figure 3. The physical augmentation factor of precision (PAFP) proposed by Ai et al. is to improve the spatial layout equivalently through multi-frequency navigation signals [22]. At this time, the PDOP value of the constellation composed of satellites transmitting multi-frequency navigation signals is significantly lower than that of the original system.

2.1. Virtual Clock Technology

The virtual clock has a pseudo-range differential effect; in the implementation of virtual clock technology, the integrated baseband of the ground station of the CAPS system is synchronized with the atomic clock of CAPS; the navigation signals are sent from the integrated baseband of the ground station and transmitted to the satellites through the radio frequency transmission channel [21]. The satellite transponder forwards the signals and broadcasts them to the CAPS coverage area, and the user receiver receives the navigation signals to achieve positioning, and the ground station itself also receives the navigation signals at the same time, as shown in Figure 4.
The observed pseudo-range of the user receiver is:
ρ = c · ( t c t v )
where c is the speed of light; t c is the total delay of the signal sent from the integrated baseband transmitting terminal of the ground station, through the radio frequency transmission channel, space uplink, satellite forwarding, and space downlink to subscriber receiver; t v is the time correction of the virtual clock.
The time correction of the virtual clock is obtained through the virtual clock model, as follows [24]:
t v = a 4 ( t t 0 ) 4 + a 3 ( t t 0 ) 3 + a 2 ( t t 0 ) 2 + a 1 ( t t 0 ) + a 0
where t is the time when the subscriber receiver receives the signals; t 0 is the start time of the virtual clock model; a 0 , a 1 , a 2 , a 3 , and a 4 are detected by the ground monitoring station.
The ground station forwards the model parameters to the subscriber receiver through the navigation message, and the CAPS navigation information is generated by the ground station, including time information, synchronization information, satellite orbit parameters, virtual clock model parameters, wide area differential information, local pseud-range differential information, local position difference data, etc.
Regarding the virtual clock deviation caused by the satellite position deviation, Li et al. also made theoretical analysis and actual measurement [21]. When there is 1 km difference between the actual position of the satellite and the predicted position of the satellite by ephemeris data, the receiver uses virtual clock technology to correct the ephemeris error, and the result of the positioning solution is almost the same as the result calculated by using the actual position of the satellite, which shows that the virtual clock technology can cancel or greatly weaken the ephemeris error.

2.2. Physical Augmentation Factor of Precision (PAFP)

Combining the resource advantages of the communication satellites of the CAPS system, Ai et al. proposed a method of using multiple frequencies to broadcast navigation signals separately, which can equivalently improve the spatial layout of the constellation [22]. The CAPS system selects communication band (including C-band and Ku-band) to transmit satellite navigation signals. There are more than 20 communication transponders on numerous communication satellite to transmit CAPS navigation signals. The working frequency bands and available frequency points are shown in Table 1.
Assuming that the number of visible satellites is n, each satellite sends down one navigation signal, and the coefficient matrix obtained by linearizing the pseudo-range observation equation is:
H = [ hx 1 hx 2 hx 3 hy 1 hy 2 hy 3 hz 1 hz 2 hz 3 hx i hy i hz i 1 1 1 1 ]
where [ hx i , hy i , hz i ] (i = 1, 2, 3, …, n) is the cosine vector from the receiver position to each satellite, and n is the number of satellites involved in the calculation. Then, the weight coefficient matrix is:
D = ( H T H ) 1 = [ i = 1 n hx i 2 i = 1 n hx i hy i i = 1 n hx i hz i i = 1 n hx i i = 1 n hx i hy i i = 1 n hy i 2 i = 1 n hy i hz i i = 1 n hy i i = 1 n hx i hz i i = 1 n hy i hz i i = 1 n hz i 2 i = 1 n hz i i = 1 n hx i i = 1 n hy i i = 1 n hz i n ] 1 = [ P 11 P 12 P 13 P 14 P 21 P 22 P 23 P 24 P 31 P 41 P 32 P 42 P 33 P 43 P 34 P 44 ]
Then, the PDOP value can be expressed as:
PDOP = P 11 + P 22 + P 33
where P ii (i = 1, 2, 3, 4) is the diagonal element of the weight coefficient matrix of the pseudo-range observation equation.
In a satellite positioning system, the variance of the pseudo-range observation error is amplified by the weight coefficient matrix to obtain the variance of the positioning error [1]. The relationship can be expressed as:
[ σ x 2 σ y 2 σ z 2 σ t 2 ] = [ P 11 P 22 P 33 P 44 ] σ e 2
where σ x 2 , σ y 2 , σ z 2 , and σ t 2 are the variance of the positioning errors in the X-axis, Y-axis, Z-axis directions, and clock difference in the ECEF coordinate system, respectively; P ii ( i   = 1, 2, 3, 4) is the diagonal element of the weight matrix of the pseudo-range observation equation.
Then, the standard deviation of the positioning error in three-dimensional space is:
σ = σ x 2 + σ y 2 + σ z 2 = P 11 + P 22 + P 33 σ e
When the number of visible satellites is n, each visible satellite sends down multi-frequency navigation signals, the number of downlink frequency is k, and its coefficient matrix is:
H B = [ hx 11 hy 11 hz 11 1 hx 21 hy 21 hz 21 1 hx i1 hy i1 hz i1 1 hx 12 hy 12 hz 12 1 hx 22 hy 22 hz 22 1 hx i2 hy i2 hz i2 1 hx ij hy ij hz ij 1 ]
where [ hx ij , hy ij , hz ij ] (i = 1, 2, …, n; j = 1, 2, …, k) is the cosine vector of the direction from the receiver to the i-th satellite that transmits the j-th signal, so:
hx i 1 = hx i 2 = = hx ij
hy i 1 = hy i 2 = = hy ij
hz i 1 = hz ij = = hz ij
According to Equations (6)–(9), the weight coefficient matrix of the pseudo-range observation equation can be obtained as:
D B = ( H B T H B ) 1 = [ k i = 1 n hx i 1 2 k i = 1 n hx i 1 hy i 1 k i = 1 n hx i 1 hz i 1 k i = 1 n hx i 1 k i = 1 n hx i 1 hy i 1 k i = 1 n hy i 1 2 k i = 1 n hy i 1 hz i 1 k i = 1 n hy i 1 k i = 1 n hx i 1 hz i 1 k i = 1 n hy i 1 hz i 1 k i = 1 n hz i 1 2 k i = 1 n hz i 1 k i = 1 n hx i 1 k i = 1 n hy i 1 k i = 1 n hz i 1 k · n ] 1 = 1 k [ P 11 P 12 P 13 P 14 P 21 P 22 P 23 P 24 P 31 P 41 P 32 P 42 P 33 P 43 P 34 P 44 ]
Then, its PDOP value can be expressed as:
PDOP B = P 11 + P 22 + P 33 k = PDOP k
It is seen that the ratio of the improved PDOP value to the original PDOP value is 1 k . When the original PDOP value is constant, the improved PDOP value is inversely proportional to the square root of the number of downlink frequencies. This ratio is called the physical augmentation factor of precision (PFAP).
Then, when each visible satellite sends down multiple-frequency navigation signals, the standard deviation of the positioning error in three-dimensional space is:
σ B = σ x 2 + σ y 2 + σ z 2 = P 11 + P 22 + P 33 k × σ e _ B

2.3. Positioning Solution of Multi-Frequency Receiver

The baseband system of the C/Ku dual-band navigation receiver in the user segment in the CAPS has 48 signal channels, which can simultaneously receive 48 downlink navigation signals, as shown in Figure 5. It uses virtual clock technology to correct satellite orbit errors, clock errors, and errors caused by atmospheric delay in the measurement pseudo-range and solves the receiver position by establishing pseudo-range equations to solve satellite positioning.
The process of satellite positioning using multi-frequency pseudo-range observations is as follows.
The pseudo-range observation equation of the navigation satellite can be expressed as:
ρ ij = c · ( t c t v ) = ( X i X ) 2 + ( Y i Y ) 2 + ( Z i Z ) 2 + ε ij
where ρ ij is the distance from the receiver to the i -th satellite that transmits the j -th signal, that is, the pseudo-range from the satellite to the receiver; c is the speed of light; t c is the total delay of the signal sent from the integrated baseband transmitting terminal of the ground station, through the radio frequency transmission channel, space uplink, satellite forwarding, and space downlink to the integrated baseband accept terminal of the subscriber receiver; t v is the time correction of the virtual clock; ( X i , Y i , Z i ) is the position coordinate of the i -th satellite; ( X , Y , Z ) is the receiver position coordinates; and ε ij is gaussian-distributed random errors due to multipath effects and receiver thermal noise, etc.
The pseudo-range observation equation is expanded by the first-order Taylor series at the initial point O ( X 0 , Y 0 , Z 0 ) . After omitting the higher-order terms, we can get:
ρ ij ρ i ( 0 ) = H · dX + ε ij
Its coefficient matrix is:
H = [ hx 11 hy 11 hz 11 1 hx 21 hy 21 hz 21 1 hx i1 hy i1 hz i1 1 hx 12 hy 12 hz 12 1 hx 22 hy 22 hz 22 1 hx i2 hy i2 hz i2 1 hx ij hy ij hz ij 1 ]
where [ hx ij , hy ij , hz ij ] (i = 1, 2, …, n; j = 1, 2, …, k) is the cosine vector of the direction from the receiver to the i-th satellite that transmits the j-th signal.
The difference between the observed pseudo-range and the predicted pseudo-range value at the initial point can be expressed as:
B = [ ρ 11 ρ 1 ( 0 ) ρ 21 ρ 2 ( 0 ) ρ i 1 ρ i ( 0 ) ρ 12 ρ 1 ( 0 ) ρ 22 ρ 2 ( 0 ) ρ i 2 ρ i ( 0 ) ρ ij ρ i ( 0 ) ]
where ρ ij is the pseudo-range observation value; ρ i ( 0 ) is the pseudo-range value predicted by the initial point, which can be expressed as:
ρ i ( O ) = ( X 0 X i ) 2 + ( Y 0 Y i ) 2 + ( Z 0 Z i ) 2
where ( X 0 , Y 0 , Z 0 ) are the coordinates of the initial point, and ( X i , Y i , Z i ) are the position coordinates of the i -th satellite.
The estimated value of the offset calculated by the least square algorithm is:
dX = [ Δ X Δ Y Δ Z Δ δ t ] = ( H T H ) 1 H T B
Add the current offset to the initial value to get the initial value of the next iteration ( X 0 + Δ X , Y 0 + Δ Y , Z 0 + Δ Z ). After reaching the threshold, stop the iteration and get the current receiver position estimate.

3. Experimental Data and Analysis Indicators

3.1. Observation Data of CAPS System

The CAPS system was in the test period in April 2021. The system uses two GEO satellites and a SIGSO satellite, namely APSTAR-7 (76.5° E), CHINASAT-10 (110.5° E), and ASIASAT-3S (147.5° E), to form a navigation constellation. The measured data of the CAPS system from 17:00 to 18:00 on 1 April 2021, with a sampling interval of one second, were collected at the Wuqing station in Tianjin, China. The measured pseudo-ranges in this data are corrected by virtual clock technology, which eliminates or greatly reduces correlated errors, including ephemeris errors, clock errors, errors caused by atmospheric delays, etc. We simulate the pseudo-range of multi-frequency measurement by adding random noise that satisfies the Gaussian distribution on the basis of the single-frequency measured pseudo-range. The position coordinates of the observation station are obtained by using the long-term static relative positioning method of commercial RTK.

3.2. Observational Data from GNSS Systems

To verify the degree of improvement of the positioning precision of this method on the GNSS system, the GPS satellite ephemeris data provided by the IGS station on 8 September 2021 at the Fangshan in Beijing, China for 24 h with a sampling interval of 30 s was selected as the satellite constellation of the simulation experiment. The visible number of satellites during the observation time is shown in Figure 6.
At present, the GNSS system can correct the ephemeris error and clock error by using the precise ephemeris and precise clock error products provided by IGS, and use the dual-frequency ionospheric elimination combined model to weaken the ionospheric error; use the Saastamonien model to correct the tropospheric delay. In the simulation experiment of GNSS, on the premise that the correlation models are used to weaken the correlation system error, a random error with a variance of 1 and a mean of 0 that satisfies the Gaussian distribution is added to the geometric distance from the position coordinates of the satellite to the position coordinates of the Fangshan station to simulate the measurement pseudo-ranges of multiple frequency points. The position coordinates of the observation station are the coordinates of the IGS annual solution.

3.3. Analysis Indicators

3.3.1. Root Mean Square Error (RMSE)

The root mean square error is used to measure the positioning precision of the positioning system in the three directions of east, north, zenith, and three-dimensional space. That is, the difference between the calculated estimated value and the reference value is calculated as:
RMSE = i = 1 n ( x i x 0 ) 2 n
3 DRMSE = i = 1 n ( x i x 0 ) 2 + ( y i y 0 ) 2 + ( z i z 0 ) 2 n
where i is the epoch time; n is the total epoch; x i ,   y i ,   z i are the estimated value scalculated in the three directions of the north, east, and zenith, respectively, at the i -th epoch time; and x 0 ,   y 0 ,   z 0 are the reference values of the observation station.

3.3.2. Improvement Percentage (IP)

The improvement percentage is used to reflect the improvement of the positioning precision; that is, when the number of downlink frequencies is different, the calculated positioning precision is improved relative to the single frequency. The calculation formula is:
IP = r i r 0 r 0 × 100 %
where r i is the positioning precision when the frequency number is i , i = 1, 2, …, n, n is the frequency number; and r 0 is the positioning precision when the frequency is single.

3.3.3. Standard Deviation (STD)

The standard deviation is used to reflect the degree of dispersion of the positioning results, that is, the difference between the calculated estimated value and the mean of the estimated value. The calculation formula is:
STD = i = 1 n ( x i x ¯ ) 2 n
where i is the epoch time; n is the total epoch; x i is the estimated value calculated at the i -th epoch time; and x ¯ is the average of the estimates computed at all epochs.

4. Results

The navigation receiver of the CAPS system has 48 signal channels and can receive at least 12 channels of C-band navigation signals at the same time. Assuming that the number of downlink frequencies of each visible satellite is at most 12, when the downlink frequencies are different, the maximum, minimum, and average values of PDOP in the CAPS system and GNSS system are counted respectively. The results are shown in Table 2.
To more intuitively represent the relationship between the PDOP value and the number of downlink frequencies, the changes of the PDOP value of CAPS and GNSS with the number of frequencies are plotted in Figure 7.
As shown in Table 2 and Figure 7, the use of multi-frequency navigation signals for satellite positioning can significantly reduce the PDOP value of the satellite navigation system. When the number of downlink frequencies is k, the PDOP value is about 1 / k of the original PDOP value. In the CAPS system, when the number of downlink frequency of each satellite is 1, that is, when multi-frequency is not used, the maximum value of the PDOP value is 6.85, the minimum value is 6.82, and the average value is 6.83; when the downlink frequencies of each satellite is 12, the maximum PDOP value drops to 1.98, the minimum value drops to 1.97, and the average value drops to 1.97. In the GNSS system, when the number of the downlink frequency of each satellite is 1, the maximum PDOP value is 2.91, the minimum value is 1.17, and the average value is 1.76; when the number is 12, its maximum PDOP value drops to 0.84, its minimum value drops to 0.34, and its average value drops to 0.51.
Based on the above satellite constellation space layout, take the observation station as the coordinate origin, and calculate the position coordinates of the receiver in the three directions of the north, east, and zenith at each moment. When the number of downlink frequencies is different, the root means square error, improvement percentage, and standard deviation of the three directions are calculated separately. The statistical results of the CAPS system are shown in Table 3, and the statistical results of the GNSS system are shown in Table 4; the variation of statistical results with the number of downlink frequencies are shown in Figure 8.
As shown in Table 3 and Table 4 and Figure 8, when the number of downlink frequencies is k, the positioning precision of the system is about 1/ k of the original. In the CAPS system, the elevation-assisted method is used for the positioning calculation, and the positioning error in the height is less than 0.01 m, so it is not counted here. When the number of downlink frequencies is 12, in the CAPS system, the positioning precision in the east and north directions are 0.57 and 2.18, respectively, which are 52% and 68% improved compared with the single frequency; the positioning precision of the two-dimensional space is improved by 67% compared with the single frequency. In the GNSS system, the positioning precision in the east, north, and zenith directions are 0.17, 0.23, 0.43, respectively, which are improved by 72%, 71%, and 71% compared with the single frequency; the positioning precision of the three-dimensional space is improved by 71% compared with the single frequency. The system standard deviation decreases as the number of frequencies increases, and it is seen that the accuracy and stability of the calculation results are improved.

5. Discussion

This paper introduces a method of using multi-frequency navigation signals to improve positioning precision and uses the observation data of the CAPS system and GNSS system to conduct simulation experiments. Statistical analysis was carried out on the experimental results in terms of root mean square error, improvement percentage, and standard deviation. It is verified that the method can improve the positioning precision, improve the stability of the positioning results, and can effectively reduce the PDOP value. The method of improving the positioning precision in this paper is different from the study of multi-scale ionospheric irregularity proposed in the paper [5]; it is also different from the improvement of ionospheric model in the paper [6]; it is also different from the improved solution of the whole week ambiguity proposed in paper [7]. Instead, the multi-channel navigation signal is received by the multi-frequency receiver, the distance from the satellite to the receiver is measured, and multiple pseudo-distance equations are established to achieve the purpose of positioning precision. This paper analyzes the improvement degree of multi-frequency on positioning precision from the time domain, and quantitatively analyzes the positioning results according to the number of downlink frequencies, which is different from the analysis of the influence of multi-frequency on the positioning results from the frequency domain in [25]. It is different from the one proposed by paper [8] to optimize satellite constellation from the angle of altitude angle to achieve the purpose of improving the PDOP value. It is also different from the weighted PDOP value method proposed in the paper [26] as the basis for selecting the best constellation. By increasing the number of downlink navigation signals of a single satellite without increasing the number of satellites, the method effectively improves the spatial layout of satellites and reduces the PDOP value. Through theoretical research and simulation analysis, it is found that when each satellite sends down multiple navigation signals, the ratio of the improved PDOP to the original PDOP is inversely proportional to the square root of the number of downlink navigation signals.

6. Conclusions

Based on the existing equipment and technology of the CAPS system, this paper introduces a method to improve the positioning precision through multi-frequency navigation signals, which can effectively reduce the PDOP value and improve the positioning precision. In the simulation results of this paper, in the CAPS system, when the number of downlink multi-frequency navigation signals from each satellite is 12, compared with single-frequency navigation signals, the average PDOP value drops from 6.83 to 1.97, the positioning precision in two-dimensional space is improved by 67%. Because the test data is collected in the CAPS test phase, the positioning error is large and does not represent the final positioning precision of the system. In the GNSS system, when the number of downlink frequencies is 12, and when the daily average PDOP of the satellite constellation is basically kept at about 2, the PDOP value can be reduced to less than 1 through multi-frequency navigation signals, and the positioning precision in 3D space is improved by 71%; the standard deviation decreases with the increase of the number of frequencies, and the improvement effect of positioning precision and stability is obvious. This research provides theoretical support and data reference for future satellite positioning technology combining communication satellites and navigation and has theoretical research and application value.

Author Contributions

T.K. performed the data collections, software, formal analysis, visualization and writing the manuscript. L.M. performed the methodology, formal analysis, investigation and writing the manuscript. G.A. performed the conceptualization, investigation and writing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number U1931125.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xie, G. Principles of GPS and Receiver Design, 1st ed.; Publishing House of Electronics Industry: Beijing, China, 2009; pp. 110–115. [Google Scholar]
  2. Chen, J.Y.; Liu, J.N.; Zhang, Y.P.; Hu, J.G.; Ge, M.R. On Distrubuted Wide Area Differential GPS Time Kinematic Positioning System. Acta Geod. Cartogr. Sin. 1998, 27, 1–8. [Google Scholar]
  3. Ashkenazi, V.; Hill, C.J.; Ochieng, W.Y.; Nagle, J. Wide-Area Differential GPS: A Performance Study. Navigation 1993, 40, 297–319. [Google Scholar] [CrossRef]
  4. Liu, X.L.; Wu, Y.; Liu, W.K. Correction of the Ephemeris and Satellite Clock of GPS Wide Area Augmentation System and User Performance Analysis. J. Geomat. 2019, 44, 79–84. [Google Scholar]
  5. Peng, Y.; Scales, W.A.; Hartinger, M.D.; Xu, Z.; Coyle, S. Characterization of multi-scale ionospheric irregularities using ground-based and space-based GNSS observations. Satell. Navig. 2021, 2, a14. [Google Scholar] [CrossRef]
  6. Sivavaraprasad, G.; Venkata Ratnam, D. Performance evaluation of ionospheric time delay forecasting models using GPS observations at a low-latitude station. Adv. Space Res. 2017, 60, 475–490. [Google Scholar] [CrossRef]
  7. Guo, J.M.; Zhou, M.D.; Shi, J.B.; Huang, C.J. FARSE scheme for single epoch GPS solution based on DUFCOM and DC algorithm and its performance analysis. Surv. Rev. 2014, 46, 426–431. [Google Scholar] [CrossRef]
  8. Wang, M.; Tao, S.; Zhang, W.W.; Hao, H. Analysis of BDS/GPS Signals’ Characteristics and Navigation Accuracy for a Geostationary Satellite. Remote Sens. 2021, 13, 1967. [Google Scholar] [CrossRef]
  9. Januszewski, J. Sources of Error in Satellite Navigation Positioning. Trans. Nav. Int. J. Mar. Navig. Saf. Sea Transp. 2017, 11, 419–423. [Google Scholar] [CrossRef]
  10. Wu, Q.; Wang, T.; Huang, Y.; Li, J. Analysis of Optimal Constellation Selection on Based on Satellite Positioning Error. J. Xi’an Jiaotong Univ. 2007, 10, 86–89. [Google Scholar]
  11. Xu, X.; Ma, L.; Ai, G. Satellite Selection with Muli-Objective Genetic Algorithm for Multi-GNSS Constellations. J. Shanghai Jiaotong Univ. 2017, 51, 281. [Google Scholar]
  12. Abedi, A.A.; Mosavi, M.R.; Mohammadi, K. A new recursive satellite selection method for multi-constellation GNSS. Surv. Rev. 2020, 52, 330–340. [Google Scholar] [CrossRef]
  13. Deng, C.; Zhuo, Y.; Wu, T. Constellation of a Regional Satellite Navigation System. J. Univ. Electron. Sci. Technol. China 2006, 5, 725–728. [Google Scholar]
  14. Yang, D.; Zeng, Q.; Liu, J.; Hu, Q. Analysis of performance of multi-constellation modernized satellite navigation system. Comput. Simul. 2019, 29, 65–68. [Google Scholar]
  15. Zhang, Y.; Li, Z.; Wang, Z.; Li, R.; Yuan, H. The improvement of BDS Observation Geometry with LEO constellations in Orbit Determination. Measurement 2021, 177, 109228. [Google Scholar] [CrossRef]
  16. Bahadur, B.; Nohutcu, M. PPPH: A MATLAB-based software for multi-GNSS precise point positioning analysis. GPS Solut. 2018, 22, 113. [Google Scholar] [CrossRef]
  17. Han, S. Quality-control issues relating to instantaneous ambiguity resolution for real-time GPS kinematic positioning. J. Geod. 1997, 71, 351–361. [Google Scholar] [CrossRef]
  18. Ai, G.X.; Shi, H.L.; Wu, H.T.; Yan, Y.H.; Bian, Y.J.; Hu, Y.H.; Li, Z.G.; Guo, J.; Cai, X.D. A positioning system based on communication satellites and the Chinese Area Positioning System (CAPS). Chin. J. Astron. Astrophys. 2008, 8, 611–630. [Google Scholar] [CrossRef] [Green Version]
  19. Jia, G.; Li, G.; Zhu, Q. Positionning Algorithm with Three Satellite Based on Chinese Area Position System. Comput. Appl. 2008, 28, 384–386. [Google Scholar]
  20. Wang, M.; Ma, L.; Zhang, L.; Ji, H.; Shi, H. Three-satellite Positioning Algorithm with Altitude Aiding for Regional Navigation Satellite System. J. Shanghai Jiaotong Univ. 2012, 46, 1647. [Google Scholar]
  21. Li, X.; Wu, H.; Bian, Y.; Wang, D. Virtual Satellite Atomic Clock with Pseudo-range Differential Function. Sci. China Ser. G Phys. Mech. Astron. 2008, 52, 353–359. [Google Scholar] [CrossRef]
  22. Ai, G.; Ma, L.; Shi, H.; Ji, H. A Method of Doubling Positioning Precision in Satellite Navigation and Positioning System. PRC Patent No. 201210090864.8. 2012.
  23. Parkinson, B.; Spilker, J., Jr.; Axelrad, P.; Enge, P. Global Positioning System Theory and Applications; American Institute of Aeronautics: Washington, DC, USA, 1996; pp. 202–207. [Google Scholar]
  24. Zhu, F. Development of Real-time Error Correction System of The Satellite Virtual Atomic Clock; University of Chinese Academy of Sciences (National Time Center, Chinese Academy of Sciences): Xi’an, China, 2011. [Google Scholar]
  25. Felski, A.; Jaskólski, K.; Zwolak, K.; Piskur, P. Analysis of Satellite Compass Error’s Spectrum. Sensors 2020, 20, 4067. [Google Scholar] [CrossRef] [PubMed]
  26. Wu, T.; Zhang, Y.; Liu, Y.; Yuan, G. BeiDou/GPS combination positioning methodology. J. Remote Sens. 2014, 18, 1087–1097. [Google Scholar]
Figure 1. Schematic diagram of the principle of the CAPS system.
Figure 1. Schematic diagram of the principle of the CAPS system.
Sensors 22 04210 g001
Figure 2. Layout of the six satellites with the smallest GDOP value.
Figure 2. Layout of the six satellites with the smallest GDOP value.
Sensors 22 04210 g002
Figure 3. Comparison of downlink single-frequency and multi-frequency signals from in-orbit satellites. (a) Each satellite sends down one navigation signal; (b) each satellite sends down multiple navigation signals with multi-frequencies.
Figure 3. Comparison of downlink single-frequency and multi-frequency signals from in-orbit satellites. (a) Each satellite sends down one navigation signal; (b) each satellite sends down multiple navigation signals with multi-frequencies.
Sensors 22 04210 g003
Figure 4. The realization process of the ground station transmitting and receiving the navigation signals.
Figure 4. The realization process of the ground station transmitting and receiving the navigation signals.
Sensors 22 04210 g004
Figure 5. CAPS navigation receiver: (a) is the physical diagram of the CAPS receiver, including multi-beam antenna, frequency converter, radio frequency device, signal channel, intelligent control terminal (including memory, CPU, display), power supply, etc.; (b) is a variable frequency navigation radio frequency device with 48 channels.
Figure 5. CAPS navigation receiver: (a) is the physical diagram of the CAPS receiver, including multi-beam antenna, frequency converter, radio frequency device, signal channel, intelligent control terminal (including memory, CPU, display), power supply, etc.; (b) is a variable frequency navigation radio frequency device with 48 channels.
Sensors 22 04210 g005
Figure 6. The situation of the visible satellites at the Fangshan station in Beijing on 8 September 2021.
Figure 6. The situation of the visible satellites at the Fangshan station in Beijing on 8 September 2021.
Sensors 22 04210 g006
Figure 7. The change of the PDOP value with the increase of the number of frequencies. Subfigures (a,b) are corresonding to the CAPS system and the GNSS system, respectively.
Figure 7. The change of the PDOP value with the increase of the number of frequencies. Subfigures (a,b) are corresonding to the CAPS system and the GNSS system, respectively.
Sensors 22 04210 g007
Figure 8. Positioning precision and improvement of multi-frequency navigation signals in CAPS system and GNSS system. Subfigures (a) at left column are standard deviation and the root mean square error in three directions and root mean square error of three-dimensional space change with the number of downlink frequencies in CAPS system; abd subfigures (b) at right column are standard deviation and the root mean square error in three directions and root mean square error of three-dimensional space change with the number of downlink frequencies in GNSS system.
Figure 8. Positioning precision and improvement of multi-frequency navigation signals in CAPS system and GNSS system. Subfigures (a) at left column are standard deviation and the root mean square error in three directions and root mean square error of three-dimensional space change with the number of downlink frequencies in CAPS system; abd subfigures (b) at right column are standard deviation and the root mean square error in three directions and root mean square error of three-dimensional space change with the number of downlink frequencies in GNSS system.
Sensors 22 04210 g008aSensors 22 04210 g008b
Table 1. Frequency band and the number of available frequency points.
Table 1. Frequency band and the number of available frequency points.
The Name of the BandUp-Band (MHz)Down-Band (MHz)Available Bandwidth (MHz)The Number of Available Frequency
Standard C-band5925~64253700~420050078
Extended C-band6425~67253400~3700500
Standard Ku-band14,000~14,50012,250~12,750500139
Extended Ku-band13,750~14,00010,950~11,200250
Planning Ku-Band12,750~13,25010,700~10,950, 11,200~11,450500
Broadcast Ku band17,300~17,80011,700~12,200500
Table 2. Statistics of PDOP values of multi-frequency navigation signals in CAPS system and GNSS system.
Table 2. Statistics of PDOP values of multi-frequency navigation signals in CAPS system and GNSS system.
The Number of Downlink FrequenciesCAPSGNSS
MaxMinAverageMaxMinAverage
16.856.826.832.911.171.76
24.844.824.832.060.831.24
33.963.943.941.680.681.01
43.433.413.411.450.580.88
53.063.053.051.30.520.79
62.82.782.791.190.480.72
72.592.582.581.10.440.66
82.422.412.411.030.410.62
92.282.272.280.970.390.59
102.172.162.160.920.370.56
112.072.062.060.880.350.53
121.981.971.970.840.340.51
Table 3. Statistical results of the CAPS system.
Table 3. Statistical results of the CAPS system.
The Number of Downlink FrequenciesRMS (m)IP (%)STD (m)
EN2DEN2DEN
11.26.836.93---2.243.91
20.934.794.882330301.562.77
30.813.974.053342421.272.32
40.753.483.563849491.12.05
50.693.163.234254530.991.86
60.662.943.024557560.91.76
70.642.752.824760590.841.64
80.622.592.664962620.781.56
90.612.512.584963630.751.52
100.62.432.55164640.691.49
110.592.332.415166650.671.43
120.572.182.265268670.621.35
Table 4. Statistical results of the GNSS system.
Table 4. Statistical results of the GNSS system.
The Number of Downlink FrequenciesRMS (m)IP (%)STD (m)
ENU3DENU3DENU
10.620.81.491.8----0.841.161.08
20.440.561.041.26292930300.580.810.76
30.360.460.841.03414343430.490.660.61
40.310.40.750.9505050500.420.580.54
50.280.360.650.79555556560.380.510.47
60.250.330.590.72595960600.330.470.43
70.230.30.570.68626262620.320.440.4
80.220.280.530.64656464640.290.410.39
90.210.260.490.59666767670.280.380.35
100.190.250.470.57696868680.260.360.34
110.180.250.440.54716970700.240.350.32
120.170.230.430.52727171710.240.330.31
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kong, T.; Ma, L.; Ai, G. Research on Improving Satellite Positioning Precision Based on Multi-Frequency Navigation Signals. Sensors 2022, 22, 4210. https://doi.org/10.3390/s22114210

AMA Style

Kong T, Ma L, Ai G. Research on Improving Satellite Positioning Precision Based on Multi-Frequency Navigation Signals. Sensors. 2022; 22(11):4210. https://doi.org/10.3390/s22114210

Chicago/Turabian Style

Kong, Ting, Lihua Ma, and Guoxiang Ai. 2022. "Research on Improving Satellite Positioning Precision Based on Multi-Frequency Navigation Signals" Sensors 22, no. 11: 4210. https://doi.org/10.3390/s22114210

APA Style

Kong, T., Ma, L., & Ai, G. (2022). Research on Improving Satellite Positioning Precision Based on Multi-Frequency Navigation Signals. Sensors, 22(11), 4210. https://doi.org/10.3390/s22114210

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