Next Article in Journal
Development and Experimental Validation of a Dry Non-Invasive Multi-Channel Mouse Scalp EEG Sensor through Visual Evoked Potential Recordings
Previous Article in Journal
Wind Tunnel Measurements for Flutter of a Long-Afterbody Bridge Deck
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

BeiDou Signal Acquisition with Neumann–Hoffman Code Modulation in a Degraded Channel

Harbin Engineering University, College of Automation, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(2), 323; https://doi.org/10.3390/s17020323
Submission received: 7 December 2016 / Revised: 19 January 2017 / Accepted: 19 January 2017 / Published: 9 February 2017
(This article belongs to the Section Remote Sensors)

Abstract

:
With the modernization of global navigation satellite systems (GNSS), secondary codes, also known as the Neumann–Hoffman (NH) codes, are modulated on the satellite signal to obtain a better positioning performance. However, this leads to an attenuation of the acquisition sensitivity of classic integration algorithms because of the frequent bit transitions that refer to the NH codes. Taking weak BeiDou navigation satellite system (BDS) signals as objects, the present study analyzes the side effect of NH codes on acquisition in detail and derives a straightforward formula, which indicates that bit transitions decrease the frequency accuracy. To meet the requirement of carrier-tracking loop initialization, a frequency recalculation algorithm is proposed based on verified fast Fourier transform (FFT) to mitigate the effect, meanwhile, the starting point of NH codes is found. Then, a differential correction is utilized to improve the acquisition accuracy of code phase. Monte Carlo simulations and real BDS data tests demonstrate that the new structure is superior to the conventional algorithms both in detection probability and frequency accuracy in a degraded channel.

1. Introduction

The application of global navigation satellite systems (GNSS) in harsh environments is drawing increasing attention among civilian users [1]. As one of the world’s major navigation system, BDS is scheduled to be completed in 2020, and will eventually operate in ocean, urban, and indoor environments. To improve the capability for anti-narrowband interference and provide more precise ranging information, a secondary modulation of the NH code is introduced [2,3]. However, the periodicity of the transmitted sequence was broken by the NH codes, leading to more polarity changes, also known as bit transitions.
For typical receiver modules, the acquisition process occurs first and determines the sensitivity of a receiver [4]. The maximum sensitivity gain is achievable by the integration of consecutive correlations of the received signal and the local code replica [5]. The integration time becomes one of the main concerns that should be addressed in challenging environments affected by bit transitions, square loss, computational burden, and so on [6].
Coherent combination (CC) is the most classical signal acquisition algorithm which can obtain the maximum signal-to-noise ratio (SNR), but the integration time is limited by bit transitions [7]. In order to solve this problem, some improved CC algorithms are proposed. The half-bit alternation method that combines CC can be adopted to skip bit transitions [8], and a block acquisition scheme using the CC method can also alleviate this effect [9], meanwhile, a number of weighted signal blocks can be coherently accumulated to obtain an improved power [10]. However, these algorithms based on the foundation that the probability of bit transitions is 5% for the global positioning system (GPS) L1 signal. However, the probability of bit transitions is 50% for the BDS B1 signal caused by the NH codes, which can significantly influence the aforementioned integration algorithms [11,12].
Apart from the CC method and its modified version, two other well-known techniques, namely, the non-coherent combination (NC) algorithm and the differential combination (DC) algorithm, are also applied in weak signal environments [13,14]. The NC algorithm mitigates the undesirable effect of polarity changes through a squaring operation, but has to suffer from square loss, which reduces the SNR [15,16]. The DC method is notable for multiplying the present integration result with the delay conjugate of the last one, which reduces the square loss and the effect of polarity changes to a certain extent, thereby making this method feasible in weak GPS signal processing [17]. Different forms of detection schemes can be found in literature [17,18,19].
Regardless of the calculation burden, bit transition is a major problem that hinders long integration time, it can degrade the probability of acquisition in terms of Galileo E1 OS, GPS L5, and GPS L1C signal in different degrees [20], as well as the BDS B1 signal with a similar modulation. To achieve high frequency accuracy and acquisition probability in weak environments, an aiding Doppler information from strap-down inertial navigation system can be added into the acquisition process to reduce the search band of frequency. Meanwhile, the Viterbi algorithm is used for bit estimation to obtain a long differential integration time [6]. To demodulate the NH code without aiding information, a search tree for secondary codes is proposed combined with the CC method [21]; however, 220 combinations for 20 ms integration time will consume a very large amount of storage. Simply adopting a whole tiered code of pseudo-random noise (PRN) codes and NH codes as the local code for acquisition is an easier alternative [22], however, conducting that many FFTs and inverse FFTs of a 20 ms sequence in each frequency unit is difficult. On the basis of the same principle, a modified CC method, that uses the top 5 ms NH chips with the same polarity as the local code, is utilized to acquire BDS signal [23], but this technique requires 19 instances of redundant computation in each frequency unit.
In the above algorithms, FFT is used for transformation from the time domain to the frequency domain, or from the frequency domain to the time domain, to reduce the computation burden. Additionally, the basic function of FFT is to analyze the spectrum with a simplified search method, and the resolution of the spectrum is proportional to the length of the data [24]. Based on the theory and the analysis of paper [25], the present study derives a formula that shows exactly how bit transitions can bring in an error Doppler estimation, and related analysis is addressed. A frequency recalculation algorithm is proposed which exploits the efficiency of FFT-based algorithms, the start point of NH codes can be found, that is to say the symbol combination is found, and long integration time becomes possible. Meanwhile, a more accurate Doppler bins can be estimated. This result is helpful for speeding up the tracking process. The developed algorithm is analyzed in terms of cross-ambiguity function (CAF). Monte Carlo simulations are used to analysis the results, and the algorithms are demonstrated through real data.
The remainder of this paper is organized as follows: Section 2 introduces the BDS signal model and reviews the acquisition process. Section 3 analyzes the side effect of bit transitions on acquisition and derives a theoretical formula in the presence of bit transitions. In Section 4, a FFT-based acquisition algorithm is proposed to get the fine frequency. In Section 5 the different integration techniques are further analyzed and compared by means of Monte Carlo simulations and real data. Finally, conclusions are drawn in Section 6.

2. Signal and System Model

Unlike conventional GPS signals, BDS B1 signal adopts secondary codes modulation by the primary purpose of speeding up bit synchronization. Thus, the period of NH codes is 20 ms, corresponding to a navigation bit interval, besides, the secondary codes and the navigation bits are time aligned. Figure 1 shows the structure of BDS B1 signal. The received signal of the B1 in-phase branch in a one-path additive white Gaussian noise channel can be modeled as [26]:
y ( t ) = k = 1 K A k d k ( t τ 0 k ) c k ( t τ 0 k ) s k ( t τ 0 k ) cos ( 2 π ( f RF + f 0 k ) t + φ 0 k ) + η ( t )
where, superscript k stands for satellite number depending on the modulated PRN code; A is the amplitude; and K is the number of visible satellites. c, d, and s represent the PRN codes, the navigation data and the NH codes, respectively. f RF is the carrier frequency, while f 0 , τ 0 , and φ 0 are the Doppler frequency, the delay and the carrier phase introduced when transmission. η is a stationary additive white Gaussian noise with variance σ η 2 .
The received signal is filtered, down-converted, sampled at the frequency f s = 1 / T s , and then digitalized. Thus, the digital signal can be expressed as:
y [ n ] = k = 1 K A k d k ( n T s τ 0 k ) c k ( n T s τ 0 k ) s k ( n T s τ 0 k ) cos ( 2 π ( f IF + f 0 k ) n T s + φ 0 k ) + η ( n T s )
where f IF is the receiver intermediate frequency.
Owing to the quasi-orthogonality of the PRN codes, different signals can be analyzed separately. Thus, the following signal model can be adopted for simplification:
y [ n ] = A d ( n T s τ 0 ) c ( n T s τ 0 ) s ( n T s τ 0 ) cos ( 2 π F 0 n + φ 0 ) + η ( n T s )
where the index k has been dropped, and F 0 = ( f IF + f 0 ) T s is introduced. It should be noted that the input signal of conventional GPS L1 signal has a similar form, but there are no secondary codes, that is, there is no s ( n T s τ 0 ) term.
The basic concept of signal acquisition is to obtain available satellite PRN numbers, estimated Doppler shifts, and code phases through correlation between signal of Equation (3) and the local duplicates, then a random variable R ( F D , τ ) is produced, where, F D = ( f IF + f d ) T s and τ stands for different carrier frequencies and delays to be tested. When the GPS L1 signal is present [26]:
R ( F D , τ ) = A 2 d sin ( π N Δ F ) π N Δ F R ( Δ τ ) exp ( j Δ φ ) + η
R ( ) is the cross correlation between the input and the local code, N is number of samples, and T c = N T s defines the coherent integration time in seconds, T c is set to 0.001s. Δ F = F 0 F D is the frequency difference between the input and local signal, while Δ τ = ( τ 0 τ ) / T s represents the difference between the input and local code delay, normalized by sampling interval T s , Δ φ is the difference of phases between the input and local carriers, d is 1 or −1 that represents the navigation bits. Equation (4) is derived based on the assumption that the navigation bits will not reverse during the integration process. It is reasonable because there exists, at most, one bit transition in every 20 ms, though the one variable affected by bit transition is no longer in the same form with Equation (4), it hardly affect the final result. However, when it comes to the BDS B1 signal, the equation is not applicable, and the general Equation (4) must be revised attributing to the NH codes. The next section will derive a new variable that fits for the BDS B1 signal and a detailed analysis will be drawn.

3. Analysis of Bit Transitions

The CC acquisition is realized by averaging the in-phase (I) and quadrature (Q) components through correlation [27]. The derivation of variable R ( F D , τ ) can easily be found in the related literature, thus, part of the derivation in this paper is omitted. The same goes for the derivation of decision variable for GPS signals, the assumption that navigation bits will not reverse during the integration process is made. When the local code replica perfectly matches the input signal, that is τ = τ 0 , the CC results can be expressed as:
R ( F D , τ ) = n = 0 N 1 A d ( n T s τ 0 ) c ( n T s τ 0 ) s ( n T s τ 0 ) cos ( 2 π F 0 n + φ 0 ) c l o c ( n T s τ ) exp ( j 2 π F D n ) + η = A 2 n = 0 N 1 s ( n T s τ 0 ) { exp [ j ( 2 π Δ F n + φ ) ] + exp [ j ( 2 π ( F 0 + F D ) n + φ ) ] } + η
Suppose that bit transitions occur in N τ ( 0 N τ N 1 ), namely:
s ( n T s τ 0 ) = { 1 f o r 0 n < N τ 1 f o r N τ n N 1
Considering that high frequency components can be mitigated by integration, then,
R ( F D , τ ) A 2 exp ( j φ ) { n = 0 N τ 1 exp [ j ( 2 π Δ F n ) ] n = N τ N 1 exp [ j ( 2 π Δ F n ) ] } + η
The two terms in Equation (7) are two truncated geometrical series, which can be summed easily as [25]:
R ( F D , τ ) A 2 exp ( j φ ) n = 0 N 1 s ( n T s τ 0 ) exp ( j 2 π Δ F n ) = A 2 exp ( j φ ) { exp ( j α 1 ) sin ( π Δ F N τ ) sin ( π Δ F ) exp ( j α 2 ) sin [ π Δ F ( N N τ ) ] sin ( π Δ F ) }
where α 1 = π Δ F ( N τ 1 ) and α 2 = π Δ F ( N + N τ 1 ) . When N τ = N / 2 , R ( F D , τ ) can be rewritten as:
R ( F D , τ ) = A 2 exp ( j φ ) { exp [ j π Δ F ( N 2 1 ) ] exp [ j π Δ F ( 3 N 2 1 ) ] } sin ( 1 2 π Δ F N ) sin ( π Δ F ) = A 2 exp ( j φ ) exp [ j π Δ F ( N 2 1 ) ] [ 1 exp ( j π Δ F N ) ] Φ ( Δ F ) = A 2 exp ( j φ ) Γ ( Δ F ) Φ ( Δ F )
where, Φ ( Δ F ) = sin ( 1 2 π Δ F N ) / sin ( π Δ F ) , and Γ ( Δ F ) = exp [ j π Δ F ( N 2 1 ) ] [ 1 exp ( j π Δ F N ) ] .
The Doppler frequency is within 10 kHz under the existing vehicle and satellite conditions [28], whereas f s is greater than twice the intermediate frequency (IF). That is, Δ F , which denotes the ratio of the frequency error, and f s , is extremely small, i.e., smaller than 0.0025. Thus, sin ( π Δ F ) can be replaced with π Δ F as follows:
R ( F D , τ ) = A N 4 exp ( j φ ) sin ( π Δ F N / 2 ) π Δ F N / 2 Γ ( Δ F ) = A N 4 exp ( j φ ) sinc ( Δ F N / 2 ) Γ ( Δ F )
Then, the decision variable of CC is determined as:
S ( F D , τ ) = | R ( F D , τ ) | 2 = [ A N 4 sinc ( 1 2 Δ F N ) ] 2 | 1 exp ( j π Δ F N ) | 2 = 2 [ A N 4 sinc ( 1 2 Δ F N ) ] 2 [ 1 cos ( π Δ F N ) ]
The NC technique is realized by summing the squared coherent correlation values. This method can be written as [13]:
N ( F D , τ ) = n = 1 N NC | R n ( F D , τ ) | 2
where N NC is the time of non-coherent integration. However, a famous square loss will be experienced because the noise is also squared [15].
The DC method will be achieved via two consecutive coherent integration periods. The output can be presented as follows [19]:
Y ( F D , τ ) = | n = 1 N DC R n ( F D , τ ) R n 1 ( F D , τ ) |
where N DC is the time of differential integration; R n ( F D , τ ) and R n 1 ( F D , τ ) are the nth output and the complex conjugate of the (n−1)-th output, respectively.
In Equation (11), sinc 2 ( Δ F N / 2 ) and [ 1 cos ( π Δ F N ) ] will determine the position of the peak in the frequency bins. sinc 2 ( Δ F N / 2 ) will monotonically decrease to zero in | Δ F | [ 0 , 2 / N ] , in other words, the peak decreases as the frequency error increases in the range | Δ f | = | Δ F | f s [ 0 , 2000 ] Hz. However, when Δ F is equal to zero, [ 1 cos ( π Δ F N ) ] is equal to zero too. That is to say, the peak is no longer consistent with the right Doppler bins. When Δ F is equal to ± 1 / N , i.e., the frequency difference is equal to ± 1000 Hz, a maximum of [ 1 cos ( π Δ F N ) ] will be achieved and sinc 2 ( Δ F N / 2 ) [ 1 cos ( π Δ F N ) ] will reach its maximum, which will bring in a large frequency error. The curve of sinc 2 ( Δ F N / 2 ) [ 1 cos ( π Δ F N ) ] is illustrated in Figure 2.
Figure 3 presents the acquisition results with simulated BDS signal of satellite 12 when T c = 0.001 s, f s = 16.368 MHz, and f IF = 4.092 MHz, thus, N = 1 , 6368 , the code phase, and also N τ , is set as 0.5 N . The starting point of the NH codes are set to 11 NH code chips, and the Doppler bins are set to 3000 Hz. The bin width of the frequency dimension and code dimension are 500 Hz and 0.5 chips, respectively, for the full paper. As shown in Figure 3, when the integration process crosses bit transitions, the peak in the frequency axis splits into two components. The corresponding Doppler bins of the peak is 2500 Hz, which indicates a failed acquisition. Considering that the CC method serves as the base of the NC and DC methods, the outputs of CC are accumulated by means of square or delay conjugate multiplication to reduce the white Gaussian noise. However, such operations cannot eliminate the frequency error, and it will be transmitted to the final decision variable. Moreover, Equation (11) indicates that the code phase is still correct, except for a loss of amplitude. Consequently, a FFT-based frequency recalculation method is designed to mitigate the undesirable effect of bit transitions.

4. FFT Based Frequency Recalculation

Considering that the frequency of conventional acquisition algorithms deviates from the desired value because of the secondary codes, a more accurate frequency must be obtained. Generally, the frequency spectrum of a digital sine wave can be obtained via discrete Fourier transform, and the revolution of the spectrum is proportional to the ratio of the sampling frequency and the samples. To apply the theory to GNSS signal acquisition, the PRN codes, NH codes, and navigation bits must be demodulated from the received signal to obtain the spectrum of the carrier. The code phase estimated from conventional acquisition algorithms can be used to remove the PRN code. Thus, there remains the NH codes and navigation bits.
In this study, a frequency search structure is adopted to remove the NH codes and the navigation bits. First, we conduct a conventional CC acquisition to obtain the coarse code phase of PRN codes. Then, we make the local PRN codes align with the received signal at the first bit edge. Next, wew correlate the signal of step two with the whole tiered codes of the PRN codes and NH codes 20 times, of which the delay is 1 ms. Finally, we convert the correlations into frequency domain. There will be 20 outputs, and each output has a maximum value. The largest of these maximums means that both the NH codes and the navigation bits are removed, and a more precise frequency of the received signal is obtained. The diagram of the frequency search structure is shown in Figure 4.
In open areas, the length of the FFT can be shortened. More specifically, any consecutive chips of the PRN-NH modulated codes can be used for the proposed algorithm. However, the sensitivity differs significantly, as shown in the next section. To satisfy the requirement of the tracking loop initialization in weak signal environments, the whole tiered codes of the PRN and NH codes are a better choice.
Figure 5 and Figure 6 present the results of the new scheme of the same signal used in Figure 3. These parameters will not repeat here. In Figure 5, the length of the FFT is 524,288 samples, That is, the revolution of the spectrum is f s / 524288 = 31.2 Hz. Figure 5 shows that the index of the maximum value is 11, which is consistent with the default value. Meanwhile, in the carrier frequency dimension, the peak appears at 4,095,018 Hz, which represents a slight error of 18 Hz. Figure 6 shows the frequency error against different NH code phases in a more direct way. Only when the local NH codes is aligned with the received signal can the peak reach the maximum value. In other cases, the frequency error of tens to hundreds of hertz occurs.
When frequent bit transitions are considered in urban environments, the NC algorithm should be adopted before using the FFT-based frequency recalculation algorithm to increase the peak. In addition to the accuracy frequency, the start point of NH codes can also be found through the proposed scheme, thus, the NH codes can be removed, and there remains only the navigation bits. Then, a conventional DC algorithm will easily be conducted with the fine frequency to correct the code phase of the DC algorithm. The final code phase τ ¯ of the scheme is determined by the following criterion:
{ τ ¯ = τ NC f o r   r NC > r DC τ ¯ = τ NC Δ τ DC f o r   r NC r DC
where, r NC and r DC represent the peak to second noise ratio of the NC and DC algorithms, respectively. The peak to second noise ratio is usually used for threshold detection, here, it indicates the reliability of the conventional NC and DC algorithms as well. The outputs with high peak to second noise ratio are more reliable. τ NC and Δ τ DC denote the code phases obtained from the conventional NC and DC algorithms, respectively. More specifically, Δ τ DC stands for the code phase error of the NC algorithm owing to step 2 of the frequency recalculation algorithm. If the code phase estimation τ NC is correct, Δ τ DC equals 0; if not, τ NC is corrected with Δ τ DC .

5. Performance Assessment

To achieve a comprehensive assessment of the proposed acquisition scheme, Monte Carlo simulations were conducted to compare the proposed scheme with the conventional NC and DC algorithms. A schematic simulation platform was designed. Furthermore, BDS signals received by a GNSS receiver and other auxiliary equipment were also used to implement the expected scheme.

5.1. Monte Carlo Simulations

Monte Carlo simulations were used for comprehensive evaluation of different algorithms. To prove the necessity and reliability of the proposed algorithm, comparisons were done to study the sensitivity of different signal lengths, and the influence of the coarse code phase obtained from the conventional NC algorithm was analyzed. Then, the detection probability and frequency error were plotted as a function of the C/N0s and the bit transition positions for comparison with the conventional algorithms.
The platform is shown in Figure 7, and the signal channel could be achieved by adjusting the amplitude and delay of the IF signal. The following basic settings were used to obtain statistic results: A software simulator generates the BDS signal with an intermediate frequency of 4.092 MHz and a sampling frequency of 16.368 MHz. The Doppler frequency distributes uniformly in [−10000, +10000] Hz, and the code phase distributes uniformly in [0, 1] chips. The front-end filter bandwidth is 4.5 MHz. The coherent integration time was set to 1 ms, considering the secondary code. Furthermore, additional white Gaussian noise was generated for each trial, and 1000 trials were used for each probability and frequency value.
The first simulation tests the effect of signal lengths on the proposed algorithm. The total integration time was set to 100 ms, and different C/N0s were tested. The detection probability and average frequency error are shown in Figure 8. With the increase of the length of FFT, the detection probability increases as well, while the frequency error decreases. It is clear that, with the increase of C/N0, the performance improves greatly. When the C/N0 is greater than 31 dB-Hz and the data length is greater than 8 ms, the performance differs slightly only in the error of the frequency. However, different data lengths show a large difference in the frequency error. In high C/N0 conditions, the FFT of 10ms length or shorter is enough for the tracking process. Considering the conditions in the following analysis, 20 ms-length FFT is used for better performance if there is no special notation.
The second simulation analyzes the reliability of the proposed algorithm against the accuracy of the coarse code phase. The total integration time was set to 100 ms. The error of code phase obtained from the conventional NC algorithm was set to [−0.5, +0.5] PRN code chips, which are equivalent to [−4, +4] samples. The results are shown in Figure 9. The top panel shows the distribution of the code phase error of successful NC acquisition, and the code phase error have been rounded up and down to an integral sample. It is clear that the error of the code phase distributes mainly within one sample. The middle and bottom panels show the detection probability and frequency error of the proposed algorithm, respectively. When the error of code phase is less than or equal to 0.125 chips, there is no loss of detection probability. When the code phase error increases further, the detection probability decreases. While the frequency error of the proposed algorithm remains stable.
The third simulation verifies a comparable acquisition performance among the three algorithms against the position of bit transitions. To better observe the side effect of bit transitions, the total integration time was set to 10 ms, and a high C/N0 of 45 dB-Hz was adopted. The detection probability and average frequency error are shown in Figure 10 and Figure 11, respectively, which clearly indicate that the detection probability and frequency error of the conventional NC and DC algorithms are functions of bit transition positions. More specifically, there shows a valley in the curve of detection probability when bit transitions appear on 0.5 chips, and the maximums show up on both sides while the frequency error presents the opposite trend. This phenomenon is consistent with the analysis of Section 3 that bit transitions do affect the acquisition performance. As for the reason why the performance of the NC algorithm is better than that of the DC algorithm, the next paragraph will give a detailed analysis. In contrast, the detection probability of the proposed algorithm is hardly affected by bit transitions, and the frequency error is much smaller than the conventional algorithms.
The fourth simulation verifies the acquisition performance against C/N0 ranging from 25 dB-Hz to 35 dB-Hz, the total integration time was set to 100 ms. It is clear that the proposed algorithm demonstrates considerable advantage in Figure 12. When the C/N0 is 30.3 dB-Hz, the proposed algorithm can achieve a probability of 90%, which is 22% higher than that of the NC algorithm, and 6.5% more than the DC algorithm. As is known, when the GPS L1 signal is present, the bit transition loss of the DC algorithm can be negligible when long integration time is experienced, and the acquisition probability is higher than that of the NC algorithm. However, Figure 12 shows that the NC algorithm exhibits a slightly better performance than the DC algorithm although it has to endure squared losses, it is consistent with the analysis that frequent bit transitions of NH codes can deteriorate the acquisition performance, especially of algorithms that are affected by the polarity of the variable. Note that the DC algorithm is not affected directly by the signal polarity, but the product of two adjacent units. Considering the polarity of the secondary codes, bit transitions can decrease the performance of the DC algorithm, too. The analysis also applies to the results of the frequency in Figure 13. Moreover, it is shown that the detection probabilities of the NC and DC algorithms increase slowly when the C/N0 is higher than 32 dB-Hz, because the performance is mainly affected by frequent bit transitions in this case. Large frequency errors caused by bit transitions are the main reason, while C/N0 is only the secondary cause.

5.2. Real Data Tests

To confirm that the new approach exhibits better behavior, real data tests were conducted. The comprehensive test platform is shown in Figure 14. It should be noted that only the GNSS antennas, the intermediate frequency signal sampler, and the ProPak6 receiver were used, the rest of the devices in the figure were used for other needs, and there will be no description here. The intermediate frequency signal sampler can convert the high frequency BDS B1 signal into a 3.996875 MHz digital IF signal with a sampling frequency of 16.369 MHz. The ProPak6 receiver is the most advanced receiver made by NovAtel Company seated in Calgary, Alberta, Canada, and the supplier is Beijing BDStar Navigation Co., Ltd. The ProPak6 receiver can provide meter to centimeter level positioning precision, thus it was used to judge the frequency accuracy of the algorithms.
The parameters are as follows: the non-coherent integration time is 100 ms combined with 1 ms coherent integration, the data length of the proposed algorithm is 20 ms. Satellite 9 was visible and the Doppler bins obtained by the ProPak6 receiver is 739.355 Hz. The results of the conventional NC algorithm and the proposed scheme are shown in Figure 15 and Figure 16, respectively.
Figure 15 shows a slight peak split in the Doppler frequency axis caused by bit transitions, which results in a frequency error of 239.355 Hz compared with the reference receiver. The peak split is not serious because several integration results without bit transitions occur as a result of the long integration time. Further increasing the non-coherent integration time will hardly help mitigating the undesirable effect of bit transitions because the ratio of the results with, and without, bit transitions are determined by the structure of the NH codes. While Figure 16 mitigates this phenomenon by using the proposed algorithm, and shows that the start point of NH codes is eight, the frequency error is about 38.7 Hz compared with the reference Doppler, much smaller than the conventional one.

6. Conclusions

This study selects BDS signal as the research object. It addresses the acquisition of a weak signal channel, and derives the formula of CAF affected by bit transitions of conventional acquisition methods. The results indicate that the frequent bit transitions in BDS signals can cause a large frequency error up to 1000 Hz, which seriously affects the acquisition performance. To address this issue, a frequency recalculation algorithm is proposed. This algorithm uses an additional FFT operation based on the conventional NC algorithm and the coarse code phase, and can obtain the starting point of NH codes at the same time. Then, a fine DC acquisition process is performed to further correct the coarse code phase obtained from the conventional NC algorithm.
Simulations show that the new acquisition technique is reliable when the conventional NC is successful, and the new acquisition scheme exhibits excellent performance in terms of frequency accuracy and detection probability compared with conventional algorithms. More specifically, the proposed acquisition algorithm can acquire BDS signal of 30.3 dB-Hz effectively within 100 ms even if there exists frequent bit transitions, about 22% higher than conventional algorithms, meanwhile, the acquisition accuracy is higher than 20 Hz. When the delay of PRN codes is 0.5 chips in 45 dB-Hz, the detection probabilities of conventional acquisition algorithms is lower than 60%, while the new algorithm can reach 98% with an accuracy of about 25 Hz. In addition, the start point of NH codes can also be obtained, which is capable of long differential integration time. With the promotion of GNSS modernization, there will be a high demand of the signal acquisition algorithms.

Acknowledgments

This work is co-funded by the National Natural Science Foundation of China under grant nos. 61304234, 61304235, 61273081, the Fundamental Research Funds for the Central Universities under grant nos. HEUCFX041403, and the Post-Doctoral Scientific Research Foundation, Heilongjiang Province under grant nos. LBH-Q15033 and LBH-Q14054.

Author Contributions

Jicheng Ding and Aimeng Liu conceived and designed the experiments; Aimeng Liu performed the experiments; Lin Zhao contributed experiment tools and analyzed the data; Aimeng Liu, Jicheng Ding and Jing Wang wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Quan, Y.; Lau, L.; Roberts, G.W.; Meng, X. Measurement Signal Quality Assessment on All Available and New Signals of Multi-GNSS (GPS, GLONASS, Galileo, BDS, and QZSS) with Real Data. J. Navig. 2015, 69, 313–334. [Google Scholar] [CrossRef]
  2. Corazza, G.E.; Palestini, C.; Pedone, R.; Villanti, M. Galileo Primary Code Acquisition Based on Multi-hypothesis Secondary Code Ambiguity Elimination. In Proceedings of the 20th International Technical Meeting of the Satellite Division of The Institute of Navigation, Fort Worth, TX, USA, 25–28 September 2007; pp. 2459–2465.
  3. China Satellite Navigation Office, BeiDou Navigation Satellite System Signal in Space Interface Control Document Open Service Signal (Version 2.1). Available online: http://en.beidou.gov.cn/index.html (accessed on 18 November 2016).
  4. Borre, K.; Akos, D.M.; Bertelsen, N.; Rinder, P.; Jensen, S.H. A Software-Defined GPS and Galileo Receiver—A Single Frequency Approach; Birkhäuser Boston: Berlin, Germany, 2007; pp. 70–71. [Google Scholar]
  5. Esteves, P.; Sahmoudi, M.; Boucheret, M.L. Sensitivity Characterization of Differential Detectors for Acquisition of Weak GNSS Signals. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 20–37. [Google Scholar] [CrossRef]
  6. Wang, X.; Ji, X.; Feng, S. A Scheme for Weak GPS Signal Acquisition Aided by SINS Information. GPS Solut. 2014, 18, 243–252. [Google Scholar] [CrossRef]
  7. Yang, C.; Miller, M.; Blasch, E.; Nguyen, T. Comparative Study of Coherent, Non-Coherent, and Semi-Coherent Integration Schemes for GNSS Receivers. In Proceedings of the 20th International Technical Meeting of the Satellite Division of The Institute of Navigation, Cambridge, MA, USA, 25–28 September 2007; pp. 572–588.
  8. Ma, Y.; Zhang, Y.; Zhang, Z.; Ma, G. Modified Method of High Dynamic & High Sensitivity GPS Signal Acquisition Systems. J. Syst. Eng. Electron. 2009, 31, 265–269. [Google Scholar]
  9. Psiaki, M.L. Block Acquisition of Weak GPS Signals in A Software Receiver. In Proceedings of the 14th International Technical Meeting of the Satellite Division of the Institute of Navigation, Salt Lake City, UT, USA, 11–14 September 2001; pp. 1–13.
  10. Sahmoudi, M.; Amin, M.; Landry, R.J. Acquisition of Weak GNSS Signals Using a New Block Averaging Pre-Processing. In Proceedings of the IEEE/ION Position Location And Navigation Symposium, Monterey, CA, USA, 6–8 May 2008; pp. 1362–1372.
  11. Yan, K.; Zhang, H.; Zhang, T.; Xu, L.; Niu, X. Analysis and Verification to the Effects of NH Code for Beidou Signal Acquisition and Tracking. In Proceedings of the 26th International Technical Meeting of The Satellite Division of the Institute of Navigation, Nashville, TN, USA, 16–20 September 2013; pp. 107–113.
  12. Zou, D.; Deng, Z.; Huang, J.; Liu, H.; Yang, L. A study of Neuman Hoffman Codes for GNSS Application. In Proceedings of the 5th International Conference on Wireless Communication, Beijing, China, 24–26 September 2009; pp. 1–4.
  13. Choi, H.; Park, S.H.; Cho, D.J.; Yun, S.J.; Kim, Y.B.; Lee, S.J. A Novel Weak Signal Acquisition Scheme for Assisted GPS. In Proceedings of the 15th International Technical Meeting of the Satellite Division of The Institute of Navigation, Portland, OR, USA, 24–27 September 2002; pp. 177–183.
  14. Borio, D.; O’Driscoll, C.; Lachapelle, G. Coherent, Noncoherent, and Differentially Coherent Combining Techniques for Acquisition of New Composite GNSS Signals. IEEE Trans. Aerosp. Electron. Syst. 2009, 45, 1227–1240. [Google Scholar] [CrossRef]
  15. Strässle, C.; Megnet, D.; Mathis, H.; Bürgi, C. The Squaring-loss Paradox. In Proceedings of the 20th International Technical Meeting of the Satellite Division of The Institute of Navigation, Fort Worth, TX, USA, 25–28 September 2007; pp. 2715–2722.
  16. Borio, D.; Akos, D. Noncoherent Integrations for GNSS Detection: Analysis and Comparisons. IEEE Trans. Aerosp. Electron. Syst. 2009, 45, 360–375. [Google Scholar] [CrossRef]
  17. Villanti, M.; Salmi, P.; Corazza, G.E. Differential Post Detection Integration Techniques for Robust Code Acquisition. IEEE Trans. Commun. 2007, 55, 2172–2183. [Google Scholar] [CrossRef]
  18. Lank, G.W.; Reed, I.S.; Pollon, G.E. A Semicoherent Detection and Doppler Estimation Statistic. IEEE Trans. Aerosp. Electron. Syst. 1973, 9, 151–165. [Google Scholar] [CrossRef]
  19. Zarrabizadeh, M.H.; Sousa, E.S. A Differentially Coherent PN Code Acquisition Receiver for CDMA Systems. IEEE Trans. Commun. 1997, 45, 1456–1465. [Google Scholar] [CrossRef]
  20. Foucras, M.; Ekambi, B.; Bacard, F.; Julien, O.; Macabiau, C. Optimal GNSS Acquisition Parameters when Considering Bit Transitions. In Proceedings of the IEEE/ION Position Location And Navigation Symposium, Monterey, CA, USA, 5–8 May 2014; pp. 804–817.
  21. Shivaramaiah, N.C.; Dempster, A.G.; Rizos, C. Exploiting the Secondary Codes to Improve Signal Acquisition Performance in Galileo Receivers. In Proceedings of the 21st International Technical Meeting of the Satellite Division of The Institute of Navigation, Savannah, GA, USA, 16–19 September 2008; pp. 1497–1506.
  22. Shi, M.; Peng, A.; Guang, O. Analysis to the Effects of NH Code for Beidou MEO/IGSO Satellite Signal Acquisition. In Proceedings of the 2014 IEEE 9th Conference on Industrial Electronics and Applications (ICIEA), Hangzhou, China, 9–11 June 2014; pp. 2075–2080.
  23. Bhuiyan, M.Z.H.; Söderholm, S.; Thombre, S.; Ruotsalainen, L.; Kuusniemi, H. Overcoming the Challenges of BeiDou Receiver Implementation. Sensors 2014, 14, 22082–22098. [Google Scholar] [CrossRef] [PubMed]
  24. Ma, C.; Zhou, Q.; Qin, J.; Xie, W.L.; Dong, Y.; Hu, W.S. Fast Spectrum Analysis for an OFDR Using the FFT and SCZT Combination Approach. IEEE Photonics Technol. Lett. 2016, 28, 657–660. [Google Scholar] [CrossRef]
  25. Sun, K.; Lo Presti, L. Bit Sign Transition Cancellation Method for GNSS Signal Acquisition. J. Navig. 2012, 65, 73–97. [Google Scholar] [CrossRef]
  26. Borio, D. M-Sequence and Secondary Code Constraints for GNSS Signal Acquisition. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 928–945. [Google Scholar] [CrossRef]
  27. El-Natour, H.A.; Escher, A.C.; Macabiau, C.; Boucheret, M.L. Impact of Multipath and Cross-correlation on GPS Acquisition in Indoor Environments. In Proceedings of the 2005 National Technical Meeting of The Institute of Navigation, San Diego, CA, USA, 24–26 January 2005; pp. 1062–1070.
  28. Tsui, J.B.Y. Fundamentals of Global Positioning System Receivers: A Software Approach; John Wiley & Sons Incorporate: New York, NY, USA, 2000; pp. 36–43. [Google Scholar]
Figure 1. BDS signal adopting NH code modulation
Figure 1. BDS signal adopting NH code modulation
Sensors 17 00323 g001
Figure 2. The curve of sinc 2 ( Δ F N / 2 ) [ 1 cos ( π Δ F N ) ] .
Figure 2. The curve of sinc 2 ( Δ F N / 2 ) [ 1 cos ( π Δ F N ) ] .
Sensors 17 00323 g002
Figure 3. Acquisition results with simulated signal of satellite 12.
Figure 3. Acquisition results with simulated signal of satellite 12.
Sensors 17 00323 g003
Figure 4. Diagram of the frequency recalculation algorithm.
Figure 4. Diagram of the frequency recalculation algorithm.
Sensors 17 00323 g004
Figure 5. Frequency recalculation result against the NH code phase and carrier frequency.
Figure 5. Frequency recalculation result against the NH code phase and carrier frequency.
Sensors 17 00323 g005
Figure 6. Frequency errors against different NH code phases.
Figure 6. Frequency errors against different NH code phases.
Sensors 17 00323 g006
Figure 7. The diagram of the simulation platform.
Figure 7. The diagram of the simulation platform.
Sensors 17 00323 g007
Figure 8. Influence of the FFT length on the proposed algorithm.
Figure 8. Influence of the FFT length on the proposed algorithm.
Sensors 17 00323 g008
Figure 9. Performance of the proposed algorithm against the code phase error.
Figure 9. Performance of the proposed algorithm against the code phase error.
Sensors 17 00323 g009
Figure 10. Detection probabilities of the NC, DC, and the proposed algorithms against the code phase.
Figure 10. Detection probabilities of the NC, DC, and the proposed algorithms against the code phase.
Sensors 17 00323 g010
Figure 11. Average frequency errors of the NC, DC, and the proposed algorithms against the code phase.
Figure 11. Average frequency errors of the NC, DC, and the proposed algorithms against the code phase.
Sensors 17 00323 g011
Figure 12. Detection probabilities of the NC, DC, and the proposed algorithms against C/N0.
Figure 12. Detection probabilities of the NC, DC, and the proposed algorithms against C/N0.
Sensors 17 00323 g012
Figure 13. Average frequency errors of the NC, DC, and the proposed algorithms against C/N0.
Figure 13. Average frequency errors of the NC, DC, and the proposed algorithms against C/N0.
Sensors 17 00323 g013
Figure 14. GNSS comprehensive test platform, (a) stand for GNSS antennas, (b) is the intermediate frequency signal sampler, and (c) is ProPak6 receiver.
Figure 14. GNSS comprehensive test platform, (a) stand for GNSS antennas, (b) is the intermediate frequency signal sampler, and (c) is ProPak6 receiver.
Sensors 17 00323 g014
Figure 15. Real data test results of the conventional NC algorithm.
Figure 15. Real data test results of the conventional NC algorithm.
Sensors 17 00323 g015
Figure 16. Real data test results of the proposed algorithm.
Figure 16. Real data test results of the proposed algorithm.
Sensors 17 00323 g016

Share and Cite

MDPI and ACS Style

Zhao, L.; Liu, A.; Ding, J.; Wang, J. BeiDou Signal Acquisition with Neumann–Hoffman Code Modulation in a Degraded Channel. Sensors 2017, 17, 323. https://doi.org/10.3390/s17020323

AMA Style

Zhao L, Liu A, Ding J, Wang J. BeiDou Signal Acquisition with Neumann–Hoffman Code Modulation in a Degraded Channel. Sensors. 2017; 17(2):323. https://doi.org/10.3390/s17020323

Chicago/Turabian Style

Zhao, Lin, Aimeng Liu, Jicheng Ding, and Jing Wang. 2017. "BeiDou Signal Acquisition with Neumann–Hoffman Code Modulation in a Degraded Channel" Sensors 17, no. 2: 323. https://doi.org/10.3390/s17020323

APA Style

Zhao, L., Liu, A., Ding, J., & Wang, J. (2017). BeiDou Signal Acquisition with Neumann–Hoffman Code Modulation in a Degraded Channel. Sensors, 17(2), 323. https://doi.org/10.3390/s17020323

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