Next Article in Journal
Temporal Evolution of Refractive Index Induced by Short Laser Pulses Accounting for Both Photoacoustic and Photothermal Effects
Previous Article in Journal
Discovery and Characterization of Chemical Compounds That Inhibit the Function of Bacterial Neuraminidase from Codonopsis ussuriensis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

GNSS Signal Acquisition Algorithm Based on Two-Stage Compression of Code-Frequency Domain

1
Innovation Academy for Microsatellites, Chinese Academy of Sciences, Shanghai 201210, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
College of Electrical and Electronic Engineering, Wenzhou University, Wenzhou 325035, China
*
Author to whom correspondence should be addressed.
Submission received: 5 May 2022 / Revised: 3 June 2022 / Accepted: 9 June 2022 / Published: 20 June 2022

Abstract

:
The recently-emerging compressed sensing (CS) theory makes GNSS signal processing at a sub-Nyquist rate possible if it has a sparse representation in certain domain. The previously proposed code-domain compression acquisition algorithms have high computational complexity and low acquisition accuracy under high dynamic conditions. In this paper, a GNSS signal acquisition algorithm based on two-stage compression of the code-frequency domain is proposed. The algorithm maps the incoming signal of the same interval to multiple carrier frequency bins and overlaps the mapped signal that belongs to the same code phase. Meanwhile, the code domain compression is introduced to the preprocessed signal, replacing circular correlation with compressed reconstruction to obtain Doppler frequency and code phase. Theoretical analyses and simulation results show that the proposed algorithm can improve the frequency search accuracy and reduce the computational complexity by about 50% in high dynamics.

1. Introduction

The global navigation satellite system (GNSS) includes constellations of satellites orbiting over the earth’s surface that provide continuous positioning, navigation, and timing (PNT) services anywhere in the world or near-earth space, under any weather conditions [1,2,3]. The first stage of a GNSS receiver is the acquisition, which plays a very important role in signal synchronization, whose aim is to detect the signal and roughly estimate the code phase and the Doppler frequency of the satellite in the field of view.
GNSS signal acquisition is a two-dimensional searching process, which needs to search all the likely spreading code phases and Doppler frequencies, including serial acquisition and parallel acquisition. The serial acquisition is a common method, which needs to reconstruct all possible code phases and Doppler frequencies and correlate with the received satellite signal one by one, and is a very time-consuming search procedure. Thus, the parallel acquisition methods are proposed, which use fast Fourier transforms (FFT) to convert the signal from time domain into frequency domain. These methods can be summarized into three categories: the parallel frequency search, which uses a FFT to search simultaneously all or a part of the Doppler frequency bins [4,5]; the parallel code search, which uses FFTs to compute the code correlation to search all the code delay bins simultaneously [6,7]; and the two-dimensional search where FFTs are used for both frequency and code domains [8,9,10,11,12]. The parallel acquisition methods can improve the acquisition efficiency, but they increase the complexity considerably. Furthermore, the above acquisition methods are based on the Nyquist sampling theory, which is a redundant sampling method.
In recent years, as a promising theory of signal sampling, compressed sensing (CS) provides a new method to solve this dilemma [13], which can greatly reduce the sampling rates of sparse signals and faithfully reconstruct the original signals back from fewer compressive measurements. It has been applied in many fields, including speech enhancement [14], magnetic resonance imaging [15], wireless sensor network [16], etc. Moreover, several CS-based methods have been proposed for GNSS signal acquisition. He et al. [17] used the sparse matrix to compress the spread spectrum signal and reconstructed the compressed signal based on the greedy algorithm. Kong et al. [18] proposed two-stage compression of GNSS signal by the Walsh–Hadamard matrix, which achieved fast acquisition with fewer correlators than traditional acquisition methods. Chang et al. [19] combined CS with subspace to enhance the acquisition performance of GNSS signal when there are a lot of interferences. He et al. [20] designed a sparse dictionary based on the features of the GLONASS navigation signal and modified the greedy reconstruction algorithm to achieve fast acquisition. However, the above algorithms are based on the conventional frequency search method with a fixed bandwidth, which results in low acquisition probabilities under high dynamic conditions.
In this paper, a GNSS signal acquisition algorithm based on two-stage compression of the code-frequency domain is proposed, which divides the Doppler frequency offset into multiple intervals, and the multiple frequency points are compressed and searched in every interval at the same time to increase the correlation bandwidth. This method can overcome the limitations of the conventional Doppler frequency search and is suitable for fast acquisition in high dynamics.
The rest of this paper is organized as follows. Section 2 introduces the compression acquisition theory and the problems of the compression acquisition algorithm in a highly dynamic environment. Section 3 proposes the model of two-stage compression and theoretically analyzes the advantages of the proposed method, including computational complexity and detection probability. Section 4 compares the performance of the proposed algorithm with other acquisition algorithms and discusses the performance with different key parameters. Section 5 concludes this paper.

2. Code-Domain Compression Acquisition Algorithm

2.1. Sparsification of GNSS Signal

The incoming digital down-converted GNSS intermediate frequency (IF) signal can be expressed as [21]
r ( n ) = A d ( n ) c ( n τ ) exp [ j ( 2 π ( f IF + f d ) n T s + φ ) ] + v ( n )
where A is the signal amplitude, d is the signal data, c is the spreading code, f IF is the intermediate frequency, f d is the Doppler frequency offset, T s is the sampling interval, τ is the code phase, v is the noise, and  φ is the unknown carrier phase.
From the previous work, we know GNSS signal acquisition is possible at a low sampling rate by employing CS theory to compress the raw GNSS signal data using a measurement matrix [22], which can be expressed as
y = Φ · x
where x is the original signal, y is the measurement vector, Φ R M × N is the measurement matrix, N is the signal length, M is the measurement length, and  M < < N .
The premise of CS is that the signal is sparse. However, most of the real signals are not sparse but can be expressed as
x = Ψ · a
where Ψ is called the sparse basis of the signal, a is a sparse vector.
Shifting the spreading code vector c = [ c ( 0 ) , c ( 1 ) , , c ( N K 1 ) ] R N K × 1 to get the code matrix C
C = [ c 0 , c 1 , , c N 1 ] = c ( 0 ) c ( 1 ) c ( 2 K 1 ) c ( K ) c ( K + 1 ) c ( 3 K 1 ) c [ ( N 1 ) K ] c [ ( N 1 ) K + 1 ] c [ ( N + 1 ) K 1 ]
where N is the amount of points in the FFT, and N = 2 K .
Padding zeros to the end of the incoming signal r ( n ) to increase its length to 2 K , which can solve the problem of non-periodicity after block processing of spreading code and reduce computation effectively [23]. It can be described as
r = [ r ( 0 ) , r ( 1 ) , , r ( K 1 ) , 0 1 × K ] T
Then, correlating the padded signal vector r with the local spreading code matrix C , we can obtain the correlation vector X
X = C · r
The acorrelation function (CF) of the incoming signal and the spreading code can be described as
X i = < c n , r > = i = 0 2 K 1 c ( i + n ) r ( i ) = v + 2 A K 1 | τ | T c , | τ | T c v , | τ | > T c
where X i is the ith element of the correlation vector X , c n is the nth row of the code matrix C , < · > is the discrete-time convolution operation, T c is the chip time, 2 A K is the maximum peak value.
As shown in Figure 1, the correlation result X can be considered sparse if 2 A K > > v in Equation (7).

2.2. Code-Domain Compression Acquisition

According to Equations (2) and (6), we can use the measurement matrix Φ to compresses the N × 1 dimensional correlation result X to obtain the M × 1 dimensional measurement data Y
Y M × 1 = Φ M × N · X N × 1 = Φ M × N · C N × 2 K · r 2 K × 1
The measurement matrix Φ will be helpful for signal reconstruction when restricted isometry property (RIP) criteria is satisfied. In other words, each column is a unit vector, and the sum of squares of all columns meets the Welch bound [24,25]
i = 1 N j = 1 N ϕ i , ϕ j 2 N 2 M
The FFT-based parallel search engine is adopted to acquire the correlation peak, which converts time-domain convolution operation into frequency-domain multiplication. The FFT of Φ · C are taken and multiplied by the FFT of the r , and then taken the inverse FFT (IFFT) to obtain Y .
The correlation result X can be reconstructed from Y through some reconstruction algorithms, which can be described as follows
min X 0 s . t . Y = Φ · X
where · 0 means the l 0 norm which counts the number of nonzero elements.
The l 0 norm is a NP-hard problem, which can be converted into l 2 norm if the measurement matrix Φ meets the RIP condition.
min X 2 2 s . t . Y = Φ · X
where · 2 is the l 2 norm, also known as the Euclidean norm, which is used as a standard quantity for measuring a vector difference [26].
The least square method (LSM) is used to reconstruct the elements from the measurement data Y [27].
X = Φ T ( Φ Φ T ) 1 Y
The correct code phase and Doppler frequency offset are determined by comparing the first K elements from the correlation result X . Replacing the circular correlation with a compressed reconstruction in this method decreases the number of correlators from N + 1 to M + 1 [28].
As we know, the correlation function R ( τ , f d ) of the spreading code can be described as follows [29]
R ( τ , f d ) = 2 A K 1 τ T c sin [ π f d f i 2 K T s ] sin [ π f d f i T s ] exp j γ ( φ )
where γ ( φ ) = π f d f i 2 K 1 T s + φ , f i is the Doppler frequency hypothesis.
We can find that R ( τ , f d ) is affected by the Doppler frequency as a sinc ( · ) function [29], and the zero of its main lobe is located at
f d = ± 1 K T s + f i
The correlation bandwidth B c denotes
B c = 2 K T s
Hence the Doppler frequency offset range in each search is obtained as
f d 1 K T s + f i , 1 K T s + f i
where the search accuracy is 1 / K T s .
In above algorithm, the acquisition is serial in the frequency domain, and it will take a long time to acquire the Doppler frequency in high dynamics.

3. Code-Frequency Two-Stage Compression Acquisition

In this section, we introduce a new frequency domain compression to preprocess the incoming signal to reduce the influence of large Doppler frequency offset.

3.1. Frequency Domain Compression Preprocessing

The incoming signal is multiplied with exp [ j 2 π ( f o i + f c l ) ] to produce the mapped signal r l ( n )
r l ( n ) = r ( n ) exp [ j 2 π ( f o i + f c l ) n T s ]
where f c l is the multiple carrier frequency bins, l { 0 , 1 , 2 , , L 1 } , L is the frequency compression ratio (FCR), the frequency search step f step = f c l + 1 f c l , f o i is the starting frequency in once frequency compression search, i { 0 , 1 , 2 , , I 1 } , I is the number of sub-carrier bins, the frequency compression steps f o i + 1 f o i = L · f step .
Hence the Doppler frequency offset can be expressed as
f d = I · L · f step
Figure 2 shows the process of the mapping and overlapping-based carrier frequency searching technique. We divide the Doppler frequency offset f d into I sub-carrier frequency bins and search for L sub-frequency points separated by f step simultaneously.
From Equation (13), we can see that if the Doppler frequency is greater than the correlation bandwidth B c = 2 / K T s , the correlation value is attenuated seriously. Accordingly, the frequency search step should satisfy f step B c / 2 , making sure at least two neighboring acquisition frequency points are located within the correlation bandwidth B c , which is beneficial to enhancing the acquisition performance, as shown in Figure 3a. However, if  f step > B c , the neighbor acquisition frequency exceeds the correlation bandwidth B c , causing the correlation value decrease and the acquisition failure, as illustrated in Figure 3b.
Then, the mapped signals which correspond to the same code phase are overlapped together. Consequently, we obtain the preprocessed incoming signal
r l ( n ) = l = 0 L 1 r l ( n )
The correlation bandwidth of the comprssion acquisition is improved by frequency domain preprocessing, which can be expressed as
B c = ( L 1 ) f step + B c
Compared with the traditional Doppler frequency search method, the L frequency bins are searched simultaneously in this method, which can reduce search number and acquisition time.

3.2. Algorithm Process

The detailed and straightforward processes of proposed GNSS acquisition algorithm are shown in Algorithm 1 and Figure 4.
Algorithm 1 The proposed GNSS acquisition algorithm.
Input: Incoming signal, r ( n ) ; Spreading code vector, c 0 ; Measurement matrix, Φ ; Decision
threshold, V t ;
Output: Code phase, τ ; Doppler frequency offset, f d ;
1:Mapping incoming signal r ( n ) to different carrier frequencies f o i ;
2:Overlapping the mapped signals to the same code phase to obtain r l ( n ) ;
3:Padding zeros to the end of the mapped signals r l ( n ) to form a signal vector r ;
4:Shifting spreading code vector c 0 to obtain the code matrix C ;
5:Performing FFT operation to find out the measurement data Y ;
6:Using the greedy algorithm to obtain the correlation result X ;
7:Comparing the first K elements from the correlation result X with the threshold V t to
obtain the code phase τ and Doppler shift f d ;
8:return τ , f d ;

3.3. Detection Probability Analysis

Assuming a partial Hadamard matrix is adopted as measurement matrix to assess the performance of the proposed algorithm, the following two hypotheses determine the test statistic for signal detection
y ˜ = Φ 1 w H 0 Φ 1 ( x + w ) H 1
where Φ 1 R M T p / α c × N denotes a random partial Hadamard matrix, α c is the maximum subsampling factor, T p is the spreading spectrum period, H 1 denotes the signal is present and correctly aligned with the local replica, H 0 denotes the signal is absent or not correctly aligned with local replica.
Consequently, the probability density function of y ˜ under hypothesis H 1 is shown as follows:
p y ˜ x ¯ ^ , H 1 = 1 ( 2 π ) M / α c / 2 σ 2 Φ 1 Φ 1 T 1 / 2 exp 1 2 ( W ) T σ 2 Φ 1 Φ 1 T 1 ( W )
where W = y ˜ D x ^ , D = Φ 1 · Ψ , and x ^ denotes
x ^ = ( Φ 1 Φ 1 T ) 1 Φ 1 T y ˜
The x ¯ ^ denotes the maximum likelihood estimation (MLE) of x under H 1 , which is shown as
x ¯ ^ = D T ( σ 2 Φ 1 Φ 1 T ) 1 D 1 D T ( σ 2 Φ 1 Φ 1 T ) 1 y ˜
If the threshold η is selected to attain the required false alarm probability, the detection performance can be obtained as
P d = Pr Γ y ˜ > η H 1 = Q χ ¯ M T p / α c 2 λ ( λ ) ( η )
where Γ y ˜ is a noncentral chi-squared distribution with M T p / α c degrees of freedom, which is expressed as
Γ y ˜ = y ˜ T ( σ 2 Φ 1 Φ 1 T ) 1 D D T ( σ 2 Φ 1 Φ 1 T ) 1 D 1 D T ( σ 2 Φ 1 Φ 1 T ) 1 y ˜
where λ is a non-centrality parameter specified by
λ = σ 2 ( x T ) 1 Φ 1 T Φ 1 Φ 1 T 1 Φ 1 x

3.4. Complexity Analysis

The complexity of the proposed algorithm consists of three parts: frequency domain compression preprocessing, code domain compression, and signal reconstruction, of which the first and third are the main parts, and the second part can be ignored as the elements of the measurement matrix are −1 or 1. In this paper, we chose the orthogonal matching pursuit (OMP) algorithm for signal reconstruction, whose complexity is O ( N · ( log 2 ( N ) ) 2 ) [30,31]. Table 1 gives the complexities of the serial acquisition algorithm, the PMF-FFT acquisition algorithm, the code-compression (CC) acquisition algorithm [26] and the proposed code-frequency compression (CFC) acquisition algorithm.

4. Simulation Results

The GPS L1 signal is adopted to evaluate the performance of the proposed method. Set the data rate R b = 50 bps, the code rate R c = 1.023 Mcps, the sampling rate f s = 2 R c , the code phase τ = 400, the code length N = 1024, the measured length M = 768, and the Doppler frequency offset range is ±80 KHz.

4.1. Results on Detection Probability

Figure 5 shows the normalized acquisition results of the proposed algorithm with different F C R when S N R is −12 dB. We can find that with the decrease of F C R , the ratio of the highest correlation peak to the second correlation peak increases, which means the detection probability increases.
Figure 6 shows the detection probabilities of the proposed algorithm with different S N R . We can compute that when the S N R = −12 dB, the detection probabilities of F C R = 2:1, 4:1, and 8:1 are 96.3%, 67.8%, and 22.4%, respectively. With the increase of F C R , the number of overlaps of the mapped signals with the same code phase increases, which makes it difficult to distinguish the correlation peaks and reduces the detection probability.

4.2. Results on Dynamic Acquisition

Figure 7 shows the detection probabilities of the proposed algorithm with different Doppler frequencies. When the S N R = −12 dB, F C R = 2:1, the Doppler frequency offsets are 0kHz, 26 kHz, 52 kHz and 78 kHz, the detection probabilities are 96.3%, 93.2%, 98.5%, and 92.3%, respectively, and their values are very close to each other. When the F C R is changed to 4:1 or 8:1, we can get a similar conclusion. As analyzed in Section 3.1, the frequency-domain preprocessing improves the correlation bandwidth, keeping more neighboring acquisition frequency points within the correlation bandwidth, which greatly reduces its sensitivity to Doppler frequency.
Figure 8 compares the detection probabilities of the proposed algorithm with other traditional acquisition algorithms. With the increase of Doppler frequency, the acquisition probabilities of the PMF-FFT and the CC acquisition algorithms decrease. A large Doppler frequency restricts the matching filter length of the PMF-FFT acquisition algorithm [32] and the correlation bandwidth of the CC acquisition algorithm, which limits their application in highly dynamic situations. As we can see, when the Doppler frequency is greater than about 52 kHz, the acquisition performance of the CFC acquisition algorithm begins to surpass that of the other two algorithms.

4.3. Results on Algorithm Complexity

Figure 9 gives the complexities of the four acquisition algorithms with different F C R and N. As we can see, the complexity of the proposed CFC acquisition algorithm is less than that of the other three algorithms. This advantage becomes more obvious with the increase of parameters F C R and N. Specifically, as shown in Figure 9a, when the F C R = 2, the code length N = 1024, the complexities of the serial, the PMF-FFT, the CC, and the CFC acquisition algorithms are 1.049 × 10 7 , 1.114 × 10 6 , 1.638 × 10 6 , 8.192 × 10 5 , respectively, and the proposed CFC acquisition algorithm is 92.19%, 26.46%, and 50% lower than the other three algorithms. As shown in Figure 9b, when the code length N = 2048, the F C R = 1:2, the complexities of the four acquisition algorithms are 6.711 × 10 7 , 4.325 × 10 6 , 3.965 × 10 6 , 1.982 × 10 6 , respectively. Compared with the other algorithms, the CFC acquisition algorithm reduces 99.7%, 54.17%, and 50%. Thus, the frequency domain compression preprocessing can reduce the resource overhead, which makes the proposed CFC acquisition algorithm better adapt to high dynamic conditions.

5. Conclusions

In this paper, we proposed a GNSS signal acquisition algorithm based on two-stage compression of the code-frequency domain. The preprocessing of frequency compression is employed to reduce the computational complexities of the CC acquisition algorithm, which maps the incoming signal into different sub-carrier frequency bins, and the mapped signals are then overlapped with the same code phase, which can reduce the number of bins that should be searched. The simulation results show that when the Doppler frequency offset is greater than 52 kHz, the detection probability of the proposed algorithm is better than that of the only code-domain compression acquisition algorithms, and the computational complexity is reduced by 50% at the same time. In addition, compared with the traditional PMF-FFT and CC acquisition algorithms, the detection probability of the proposed CFC algorithm is insensitive to the Doppler frequency offsets in highly dynamic environments. In the future, the measurement matrix Φ could be optimized and the reconstruction method could be explored to further improve the acquisition performance. Moreover, the impacts of multipath and non-Gaussian noise on compression acquisition are also very interesting research topics.

Author Contributions

Conceptualization, F.Z. proposed the main idea, finished the draft manuscript; F.Z., L.L. conceived of the experiments and drew the figures and tables; methodology, F.Z.; L.Z. and L.L. analyzed the data; F.Z. and L.L.wrote the paper; X.J.; Y.H.; J.Y. and G.L. reviewed the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Chinese Academy of Sciences Youth Innovation Promotion Association, grant numbers Y202068, Shanghai Industrial Collaborative Innovation Project, grant numbers No.2021-CYXT2-KJ03, Wenzhou Major Scientific and Technological Innovation Projects of China, grant numbers ZG2021029, Scientific Research Project of Zhejiang Provincial Department of Education, grant numbers Y202146796, Design and Key Technology of Electro-Optic Switching for Satellite Internet Project, grant numbers KJ002-02, National Natural Science Foundation of China, grant numbers U21A20443.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
GNSSGlobal Navigation Satellite System
CSCompressed sensing
FFTFast Fourier Transform
IFFTInverse Fast Fourier Transform
RIPRestricted Isometry Property
CFCCode-frequency compression
CCCode-compression
FCRFrequency Compression Ratio
SNRSignal-to-Noise Ratio

References

  1. Lu, K.; Wang, X.; Shen, L.; Chen, D. A GPS signal acquisition algorithm for the high orbit space. GPS Solut. 2021, 25, 92. [Google Scholar] [CrossRef]
  2. Xiao, Y.; Zhou, X.; Wang, J.; He, Z.; Zhou, H. Observability Analysis and Navigation Filter Optimization of High-Orbit Satellite Navigation System Based on GNSS. Appl. Sci. 2020, 10, 7513. [Google Scholar] [CrossRef]
  3. Carrasco, P.; Cuesta, F.; Caballero, R.; Perez-Grau, F.J.; Viguria, A. Multi-Sensor Fusion for Aerial Robots in Industrial GNSS-Denied Environments. Appl. Sci. 2021, 11, 3921. [Google Scholar] [CrossRef]
  4. Le, W.; Wenzhao, T.; Mi, H. High Dynamic Spread Spectrum Signal Acquisition Algorithm based on Delay Autocorrelation. In Proceedings of the 2018 International Conference on Sensor Networks and Signal Processing (SNSP), Xi’an, China, 28–31 October 2018; pp. 382–385. [Google Scholar] [CrossRef]
  5. Leclère, J.; Landry, R.; Botteron, C. Comparison of L1 and L5 Bands GNSS Signals Acquisition. Sensors 2018, 18, 2779. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Borre, K.; Akos, D.M.; Bertelsen, N.; Rinder, P.; Jensen, S.H. A Software-Defined GPS and Galileo Receiver: A Single-Frequency Approach; Springer Science & Business Media: Berlin, Germany, 2007. [Google Scholar]
  7. Kim, B.; Kong, S.H. Design of FFT-Based TDCC for GNSS Acquisition. IEEE Trans. Wirel. Commun. 2014, 13, 2798–2808. [Google Scholar] [CrossRef]
  8. Foucras, M. Performance Analysis of the Modernized GNSS Signal Acquisition. Ph.D. Thesis, Institut National Polytechnique de Toulouse, Toulouse, France, 2015. [Google Scholar]
  9. Xu, Y.; Xu, L.; Yuan, H.; Luo, R. Direct P-code acquisition algorithm based on bidirectional overlap technique. J. Syst. Eng. Electron. 2014, 25, 538–546. [Google Scholar]
  10. Zhang, Y.; Li, Q. Fast Acquisition Algorithm for GPS L5 Signal Based on Folding. In Proceedings of the 2017 International Conference on Computer Technology, Electronics and Communication (ICCTEC), Dalian, China, 19–21 December 2017; pp. 439–442. [Google Scholar]
  11. Zeng, Q.; Tang, L.; Zhang, P.; Pei, L. Fast acquisition of L2C CL codes based on combination of hyper codes and averaging correlation. J. Syst. Eng. Electron. 2016, 27, 308–318. [Google Scholar] [CrossRef]
  12. Li, H.; Lu, M.; Feng, Z. Mapping and overlapping based carrier frequency searching technique for rapid GNSS long PN-code acquisition. Sci. China Inf. Sci. 2010, 53, 2642–2652. [Google Scholar] [CrossRef]
  13. Holger, B.; Robert, C.; Gitta, K. Compressed Sensing and Its Applications; Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  14. Haneche, H.; Boudraa, B.; Ouahabi, A. A new way to enhance speech signal based on compressed sensing. Measurement 2020, 151, 107–117. [Google Scholar] [CrossRef]
  15. Sandino, C.M.; Cheng, J.Y.; Chen, F.; Mardani, M.; Pauly, J.M.; Vasanawala, S.S. Compressed Sensing: From Research to Clinical Practice with Deep Neural Networks: Shortening Scan Times for Magnetic Resonance Imaging. IEEE Signal Processing 2020, 37, 117–127. [Google Scholar] [CrossRef]
  16. Liang, J.; Li, L.; Zhao, C. A Transfer Learning Approach for Compressed Sensing in 6G-IoT. IEEE Internet Things J. 2021, 8, 15276–15283. [Google Scholar] [CrossRef]
  17. He, G.; Song, M.; He, X.; Hu, Y. GPS signal acquisition based on compressive sensing and modified greedy acquisition algorithm. IEEE Access. 2019, 7, 40445–40453. [Google Scholar] [CrossRef]
  18. Kong, S.H. A Deterministic Compressed GNSS Acquisition Technique. IEEE Trans. Veh. Technol. 2012, 62, 511–521. [Google Scholar] [CrossRef]
  19. Chang, C.L. Modified compressive sensing approach for GNSS signal reception in the presence of interference. GPS Solut. 2016, 20, 201–213. [Google Scholar] [CrossRef]
  20. He, G.; Song, M.; Zhang, S.; Song, P.; Shu, X. Sparse GLONASS Signal Acquisition Based on Compressive Sensing and Multiple Measurement Vectors. Math. Probl. Eng. 2020, 2020, 9654120. [Google Scholar] [CrossRef]
  21. Chao, W.; Erxiao, L.; Zhihua, J. Two-step compressed acquisition method for Doppler frequency and Doppler rate estimation in high-dynamic and weak signal environments. J. Syst. Eng. Electron. 2021, 32, 831–840. [Google Scholar] [CrossRef]
  22. Bi, X.; Leng, L.; Kim, C.; Liu, X.; Du, Y.; Liu, F. Constrained Backtracking Matching Pursuit Algorithm for Image Reconstruction in Compressed Sensing. Appl. Sci. 2021, 11, 1435. [Google Scholar] [CrossRef]
  23. Han, L.; Meng, Y.; Wang, Y.; Han, X. A Fast Algorithm of GNSS-R Signal Processing Based on DBZP. In Proceedings of the China Satellite Navigation Conference (CSNC) 2017 Proceedings: Volume I; Sun, J., Liu, J., Yang, Y., Fan, S., Yu, W., Eds.; Springer: Singapore, 2017; pp. 187–197. [Google Scholar]
  24. Datta, S. Welch bounds for cross correlation of subspaces and generalizations. Linear Multilinear Algebra. 2016, 64, 1484–1497. [Google Scholar] [CrossRef]
  25. Waldron, S. A Sharpening of the Welch Bounds and the Existence of Real and Complex Spherical t–Designs. IEEE Trans. Inf. Theory 2017, 63, 6849–6857. [Google Scholar] [CrossRef]
  26. Dumitrescu, B.; Irofti, P. Dictionary Learning Algorithms and Applications; Springer: Berlin/Heidelberg, Germany, 2018. [Google Scholar]
  27. Zhang, J.; Liu, J.; Fan, W.; Qiu, W.; Luo, J. Partial Hadamard Encoded Synthetic Transmit Aperture for High Frame Rate Imaging with Minimal l2 norm Least Squares Method. Phys. Med. Biol. 2022, 67, 105002. [Google Scholar] [CrossRef]
  28. Wang, J.; Wang, W.; Chen, J. Adaptive Rate Block Compressive Sensing Based on Statistical Characteristics Estimation. IEEE Trans. Image Processing 2021, 31, 734–747. [Google Scholar] [CrossRef] [PubMed]
  29. Kong, S.H.; Kim, B. Two-Dimensional Compressed Correlator for Fast PN Code Acquisition. IEEE Trans. Wirel. Commun. 2013, 12, 5859–5867. [Google Scholar] [CrossRef]
  30. Zhu, H.; Chen, W.; Wu, Y. Efficient Implementations for Orthogonal Matching Pursuit. Electronics 2020, 9, 1507. [Google Scholar] [CrossRef]
  31. Lee, J.; Gil, G.T.; Lee, Y.H. Channel Estimation via Orthogonal Matching Pursuit for Hybrid MIMO Systems in Millimeter Wave Communications. IEEE Trans. Commun. 2016, 64, 2370–2386. [Google Scholar] [CrossRef]
  32. Deng, Z.; Jia, B.; Tang, S.; Fu, X.; Mo, J. Fine Frequency Acquisition Scheme in Weak Signal Environment for a Communication and Navigation Fusion System. Electronics 2019, 8, 829. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Model of correlation output.
Figure 1. Model of correlation output.
Applsci 12 06255 g001
Figure 2. Schematic diagram of frequency domain compression processing.
Figure 2. Schematic diagram of frequency domain compression processing.
Applsci 12 06255 g002
Figure 3. The relationship between frequency search step and correlation bandwidth: (a) f step B c / 2 ; (b) f step > B c .
Figure 3. The relationship between frequency search step and correlation bandwidth: (a) f step B c / 2 ; (b) f step > B c .
Applsci 12 06255 g003
Figure 4. GNSS signal acquisition based on code-frequency compression.
Figure 4. GNSS signal acquisition based on code-frequency compression.
Applsci 12 06255 g004
Figure 5. The acquisition results with different F C R : (a) F C R = 8:1; (b) F C R = 4:1; (c) F C R = 2:1.
Figure 5. The acquisition results with different F C R : (a) F C R = 8:1; (b) F C R = 4:1; (c) F C R = 2:1.
Applsci 12 06255 g005
Figure 6. The detection probabilities with different S N R .
Figure 6. The detection probabilities with different S N R .
Applsci 12 06255 g006
Figure 7. The detection probabilities with different Doppler frequencies.
Figure 7. The detection probabilities with different Doppler frequencies.
Applsci 12 06255 g007
Figure 8. The detection probabilities of several algorithms with different S N R : (a) Doppler = 0 kHz; (b) Doppler = 26 kHz; (c) Doppler = 52 kHz; (d) Doppler = 78 kHz.
Figure 8. The detection probabilities of several algorithms with different S N R : (a) Doppler = 0 kHz; (b) Doppler = 26 kHz; (c) Doppler = 52 kHz; (d) Doppler = 78 kHz.
Applsci 12 06255 g008
Figure 9. Algorithm complexities with different F C R and N: (a) F C R varies from 2 to 10; (b) N varies from 256 to 4096.
Figure 9. Algorithm complexities with different F C R and N: (a) F C R varies from 2 to 10; (b) N varies from 256 to 4096.
Applsci 12 06255 g009
Table 1. Table of algorithm complexities.
Table 1. Table of algorithm complexities.
Acquisition AlgorithmComplexity
Serial acquisition algorithm O ( N 2 · I · L )
PMF-FFT acquisition algorithm O ( N · ( L · log 2 ( L ) + N ) )
CC acquisition algorithm O ( N · ( log 2 ( N ) ) 2 · I · L )
CFC acquisition algorithm O ( N · ( log 2 ( N ) ) 2 · I )
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, F.; Zhao, L.; Li, L.; Hu, Y.; Jiang, X.; Yu, J.; Liang, G. GNSS Signal Acquisition Algorithm Based on Two-Stage Compression of Code-Frequency Domain. Appl. Sci. 2022, 12, 6255. https://0-doi-org.brum.beds.ac.uk/10.3390/app12126255

AMA Style

Zhou F, Zhao L, Li L, Hu Y, Jiang X, Yu J, Liang G. GNSS Signal Acquisition Algorithm Based on Two-Stage Compression of Code-Frequency Domain. Applied Sciences. 2022; 12(12):6255. https://0-doi-org.brum.beds.ac.uk/10.3390/app12126255

Chicago/Turabian Style

Zhou, Fangming, Lulu Zhao, Limin Li, Yifei Hu, Xinglong Jiang, Jinpei Yu, and Guang Liang. 2022. "GNSS Signal Acquisition Algorithm Based on Two-Stage Compression of Code-Frequency Domain" Applied Sciences 12, no. 12: 6255. https://0-doi-org.brum.beds.ac.uk/10.3390/app12126255

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