Next Article in Journal
A Stacked Back Side-Illuminated Voltage Domain Global Shutter CMOS Image Sensor with a 4.0 μm Multiple Gain Readout Pixel
Next Article in Special Issue
A New DGNSS Positioning Infrastructure for Android Smartphones
Previous Article in Journal
Optimal Energy Resources Allocation Method of Wireless Sensor Networks for Intelligent Railway Systems
Previous Article in Special Issue
A New Cycle Slip Detection and Repair Method Using a Single Receiver’s Single Station B1 and L1 Frequencies in Ground-Based Positioning Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unambiguous Acquisition/Tracking Technique Based on Sub-Correlation Functions for GNSS Sine-BOC Signals

1
College of Automation, Harbin Engineering University, No.145, Nantong Street, Harbin 150001, China
2
State Key Laboratory of Satellite Navigation System and Equipment Technology, Shijiazhuang 050081, China
3
The 54th Research Institute of China Electronics Technology Group Corporation, Shijiazhuang 050081, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(2), 485; https://doi.org/10.3390/s20020485
Submission received: 8 November 2019 / Revised: 9 January 2020 / Accepted: 10 January 2020 / Published: 15 January 2020

Abstract

:
The autocorrelation function (ACF) of the Binary Offset Carrier modulation (BOC) signal for Global Navigation Satellite System (GNSS) has multiple peaks, ambiguity is easily generated during the synchronization of the baseband signal. Some methods have been proposed to remove the ambiguity, but the performance is not suitable for high-order BOC signals or does not maintain narrow correlation characteristics. This paper proposes a sub-function reconstruction synchronization algorithm to solve this problem, of which the key is to design a new local auxiliary code: the local Pseudo-Random Noise (PRN) code is divided into several new codes with different delays. The auxiliary code performs a coherent integration operation with the received signal. Then, a correlation function without any positive side peaks is obtained by multiplying the two correlation results to make the acquisition/tracking completely unambiguous. The paper gives a design scheme of navigation signal acquisition/tracking and deduces the theoretical analysis of detection performance. The phase discrimination function is provided. The performance of the method is analyzed from both theoretical and simulation aspects. Compared with the Binary phase shift keying-like (BPSK-LIKE) method, Subcarrier Phase Cancellation (SCPC) method and the Autocorrelation Side-Peak Cancellation Technique (ASPeCT) method, the proposed method has the best detection probability for the acquisition, which is 0.5 dB-Hz better than ASPeCT. For tracking, the proposed method performs best in terms of phase-detection curve, anti-multipath performance, and anti-noise performance. For high-order BOC signals, the SRSA technique successfully removes the false lock points, and there is only one multipath error envelope, and the code tracking error is almost the same as the ASPeCT method.

1. Introduction

Compared to traditional the Global Navigation Satellite System (GNSS), the addition of BOC modulation signals marks a breakthrough in the development of modern GNSS. The frequency spectrum splitting characteristic of the BOC signal helps different GNSS signals efficiently utilize limited band resources. The ACF of the BOC modulation signal also has a narrow correlation main peak, and a large main peak slope allows the signal to have better confidentiality and location accuracy. However, there is a severe ambiguity in the synchronization processing of BOC signal, because the BOC signal has not only the spectral splitting characteristics, but is also characterized by correlation function multi-peak. Only the main peak is the correct sync position of the corresponding signal, while other side peaks will cause false acquisition and mis-locking of the BOC signal. Although the BOC modulation signal possesses stronger anti-multipath performance compared with the conventional BPSK modulation signal, the multipath error is still the primary source of localization error of the BOC signal [1]. Therefore, how to remove the ambiguity of the BOC signal and enhance its anti-multipath performance has become a research hotspot in the field of navigation.
At present, there are three main types of processing methods for acquisition and tracking the BOC signal: The first type is to avoid the ambiguity detection method. These methods, such as the double estimation technique (DET) method [2,3] and the Bump-jump (BJ) method [4], detect the occurrence of error acquisition or mis-locking during the synchronization process of the BOC signal by adding a correlator. Because such methods require more extended detection and recovery times, they are only suitable for strong signal situations. The second category is the idea of degraded processing. This method uses the spectral relationship between BOC and BPSK to convert the BOC correlation function into a BPSK- LIKE correlation function with only a single peak, such as the BPSK-LIKE method [5]. The correlation main peak width obtained by this method is the same as for BPSK, which sacrifices the advantages of anti-multipath and high precision brought by the narrow correlation characteristics of the BOC signal. The third type is the idea of combining a single narrow main peak. These methods achieve the elimination or suppression of the sub-peak through a particular combination algorithm by locally generating a unique auxiliary code, and retains or generates a new single narrow main peak, such as ASPeCT [6], SCPC [7], GRASS [8] and PUDLL [9]. The ASPeCT method combines an unambiguous new function by the ACF of the BOC signal and its cross-correlation function (CCF) with the PRN code. But it only applies to the BOC(n, n) signal. The SCPC method performs coherent integration by receiving a signal with a local BOC signal modulated by two orthogonal subcarriers to obtain an unambiguous new function. Although this method is applicable to BOC(m, n) signals, it does not retain the narrow correlation characteristics of BOC signals. PUDLL [9] and GRASS [8] design some local auxiliary signals and generated a new correlation function with no side-peaks. However, as the PUDELL method increases with the order of BOC modulation, the code tracking performance decreases rapidly, and the GRASS method sacrifices performance. These methods take advantage of the BOC signal and become a vital research program.
This paper proposes a sub-function reconstruction synchronization algorithm (SRSA), which can solve the ambiguity. The structure of the acquisition and tracking algorithm is designed. The corresponding unambiguous non-coherent discrimination function for tracking process and the principle of combination are presented. The acquisition/tracking performance of this algorithm is analyzed from both theoretical and simulation aspects. Comparing BPSK-LIKE, SCPC, and ASPeCT with the proposed algorithm, simulation experiments show that the proposed algorithm can remove the interference of side peaks while retaining the narrow main peak characteristics of the correlation function. And code tracking accuracy and acquisition sensitivity have been significantly improved.
The organization of this paper is as follows. Section 2 presents the ambiguity problem, defines a new local auxiliary signal, and puts forward an algorithm for solving the ambiguity problem. In Section 3, an unambiguous acquisition scheme is provided, and the acquisition performance of the proposed algorithm is theoretical and simulated. Section 4 demonstrates tracking methods and the phase discrimination function. The tracking performance of the proposed algorithm is analyzed from three aspects of phase discrimination, anti-multipath, and anti-noise. And the final section draws a conclusion.

2. The Combination Principle of SRSA

2.1. The Characteristics of BOC(Kn,n) Correlation Function and Problems

According to the BOC modulation principle, BOC(kn, n) can be expressed as the following:
S B O C ( k n , n ) ( t ) = A · D ( t ) · c ( t ) · s i g n ( sin ( 2 π f s c t ) )
In the above formula, A represents the signal amplitude; D ( t ) represents the navigation signal data code; c ( t ) is PRN code; Subcarriers are represented by s i g n ( sin ( 2 π f s c t ) ) . Using f 0 = 1.023   MHz as the reference frequency, BOC(kn, n) has the characteristics of PRN code rate f c = n f 0 and subcarrier frequency f s c = k n f 0 . N = 2 k n / n = 2 k is defined as the modulation factor of BOC(kn, n). The BOC(kn, n) code after PRN code and subcarrier modulation can be expressed as:
S B O C ( t ) = c ( t ) · s i g n ( sin ( 2 π f s c t ) )
Since the PRN code and subcarrier have specific correlation characteristics, the correlation function of the BOC(kn, n) signal has multiple peaks. The mathematical formula of the autocorrelation function R B ( τ ) and the cross-correlation function R B / P ( τ ) of the BOC(kn, n) is [10]
R B ( τ ) = i = 1 N 1 { [ ( 1 ) i t r i ( τ + T c N ) + ( 1 ) i t r i ( τ T c N ) ] × N i N } + t r i 0 ( τ )
R B / P ( τ ) = i = 1 N / 2 { [ t r i ( τ + ( 2 i 1 ) T c N ) + ( 1 ) t r i ( τ ( 2 i 1 ) T c N ) ] × 1 N }
where t r i ( τ a T c ) represents a triangular function with peak at a · T c . The triangle duration is 2 T s c , and T s c = T c / N . T c = 1 / f c is a PRN code chip duration. τ indicates the code phase delay. Figure 1 is a comparison of autocorrelation functions of BPSK, BOC(n, n) and BOC(2n, n). Figure 2 is a comparison of the cross-correlation functions between the BOC(n, n)/BOC(2n, n) signals and the PRN code. As can be seen from the figures, the larger the modulation factor, the greater the number of correlation peaks.
However, for BOC signals, once the side peaks are acquired, the false locks will occur during the subsequent tracking phase [11]. Especially for high-order BOC signals, it is difficult to distinguish the main peak and the sub-peak of the correlation function, which is more likely to cause the false lock phenomenon. The relationship between the false acquisition probability of different factors N and the signal-to-noise ratio (SNR) is given in Figure 3. SNR is defined as the ratio of signal to noise in an electronic device or electronic system. It can be seen that in the same SNR environment, the larger the modulation factor, the greater the probability of false acquisition. At −27 dB, the false acquisition probability of N = 4 and N = 2 are the same, close to zero, and the false acquisition probability of N = 8 is still 0.05.
In the code tracking stage, the classic EMLP (Early-Minus-Late-Power) phase discriminator is used, and the phase discrimination formula is
V ( Δ τ ) = R B 2 ( Δ τ d / 2 ) R B 2 ( Δ τ + d / 2 )
Δ τ represents the estimation error of the code phase. d is the delay interval between the early and late correlators. Figure 4 shows the discrimination curves for the BPSK signal, the BOC(1,1) signal, and the BOC(10, 5) signal in an infinite bandwidth. False-lock point corresponds to a zero of the discriminator where the S-curve presents a positive slope. Obviously, the larger the modulation factor, the more false-lock points, which are generated by the side peaks of ACF. This situation will result in intolerable code tracking errors.

2.2. Generation of Sub-PRN Code Sequences

By definition, the BOC code and the PRN code are expressed in the form of pulses as follows.
S P R N ( t ) = i = C i P T C ( t i T C ) = i = 0 M S P R N , i ( t )
s B O C ( t ) = j = 0 M 1 i = C i ( 1 ) j P T s c ( t i T C j T s c )
P T C is a rectangular pulse with a period of T C and an amplitude of 1. P T s c is a rectangular pulse with a duration of T s c and an amplitude of 1, and T s c = T c / N . C i is the symbol of the chip, C i ( 1 , 1 ) . M is the total number of pulses in a chip interval, which means that each PRN code chip is divided into M segments. Here M = N. This paper uses S P R N , i ( t ) to represent the i-th sub-PRN code. Then, separate the local PRN code, retain the i-th segment, and set the rest to zero. When M = 2, as shown in Figure 5, when M = 4, as shown in Figure 6.
It is supposed that the coherent integration time is T c o h , then the cross-correlation function of BOC code and PRN code is:
R B / P ( τ ) = 1 T c o h 0 T c o h S B O C ( t ) · S P R N ( t + τ ) d t = 1 T c o h T c o h T c τ T c S B O C ( t ) · S P R N ( t + τ ) d t = i = 1 M R B / P i ( τ )
R B / P i ( τ ) = 1 T c o h 0 T c o h S B O C ( t ) · S P R N , i ( t + τ ) d t
Figure 7 is a schematic diagram of sub-code separation for M = 2. And Figure 8 is a schematic diagram of sub-code separation for M = 4.
And it can also be represented by triangular functions. When M = 2, it can be expressed as
R B ( n , n ) / p 1 = 1 2 i = 1 1 t r i ( τ + ( 2 i 1 ) T s c T s c ) t r i ( τ ( 2 i 1 ) T s c )
R B ( n , n ) / p 2 = 1 2 i = 1 1 t r i ( τ + ( 2 i 1 ) T s c ) t r i ( τ ( 2 i 1 ) T s c + T s c )
And when M = 4 and the triangle duration is 2 T s c , it can be expressed as
R B ( 2 n , n ) / p 1 = 1 4 i = 1 2 t r i ( τ + ( 2 i 1 ) T s c 3 T s c ) t r i ( τ ( 2 i 1 ) T s c )
R B ( 2 n , n ) / p 2 = 1 4 i = 1 2 t r i ( τ + ( 2 i 1 ) T s c 2 T s c ) t r i ( τ ( 2 i 1 ) T s c + T s c )
R B ( 2 n , n ) / p 3 = 1 4 i = 1 2 t r i ( τ + ( 2 i 1 ) T s c T s c ) t r i ( τ ( 2 i 1 ) T s c + 2 T s c )
R B ( 2 n , n ) / p 4 = 1 4 i = 1 2 t r i ( τ + ( 2 i 1 ) T s c ) t r i ( τ ( 2 i 1 ) T s c + 3 T s c )

2.3. Correlation Function Combination

By analyzing the sub-correlation functions of the BOC(kn, n) signal with each modulation factor, it is found that the main peak position of the correlation function of the first segment corresponds to the main peak position of the correlation function of the M-th segment, and the side peak corresponds to the zero points. No matter what value M is, the product of R B / P 1 and R B / P M can yield non-zero value only near the (0,0) point. And R B / P M can be obtained by R B / P 1 shift ( M 1 ) T c / M and multiply by −1. Thus, combining correlation functions, such as formulas
R 0 = R B / P 1 · R B / P M
It can be found that R 0 has two small sub-peaks with a negative value and a positive main peak. To achieve the purpose of eliminating sub-peaks, the combination algorithm is
R S R S A = | | R B / P 1 · R B / P M | + R B / P 1 · R B / P M |
Figure 9 is a simulation diagram of the combination function for M = 2 and M = 4. It can be seen from the graph that the larger M is, the smaller the main peak width of the combination function is.

3. Acquisition Structure Design and Performance Analysis Based on SRSA Principle

3.1. The Acquisition Structure Diagram

When the receiver generates a local signal, the local PRN sequence is separated to obtain a sub-signal. The sub-signal is correlated with the received signal to obtain a sub-correlation function, and the sub-correlation function is combined according to Equation (17) to get a correlation function of the no secondary peak. Figure 10 shows the SRSA acquisition schematic. It is worth noting that the sub-correlation functions can be obtained from each other. According to this property, only one sub-correlation function can be generated, and the generated sub-correlation function can be delayed to get other functions.
The steps for signal acquisition are as follows:
Step 1: The carrier NCO generates two orthogonal local carriers, and mixes with the received intermediate frequency (IF) BOC(kn, n) signal with the purpose of cancelling the carrier.
Step 2: Generate a local PRN code. And then generate a local auxiliary code according to the method in Section 2.2.
Step 3: The auxiliary code is coherently integrated with the input signal to obtain a sub-correlation function R B / P 1 .
Step 4: Delay R B / P 1 by sample points of ( M 1 ) / M PRN chip length, and a new sub-function R B / P M is obtained.
Step 5: According to the Formula (15), the final algorithm operation is performed to obtain a correlation function without a sub-peak, that is, the detection value V.

3.2. Detection Performance Analysis

The input BOC signal is expressed as follows:
S ( t ) = 2 A × C ( t τ ) × D ( t τ ) × S C ( t τ ) × cos ( 2 π ( f I F + Δ f D ) t + φ ) + N n
2 A is the amplitude of the input signal, C ( t ) is the PRN code, D ( t ) is the navigation data, τ is the code delay of the input signal, Δ f D is the Doppler frequency of the input signal, f I F is the intermediate frequency, S c ( t ) is the subcarrier, φ is the initial phase of carrier and N n is the noise term.
Mixing input signal with local carrier. For simplicity, we exploit the baseband complex envelope expression of S O ( t ) .
S O ( t ) = 2 A × C ( t τ ) × D ( t τ ) × S C ( t τ ) × e j ( 2 π Δ f D t + φ ) + N c + j N s
N c and N s are independent Gaussian random processes. They have the same zero mean and double-sided power spectrum density N 0 . Then, multiply the local auxiliary signal according to the SRSA principle.
S 1 ( t )   = 2 A × C ( t τ ) × D ( t τ ) × S C ( t τ ) × S P R N , 1 ( t ) × e j ( 2 π Δ f D t + φ ) + N c + j N s
S M ( t )   = 2 A × C ( t τ ) × D ( t τ ) × S C ( t τ ) × S P R N , M ( t ) × e j ( 2 π Δ f D t + φ ) + N c + j N s
Then perform coherent integration, D ( t ) is a constant, does not affect the overall process, omitted here, no need to consider the change of navigation data bits.
S ¯ 1 = A × R B / P , 1 ( τ ) × sin c ( π Δ f D T c o h ) × e j ( 2 π Δ f D T c o h + φ ) + N 1                     A × R B / P , 1 ( τ ) × sin c ( π Δ f D T c o h ) + N 1                     = S ¯ 1 ( τ , Δ f D ) + N 1
S ¯ M = A × R B / P , M ( τ ) × sin c ( π Δ f D T c o h ) × e j ( 2 π Δ f D T c o h + φ ) + N m                           A × R B / P , M ( τ ) × sin c ( π Δ f D T c o h ) + N m                           = S ¯ M ( τ , Δ f D ) + N m
When the carrier and code loop tracking is successful, the complex exponential is almost 1. S ¯ 1 ( τ , Δ f D ) is the signal part of S ¯ 1 , S ¯ M ( τ , Δ f D ) is the signal part of S ¯ M , and the non-coherent accumulation time is M n o h . In digital signal processing and communication theory, the normalized sin c function is usually defined as sin c ( x ) = sin ( π x ) / π x . According to the proposed combination principle, the final detection value is
S f i n a l = i = 1 M n o h ( | S 1 ¯ × S M ¯ | + S 1 ¯ × S M ¯ ) 2 i = 1 M n o h ( | S 1 ¯ × S M ¯ | ) = 2 i = 1 M n o h ( | S 1 ¯ ( τ , Δ f D ) + N 1 | × | S M ¯ ( τ , Δ f D ) + N m | ) = 2 i = 1 M n o h ( | S 1 ¯ ( τ , Δ f D ) | × | S M ¯ ( τ , Δ f D ) | ) + 2 i = 1 M n o h ( | S 1 ¯ ( τ , Δ f D ) | × N m ) + 2 i = 1 M n o h ( | S M ¯ ( τ , Δ f D ) | × N 1 ) +   2 i = 1 M n o h ( N 1 × N m )
The main peak of the correlation function is much higher than the side peak. | S 1 ¯ × S M ¯ | + S 1 ¯ × S M ¯ can be approximated as 2 | S 1 ¯ × S M ¯ | . Assuming that X 0 is the case that only noise presents in the received signal, X 1 is the other case that both signal and noise exist, here the two cases are analysed [12]. When the target signal does not exist, it can also be said that the received signal is phase-aligned with the local signal code outside of one chip. At this time, the detection statistic value V 1 obeys Rayleigh distributions. Suppose the probability density function (PDF) at this time is P V 1 , and the false alarm probability can be expressed as
P f a = P [ V 1 > V ] = V + P V 1 ( x ) d x = V + V 1 ( σ n ) 2 · exp { V 1 2 2 ( σ n ) 2 } d V 1
Therefore, according to Equation (25), the threshold value V can be reversed after the false alarm probability and noise power are given. σ n 2 = C / ( 10 0.1 dB ) / 0.004 . C is carrier power of the input (IF) signal. d B is the carrier-to-noise ratio range.
When the received signal is aligned with the locally generated chip, the detected value obeys the Rice distribution, and its mean and variance are
E ( V 2 ) = 4 M n o h A 2 T c o h 2 | R B / P 1 ( τ ) | | R B / P M ( τ ) |
D ( V 2 ) = 16 M n o h σ 4 + l = 0 M n o h [ 4 A 2 T c o h 2 ( | R B / P 1 ( τ ) | l × | R B / P M ( τ ) | l ) ]
when M n o h is large, the detection statistic V 2 is approximately obeying the Rice distribution.
P V 2 ( x ) is assumed to be the probability density at this time, and the detection probability P d is:
P d = V + P V 2 ( x ) d x
P V 2 ( x ) = V 2 σ 2 exp { V 2 2 + a 2 2 σ 2 } I 0 ( a V 2 σ 2 )
a 2 / σ 2 is the signal-to-noise ratio (SNR), I 0 ( x ) is the Modified Bessel function of the first kind. The detection statistics can only be maximized when the code phase error τ is close to zero.

3.3. Performance Simulation and Analysis

Based on the Matlab platform, the intermediate frequency of the input signal is 4.092 MHz, the sampling rate is set to 40.92 MHz, the code phase offset is in the 601th sampling point, the Doppler is set to 2000 Hz. The coherent integration time is T c o h = 1   ms . The four acquisition methods of SRSA, SCPC, ASPeCT and BPSK-LIKE are compared to the performance.

3.3.1. Disambiguation Performance Analysis

Figure 11 and Figure 12 shows a normalized two-dimensional acquisition comparison of the four methods. When M = 2, it can be found that the SRSA method is better than the SCPC, ASPeCT, and BPSK-LIKE methods. In terms of main peak span, the SRSA main peak spans 20 sampling points (half chip points) and is the same as the ASPeCT method, and is better than the BPSK-LIKE method and the SCPC method 80 sampling points (i.e., two-chip points). For the peak value, the SRSA method is the same as the ASPeCT and SCPC methods and 10% higher than the BPSK-LIKE method. And the SRSA method removes the effects of side peaks while the ASPeCT still has two smaller sub-peaks. When M = 4, the SCPC and BPSK-LIKE methods remove the sub-peaks, but still do not retain the narrow correlation characteristic of BOC signal. The ASPeCT method has more sub-peaks of the correlation function when M = 4, which is not applicable to BOC signals above the 2nd order. And the SRSA can still completely remove all the sub-peaks of the 4th-order BOC, and the main peak span is only ten sampling points.

3.3.2. Computational Complexity Analysis

The one millisecond signal is received in the pre-set simulation environment, and the number of frequency searches is 41 times. Using the FFT parallel code phase search method, the number of complex additions required to complete a 40,920 points FFT is:
n log 2 n = 40 , 920 log 2 40 , 920 = 626 , 916
The number of complex multiplications required is
n 2 log 2 n = 40 , 920 2 log 2 40 , 920 = 313 , 458
One complex multiplication is equal to four real multiplications and two real additions, and a complex addition equals two real additions [13,14]. Table 1 is the computational complexity of the four algorithms in the acquisition. It shows from the table that the computational complexity of SRSA is the same as ASPeCT, and is better than BPSK-Like and SCPC.

3.3.3. Acquisition Probability Analysis

The Monte Carlo method is used to simulate the detection probability of BOC(1, 1) and BOC(10, 5) under different acquisition methods. Set the carrier-to-noise ratio range of the simulation environment to [20-50], set the false alarm probability, the non-coherent accumulation number is L = 10, and perform 2000 statistics as shown in Figure 13 and Figure 14.
As can be seen from the above figure, the SRSA method has obvious advantages when the acquisition probability is higher than 0.9. For 2nd order BOC signals, the SRSA acquisition probability is better than the BPSK-LIKE, SCPC and ASPeCT methods 2.5 dB-Hz, 1 dB-Hz, and 0.5 dB-Hz. For high-order BOC signals, the ASPeCT method cannot suppress the sub-peaks of the correlation function, and the effect is not satisfactory. There is no comparison here. The SRSA acquisition probability is better than the BPSK-LIKE and SCPC methods 0.5 dB-Hz and 0.25 dB-Hz when the acquisition probability is higher than 0.9.

3.3.4. Peak-to-Average Ratio Analysis

The peak-to-average ratio is also an important indicator for evaluating the performance of an algorithm, and is defined as the ratio of the peak value to the average value of the correlation result obtained by the acquisition. In the preset simulation environment, the BPSK-LIKE, SCPC, ASPeCT and SRSA methods were compared by repeating 3000 experiments, and the peak-to-average ratios of the four methods at different signal-to-noise ratios were obtained, as shown in Figure 15 and Figure 16.
When M = 2, it can be seen that the performance of the SRSA method is optimal. At −20 dB, the value of the SRSA method is 3.5 times that of the BPSK-LIKE method, three times that of the SCPC method, and twice that of the ASPeCT method. The values of the four methods SRSA, ASPeCT, SCPC, and BPSK-LIKE are 280.6, 146.6, 93.27 and 81.17 at −20 dB.
When M = 4, the curve of the SRSA method intersects the curve of the ASPeCT method at −24 dB. The curve of the SCPC method and the BPSK-LIKE method almost coincide. At −20 dB, the values of the four methods SRSA, ASPeCT, SCPC, and BPSK-LIKE are 150.6, 150.4, 90.14 and 90.14.

4. Tracking Loop Structure Design and Performance Analysis Based on SRSA Principle

4.1. Discriminator Function and Tracking Loop Structure

Using the first and the M-th sub-PRN codes as two local auxiliary codes in the tracking loop. Without considering multipath and interference, the BOC IF signal received by a GNSS receiver from a satellite can be expressed as [15]:
r ( t ) = 2 P · D ( t ) · S B O C ( t ) · cos ( 2 π f I F t + φ 0 ) + n ( t )
Among them, P is the total power of the received signal, D ( t ) is the navigation data code, S B O C ( t ) is defined in Equation (6), f I F is the intermediate frequency after down-converting the received signal, φ 0 is the initial phase of the carrier. n ( t ) is Band-limited white noise [16], and it has the zero mean and double-sided power spectrum density N 0 .
After the carrier is cancelled, the I and Q signals perform coherent integration operations with the Early and Late local auxiliary codes, and the following is obtained:
{ I E 1 + j Q E 1 = 2 P R B / P 1 ( Δ τ d 2 ) e j Δ φ + n I E 1 + j n Q E 1 I E M + j Q E M = 2 P R B / P M ( Δ τ d 2 ) e j Δ φ + n I E M + j n Q E M I L 1 + j Q L 1 = 2 P R B / P 1 ( Δ τ d 2 ) e j Δ φ + n I L 1 + j n Q L 1 I L M + j Q L M = 2 P R B / P M ( Δ τ d 2 ) e j Δ φ + n I L M + j n Q L M
In the formula, the subscripts 1 and M indicate the local code number, and the subscripts E and L respectively indicate the Early and Late branches, and other branches are similarly represented. Δ τ and Δ φ respectively represent the estimation error of the code phase and the initial phase of the carrier, d is the delay spacing between the early and late correlators. All n are noise terms obeying the Gaussian distribution [16,17].
Combined with the Formulas (26) and (27), The distribution function of the outputs I E 1 , I E M , I L 1 , I L M at Δ τ = 0 is
( I E 1 , I E M , I L 1 , I L M ) T ~ N ( μ cos ( Δ φ ) , ω ) ( Q E 1 , Q E M , Q L 1 , Q L M ) T ~ N ( μ sin ( Δ φ ) , ω )
with
μ = 2 P · sin c ( π Δ f T c o h ) · [ R B / P 1 ( d / 2 )           R B / P 1 ( d / 2 )           R B / P M ( d / 2 )           R B / P M ( d / 2 ) ] T
ω = N 0 T c o h [ R B ( 0 ) R B ( d ) R B / P 1 ( 0 ) R B / P 1 ( d ) R B ( d ) R B ( 0 ) R B / P 1 ( d ) R B / P 1 ( 0 ) R B / P M ( 0 ) R B / P M ( d ) R P 1 / P M ( 0 ) R P 1 / P M ( d ) R B / P M ( d ) R B / P M ( 0 ) R P 1 / P M ( d ) R P 1 / P M ( 0 ) ]
The output of Q branch and I branch are independent of each other. The distribution function of the outputs Q E 1 , Q E M , Q L 1 , Q L M at Δ τ = 0 is N ( 0 , ω ) . R P 1 / P M is the cross-correlation function of the auxiliary signal.
According to Equations (17), the formula for the incoherent discriminator function is:
D S R S A ( Δ τ ) = ( | ( I E 1 + Q E 1 ) · ( I E M + Q E M ) | + ( I E 1 + Q E 1 ) · ( I E M + Q E M ) ) ( | ( I L 1 + Q L 1 ) · ( I L M + Q L M ) | + ( I L 1 + Q L 1 ) · ( I L M + Q L M ) ) = 2 P ( R S R S A 2 ( Δ τ d / 2 ) R S R S A 2 ( Δ τ + d / 2 ) )
Figure 17 shows the discriminator curve for the proposed method at different delay spacings. Compared with Figure 4, the proposed method successfully removes the mis-lock point.
Figure 18 shows the proposed code tracking loop structure. The received intermediate frequency signal is first subjected to carrier cancellation operation. At the same time, the code loop generates two local auxiliary codes and performs Early and Late delays. The received signal is correlated with the local code. Then the local code phase is adjusted by the phase detector and the NCO (Numerically Controlled Oscillator) to complete the unambiguous tracking of the BOC signal.

4.2. The Code Tracking Error Analysis

The code tracking error formula is given in [18].
σ c t e 2 = 2 B L ( 1 0.5 B L T c o h ) T c o h σ v 2 K v 2
where B L is the code loop filter bandwidth, σ v is the discriminator output standard deviation, and K v is the discriminator gain [19].
K V = d D S R S A d Δ τ | Δ τ = 0 = 16 π P [ f H ( f ) G B / P 1 ( f ) sin ( π f d ) d f f H ( f ) G B / P M ( f ) cos ( π f d ) d f f H ( f ) G B / P M ( f ) sin ( π f d ) d f f H ( f ) G B / P 1 ( f ) cos ( π f d ) d f ]
In Equation (39), H ( f ) is the power spectral density of the front-end filter, and G B / P represents the cross-power spectral density between the BOC signal and the auxiliary signal. For convenience, let
{ V I E = | I E 1 I E M | + I E 1 I E M V I L = | I L 1 I L M | + I L 1 I L M V Q E = | Q E 1 Q E M | + Q E 1 Q E M V Q L = | Q L 1 Q L M | + Q L 1 Q L M
then
D S R S A = V I E V I L + V Q E V Q L
when the tracking process is in a steady state, considering that Δ f 0 , Δ φ 0 and Δ τ 0 , we have
E ( D S R S A ) = 0 E [ ( V I E V I L ) · ( V Q E V Q L ) ] = 0 σ V 2 = E ( D S R S A 2 ) = E [ ( V I E V I L ) 2 ] + E [ ( V Q E V Q L ) 2 ]
As long as E [ ( V I E V I L ) 2 ] and E [ ( V Q E V Q L ) 2 ] are solved, σ V 2 can be obtained. It can be calculated using the eigenfunction of the joint Gaussian distribution.
E [ ( V I E V I L ) 2 ] = 2 E [ ( I E 1 I E M ) 2 ] + 2 E [ ( I L 1 I L M ) 2 ] 2 E [ | I E 1 I E M | · ( I E 1 I E M ) ]   2 E [ | I L 1 I L M | · ( I L 1 I L M ) ] 2 E [ | I E 1 I E M | · | I L 1 I L M | ]   2 E [ ( I E 1 I E M ) · ( I L 1 I L M ) ]         + 2 E [ ( I E 1 I E M ) · | I L 1 I L M | ] + 2 E [ | I E 1 I E M | · ( I L 1 I L M ) ]       2 E [ ( I E 1 I E M ) 2 ] + 2 E [ ( I E 1 I E M ) ]   2       = 2 [ ( 1 + ρ 0 2 ) σ 0 4 + cos 2 ( Δ φ ) ( μ 2 ( 1 ) + μ 2 ( 3 ) + 2 μ ( 1 ) μ ( 3 ) ρ 0 ) σ 0 2 ]                              
E [ ( V Q E V Q L ) 2 ] 2 E [ ( Q E 1 Q E M ) 2 ] + 2 E [ ( Q E 1 Q E M ) ]   2       = 2 [ ( 1 + ρ 0 2 ) σ 0 4 + sin 2 ( Δ φ ) ( μ 2 ( 1 ) + μ 2 ( 3 ) + 2 μ ( 1 ) μ ( 3 ) ρ 0 ) σ 0 2 ]
In Equations (43) and (44), in the case where R ^ B ( 0 ) is approximately equal to R ^ P 1 / P M ( 0 ) , ρ 0 = R ^ B / P M 2 ( 0 ) / R ^ B / P 1 2 ( 0 ) , σ 2 = N 0 T c o h R ^ B ( 0 ) .
Substituting Equations (39), (43), and (44) into Equation (38), the code tracking error variance of the proposed algorithm can be derived.

4.3. Performance Simulation and Analysis

Based on the Matlab platform, the intermediate frequency of the input signal is 4.092 MHz, the sampling rate is set to 40.92 MHz, the code phase offset is in the 609th sampling point, the Doppler search range is [−10 kHz, 10 kHz], and the step is 500 Hz. The four methods of SRSA, SCPC, ASPeCT and BPSK-LIKE are compared to the performance.

4.3.1. Phase Discrimination Curve Analysis

The classical EMLP (Early-Minus-Late-Power) phase detector is often used to generate the phase discrimination curve to evaluate the tracking performance of each method [11]. Assuming that the front-end bandwidth is infinite, compare the traditional Delay Lock Loop (DLL), BPSK-like, ASPeCT, and SRSA methods. The results show that for the second-order BOC signal (Figure 19), the phase discrimination curve of the traditional DLL has two mis-locking points, and the other three algorithms can remove the mis-locking point. The stable region is [−0.1 Tc, +0.1 Tc]. However, the slope gain of the phase discrimination curve of the DLL, ASPeCT, and SRSA methods relative to the linear region of the phase discrimination curve of BPSK-like is 5.2 dB. For the 4th-order BOC signal (Figure 20), the phase discrimination curve of the traditional DLL has six mis-locking points. The stability of the other three algorithms is [−0.05 Tc, +0.05 Tc], but the ASPeCT method has four mis-locking points, which cannot remove the tracking ambiguity. The proposed SRSA method can completely remove the mis-locking point, and the slope gain is 7.2 dB compare with BPSK-like, and the slope of the phase discrimination curve can be kept large. The phase discrimination error output is proportional to the slope of the phase discrimination curve at the center zero-crossing point, and the performance such as anti-noise and tracking jitter accuracy is closely related to this.

4.3.2. Anti-Multipath Performance Analysis

Multipath Error Envelope (MEE) is a typical indicator for evaluating the multipath performance of tracking loops, reflecting the sensitivity of a code tracking loop to different parameters of multipath signals [20,21]. Figure 21 and Figure 22 show the multipath envelope comparison of BOC(1, 1) and BOC(10, 5) obtained by the four tracking algorithm when the correlator interval is set to 0.1 Tc and 0.05 Tc. The envelope extreme value (the absolute maximum value of the MEE), the length of the envelope interval length (the sum of the abscissa intervals when the MEE takes a non-zero error), and the envelope area (the area enclosed by the MEE) are three measures of anti-multipath performance [22]. The smaller the three indicators, the better the anti-multipath performance. It can be seen from the Figure 21 and Figure 22 that the SRSA and ASPeCT methods have better anti-multipath performance under short multipath delay conditions. Still, ASPeCT does not completely remove the sub-peaks and causes errors in a certain multipath delay range. Therefore, the method proposed in this paper has excellent anti-multipath performance.

4.3.3. Anti-Noise Performance Analysis

Thermal noise is another important cause of tracking error, and loop code tracking error is an important indicator to measure the anti-noise performance of tracking methods [23]. As shown in Figure 23, under different carrier-to-noise ratios, the single-sideband loop bandwidth is set to B = 2   Hz , the coherent integration time is T c o h = 1   ms , and the tracking interval is 0.1   Tc . The standard deviation of the code tracking error is obtained by tracking the BOC(1,1) using several tracking methods. The code tracking error of ASPeCT and SRSA is improved compared with the traditional EMLP phase discriminator at a low carrier-to-noise ratio, and the BPSK-LIKE method has poor code tracking error. And compared with BPSK-LIKE, the standard deviation of code tracking error is reduced by 0.045   Tc . As shown in Figure 24, under different carrier-to-noise ratios, the single-sideband loop bandwidth is set to B = 2   Hz , the coherent integration time is T coh = 1   ms , and the tracking interval is 0.05   Tc . The code tracking error standard deviation is obtained by tracking the BOC(10,5) by several tracking methods. Because ASPeCT is only applicable to 2nd-order BOC signals, only BPSK-LIKE, EMLP, and SRSA are analyzed here. The performance of the SRSA method is close to that obtained by EMLP, compared with BPSK-LIKE, the standard deviation of code tracking error is reduced by 0.03   Tc . The BPSK-LIKE method is still poor.

5. Conclusions

Aiming at the ambiguity of BOC signal acquisition/tracking, many methods have been proposed [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31]. Sub-function reconstruction synchronization algorithm is proposed in this paper. By designing two auxiliary codes and performing cross-correlation operation with BOC signal, and then obtaining a single peak correlation function by Formula (17), we achieved the purpose of eliminating ambiguity. In the acquisition stage, the acquisition probability, computational complexity, and peak-to-average ratio were evaluated. Compared with ASPeCT, BPSK-LIKE, and SCPC, the SRSA method has excellent performance in the acquisition process. The SRSA method preserves the narrow correlation of the BOC signal correlation function and also removes the interference of the sub-peak. In the tracking stage, the tracking loop uses incoherent discriminator function (36) to eliminate ambiguity. The SRSA method removes the mis-locking point. In terms of anti-multipath performance, the SRSA method has the smallest envelope interval length and envelope area compared to ASPeCT and BPSK-LIKE. In terms of anti-noise performance, the SRSA method has almost the same code tracking error as the ASPeCT method and is superior to the BPSK-LIKE method. In summary, the proposed method can completely remove the threat of ambiguity while achieving significant performance advantages in terms of acquisition and tracking.
The main contributions of this paper are as follows:
(1)
A sub-function reconstruction synchronization algorithm is proposed. Design auxiliary code and perform cross-correlation operation with BOC signal to obtain cross-correlation functions, as shown in Figure 7 or Figure 8. The first sub-function is at the same position as the main peak of the last sub-function and the sub-peak corresponds to 0 point. Based on this characteristic, combining Equation (17), we can obtain a new function with a single peak.
(2)
The proposed method can be applied to BOC signals of arbitrary modulation order.
(3)
Reduced computational complexity and resource consumption.

Author Contributions

All authors contributed to the manuscript and discussed the results. All authors together developed the idea that led to this paper. X.G. and F.H. conceived the experiments and analyzed the data. L.H. and R.J. performed data processing and contributed to the manuscript. B.Y. and H.Z. provided critical comments and contributed to the final revision of the paper. B.W. contributed to the expression and the design of programs. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Key Research and Development Plan of China (project: Indoor Hybrid Intelligent Positioning and Indoor GIS Technology (No. 2016YFB0502100, 2016YFB0502101)).

Acknowledgments

This work was supported in part by the National Key Research and Development Plan of China (project: Indoor Hybrid Intelligent Positioning and Indoor GIS Technology (No. 2016YFB0502100, 2016YFB0502101)).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shen, F.; Xu, G.; Li, Q. Non-Coherent Unambiguous Tracking Method for Cosine-BOC Signals Based on an S-Curve Shaping Technique. IEEE Signal Process. Lett. 2015, 22, 752–756. [Google Scholar] [CrossRef]
  2. Hodgart, M.S.; Blunt, P.D.; Unwin, M. The optimal dual estimate solution for robust tracking of binary offset carrier (BOC) modulation. In Proceedings of the ION GNSS 2007, Fort Worth, TX, USA, 25–28 September 2007; pp. 1017–1027. [Google Scholar]
  3. Hodgart, M.S.; Simons, E. Improvements and additions to the double estimation technique. In Proceedings of the Satellite Navigation Technologies and European Workshop on GNSS Signals and Signal Processing (NAVITEC), Noordwijk, The Netherlands, 5–7 December 2012. [Google Scholar]
  4. Fine, P.; Wilson, W. Tracking algorithm for GPS offset carrier signals. In Proceedings of the ION NTM 1999, San Diego, CA, USA, 25–27 January 1999; pp. 671–676. [Google Scholar]
  5. Burian, A.; Lohan, E.S.; Renfors, M. BPSK-like Methods for Hybrid-Search Acquisition of Galileo Signals. In Proceedings of the IEEE International Conference on Communications, Istanbul, Turkey, 11–15 June 2006; pp. 5211–5216. [Google Scholar]
  6. Julien, O.; Macabiau, C.; Cannon, M.E.; Lachapelle, G. ASPeCT: Unambiguous sine-BOC(n,n) acquisition/tracking technique for navigation applications. IEEE Trans. Aerosp. Electron. Syst. 2011, 43, 150–162. [Google Scholar] [CrossRef] [Green Version]
  7. Liu, Z.; Xu, B.; Tang, X. Multipath mitigating technique for BOC(kn, n) signal in GNSS. In Proceedings of the International Conference on Wireless Communications and Signal Processing, Nanjing, China, 15–17 October 2015; pp. 1–4. [Google Scholar]
  8. Yao, Z.; Lu, M.; Feng, Z. Unambiguous sine-phased binary offset carrier modulated signal acquisition technique. IEEE Trans. Wirel. Commun. 2010, 9, 577–580. [Google Scholar] [CrossRef]
  9. Yao, Z.; Cui, X.; Lu, M.; Feng, Z. Pseudo-correlation-function-based unambiguous tracking technique for sine-BOC signals. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 1782–1796. [Google Scholar] [CrossRef]
  10. Feng, S.; Xu, G.H.; Fang, H.Y. Synthesized correlation function based unambiguous acquisition technique for sin-BOC/MBOC modulated signals. J. Syst. Eng. Electron. 2015, 37, 1980–1986. [Google Scholar]
  11. Yan, T.; Wei, J.; Tang, Z.; Qu, B.; Zhou, Z. Unambiguous Acquisition/Tracking Technique for High-Order Sine-Phased Binary Offset Carrier Modulated Signal. Wirel. Pers. Commun. 2015, 84, 2835–2857. [Google Scholar] [CrossRef]
  12. Ji, Y.-F.; Liu, Y.; Zhen, W.-M.; Sun, X.-Y.; Yu, B.-G. An unambiguous acquisition algorithm based on unit correlation for BOC(n,n) signal. IEICE Trans. Commun. 2017, 8, 1507–1513. [Google Scholar] [CrossRef]
  13. Wang, C.; Yang, C.; Ku, S.; Liu, C. A Low-Complexity Peak-to-Average Power Ratio Estimation Method for OFDM Signals. In Proceedings of the IEEE Vehicular Technology Conference, Montreal, QC, Canada, 25–28 September 2006. [Google Scholar] [CrossRef]
  14. Wu, J.; Jiang, L.; Chen, Y.; Senhadji, L.; Shu, H. A generalized modified split-radix FFT algorithm for N = q × 2m and its applications. In Proceedings of the IEEE 2014 7th International Congress on Image and Signal Processing, Dalian, China, 14–16 October 2014. [Google Scholar] [CrossRef]
  15. Lei, D.; Lu, W.; Yu, D. Resource-efficient acquisition architecture for BOC-modulated signals. IEICE Electron. Express 2015, 11, 20140358. [Google Scholar] [CrossRef] [Green Version]
  16. Zhou, Y.; Hu, X.; Ke, T.; Tang, Z. Ambiguity Mitigating Technique for Cosine-Phased Binary Offset Carrier Signal. IEEE Trans. Wirel. Commun. 2012, 11, 1981–1984. [Google Scholar] [CrossRef]
  17. Yan, T.; Wei, J.; Tang, Z.; Qu, B.; Zhou, Z. Unambiguous combined correlation functions for sine-BOC signal tracking. GPS Solut. 2015, 19, 623–638. [Google Scholar] [CrossRef]
  18. Kao, T.L.; Juang, J.C. Weighted discriminators for GNSS BOC signal tracking. GPS Solut. 2012, 16, 339–351. [Google Scholar] [CrossRef]
  19. Li, T.; Tang, Z.; Wei, J.; Zhou, Z.; Wang, B. Unambiguous Tracking Technique Based on Combined Correlation Functions for Sine BOC Signals. J. Navig. 2019, 72, 140–154. [Google Scholar] [CrossRef]
  20. Mao, H.; Wu, D. Evaluation of Anti-multipath Performance to Satellite Navigation Signal Tracking on L1/E1/B1 Frequency Band. J. Syst. Simul. 2016, 28, 106–113. [Google Scholar]
  21. Deng, Z.; Hu, E.; Yin, L.; Liu, W. An Unambiguous Tracking Technique for Sine-BOC(kn,n) Modulated GNSS Signals. Wirel. Pers. Commun. 2018, 103, 1101–1112. [Google Scholar] [CrossRef]
  22. Liu, H.C.; Cheng, X.; Ni, S.J. Evaluation of Multipath Mitigation Performances Based on Error Envelope. J. Natl. Univ. Def. Technol. 2011, 33, 72–75. [Google Scholar]
  23. Wu, J.; Dempster, A.G. Code Tracking Variance Analysis for GNSS Receivers with “Strobe Correlators”. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 2760–2782. [Google Scholar] [CrossRef]
  24. Borio, D. Coherent side-band BOC processing. IET Radar Sonar Navig. 2017, 11, 1455–1466. [Google Scholar] [CrossRef]
  25. Feng, T.; Kai, Z.; Liang, C. Unambiguous Tracking of BOC Signals Using Coherent Combination of Dual Sidebands. Comput. Educ. 2016, 20, 1555–1558. [Google Scholar] [CrossRef]
  26. Li, T.; Wei, J.L.; Tang, Z.P.; Zhou, Z.H.; Wang, B.Y. An Optimizing Combined Unambiguous Correlation Functions for BOC Signals Tracking. In Proceedings of the 2017 International Technical Meeting of the Institute of Navigation, Washington, DC, USA, 30 January–2 February 2017; pp. 388–400. [Google Scholar]
  27. Zhang, H.; Ba, X.; Chen, J.; Zhou, H. Unambiguous acquisition technique for BOC(m,n) modulated signals. Acta Aeronaut. Astronaut. Sin. 2017, 38, 320394. [Google Scholar]
  28. Shim, D.-S.; Jeon, J.-S. An Unambiguous Delay-And-Multiply Acquisition Scheme for GPS L1C Signals. Sensors 2018, 18, 1739. [Google Scholar] [CrossRef] [Green Version]
  29. Xu, C.; Liu, Z.; Tang, X.; Wang, F. A Design Method of Code Correlation Reference Waveform in GNSS Based on Least-Squares Fitting. Sensors 2016, 16, 1194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Zhen, L.; Jie, H.; Wang, J.; Zhao, Y.; Chen, S.W. Generalized unambiguous tracking method based on pseudo correlation function for multi-level coded symbol modulated signals. Acta Phys. Sin. 2017, 66, 139101. [Google Scholar]
  31. Zhen, L.; Zhang, J.Y.; Lu, M.Q.; Jie, H.; Zhao, Y.J. Universal evaluation criteria for code delay estimation error of satellite navigation signals. Acta Phys. Sin. 2017, 66, 129101. [Google Scholar]
Figure 1. Auto-correlation function.
Figure 1. Auto-correlation function.
Sensors 20 00485 g001
Figure 2. Cross-correlation functions.
Figure 2. Cross-correlation functions.
Sensors 20 00485 g002
Figure 3. The false acquisition probability of three different factors.
Figure 3. The false acquisition probability of three different factors.
Sensors 20 00485 g003
Figure 4. The discrimination curves for the three methods.
Figure 4. The discrimination curves for the three methods.
Sensors 20 00485 g004
Figure 5. The Separation of PRN at M = 2.
Figure 5. The Separation of PRN at M = 2.
Sensors 20 00485 g005
Figure 6. The Separation of PRN at M = 4.
Figure 6. The Separation of PRN at M = 4.
Sensors 20 00485 g006
Figure 7. The normalized sub cross-correlation functions at M = 2.
Figure 7. The normalized sub cross-correlation functions at M = 2.
Sensors 20 00485 g007
Figure 8. The normalized sub cross-correlation functions at M = 4.
Figure 8. The normalized sub cross-correlation functions at M = 4.
Sensors 20 00485 g008
Figure 9. The combination results of correlation functions by SRSA.
Figure 9. The combination results of correlation functions by SRSA.
Sensors 20 00485 g009
Figure 10. The SRSA acquisition schematic.
Figure 10. The SRSA acquisition schematic.
Sensors 20 00485 g010
Figure 11. The 2-D acquisition result at M = 2.
Figure 11. The 2-D acquisition result at M = 2.
Sensors 20 00485 g011
Figure 12. The 2-D acquisition result at M = 4.
Figure 12. The 2-D acquisition result at M = 4.
Sensors 20 00485 g012
Figure 13. The acquisition probability of BOC(1, 1).
Figure 13. The acquisition probability of BOC(1, 1).
Sensors 20 00485 g013
Figure 14. The acquisition probability of BOC(10, 5).
Figure 14. The acquisition probability of BOC(10, 5).
Sensors 20 00485 g014
Figure 15. The peak to average ratio at M = 2.
Figure 15. The peak to average ratio at M = 2.
Sensors 20 00485 g015
Figure 16. The peak to average ratio at M = 4.
Figure 16. The peak to average ratio at M = 4.
Sensors 20 00485 g016
Figure 17. The discriminator curve at different delay spacings of BOC(1, 1).
Figure 17. The discriminator curve at different delay spacings of BOC(1, 1).
Sensors 20 00485 g017
Figure 18. The proposed code tracking loop.
Figure 18. The proposed code tracking loop.
Sensors 20 00485 g018
Figure 19. Discriminator output curve of BOC(1, 1).
Figure 19. Discriminator output curve of BOC(1, 1).
Sensors 20 00485 g019
Figure 20. Discriminator output curve of BOC(10,5).
Figure 20. Discriminator output curve of BOC(10,5).
Sensors 20 00485 g020
Figure 21. Multipath error envelope of BOC(1, 1).
Figure 21. Multipath error envelope of BOC(1, 1).
Sensors 20 00485 g021
Figure 22. Multipath error envelope of BOC(10,5).
Figure 22. Multipath error envelope of BOC(10,5).
Sensors 20 00485 g022
Figure 23. Code tracking error standard deviation of BOC(1,1).
Figure 23. Code tracking error standard deviation of BOC(1,1).
Sensors 20 00485 g023
Figure 24. Code tracking error standard deviation of BOC(10, 5).
Figure 24. Code tracking error standard deviation of BOC(10, 5).
Sensors 20 00485 g024
Table 1. Comparison calculation between four algorithms.
Table 1. Comparison calculation between four algorithms.
NumberAlgorithm NameMultiplication/TimesAddition/Times
1BPSK-Like10Nlog2N + 10N17.5Nlog2N + 7N
2ASPeCT10Nlog2N + 9N17.5Nlog2N + 7N
3SCPC10Nlog2N + 10N17.5Nlog2N + 7N
4SRSA10Nlog2N + 9N17.5Nlog2N + 7N

Share and Cite

MDPI and ACS Style

Hao, F.; Yu, B.; Gan, X.; Jia, R.; Zhang, H.; Huang, L.; Wang, B. Unambiguous Acquisition/Tracking Technique Based on Sub-Correlation Functions for GNSS Sine-BOC Signals. Sensors 2020, 20, 485. https://doi.org/10.3390/s20020485

AMA Style

Hao F, Yu B, Gan X, Jia R, Zhang H, Huang L, Wang B. Unambiguous Acquisition/Tracking Technique Based on Sub-Correlation Functions for GNSS Sine-BOC Signals. Sensors. 2020; 20(2):485. https://doi.org/10.3390/s20020485

Chicago/Turabian Style

Hao, Fang, Baoguo Yu, Xingli Gan, Ruicai Jia, Heng Zhang, Lu Huang, and Boyuan Wang. 2020. "Unambiguous Acquisition/Tracking Technique Based on Sub-Correlation Functions for GNSS Sine-BOC Signals" Sensors 20, no. 2: 485. https://doi.org/10.3390/s20020485

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