Research Article  Open Access
Adaptive Gain Control Method of a PhaseLocked Loop for GNSS Carrier Signal Tracking
Abstract
The global navigation satellite system (GNSS) has been widely used in both military and civil fields. This study focuses on enhancing the carrier tracking ability of the phaselocked loop (PLL) in GNSS receivers for highdynamic application. The PLL is a very popular and practical approach for tracking the GNSS carrier signal which propagates in the form of electromagnetic wave. However, a PLL with constant coefficient would be suboptimal. Adaptive loop noise bandwidth techniques proposed by previous researches can improve PLL tracking behavior to some extent. This paper presents a novel PLL with an adaptive loop gain control filter (AGCFPLL) that can provide an alternative. The mathematical model based on second and thirdorder PLL was derived. The error characteristics of the AGCFPLL were also derived and analyzed under different signal conditions, which mainly refers to the different combinations of carrier phase dynamic and signal strength. Based on error characteristic curves, the optimal loop gain control method has been achieved to minimize tracking error. Finally, the completely adaptive loop gain control algorithm was designed. Comparable test results and analysis using the new method, conventional PLL, FLLassisted PLL, and FABLL demonstrate that the AGCFPLL has stronger adaptability to high target movement dynamic.
1. Introduction
The global navigation satellite system (GNSS) is a system that uses satellites to provide autonomous positioning, which has been widely used in both military and civil fields. It allows electronic receivers that are of small volume to determine their navigation parameters (i.e., position, velocity, and time (PVT)) with high accuracy using signals transmitted by satellites anytime and anywhere around the globe. The GNSS carrier signal transmitted by satellites propagates in the form of electromagnetic wave with binary navigation data and pseudorandom noise (PRN) code modulated on it. After being transmitted by GNSS satellites, the GNSS signal would propagate through the atmosphere, subsequently being received by an antenna of the GNSS receiver near the ground and further processed by the receiver itself. So far, there exist four GNSS in the world. They are the Chinese Beidou System (BDS), the American Global Positioning System (GPS), the European Union’s Galileo Navigation Satellite System, and the Russian GLONASS.
In the process of GNSS signal propagation, the frequency of the carrier signal received by the GNSS receiver differs from the frequency of the carrier transmitted by a satellite. This phenomenon is referred to as the Doppler effect which is caused by relative motion between the GNSS satellites and GNSS receiver. Fortunately, it is the Doppler effect that provides us a method to measure the relative velocity between the receiver and satellites, which is the precondition of receiver velocity determination. That is, the technique to compute the velocity of a user is based on the received carrier Doppler frequency or phase measurements between several satellites and the receiver. Specifically, the receiver has to synchronize its own local carrier with the incoming signals to extract the navigation data and perform the Doppler measurements. Traditionally, carrier synchronization is achieved by a carrier tracking loop, which is often implemented as PLL. Apart from PLL, a frequencylocked loop (FLL) and its combination with PLL are also common carrier tracking loop structures to improve receivers’ dynamic performance [1]. Besides, vector tracking algorithms and sensoraided GNSS carrier tracking techniques are very popular in robust receiver design [2]. The carrier tracking loop is designed to track the incoming carrier’s phase dynamic or its Doppler frequency caused by the relative motion between the receiver and satellites. The Doppler frequency estimated by the carrier tracking loop mainly contains three terms: the true Doppler frequency, the receiver clock drift, and the noise [2]. The second term could be formulated after obtaining the receiver velocity. However, the noise term, which results from signal propagation and receiver processing, cannot be estimated exactly. Therefore, the carrier tracking loop is also designed to filter the noise in the incoming carrier so that more precise carrier phase or Doppler frequency measurements can be obtained, which will bring about more precise velocity measurements.
Generally, the function of PLL in GNSS receivers is to try to reduce the phase difference between the incoming carrier and local generated carrier. When PLL is in a locking state, the local carrier can be considered a precise replica of the received carrier. Inevitably, the measurements for the incoming carrier phase contain error. The dominant sources of phase error are phase jitter which is mainly caused by thermal noise and dynamic stress error which is caused by relative motion between the receiver and satellites [3]. For PLL, dynamic stress error could be relieved by increasing the PLL loop noise bandwidth whereas phase jitter error can be reduced by decreasing the PLL loop noise bandwidth. That means that improving receiver noise performance would degrade its dynamic performance. Hence, the confliction between PLL dynamic performance and noise performance necessitates some compromises in PLL design, and many methods have been developed to address this typical issue [4–8].
A method to determine appropriate loop noise bandwidth is considering the worst carriertonoise power density ratio (CNR) and max relative dynamic (i.e., max relative acceleration for the 2ndorder PLL and max relative jerk for the 3rdorder PLL) [4]. Unfortunately, it will lead to a suboptimal loop for low relative dynamic. A PLL with two sets of loop coefficients working for pullin mode and tracking mode, respectively, can decrease pullin time [5]. However, it is hard to improve carrier tracking capability. Adapting loop noise bandwidth depending on the realtime estimations of the relative dynamics and CNR, known as fast adaptive bandwidth lock loops (FABLL), is a realtime optimal PLL through minimizing the thermal noise within the total tracking error [6, 9]. The continuation of FABLL has been also proposed in [10]. The second version of FABLL is able to automatically and dynamically set loop noise filter parameters of PLL and DLL in order to minimize the total tracking error. However, the noisy and timevarying CNR and relative dynamic estimations which might cause intense and frequent loop noise bandwidth adjustment limit its application. The PLL does not allow its coefficient to change quickly; otherwise, the loop filter will output impulses resulting in PLL losing lock [6]. To solve this problem, a simple method is to insert a lowpass filter to delay the bandwidth change rate [3].
Besides, PLL behavior can be affected significantly by a loop gain, which is not discussed in detail in previous research. Researches have demonstrated that a reasonable loop gain could promote carrier tracking performance [11, 12]. By contrast, an adaptive loop gain control method assisted by FLL has been designed for fast convergence from pullin to tracking mode and decreasing tracking error when PLL has to deal with highdynamic movement [7]. This method establishes an empirical relation between the loop gain and frequency difference between the received carrier and local carrier. The loop gain could be adaptively adjusted according to the frequency discriminator outputs. However, the use of empirical relation without a mathematical model accurately modeling and analyzing the system leads to a suboptimal PLL. The performance of this kind of method has not been compared with that of FABLL. A mathematical model describing the 2ndorder PLL with a loop gain control factor has been derived in [8]. It has been shown that a loop gain can control loop noise bandwidth and characteristic frequency. According to its results, an optimal loop gain can be calculated to minimize the dominating tracking error in different signal conditions. It also shows the possibility that the loop gain could be adjusted adaptively based on CNR and relative dynamic estimations. Furthermore, the 3rdorder PLL can track the phase acceleration input and has advantage over the 2ndorder PLL in tracking carrier phase dynamic. From the comprehensive theoretical research perspective, the modeled relationship should also be analyzed.
The focus of this paper is the development and analysis of an adaptive loop gain control technique for 2nd and 3rdorder PLL to enhance PLL highdynamic adaptability. The structure of the paper is shown as follows: In Section 2, the mathematical model of AGCFPLL is derived. A loop gain control factor is inserted into the traditional PLL to adjust the loop gain conveniently. In order to avoid intense adjustment of a loop gain, an inertial element is added to make the loop gain adjustment smoother. In Section 3, the relationship between the behavior of AGCFPLL and its loop gain is analyzed in detail in time and frequency domains, respectively. In Section 4, an optimal loop gain adjustment strategy is designed to minimize the dominating tracking error according to estimated CNR and relative dynamic between the receiver and satellites. A simple averaging operation can be used to reduce the noise in the relative acceleration estimations and relative jerk estimations. Comparable tests and analysis are conducted between the novel adaptive loop gain control PLL, the conventional PLL, the FLLassisted PLL, and the FABLL, to show the loop dynamic adaptability of the novel algorithms in Section 5. Finally, some concluding remarks are given in Section 6.
2. GNSS Carrier Tracking Loop with an Adaptive Loop Gain Control Filter (AGCFPLL)
This section reviews the basic concepts of the traditional GNSS carrier tracking loop. The idea to incorporate the loop gain control filter to change the loop gain conveniently and smoothly will build upon the traditional loop structure, and the mathematical model of this new PLL is derived. Then, the digital implementation of this kind of PLL will be discussed. The traditional PLL linear model was discussed in more detail in [3, 5, 12–16].
2.1. Traditional Linear Model of PLL
The function of PLL within GNSS receivers is to try to generate a local sinusoidal carrier whose phase is the same with the received carrier phase. When PLL is in a locking state, the received carrier phase can be measured by observing the local generated carrier. Meantime, the carrier is wiped off from the received signal. Figure 1 shows the traditional linear PLL model. There are three components: phase discriminator (PD), loop filter (LF), and numerically controlled oscillator (NCO). The phase discriminator calculates the phase difference between the received carrier and local carrier. The phase difference is sent into the loop filter which not only filters the noise but also generates the control signals to NCO. NCO transforms the control signals into local carrier phase adjustment so that the local carrier phase is precisely aligned with the received carrier phase.
The PLL transfer function and its error transfer function can be derived from Figure 1 as where is the transfer function of the loop filter and is the loop gain that is defined as
1storder LF and 2ndorder LF are often utilized to track certain phase dynamic in a received signal [17]. The 2ndorder PLL with a 1storder loop filter is able to track phase acceleration input without steadystate error, which is caused by velocity in relative movement. The 3rdorder PLL with a 2ndorder loop filter is able to track phase jerk input caused by acceleration in relative movement without steadystate error. The transfer functions of 2ndorder and 3rdorder loops are expressed as follows: where is the damping factor of the 2ndorder loop, and are the loop coefficients of the 3rdorder loop, and and are the natural frequencies of the 2ndorder and 3rdorder PLL, respectively.
Substituting (3) into (1), one can rewrite the 2ndorder PLL system transfer function and error transfer function as
Similarly, system transfer function and error transfer function for the 3rdorder PLL are given by
2.2. Model of AGCFPLL
In order to control the loop gain conveniently and smoothly according to different relative dynamics and signal strengths, a loop gain control filter can be inserted into the traditional PLL structure, including a gain control factor () and inertial element () as shown in Figure 2. Consequently, the actual loop gain could be expressed as
The loop gain could then be adjusted by . Generally, and are set to one when arc tangent phase discriminator is employed. Thus, is the actual loop gain. The transfer function of the loop gain control filter is where is the time constant of the inertial element in the unit of seconds and is the optimal loop gain which can minimize the dominating tracking error. The inertial element can delay the change rate of the actual loop gain, as shown in Figure 3. It takes about seconds before the actual loop gain reaches its steady value when the optimal loop gain changes.
Figure 4 shows the distribution rejection ability of the inertial element in loop coefficient adjustment. Here, we create a scene where there are some step changes in loop gain adjustment. There are two options in loop gain adjustment. The loop gain is adjusted directly in option 1, while the loop gain passed the inertial element in option 2. As shown in the above figure of Figure 4. The ability of the inertial element rests in the fact that it would delay the change rate of the loop gain. The below figure in Figure 4 shows the carrier phase error discriminator output in two options. The discriminator output represents the carrier phase tracking error and reflects the locking condition of PLL for the carrier phase. It is obvious that the direct loop gain adjustment will introduce disturbance to the tracking loop. It takes some time before the PLL recovers to a locking state from this disturbance. As for option 2, the inertial element adjusts the loop gain in a smooth way so that the disturbance is eliminated. Another function of the inertial element lies in decreasing the noisy optimal K_{c} calculation, which will be explained in detail in Section 4.
According to Figure 2, the 2nd and 3rdorder AGCFPLL transfer function and error transfer function can be derived to be their general forms as follows:
The inertial element is excluded when deriving these equations because of its delay function. It should be noted that actual would follow optimal . Furthermore, the inertial element can be ignored when reaches its steady value.
To normalize the AGCFPLL models for the sake of finding the relationship of loop coefficients between the ordinary PLL and AGCFPLL, two sets of new loop filter coefficients are introduced here. The first set of new coefficients is used to normalize the 2ndorder AGCFPLL and defined as follows:
Substituting (12) and (13) into (8) and (9), one can rewrite the 2ndorder AGCFPLL transfer function and error transfer function to be their normalized forms as
Obviously, these two equations are similar to the model of the ordinary 2ndorder PLL, as expressed in (4).
The second set of coefficients is used to normalize the 3rdorder AGCFPLL and stated as follows:
Substituting (15) and (16) into (10) and (11), one can rewrite the 3rdorder AGCFPLL transfer function and error transfer function to be their normalized forms as
It is obvious that these two equations are similar to the model of the ordinary 3rdorder PLL, as expressed in (5).
By comparing the ordinary PLL model with the normalized AGCFPLL model, it can be found that they have similar form with different loop coefficients. The change of in AGCFPLL is equivalent to change loop coefficients in the ordinary PLL. This relation is shown in (12), (13), (15), (16), and (17). It is easy to find out that the AGCFPLL is equivalent to the ordinary PLL if the loop gain is set to one.
2.3. Digital Implement of AGCFPLL
In GNSS receivers, the PLL is implemented in the form of a Costas loop which is insensitive to phase transitions due to navigation bits [15]. Figure 5 shows the digital implementation of AGCFPLL including components like mixers and integration and dump operations.
The digital intermediate frequency (IF) signal is first mixed with the local generated carrier and then correlated with the local pseudorandom noise (PRN) code. The correlation results are used in a phase discriminator to obtain the phase error that contains noise. Subsequently, the LF filters noisy phase error and generates a NCO control signal to synchronize the local carrier with the received signals. Relative dynamic and signal C/N_{0} estimations are used to calculate the optimal loop gain in order to minimize the dominating tracking error.
The 1st and 2ndorder loop filter and inertial element are all discredited to their digital forms by bilinear transformation, which are derived as follows: where is the coherent integration time in the unit of seconds. Arc tangent function was chosen as the phase discriminator in this paper.
3. Time Domain and Frequency Domain Analyses of AGCFPLL
In this section, AGCFPLL’s behavior is analyzed in the time domain and frequency domain using the models derived in Section 2. The relationship between the AGCFPLL’s behavior and will be discussed in detail without considering the inertial element. The reason has been explained in Section 2. An averaging method of estimating relative acceleration and jerk will be derived based on time domain analysis. In the following discussion, the initial values of loop coefficients of the 2nd and 3rdorder AGCFPLL are chosen to be their common values as follows:
These loop coefficients can be adjusted by changing , as discussed in Section 2.
3.1. Time Domain Analysis
The primary function of PLL is to make the local replicated carrier phase track the received one [16]. Therefore, the PLL could be regarded as a control system, or more specifically, a servo control system. Furthermore, the PLL system is a nonlinear feedback control system in the classical control theory point of views, where the loop filter is the controller and NCO is the plant. Its nonlinearity results from the nonlinear phase discriminator. Fortunately, the PLL’s behavior in a locking state can be predicted by its linearized model, as derived in Section 2. Therefore, some classical analysis methods of the control system can be employed to analyze the AGCFPLL. According to (8) and (10), the characteristic equation of the second and thirdorder AGCFPLL can be found as
Figures 6 and 7 show the root locus of the 2nd and 3rdorder AGCFPLL for different . It can be seen that with an increase in , the closedloop poles move towards minus infinity or openloop zero points from openloop poles. With the help of the classical control theory, it can be concluded that the response of the AGCFPLL would be fast with less sinusoidal oscillation and less overshoot with increment. It is worth noting that the 2ndorder AGCFPLL is unconditionally stable while the 3rdorder AGCFPLL is not. The stable condition of the thirdorder AGCFPLL is
Figures 8 and 9 show the step response of the AGCFPLL for different values of the loop gain. The results of step response simulation are consistent with those of the previous analysis and more visualized.
The PLL tracking error caused by relative dynamics could be written as [14] where is the order of PLL, is the tracking error in the unit of meters, and is the carrier wavelength in the unit of meters. According to (13), (15), and (24), increasing can decrease dynamic stress error. In conclusion, large brings about high convergence speed and less dynamic stress error. Thus, large is expected in the time domain without regard to the noise in the input carrier.
3.2. Frequency Domain Analysis
The PLL tracks the noisy received carrier phase so that the local replicated carrier contains noise unavoidably. Therefore, the secondary function of PLL is to filter the noise in the received carrier phase in order to obtain an accurate carrier replica. As one of the important parameters of the PLL, noise bandwidth controls the amount of noise allowed in the PLL. This parameter is defined as follows [15]:
The relationship between the noise bandwidths and loop coefficients for the 2nd and 3rdorder PLL can be expressed as follows: where and are the loop noise bandwidths of 2nd and 3rdorder PLL, respectively.
The relationship between the noise bandwidth and of the 2ndorder AFCFPLL can be derived by inserting (12) and (13) into (26). Similarly, this relationship for the 3rdorder AGCFPLL can be derived by inserting (15), (16), and (17) into (27). Finally, those two relationships are expressed as follows:
Here, and are the loop noise bandwidths of the 2nd and 3rdorder AGCFPLL, respectively.
Figures 10 and 11 show the relationship between the loop noise bandwidth and for the AGCFPLL in graphic way. It can be concluded that the loop noise bandwidth of the 2nd order AGCFPLL is proportional to the while the loop noise bandwidth of the 3rd order AGCFPLL is positively correlated with the .
3.3. The Estimation of Relative Dynamic and C/N_{0}
As the first step to calculate optimal , relative acceleration and jerk estimations play important roles in AGCFPLL. Actually, (24) provides us a method to calculate relative acceleration or jerk according to phase error. As long as the AGCFPLL is in a locking state, the relative acceleration can be obtained from the 2ndorder AGCFPLL and the relative jerk can be obtained from the 3rdorder AGCFPLL. This relation is shown as follows: where is the nth phase error value estimated by the 2ndorder AGCFPLL and is the nth phase error value estimated by the 3rdorder AGCFPLL. However, the estimated phase error contains much noise, which results in noisy acceleration and jerk estimations. Therefore, a simple averaging operation is employed to filter the noise. Thus, (29) can be rewritten as follows: where M is the number of phase error values from a phase discriminator. M should not be too small considering the noise influence. In this paper, N is set as 500.
C/N_{0} estimations are another coefficient that the AGCFPLL algorithms need to calculate the optimal . The correlator comparison method is used to measure C/N_{0} in this paper [18]. The measurement, Z, is first formed as follows: where and are the accumulated outputs of the prompt correlator, is the unsynchronized accumulated correlator output, and denotes the number of iterations averaged over to smooth noise. Then, C/N_{0} estimations could be calculated employing the following equation: where is in the unit of dBHz while is in the unit of Hz.
4. Optimized Loop Gain of AGCFPLL
In this section, the AGCFPLL measurement error will be derived and discussed in different conditions first. Then, the optimal is calculated according to the estimated and relative dynamic to obtain the minimum tracking error. Finally, the loop gain control method will be proposed.
4.1. Error Characteristic of AGCFPLL
The dominant sources of the phase error in the PLL are phase jitter and dynamic stress error. The phase jitter contains thermal noise and oscillator noise (i.e., vibrationinduced oscillator jitter and Allan varianceinduced oscillator jitter) [15]. The 1σ tracking error of the PLL is given by where R is the lineofsight (LOS) distance between the satellite and receiver in the unit of meters. is chosen to be 1 millisecond in this paper. Vibrationinduced oscillator jitter (approximately ) and Allen varianceinduced oscillator jitter (approximately ) are also considered but not expressed in the equation for the sake of simplicity.
The 1σ tracking error of the 2ndorder AGCFPLL could be derived by inserting (12) and (13) into (33):
Similarly, the 1σ tracking error of the 3rdorder AGCFPLL could be derived by inserting (15), (16), and (17) into (33):
The error characteristic with an acceleration of 4 and different C/N_{0} is shown in Figure 12, where is the acceleration of the gravity. It can be observed that the 1σ tracking error decreases to minimum, then increases to positive infinity with the increase in . So there exists minimum, maximum, and optimum . The maximal and minimal values were located at the crossing point of the error characteristic curve and tracking threshold curve (). The optimum is located on the point that has the minimal phase error. It could also be observed that the high would extend the range of C/N_{0}. The error characteristic of the 3rdorder AGCFPLL is similar to that of the 2ndorder AGCFPLL and sensitive to relative jerk, which is not shown here.
4.2. Optimized Loop Gain of AGCFPLL
According to the analysis above, the optimal value of could be calculated by (36) while the maximum and minimum values of could be calculated by (37). Because high would degrade the loop performance when the transformation (bilinear, boxcar transform) method is used to implement the loop filter [19, 20], another constraint () is added.
Figures 13 and 14 show the maximum, minimum, and optimal of the 2nd and 3rdorder AGCFPLL in different conditions. Figures 15 and 16 show their minimal tracking errors in different conditions when optimal is employed.
4.3. Adaptive Loop Gain Control Strategy
Based on the above analysis, the maximal K_{c} could be used initially to make the tracking loop converge from the pullin state to the locking state. The reason has been discussed in Section 3. After the loop goes into the locking state, K_{c} can be adjusted adaptively according to the C/N_{0} and relative dynamic estimations. Considering the complexity to solve (36) and (37) that have a radical sign, a lookup table is used to store the relationship among C/N_{0} estimations, relative dynamic estimations, and optimal K_{c} in a discrete form. The C/N_{0} and relative dynamic estimations discussed in Section 3 contain accumulation that could be considered to be the first lowpass filtering. The estimated C/N_{0} and relative dynamic will be the input of the lookup table. The twodimensional linear interpolation method is employed to calculate the optimal K_{c} subsequently. However, the estimated C/N_{0} and relative dynamic still contain noise. Thus, the optimal K_{c} is noisy too. The inertial element introduced here would further filter the noise in optimal K_{c}. Therefore, the adjustment frequency decreases. A part of the lookup table with specific C/N_{0} and dynamic for AGCFPLL is shown in Tables 1 and 2.


5. Test, Results, and Analysis
The algorithms discussed in this paper were implemented in the software GPS receiver and tested using simulated GPS signals from a navigation satellite signal simulator. Because the 2ndorder AGCFPLL is sensitive to relative acceleration while the 3rdorder AGCFPLL is sensitive to relative jerk, two scenarios were designed to demonstrate the tracking performance under highdynamic application—highacceleration scenario and highjerk scenario. The AGCFPLL is compared to their constant loop gain versions (conventional PLL), FLLassisted PLL, and FABLL, to show the improvements.
5.1. Experimental Setup
In Figure 17, the navigation satellite signal simulator (HWARNSS7300 manufactured by HWA Create Corporation Ltd.) controlled by a laptop can provide L1 C/A code signals under various scenarios. The complicated straight motion scenario is employed to generate highdynamic motion. The power of the carrier signal can be adjusted by the laptop control signal, which can be used to change the C/N_{0} of C/A code signals. A lownoise amplifier (LNA; HLLNA012102B) is used to play the same role as the LNA in a real GPS antenna. Connected with LNA is the GPS IF signal collector (NewStar210M made by OlinkStar Corporation) which can transfer digital IF signals to the laptop. At last, the postprocessing including acquisition and tracking for the IF signal is done by software.
5.2. Metrics
The metrics used to assess the algorithm performance is the cosine value of the phase difference (), which is given by
When PLL is locked, the phase error should keep close to 0 and therefore should be approximately equal to one. Moreover, deviating from one indicates that PLL performance degrades, while the irregular alterations of mean that the loop is losing lock. It is safe to assume that if the value of the lock indicator is greater than 0.8, the PLL is locked [4]. Although this value is corrupted with noise, it is this value that indicated that the loop locking status and its mean value can be used to assessing carrier tracking performance of the loop.
5.3. Test Results and Analysis for the 2ndOrder AGCFPLL
The highacceleration scenario (test scenario A) is simulated and used to test the 2ndorder AGCFPLL algorithms, and its performance is compared with the traditional 2ndorder ordinary PLL, 1storder FLLassisted 2ndorder PLL, and FABLL. At the beginning, the vehicle stays for 10 seconds and then a highacceleration motion is simulated. The bandwidth of the ordinary 2ndorder PLL is set to 20 Hz. The bandwidth of FLL and PLL for FLLassisted PLL is set to 5 Hz and 20 Hz, respectively. The test GPS signal came from satellite 21.
In order to adjust the K_{c}, the C/N_{0} and relative acceleration are estimated firstly, as shown in Figure 18. Compared with estimations given by [4], the averaging method makes the estimated acceleration less noisy, which will decrease the regulating frequency of K_{c}. According to the estimated signal conditions, the optimal K_{c} can be found through the lookup table and twodimensional interpolation, as discussed in Section 4. Figure 19 shows the optimal and actual K_{c}. We can find that the actual K_{c} would track the optimal K_{c} in a slow and smooth way, avoiding PLL losing lock which is caused by sharp K_{c} adjustment.
Figure 20 shows the arc tangent discriminator outputs of the four kinds of PLL. The arc tangent discriminator outputs represent the carrier phase tracking error. It can be found that the ordinary 2ndorder PLL cannot bear the highacceleration motion and lost lock for the carrier phase eventually. The 1storder FLLassisted 2ndorder PLL is able to track the carrier phase dynamic. However, there is a steadystate error. The carrier phase tracking error is very close to the effective working range (−0.25 cycle, 0.25 cycle) of the arc tangent discriminator. Thus, this PLL is easy to lose lock, especially when relative acceleration or noise increases. As for the FABLL and AGCFPLL, they can keep the carrier phase lock in the highacceleration scenario with less carrier phase tracking error. The AGCFPLL has comparable performance with FABLL according to this figure. In order to quantitatively evaluate their performance, further evaluation is done and shown as follows.
Figure 21 shows the PLL lock status of the four kinds of PLL in test scenario A, while Table 3 presents the mean value of the lock indicator. Because the ordinary PLL loses lock for the carrier phase in test scenario A, the mean value of the lock indicator is invalid and not shown in this table. The mean value of the lock indicator of the AGCFLL during highacceleration motion (10~40 s) is about 0.9883, which is closer to 1 than that of the FABLL. This suggests that the 2ndorder AGCFPLL outperforms the other kinds of 2ndorder PLL in highacceleration scenarios.

5.4. Test Results and Analysis for the 3rdOrder AGCFPLL
Similarly, a highjerk scenario (test scenario B) is simulated and employed to test the performance of the ordinary 3rdorder PLL, 2ndorder FLLassisted 3rdorder PLL, 3rdorder FABLL, and 3rdorder AGCFPLL. The bandwidth of the ordinary 3rdorder PLL is set to 18 Hz. The bandwidth of FLL and PLL for FLLassisted PLL is set to 5 Hz and 18 Hz, respectively. The test GPS signal came from satellite 31.
Figure 22 shows the estimated C/N_{0} and relative jerk of AGCFPLL in test scenario B. These two data determine the optimal K_{c}, which is able to minimize the carrier phase tracking error. The inertial element is employed to make the actual K_{c} track the optimal K_{c} in a smooth way, thus decreasing the effect of noisy C/N_{0} estimation and relative jerk estimation. This relation is shown in Figure 23.
Figure 24 shows the arc tangent discriminator outputs of the four kinds of PLL in test scenario B. The four kinds of PLL are all able to track the carrier phase dynamic. However, the carrier phase tracking errors of the ordinary 3rdorder PLL and 2rdorder FLLassisted 3rdorder PLL are bigger than those of AGCFPLL and FABLL. The FABLL and AGCFPLL can keep the carrier phase lock in the highjerk scenario with less carrier phase tracking error. The AGCFPLL has comparable performance with FABLL according to this figure. In order to quantitatively evaluate their performance, we did further evaluation and the results are shown as the following parts.
Figure 25 shows the PLL lock status of the four kinds of PLL in test scenario B, while Table 4 presents the mean value of the lock indicator. The mean value of the lock indicator of 3rdorder AGCFPLL during highjerk motion (18~40 s) is about 0.9916, which is closest to 1, than that of other kinds of 3rdorder PLL. This suggests that 3rdorder AGCFPLL has the least carrier phase tracking error in these four kinds of PLL.

6. Conclusions
This paper presents a novel adaptive loop gain control approach for second and thirdorder PLL in GNSS receivers, and its performance was tested and verified compared with that of the traditional PLL. The new system model, time and frequency domain analyses, and error characteristic of the new PLL are derived and displayed in detail, which could be regarded as the first contribution of this paper. Although only GPS L1 C/A code signals are tested, it is easy to generalize this method to other GNSS signals. The loop gain can be adjusted to its optimal value according to C/N_{0} and relative dynamic estimations to obtain the minimal tracking error. Compared with the conventional PLL, FLLassisted PLL, FABLL, the proposed adaptive loop gain algorithms have a stronger tracking ability to deal with highdynamic stress.
An inertial element is inserted in the new adaptive PLL in order to decrease the change rate of the loop gain, whose intense adjustment may cause a disturbance to PLL or results in PLL losing lock for the carrier phase.
The third contribution of this paper is the lookup table which stores the relationships among the signal conditions and optimal loop gain. Using the method of the lookup table to calculate the loop gain decreases the computation cost depending on calculating equation with radicals so that a realtime AGCFPLL can be realized.
The predetection integration time and loop order are another two coefficients of PLL that can be made adaptive. These two coefficients can also be adjusted according to signal condition to increase the tracking performance of PLL. But the predetection integration time adjustment and loop order adjustment are not discussed in this paper and would become the next research orientation.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work is funded by the National Natural Science Foundation of China under Grant nos. 61633008, 61304234, 61773132, and 61374007.
References
 P. W. Ward, Performance Comparisons between FLL, PLL and a Novel FLLAssistedPLL Carrier Tracking Loop under RF Interference Conditions, Proceedings of the 11th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GPS 1998), Nashville, USA, 1998.
 S. Alban, D. M. Akos, S. M. Rock, and D. GebreEgziabher, Performance Analysis and Architectures for INSAided GPS Tracking Loops, Institute of Navigation–ITM (2003), San Diego, California, USA, 2003.
 P. Misra and P. Enge, Global Positioning System: Signals, Measurements and Performance, GangaJamuna Press, Lincoln, MA, Second edition, 2011.
 P. Lian, Improving tracking performance of PLL in high dynamic applications, [Ph.D. thesis], Department of Geomatics Engineering University of Calgary, Calgary, Canada, 2004.
 P. Rinder, Design of a single frequency GPS software receiver, [Ph.D. thesis], Institute of Electronic Systems, Aalborg University, Aalborg, Denmark, 2004.
 F. Legrand, C. Macabiau, J. L. Issler, L. Lestarquit, and C. Mehlen, Improvement of Pseudorange Measurements Accuracy by Using Fast Adaptive Bandwidth Lock Loops, ION GPS 2000 13th International Technical Meeting of the Satellite Division of the Institute of Navigation, Salt Lake City, USA, 2000.
 J. B. Zhang, S. F. Zhang, Q. Hu, Y. Jiang, and X. Sun, Design of GNSS Carrier Tracking Loop Base on Adaptive Gain Control, Cross Strait QuadReginal Radio Science and Wireless Technology Conference, 2011, Harbin, China, 2011.
 J. C. Ding, Z. B. Luo, and L. Zhao, GNSS Carrier Tracking Loop with Loop Gain Control Factor, International Conference about Information and System, 2015, Harbin, China, 2015.
 F. Legrand and C. Macabiau, Results of the Implementation of the Fast Adaptive Bandwidth Lock Loops on a Real GPS Receiver in a High Dynamics Context, GNSS 2001, International Conference on Global Navigation Satellite System, Sevilla, Spain, 2001.
 F. Legrand and C. Macabiau, “Realtime minimization of the total tracking error in phase and delay lock loopsa second approach of the fast adaptive bandwidth algorithm,” in Proceedings of the 57th Annual Meeting of the Institute of Navigation, Albuquerque, NM, USA, 2001. View at: Google Scholar
 K. Borre, D. M. Akos, N. Bertelsen, P. Rinder, and S. H. Jensen, A SoftwareDefined GPS and Galileo Receiver: A SingleFrequency Approach, Springer Science & Business Media, 2007.
 J. B.Y. Tsui, Fundamentals of Global Positioning System Receivers: A Software Approach, Jone Wiley & Sons, 2000. View at: Publisher Site
 J.H. Won, D. Dötterböck, and B. Eissfeller, “Performance comparison of different forms of Kalman filter approaches for a vectorbased GNSS signal tracking loop,” Navigation, vol. 57, no. 3, pp. 185–199, 2010. View at: Publisher Site  Google Scholar
 G. Xie, Principles of GPS and Receiver Design, Publishing House of Electronics Industry, Beijing, 2009.
 E. D. Kaplan and C. Hegarty, Understanding GPS: Principles and Applications, Artech house, 2005.
 F. M. Gardner, Phase Lock Techniques, John Wiley & Sons, 2005. View at: Publisher Site
 C. R. Hamm, W. S. Flenniken, D. M. Bevly, and D. Lawrence, “Comparative performance analysis of aided carrier tracking loop algorithms in high noise/high dynamic environments,” in Proceedings of the 2004 ION GNSS Conference, Long Beach, California, USA, 2004. View at: Google Scholar
 P. D. Groves, “GPS signaltonoise measurement in weak signal and highinterference environments,” Navigation, vol. 52, no. 2, pp. 83–94, 2005. View at: Publisher Site  Google Scholar
 P. L. Kazemi, “Optimum digital filters for GNSS tracking loops,” in Proceedings of ION GNSS, 2008, Savannah, GA, USA, 2008. View at: Google Scholar
 P. A. R. Roncagliolo, C. E. De Blasis, and C. H. Muravchik, “GPS digital tracking loops design for high dynamic launching vehicles,” in 2006 IEEE Ninth International Symposium on Spread Spectrum Techniques and Applications, ManausAmazon, Brazil, August 2006. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Zhibin Luo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.