Unambiguous Acquisition and Tracking ... - Radioengineering

1 downloads 0 Views 335KB Size Report
BOC if a classical delay lock loop is used in signal tracking. Therefore, an ...... thermal noise level, and the code tracking error standard de- viation caused by ...
FENG SHEN, GUANGHUI XU, JOON WAYN CHEONG, ET AL., UNAMBIGUOUS ACQUISITION AND TRACKING . . .

840

Unambiguous Unambiguous Acquisition Acquisition and and Tracking Tracking Technique Technique for for General BOC BOC Signals signals General Feng SHEN 1 , Guanghui XU 1 , Joon Wayn CHEONG 2 , Haiyu FENG 1 1 College

2 School

of Automation, Harbin Engineering University, Street Nantong, Harbin, China of Electrical Engineering and telecommunications, University of New South Wales, Sydney, NSW, 2052, Australia {sf407, xghdh407, fhydh407}@126.com, [email protected]

Abstract. This article presents a new unambiguous acquisition and tracking technique for general Binary Offset Carrier (BOC) ranging signals, which will be used in modern GPS, European Galileo system and Chinese BeiDou system. The test criterion employed in this technique is based on a synthesized correlation function which completely removes positive side peaks while keeping the sharp main peak. Simulation results indicate that the proposed technique completely removes the ambiguity threat in the acquisition process while maintaining relatively higher acquisition performance for low order BOC signals. The potential false lock points in the tracking phase for any order BOC signals are avoided by using the proposed method. Impacts of thermal noise and multipath on the proposed technique are investigated; the simulation results show that the new method allows the removal of false lock points with slightly degraded tracking performance. In addition, this method is convenient to implement via logic circuits.

Keywords Modern GPS, Galileo, BeiDou, BOC, unambiguous acquisition and tracking, multipath, thermal noise.

1. Introduction In recent years, applications of Global Navigation Satellite System (GNSS) are developing rapidly and attracting a great interest. However, GNSSs have to coexist in a limited number of frequency bands [1]. The traditional GPS L1 signal uses a spread spectrum technique with a Binary-Phase Shift Keying (BPSK) modulation. With the upcoming GNSS generation some civil signals will implement BOC modulation that allows improved code delay tracking while offering a spectral separation from BPSK signals due to its split spectrum [2]. However, the main drawback of every BOC modulated signal is its multi-peak AutoCorrelation Function (ACF) that implies possible false acquisition or biased tracking which will result in biased code delay measurements of up to hundreds of meters if no special care is taken.

DOI: 10.13164/re.2015.0840

Several techniques have been investigated to remove this biased threat [3]-[16]. Sideband filtering techniques in [3]-[5] can mitigate the ambiguity threat completely, but the filtering and the dual sideband processing increases the implementation complexity and causes correlation loss due to the mismatch between the received code (BOC) and the local code (shifted BPSK). Besides, these methods produce a wider ACF peak than the main ACF peak of a sine-phased BOC if a classical delay lock loop is used in signal tracking. Therefore, an additional transition from acquisition to tracking is necessary to avoid falsely locking onto an ACF sidepeak. [6], [7] are autocorrelation function based acquisition and tracking techniques. Bump-Jumping (BJ) technique aims at determining whether or not the peak being tracked is the correct one [7]. It exhibits high tracking accuracy when it is locked on the main peak, but the acquisition may have a high probability of false alarm when the signal-to-noise ratio is low. In addition, it needs time to recover from false lock, which makes it inapplicable in some time critical applications. Sub-carrier processing method such as Double Estimator Technique in [8] can completely solve the ambiguous problem while maintaining the tracking accuracy, but this method needs additional tracking loop for Sub-carrier. Sidepeak cancellation techniques [9]-[16] are based on the combination of constitutive sub-correlations of the BOC autocorrelation function. A General Removing Ambiguity via Side peak Suppression (GRASS) technique is proposed in [9] for sine-BOC(kn, n) signal and in [10] for cosine-BOC(kn, n) signal, which can suppress all of the undesired side peaks by combining the cross-correlation function between the received signal and the local auxiliary signal with the ACF of the received signal. However, the method in [10] presents high implementation complexity because of its residual two positive side peaks, which need to be calculated previously under different frond-end bandwidth. A Pseudo-CorrelationFunction (PCF)-based method presented in [12] does not have any side peaks and can completely remove all of the false lock points on the discriminator output. However, it is only work with sine phased BOC signals. A method referred to Autocorrelation Side-Peak Cancellation Technique (ASPeCT) was proposed in [16] which can realize both unambiguous acquisition and tracking. However, it is only dedicated to sine-BOC(n, n) signal. As a summation of the methods mentioned above, most of the techniques can not be

SIGNALS

RADIOENGINEERING, VOL. 24, NO. 3, SEPTEMBER 2015

841

used for both BOC signals acquisition and tracking process simultaneously, and most of them are also only suitable to a dedicated kind of BOC modulated signal. In this paper, an ambiguity removal technique suitable for generic BOC(x, y) signals is developed. Different from aforementioned side-peak cancellation techniques, the proposed method can realize both unambiguous acquisition and tracking process for generic BOC(x, y) signals. By employing two symmetrical local auxiliary signals, and combining the first correlation function between the received signal and one local auxiliary signal with the second correlation function between the received signal and another local auxiliary signal, the proposed technique constructs a test criterion which removes all of the undesired positive side peaks. Meanwhile, the false lock points can also be avoided by using modified Early-Minus- Late Power discriminator. The remainder of this paper is organized as follows. Firstly, we discuss the generic solution of the local modulated symbols for both sine-phased and cosine-phased BOC modulations. Then, the acquisition performance including detection probability of the proposed method is analyzed and evaluated. Meanwhile, the tracking accuracy including multipath and thermal noise resistance performance are discussed in detail. Followed by the implementation complexity analysis and structure simplification. Finally, the contributions of this paper and future steps are summarized.

2. BOC Modulated Signals and Method Description 2.1 BOC Signal Model Definition In this paper, a BOC modulation is denoted as BOC(x, y), where x means the ratio of the square wave frequency fs to 1.023 MHz, and y denotes the ratio of the spreading code rate fc to 1.023 MHz. The ratio M = 2x/y is referred as BOC-modulation order and it is constrained to a positive integer value. Since M is even in most of BOC signals employed in GNSS, for simplicity, we will discuss only the even case. The baseband direct sequence spread spectrum signal can be expressed as [12]

 ψk (t) =

1 t ∈ [kTc /M, (k + 1)Tc /M], 0 others.

(3)

For cosine-BOC signal, the modulated symbol is divided into 2M segments with equal length Ts = Tc /2M, then 2M−1

p(t) =



dk ψk (t)

(4)

k=0

where dk = (−1)k/2 , k = 0, 1, ..., 2M − 1, · represents the ceiling operation; and  1 t ∈ [kTs , (k + 1)Ts ], ψk (t) = (5) 0 others. Specially, for Multiplexed Binary Offset Carrier (MBOC) modulated signal M = 12 which has been widely (−1)k γ2 [13]; used in modern GNSS, dk = (−1) k/6 γ1 +  where · represents the floor operator, γ1 = 10/11 and γ2 = 1/11 are the amplitude weighting factors. The basic shape ψk (t) is the same as (3).

2.2 Proposed Method Description In [12], Yao et al. presented the Cross-Correlation Function (CCF) between the spreading signal and the Stepshape Code Symbol signal (SCS) which have the same fc and M. In fact, the step-shape modulated symbol p(t) can be identified by vector dscs = [d0 , d1 , ..., dM−1 ]T and spreading sequence rate fc = 1/Tc . So a step-shape modulated symbol can be denoted by p(t, dscs , fc ). Some examples of stepshape modulated symbols are shown in Fig. 1. In this paper, the CCF will be further generalized to cosine-BOC signals.                             Fig. 1. Step-shape modulated symbols.

(2)

Under the assumption that {ci } have ideal autocorrelation characteristics, the CCF between cosine-BOC and SCS signals using their waveforms pBOCc (t; dBOCC ; fc ) and pSCS (t; dSCS ; fc )(dSCS = [d0 , d1 , ..., d2M−1 ]T ) can be expressed as ⎧ (τ/Ts − k)(rk+1 − rk ) + rk ⎪ ⎪ ⎪ ⎪ τ ∈ [kTs , (k + 1)Ts ], ⎨ (τ/Ts − k + 2M)(rk−2M+1 − rk−2M ) + rk−2M RB/L (τ) = ⎪ ⎪ τ ∈ [(k − 2M)Ts , (k − 2M + 1)Ts ], ⎪ ⎪ ⎩ 0 others (6)

The shape vector can be described as dk = (−1)k , k = 0, 1, ..., M − 1; and the basic shape ψk (t) can be expressed as

where k = 0, 1, ..., 2M − 1, and rk is the Aperiodic CrossCorrelation Function (ACCF) between dBOCC and dSCS , i.e.,

S(t) = ∑(−1)ci p(t − iTc )

(1)

i

where ci is the spreading sequence of the binary digits {0, 1}, Tc is the symbol duration, and p(t) is the modulated symbol which is further divided into M segments with equal length Ts = Tc /M for sine-BOC signal. Then the modulated symbol can be represented as M−1

p(t) =

∑ dk ψk (t).

k=0

FENG SHEN, GUANGHUI XU, JOON WAYN CHEONG, ET AL., UNAMBIGUOUS ACQUISITION AND TRACKING . . .

842

⎧ ⎪ ⎪ ⎪ ⎪ ⎨ rk (dSCS ) =

⎪ ⎪ ⎪ ⎪ ⎩

⎧ (1) (2) ⎪ ⎨ r0 = r0 , (1) (2) ri · ri ≥ 0, i = 0, ⎪ ⎩ (1) (2) (1) (2) |ri + ri | = ( 1+k 1−k )|ri − ri |,

(−1)i/2 di+k , ∑2M−1−k i=0 0 ≤ k ≤ 2M − 1, 2M−1+k 1 (i−k)/2 d , (−1) ∑ i 2M i=0 1 − 2M ≤ k < 0, 0 |k| ≥ 2M. 1 2M

(7) kTc ) RB/L ( 2M

= rk From (6) and (7), it can be noticed that kTc and that RB/L is piecewise linear between the values at 2M . Therefore, the shape of RB/L can be changed by changing the shape vectors of those two local step-shape modulated symbols The proposed method in this paper uses two local signals with step-shape modulated symbols which noncoherently combine the outputs of the correlators to create a no-side-peak function instead of the original ambiguous ACF which contains 2 or more side-peaks. The CCF between a cosine-BOC(kn, n) signal and two local step-shape modulated signals will be calculated as following.

(2)

(1)

(1) The CCFs between received BOC signal and local SCS signals are piecewise linear, and they are symmetrical relative to the perpendicular axis. (2) The CCFs have the same zero crossing points and the same sign when |τ| > Ts as shown in Fig. 2; and the (1) (2) ACCF rk , rk follow the terms in (12). 

(2)

rk

1 = 2M (−1)i/2 di+k ∑2M−1−k i=0 (1) 1 = 2M (−1)(2M−i−1)/2 d2M−k−i−1 ∑2M−1−k i=0 (1) 1 = 2M (−1)i/2 di−k ∑2M−1−k i=0 (1) = r−k . (2)

(9)

Substituting (9) into (6), we have RB/L2 (τ) = RB/L1 (−τ).

(10)

Just as the case in [12], a generic solution (i.e. the specific shape vectors of those two local step-shape modulated symbols) for high order BOC(kn, n) signal is hard to be given. For simplicity, a specific family of modulated symbol of local signal is dedicated to thereinafter. For cosine-BOC(kn, n) signal, the shape vectors of two local step-shape modulated symbols can be expressed as:   2M 2M T d1 (k) = [ 1+k 2 , 0, ..., 0, k 2] ,   1+k (11) 2M 2M T ] . d2 (k) = [k 1+k 2 , 0, ..., 0, 1+k2 Note that these two step-shape symbols have been energy normalized and the shapes are identified by the shape (1) (1) factor k = d2M−1 /d0 , k ∈ [0, 1). After substituting (11) into (7), it is interesting to find that

 W





 W

   

 



    



         

(8)

The correlation functions between the received cosineBOC(kn, n) signal and these two local signals can be denoted as RB/L1 (τ) and RB/L2 (τ), respectively. Note that (−1)i/2 = (−1)(2M−i−1)/2 , i = 0, 1, ..., 2M − 1. Substituting it into (7) yields

i = 0.

As analyzed above, we can draw the following conclusions.

The modulated symbol of one local signal is denoted (1) (1) (1) as pSCSc (t, d1 , fc ), where d1 = [d0 , d1 , ..., d2M−1 ]T , and the modulated symbol of another local signal, denoted as pSCSc (t, d2 , fc ), is the mirror image of pSCSc (t, d1 , fc ), i.e., dk = d2M−k−1 .

(12)









    



    







  

 





W

Fig. 2. Schematic diagram of correlation function RB/Li (τ).

Utilizing the characteristics of CCFs and the relationships shown in (10) and (12), it is possible to obtain a Synthesized Correlation Function (SCF) that completely removes any ACF side-peak. This eliminates the ambiguity threat in the acquisition and tracking process. Similar to [12], the SCF can be chosen as  2  1+k R(τ) = (RB/L1 + RB/L2 )2 − · (RB/L1 − RB/L2 ) . 1−k (13) In order to realize the acquisition and tracking process for both sine- and cosine-BOC(kn, n) signals in the same loop without changing the receiver architecture, the SCF for sine-BOC(kn, n) signal cannot be changed. Therefore, the modulated symbols of local signals for sine-BOC(kn, n) signal are chosen as pSCSs (t; d1 ; fc ) and pSCSs (t; d2 ; fc ), wherein (2)

(1)

dk = −dM−k−1 .

(14)

Substituting (14) into the ACCF in [12], it can be seen that (2)

(2)

i = M1 ∑M−1 i=0 (−1) di+k 1 M−1 1−M+2i M−i−1 (−d (1) = (−1) M−i−k−1 ) M ∑i=0 (−1) (1) 1 M−1 i = M ∑i=0 (−1) di−k (1) = r−k . (15) Corresponding to cosine-BOC(kn, n) signals, for sineBOC(kn, n) signal, we have RB/L2 (τ) = RB/L1 (−τ); and the shape vectors of two local step-shape modulated symbols, similar with the shape vectors given in (11), are chosen as

rk

RADIOENGINEERING, VOL. 24, NO. 3, SEPTEMBER 2015





M M T 2 , 0, ..., 0, −k 2] , 1+k  1+k M M T [k 1+k ] . 2 , 0, ..., 0, − 1+k2

d1 (k) = [ d2 (k) =

843

(16)

3.5

SCF,k=0.1 R (τ),k=0.1

3

B/L

1

Considering the fact that MBOC signal is the combination of two sine phased BOC signals, and most of the energy are located at sine-BOC(1,1) signal, so that the ACF of the MBOC signal is mainly decided by sine-BOC(1,1) signal. Therefore, the shape vectors of the two local step-shape modulated symbols are chosen like sine-BOC(1,1) signal, i.e., d1 (k) = ⎡

⎤T

⎢ 2 ⎢ ⎣ 1 + k2 , ..., ! "#

2 , −k 1 + k2 $!

6

2 , ..., −k 1 + k2 "#

2

2.5

SCF,k=0.5 RB/L (τ),k=0.5 1

2

R

B/L

1.5

1

0.5

0

−0.5

−1 −1.5

−1

−0.5

0

0.5

1

Fig. 3. RB/Li and SCF for MBOC signal, k = 0.1 and k = 0.5, with 6 MHz frond-end filter (double-sided).

2 ⎥ ⎥ 1 + k2 ⎦ $

1.2

SCF,6MHz RB/L (τ),6MHz

1

1

R

2

⎤T

SCF,Infinite R (τ),Infinite B/L

0.6

6

6

(17) Figures 3 to 5 show the SCFs of the proposed method and the impact of a front-end filter with a 6 MHz doubleside bandwidth on the SCFs for MBOC, sine-BOC and cosine-BOC, respectively. As can be seen, the low-pass filter changes the shapes of the correlation functions RB/L1 and RB/L2 . Note that the SCFs do not have any positive side peak. Therefore, the proposed method can completely remove the ambiguity threat for various BOC signals. Meanwhile, the main peak of the SCFs maintains its sharpness which translates into good tracking performance. Two negative side-peaks adjacent to the main peak are an artefact of unmatched slopes (refer to (12)) between the two correlation functions when |τ| < Ts as shown in Fig. 2. However, by using the special designed EMLP discriminator for the proposed method which will be illustrated in Section 4, these side-peaks will not maintain a tracking lock as they are negative valued [16].

1

(τ),Infinite

B/L

2

Correlation Function

2 , ..., − 1 + k2 "#

R 0.4

Sine−BOC(2n,n),6MHz

0.2

0

−0.2

−0.4

−0.6

−0.8 −1.5

−1

−0.5

0

0.5

1

1.5

Code Delay (chips)

Fig. 4. RB/Li SCF and ACF for sine-BOC(2, 1) signal, k = 0.1, with front-end filter 6 MHz and infinite bandwidth. 0.8

SCF,6MHz RB/L (τ),6MHz 1

R

B/L

0.6

(τ),6MHz

2

SCF,Infinite RB/L (τ),Infinite 1

R

B/L

0.4

Correlation Function

2 ,− 1 + k2 $!

2 ⎥ ⎥ 1 + k2 ⎦ $

(τ),6MHz

B/L

0.8

2 , ..., k 1 + k2 "#

1.5

Code Delay (chips)

6

d2 (k) = ⎡ ⎢ ⎢k ⎣ !

(τ),k=0.5

2

Correlation Function

Substituting (16) into the ACCF in [12], the same conclusions as (12) can be drawn for sine-BOC(kn, n) signals. Therefore, the SCF R(τ) is also suitable for sine-BOC(kn, n) signal when the shape vectors (16) are adopted for sinephased BOC signals.

RB/L (τ),k=0.1

(τ),Infinite

2

0.2

0

3. Acquisition using Proposed Method −0.2

3.1 Detection Probability Analysis. Figure 6 illustrates the receiver acquisition structure of the proposed method, wherein (+/−)2 represents (I1p ± I2p )2 or (Q1p ± Q2p )2 , I1p , I2p and Q1p , Q2p are the in-phase and quadra-phase integrator outputs, respectively.

−0.4 −1.5

−1

−0.5

0

0.5

1

1.5

Code Delay (chips)

Fig. 5. RB/Li and SCF for cosine-BOC(2, 1) signal, k = 0.1, with front-end filter 6 MHz and infinite bandwidth.

FENG SHEN, GUANGHUI XU, JOON WAYN CHEONG, ET AL., UNAMBIGUOUS ACQUISITION AND TRACKING . . .

844



  

!



 

   #$



⎤T 0 0 R11 (0) R12 (0) 0 0 ⎥ N0 ⎢ ⎥ ⎢R21 (0) R22 (0) (23) σ2 = ⎥ ⎢ 4Tcoh ⎣ 0 0 R11 (0) R12 (0)⎦ 0 0 R21 (0) R22 (0) where Ri j is the autocorrelation (i = j) or cross-correlation (i = j) function of the local auxiliary signals, i.e., Ri j = E(Li ⊗L j ). Here, we use Monte Carlo simulations to analyze the acquisition performance and assume that the Doppler frequency error and code delay error are small [10]. Meanwhile, the noise is generated using (20). ⎡

   D  u

 

 

 

 #$

D



 

!   



 D  u   

 

Fig. 6. Block diagram of the proposed acquisition structure.

Since the discriminator is non-coherent, we can ignore the data modulation. Therefore, the received Intermediate Frequency (IF) signal can be modeled as r(t) = Aexp(2π fIF t + θ)s(t − τ) + n(t)

(18)

where A is the signal amplitude, fIF is the intermediate frequency of the receiver, θ is the phase of the carrier, τ is the propagation code delay, s(t) is the baseband direct sequence spread spectrum signal, n(t) is the band-limited white noise with zero mean and double-sided power spectral density N0 , respectively. As shown in Fig. 6, the test criterion of the proposed method can be expressed as %

& N (IB/L1 ,i + IB/L2 ,i )2 + (QB/L1 ,i + QB/L2 ,i )2 − Tproposal = ∑

2 2 i=1 α (IB/L1 ,i − IB/L2 ,i ) + (QB/L1 ,i − QB/L2 ,i ) (19) 1+k 2 where α = ( 1−k ) , k is the shape factor. IB/L and QB/L are the in-phase and quadrature CCF results between the incoming BOC signal and the local auxiliary signal, respectively, which are expressed as

IB/L j ,i QB/L j ,i

= (A/2)sinc(πΔ f Tcoh )RB/L j (Δτ)cos(Δϕ) + nI j ,i = (A/2)sinc(πΔ f Tcoh )RB/L j (Δτ)sin(Δϕ) + nQ j ,i (20)

where j = 1, 2, sinc(x) = sin(x)/x. Tcoh is the integration time. Δ f , Δτ, Δϕ refer the frequency wipe-off error, the code delay and carrier phase estimation errors, respectively. Since nI j ,i and nQ j ,i are generated by the same noise process, its probability density function (pdf) follows a Gaussian distribution [10]. Therefore, at a code delay of Δτ, the joint distribution of the four branches is (IB/L1 (Δτ), IB/L2 (Δτ), QB/L1 (Δτ), QB/L2 (Δτ)) ∼ N(μ, σ2 ) (21) where ⎡ ⎤T RB/L1 (Δτ)cos(Δϕ) ⎢ ⎥ ⎢ RB/L2 (Δτ)cos(Δϕ)⎥ ⎥ , (22) μ = (A/2)sinc(πΔ f Tcoh ) ⎢ ⎢R ⎥ ⎣ B/L1 (Δτ)sin(Δϕ) ⎦ RB/L2 (Δτ)sin(Δϕ)

The detection probability of MBOC, sine-BOC(2n,n) and cosine-BOC(2n,n) versus CNR with a fixed false alarm probability Pf a = 10−5 for different acquisition methods are shown in Fig. 7, Fig. 8 and Fig. 9, respectively. The test threshold can be asymptotically tuned from a straightforward evaluation of the Gaussian integral for fixed Pf a , under the H0 hypothesis (i.e., absence of the tested code with the considered code offset) [17]: V = E[Zk (τ)|H0 ] + (2 · var[Zk (τ)|H0 ])1/2 · er f −1 (1 − 2Pf a ) (24) where

1 W (25) ∑ R(τ) W w=1 where W is the accumulation blocks number of the power detector; R(τ) is the synthesized correlation function of the proposed method. Zk (τ) =

In the simulation, the detection probability of traditional ambiguous acquisition is under the assumption that no false acquisition occurs and it is reviewed here for comparison. From Fig. 7, Fig. 8 and Fig. 9, it can be observed that with the increase of non-coherent summation number N and the integral time Tcoh , the detection probability increase significantly. However, the non-coherent summations can widen the performance gap between the traditional acquisition scheme and the proposed method as shown in Fig. 8 and Fig. 9. The performance degradation is significant with a larger shape factor k which is proportional to the weight coefficient α. Since the weighted sum operation in (19) introduces more noise relative to the traditional method, the performance degradation seems to increase with the increasing of weight coefficient α. However, it must be kept in mind that the comparison is based on the assumption that no false acquisition occurred for traditional acquisition scheme. From Fig. 7, it also can be seen that for MBOC signal, the proposed method outperforms the traditional acquisition scheme in terms of detection probability when the non-coherent summation number is low. That is mainly because the main peak of the SCF of the proposed method is the sum of two CCFs between the incoming signal and two local auxiliary signals. For MBOC with k = 0.1, the signal to noise ratio (SNR) of the proposed method is slightly higher than the traditional method. However, with the increase of modulation order M of BOC signals, the amplitude of the SCF peak decreases, resulting in degradation in detection performance; this degradation is more obvious

RADIOENGINEERING, VOL. 24, NO. 3, SEPTEMBER 2015

for cosine phased BOC signals comparing to Fig. 3, Fig. 4 and Fig. 5. However, it must be kept in mind that the comparison is based on the assumption that no false acquisition occurred. In practical applications, the multi-peak characteristic makes it easy for false detection to occur in the traditional acquisition scheme, especially in low C/N0 environment. But the SCF of the proposed method does not contain any positive side peak, which makes the proposed technique more robust against false detection than traditional acquisition scheme. Meanwhile, considering that longer integral time leads to higher detection probability, the proposed method can be used for acquisition on data-less channels where the coherent integral time can be prolonged to compensate for its lack of detection performance.

1

Proposal Traditional 0.9

0.8

Detection Probability

0.7

0.6

0.5

Tcoh=10ms N=1

Tcoh=1ms N=1

0.4

845

Tcoh=1ms N=20 0.3

0.2

0.1

25

30

35

40

45

50

C/N0 (dB−Hz)

Fig. 7. The detection probability for MBOC, Pf a = 10−5 , k = 0.1. 1

0.9

Traditional Proposal,k=0.1 Proposal,k=0.5

It also can be seen form Fig. 3, Fig. 4 and Fig. 5 that the test criterion of the proposed method has a narrower peak comparing with BPSK-Like [18] [19] methods, which can provide sufficiently small code delay uncertainty. Therefore the additional transition from acquisition to tracking can be avoided. However, it also means that the proposed method needs smaller code delay bins which will increase mean acquisition time.

0.8

Detection Probability

0.7

4. Tracking Performance Analysis of the Proposed Method

0.6

0.5

Tcoh=1ms N=1

Tcoh=1ms N=20

4.1 S-curve Simulation of the Proposed Technique

0.4

0.3

0.2

0.1

0 25

30

35

40

45

50

55

C/N0 (dB−Hz)

Fig. 8. The detection probability for sine-BOC(2n,n), Pf a = 10−5 , k = 0.1 and k = 0.5. 1

0.9

Traditional Proposal,k=0.1 Proposal,k=0.3

0.8

(26)

Detection Probability

0.7

Tcoh=1ms N=20

0.6

0.5

Tcoh=1ms N=1

0.4

0.3

0.2

0.1

0 25

In order to evaluate the tracking performance, the Early-Minus-Late-Power (EMLP) discriminator is used hereafter. For the proposed method, the discriminator can be obtained as ⎡' ( ⎤ L L 2 L L 2 + I ) + (Q + Q ) (IB/L B/L B/L B/L 1 2 1 2 ⎥ ⎢ ' (⎦ DELMP proposal = ⎣ E E 2 E E 2 − (IB/L1 + IB/L2 ) + (QB/L1 + QB/L2 ) ⎡' ( ⎤ L L 2 L L 2 − I ) + (Q − Q ) (IB/L B/L2 B/L1 B/L2 1 ⎥ ⎢ (⎦ −α⎣ ' E E 2 E E 2 − (IB/L1 − IB/L2 ) + (QB/L1 − QB/L2 )

30

35

40

45

50

55

C/N0 (dB−Hz)

Fig. 9. The detection probability for cosine-BOC(2n,n), Pf a = 10−5 , k = 0.1and k = 0.3.

Figure 10 shows EMLP discriminator output of sineBOC(2n, n) signal for the traditional [21-23], PCF [12] and proposed method without front-end filter and Early-Late spacing d = 0.1 chips, respectively,.. It can be seen that the proposed method clearly removes any potential false lock points [16] [24]. Moreover, the discriminator stability domain is slightly larger than PCF and the traditional method: from [−0.16; 0.16] chips for PCF and [−0.15; 0.15] chips for the traditional method to [−0.188; 0.188] chips for the proposed method. In addition, with the increase of shape factor k, the slope of the proposed method at τ ≈ 0 increased, but the discriminator stability domain remains unchanged.

FENG SHEN, GUANGHUI XU, JOON WAYN CHEONG, ET AL., UNAMBIGUOUS ACQUISITION AND TRACKING . . .

846

2

0.06 PCF,k=0 Traditional sine−BOC(2n,n) Proposal,k=0.1 Proposal,k=0.5 Traditional BPSK

Proposal,k=0.1 Proposal,k=0.5 Traditional PCF,k=0.1

1.5

0.05

Code Tracking Error Standard Deviation (chips)

Discriminator Output

1

0.5

0

−0.5

0.04

0.03

0.02

−1

0.01 −1.5

−2 −1

−0.8

−0.6

−0.4

−0.2

0

0.2

0.4

0.6

0.8

1

Input Code Delay (chips)

Fig. 10. S-curve of sine-BOC(2n,n) for traditional method, PCF with k = 0.1, and proposed method with k = 0.1 and k = 0.5, without front-end filter, and Early-Late spacing d = 0.1 chips.

4.2 Impact of Thermal Noise The BOC signals tracking accuracy depends on the thermal noise level, and the code tracking error standard deviation caused by thermal noise is defined as [25] √ 2BL Tcoh σD (27) στ = G where σD is the measured standard deviation of the discriminator output, BL is the single-sided closed loop noise bandwidth in Hz, Tcoh is the coherent integration time, G = dD(τ) dτ |τ=0 is the normalized discriminator gain. Figure 11 shows the standard deviation of code tracking error for different tracking technique with BL = 2 Hz, Tcoh = 1 ms, d = 0.05 chips. From Fig. 11, it can be seen that the proposed method has a worse performance than the traditional DLL and PCF (k = 0). However, there is only slight absolute gap between the proposed method and the traditional DLL with the CNR increasing. Although the tracking accuracy of the proposed method degraded comparing with the traditional DLL for BOC signals in the presence of thermal noise, it still outperforms the traditional DLL for BPSK signals. It also can be seen that the tracking accuracy curves of the proposed method for k = 0.1 and k = 0.5 nearly overlap, which means that the shape factor k has no impact on the thermal noise caused code tracking error.

0 20

22

24

26

28

30 C/N0 (dB−Hz)

32

34

36

38

40

Fig. 11. Code tracking error standard deviation versus C/N0 with BL = 2 Hz, Tcoh = 1 ms, d = 0.05 chips for sineBOC(2n,n) signal.

Figure 12 and Figure 13 show the multipath-induced error envelop for the traditional MBOC and BPSK tracking method, as well as for the proposed method with k = 0.1 and k = 0.5, respectively. In this simulation, the reflection path is 6 dB weaker than the direct path, and the Early-Late spacing, d, is 0.1 chips, and the frond-end bandwidth is infinite. Figure 14 and Figure 15 show the case for sine-BOC(2n,n) and cosine-BOC(2n,n) respectively with infinite front-end bandwidth, d = 0.1 chips, and several values of k. As shown in Fig. 12 and Fig. 13, the proposed method exhibits inferior multipath resistance performance in short and medium delay for MBOC signals. For delay [0.08; 0.51] chips, the proposed method has inferior performance than traditional BPSK signal. But for delay [0.56; 0.85] chips, it is superior to the traditional DLL for MBOC signals. For medium and long delay multipath, it exhibits familiar multipath mitigation performance with the traditional DLL for MBOC signals. However, it must be kept in mind that the proposed method completely removes ambiguity threat for both MBOC signals tracking and acquisition process.

4.3 Impact of Multipath

From Fig. 14 and Fig. 15, it can be seen that, compared with the traditional DLL, the multipath mitigation performance of proposed method decreased slightly for sineBOC(2n,n) and cosine-BOC(2n,n) signals. There is only slight absolute gap between the proposed method and the traditional DLL for cosine-BOC(2n,n) signal with the multipath delay increasing. But the proposed method outperforms the traditional DLL for BPSK signals significantly, especially for long delay multipath.

To investigate the effect of multipath on the code tracking, in this paper, we consider a simple model of multipath with only one-path specular reflection, Meanwhile, we assume the amplitude and phase delay of the multipath are time-invariant over the time period of interest.

It should be noted that the multipath induced error envelops of the proposed method for k = 0.1 and k = 0.5 are nearly overlap as shown from Fig. 12 to Fig. 15, which means the changing of shape factor k cannot improve the multipath resistance performance. The conclusion that the

RADIOENGINEERING, VOL. 24, NO. 3, SEPTEMBER 2015

847

0.04

0.025

BPSK MBOC Proposal,k=0.5 Proposal,k=0.1

0.03

BPSK sine−BOC(2n,n) Proposal, k=0.1 Proposal, k=0.5

0.02

Running Everage Error (chips)

Multipath Error Envelop (chips)

0.02

0.01

0

−0.01

0.015

0.01

−0.02 0.005

−0.03

−0.04

0

0.5

1

0

1.5

0

0.5

1

1.5

Multipath Delay (chips)

Multipath Delay (chips)

Fig. 12. Code tracking multipath envelop for MBOC, BPSK, and proposed method (k = 0.1 and k = 0.5) for EarlyLate spacing of 0.1 chips, and infinite front-end filter.

Fig. 14. Running average errors of BPSK, sine-BOC(2n,n) and proposed method (k = 0.1 and k = 0.5), Early-Late spacing of 0.1 chips, and infinite front-end filter.

0.03

0.025

Proposal, k=0.1 Proposal, k=0.5 BPSK MBOC

BPSK cosine−BOC(2n,n) Proposal, k=0.5 Proposal, k=0.1

0.025

0.02

Running Everage Errors (chips)

Running Everage Error (chips)

0.02

0.015

0.01

0.015

0.01

0.005 0.005

0

0

0.5

1

0

1.5

0

0.5

1

1.5

Multipath Delay (chips)

Multipath Delay (chips)

Fig. 13. Running average errors of BPSK, MBOC and proposed method (k = 0.1 and k = 0.5), Early-Late spacing of 0.1 chips, and infinite front-end filter.

Fig. 15. Running average errors of BPSK, cosine-BOC(2n,n) and proposed method (k = 0.1 and k = 0.5), Early-Late spacing of 0.1 chips, and infinite front-end filter.

shape factor k has no influence on multipath resistance is similar with the thermal noise resistance performance mentioned above. However, a lower value of k remains a relative higher detection probability for the proposed method. Therefore, a lower k, such as k = 0.1, is suggested in the practical applications.

In fact, the shape vector of local signal can be divided (1) (2) into two parts, i.e. dL = αdL ± βdL . Then the ACCF between BOC signal and local auxiliary signal can be expressed as rk

= = =

5. Implementation Issues In the proposed method, the local modulated symbols are step-shape and k is a real tunable variable within [0, 1), so there are some difficulties to implementing the correlators and proposed modulators using logic circuits. To solve the problem, we herein reanalysis the SCF of the proposed method.

∑M−1 i=0 dBOC,i dL,k+i (1) (2) ∑M−1 i=0 dBOC,i (αdL,k+i ± βdL,k+i ) (1) (2) αrk ± βrk . 1 M 1 M

(28)

Note that the CCF between received BOC signal and local auxiliary signal is a linear combination of ACCF, therefore, after substituting (28) into the CCF, for example (6), then (1) (2) (29) RB/L (τ) = αRB/L (τ) ± βRB/L (τ) (1)

(2)

where RB/L (τ) and RB/L (τ) are the CCFs between received BOC signal and local auxiliary signals with shape vectors

FENG SHEN, GUANGHUI XU, JOON WAYN CHEONG, ET AL., UNAMBIGUOUS ACQUISITION AND TRACKING . . .

848

(1)

(2)

dL , dL respectively. In (11), the shape vectors of local auxiliary signal d1 (k) and d2 (k) can be rephrased as d1 (k) = d2 (k) = where

(30)



2M , 0, ..., 0]T 1+k2  2M T [0, ..., 0, 1+k 2]

d (1) (k) = [ d (2) (k) =

then

d (1) (k) + kd (2) (k) kd (1) (k) + d (2) (k)

(1)

(2)

RB/L1 (τ) = RB/L (τ) + kRB/L (τ) (1)

(31)

(2)

RB/L2 (τ) = kRB/L (τ) + RB/L (τ).

(32)

It can be proved that the shape of the CCF between received BOC signal and local auxiliary signal keeps unchanged. However, the spreading symbols in this form are easier to implement using logic circuits. In this form, after the carrier stripping process, the received signal does not need to be correlated with two stepshape symbol modulated local signals. It just multiplies a replica of the spreading code and then does gating integral and dump. RB/L1 (τ) and RB/L2 (τ) can be obtained as a linear combination of one integrator output which only works during the first Ts interval of each Tc and another integrator output which only works during the last Ts interval of each Tc . The integration time is T . Such a high speed but simple process is easily implemented in an ASIC or FPGA, and this process provides filtering and resampling.

can maintain the DLL stability in low C/N0 environment. In term of thermal noise resistance performance, the proposed method has very limited degradation compared with PCF and traditional tracking method with the C/N0 increasing. The multipath rejection performance of the proposed method is decreased compared with the traditional tracking method, especially for MBOC signals, meanwhile, the degradation decreases with the increase of modulation order, especially for cosine phased BOC signals. It must be kept in mind that the comparison is based on the assumption that no false acquisition and tracking occurred for the traditional method. But in practical applications, the multi-peak characteristic makes the traditional acquisition and tracking schemes easy to acquire or track false peak. Concerning the shape factor k, it is deserved to note that the shape factor k has an obvious effect on BOC signals acquisition process but little impact on signals tracking accuracy. A large value of k will amplify the noise, resulting in decreased detection probability. Therefore, a low value of k is suggested to be chosen if the proposed method is used in the receiver. In addition, the proposed method is practicable and requires the same number of correlators as ASPeCT and the bump-jumping method. Moreover, it has an equivalent form to the proposed method whose correlators and modulators are easily implemented via logic circuits.

Note that the number of correlators the proposed technique employs is the same as the bump-jumping method and ASPeCT, which is small compared with those which require dozens of correlators to obtain an unambiguous discriminator S-curve, such as S-curve shaping technique proposed in [26].

This work was supported by National Nature Science Foundation of China under grants 61102107 and 61374208, and by the Fundamental Research Funds for the Central Universities under grants HEUCFX41310.

6. Conclusions

References

In this paper, based on the combination of constitutive sub-correlations of the BOC autocorrelation, an unambiguous acquisition and tracking technique for generic BOC signal is proposed, which is suitable for MBOC, sine phased and cosine phased BOC signals. Theoretic analysis and simulation results show that the proposed method completely removes ambiguity problem in both BOC signals acquisition and tracking process. In the acquisition process, for MBOC signals, the detection probability of the proposed method keeps higher than the traditional acquisition method if a low value of shape factor k is adopted. Whereas for higher modulation order sine and cosine phased BOC signals, the acquisition sensitivity of the proposed method decreases with the increase of shape factor k, but the sensitivity degradation of the proposed method is acceptable when the shape factor k keeps small. In the tracking process, the discriminator output of the proposed method keeps a slightly larger stability domain compared with the traditional tracking and PCF methods, which

Acknowledgments

[1] GODET, J., DE MATEO, J. C., ERHARD, P., et al. Assessing the radio frequency compatibility between GPS and Galileo. In Proceeding of U.S. Institute of Navigation GPS Conference. Portland (OR, USA), 2002, p. 1260–1269. [2] BETZ, J. W. Binary offset carrier modulations for radionavigation. Navigation, Journal of the Institute of Navigation, 2002, vol. 48, no. 4, p. 227–246. DOI: 10.1002/j.2161-4296.2001.tb00247.x [3] MARTIN, N., LEBLOND, V., GUILLOTEL, G., HEIRIES, V. BOC(x, y) signal acquisition technique and performances. In Proceedings of the 16th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GPS/GNSS). Portland (OR, USA), 2003, p. 188–198. [4] FISHMAN, P., BETZ, J. W. Predicting performance of direct acquisition for the M code signal. In Proceedings of the 2000 National Technical Meeting of The Institute of Navigation. Anaheim (CA, USA), 2000, p. 574–582. [5] MAO, W.-L., HWANG., C.-S., HUNG, C.-W., et al. Unambiguous BPSK-like CSC method for Galileo acquisition. In 18th International Conference on Methods and Models in Automation and Robotics (MMAR). Miedzyzdroje (Poland), 2013, p. 627–632. DOI: 10.1109/MMAR.2013.6669983

RADIOENGINEERING, VOL. 24, NO. 3, SEPTEMBER 2015

[6] BENEDETTO, F., GIUNTA, G., LOHAN, E. S., RENFORS, M. A fast unambiguous acquisition algorithm for BOC-modulated signals. IEEE Transactions on Vehicular Technology, 2013, vol. 62, no. 3, p. 1350–1355. DOI: 10.1109/TVT.2012.2228681 [7] FINE, P., WILSON, W. Tracking algorithm for GPS offset carrier signals. In Proceedings of the 1999 National Technical Meeting of The Institute of Navigation. San Diego (CA, USA), 1999, p. 671–676. [8] HODGART, M., BLUNT, P., UNWIN, M. Double estimator – a new receiver principle for tracking BOC signals. Inside GNSS, Spring 2008, p. 26–36. [9] YAO, Z., LU, M., FENG, Z. Unambiguous sine-phased binary offset carrier modulated signal acquisition technique. IEEE Transactions on Wireless Communications, 2010, vol. 9, no. 2, p. 577–580. DOI: 10.1109/TWC.2010.02.091066 [10] ZHOU, Y. Ambiguity mitigating technique for cosine-phased binary offset carrier signal. IEEE Transactions on Wireless Communications, 2012, vol. 46, no. 4, p. 1981–1984. DOI: 10.1109/TWC.2012.041612.110531 [11] KIM, H., LEE, Y., YOON, S. An unambiguous acquisition scheme for binary offset carrier signals. In 2013 International Symposium on Intelligent Signal Processing and Communications Systems (ISPACS). Naha (Japan), 2013, p. 340–344. DOI: 10.1109/ISPACS.2013.6704573 [12] YAO, Z., et al. Pseudo-correlation function-based unambiguous tracking technique for sine-BOC signals. IEEE Transactions on Aerospace and Electronic Systems, 2010, vol. 46, no. 4, p. 1782–1796. DOI: 10.1109/TAES.2010.5595594 [13] XU, D., LIU, M., SHEN, F. Ambiguity mitigating technique for multiplexed binary offset carrier signal tracking. IEEE Communications Letters, 2013, vol. 17, no. 11, p. 2021–2024. DOI: 10.1109/LCOMM.2013.091113.131685 [14] WARD, P. W. A design technique to remove the correlation ambiguity in binary offset carrier (BOC) spread spectrum signals. In Proceedings of the 59th Annual Meeting of The Institute of Navigation and CIGTF 22nd Guidance Test Symposium. Albuquerque (NM, USA), 2003, p. 146–155. [15] PHAM-VIET, H., DAO-NGOC, C., NGUYEN-VAN, K. A scheme of cosine-BOC(n,n) side peaks cancellation for navigation applications. In 2013 International Conference on Advanced Technologies for Communications (ATC’13). Ho Chi Minh City (Vietnam), 2013, p. 180–184. DOI: 10.1109/ATC.2013.6698101 [16] JULIEN, O., MACABIAU, C., CANNON, M. E., LACHAPELLE, G. ASPeCT: unambiguous sine-BOC(n,n) acquisition/tracking technique for navigation applications. IEEE Transactions on Aerospace and Electronic Systems, 2007, vol. 43, no. 1, p. 150–162. DOI: 10.1109/TAES.2007.357123 [17] BENEDETTO, F., GIUNTA, G. A self-synchronizing method for asynchronous code acquisition in band-limited spread spectrum communications. IEEE Transactions on Communications, 2006, vol. 57, no. 8, p. 2410–2419. DOI: 10.1109/TCOMM.2008.08.070030 [18] HEIRISES, V., ROVIRAS, D., RIES, L. Analysis of non-ambiguous BOC signal acquisition performance. In Proceedings of the 17th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS). Long Beach (CA, USA), 2004, p. 671–676. [19] LOHAN, E. S., BURIAN, A., RENFORS, M. Low-complexity unambiguous acquisition methods for BOC modulated CDMA signals. International Journal of Satellite Communications and Networking, 2008, vol. 26, p. 503–522. DOI: 10.1002/sat.922 [20] HOLMES, J. K. Coherent Spread Spectrum Communications. Melbourne (FL, USA): Krieger Publishing Co., 1990. ISBN: 9780071382151 [21] VAN DIERENDONCK, A. J. GPS receivers. In Global Positioning System: Theory and Applications. New York (USA): American Institute of Astronautics and Aeronautics, 1996, p. 329–408. ISBN: 978-1563471063

849

[22] WARD, P. Satellite signal acquisition and tracking. In Understanding GPS-Principles and applications. Norwood (MA, USA): Artech House, 1997, p. 119–218. ISBN: 978-1580538947 [23] HOLMES, J. K., RAGHAVAN, S. H., LAZAR, S. Acquisition and tracking performance of NRZ and square-wave modulated symbols for use in GPS. In Proceedings of the 54th Annual Meeting of The Institute of Navigation. Denver (CO, USA), 1998, p. 611–626. [24] JULIEN, O., MACABIAU, C., CANNON, M. E., et al. New unambiguous BOC(n,n) tracking technique. In Proceeding of NaviTec Conference. Noordwidjk (The Netherlands), 2004, 8 pages. [25] NUNES, F. D., SOUSA, F. M. G., LEITAO, J. M. N. Gating functions for multipath mitigation in GNSS BOC signals. IEEE Transactions on Aerospace and Electronic Systems, 2007, vol. 43, no. 3, p. 951–964. DOI: 10.1109/TAES.2007.4383585 [26] PANY, T., IRSIGLER, M., EISSFELLER, B. S-curve shaping: A new method for optimum discriminator based code multipath mitigation. In Proceedings of the 18th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS). Long Beach (CA, USA), 2005, p. 2139–2154.

About the Authors . . . Feng SHEN was born in Jiangsu,China in 1981. He received his Ph.D. from Harbin Engineering University in 2009. He is now an associate professor at the College of Automation, Harbin Engineering University, Harbin, China. He developed the signal acquisition and tracking of the Spread spectrum receiver for the GPS ground-based augmentation systems. From 2014 to 2015, he is a visiting academic with the School of Electrical Engineering, University of New South Wales. His research interests include weak signal acquisition and tracking, GNSS/INS integrated navigation system. Guanghui XU (corresponding author) was born in Shandong, China in 1987. He received his B.S. from Weifang University, Shandong, China in 2011. He is currently working toward the Ph.D. degree at Harbin Engineering University. His research interests include satellite navigation, signal structure design and signal processing. Joon Wayn CHEONG is a research associate at the School of Electrical Engineering, University of New South Wales (UNSW). He received his Ph.D. (2012) and B.E. (2008) from UNSW. He developed the firmware for the space-qualified Namuru family of GPS/GALILEO integrated receivers under the Garada, Biarri and QB50 project, where he cracked the Locata pseudolite positioning system’s code and derived a hybrid Collective Detection algorithm. His other research interests in the GNSS field include weak signal acquisition, A-GPS, GNSS/pseudolite integrated signal processing, GNSS interference and spoofing. He authored more than 15 peer-reviewed publications. Haiyu FENG was born in Shandong, China in 1989. She received her B.S. from Binhai University, Shandong, China in 2011. She is currently working toward the M.Sc. degree at Harbin Engineering University. Her research interests include satellite navigation, target detection under sea clutter.