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

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.


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].

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]: 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: where f IF is the receiver intermediate frequency.
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. RF f is the carrier frequency, while 0 f ,  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  ss 1/ fT , and then digitalized. Thus, the digital signal can be expressed as :  (2) where IF f is the receiver intermediate frequency.

Figure 1. BDS signal adopting NH code modulation
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 :  () R is the cross correlation between the input and the local code, N is number of samples, and  cs T NT defines the coherent integration time in seconds, c T is set to 0.001s.    0D F F F is the 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: 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(nT 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(·) is the cross correlation between the input and the local code, N is number of samples, and T c = NT 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.

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: Considering that high frequency components can be mitigated by integration, then, The two terms in Equation (7) are two truncated geometrical series, which can be summed easily as [25]: where α 1 = π∆F(N τ − 1) and α 2 = π∆F(N + N τ − 1). When N τ = N/2, R(F D , τ) can be rewritten as: where, Φ(∆F) = sin 1 2 π∆FN /sin(π∆F) , and Γ(∆F) = exp jπ∆F N 2 − 1 · [1 − exp(jπ∆FN)]. 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: Then, the decision variable of CC is determined as: The NC technique is realized by summing the squared coherent correlation values. This method can be written as [13]: (12) 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]: 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.
The curve of FN is illustrated in Figure 2.  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 los s of amplitude. Consequently, a FFT-based frequency recalculation method is designed to mitigate the undesirable effect of bit transitions.   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.
FN is illustrated in Figure 2.  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 los s of amplitude. Consequently, a FFT-based frequency recalculation method is designed to mitigate the undesirable effect of bit transitions.

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. Figures 5 and 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.

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 larges t 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  / 524288 31.2 s f 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: where, NC r and DC r 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 wi th high peak to second noise ratio are more reliable. NC   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: where, NC r and DC r 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 wi th high peak to second noise ratio are more reliable. NC  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: 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 .

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.

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/N 0 s 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

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.

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 chan nel could be achieved by adjusting the amplitude and delay of the IF signal. The following basic settings were used to obtain statistic results:  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 accur acy 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 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/N 0 s 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/N 0, the performance improves greatly. When the C/N 0 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/N 0 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 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 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/N 0 of 45 dB-Hz was adopted. The detection probability and average frequency error are shown in Figures 10 and 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 ex perienced, 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.  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 ex perienced, 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. The fourth simulation verifies the acquisition performance against C/N 0 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/N 0 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/N 0 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/N 0 is only the secondary cause.

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 Figures 15 and 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

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 Figures 15 and 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

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 Figures 15 and 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. 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.    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.    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.    Figure 16. Real data test results of the proposed algorithm.

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.