A New Acquisition Algorithm with Elimination Side Peak for All BOC Signals

A new inhibition side peak acquisition (ISPA) algorithm is proposed for binary offset carrier (BOC) modulated signals, which will be utilized in global navigation satellite systems (GNSS). We eliminate all side peaks of the BOC correlation function (CF) by structuring special sequences composed of PRN code and cycle rectangular sequences. The new algorithm can be applied to both generic sineand cosine-phased BOC signals, as well as to all modulation orders. Theoretical and simulation results demonstrate that the new algorithm can completely eliminate the ambiguity threat in the acquisition process, and it can adapt to lower SNR. In addition, this algorithm is better than the traditional algorithms in acquisition performance and inhibition side peak ability.


Introduction
With the development and application of global navigation satellite systems (GNSS) [1], GNSS signal receiving methods have become highly valued.Because the acquisition technology is the core of receiving; therefore, it also becomes a focus problem.Thus, mass acquisition algorithms [2,3] are proposed for GNSS signals to improve receiving performance.However, modern GNSS has provided new signals with longer PRN (pseudo random noise) codes and newer modulation methods, which aim to improve the positioning performance.Binary offset carrier (BOC) [4] modulated signals are the most widely used signal families in GNSS, and their side peak characteristics also require the highest technique complexity from GNSS receivers.
BOC modulation signal acquisition techniques focus on recovering the main correlation peak or eliminating ambiguities in the form of side peaks.At present, various techniques are proposed for side peak cancellation and are built on the basis of the correlation function (CF) of the BOC signals.Thus, the side band processing method originated from BPSK-like method [5,6], and then some improved methods [7][8][9] are proposed.The partial band is obtained by filtering or frequency domain processing in these kind methods, and then the main peak was estimated using similar BPSK characteristic.These kind methods can reduce the influence of subcarrier, but the energy and the necessary information are lost.Thus, the auxiliary signal methods [10,11] are mainly through the local auxiliary signal establishment to reach the purpose of removing side peaks.These kind methods can remove the side peak, but they lack universality.Thus, some effective methods [12][13][14] are proposed to improve the processing performance.However, these techniques apply only to sine-phased BOC signals.Thus, in [15], a mitigating ambiguity acquisition method is proposed.This technique can counterbalance the undesired side peaks, but it applies only to cosine-phased BOC signals.
In this paper, considering filter restriction and generic deficiency problems in traditional algorithms, we propose an inhibition side peak acquisition algorithm, which is applicable to all orders and to both generic sine-and cosine-phased BOC signals.where   is the modulated PRN code,    () is the subcarrier, 2  is the subcarrier cycle,    () is the spread spectrum symbol,  is the modulation order, and  and  0 , respectively, express the phase and time offset.
The BOC signal is usually expressed as BOC(  ,   ), the frequency of the subcarrier is   times the benchmark frequency, and the frequency of the PRN code is   times the benchmark frequency.The benchmark frequency is 1.023 MHz.The autocorrelation function of the BOC signal has multiple peaks and passes through zero many times.Its autocorrelation function consists of the positive peaks and the negative peaks, and the number of peaks is 2 − 1.The distance between peaks is   , and each peak height is (−1)  (− ||)/, where  is the serial number of the peaks.

The Acquisition Analysis.
From the perspective of algorithm generality, the acquisition algorithm for BOC modulation signal is usually divided into three categories, namely, the full band acquisition (FBA) algorithm [16], the peak optimization acquisition (POA) algorithm [17], and the single peak recovery acquisition (SPRA) algorithm [18].Their principles are shown in Figures 1, 2, and 3, respectively.FBA is a class of traditional algorithms, in which the correlation arithmetic is executed between the received signal and the original PRN code modulated by a square wave.POA is a class of improved algorithms, in which multiple correlations are executed to improve the main peak.SPRA is a class of new methods, in which a partial signal is separated from the received signal by the corresponding operations to inhibit the square wave.

ISPA Algorithm Structure
Let   be the sampling frequency of the BOC signal, and the frequency of the subcarrier and PRN code are   times and   times the benchmark frequency, respectively.Considering square wave modulation characteristics, the product model of the spread spectrum sequence and a series of rectangular sequences is structured, which can be approximately expressed as the BOC base-band signal model.Hence, the base-band signal may be represented by the following equation: where () is the message, () is the PRN code,  0 () is the mixed noise function caused by the discarded samples,  Δ () is the frequency error function cause by the front processing,  is the sequence position,  is the number of chips in accumulation time, and  is both the modulation order and

One chip
One chip The structured process of the rectangular sequence model.
the number of rectangular sequences in one chip, which is expressed as (3). is both the number of sampling points and the rectangular sequence width, which is expressed as (4).  ( +  − ) is a shifting rectangular sequence, which is expressed in (5), where () is the step sequence: Considering the represented model of the BOC baseband signal, the local rectangular sequence model is structured to inhibit the acquisition of side peaks.The structured process is shown in Figure 4, in which  is an odd number.The square wave sequence is shown in Figure 4(a), and the two structured cycle rectangular sequences are shown in Figures 4(b) and 4(c).The cycle rectangular sequences can also be structured for an even number  using the same principle.Further, the th cycle of two local channel rectangular sequences can be expressed as The original PRN code is, respectively, multiplied by the two-channel cycle rectangular sequences to structure the two new local channel sequences, which are expressed as where ( + ) is the delay PRN code and  is the time delay.
The beforehand processing received signal is executed by the correlation circumferential arithmetic with  1 () and  2 (), respectively, which are expressed as where ∧   () is the trigonometry sequence of width   and is expressed as When  is five, the autocorrelation result of the BOC base-band signal is shown in Figure 5(a), and the two structured correlation results are shown in Figures 5(b) and 5(c), respectively.The results show that the positions of the two channel main peaks exactly coincide with the position of the autocorrelation main peak, and the numbers of peaks are the same in both channels.In addition, the positions of the two channel peaks are symmetrical about the main peak position of the autocorrelation function.In view of these characteristics and combining ( 8) and ( 9), the addition and subtraction operations are performed by using the two structured correlation results, expressed as Thus, the new correlation function is structured to eliminate side peaks, and the processing is expressed as When the impacts of the frequency error and noise function are likely to be relatively weak, the relationship of the main peak value  1 in the Δ() function and the BOC autocorrelation function value  2 is expressed as To improve the peak, the result of Δ() is multiplied by a coefficient of /2 to obtain the final expression as

Performance Analysis
The Δ 1 () and Δ 2 () may be approximately represented by Hence, Δ 1 () satisfies a Gaussian distribution whose mean is (/4) and whose variance is  2 /2, and Δ 2 () satisfies Gaussian distribution whose mean is 0 and whose variance is Where  is the signal amplitude and  2 is the noise variance, the probability density function |Δ 1 ()| is expressed as and the probability density function |Δ 2 ()| is expressed as Thus, the Δ() probability density function is expressed as The false alarm probability of the ISPA algorithm is expressed as The acquisition detection probability of the ISPA algorithm is expressed as where  is the acquisition threshold.

Analysis and Simulation
5.1.Side Peak Inhibition Analysis.Equations ( 9) and (13) show that the final correlation result has a single peak whose main waveform is a triangular peak.Thus, the ISPA algorithm can achieve the goal of side peak inhibition.The new algorithm is then simulated using the following parameters: 10.23 MHz PRN code frequency, 15.345 MHz square wave frequency, and 122.76 MHz sampling frequency, modulation order of 3, and the sine-phased BOC signal for these parameters is expressed as sinBOC (15,10).
The ISPA result for sinBOC(15, 10) is shown in Figure 6.The ISPA results for sinBOC (10,2), cosBOC (10,5), and cosBOC(6, 1) are shown in Figures 7, 8, and 9, respectively.The simulation results show that the ISPA algorithm can clearly recover the main peak, whose position is the same as the main peak position of the autocorrelation function.
In particular, the ISPA algorithm can effectively inhibit side peaks.

Adaptability Analysis.
The new algorithm result is influenced by the frequency error and the mixed noise, according to (13).The algorithm result approximately conforms to the cycle equation because of the frequency error function cycle characteristics.When the relationship of the frequency error and accumulation time  satisfies (21), the algorithm result error reaches its maximum.We also find that the ISPA algorithm result decreases gradually along with the increase of mixed noise, according to where  is the positive integer.Furthermore, the ISPA algorithm's adaptability is simulated with the following parameters: 15.345 MHz square wave frequency and 122.76 MHz sampling frequency, modulation mode is sine mode, and modulation orders are 15, 10, 6, and 3, respectively.
The relationship between the relative main peak and frequency error is shown in Figure 10, which shows that the results satisfy (21).The relationship between the relative main peak and SNR is shown in Figure 11, revealing that the relative main peak decreases gradually with decreasing SNR.And the ISPA algorithm's adaptability to the SNR environment is more than −25 dB, according to (13) and Figure 11.

Superiority Analysis.
To verify the superiority of the ISPA algorithm, this ISPA algorithm is compared with other algorithms, namely, the FBA algorithm, POA algorithm, and SPRA algorithm.The simulation parameters are as follows: 2.046 MHz PRN code frequency, and the modulation mode is sine mode.
With changing modulation order, the main peak width changes and the main peak relative changes are shown in Figures 12 and 13.The results show that the ISPA algorithm's main peak width is the smallest, and its main peak relative result is the greatest, demonstrating that this algorithm's acquisition and tracking performance is the best.
The side peak relative changes and the main/side peak ratio changes with changing modulation order are shown in Figures 14 and 15.The results show that the ISPA algorithm side peak relative result is the smallest, and the main/side  peak ratio is the greatest, demonstrating that this algorithm's side peak inhibition ability is best.
The main peak relative changes with changing SNR are shown in Figure 16.The results show that the adaptability of the ISPA algorithm is better than the FBA algorithm and POA algorithm, but there are no significant differences between the ISPA algorithm and the SPRA algorithm.

Conclusions
In this paper, the principle and characteristics of BOC modulation signals have been studied.To implement the BOC modulated signal acquisition, effective algorithms have been studied, including the full band acquisition (FBA) algorithm, the peak optimization acquisition (POA) algorithm, and the single peak recovery acquisition (SPRA) algorithm.Considering the filter restriction and generic deficiency problems in traditional algorithms, we propose the ISPA algorithm.We eliminate all side peaks of the BOC correlation function (CF) by structuring special sequences composed of PRN code and cycle rectangular sequences.The ISPA algorithm can be applied to both generic sine-and cosine-phased BOC signals and to all modulation orders.In addition, it outperforms the traditional algorithms in acquisition, inhibition side peak ability, and adaptability to lower SNR.

Figure 3 :
Figure 3: The single peak recovery acquisition algorithm principle.

Figure 5 :
Figure 5: The correlation result of BOC.

Figure 10 :
Figure 10: The relationship between the relative main peak and frequency error.

Figure 11 :
Figure 11:  The relationship between the relative main peak and SNR.

Figure 12 :
Figure 12:  The relationship between the main peak width and the modulation order.

Figure 13 :Figure 14 :
Figure 13:  The relationship between the relative main peak and the modulation order.

Figure 15 :Figure 16 :
Figure 15:  The relationship between the main/side peak ratio and the modulation order.