Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (124)

Search Parameters:
Keywords = LEO navigation

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
18 pages, 6204 KiB  
Article
An Integrity Monitoring Method for Navigation Satellites Based on Multi-Source Observation Links
by Jie Xin, Dongxia Wang and Kai Li
Remote Sens. 2024, 16(23), 4574; https://doi.org/10.3390/rs16234574 - 6 Dec 2024
Viewed by 298
Abstract
The BeiDou-3 navigation satellite system (BDS-3) has officially provided positioning, navigation, and timing (PNT) services to global users since 31 July 2020. With the application of inter-satellite link technology, global integrity monitoring becomes possible. Nevertheless, the content of integrity monitoring is still limited [...] Read more.
The BeiDou-3 navigation satellite system (BDS-3) has officially provided positioning, navigation, and timing (PNT) services to global users since 31 July 2020. With the application of inter-satellite link technology, global integrity monitoring becomes possible. Nevertheless, the content of integrity monitoring is still limited by the communication capacity of inter-satellite links and the layout of ground monitoring stations. Low earth orbit (LEO) satellites have advantages in information-carrying rate and kinematic velocity and can be used as satellite-based monitoring stations for navigation satellites. Large numbers of LEO satellites can provide more monitoring data than ground monitoring stations and make it easier to obtain full-arc observation data. A new challenge of redundant data also arises. This study constructs multi-source observation links with satellite-to-ground, inter-satellite, and satellite-based observation data, proposes an integrity monitoring method with optimization of observation links, and verifies the performance of integrity monitoring with different observation links. The experimental results show four findings. (1) Based on the integrity status of BDS-3, the proposed system-level integrity mode can realize full-arc anomaly diagnosis in information and signals according to the observation conditions of the target satellite. Apart from basic navigation messages and satellite-based augmentation messages, autonomous messages and inter-satellite ranging data can be used to evaluate the state of the target satellite. (2) For a giant LEO constellation, only a small number of LEO satellites need to be selected to construct a minimum satellite-based observation unit that can realize multiple returns of navigation messages and reduce the redundancy of observation data. With the support of 12 and 30 LEO satellites, the minimum number of satellite-based observation links is 1 and 4, respectively, verifying that a small amount of LEO satellites could be used to construct a minimum satellite-based observation unit. (3) A small number of LEO satellites can effectively improve the observation geometry of the target satellite. An orbit determination observation unit, which consists of chosen satellite-to-ground and/or satellite-based observation links based on observation geometry, is proposed to carry out fast calculations of satellite orbit. If the orbit determination observation unit contains 6 satellite-to-ground monitoring links and 6/12/60 LEO satellites, the value of satellite position dilution of precision (SPDOP) is 38.37, 24.60, and 15.71, respectively, with a 92.95%, 95.49%, and 97.12% improvement than the results using 6 satellite-to-ground monitoring links only. (4) LEO satellites could not only expand the resolution of integrity parameters in real time but also augment the service accuracy of the navigation satellite system. As the number of LEO satellites increases, the area where UDRE parameters can be solved in real time is constantly expanding to a global area. The service accuracy is 0.93 m, 0.88 m, and 0.65 m, respectively, with augmentation of 6, 12, and 60 LEO satellites, which is an 8.9%, 13.7%, and 36.3% improvement compared with the results of regional service. LEO satellites have practical application values by improving the integrity monitoring of navigation satellites. Full article
Show Figures

Figure 1

Figure 1
<p>Design of integrity monitoring system for the BDS-3 satellite.</p>
Full article ">Figure 2
<p>Design of integrity monitoring system with support of multi-source observation links.</p>
Full article ">Figure 3
<p>The chosen valid satellite-based observation links (<math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>≤</mo> <msub> <mi>γ</mi> <mn>0</mn> </msub> </mrow> </semantics></math>).</p>
Full article ">Figure 4
<p>The chosen valid satellite-based observation links (<math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>&gt;</mo> <msub> <mi>γ</mi> <mn>0</mn> </msub> </mrow> </semantics></math>).</p>
Full article ">Figure 5
<p>The chosen valid satellite-to-ground observation links.</p>
Full article ">Figure 6
<p>Processing of integrity monitoring system.</p>
Full article ">Figure 7
<p>Divided grid points and chosen ground monitoring stations.</p>
Full article ">Figure 8
<p>Multiple numbers in scenario 3.</p>
Full article ">Figure 9
<p>Coverage rate of four-multiple numbers in scenario 3.</p>
Full article ">Figure 10
<p>Coverage rate of more than four-multiple numbers in scenario 3. (The grids in the blue area can be monitored by no less than four LEO satellites; the grids in the yellow area can be monitored by less than four LEO satellites).</p>
Full article ">Figure 11
<p>Number of available links and SPDOP values in scenario 1.</p>
Full article ">Figure 12
<p>Number of available links and SPDOP values in scenario 2.</p>
Full article ">Figure 13
<p>Number of available links and SPDOP values in scenario 3.</p>
Full article ">Figure 14
<p>Number of available links and SPDOP values in scenario 4.</p>
Full article ">Figure 15
<p>Number of available links and PDOP values of BJFS in scenario 4.</p>
Full article ">
21 pages, 10795 KiB  
Article
COSMIC-2 RFI Prediction Model Based on CNN-BiLSTM-Attention for Interference Detection and Location
by Cheng-Long Song, Rui-Min Jin, Chao Han, Dan-Dan Wang, Ya-Ping Guo, Xiang Cui, Xiao-Ni Wang, Pei-Rui Bai and Wei-Min Zhen
Sensors 2024, 24(23), 7745; https://doi.org/10.3390/s24237745 - 4 Dec 2024
Viewed by 422
Abstract
As the application of the Global Navigation Satellite System (GNSS) continues to expand, its stability and safety issues are receiving more and more attention, especially the interference problem. Interference reduces the signal reception quality of ground terminals and may even lead to the [...] Read more.
As the application of the Global Navigation Satellite System (GNSS) continues to expand, its stability and safety issues are receiving more and more attention, especially the interference problem. Interference reduces the signal reception quality of ground terminals and may even lead to the paralysis of GNSS function in severe cases. In recent years, Low Earth Orbit (LEO) satellites have been highly emphasized for their unique advantages in GNSS interference detection, and related commercial and academic activities have increased rapidly. In this context, based on the signal-to-noise ratio (SNR) and radio-frequency interference (RFI) measurements data from COSMIC-2 satellites, this paper explores a method of predicting RFI measurements using SNR correlation variations in different GNSS signal channels for application to the detection and localization of civil terrestrial GNSS interference signals. Research shows that the SNR in different GNSS signal channels shows a correlated change under the influence of RFI. To this end, a CNN-BiLSTM-Attention model combining a convolutional neural network (CNN), bi-directional long and short-term memory network (BiLSTM), and attention mechanism is proposed in this paper, and the model takes the multi-channel SNR time series of the GNSS as the input and outputs the maximum measured value of RFI in the multi-channels. The experimental results show that compared with the traditional band-pass filtering inter-correlation method and other deep learning models, the model in this paper has a root mean square error (RMSE), mean absolute error (MAE), and correlation coefficient (R2) of 1.0185, 1.8567, and 0.9693, respectively, in RFI prediction, which demonstrates a higher RFI detection accuracy and a wide range of rough localization capabilities, showing significant competitiveness. Since the correlation changes in the SNR can be processed to decouple the signal strength, this model is also suitable for future GNSS-RO missions (such as COSMIC-1, CHAMP, GRACE, and Spire) for which no RFI measurements have yet been made. Full article
Show Figures

Figure 1

Figure 1
<p>The different signal transmission paths between the RFI source, the GNSS satellites, and the GNSS RO satellites (not to scale).</p>
Full article ">Figure 2
<p>SNR, S4 scintillation index, elevation angle, and RFI measurements for the POD 01 antenna of C2E1 satellite near 1:35 UTC on 1 January 2023: (<b>a</b>) SNR sequence of CA code L1 band for different channels. (<b>b</b>) S4 scintillation index. (<b>c</b>) Elevation angle of the LEO-GPS link. (<b>d</b>) Maximum value of RFI measurements in multiple channels.</p>
Full article ">Figure 3
<p>Changes in SNR, S4 scintillation index, and RFI measurements of the C2E1 satellite POD 01 antenna when scintillation occurs on 1 January 2023 near 1:10 UTC. (<b>a</b>) SNR sequence of CA code L1 band for different channels. (<b>b</b>) S4 scintillation index. (<b>c</b>) Maximum value of RFI measurements in multiple channels.</p>
Full article ">Figure 4
<p>Results of the interference detection algorithm: (<b>a</b>) The result of band-pass filtering and normalization of the multi-channel SNR sequence in <a href="#sensors-24-07745-f002" class="html-fig">Figure 2</a>a. (<b>b</b>) The calculated cross-correlation sequence after moving window normalization and filtering, as well as interference, can be detected by setting a threshold (set to 0.01 in this example).</p>
Full article ">Figure 5
<p>Spatial distribution of the orbits of GPS and COSMIC-2 satellites during the SNR duration in <a href="#sensors-24-07745-f002" class="html-fig">Figure 2</a>a.</p>
Full article ">Figure 6
<p>Calculation results of the RFI measurement sequence and cross-correlation sequence output by the C2E1 satellite 01 antenna on January 1, 2023 UTC. The two sequences also show a certain degree of correlation over the day: (<b>a</b>) RFI measurement sequence; (<b>b</b>) Calculated normalized cross-correlation sequence after band-pass filtering.</p>
Full article ">Figure 7
<p>Basic structure of a CNN model.</p>
Full article ">Figure 8
<p>LSTM model schematic.</p>
Full article ">Figure 9
<p>Schematic diagram of the BiLSTM model.</p>
Full article ">Figure 10
<p>Basic structure of the CNN-BiLSTM-Attention model used in this paper.</p>
Full article ">Figure 11
<p>Flowchart of the algorithm.</p>
Full article ">Figure 12
<p>Predicted RFI measurements and dRFI for the training set on the COSMIC-2 C2E1 satellite for the CNN-BiLSTM-Attention model. The dashed line in the figure indicates the interference threshold y = ±0.001, with the time resolution downsampled to 3 h: (<b>a</b>) the blue line is the true value of the RFI measurement, and the red line is the prediction result; (<b>b</b>) dRFI.</p>
Full article ">Figure 13
<p>Predicted RFI measurements and dRFI for the training set on the COSMIC-2 C2E1 satellite for the BiLSTM-Attention model. The dashed line in the figure indicates the interference threshold y = ±0.001, with the time resolution downsampled to 3 h: (<b>a</b>) the blue line is the true value of the RFI measurement, and the red line is the prediction result; (<b>b</b>) dRFI.</p>
Full article ">Figure 14
<p>Predicted RFI measurements and dRFI for the training set on the COSMIC-2 C2E1 satellite for the LSTM model. The dashed line in the figure indicates the interference threshold y = ±0.001, with the time resolution downsampled to 3 h: (<b>a</b>) the blue line is the true value of the RFI measurement, and the red line is the prediction result; (<b>b</b>) dRFI.</p>
Full article ">Figure 15
<p>Predicted RFI measurements and dRFI for the test set on the COSMIC-2 C2E1 satellite for the CNN-BiLSTM-Attention model. The dashed line in the figure indicates the interference threshold y = ±0.001, with the time resolution downsampled to 3 h: (<b>a</b>) the blue line is the true value of the RFI measurement, and the red line is the prediction result; (<b>b</b>) dRFI.</p>
Full article ">Figure 16
<p>Predicted RFI measurements and dRFI for the test set on the COSMIC-2 C2E1 satellite for the BiLSTM-Attention model. The dashed line in the figure indicates the interference threshold y = ±0.001, with the time resolution downsampled to 3 h: (<b>a</b>) the blue line is the true value of the RFI measurement, and the red line is the prediction result; (<b>b</b>) dRFI.</p>
Full article ">Figure 17
<p>Predicted RFI measurements and dRFI for the test set on the COSMIC-2 C2E1 satellite for the LSTM model. The dashed line in the figure indicates the interference threshold y = ±0.001, with the time resolution downsampled to 3 h: (<b>a</b>) the blue line is the true value of the RFI measurement, and the red line is the prediction result; (<b>b</b>) dRFI.</p>
Full article ">Figure 18
<p>Global RFI situation map of the six COSMIC-2 satellites’ 01 and 02 antenna superposition cases using the four methods mentioned in this paper and the actual measured values: (<b>a</b>) real measured values, and the three green dotted rectangles inside the marker are the primary sources of prediction error for various algorithms (<b>b</b>) CNN-BiLSTM-Attention, (<b>c</b>) BiLSTM-Attention, (<b>d</b>) LSTM, (<b>e</b>) normalized cross-correlation method with band-pass filtering.</p>
Full article ">Figure 18 Cont.
<p>Global RFI situation map of the six COSMIC-2 satellites’ 01 and 02 antenna superposition cases using the four methods mentioned in this paper and the actual measured values: (<b>a</b>) real measured values, and the three green dotted rectangles inside the marker are the primary sources of prediction error for various algorithms (<b>b</b>) CNN-BiLSTM-Attention, (<b>c</b>) BiLSTM-Attention, (<b>d</b>) LSTM, (<b>e</b>) normalized cross-correlation method with band-pass filtering.</p>
Full article ">
23 pages, 2184 KiB  
Article
Research on High-Dynamic Tracking Algorithms for FH-BOC Signals
by Xue Li, Shun Zhao, Xinyue Hou, Lulu Wang and Yinsen Zhang
Aerospace 2024, 11(12), 987; https://doi.org/10.3390/aerospace11120987 - 28 Nov 2024
Viewed by 324
Abstract
The rapid development of Low Earth Orbit (LEO) satellite navigation systems requires modulation schemes with strong anti-jamming capabilities, high spectral efficiency, and the ability to achieve precise tracking in high-dynamic environments. Traditional Binary Offset Carrier (BOC) modulation suffers from multi-peak ambiguity, leading to [...] Read more.
The rapid development of Low Earth Orbit (LEO) satellite navigation systems requires modulation schemes with strong anti-jamming capabilities, high spectral efficiency, and the ability to achieve precise tracking in high-dynamic environments. Traditional Binary Offset Carrier (BOC) modulation suffers from multi-peak ambiguity, leading to false lock issues. To address this, FH-BOC modulation, which integrates BOC modulationand frequency hopping, significantly improves both anti-jamming performance and spectral efficiency. Against this background, this paper proposes a comprehensive high-dynamic tracking algorithm for FH-BOC signals. (1) Based on the adaptive Kalman filter algorithm, high-precision carrier tracking was achieved in high-dynamic environments. (2) By leveraging the correlation between the ranging code and frequency-hopping offset carrier, a composite pseudo-code is generated through the XOR operation, and a corresponding composite code-tracking loop is introduced. (3) Based on code loop tracking results, the frequency-hopping moments of the subcarrier are detected, and a phase-locked loop for the frequency-hopping subcarrier is established. Simulation results indicate that the algorithm achieves centimeter-level pseudorange measurement accuracy for FH-BOC navigation signals under the JPL high-dynamic model. Full article
(This article belongs to the Section Astronautics & Space Science)
Show Figures

Figure 1

Figure 1
<p>FH-BOC transmitted signal generation process.</p>
Full article ">Figure 2
<p>The simulation results of the autocorrelation function and power spectral density of the FH-BOC signal.</p>
Full article ">Figure 3
<p>JPL high-dynamic signal model.</p>
Full article ">Figure 4
<p>Linear Kalman carrier tracking model.</p>
Full article ">Figure 5
<p>Flow chart of the linear Kalman carrier tracking algorithm.</p>
Full article ">Figure 6
<p>Linear Kalman tracking simulation results: (<b>a1</b>,<b>a2</b>) carrier tracking error under different CNRs; (<b>b1</b>,<b>b2</b>) carrier Doppler tracking error under different CNRs.</p>
Full article ">Figure 7
<p>Flow chart of the adaptive Kalman carrier tracking algorithm.</p>
Full article ">Figure 8
<p>Adaptive Kalman tracking simulation results: (<b>a1,a2</b>) carrier tracking error under different CNRs; (<b>b1,b2</b>) carrier Doppler tracking error under different CNRs.</p>
Full article ">Figure 8 Cont.
<p>Adaptive Kalman tracking simulation results: (<b>a1,a2</b>) carrier tracking error under different CNRs; (<b>b1,b2</b>) carrier Doppler tracking error under different CNRs.</p>
Full article ">Figure 9
<p>Simulation results of carrier tracking algorithms under different CNRs: (<b>a</b>) Carrier phase tracking and standard error; (<b>b</b>) Carrier Doppler tracking and standard error.</p>
Full article ">Figure 10
<p>FH-BOC modulation composite code-tracking loop.</p>
Full article ">Figure 11
<p>The second-order DLL loop filter.</p>
Full article ">Figure 12
<p>Composite code-tracking loop pseudorange measurement error under different CNRs: (<b>a</b>) CNR = 45 dB·Hz; (<b>b</b>) CNR = 50 dB·Hz; (<b>c</b>) CNR = 55 dB·Hz; (<b>d</b>) CNR = 60 dB·Hz.</p>
Full article ">Figure 13
<p>Frequency-hopping subcarrier phase-locked loop.</p>
Full article ">Figure 14
<p>Simulation results of the frequency-hopping subcarrier tracking loop: (<b>a</b>) pseudorange measurement error with CNR = 50 dB·Hz; (<b>b</b>) pseudorange measurement error with CNR = 60 dB·Hz; (<b>c</b>) FH subcarrier phase tracking error with CNR = 50 dB·Hz r. (<b>d</b>) FH-subcarrier phase tracking error with CNR = 60 dB·Hz.</p>
Full article ">Figure 14 Cont.
<p>Simulation results of the frequency-hopping subcarrier tracking loop: (<b>a</b>) pseudorange measurement error with CNR = 50 dB·Hz; (<b>b</b>) pseudorange measurement error with CNR = 60 dB·Hz; (<b>c</b>) FH subcarrier phase tracking error with CNR = 50 dB·Hz r. (<b>d</b>) FH-subcarrier phase tracking error with CNR = 60 dB·Hz.</p>
Full article ">Figure 15
<p>The standard deviation of pseudorange measurements under different CNRs.</p>
Full article ">
25 pages, 5533 KiB  
Article
Pulsed Orthogonal Time Frequency Space: A Fast Acquisition and High-Precision Measurement Signal for Low Earth Orbit Position, Navigation, and Timing
by Dong Fu, Honglei Lin, Ming Ma, Muzi Yuan and Gang Ou
Remote Sens. 2024, 16(23), 4432; https://doi.org/10.3390/rs16234432 - 27 Nov 2024
Viewed by 374
Abstract
The recent rapid development of low Earth orbit (LEO) constellation-based navigation techniques has enhanced the ability of position, navigation, and timing (PNT) services in deep attenuation and interference environments. However, existing navigation modulations face the challenges of high acquisition complexity and do not [...] Read more.
The recent rapid development of low Earth orbit (LEO) constellation-based navigation techniques has enhanced the ability of position, navigation, and timing (PNT) services in deep attenuation and interference environments. However, existing navigation modulations face the challenges of high acquisition complexity and do not improve measurement precision at the same signal strength. We propose a pulsed orthogonal time frequency space (Pulse-OTFS) signal, which naturally converts continuous signals into pulses through a special delay-Doppler domain pseudorandom noise (PRN) code sequence arrangement. The performance evaluation indicates that the proposed signal reduces at least 89.4% of the acquisition complexity. The delay measurement accuracy is about 8 dB better than that of the traditional binary phase shift keying (BPSK) signals with the same bandwidth. It also provides superior compatibility and anti-multipath performance. The advantages of fast acquisition and high-precision measurement are verified by processing the real signal in the developed software receiver. As Pulse-OTFS occupies only one time slot of a signal period, it can be easily integrated with OTFS-modulated communication signals and used as a navigation signal from broadband LEO satellites as an effective complement to the global navigation satellite system (GNSS). Full article
(This article belongs to the Special Issue LEO-Augmented PNT Service)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Block diagram of OTFS modulation.</p>
Full article ">Figure 2
<p>PSDs of the Pulse-OTFS signal.</p>
Full article ">Figure 3
<p>Comparison of PSDs for Pulse-OTFS signals with BPSK(5), and Pulse-BPSK(5).</p>
Full article ">Figure 4
<p>Comparison of ACFs for Pulse-OTFS signals with BPSK(5), and Pulse-BPSK(5).</p>
Full article ">Figure 5
<p>Simulation PSDs with different quantization bits for OTFS(10,1023).</p>
Full article ">Figure 6
<p>Ambiguity function envelope of Pulse-OTFS(10,1023).</p>
Full article ">Figure 7
<p>Multiplication and addition complexity ratio of Pulse-OTFS/Pulse-BPSK to OTFS/BPSK.</p>
Full article ">Figure 8
<p>Detection probability of different signals.</p>
Full article ">Figure 9
<p>Comparison of code-tracking performance for different signals. (<b>a</b>) Gabor bandwidth and (<b>b</b>) NELP DLL code-tracking error with a 20 MHz pre-filtering bandwidth.</p>
Full article ">Figure 10
<p>Comparison of S-curves for different signals with a 20 MHz pre-filtering bandwidth.</p>
Full article ">Figure 11
<p>Comparison of anti-multipath performance for different signals. (<b>a</b>) Multipath error envelope and (<b>b</b>) average multipath error envelope. <span class="html-italic">a</span> = −6 dB, <span class="html-italic">d</span> = 1 chip.</p>
Full article ">Figure 12
<p>Simulated detection probabilities for different signals.</p>
Full article ">Figure 13
<p>Flow chart of the experimental platform.</p>
Full article ">Figure 14
<p>Photograph of the experimental platform.</p>
Full article ">Figure 15
<p>Experimental results of the NELP DLL code-tracking error. (<b>a</b>) <span class="html-italic">T<sub>coh</sub></span> = 1 ms and (<b>b</b>) <span class="html-italic">T<sub>coh</sub></span> = 10 ms.</p>
Full article ">Figure 16
<p>Correlator peak magnitude for different signals.</p>
Full article ">
25 pages, 5692 KiB  
Article
Initial Design for Next-Generation BeiDou Integrity Subsystem: Space–Ground Integrated Integrity Monitoring
by Weiguang Gao, Lei Chen, Feiren Lv, Xingqun Zhan, Lin Chen, Yuqi Liu, Yongshan Dai and Yundi Jin
Remote Sens. 2024, 16(22), 4333; https://doi.org/10.3390/rs16224333 - 20 Nov 2024
Viewed by 554
Abstract
It is essential to provide high-integrity navigation information for safety-critical applications. Global navigation satellite systems (GNSSs) play an important role in these applications because they can provide global, high-accuracy, all-weather navigation services. Therefore, it has been a hot topic to improve GNSS integrity [...] Read more.
It is essential to provide high-integrity navigation information for safety-critical applications. Global navigation satellite systems (GNSSs) play an important role in these applications because they can provide global, high-accuracy, all-weather navigation services. Therefore, it has been a hot topic to improve GNSS integrity performance. This paper focuses on an initial proposal of the next-generation BeiDou Navigation Satellite System (BDS) integrity subsystem, with the aim of providing high-quality and global integrity services for the BDS. This paper first reviews the current status of the third-generation BDS integrity service. Following this, this paper proposes a space–ground integrated integrity monitoring design for the BDS that integrates the traditional ground-based integrity monitoring method, the advanced satellite autonomous integrity monitoring (A-SAIM) method, and the augmentation from low-earth-orbit (LEO) satellites. Specifically, this work offers an initial design of the A-SAIM method, which considers both single-satellite autonomous integrity monitoring and multi-satellite joint integrity monitoring. In addition, this work describes two different ways to augment BDS integrity with LEO satellites, i.e., (a) LEO satellites act as space monitoring stations and (b) LEO satellites act as navigation satellites. Simulations are carried out to validate the proposed design using CAT-I operation in civil aviation as an example. Simulation results indicate the effectiveness of the proposed design. In addition, simulation results suggest that if the fault probability of LEO satellites is worse than 1 × 10−4, LEO satellites can contribute more to BDS integrity performance improvement by acting as space monitoring stations; otherwise, it would be better to employ LEO satellites to broadcast navigation signals. The results also suggest that after taking LEO satellites into account, the global coverage of CAT-I can be potentially improved from 67% to 99%. This work is beneficial to the design of the next-generation BDS integrity subsystem. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The integrity architecture of the BDS-3, comprising fundamental integrity services and SBAS integrity services.</p>
Full article ">Figure 2
<p>Architecture of satellite-based augmentation system (SBAS).</p>
Full article ">Figure 3
<p>The architecture of the integrity monitoring method of space–ground integration.</p>
Full article ">Figure 4
<p>Single-satellite autonomous integrity monitoring.</p>
Full article ">Figure 5
<p>Communication between LEO and BeiDou satellites with Beidou short messages.</p>
Full article ">Figure 6
<p>Communication between LEO and BeiDou satellites based on space–ground connections.</p>
Full article ">Figure 7
<p>LEO satellites augment navigation integrity by broadcasting navigation signals.</p>
Full article ">Figure 8
<p>Number of visible satellites (BDS-3 only).</p>
Full article ">Figure 9
<p>Number of visible satellites (LEO + BDS-3).</p>
Full article ">Figure 10
<p>PDOP world map for BDS-3.</p>
Full article ">Figure 11
<p>PDOP world map for BDS-3 and LEO satellites.</p>
Full article ">Figure 12
<p>Pseudocode for integrity performance evaluation under CAT-I.</p>
Full article ">Figure 13
<p>Protection level world map with BDS-3 and GPS.</p>
Full article ">Figure 14
<p>Protection level world map with BDS-3, GPS, and LEO satellites, while LEO satellites act as space monitoring stations.</p>
Full article ">Figure 15
<p>Protection level world map with BDS-3, GPS, and LEO satellites, while LEO satellites act as navigation satellites. (LEO satellite fault probability is 10<sup>−5</sup>; LEO constellation fault probability is 10<sup>−5</sup>).</p>
Full article ">Figure 16
<p>Protection level world map with BDS-3, GPS, and LEO satellites, while LEO satellites act as navigation satellites. (LEO satellite fault probability is 10<sup>−4</sup>; LEO constellation fault probability is 10<sup>−4</sup>).</p>
Full article ">Figure 17
<p>Protection level world map with BDS-3, GPS, and LEO satellites, while LEO satellites act as navigation satellites. (LEO satellite fault probability is 10<sup>−3</sup>; LEO constellation fault probability is 10<sup>−3</sup>).</p>
Full article ">Figure 18
<p>Protection level world map with BDS-3, GPS, and LEO satellites, while LEO satellites act as navigation satellites. (LEO satellite fault probability is 10<sup>−3</sup>; LEO constellation fault probability is 10<sup>−2</sup>).</p>
Full article ">Figure 19
<p>Global average vertical protection level comparison.</p>
Full article ">Figure 20
<p>Global average horizontal protection level comparison.</p>
Full article ">Figure 21
<p>Global coverage rate comparison under Category I precision approach (CAT-I).</p>
Full article ">
18 pages, 6209 KiB  
Article
Impact of Latency and Continuity of GNSS Products on Filter-Based Real-Time LEO Satellite Clock Determination
by Meifang Wu, Kan Wang, Jinqian Wang, Wei Xie, Jiawei Liu, Beixi Chen, Yulong Ge, Ahmed El-Mowafy and Xuhai Yang
Remote Sens. 2024, 16(22), 4315; https://doi.org/10.3390/rs16224315 - 19 Nov 2024
Viewed by 500
Abstract
High-precision Low Earth Orbit (LEO) satellite clocks are essential for LEO-augmented Positioning, Navigation, and Timing (PNT) services. Nowadays, high-precision LEO satellite clocks can be determined in real-time using a Kalman filter either onboard or on the ground, as long as the GNSS observations [...] Read more.
High-precision Low Earth Orbit (LEO) satellite clocks are essential for LEO-augmented Positioning, Navigation, and Timing (PNT) services. Nowadays, high-precision LEO satellite clocks can be determined in real-time using a Kalman filter either onboard or on the ground, as long as the GNSS observations collected onboard LEO satellites can be transmitted to the ground in real-time. While various real-time and high-precision GNSS products are available nowadays in the latter case, their continuity and latencies in engineering reality are not as perfect as expected and will lead to unignorable impacts on the precision of the real-time LEO satellite clocks. In this study, based on real observations of Sentinel-3B, the impacts of different latencies and continuity of the real-time GNSS products on LEO real-time clocks are determined and discussed for two scenarios, namely the “epoch estimation” and “arc estimation” scenarios. The former case refers to the traditional filter-based processing epoch-by-epoch, and the latter case connects LEO satellite clocks from different rounds of filter-based processing under a certain arc length. The two scenarios lead to the “end-loss” and “mid-gap” situations. Latencies of the real-time GNSS products are discussed for the cases of orbit-only latency, clock-only latency, and combined forms, and different handling methods for the missing GNSS satellite clocks are discussed and compared. Results show that the real-time LEO satellite clock precision is very sensitive to the precision of real-time GNSS satellite clocks, and prediction of the latter becomes essential in case of their latencies. For the “end-loss” situation, with a latency of 30 to 120 s for the GNSS real-time clocks, the LEO satellite clock precision is reduced from about 0.2 to 0.28–0.57 ns. Waiting for the GNSS products in case of their short latencies and predicting the LEO satellite clocks instead could be a better option. For “arc-estimation”, when the gap of GNSS real-time products increases from 5 to 60 min, the real-time LEO clock precision decreases from 0.26 to 0.32 ns. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The continuity of the received real-time GPS satellite clocks from four analysis centers.</p>
Full article ">Figure 2
<p>Epoch-processing (<b>left</b>) and arc-processing (<b>right</b>) of real-time LEO satellite clock determination based on the kinematic model using the Kalman filter.</p>
Full article ">Figure 3
<p>Processing timelines of the two methods for real-time LEO satellite clocks determination based on the kinematic model and using the Kalman filter.</p>
Full article ">Figure 4
<p>(<b>Top</b>): Sentinel-3B satellite clocks determined using the Kalman filter-based kinematic model with the CNES real-time products without latency (blue line) and using the BLS-based reduced-dynamic model with the CODE final products (red line); (<b>Bottom</b>): differences between the two sets of clocks (real-time and final) in the top panel after unifying the time reference.</p>
Full article ">Figure 5
<p>Sentinel-3B clocks determined using the Kalman filter-based kinematic model and CNES real-time products with different latencies for Scenario A. The clocks using CNES real-time products without loss at the end are used as the reference.</p>
Full article ">Figure 6
<p>Sentinel-3B clocks determined using the Kalman filter-based kinematic model and CNES real-time products with different latencies for Scenario A. The BLS clocks using the CODE final products are used as the reference.</p>
Full article ">Figure 7
<p>Sentinel-3B clocks in the five scenarios (see <a href="#remotesensing-16-04315-t003" class="html-table">Table 3</a>) with a latency of 120 s compared with the clocks estimated using the CNES real-time products without latencies.</p>
Full article ">Figure 8
<p>Sentinel-3B clocks in the five scenarios (with a latency of 120 s, see <a href="#remotesensing-16-04315-t003" class="html-table">Table 3</a>) and the clocks without latencies in the CNES products, compared with the BLS clocks using the CODE final products.</p>
Full article ">Figure 9
<p>The MDEVs of the satellite clocks (<b>left</b>) and the clock estimation errors and prediction errors over 120 s (<b>right</b>) for G16 and G30.</p>
Full article ">Figure 10
<p>Sentinel-3B satellite clocks determined using the Kalman filter-based kinematic model with the CNES real-time products without latency (black line); using the BLS-based reduced-dynamic model with the CODE final products (red line); the connected near-real-time clocks estimated in a Kalman filter using the CNES real-time products without gaps (blue line).</p>
Full article ">Figure 11
<p>Differences between the C and K clocks (<b>top</b>) and C and R clocks (<b>bottom</b>) for Sentinel-3B. The explanations of C, K, and R clocks are given in <a href="#remotesensing-16-04315-f010" class="html-fig">Figure 10</a>.</p>
Full article ">Figure 12
<p>Near-real-time Sentinel-3B clocks based on a Kalman filter-based kinematic model, using CNES products with different “mid-gap” durations. (<b>Top</b>): Sentinel-3B clocks for the entire day; (<b>Bottom</b>): Sentinel-3B clocks from the first epoch affected to the last epoch affected.</p>
Full article ">Figure 13
<p>Near-real-time Sentinel-3B clock errors using the Kalman filter-based kinematic model with the CNES real-time products suffering from different “mid-gap” durations, compared with those using continuous CNES real-time products. (<b>Top</b>): Sentinel-3B clocks for the entire day; (<b>Bottom</b>): Sentinel-3B clocks from the first epoch affected to the last epoch affected.</p>
Full article ">Figure 14
<p>Near-real-time Sentinel-3B clock errors determined using the Kalman filter-based kinematic model with CNES products of different “mid-gap” durations, compared with the BLS clocks determined using the CODE final products. (<b>Top</b>): Sentinel-3B clocks for the entire day; (<b>Bottom</b>): Sentinel-3B clocks from the first epoch affected to the last epoch affected.</p>
Full article ">
15 pages, 5216 KiB  
Article
Deep Learning Evidence for Global Optimality of Gerver’s Sofa
by Kuangdai Leng, Jia Bi, Jaehoon Cha, Samuel Pinilla and Jeyan Thiyagalingam
Symmetry 2024, 16(10), 1388; https://doi.org/10.3390/sym16101388 - 18 Oct 2024
Viewed by 1307
Abstract
The moving sofa problem, introduced by Leo Moser in 1966, seeks to determine the maximal area of a 2D shape that can navigate an L-shaped corridor of unit width. Joseph Gerver’s 1992 solution, providing a lower bound of approximately 2.2195, is the [...] Read more.
The moving sofa problem, introduced by Leo Moser in 1966, seeks to determine the maximal area of a 2D shape that can navigate an L-shaped corridor of unit width. Joseph Gerver’s 1992 solution, providing a lower bound of approximately 2.2195, is the best known, though its global optimality remains unproven. This paper leverages neural networks’ approximation power and recent advances in invexity optimization to explore global optimality. We propose two approaches supporting Gerver’s conjecture that his sofa is the unique global maximum. The first approach uses continuous function learning, discarding assumptions about the monotonicity, symmetry, and differentiability of sofa movements. The sofa area is computed as a differentiable function using our “waterfall” algorithm, with the loss function incorporating both differential terms and initial conditions based on physics-informed machine learning. Extensive training with diverse network initialization consistently converges to Gerver’s solution. The second approach applies discrete optimization to the Kallus–Romik upper bound, improving it from 2.37 to 2.3337 for five rotation angles. As the number of angles increases, our model asymptotically converges to Gerver’s area from above, indicating that no larger sofa exists. Full article
(This article belongs to the Section Computer)
Show Figures

Figure 1

Figure 1
<p>The moving sofa problem proposed by Leo Moser [<a href="#B1-symmetry-16-01388" class="html-bibr">1</a>] in 1966 and some lower bounds. The best-known two are found, respectively, by John Hammersley [<a href="#B2-symmetry-16-01388" class="html-bibr">2</a>] and Joseph Gerver [<a href="#B3-symmetry-16-01388" class="html-bibr">3</a>]. The point markers separate the sections formed by different contact mechanisms (i.e., with either the inner corner, the walls, or the wall envelopes of the corridor).</p>
Full article ">Figure 2
<p>Geometry of the moving sofa problem. The movement of the corridor is described by the trajectory of its inner corner <span class="html-italic">P</span>, as denoted by <span class="html-italic">p</span>, and the rotation angle <math display="inline"><semantics> <mi>α</mi> </semantics></math>. The four walls form the four families of lines: <math display="inline"><semantics> <msub> <mi>l</mi> <mi>ih</mi> </msub> </semantics></math>, <math display="inline"><semantics> <msub> <mi>l</mi> <mi>iv</mi> </msub> </semantics></math>, <math display="inline"><semantics> <msub> <mi>l</mi> <mi>oh</mi> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>l</mi> <mi>ov</mi> </msub> </semantics></math>, with the subscripts showing their initial positions (i for inner, o for outer, h for horizontal and v for vertical). Their envelopes are, respectively, <math display="inline"><semantics> <msub> <mi>e</mi> <mi>ih</mi> </msub> </semantics></math>, <math display="inline"><semantics> <msub> <mi>e</mi> <mi>iv</mi> </msub> </semantics></math>, <math display="inline"><semantics> <msub> <mi>e</mi> <mi>oh</mi> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>e</mi> <mi>ov</mi> </msub> </semantics></math>. In this example, <span class="html-italic">p</span> is a semi-ellipse with its major and minor lengths being, respectively, 1.8 and 1.1, leading to some complexities in the envelopes, as highlighted by the magnified windows.</p>
Full article ">Figure 3
<p>Wall envelopes and sofa shape generated by an irregular corridor movement (red curve). Refer to <a href="#symmetry-16-01388-f002" class="html-fig">Figure 2</a> for detailed descriptions of the geometric elements. The envelope’s complex geometry requires a robust algorithm to ensure precise area calculation.</p>
Full article ">Figure 4
<p>The <span class="html-italic">waterfall</span> algorithm for area calculation. Left: The <span class="html-italic">waterfall</span> is tested on two hand-drawn curves from both below and above. Right: The <span class="html-italic">waterfall</span> is applied to the curves formed by the corridor movement in <a href="#symmetry-16-01388-f002" class="html-fig">Figure 2</a>. Water sources are placed along the horizontal line in the middle, with the arrows indicating the falling direction.</p>
Full article ">Figure 5
<p>Sets defined for Kallus–Romik upper bound. In this example, the angles are <math display="inline"><semantics> <mrow> <msub> <mi>α</mi> <mn>1</mn> </msub> <mo>≈</mo> <mn>16.26</mn> <mo>°</mo> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>α</mi> <mn>2</mn> </msub> <mo>≈</mo> <mn>30.51</mn> <mo>°</mo> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>α</mi> <mn>3</mn> </msub> <mo>≈</mo> <mn>44.76</mn> <mo>°</mo> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>β</mi> <mn>1</mn> </msub> <mo>≈</mo> <mn>73.74</mn> <mo>°</mo> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi>β</mi> <mn>2</mn> </msub> <mo>≈</mo> <mn>79.61</mn> <mo>°</mo> </mrow> </semantics></math>, corresponding to Equation (28) in [<a href="#B6-symmetry-16-01388" class="html-bibr">6</a>].</p>
Full article ">Figure 6
<p>Function initialization with network weights sampled from scaled uniform distributions. Each function is plotted in a different color to distinguish the individual curves. We sample weights and biases from <math display="inline"><semantics> <mrow> <mi mathvariant="script">U</mi> <mfenced separators="" open="(" close=")"> <mo>−</mo> <mi>s</mi> <msqrt> <mi>k</mi> </msqrt> <mo>,</mo> <mi>s</mi> <msqrt> <mi>k</mi> </msqrt> </mfenced> </mrow> </semantics></math>, where <span class="html-italic">k</span> is the reciprocal of input size and <span class="html-italic">s</span> a number we obtain for each target function by trial and error until its admissible range is safely covered from above, as indicated by the dashed lines. Note that <math display="inline"><semantics> <mrow> <mi>s</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math> is PyTorch default.</p>
Full article ">Figure 7
<p>Schematic of our physics-informed network architecture and loss function. The <math display="inline"><semantics> <mi mathvariant="script">F</mi> </semantics></math>’s represent unconstrained ReLU-based FCNs. Once <math display="inline"><semantics> <mrow> <msub> <mi>x</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>y</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math> are computed, automatic differentiation (AD) is applied to obtain their derivatives, which are then used in Equation (<a href="#FD4-symmetry-16-01388" class="html-disp-formula">4</a>) to calculate the envelopes, and subsequently, the sofa area via our <span class="html-italic">waterfall</span> algorithm.</p>
Full article ">Figure 8
<p>Landscape of sofa area. In this visualized model, the FCNs have three hidden layers of size 128, trained with a learning rate of <math display="inline"><semantics> <msup> <mn>10</mn> <mrow> <mo>−</mo> <mn>4</mn> </mrow> </msup> </semantics></math>. The landscape is anchored by the dominant eigenvectors of the Hessian matrix with respect to model weights [<a href="#B32-symmetry-16-01388" class="html-bibr">32</a>,<a href="#B33-symmetry-16-01388" class="html-bibr">33</a>], normalized to unit length. The range of the plot is <math display="inline"><semantics> <mrow> <mo>[</mo> <mo>−</mo> <mn>1.5</mn> <mo>,</mo> <mn>1.5</mn> <mo>]</mo> <mo>×</mo> <mo>[</mo> <mo>−</mo> <mn>1.5</mn> <mo>,</mo> <mn>1.5</mn> <mo>]</mo> </mrow> </semantics></math>, centered at Gerver’s area (the summit). There are no other peaks in the landscape across <math display="inline"><semantics> <mrow> <mfenced separators="" open="[" close="]"> <mo>−</mo> <mn>10</mn> <mo>,</mo> <mn>10</mn> </mfenced> <mo>×</mo> <mfenced separators="" open="[" close="]"> <mo>−</mo> <mn>10</mn> <mo>,</mo> <mn>10</mn> </mfenced> </mrow> </semantics></math>. Point <span class="html-italic">G</span> represents Gerver’s solution, the global maximum. Points <span class="html-italic">B</span> and <span class="html-italic">C</span> are local maxima, marked in each viewport to help associate perspectives across the different views.</p>
Full article ">Figure 9
<p>Convergence of <math display="inline"><semantics> <msubsup> <mi>G</mi> <mrow> <msub> <mi>γ</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>γ</mi> <mn>2</mn> </msub> <mo>,</mo> <mo>⋯</mo> <mo>,</mo> <msub> <mi>γ</mi> <mrow> <mi>n</mi> <mo>−</mo> <mn>2</mn> </mrow> </msub> </mrow> <mrow> <msub> <mi>γ</mi> <mrow> <mi>n</mi> <mo>−</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>γ</mi> <mi>n</mi> </msub> </mrow> </msubsup> </semantics></math> to Gerver’s area, where <math display="inline"><semantics> <mrow> <msub> <mi>γ</mi> <mi>j</mi> </msub> <mo>=</mo> <mstyle scriptlevel="0" displaystyle="true"> <mfrac> <mi>j</mi> <mi>n</mi> </mfrac> </mstyle> <mstyle scriptlevel="0" displaystyle="true"> <mfrac> <mi>π</mi> <mn>2</mn> </mfrac> </mstyle> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>⋯</mo> <mo>,</mo> <mi>n</mi> </mrow> </semantics></math>. (<b>Left</b>) <math display="inline"><semantics> <msubsup> <mi>G</mi> <mrow> <msub> <mi>γ</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>γ</mi> <mn>2</mn> </msub> <mo>,</mo> <mo>⋯</mo> <mo>,</mo> <msub> <mi>γ</mi> <mrow> <mi>n</mi> <mo>−</mo> <mn>2</mn> </mrow> </msub> </mrow> <mrow> <msub> <mi>γ</mi> <mrow> <mi>n</mi> <mo>−</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>γ</mi> <mi>n</mi> </msub> </mrow> </msubsup> </semantics></math> in linear scale. (<b>Right</b>) <math display="inline"><semantics> <mrow> <msubsup> <mi>G</mi> <mrow> <msub> <mi>γ</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>γ</mi> <mn>2</mn> </msub> <mo>,</mo> <mo>⋯</mo> <mo>,</mo> <msub> <mi>γ</mi> <mrow> <mi>n</mi> <mo>−</mo> <mn>2</mn> </mrow> </msub> </mrow> <mrow> <msub> <mi>γ</mi> <mrow> <mi>n</mi> <mo>−</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>γ</mi> <mi>n</mi> </msub> </mrow> </msubsup> <mo>−</mo> <mn>2.219532</mn> </mrow> </semantics></math> in logarithmic scale.</p>
Full article ">Figure 10
<p>Set determining <math display="inline"><semantics> <msub> <mi>G</mi> <mrow> <msub> <mi>α</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>α</mi> <mn>2</mn> </msub> <mo>,</mo> <msub> <mi>α</mi> <mn>3</mn> </msub> <mo>,</mo> <msub> <mi>α</mi> <mn>4</mn> </msub> <mo>,</mo> <msub> <mi>α</mi> <mn>5</mn> </msub> </mrow> </msub> </semantics></math> (<math display="inline"><semantics> <mrow> <mo>≈</mo> <mn>2.3337</mn> </mrow> </semantics></math>) and landscape of <math display="inline"><semantics> <msub> <mi>g</mi> <mrow> <msub> <mi>α</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>α</mi> <mn>2</mn> </msub> <mo>,</mo> <msub> <mi>α</mi> <mn>3</mn> </msub> <mo>,</mo> <msub> <mi>α</mi> <mn>4</mn> </msub> <mo>,</mo> <msub> <mi>α</mi> <mn>5</mn> </msub> </mrow> </msub> </semantics></math>. The landscape uses the rotation centers <math display="inline"><semantics> <mfenced separators="" open="(" close=")"> <msub> <mi mathvariant="bold-italic">u</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi mathvariant="bold-italic">u</mi> <mn>2</mn> </msub> <mo>,</mo> <mo>⋯</mo> <mo>,</mo> <msub> <mi mathvariant="bold-italic">u</mi> <mn>5</mn> </msub> </mfenced> </semantics></math> (instead of the network weights) as the high-dimensional variable for directional perturbation, based on the discrete nature of the Kallus–Romik upper bound (namely, the weights are less relevant since the <math display="inline"><semantics> <mi>α</mi> </semantics></math>-sequence is prescribed). The two surface plots show the same peak from different viewpoints, with a range of <math display="inline"><semantics> <mrow> <mfenced separators="" open="[" close="]"> <mo>−</mo> <mn>1</mn> <mo>,</mo> <mn>1</mn> </mfenced> <mo>×</mo> <mfenced separators="" open="[" close="]"> <mo>−</mo> <mn>1</mn> <mo>,</mo> <mn>1</mn> </mfenced> </mrow> </semantics></math> centered at <math display="inline"><semantics> <msub> <mi>G</mi> <mrow> <msub> <mi>α</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>α</mi> <mn>2</mn> </msub> <mo>,</mo> <msub> <mi>α</mi> <mn>3</mn> </msub> <mo>,</mo> <msub> <mi>α</mi> <mn>4</mn> </msub> <mo>,</mo> <msub> <mi>α</mi> <mn>5</mn> </msub> </mrow> </msub> </semantics></math> (the summit). There are no other peaks in the landscape across <math display="inline"><semantics> <mrow> <mfenced separators="" open="[" close="]"> <mo>−</mo> <mn>10</mn> <mo>,</mo> <mn>10</mn> </mfenced> <mo>×</mo> <mfenced separators="" open="[" close="]"> <mo>−</mo> <mn>10</mn> <mo>,</mo> <mn>10</mn> </mfenced> </mrow> </semantics></math>.</p>
Full article ">
17 pages, 5448 KiB  
Article
Orbit Determination Method for BDS-3 MEO Satellites Based on Multi-Source Observation Links
by Jie Xin and Kai Li
Remote Sens. 2024, 16(19), 3702; https://doi.org/10.3390/rs16193702 - 4 Oct 2024
Viewed by 722
Abstract
Research on augmentation and supplement systems for navigation systems has become a significant aspect in comprehensive positioning, navigation and timing (PNT) studies. The BeiDou-3 navigation satellite system (BDS-3) has constructed a dynamic inter-satellite network to gain more observation data than ground monitoring stations. [...] Read more.
Research on augmentation and supplement systems for navigation systems has become a significant aspect in comprehensive positioning, navigation and timing (PNT) studies. The BeiDou-3 navigation satellite system (BDS-3) has constructed a dynamic inter-satellite network to gain more observation data than ground monitoring stations. Low Earth orbit (LEO) satellites have advantages in their kinematic velocity and information carrying rate and can be used as satellite-based monitoring stations for navigation satellites to make up for the distribution limitation of ground monitoring stations. This study constructs multi-source observation links with satellite-to-ground, inter-satellite and satellite-based observation data, proposes an orbit synchronization method for navigation satellites and LEO satellites and verifies the influence thereof on orbit accuracy with different observation data. The experimental results under conditions of real and simulated observation data showed the following: (1) With the support of satellite-based observation links, the orbit accuracy of the BDS-3 MEO satellites could be improved significantly, with a 78% improvement with the simulation data and a 76% improvement with the real data. When the navigation satellites leave the monitoring area of the ground monitoring stations, the accuracy reduction tendency of the orbit prediction could also be slowed down with the support of the LEO satellites and the accuracy could be maintained within centimeters. (2) Comparing the orbit accuracy with the support of the satellite-to-ground observation links, the orbit accuracy of the MEO satellites could be improved by 65.5%, 73.7% and 79.4% with the support of the 6, 12 and 60 LEO satellites, respectively. When the observation geometry and the covering multiplicity meet the basic requirement of orbit determination, the improvements to the orbit accuracy decrease with the growth of LEO satellite numbers. (3) The accuracy of orbit determination with the support of the LEO satellites or the inter-satellite links was at the centimeter level for both, verifying that inter-satellite links and satellite-based links can be used as each other’s backups for navigation satellites. (4) The accuracy of orbit determination with the multi-source observation links was also at the centimeter level, which was not better than the results with the support of the satellite-to-ground and inter-satellite links or the satellite-to-ground and satellite-based links. Full article
Show Figures

Figure 1

Figure 1
<p>Multi-source observation network.</p>
Full article ">Figure 2
<p>Processing chart of multi-source observation equations.</p>
Full article ">Figure 3
<p>Calculating chart for the normal equations.</p>
Full article ">Figure 4
<p>Designed constellation with BDS-3 MEO satellites and 60 LEO satellites. (The satellites in blue, yellow and pink lines are BDS-3 MEO satellites and the satellites in low orbit lines are LEO satellites).</p>
Full article ">Figure 5
<p>LEO constellation with 5 LEO satellites.</p>
Full article ">Figure 6
<p>Distribution of the chosen stations (simulated data).</p>
Full article ">Figure 7
<p>Distribution of the chosen stations (real data).</p>
Full article ">Figure 8
<p>Number of satellite-to-ground observation links.</p>
Full article ">Figure 9
<p>Number of satellite-to-ground and satellite-based observation links, with the support of 6 LEO satellites.</p>
Full article ">Figure 10
<p>Number of satellite-to-ground and satellite-based observation links, with the support of 12 LEO satellites.</p>
Full article ">Figure 11
<p>Compared results of orbit determination under conditions of different elevation angles in Scenario 1 and Scenario 2 (the colorful lines are the observation angles related to different monitoring stations; the blue dots and the red dots are, respectively, the results for orbit determination with different observation links in Scenario 1 and Scenario 2).</p>
Full article ">Figure 12
<p>Compared results of orbit determination under conditions of different elevation angles in Scenario 1 and Scenario 3.</p>
Full article ">Figure 13
<p>Statistical results of the three-dimensional position errors in Scenarios 1 to 7 (one day).</p>
Full article ">Figure 14
<p>Number of satellite-to-ground and satellite-based receivers.</p>
Full article ">Figure 15
<p>Compared results of orbit determination under conditions of different elevation angles in Scenario 1 and Scenario 2 (real data).</p>
Full article ">
21 pages, 6869 KiB  
Article
Research on Design and Staged Deployment of LEO Navigation Constellation for MEO Navigation Satellite Failure
by Wen Xue, Min Hu, Yongjing Ruan, Xun Wang and Moyao Yu
Remote Sens. 2024, 16(19), 3667; https://doi.org/10.3390/rs16193667 - 1 Oct 2024
Viewed by 1035
Abstract
Low Earth orbit (LEO) satellites have unique advantages in navigation because of their high signal intensity and rapid geometric changes in a short period. In order to solve the problem of constellation performance degradation after a potential failure pertaining one or more medium [...] Read more.
Low Earth orbit (LEO) satellites have unique advantages in navigation because of their high signal intensity and rapid geometric changes in a short period. In order to solve the problem of constellation performance degradation after a potential failure pertaining one or more medium Earth orbit (MEO) navigation satellites, this paper designs the LEO navigation constellation and considers the task requirements of different stages of constellation deployment. Firstly, the LEO navigation constellation is designed by a non-dominated sorting genetic algorithm II (NSGA-II). The average position dilution of precision (PDOP) is 1.676, which is an improvement compared to the average PDOP offered by the four traditional GNSS. Secondly, the staged deployment of constellation takes into account the degradation of constellation performance caused by the failure of MEO navigation satellites, and the Monte Carlo method is used to analyze the case of three simultaneous satellite failures. The results show that a single satellite failure within each orbital plane and adjacent satellites with close phase separation has a great impact on the performance of the MEO navigation constellation. On this basis, a staged deployment strategy was adopted in order to balance cost, risk, and performance. The three phases deploy 66, 156, and 288 satellites, respectively; as a make-up constellation under contingencies, a navigation enhancement constellation, and an independent navigation constellation, the deployment of the staged sub-constellations meets the mission requirements. The constellation design and staged deployment method proposed in this paper can provide reference for the future study of LEO navigation constellations. Full article
(This article belongs to the Special Issue Advances in the Study of Intelligent Aerospace)
Show Figures

Figure 1

Figure 1
<p>Diagram of satellite orbital plane.</p>
Full article ">Figure 2
<p>Constellation optimization design model based on NSGA-II algorithm.</p>
Full article ">Figure 3
<p>Analysis of three-satellite failure based on Monte Carlo method.</p>
Full article ">Figure 4
<p>The staged deployment strategy of constellations.</p>
Full article ">Figure 5
<p>Analysis of altitude parameters for orbital constellation: (<b>a</b>) Distribution of ionospheric, Van Allen radiation belts, and the satellite orbital altitudes; (<b>b</b>) the relationship between satellite orbital altitude and inclination and the period of regression period.</p>
Full article ">Figure 6
<p>Average number of visible satellites.</p>
Full article ">Figure 7
<p>Optimization results of LEO tilt constellation: (<b>a</b>) optimization results and (<b>b</b>) evaluation index.</p>
Full article ">Figure 8
<p>Inclination optimization range of polar constellation.</p>
Full article ">Figure 9
<p>Optimization results of LEO polar constellation: (<b>a</b>) optimization results and (<b>b</b>) evaluation index.</p>
Full article ">Figure 10
<p>Variation in CV following the failure of a single and double satellite: (<b>a</b>) single-satellite failure and (<b>b</b>) double-satellite failure.</p>
Full article ">Figure 11
<p>Simulation results of three-satellite failure: (<b>a</b>) proportion of different failure modes and (<b>b</b>) simulation times and CV.</p>
Full article ">Figure 12
<p>Ten modes of three-satellite failure.</p>
Full article ">Figure 13
<p>Staged deployment spatial configuration of the LEO hybrid navigation constellation. (<b>a</b>) Stage I compensation. (<b>b</b>) Stage II global navigation augmentation. (<b>c</b>) Stage III global BDS-3 level.</p>
Full article ">Figure 14
<p>Five optimization schemes are compared with the PDOP of GNSS constellations: (<b>a</b>) LEO tilt constellation scheme and (<b>b</b>) LEO hybrid constellation scheme.</p>
Full article ">Figure 15
<p>The space phase of a single failure within each orbital plane: (<b>a</b>) failure modes of 14-22-35; (<b>b</b>) failure modes of 12-28-33; (<b>c</b>) failure modes of 11-24-32; (<b>d</b>) failure modes of 15-23-31.</p>
Full article ">Figure 16
<p>Staged deployment of PDOP in various scenarios for LEO hybrid navigation constellations.</p>
Full article ">
13 pages, 2715 KiB  
Technical Note
Laser Observations of GALILEO Satellites at the CBK PAN Astrogeodynamic Observatory in Borowiec
by Paweł Lejba, Piotr Michałek, Tomasz Suchodolski, Adrian Smagło, Mateusz Matyszewski and Stanisław Zapaśnik
Remote Sens. 2024, 16(15), 2862; https://doi.org/10.3390/rs16152862 - 5 Aug 2024
Viewed by 878
Abstract
The laser station (BORL) owned by the Space Research Centre of the Polish Academy of Sciences and situated at the Astrogeodynamic Observatory in Borowiec near Poznań regularly observes more than 100 different objects in low Earth orbit (LEO) and medium Earth orbit (MEO). [...] Read more.
The laser station (BORL) owned by the Space Research Centre of the Polish Academy of Sciences and situated at the Astrogeodynamic Observatory in Borowiec near Poznań regularly observes more than 100 different objects in low Earth orbit (LEO) and medium Earth orbit (MEO). The BORL sensor’s laser observation range is from 400 km to 24,500 km. The laser measurements taken by the BORL sensor are utilized to create various products, including the geocentric positions and movements of ground stations, satellite orbits, the components of the Earth’s gravitational field and their changes over time, Earth’s orientation parameters (EOPs), and the validation of the precise Galileo orbits derived using microwave measurements, among others. These products are essential for supporting local and global geodetic and geophysics research related to time. They are crucial for the International Terrestrial Reference Frame (ITRF), which is managed by the International Earth Rotation and Reference Systems Service (IERS). In 2023, the BORL laser station expanded its list of tracked objects to include all satellites of the European satellite navigation system GALILEO, totaling 28 satellites. During that year, the BORL laser station recorded 77 successful passes of GALILEO satellites, covering a total of 21 objects. The measurements taken allowed for the registration of 7419 returns, resulting in 342 normal points. The average RMS for all successful GALILEO observations in 2023 was 13.5 mm. Full article
(This article belongs to the Special Issue Advanced Remote Sensing Technology in Modern Geodesy)
Show Figures

Figure 1

Figure 1
<p>The CBK PAN laser sensor in Borowiec consists of two independent laser modules: a smaller, nanosecond one in the foreground and a larger, picosecond one in the background.</p>
Full article ">Figure 2
<p>Number of tracked objects by BORL station in years 1992–2023.</p>
Full article ">Figure 3
<p>The constellation of GALILEO satellites [<a href="#B17-remotesensing-16-02862" class="html-bibr">17</a>].</p>
Full article ">Figure 4
<p>Artist’s view of a GALILEO Full Operational Capability (FOC). The green circle marks the panel with laser retroreflectors [<a href="#B18-remotesensing-16-02862" class="html-bibr">18</a>].</p>
Full article ">Figure 5
<p>Pass of GALILEO-212 (E03) registered by the BORL sensor on 3 May 2023, 23:12 UTC, showing slant ranges (<b>a</b>) and fit residuals (<b>b</b>) with 159 valid returns. The min. and max. distances to the satellite were 23,589 km and 23,790 km, respectively.</p>
Full article ">Figure 6
<p>Pass of GALILEO-104 (E20) registered by the BORL sensor on 21 May 2023, 21:27 UTC, showing slant ranges (<b>a</b>) and fit residuals (<b>b</b>) with 592 valid returns. The min. and max. distances to the satellite were 23,701 km and 24,115 km, respectively.</p>
Full article ">Figure 7
<p>Pass of GALILEO-208 (E08) registered by the BORL sensor on 22 May 2023, 20:59 UTC, showing slant ranges (<b>a</b>) and fit residuals (<b>b</b>) with 214 valid returns. The min. and max. distances to the satellite were 23,612 km and 23,958 km, respectively.</p>
Full article ">Figure 8
<p>Pass of GALILEO-202 (E14) registered by the BORL sensor on 21 September 2023, 22:37 UTC, showing slant ranges (<b>a</b>) and fit residuals (<b>b</b>) with 310 valid returns. The min. and max. distances to the satellite were 19,701 km and 20,292 km, respectively.</p>
Full article ">
23 pages, 6378 KiB  
Article
Navigation Resource Allocation Algorithm for LEO Constellations Based on Dynamic Programming
by Sixin Wang, Xiaomei Tang, Jingyuan Li, Xinming Huang, Jiyang Liu and Jian Liu
Remote Sens. 2024, 16(12), 2231; https://doi.org/10.3390/rs16122231 - 19 Jun 2024
Viewed by 852
Abstract
Navigation resource allocation for low-earth-orbit (LEO) constellations refers to the optimal allocation of navigational assets when the number and allocation of satellites in the LEO constellation have been determined. LEO constellations can not only transmit navigation enhancement signals but also enable space-based monitoring [...] Read more.
Navigation resource allocation for low-earth-orbit (LEO) constellations refers to the optimal allocation of navigational assets when the number and allocation of satellites in the LEO constellation have been determined. LEO constellations can not only transmit navigation enhancement signals but also enable space-based monitoring (SBM) for real-time assessment of GNSS signal quality. However, proximity in the frequencies of LEO navigation signals and SBM can lead to significant interference, necessitating isolated transmission and reception. This separation requires that SBM and navigation signal transmission be carried out by different satellites within the constellation, thus demanding a strategic allocation of satellite resources. Given the vast number of satellites and their rapid movement, the visibility among LEO, medium-earth-orbit (MEO), and geostationary orbit (GEO) satellites is highly dynamic, presenting substantial challenges in resource allocation due to the computational intensity involved. Therefore, this paper proposes an optimal allocation algorithm for LEO constellation navigation resources based on dynamic programming. In this algorithm, a network model for the allocation of navigation resources in LEO constellations is initially established. Under the constraints of visibility time windows and onboard transmission and reception isolation, the objective is set to minimize the number of LEO satellites used while achieving effective navigation signal transmission and SBM. The constraints of resource allocation and the mathematical expression of the optimization objective are derived. A dynamic programming approach is then employed to determine the optimal resource allocation scheme. Analytical results demonstrate that compared to Greedy and Divide-and-Conquer algorithms, this algorithm achieves the highest resource utilization rate and the lowest computational complexity, making it highly valuable for future resource allocation in LEO constellations. Full article
(This article belongs to the Special Issue Space-Geodetic Techniques (Third Edition))
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>LEO constellation navigation enhanced function architecture diagram.</p>
Full article ">Figure 2
<p>BDS constellation three-dimensional space allocation diagram. The red line represents the IGSO orbit, the green line represents the GEO orbit, and the blue line represents the MEO orbit.</p>
Full article ">Figure 3
<p>BDS constellation subsatellite point trajectory chart.</p>
Full article ">Figure 4
<p>Global average coverage of BDS.</p>
Full article ">Figure 5
<p>Global GDOP value distribution of BDS.</p>
Full article ">Figure 6
<p>GNSS satellites and LEO satellites.</p>
Full article ">Figure 7
<p>Dynamic programming algorithm schematic diagram.</p>
Full article ">Figure 8
<p>Flowchart of the NRAA-DP.</p>
Full article ">Figure 9
<p>The distribution diagram of the LEO constellation and BDS constellation. The blue lines represent the orbital planes of the BDS constellation, the yellow lines indicate the near-polar orbital planes of the LEO constellation, and the red lines denote the inclined orbital planes of the LEO constellation.</p>
Full article ">Figure 10
<p>Visual links between LEO satellite S0101 and BDS satellites within 4 h. Lines of different colors represent different visible links.</p>
Full article ">Figure 11
<p>Distribution map of ground stations used for evaluating constellation global coverage. The blue circles represent ground stations in the northern hemisphere along the 0° longitude from 0° to 90°N.</p>
Full article ">Figure 12
<p>Visible links between the LEO satellite I0101 and the ground stations. Lines of different colors represent different visible links.</p>
Full article ">Figure 13
<p>Constellation distribution map of the resource allocation scheme by NRAA-DP. (<b>a</b>) Inclined orbit satellites. (<b>b</b>) Near-polar orbit.</p>
Full article ">Figure 14
<p>Constellation distribution map of the resource allocation scheme by GA. (<b>a</b>) Inclined orbit satellites. (<b>b</b>) Near-polar orbit.</p>
Full article ">Figure 15
<p>Constellation distribution map of the resource allocation scheme by DCA. (<b>a</b>) Inclined orbit satellites. (<b>b</b>) Near-polar orbit.</p>
Full article ">Figure 16
<p>Coverage performance of different navigation resource allocation schemes for ground stations. (<b>a</b>) NRAA-DP. (<b>b</b>) GA. (<b>c</b>) DCA.</p>
Full article ">Figure 17
<p>Coverage performance of different navigation resource allocation schemes for BDS satellites. (<b>a</b>) NRAA-DP. (<b>b</b>) GA. (<b>c</b>) DCA.</p>
Full article ">
20 pages, 6559 KiB  
Article
Study of Fast and Reliable Time Transfer Methods Using Low Earth Orbit Enhancement
by Mingyue Liu, Rui Tu, Qiushi Chen, Qi Li, Junmei Chen, Pengfei Zhang and Xiaochun Lu
Remote Sens. 2024, 16(11), 2044; https://doi.org/10.3390/rs16112044 - 6 Jun 2024
Viewed by 742
Abstract
The Global Navigation Satellite System (GNSS) can be utilized for long-distance and high-precision time transmission. With the ongoing development of low Earth orbit (LEO) satellites and the rapidly changing geometric relationships between them, the convergence rate of ambiguity parameters in Precise Point Positioning [...] Read more.
The Global Navigation Satellite System (GNSS) can be utilized for long-distance and high-precision time transmission. With the ongoing development of low Earth orbit (LEO) satellites and the rapidly changing geometric relationships between them, the convergence rate of ambiguity parameters in Precise Point Positioning (PPP) algorithms has increased, enabling fast and reliable time transfer. In this paper, GPS is used as an experimental case, the LEO satellite constellation is designed, and simulated LEO observation data are generated. Then, using the GPS observation data provided by IGS, a LEO-enhanced PPP model is established. The LEO-augmented PPP model is employed to facilitate faster and more reliable high-precision time transfer. The application of the LEO-augmented PPP model to time transfer is examined and discussed through experimental examples. These examples show multiple types of time transfer links, and the experimental outcomes are uniform. GPS + LEO is compared with exclusive GPS time transfer schemes. The clock offset of the time transfer link for the GPS + LEO scheme converges more swiftly, meaning that the time required for the clock offset to reach a stable level is the briefest. In this paper, standard deviation is employed to assess stability, and Allan deviation is utilized to assess frequency stability. The results show that the clock offset stability and frequency stability achieved by the GPS + LEO scheme are superior within the convergence time range. Controlled experiments with different numbers of satellites for LEO enhancement indicate that time transfer performance can be improved by increasing the number of satellites. As a result, augmenting GPS tracking data with LEO observations enhances the time transfer service compared to GPS alone. Full article
(This article belongs to the Topic GNSS Measurement Technique in Aerial Navigation)
Show Figures

Figure 1

Figure 1
<p>Global distribution of stations and time transfer links.</p>
Full article ">Figure 2
<p>Average number of daily visible satellites of the LEO constellation on the first day of 2022.</p>
Full article ">Figure 3
<p>Receiver clock offset time series of stations (black denotes GPS alone; blue, GPS + LEOs.</p>
Full article ">Figure 4
<p>The first-order difference of the receiver clock offset time series of the stations (black denotes GPS alone; blue, GPS + LEOs). (<b>a</b>) The first-order difference results of the clock offset sequence for the areg station, (<b>b</b>) The first-order difference results of the clock offset sequence for the harb station, (<b>c</b>) The first-order difference results of the clock offset sequence for the ons1 station.</p>
Full article ">Figure 5
<p>Receiver clock offset convergence time.</p>
Full article ">Figure 6
<p>Time transfer link clock offset time series (black denotes GPS alone; blue GPS + LEOs). (<b>a</b>) The clock offset sequence of the areg-mcil time transfer link, (<b>b</b>) The clock offset sequence of the gold-pie1 time transfer link, (<b>c</b>) The clock offset sequence of the pie1-ons1 time transfer link, (<b>d</b>) The clock offset sequence of the ons1-bor1 time transfer link, (<b>e</b>) The clock offset sequence of the bor1-usud time transfer link, (<b>f</b>) The clock offset sequence of the kiru-dlf1 time transfer link, (<b>g</b>) The clock offset sequence of the dlf1-harb time transfer link, (<b>h</b>) The clock offset sequence of the harb-sydn time transfer link, (<b>i</b>) The clock offset sequence of the kiru-syog time transfer link, (<b>j</b>) The clock offset sequence of the ohi2-syog time transfer link.</p>
Full article ">Figure 7
<p>The first-order difference results of clock offsets time series in time transfer links (black denotes GPS alone; blue, GPS + LEOs). (<b>a</b>) The first-order difference results of the clock offset sequence for the areg-mcil time transfer link, (<b>b</b>) The first-order difference results of the clock offset sequence for the gold-pie1 time transfer link, (<b>c</b>) The first-order difference results of the clock offset sequence for the pie1-ons1 time transfer link, (<b>d</b>) The first-order difference results of the clock offset sequence for the ons1-bor1 time transfer link, (<b>e</b>) The first-order difference results of the clock offset sequence for the bor1-usud time transfer link, (<b>f</b>) The first-order difference results of the clock offset sequence for the kiru-dlf1 time transfer link, (<b>g</b>) The first-order difference results of the clock offset sequence for the dlf1-harb time transfer link, (<b>h</b>) The first-order difference results of the clock offset sequence for the harb-sydn time transfer link, (<b>i</b>)The first-order difference results of the clock offset sequence for the kiru-syog time transfer link, (<b>j</b>) The first-order difference results of the clock offset sequence for the ohi2-syog time transfer link.</p>
Full article ">Figure 8
<p>Time transfer link convergence time histogram.</p>
Full article ">Figure 9
<p>Histogram of standard deviation of time transfer link clock offsets in the convergence interval.</p>
Full article ">Figure 10
<p>The Allan deviation results of the time transfer link in the convergence interval.</p>
Full article ">Figure 11
<p>Clock offset diagram of pie1-ons1 time transfer link in the four enhancement scenarios.</p>
Full article ">Figure 12
<p>First-order difference diagram of pie1-ons1 time transfer link clock offset in the four enhancement scenarios.</p>
Full article ">Figure 13
<p>Clock offset diagram of harb-sydn time transfer link in the four enhancement scenarios.</p>
Full article ">Figure 14
<p>First-order difference diagram of harb-sydn time transfer link clock offset in the four enhancement scenarios.</p>
Full article ">Figure 15
<p>Clock offset diagram of ohi2-syog time transfer link in the four enhancement scenarios.</p>
Full article ">Figure 16
<p>First-order difference diagram of ohi2-syog time transfer link clock offset in the four enhancement scenarios.</p>
Full article ">Figure 17
<p>Convergence time histograms of the three links across the four enhancement scenarios.</p>
Full article ">
15 pages, 5146 KiB  
Article
A Fast Time Synchronization Method for Large Scale LEO Satellite Networks Based on A Bionic Algorithm
by Yue Xu, Tao Dong, Jie Yin, Ziyong Zhang, Zhihui Liu, Hao Jiang and Jing Wu
Photonics 2024, 11(5), 475; https://doi.org/10.3390/photonics11050475 - 19 May 2024
Viewed by 1282
Abstract
A fast time synchronization method for large-scale LEO satellite networks based on a bionic algorithm is proposed. Because the inter-satellite links are continuously established and interrupted due to the relative motion of the satellites, the topology of the LEO satellite networks is time [...] Read more.
A fast time synchronization method for large-scale LEO satellite networks based on a bionic algorithm is proposed. Because the inter-satellite links are continuously established and interrupted due to the relative motion of the satellites, the topology of the LEO satellite networks is time varying. Firstly, according to the ephemeris information in navigation messages, a connection table which records the connections between satellites is generated. Then, based on the connection table, the current satellite network topology is calculated and generated. Furthermore, a bionic algorithm is used to select some satellites as time source nodes and calculate the hierarchy of the clock transmission tree. By taking the minimum level of the time transmission tree as the optimization objective, the time source nodes and the clock stratums of the whole satellite networks are obtained. Finally, the onboard computational center broadcasts the time layer table to all the satellites in the LEO satellite networks and the time synchronization links can be established or recovered fast. Full article
(This article belongs to the Special Issue Novel Advances in Optical Communications)
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of space and earth integration networks.</p>
Full article ">Figure 2
<p>Time transfer trees: (<b>a</b>) the second group, (<b>b</b>) the fifth group in <a href="#photonics-11-00475-t001" class="html-table">Table 1</a>.</p>
Full article ">Figure 3
<p>The curves of the number of iterations of the PSO algorithm and the traversal algorithm with respect to the minimum sum of the weighted time in the first time slice.</p>
Full article ">Figure 4
<p>Schematic diagram of clock transmission path and topology of the first slice in LEO satellite networks.</p>
Full article ">Figure 5
<p>Network topology for another time slice.</p>
Full article ">Figure 6
<p>Time transfer trees of the fourth group in <a href="#photonics-11-00475-t004" class="html-table">Table 4</a>.</p>
Full article ">Figure 7
<p>The curves of the number of iterations of the PSO algorithm and the traversal algorithm with respect to the minimum sum of the weighted time in another time slice.</p>
Full article ">Figure A1
<p></p>
Full article ">Figure A2
<p></p>
Full article ">Figure A3
<p></p>
Full article ">
23 pages, 7491 KiB  
Article
LEO-Enhanced GNSS/INS Tightly Coupled Integration Based on Factor Graph Optimization in the Urban Environment
by Shixuan Zhang, Rui Tu, Zhouzheng Gao, Decai Zou, Siyao Wang and Xiaochun Lu
Remote Sens. 2024, 16(10), 1782; https://doi.org/10.3390/rs16101782 - 17 May 2024
Viewed by 1178
Abstract
Precision point positioning (PPP) utilizing the Global Navigation Satellite System (GNSS) is a traditional and widely employed technology. Its performance is susceptible to observation discontinuities and unfavorable geometric configurations. Consequently, the integration of the Inertial Navigation System (INS) and GNSS makes full use [...] Read more.
Precision point positioning (PPP) utilizing the Global Navigation Satellite System (GNSS) is a traditional and widely employed technology. Its performance is susceptible to observation discontinuities and unfavorable geometric configurations. Consequently, the integration of the Inertial Navigation System (INS) and GNSS makes full use of their respective advantages and effectively mitigates the limitations of GNSS positioning. However, the GNSS/INS integration faces significant challenges in complex and harsh urban environments. In recent years, the geometry between the user and the satellite has been effectively improved with the advent of lower-orbits and faster-speed Low Earth Orbit (LEO) satellites. This enhancement provides more observation data, opening up new possibilities and opportunities for high-precision positioning. Meanwhile, in contrast to the traditional extended Kalman filter (EKF) approach, the performance of the LEO-enhanced GNSS/INS tightly coupled integration (TCI) can be significantly improved by employing the factor graph optimization (FGO) method with multiple iterations to achieve stable estimation. In this study, LEO data and the FGO method were employed to enhance the GNSS/INS TCI. To validate the effectiveness of the method, vehicle data and simulated LEO observations were subjected to thorough analysis. The results suggest that the integration of LEO data significantly enhances the positioning accuracy and convergence speed of the GNSS/INS TCI. In contrast to the FGO GNSS/INS TCI without LEO enhancement, the average enhancement effect of the LEO is 22.16%, 7.58%, and 10.13% in the north, east, and vertical directions, respectively. Furthermore, the average root mean square error (RMSE) of the LEO-enhanced FGO GNSS/INS TCI is 0.63 m, 1.21 m, and 0.85 m in the north, east, and vertical directions, respectively, representing an average improvement of 41.91%, 13.66%, and 2.52% over the traditional EKF method. Meanwhile, the simulation results demonstrate that LEO data and the FGO method effectively enhance the positioning and convergence performance of GNSS/INS TCI in GNSS-challenged environments (tall buildings, viaducts, underground tunnels, and wooded areas). Full article
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of LEO-enhanced GNSS/INS TCI structure based on EKF.</p>
Full article ">Figure 2
<p>Schematic diagram of LEO-enhanced GNSS/INS TCI structure based on FGO.</p>
Full article ">Figure 3
<p>The trajectory of the land vehicle experiment in Beijing, China.</p>
Full article ">Figure 4
<p>Time series of the number of available satellites and PDOP for the GPS, GPS/LEO, BDS, and BDS/LEO.</p>
Full article ">Figure 5
<p>The position differences for the GPS/INS, LEO-enhanced GPS/INS, BDS/INS, and LEO-enhanced BDS/INS tightly coupled integration based on EKF.</p>
Full article ">Figure 6
<p>RMSE of position differences for the GPS/INS, LEO-enhanced GPS/INS, BDS/INS, and LEO-enhanced BDS/INS tightly coupled integration based on EKF.</p>
Full article ">Figure 7
<p>The position differences of the GPS/INS, LEO-enhanced GPS/INS, BDS/INS, and LEO-enhanced BDS/INS tightly coupled integration based on FGO.</p>
Full article ">Figure 8
<p>RMSE of position differences for the GPS/INS, LEO-enhanced GPS/INS, BDS/INS, and LEO-enhanced BDS/INS tightly coupled integration based on FGO.</p>
Full article ">Figure 9
<p>Cumulative distribution function of position difference for the GPS/INS, LEO-enhanced GPS/INS, BDS/INS, and LEO-enhanced BDS/INS tightly coupled integration based on FGO.</p>
Full article ">Figure 10
<p>Time series of the number of available satellites for the FGO GPS, GPS/LEO, BDS, and BDS/LEO under the low-observability environment.</p>
Full article ">Figure 11
<p>The position differences of the FGO GPS/INS, LEO-enhanced GPS/INS, BDS/INS, and LEO-enhanced BDS/INS tightly coupled integration under the low-observability environment.</p>
Full article ">Figure 12
<p>RMSE of position differences for the FGO GPS/INS, LEO-enhanced GPS/INS, BDS/INS, and LEO-enhanced BDS/INS tightly coupled integration under the low-observability environment.</p>
Full article ">Figure 13
<p>Convergence performance of the GPS/INS, LEO-enhanced GPS/INS, BDS/INS, and LEO-enhanced BDS/INS tightly coupled integration based on FGO.</p>
Full article ">Figure 14
<p>RMSE of position differences for the four tightly coupled integrations based on FGO.</p>
Full article ">
16 pages, 4218 KiB  
Article
A Study on Anti-Jamming Algorithms in Low-Earth-Orbit Satellite Signal-of-Opportunity Positioning Systems for Unmanned Aerial Vehicles
by Lihao Yao, Honglei Qin, Boyun Gu, Guangting Shi, Hai Sha, Mengli Wang, Deyong Xian, Feiqiang Chen and Zukun Lu
Drones 2024, 8(4), 164; https://doi.org/10.3390/drones8040164 - 20 Apr 2024
Cited by 1 | Viewed by 1799
Abstract
Low-Earth-Orbit (LEO) satellite Signal-of-Opportunity (SOP) positioning technology has gradually matured to meet the accuracy requirements for unmanned aerial vehicle (UAV) positioning in daily scenarios. Advancements in miniaturization technology for positioning terminals have also made this technology’s application to UAV positioning crucial for UAV [...] Read more.
Low-Earth-Orbit (LEO) satellite Signal-of-Opportunity (SOP) positioning technology has gradually matured to meet the accuracy requirements for unmanned aerial vehicle (UAV) positioning in daily scenarios. Advancements in miniaturization technology for positioning terminals have also made this technology’s application to UAV positioning crucial for UAV development. However, in the increasingly complex electromagnetic environment, there remains a significant risk of degradation in positioning performance for UAVs in LEO satellite SOP positioning due to unintentional or malicious jamming. Furthermore, there is a lack of in-depth research from scholars both domestically and internationally on the anti-jamming capabilities of LEO satellite SOP positioning technology. Due to significant differences in the downlink signal characteristics between LEO satellites and Global Navigation Satellite System (GNSS) signals based on Medium Earth Orbit (MEO) or Geostationary Earth Orbit (GEO) satellites, the anti-jamming research results of traditional satellite navigation systems cannot be directly applied. This study addresses the narrow bandwidth and high signal-to-noise ratio (SNR) characteristics of signals from LEO satellite constellations. We propose a Consecutive Iteration based on Signal Cancellation (SCCI) algorithm, which significantly reduces errors during the model fitting process. Additionally, an adaptive variable convergence factor was designed to simultaneously balance convergence speed and steady-state error during the iteration process. Compared to traditional algorithms, simulation and experimental results demonstrated that the proposed algorithm enhances the effectiveness of jamming threshold settings under narrow bandwidth and high-power conditions. In the context of LEO satellite jamming scenarios, it improves the frequency-domain anti-jamming performance significantly and holds high application value for drone positioning. Full article
(This article belongs to the Special Issue Advances of Drones in Green Internet-of-Things)
Show Figures

Figure 1

Figure 1
<p>Complex electromagnetic environment faced by UAVs.</p>
Full article ">Figure 2
<p>Iterative process flow chart.</p>
Full article ">Figure 3
<p>The relationship between the mean square error (MSE) and the number of iterations. (<b>a</b>) The impact of parameter α; (<b>b</b>) the impact of parameter λ.</p>
Full article ">Figure 4
<p>Convergence curves of each algorithm.</p>
Full article ">Figure 5
<p>Flow chart of SCCI algorithm.</p>
Full article ">Figure 6
<p>Change in power spectrum before and after anti-jamming under strong jamming for the SCCI algorithm and SCCME algorithm (JSR 30 dB). (<b>a</b>) Before anti-jamming; (<b>b</b>) after anti-jamming (SCCI); (<b>c</b>) after anti-jamming (SCCME).</p>
Full article ">Figure 7
<p>Change in power spectrum before and after anti-jamming under weak jamming for SCCI algorithm and SCCME algorithm (JSR 15 dB). (<b>a</b>) Before anti-jamming; (<b>b</b>) after anti-jamming (SCCI); (<b>c</b>) after anti-jamming (SCCME).</p>
Full article ">Figure 8
<p>The interference detection ratio (IDR). (<b>a</b>) Scene One; (<b>b</b>) Scene Two; (<b>c</b>) Scene Three.</p>
Full article ">Figure 9
<p>Verification of anti-jamming performance. (<b>a</b>) Scene One; (<b>b</b>) Scene Two; (<b>c</b>) Scene Three.</p>
Full article ">Figure 10
<p>Hardware connection diagram.</p>
Full article ">Figure 11
<p>Constellation map during the satellite visibility period.</p>
Full article ">Figure 12
<p>The anti-jamming results of the strong jamming scenario for the SCCI algorithm and SCCME algorithm (JSR 35 dB). (<b>a</b>) Before anti-jamming; (<b>b</b>) after anti-jamming (SCCI); (<b>c</b>) after anti-jamming (SCCME).</p>
Full article ">Figure 13
<p>The anti-jamming results of the weak jamming scenario for the SCCI algorithm and SCCME algorithm (JSR 15 dB). (<b>a</b>) Before anti-jamming; (<b>b</b>) after anti-jamming (SCCI); (<b>c</b>) after anti-jamming (SCCME).</p>
Full article ">
Back to TopTop