Open AccessArticle

Frequency-Locked Detector Threshold Setting Criteria Based on Mean-Time-To-Lose-Lock (MTLL) for GPS Receivers

Tian Jin

Heliang Yuan

Na Zhao

¹,

Honglei Qin

^1,*,

Kewen Sun

² and

Yuanfa Ji

School of Electronic and Information Engineering, Beihang University, Beijing 100191, China

School of Computer and Information, Hefei University of Technology, Hefei 230009, China

Guangxi Key Laboratory of Precision Navigation Technology and Application, Guilin 541004, China

Author to whom correspondence should be addressed.

Sensors 2017, 17(12), 2808; https://doi.org/10.3390/s17122808

Submission received: 17 October 2017 / Revised: 21 November 2017 / Accepted: 28 November 2017 / Published: 4 December 2017

(This article belongs to the Special Issue Sensor Signal and Information Processing)

Download

Browse Figures

Versions Notes

Abstract

Frequency-locked detector (FLD) has been widely utilized in tracking loops of Global Positioning System (GPS) receivers to indicate their locking status. The relation between FLD and lock status has been seldom discussed. The traditional PLL experience is not suitable for FLL. In this paper, the threshold setting criteria for frequency-locked detector in the GPS receiver has been proposed by analyzing statistical characteristic of FLD output. The approximate probability distribution of frequency-locked detector is theoretically derived by using a statistical approach, which reveals the relationship between probabilities of frequency-locked detector and the carrier-to-noise ratio (C/N₀) of the received GPS signal. The relationship among mean-time-to-lose-lock (MTLL), detection threshold and lock probability related to C/N₀ can be further discovered by utilizing this probability. Therefore, a theoretical basis for threshold setting criteria in frequency locked loops for GPS receivers is provided based on mean-time-to-lose-lock analysis.

Keywords:

GPS receiver; frequency locked loop; frequency-locked detector; mean-time-to-lose-lock

1. Introduction

In the Global Positioning System (GPS) receiver, the lock detector is usually used to indicate the signal tracking status by comparing with a threshold. Code, frequency and phase are tracked by delay, frequency, and phase locked loops, respectively. These loops have their own indicators [1,2]. The GPS receiver mainly relies on phase locked loop (PLL) and delay locked loop (DLL) to track signal and frequency locked loop (FLL) is a transition part to bridge acquisition and PLL + DLL. However, in high performance GPS receivers, such as high dynamic and high sensitivity receivers, the FLL has been widely used when the carrier phase could not be tracked by PLL or aids PLL to track signal [3,4]. It indicates the significance of analyzing frequency-locked detector. For example, a detailed discussion and comparison of PLL and FLL can be found in [5]. Yang and Huang [6] proposed a non-linear carrier NCO unit to track carrier precisely by selecting interpolating filter orders, which is derived from the frequency discriminator in high dynamic situation. Curran et al. [7] discussed the design and steady-state performance of one first-order, two second-order FLL loop filters, and four carrier frequency discriminators, then a new FLL design under weak signal conditions was proposed. Meanwhile, FLLs have some distinct advantages over PLL counterparts [8]. By neglecting absolute phase error and permitting relative phase rotation of the received signal and the local carrier replica, an FLL can, typically, acquire and track signals that are at higher frequency offsets than a PLL. It is easy to understand that FLL plays a key role in facilitating reliable signal tracking. The frequency-locked detector (FLD) is used to indicate FLL signal tracking state. Natali [9] presented several AFC loops and performance for different signal. Messerschmi [10] has described two simply implemented frequency detectors to aid PLL timing and carrier acquisition, even with very small loop bandwidths and large initial frequency offsets. Mileant and Hinedi [11] proposed an FLD based on square law, which experiences less degradation due to phase jitter than the absolute value detector. Linn and Peleg [12] suggested a family of PLDs for M-PSK receivers operating in additive white Gaussian noise channels and derived the statistical properties of PLD. Kratyuk et al. [13] discussed Frequency detector for fast frequency lock of digital PLLs, which provides frequency difference information at each reference cycle to guarantee fast frequency acquisition. Previous researches mainly focus on the structure of FLL and its detector, but how to set the FLL detector threshold is seldom discussed, which make it hard to indicate the tracking frequency is stabilized or not. Our main purpose is to analyze the statistical performance of FLD by mean-time-to-lose-lock approaches, and make comparison with PLD, which has been studied by Jin et al. [14].

2. Frequency-Locked Detector Output $C_{2 φ}$

In the GPS receiver, the GPS signal received by the antenna is down converted to intermediate frequency through low noise amplifier, filters, and down-converter [15], the obtained signal

S (t)

can be written as follows:

S (t) = A C (t) D (t) \cos [(ω_{I F} + ω_{D o p p l e r}) t + ϕ_{0}] + n (t)

(1)

where

A

is the signal amplitude,

C (t)

is PRN code of GPS,

D (t)

is the navigation data,

ω_{I F}

is the intermediate frequency,

ω_{D o p p l e r}

is the Doppler frequency shift,

ϕ_{0}

is the initial carrier phase,

n (t)

is the additive white Gaussian noise.

According to [4,16,17], LOS (line of sight) signal through the front-end A/D sampling, quantization and integration, the in-phase, and quadrature-phase integrations at

i

th interval are I_i and Q_i, which are shown as:

\begin{matrix} I_{i} & \approx \sqrt{2 \frac{C}{N_{0}} T_{c o h}} R (τ) \sin c (π Δ f_{i} T_{c o h}) \cos (π Δ f_{i} T_{c o h} + Δ ϕ_{i}) + n_{I, i} \\ = \sqrt{2 \frac{C}{N_{0}} T_{c o h}} R (τ) \sin c (π Δ f_{i} T_{c o h}) \cos (Δ θ_{i}) + n_{I, i} \\ = A_{i} \cos (φ_{i}) \end{matrix}

(2)

\begin{matrix} Q_{i} & \approx \sqrt{2 \frac{C}{N_{0}} T_{c o h}} R (τ) \sin c (π Δ f_{i} T_{c o h}) \sin (π Δ f_{i} T_{c o h} + Δ ϕ_{i}) + n_{Q, i} \\ = \sqrt{2 \frac{C}{N_{0}} T_{c o h}} R (τ) \sin c (π Δ f_{i} T_{c o h}) \sin (Δ θ_{i}) + n_{I, i} \\ = A_{i} \sin (φ_{i}) \end{matrix}

(3)

where

C / N_{0}

is carrier to noise ratio,

T_{c o h}

is coherent integration time,

R (τ)

is the auto-correlation function of PRN code,

τ

is the code phase delay between received and local replica,

Δ f_{i}

and

Δ ϕ_{i}

are the frequency error between received and local replica and the initial phase error between received and local replica,

Δ θ_{i}

is total phase error that includes

Δ f_{i}

and

Δ ϕ_{i}

in the ith intervals.

n_{I, i}, n_{Q, i}

are independent white Gaussian noises corresponding to the I and Q branches, respectively. For the convenience of calculation, we take normalization processing of

n_{I, i} ~ N (0, 1), n_{Q, i} ~ N (0, 1)

[14].

A_{i}

and

φ_{i}

are signal amplitude and average phase estimation, containing noise component, relative to

\sqrt{2 \frac{C}{N_{0}} T_{c o h}} R (τ) \sin c (π Δ f_{i} T_{c o h})

and

Δ θ_{i}

. Detailed discussions of the effects of signal quantization are given in [18,19].

The output of frequency-locked detector during a period of time was proposed by [20], which is shown as:

C_{2 φ} = \frac{1}{M} \sum_{i = 1}^{M} \frac{D_{i}}{P_{i}}

(4)

where

D_{i} = {(D o t)}^{2} - {(C r o s s)}^{2}

P_{i} = {(C r o s s)}^{2} + {(D o t)}^{2}

D o t + j C r o s s = (I_{i} + j Q_{i}) (\bar{I_{i - 1} + j Q_{i - 1}})

. M is the times of accumulation.

Further, according to Spiegel [21], (4) can be simplified as:

D_{i} = {(D o t)}^{2} - {(C r o s s)}^{2} = [{(D o t)}^{2} + {(C r o s s)}^{2}] \cos (2 φ_{i} - 2 φ_{i - 1})

(5)

\begin{matrix} C_{2 φ} & = \frac{1}{M} \sum_{i = 1}^{M} \frac{D_{i}}{P_{i}} = \frac{1}{M} \sum_{i = 1}^{M} \frac{{(D {o t}_{i})}^{2} - {(C r o s s_{i})}^{2}}{{(D o t_{i})}^{2} + {(C r o s s_{i})}^{2}} \\ = \frac{1}{M} \sum_{i = 1}^{M} \frac{[{(D o t_{i})}^{2} + {(C r o s s_{i})}^{2}] \cos (2 φ_{i} - 2 φ_{i - 1})}{{(D o t_{i})}^{2} + {(C r o s s_{i})}^{2}} \\ = \frac{1}{M} \sum_{i = 1}^{M} \cos (2 φ_{i} - 2 φ_{i - 1}) \\ = \frac{1}{M} \sum_{i = 1}^{M} \cos (2 Δ φ_{i}) \end{matrix}

(6)

where

φ_{i - 1}

and

φ_{i}

are average phase estimation errors in the (i − 1)th and ith intervals, respectively.

Δ φ_{i}

is the difference between

φ_{i}

and

φ_{i - 1}

3. Distribution of FLD Output $C_{2 φ}$

According to Zhuang [22], the statistical properties of

φ

obeys Rician distribution with

| φ | \leq π

, its probability density function is

f (φ) = \frac{1}{2 π} e^{- \frac{C}{N_{0}} T_{c o h}} [1 + \sqrt{2 \frac{C}{N_{0}} T_{c o h}} \cos (φ - Δ θ) \exp (\frac{C}{N_{0}} T_{c o h} \cos^{2} (φ - Δ θ)) \int_{- \infty}^{a} e^{- \frac{x^{2}}{2}} d x]

(7)

where

a = \cos (φ - Δ θ) \sqrt{2 \frac{C}{N_{0}} T_{c o h}}

. The approximate simplification of (7) shows that

φ

is nearly subjected to Gaussian distribution under high

C / N_{0}

, while it obeys uniform distribution under low

C / N_{0}

. Based on the different distribution of

φ

, under different

C / N_{0}

conditions, we can analyze the distribution characteristics of FLD output

C_{2 φ}

individually.

3.1. $C_{2 φ}$ Distribution under Frequency Lock

When frequency locked loop is in lock and

C / N_{0}

is high, we can make the following assumptions:

\exp (- \frac{C}{N_{0}} T_{c o h}) \approx 0

(8)

\cos^{2} (φ - Δ θ) = 1 - \sin^{2} (φ - Δ θ) \approx 1 - {(φ - Δ θ)}^{2}

(9)

\int_{- \infty}^{\cos (φ - Δ θ) \sqrt{2 \frac{C}{N_{0}} T_{c o h}}} e^{- \frac{x^{2}}{2}} d x \approx \sqrt{2 π}

(10)

Based on the previous assumptions, (7) can be simplified to

f (φ) \approx \frac{1}{\sqrt{2 π * (\frac{1}{2 T_{c o h} C / N_{0}})}} * \exp (- \frac{{(φ - Δ θ)}^{2}}{2 * (\frac{1}{2 T_{c o h} C / N_{0}})}) | φ | \leq π

(11)

where

φ

follows Gaussian distribution with the mathematical expectation

Δ θ

and variance

1 / (2 T_{c o h} \cdot \frac{C}{N_{0}})

(

σ^{2} = 1 / (T_{c o h} \cdot \frac{C}{N_{0}})

). Then, we can get the probability distribution function of

Δ φ

shown below. When considering noise of ith and (i − 1)th are independent with each other,

φ_{i}

and

φ_{i - 1}

are independent of each other, while

Δ θ_{i}

and

Δ θ_{i - 1}

(mathematical expectation of

φ_{i}

and

φ_{i - 1}

) are approximately equal to zero under frequency lock and interval is small. Thus, we can calculate the mathematical expectation

E {Δ φ_{i}}

and variance

V A R {Δ φ_{i}}

, and probability distribution function of

Δ φ

, which are shown as follows:

E {Δ φ_{i}} = E {φ_{i}} - E {φ_{i - 1}} = Δ θ_{i} - Δ θ_{i - 1} \approx 0

(12)

V A R {Δ φ_{i}} = V A R {φ_{i} - φ_{i - 1}} = V A R {φ_{i}} + V A R {φ_{i - 1}} = σ^{2}

(13)

f (Δ φ_{i}) = \int_{- \infty}^{+ \infty} f (Δ φ_{i} - φ_{i - 1}, φ_{i - 1}) d φ_{i - 1} = \frac{1}{\sqrt{2 π} σ} \exp (- \frac{Δ φ_{i}^{2}}{2 σ^{2}})

(14)

The distribution characterization of FLD output represents a Gaussian stochastic process through a cosine system. We can get the mathematical expectation and variance of

C_{2 φ}

in (6) as follows:

\begin{matrix} E {C_{2 φ}} & = E {\frac{1}{M} \sum_{i = 1}^{M} \cos (2 Δ φ_{i})} \\ = \frac{1}{M} \sum_{i = 1}^{M} \int_{- \infty}^{\infty} \cos (2 Δ φ_{i}) p (Δ φ_{i}) d Δ φ_{i} \\ \approx \frac{1}{\sqrt{2 π} σ} \int_{- \infty}^{\infty} \cos (2 Δ φ_{i}) \exp (- \frac{Δ φ_{i}^{2}}{2 σ^{2}}) d Δ φ_{i} \\ = \exp (- 2 σ^{2}) \end{matrix}

(15)

\begin{matrix} V A R {C_{2 φ}} & = V A R {\frac{1}{M} \sum_{i = 1}^{M} \cos (2 Δ φ_{i})} \\ = \frac{1}{M^{2}} \sum_{i = 1}^{M} V A R {\cos (2 Δ φ_{i})} \\ = \frac{E {\cos^{2} (2 Δ φ)} - E^{2} {\cos (2 Δ φ)}}{M} \\ = \frac{1}{M} \frac{1}{\sqrt{2 π} σ} \int_{- \infty}^{\infty} \cos^{2} (2 Δ φ) \exp (- \frac{Δ φ^{2}}{2 σ^{2}}) d Δ φ - \frac{1}{M} \exp (- 4 σ^{2}) \\ = \frac{1}{2 M} (1 + e^{- 8 σ^{2}}) - \frac{1}{M} e^{- 4 σ^{2}} \end{matrix}

(16)

When the signal is locked,

Δ φ_{i}

is stable around zero. Therefore,

C_{2 φ}

approximates to

C_{2 φ} = \frac{1}{M} \sum_{i = 1}^{M} \cos (2 Δ φ_{i}) \approx 1 - \frac{2}{M} \sum_{i = 1}^{M} Δ φ_{i}^{2}

(17)

For simplification, we can assume a new variable Y related with

C_{2 φ}

, which nearly follows the chi-square distribution

Y = \frac{(1 - C_{2 φ}) M}{2 σ^{2}} = \sum_{i = 1}^{M} {(\frac{Δ φ_{i}}{σ})}^{2} ~ χ^{2} (M)

(18)

According to (18), the probability distribution function of FLD output is shown as:

f_{C_{2 φ}} (c_{2 φ}) = \frac{M}{2 σ^{2}} f_{Y} (\frac{M - M c_{2 φ}}{2 σ^{2}}) | c_{2 φ} | \leq 1

(19)

where

f_{Y} (y)

is the probability distribution function of the random variable

Y

(

Y ~ χ^{2} (M)

3.2. $C_{2 φ}$ Distribution under Frequency Unlock

Under low C/N₀ condition, frequency loop may be unlocked. After simplification of Equation (7),

φ

is uniformly distributed in

(- π, π)

, and its probability distribution function is shown as:

f (φ) \approx \frac{1}{2 π} | φ | \leq π

(20)

When M = 1, the probability distribution functions of

Δ φ

and

f_{C_{2 φ}} (c_{2 φ})

are shown follows, respectively:

f_{Δ φ} (Δ φ) = {\begin{matrix} \frac{1}{2 π} + \frac{Δ φ}{4 π^{2}} & (- 2 π, 0) \\ \frac{1}{2 π} - \frac{Δ φ}{4 π^{2}} & (0, 2 π) \end{matrix}

(21)

f_{C_{2 φ}} (c_{2 φ}) = {\frac{1}{π \sqrt{(1 + c_{2 φ}) (1 - c_{2 φ})}} | c_{2 φ} | \leq 1

(22)

According to (22), the expectation and variance of

C_{2 φ}

are provided:

E {C_{2 φ}} = \int_{- \infty}^{+ \infty} \cos (2 Δ φ) f (Δ φ) d Δ φ = 0

(23)

V A R {C_{2 φ}} = E {C_{2 φ}^{2}} - E^{2} {C_{2 φ}} = \int_{- \infty}^{+ \infty} \cos^{2} (2 Δ φ) f (Δ φ) d Δ φ = 0.5

(24)

From the above section, the distribution of FLD output in unlock situation is completely different from that in lock situation.

4. Comparison between Theoretical and Simulation Results of $C_{2 φ}$

The statistics characteristic of

C_{2 φ}

has been discussed in the previous sections. It provides a threshold setting that indicates the lock/unlock condition of FLD. In order to validate the validity of the theoretical analysis, we will compare between theoretical and simulation results.

The probability density functions of

C_{2 φ}

under lock and unlock conditions can be obtained from (19) and (20), respectively. At the same time, Monte Carlo simulation experiments that software receiver [23] acquires and tracks simulated GPS L1 signal are carried out. The parameters of the simulation experiment are shown in Table 1.

The simulated GPS signal with 4.996 MHz intermediate frequency is generated by the software program. Four signals with different C/N₀ values and a noise signal are simulated. The software receiver is used to process the signal and the probability density function (PDF) of FLD output can be retrieved under different situations.

The comparison between theoretical and simulation results are shown in Figure 1, Figure 2 and Figure 3. Figure 1 and Figure 2 show the PDF of

C_{2 φ}

with different average length M when signal presented. However, both of the figures indicate that maximum probability point of the horizontal axis moves closer to one when increasing the C/N₀ values, which means that

Δ φ_{i}

become smaller. Normally,

C_{2 φ}

decreases when reducing C/N₀. Maximum probability points in probability density curves change dramatically with different C/N₀ values in Figure 1 but no evident differences in Figure 2. It means that probability density curve could not indicate

C_{2 φ}

changing with different C/N₀ because of interference by noise when M = 1 in Figure 2. So, M should be set larger. Figure 3 shows

C_{2 φ}

obeys uniform distribution when no signal presented.

4.1. Lock Probability Analysis of FLD Output

The output distribution of FLD is mainly related with the C/N₀ of the received signal. With a given

C / N_{0} = λ

, the locked threshold of the tracking loop can be determined by the false alarm probability

P_{F A}

and the lock probability

P_{D}

P_{F A} = P {C_{2 φ} > T h | n o i s e i n p u t o n l y} = \int_{T h}^{1} \frac{1}{π \sqrt{(1 + c_{2 φ}) (1 - c_{2 φ})}} d c_{2 φ}

(25)

\begin{matrix} P_{D} & = P {C_{2 φ} > T h | C / N_{0} = λ} \\ = P {Y = \frac{M (1 - C_{2 φ})}{2 σ^{2}} < \frac{M (1 - T h)}{2 σ^{2}} | C / N_{0} = λ} \\ = F_{Y} (\frac{M (1 - T h)}{2 σ^{2}} | C / N_{0} = λ) \end{matrix}

(26)

where

F_{Y} (y)

is the distribution function of the random variable

Y

(

Y ~ χ^{2} (M)

). Figure 4, Figure 5 and Figure 6 describe the relationships between the lock probability and the false alarm probability. As shown in Figure 4 and Figure 5, the lock probability increases with the increasing coherent integration time under the same C/N₀. The lock probability can be also increased by decreasing the threshold, but the false alarm probability also increases under that situation. Figure 6 shows the relation between the FLD output false alarm probability and threshold. As a result, in order to set optimal FLD threshold, C/N₀ and false alarm should be considered at same time.

4.2. Mean-Time-To-Lose-Lock (MTLL) of FLD

In order to evaluate the performance of the tracking loop, we need to analyze the mean-time-to-lose-lock (MTLL) of the frequency tracking loop. The time should be shorter when the signals are not present, while longer when the signals are present. The distinct threshold will influence how long the tracking loop can detect loop state changes. To be more specific, MTLL (

E {T_{F}}

) is used to indicate the average time required that tracking loop detect state changed from false alarm to correct detection under no signal condition. MTLL (

E {T_{D}}

) will imply that how long tracking loop change the signal detection state to missing detection state when the signal presented. Over a total time TT, MTLL (

E {T_{F}}

) under false alarm

P_{F A}

can be obtained [14] as:

\begin{matrix} E {T_{F}} |_{T T - > \infty} & = \sum_{N = 1}^{T T} N P_{F A}^{N} (1 - P_{F A}) |_{T T - > \infty} \\ = \frac{1}{1 - P_{F A}} \end{matrix}

(27)

where the false alarm probability

P_{F A}

with a certain threshold (M = 1) is shown in (25), N is time slot. MTLL (

E {T_{F}}

) indicates average time that frequency locked loop used to detect no signal state.

Based on (25) and (27), we can get the relationship between MTLL (

E {T_{F}}

) and lock threshold. As shown in Figure 7, it indicates that relation between MTLL (

E {T_{F}}

) and threshold under no signal situation. The continuous tracking time is gradually decreased when increasing

T h

. MTLL (

E {T_{F}}

) is almost same when

T h

above zero. The MTLL (

E {T_{F}}

) of false alarm in FLD shown in Figure 7 is similar with that in PLD in [14].

When the signals present, the MTLL (

E {T_{D}}

) can be obtained [14] as:

E {T_{D}} |_{T T - > \infty} = \frac{1}{1 - P_{D}}

(28)

where the lock probability

P_{D}

with a certain threshold

T h

is shown in (26). MTLL (

E {T_{D}}

) indicates the average time that the frequency locked loop takes to maintain signal detection state.

Based on (26) and (28), we can get the relationship between MTLL (

E {T_{D}}

) and C/N₀ with a certain lock threshold. The relationship curves are shown in Figure 8 and Figure 9 with M = 1 and M = 20, the continuous tracking time MTLL (

E {T_{D}}

) gradually increases when increasing C/N₀ or coherent integration time. At the same time, the verification of MTLL is conducted with the help of Monte Carlo simulation in Figure 10 and Figure 11.

The theoretical results show that the MTLL (

E {T_{F}}

) is only 2 s when the signal does not present and the threshold

T h

is equal to 0 in Figure 7. Oversized noise can cause false alarm that leads to error lock in frequency locked loop. We can judge the error lock state in is 2 s. The MTLL (

E {T_{D}}

) is about 40 s when the integration time is 1 ms under signal with 40 dB-Hz C/N₀, and 10 ms for the signal with 30 dB-Hz C/N₀ that threshold

T h

is equal to 0 (M = 1) in Figure 8 when the signal present.

The comparison with PLD [14] is shown in Figure 12. The threshold setting of PLD at 0.7 has a similar curve to the threshold setting of FLD at 0.4, while PLD at 0.6 is similar to the FLD at 0.2. It means that FLD has lower threshold than PLD under similar MTLL. Therefore, FLD can track signal more stable than PLD because of lower lock threshold. Also, the PLD threshold setting criteria experience cannot apply directly to FLD. These should be considered in the practical receiver implementation.

4.3. Setting Threshold with Actual Data

In order to set reasonable threshold of FLD, we can get the relationship curve about FLD output and C/N₀ based on (6) and thermal jitter of FLL [24]:

σ_{t F L L} = \frac{1}{2 π T_{c o h}} \sqrt{\frac{4 F B_{L}}{C / N_{0}} (1 + \frac{1}{T_{c o h} * C / N_{0}})}

(29)

where

σ_{t F L L}

is the frequency jitter caused by thermal noise,

B_{L}

is equivalent noise bandwidth of tracking loop,

F

is equal to 1 and 2 for low, and high C/N₀ individually. Normally,

3 σ_{t F L L}

is used to estimate tracking performance of FLL. The relationship between threshold of FLD and C/N₀ based on (6) and (29) is

T h_{C_{2 φ}} = \frac{1}{M} \sum_{i = 1}^{M} \cos (12 π T_{c o h} σ_{t F L L})

(30)

Assuming that M = 20, coherent integration length

T_{c o h} = 1 ms

B_{L} = 25 Hz

, we can figure (30) in Figure 13.

Then, we recorded actual real-time data samples from the intermediate frequency (IF) signal of GPS receiver. The data were recorded on building roof for a duration of 100 s at 62 million samples per second. The C/N₀ of the GPS satellites were 37 dB-Hz (SVN 9), 40 dB-Hz (SVN 6), 43 dB-Hz (SVN 29), and 47 dB-Hz (SVN 5). From Figure 11, we can know that the theoretical threshold of FLD is 0.2 (37 dB-Hz), 0.6 (40 dB-Hz), 0.8 (43 dB-Hz), and 0.9 (47 dB-Hz). Figure 14 shows the comparison of FLD outputs and theoretical thresholds with different satellites. The results show that the theoretical thresholds fit well with FLD outputs of real signals, which provide a solid support for FLD threshold setting criteria.

5. Conclusions

The approximate probability distribution of frequency-locked detector output has been theoretically derived by statistical approach. The mean-time-to-lose-lock (MTLL) of the frequency tracking loop has been analyzed based on the statistical characteristic of the FLD output to evaluate the performance of the tracking loop, as well as to compare with the PLD output. MTLL is a kind of evaluation criteria to measure FLD or PLD performance. The relationship among mean-time-to-lose-lock, detection threshold, and C/N₀ can be revealed with this method, which shows the difference between threshold settings in the phase-locked detector and the frequency-locked detector. The FLD threshold is, as demonstrated by the analysis results, much lower in comparison to the PLD threshold in the PLL with the same MTLL performance. The FLD output is an accumulation of adjacent periods, which means that FLD can be used to indicate a weaker signal to ensure the stability of the tracking loop. Therefore, a theoretical basis for the threshold setting in the frequency-locked detector of the GPS receiver is provided.

Meanwhile, the result of the experiment using the actual data indicates several suggested thresholds under different C/N₀ values, which can be used as a reference for the threshold setting criteria of FLD in receiver.

Acknowledgments

Sponsored by National Natural Science Foundation of China(61471017, 41374137) and Guangxi Key Laboratory of Precision Navigation Technology and Application, Guilin University of Electronic Technology (No DH201514).

Author Contributions

Tian Jin, Heliang Yuan, Na Zhao and Honglei Qin conceived and designed the experiments; Tian Jin, Heliang Yuan and Na Zhao developed the algorithms; Tian Jin, Heliang Yuan and Na Zhao coded and debugged the algorithms; Heliang Yuan performed the experiments; Tian Jin, Heliang Yuan, Kewen Sun and Yuanfa Ji analysed the data; Heliang Yuan provided field data; Tian Jin, Heliang Yuan contributed field surveying tools; all authors wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

Roudier, M.; Pena, A.J.G.; Julien, O.; Grelier, T.; Ries, L.; Poulliat, C.; Boucheret, M.-L. Demodulation Performance Assessment of New GNSS Signals in Urban Environments. In Proceedings of the 27th International Technical Meeting of The Satellite Division of the Institute of Navigation (ION GNSS+ 2014), Tampa, FL, USA, 8–12 September 2014; pp. 3411–3429. [Google Scholar]
Spilker, J.J. GPS signal structure and theoretical performance. Glob. Position. Syst. Theory Appl. 1996, 1, 57–119. [Google Scholar]
Roncagliolo, P.A.; De Blasis, C.E.; Muravchik, C.H. GPS digital tracking loops design for high dynamic launching vehicles. In Proceedings of the 2006 IEEE Ninth International Symposium on Spread Spectrum Techniques and Applications, Manaus-Amazon, Brazil, 28–31 August 2006; pp. 41–45. [Google Scholar]
Duan, R.; Liu, R.; Zhou, Y.; Song, Q.; Li, Z. A Carrier Acquisition and Tracking Algorithm for High-Dynamic Weak Signal. In Proceedings of the 26th Conference of Spacecraft TT&C Technology, Nanjing, China, 29 September 2013; Springer: Berlin/Heidelberg, Germany, 2013; pp. 211–219. [Google Scholar]
Curran, J.T. Weak Signal Digital GNSS Tracking Algorithms. Ph.D. Thesis, Department of Electrical and Electronic Engineering, National University of Ireland, Cork, Ireland, 2010. [Google Scholar]
Yang, Y.; Huang, Z. High performance digital carrier tracking loop design for high dynamic GPS receiver. In Proceedings of the 2009 5th International Conference on Wireless Communications, Networking and Mobile Computing, Beijing, China, 24–26 September 2009; pp. 1–4. [Google Scholar]
Curran, J.T.; Lachapelle, G.; Murphy, C.C. Improving the design of frequency lock loops for GNSS receivers. IEEE Trans. Aerosp. Electr. Syst. 2012, 48, 850–868. [Google Scholar] [CrossRef]
Ward, P.W. Performance comparisons between FLL, PLL and a novel FLL-assisted-PLL carrier tracking loop under RF interference conditions. In Proceedings of the 11th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GPS-98), Nashville, TN, USA, 15–18 September 1998; pp. 783–795. [Google Scholar]
Natali, F. AFC tracking algorithms. IEEE Trans. Commun. 1984, 32, 935–947. [Google Scholar] [CrossRef]
Messerschmitt, D. Frequency detectors for PLL acquisition in timing and carrier recovery. IEEE Trans. Commun. 1979, 27, 1288–1295. [Google Scholar] [CrossRef]
Mileant, A.; Hinedi, S. Lock detection in costas loops. IEEE Trans. Commun. 1992, 40, 480–483. [Google Scholar] [CrossRef]
Linn, Y.; Peleg, N. A family of self-normalizing carrier lock detectors and E S/N 0 estimators for M-PSK and other phase Modulation schemes. IEEE Trans. Wirel. Commun. 2004, 3, 1659–1668. [Google Scholar] [CrossRef]
Kratyuk, V.; Hanumolu, P.; Moon, U.K.; Mayaram, K. Frequency detector for fast frequency lock of digital PLLs. Electron. Lett. 2007, 43, 13–14. [Google Scholar] [CrossRef]
Jin, T.; Wang, Y.; Lv, W. Study of mean time to lose lock and lock detector threshold in GPS carrier tracking loops. Chin. J. Electron. 2013, 22, 46–50. [Google Scholar]
Kaplan, E.; Hegarty, C. Understanding GPS: Principles and Applications; Artech House: Norwood, MA, USA, 2005. [Google Scholar]
Bradford, P.W.; Spilker, J.; Enge, P. Global Positioning System: Theory and Applications; AIAA: Washington, DC, USA, 1996. [Google Scholar]
Misra, P.; Enge, P. Global Positioning System: Signals, Measurements and Performance, 2nd ed.; Ganga-Jamuna Press: Lincoln, MA, USA, 2006. [Google Scholar]
Curran, J.; Borio, D.; Murphy, C.C. Front-end filtering and quantisation effects on GNSS signal processing. In Proceedings of the 1st International Conference on Wireless Communication, Vehicular Technology, Information Theory and Aerospace & Electronic Systems Technology, Aalborg, Denmark, 17–20 May 2009; pp. 227–231. [Google Scholar]
Curran, J.T.; Borio, D.; Lachapelle, G.; Murphy, C.C. Reducing front-end bandwidth may improve digital GNSS receiver performance. IEEE Trans. Signal Process. 2010, 58, 2399–2404. [Google Scholar] [CrossRef]
Mongrédien, C.; Lachapelle, G.; Cannon, M.E. Testing GPS L5 acquisition and tracking algorithms using a hardware simulator. In Proceedings of the 19th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS), Fort Worth, TX, USA, 26–29 September 2006; pp. 2901–2913. [Google Scholar]
Spiegel, M.R. Mathematical Handbook of Formulas and Tables; McGraw-Hill: New York, NY, USA, 1968. [Google Scholar]
Zhuang, W. Performance analysis of GPS carrier phase observable. IEEE Trans. Aerosp. Electron. Syst. 1996, 32, 754–767. [Google Scholar] [CrossRef]
Juang, J.; Chen, Y. Phase/frequency tracking in a GNSS software receiver. IEEE J. Sel. Top. Signal Process. 2009, 3, 651–660. [Google Scholar] [CrossRef]
Razavi, A.; Gebre-Egziabher, D.; Akos, D.M. Carrier loop architectures for tracking weak GPS signals. IEEE Trans. Aerosp. Electron. Syst. 2008, 44, 697–710. [Google Scholar] [CrossRef]

Figure 1. PDF of theoretical and simulation under high C/N₀ (T = 1 ms, M = 20).

Figure 2. Probability density function (PDF) of theoretical and simulation under high C/N₀ (T = 1 ms, M = 1).

Figure 3. Frequency-locked detector (FLD) output probability density with no signal (M = 1).

Figure 4. Loop lock probability under different C/N₀ with M = 20.

Figure 5. Loop lock probability under different C/N₀ with M = 1.

Figure 6. Relation between FLD false alarm probability and threshold (Th).

Figure 7. Relation between MTLL (

E {T_{F}}

) and threshold

T h

under no signal presented situation.

Figure 7. Relation between MTLL (

E {T_{F}}

) and threshold

T h

under no signal presented situation.

Figure 8. Relation between mean-time-to-lose-lock (MTLL) (

E {T_{D}}

) and C/N₀ (M = 1) under signal presented situation.

Figure 8. Relation between mean-time-to-lose-lock (MTLL) (

E {T_{D}}

) and C/N₀ (M = 1) under signal presented situation.

Figure 9. Relation between MTLL (

E {T_{D}}

) and C/N₀ (M = 20) under signal presented situation.

Figure 9. Relation between MTLL (

E {T_{D}}

) and C/N₀ (M = 20) under signal presented situation.

Figure 10. MTLL (

E {T_{D}}

) of theoretical and simulation C/N₀ (M = 1).

Figure 10. MTLL (

E {T_{D}}

) of theoretical and simulation C/N₀ (M = 1).

Figure 11. MTLL (

E {T_{D}}

) of theoretical and simulation C/N₀ (M = 20).

Figure 11. MTLL (

E {T_{D}}

) of theoretical and simulation C/N₀ (M = 20).

Figure 12. MTLL comparison between FLD (M = 1) and PLD.

Figure 13. Relation between threshold of FLD and C/N₀ (M = 20).

Figure 14. FLD output of actual data with different C/N₀ (M = 20).

Table 1. Simulation parameters of software receiver.

Parameters	Value
Signal	GPS L1 Signal
Signal C/N₀	High C/N₀: 48 dB-Hz, 44 dB-Hz, 40 dB-Hz, 36 dB-Hz Low C/N₀: No Signal
Signal sampling rate	12 MHz
Coherent integration length	1 ms
Doppler frequency	0 Hz
Simulation times	5000 times, 5 s (M = 1 M = 20)

© 2017 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Jin, T.; Yuan, H.; Zhao, N.; Qin, H.; Sun, K.; Ji, Y. Frequency-Locked Detector Threshold Setting Criteria Based on Mean-Time-To-Lose-Lock (MTLL) for GPS Receivers. Sensors 2017, 17, 2808. https://doi.org/10.3390/s17122808

AMA Style

Jin T, Yuan H, Zhao N, Qin H, Sun K, Ji Y. Frequency-Locked Detector Threshold Setting Criteria Based on Mean-Time-To-Lose-Lock (MTLL) for GPS Receivers. Sensors. 2017; 17(12):2808. https://doi.org/10.3390/s17122808

Chicago/Turabian Style

Jin, Tian, Heliang Yuan, Na Zhao, Honglei Qin, Kewen Sun, and Yuanfa Ji. 2017. "Frequency-Locked Detector Threshold Setting Criteria Based on Mean-Time-To-Lose-Lock (MTLL) for GPS Receivers" Sensors 17, no. 12: 2808. https://doi.org/10.3390/s17122808

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Frequency-Locked Detector Threshold Setting Criteria Based on Mean-Time-To-Lose-Lock (MTLL) for GPS Receivers

Abstract

1. Introduction

2. Frequency-Locked Detector Output $C_{2 φ}$

3. Distribution of FLD Output $C_{2 φ}$

3.1. $C_{2 φ}$ Distribution under Frequency Lock

3.2. $C_{2 φ}$ Distribution under Frequency Unlock

4. Comparison between Theoretical and Simulation Results of $C_{2 φ}$

4.1. Lock Probability Analysis of FLD Output

4.2. Mean-Time-To-Lose-Lock (MTLL) of FLD

4.3. Setting Threshold with Actual Data

5. Conclusions

Acknowledgments

Author Contributions

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI

Article Menu

Frequency-Locked Detector Threshold Setting Criteria Based on Mean-Time-To-Lose-Lock (MTLL) for GPS Receivers

Abstract

1. Introduction

2. Frequency-Locked Detector Output C 2 φ

3. Distribution of FLD Output C 2 φ

3.1. C 2 φ Distribution under Frequency Lock

3.2. C 2 φ Distribution under Frequency Unlock

4. Comparison between Theoretical and Simulation Results of C 2 φ

4.1. Lock Probability Analysis of FLD Output

4.2. Mean-Time-To-Lose-Lock (MTLL) of FLD

4.3. Setting Threshold with Actual Data

5. Conclusions

Acknowledgments

Author Contributions

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI

2. Frequency-Locked Detector Output $C_{2 φ}$

3. Distribution of FLD Output $C_{2 φ}$

3.1. $C_{2 φ}$ Distribution under Frequency Lock

3.2. $C_{2 φ}$ Distribution under Frequency Unlock

4. Comparison between Theoretical and Simulation Results of $C_{2 φ}$