Research Article  Open Access
Enhancing WeakSignal Carrier Phase Tracking in GNSS Receivers
Abstract
Examining the performance of the GNSS PLL, this paper presents novel results describing the statistical properties of four popular phase estimators under both strong and weaksignal conditions when subject to thermal noise, deterministic dynamics, and typical pedestrian motion. Design routines are developed which employ these results to enhance weaksignal performance of the PLL in terms of transient response, steadystate errors, and cycleslips. By examining both single and datapilot signals, it is shown that appropriate design and tuning of the PLL can significantly enhance tracking performance, in particular when used for pedestrian applications.
1. Introduction
Despite the military origins of Global Navigation Satellite Systems (GNSS), the most widespread use of GNSS receivers is civilian and the single most common receiver platform is the cellular handset. Although the civilian user is, generally, less demanding in terms of position, velocity, and timing accuracy, signal processing for civilian applications is not a simple task. Severe attenuation experienced in the indoor environment, multipath propagation through urban environments, and the limitations of consumergrade receivers are all obstacles to maintaining acceptable receiver performance.
While many receivers can adequately track carrier frequency under most operating conditions, including in the indoor environment, reliable carrier phase tracking still proves challenging. Owing to a very short wavelength, when subject to any appreciable attenuation, the dynamics of pedestrian motion can induce carrier phase cycleslips or even loss of phaselock. Despite these challenges, the ability to track carrier phase is desirable for many reasons including enhanced bitsynchronization, reduced biterrorrate, enhanced range estimation, improved velocity estimation, and, ultimately, provision for carrierbased positioning.
In response to this challenge, this paper focuses on the process of carrier phase tracking in a scalar phaselockloop (PLL). The primary weakness of the PLL when operating on attenuated signals is the process of phase error estimation or phase discrimination. The performance of phase discriminator functions typically degrades rapidly with reduced signal strength and their behavior under weaksignal conditions is generally unique to each discriminator function. To best design a PLL, therefore, this behavior must be understood. This work aims to develop a thorough mathematical model for the carrier phase discriminator and, from this model, to infer best practices for GNSS PLL design. In particular, two case studies are investigated: pedestrian navigation using the GPS L1 C/A signal and datapilot tracking of the Galileo E1 B/C signal.
Two classes of phase discriminator will be examined, those which employ purePLL discriminators and those which employ Costas discriminators. PurePLL discriminators are those which are designed to capture the entire phase error on the interval and are therefore useful for synchronization with continuous wave signals or those with smooth modulation, such as frequencymodulation. They represent the earliest form of PLL, dating back to the 1930s [1] and over the last decade have seen applications in GNSS receivers for modernized signals which include a pilot signalcomponent. By the 1950s, the use of suppressedcarrier modulation required the development of PLLs which were insensitive to carriermodulation of which the most notable is the Costas PLL [2]. This type of PLL, capturing the phase error on the interval , is widely used in GNSS receivers for BPSK modulated signals, such as GPS L1 C/A or Galileo E1B. Strictly speaking, the Costas PLL is that which performs phase estimation via the product of the inphase and quadrature baseband channels; however, the term Costas PLL or Costas discriminator has become synonymous with the class of all modulationinsensitive phase discriminators.
The paper is organized as follows: Section 2 introduces the GNSS signal, the PLL architecture, and the linearized PLL model. A statistical analysis of four popular carrier phase estimators is developed in Section 3. Weaksignal effects on the transient and steadystate performance of the PLL are considered in Section 4 and Section 5 presents the application of the theory developed here to the problem of PLL design.
2. Receiver Model and PLL Architecture
To facilitate the following analysis, the PLL is modeled as a simplified linear, timeinvariant (LTI) system. A model of the received signal and the corresponding correlator values are developed and a general description of the classical PLL is introduced. A selection of discriminator functions are examined and equivalent linear models are provided, including an assessment of the operating region over which the linearization is accurate. These component models are then combined to yield a linearized system describing the PLL operation. Through these models, it is proposed that the PLL behavior under weaksignal conditions can be described as the superposition of the response of an equivalent linear model of the PLL to various stimuli, including that of thermal noise and of phase variations, where the particular linear model is a function of the prevailing signal strength.
2.1. Downconversion and IF Signal Processing
The correlation of the local replica signals with the incoming digital intermediate frequency (IF) signal over the interval can be approximated by the wellknown expressions for the inphase, , and quadrature, , values [3, 4]:where , , and denote the mean code phase, carrier frequency, and carrier phase errors, respectively, and is the spreading code autocorrelation function. The variable denotes the coherent integration period and also defines the interval between successive updates of the tracking loop. It is assumed that the coherent integration period is aligned with the data modulation symbol boundaries, such that the variable denotes the data sign, which is constant during correlation interval. Under normal PLL operation, the code phase and carrier frequency are reasonably well tracked by the receiver, such that , and so they have a negligible effect on (1). The propagation of the thermal noise to the correlator values is modeled as additive white Gaussian noise (AWGN):where represents the onesided thermal noise floor in W/Hz. An estimate of the carrier phase tracking error, , is then made by applying a carrier phase discriminator to the values and . This estimation procedure is discussed in more detail in Section 3.
2.2. The PhaseLock Loop
The standard phaselock loop is a feedback control loop which tracks the carrier phase using estimates of the carrier phase tracking error. Although all realizable phase error estimators are nonlinear, if the estimate is linearized around zero phase error and normalized such that the noisefree estimate has unity gain, the phase error estimate, denoted by , can be approximated by [4]where represents the linear region of the discriminator. That is, the phase error estimate is approximately equal to a constant times the true phase error, plus a zero mean, white noise, . The constant gain, , is referred to as the discriminator gain and depends on the chosen discriminator function and the prevailing signaltonoiseratio. The variance of is also dependent on the phase discriminator used and the received signaltonoise ratio. The twosided spectral density of is denoted here by . The linear region is defined as the interval , over which this discriminator approximation is valid. The exact details of the linearization of this phase error estimate and the values of the PSD of for various discriminators will be given in Section 3.
The remainder of the PLL is linear and can be represented by a system of domain transfer functions, where the update interval of the system is . Such a linearized loop model is useful as it facilitates the estimation of loop stability and tracking performance. Of particular interest are the transfer functions between the carrier phase, , and the carrier phase estimate, , between the carrier phase, , and the tracking error, , and between the thermal noise, , and the tracking error, . These quantities are depicted in a linearized loop model in Figure 1. The transfer functions of interest are given bywhere uppercase symbols represent the transform of the corresponding lowercase time series. The functions and represent the transform of the loop filter and the numerically controlled oscillator, respectively. The numerically controlled oscillator is defined as [4]
Generally, a proportional and integral (PI) controller is used in GNSS PLLs. A generalization of this type of controller takes the form [5]where is the order of the resultant closed loop system.
3. Carrier Phase Estimation
As discussed in Section 2.2, the performance of the PLL in the presence of AWGN can be estimated by examining the linear model and the noise performance of the carrier phase discriminator. Four popular carrier phase discriminators are examined here: the fourquadrant arctangent discriminator, the arctangent discriminator, the decisiondirected discriminator, and the simple quadrature discriminator. The discriminators are characterized in terms of gain, , and variance, , whereand, in Section 3.7, the linear region will also be considered.
In the carrier phase discriminator analysis that follows, it is assumed that the PLL is operating normally, with a mean frequency error of zero and a maximum frequency error that is reasonably small relative to the update interval, such that the phase error accrued over the update interval is less than the linear region of the discriminator, for example.
3.1. Measuring Signal Quality
It will be shown that the performance of the estimators considered in this work varies with the signaltonoise ratio (SNR) of the correlator values, specifically the coherent SNR, defined as [3, 6]
This metric represents the quality of the signals which are applied to the discriminator and is largely the same as the metric used in, for example, [7], for characterizing baseband communication systems with the distinction that need not, necessarily, correspond to a full bit period. A number of factors influence the value of the including the received signal power, the receiver’s noise floor, the coherent integration period, the frontend filter, and the quantizer configuration. In an ideal receiver, the value of can be related directly to the received carriertonoisedensity ratio (often denoted or CNR) and the coherent integration time. Similarly, if the losses induced by factors such as frontend filtering and quantization can be modeled as a single loss value, denoted here by , then the following approximation is valid [3, 8]:
The advantage of using as a signal quality metric, as opposed to , for example, is that it reflects all of the signal processing effects applied to the received signal. Therefore, the performance of various discriminators can be related to one signal metric, as opposed to the ensemble of quantities: , , , and, perhaps, others. Moreover, as will be shown in Section 5.2, to achieve a specified loop performance it may be necessary to maintain a particular value of and so, accordingly a designer may wish to adjust , given a particular .
As a numerical example, consider a typical received GPS L1 C/A signal under open sky conditions and a typical consumer grade receiver. The received signal power using a patch antenna will be approximately −160 dBW and the receiver thermal noise floor may be assumed to be −205 dBW/Hz. If the receiver employs a frontend filter with a 2 MHz bandwidth and a onebit quantizer then the combined receiver 10 processing losses will be approximately 2 dB. Tallying these figures and assuming a coherent integration period of 1 ms, the expected value of is 16 dB. Alternatively, if a coherent integration period of 20 ms was assumed then the expected value of would be 29 dB. Under weaksignal conditions, however, such as the indoor environment, can fall to 0 dB and below.
3.2. The FourQuadrant Arctangent Discriminator (Atan2)
The fourquadrant arctangent discriminator is defined as [8]and is a purePLL discriminator, appropriate for pilot signals, or when data wipeoff is employed. The mean response of this discriminator to phase error, denoted here by , is well known (see, e.g., [8–10]). For high values of , is relatively linear across a wide range of phase error values and has approximately unity gain. As the value of reduces, the gain reduces considerably and the linear region diminishes. Expressions describing the exact mean response of the discriminator and its variance are developed as follows.
If the correlator values are interpreted as a complex pair, , then the argument, , can be shown to be distributed according to the probability density function, , defined as [9]Using , it can be readily shown that the value of is given by
From (9), taking the first derivative of (15) and setting , the discriminator gain is found to bewhere is the first derivative of with respect to , given by
A plot of versus for this discriminator is shown in Figure 2. It is evident that for values below approximately 6 dB, the discriminator gain reduces rapidly. This reduction in discriminator gain has implications for the closed loop poles of (4), (5), and (6). This will be discussed further in Section 4.
In a similar fashion to the mean of the fourquadrant arctangent discriminator, the variance of this carrier phase estimate can be found viaThis variance estimate, however, is a function of . Assuming that the PLL is tracking with an approximately zero mean phase error, it is useful to linearize this estimate around a zero phase error, :
A plot of for this discriminator is shown in Figure 3. For convenience, approximate solutions to (16) and (19) are provided in the appendix. Unsurprisingly, the discriminator variance changes linearly with for high values, bearing the approximate relationship: . In this region, (14) is approximately Gaussian. At dB the discriminator variance exceeds with reducing . At this point, (14) has begun to resemble a truncated Gaussian distribution. As is further reduced, (14) falls below the curve and approaches a uniform distribution over the interval , reaching a maximum variance of . This nonlinear relationship between and has a significant impact on the performance of the PLL under weaksignal conditions and, in conjunction with the discriminator gain effects described earlier, can result in severely degraded tracking performance. These effects must, therefore, be considered in the design of the PLL and will be discussed further in Section 4.
3.3. The Arctangent Discriminator (Atan)
The arctangent discriminator is defined as [7, 8, 11]and is a Costas discriminator, suitable for use on datamodulated signals. Similar to the fourquadrant arctangent discriminator, for high values of , changes in a relatively linear fashion with changing and has approximately unity gain. As the value of reduces, the gain reduces considerably and the linear region diminishes. This occurs at a higher value for the arctangent discriminator than for the fourquadrant arctangent discriminator, owing to its smaller linear region (discussed further in Section 3.7).
Similar to Section 3.2, it can be shown that the mean response of the arctangent discriminator, after some simplification, is given bywhere the limits of integration have been manipulated such that the arctangent function and its arguments reduce to simple linear combinations of , , and .
Again, from (9), taking the first derivative of (21) and setting , the arctangent discriminator gain, , is found to beFigure 2 depicts the relationship between and for this discriminator. For values below approximately 10 dB, the discriminator gain reduces rapidly. Although the trend is similar to that of the fourquadrant arctangent discriminator, it occurs at a higher value and the reduction in with is greater.
Similar to the fourquadrant arctangent discriminator, the variance of this carrier phase estimate, linearized around a zero phase error, can be found viaFigure 3 illustrates this relationship across an appropriate range of values. Again, similar to the fourquadrant arctangent discriminator, the discriminator variance changes linearly with for high values. As is reduced, (14) approaches a uniform distribution over the interval and reaches a maximum variance of . Once again, approximate solutions to (22) and (23) are provided in the Appendix.
3.4. The Quadrature Discriminator ()
The quadrature discriminator is defined, as its name suggests, as which is a purePLL discriminator, appropriate for nondata modulated signals. It is also notable both as being the phase detector used in the earliest PLLs; and as being the only discriminator for which the resultant PLL admits tractable nonlinear analysis [12, 13]. This discriminator function is, by far, the simplest form of carrier phase estimator. Owing to its simple definition, the characteristics of this discriminator are quite easily expressed, having a mean value, , ofUnlike the arctangentbased discriminator functions, this phase estimate is not selfnormalizing; that is, the estimate is a function of the nominal received signal power. To use this discriminator, even for high values (where the arctangentbased discriminators are completely selfnormalizing), this phase estimate must be normalized by an estimate of . The discriminator gain, , is given bywhich, unlike the previous three discriminators, is independent of . The variance of this carrier phase estimate is given by
Note that when the carrier phase estimate has been correctly normalized by , then the variance of this carrier phase estimate is given by . This curve is illustrated in Figure 3, providing a comparison with the arctangent based discriminators.
3.5. The DecisionDirected Quadrature Discriminator ()
The decisiondirected discriminator is defined as [8]:
The purpose of the term in this discriminator function is to render it insensitive to data modulation. The term provides an estimate of multiplied by the data value while the term provides an estimate of . As , this discriminator is (ideally) insensitive to data modulation. The value of can be found fromwhere denotes the expectation operator. Since and are statistically independent, thenFinding is equivalent to the estimation of given the AWGN corrupted sample: . It can, thus, be readily shown that [7]
The gain, , of this discriminator can be shown, using (9) and (31), to beA plot of versus is shown in Figure 2.
Unlike the arctangent based discriminators, the variance of the decisiondirected discriminator can be readily related to . From (10) and (28), we find
As with the quadrature discriminator, once the carrier phase estimate has been correctly normalized by , then the variance of this carrier phase estimate is given by , as depicted in Figure 3.
3.6. The GaintoNoise Ratio (GNR)
Sections 4.3 and 5 will illustrate that the performance of the tracking loop can be related, amongst other things, to the gain and variance of the discriminator. In fact, it will be shown that, under steadystate conditions, it is directly related to the ratio of the square of the gain to the variance. It is useful, therefore, to consider this ratio as a metric by which the tracking capability of each discriminator can be compared. This metric, termed the gaintonoise ratio and denoted by , is defined asA plot of for each of the four discriminators is shown in Figure 4. As can be expected from the analysis presented in Section 3, under high conditions ( dB, e.g.), the for each discriminator is similar. The reason for this is that for these high values and for each of the four discriminators. For reduced conditions, however, the unique relationship between , , and for each discriminator becomes evident. As the gain and variance characteristics of each discriminator are different, the curves diverge as is reduced. Because the tracking capability of the PLL can be directly related to the (as will be shown in Section 4.3), it provides insight into the relative tracking performance of each discriminator. Interestingly this metric also provides some insight into the relative benefits of a pilot signal, as it is clear that the of the purePLL discriminators is noticeably higher than that of the Costas discriminators for values below approximately 5 dB, as will be discussed further in Section 5.4.
While the can provide valuable insight into the operation of the PLL in its linear region, it does not completely characterize the discriminator’s influence on closed loop operation, as will be discussed next.
3.7. The Discriminator Linear Region
The linearized discriminator model employed in previous sections is an optimistic performance model and, generally, is only accurate for a limited range of . This range is termed the linear region and is finite for all discriminators. Indeed, it is ultimately limited to the range , owing to the periodic nature of the sinusoid. As has been shown in the previous sections, the values of and are dependent on the discriminator function, and, with the exception of the quadrature discriminator, are also dependent on .
In general, the linear region is symmetric around the origin and so it can be defined by the single scalar LR such that the linear region is the interval . LR is defined as the value of at which the true value of and the approximation differ by a certain percentage. The percentage is chosen arbitrarily, often depending on the application, but typical values are 5% and 10%. Specifically, LR for an linear region, denoted , is defined as where the notation can be interpreted as equals values of in the set such that condition is satisfied.
It is, thus, the intersection of the mean discriminator curve and the line that defines the linear region.
A plot of versus is shown in Figure 5. Examining the purePLL discriminators, it can be seen that the fourquadrant arctangent discriminator has a significantly larger linear region than the quadrature discriminator over the entire range of interest. For the Costas case, the arctangent discriminator has a significantly larger linear region than the decisiondirected discriminator, for high values. For values below approximately 7 dB, however, the linear region of both discriminators converge. The implications of the specific LR values and their dependence on will be discussed further in Section 5.
4. Closed Loop Operation
This section examines the closed loop operation of the PLL, specifically investigating the relationship between the discriminator gain and variance and their dependence on the closed loop transient and steadystate behavior. The relationship between the tracking bandwidth and the dependent discriminator gain is examined theoretically and the resultant influence on the transient response is illustrated via simulated phase steptests. In terms of steadystate performance, the significance of the metric as a means of predicting thermal noise induced tracking error is examined.
4.1. Tracking Bandwidth
The design of PLL loop filters is often a delicate balance between a sufficiently fast loop to cope with satellitetouser dynamics, and a sufficiently slow loop to resist thermal noise induced tracking error. It is crucial, therefore, that a designer has control over the exact placement of the loop poles. In general, direct specification of the closed loop poles is not intuitive and, so, the pole placement is often specified in terms of damping coefficient and equivalent bandwidth (secondorderdominant systems are generally parametrized in terms of the system damping coefficient, , and natural frequency, , each of which a uniquely observable effect in the timedomain [5]; when considering the system in terms polezero placement, however, it is more convenient to reparametrize the system in terms of the dominant poles, given by which can be related via and .)
The effective twosided rectangular bandwidth of the closed loop transfer function, denoted here by , is defined asGenerally, will be lowpass, with a relatively smooth passband. , therefore, is indicative of the speed at which the PLL will settle.
It is important to note the presence of in both the numerator and denominator of in (4). As has been shown in Section 3, is dependent on the prevailing . For low values, the value of is less than unity and, therefore, the effective bandwidth of the PLL will be less than its design value. As will be shown in Sections 4.2 and 5.1, this can have a significant impact on the overall loop performance.
In light of this effect, it is convenient to denote the high value (or the design value) of by and to define it as
Examining once again, it is clear that the denominator can be rendered independent of by scaling or, more specifically, the gains, , by a factor . To implement this gaincompensation the receiver must estimate the prevailing value of and calculate the value of corresponding to this value and the particular discriminator employed. This approach maintains the value of at the prescribed value , regardless of the prevailing . Note that although this gaincompensation modifies the loop filter gains as the prevailing signal strength changes, it is not an adaptive loop; it merely corrects for gain degradation in order to maintain a constant loop bandwidth. The effects of and benefits of this gaincompensation are explored next.
4.2. Transient Response
This section examines the transient response of a secondorder PLL employing an arctangent discriminator and a loop update rate of ms under both high and low conditions. The loop filter, detailed in Table 1, effects a critically damped system, , with a tracking bandwidth of Hz. The system was excited by a simultaneous step in phase of −0.12 rad and a step in frequency of −1.2 rad/s. An example of a simulated response of the PLL to this excitation, for a signal received at an of 23 dB, is plotted in Figure 6(a) and labeled “High.” As can be seen, the PLL exhibits a smooth, critically damped response which settles to within 5% of its peak value within 0.5 s.

(a) One instance of the PLL response
(b) The mean response calculated over 500 trials
To illustrate the impact of the dependent on the transient response of the PLL, this simulation was repeated under “Low” conditions. The particular case of of 0 dB was chosen as it corresponds to for the arctangent discriminator (see Figure 2). An example of the PLL response in this case is shown in Figure 6(a) and labeled “Low.” It is evident, apart from the increased noise, that the response of the PLL has become slower and more oscillatory.
This transient was simulated a total of 500 times, for both the high and the low cases and the average response was calculated and is presented in Figure 6(b), labeled “High” and “Low,” respectively. Indeed, it can be seen that the reduction in has induced a slower and underdamped response. Given , this has been calculated to be .
Using (36), was calculated for both the high case and the low cases to be 10 Hz and 5.28 Hz, respectively. It is clear that the value of is significantly reduced by the reduction in , an observation which agrees with Figure 6(b), in the sense that reduced results in an increased mean time to settle.
To eliminate this effect, the gain compensation discussed in Section 4.1 was applied to the PLL and the low scenario of 0 dB was reprocessed. The parameters of the compensated loop are presented in Table 1, in the column labeled “Comp.”. An example of the compensated loop response is shown in Figure 6(a). Again, the mean value of this response is estimated over 500 trials and is plotted in Figure 6(b). It is clear that the mean response is restored to that of the high case. Restoration of the PLL transient performance does, unfortunately, come at a price. It can be seen in Figure 6(a) that the compensated loop response exhibits significantly more thermal noise induced tracking error. This has implications for the steadystate operation of the PLL and is discussed next.
4.3. Tracking Error/Jitter
Following the transient response of the PLL, once the signal parameters (phase, Doppler, and higher order effects) have been estimated, the PLL settles and tracks the carrier phase. This socalled steadystate performance is, typically, dominated by thermal noise. The performance of the PLL in the presence of thermal noise can be measured in terms of the steadystate tracking error variance or tracking jitter, denoted here by . In the case of the PLL, the noise which corrupts the estimate of the carrier phase is and has propagated through the discriminator. Similar to the thermal noise floor, , it is convenient to consider an equivalent noise floor for the tracking error estimate, . Denoted here by , the noise floor of the phase error estimate, in rad^{2}/Hz, is defined asNote that, unlike , is defined as a twosided PSD. Given the transfer function and (38), the tracking error variance can be estimated as
Although can be well approximated by for high values of , for lower values of their values diverge and alone cannot be used to predict closed loop performance.
Using (39), the impact of on the noise performance of the PLL can be examined. As discussed in Section 4.2, it is necessary to compensate for the induced reduction in by increasing the filter gains, , by a factor . Using such gaincompensation, from (6), is given by
This result implies that, given perfect compensation, is equal to a constant term, , divided by the ratio . Although the ratio appears in (40), its reciprocal is chosen as the definition of the ; this is done so that conforms with metrics such as , , and , where the numerator pertains to the signal and the denominator pertains to the noise; also, this definition of can be used as a measure of usefulness; the higher the , the more useful the discriminator estimate. As the name suggests, the constant term, , is chosen by the designer. The ratio, , is related to via a function which is particular to each discriminator, defined earlier as .
To illustrate the usefulness of the in predicting the relative closed loop performance of various discriminators in the presence of thermal noise, the tracking error variance of a simulated PLL was measured for each of the four discriminators, across a range of conditions. The loop filter configuration of Table 1 was used and perfect compensation was applied to the loop filter gains for each case. A total of 29 conditions were simulated, ranging from −5 dB to 23 dB which corresponds to a range of to dB Hz and ms, for each of the four discriminators. The results of the MonteCarlo simulations are presented in Figure 7. Using (40), the theoretically predicted variance was calculated and is also plotted in Figure 7, exhibiting good agreement with the simulation results. For values below approximately 0 dB, the simulation results for the arctangent and decisiondirected discriminator have been omitted. In these cases, the PLL has lost lock and the resulting measurements of tracking error variance are meaningless.
Examining the relative performance of the four discriminators in Figure 7, we see that, for high values, all four discriminators perform equally well. As the is reduced, however, the individual characteristics of each discriminator influence the performance. These trends compare well with those observed in Figure 4. In fact, from (40), the relative relationship is identical as the curves of Figure 7 are simply the reciprocal of the curves of Figure 4 multiplied by the constant . It is noticeable that, for dB, both of the Costas discriminators perform more poorly than the purePLL discriminators. This is to be expected and is the unavoidable cost of achieving insensitivity to data modulation.
It is worth mentioning how the tracking jitter curves presented here compare with the traditional theory, which offers two different equations, one representing the class of purePLL discriminators, both the quadrature and fourquadrant arctangent, and another representing the Costas discriminators, including the decisiondirected quadrature discriminator and the arc tangent discriminator.
Included in Figure 7, for comparison purposes, is a plot of the traditional theoretical performance estimate, representing the general class of Costas discriminators [8, 16]. While this offers a reasonable fit for high values and aligns within 15% of the measured performance of the decisiondirected quadrature discriminator, it diverges from the measured performance of the arctangent discriminator with reducing , being in error by over 50% by an of 4 dB, and further diverging below this value. Interestingly, the traditional theory describing the performance of the class purePLL discriminators coincides exactly with that presented here for the quadrature discriminator. However, it offers a very poor fit to the performance of the fourquadrant arctangent discriminator being in error by 50% by an of 7 dB.
It is reasonable, therefore, to conclude that a comparison of the relative closed loop tracking performance of PLLs which employ the same loop filter, but different discriminators, can be inferred directly by simply examining the relative GNR of the discriminators. That is, the relative linear closed loop performance of two PLLs, for any loop configuration, can be inferred by simply examining the open loop behavior of their respective discriminators. This will be discussed further in Section 5.2 and the usefulness of the in choosing a particular discriminator for a given application will be discussed.
5. Applications to Receiver Design
This section discusses applications of the theory developed in the previous section to GNSS receivers in the context of initial design choices and runtime receiver tuning. Firstly, the importance of acknowledging the dependence of the discriminator gain on the prevailing signal conditions and the benefits of compensating for this gain are considered by examining real GPS L1 C/A for a pedestrian navigation scenario. Secondly, the problem of choosing an appropriate discriminator, given a receiver configuration and received signal strength, is addressed by utilizing the and the linear region metrics. Finally, and once again employing these two metrics, the issue of optimal combining of carrier phase error estimates in datapilot systems is examined using the Galileo E1 B/C signals as a case study.
5.1. Maintaining the Design Loop Bandwidth
The impact of the discriminator gain on the performance of a GPS L1 C/A tracking loop is examined here in the context of pedestrian navigation. The experiment encompassed a range of conditions and considered both tracking loops which employ gaincompensation and those which do not. Results show that loops which compensate for exhibit significantly improved cycleslip performance.
A set of IF data was collected using a GPS1A frontend and an Antcom antenna [17, 18] which logged twobit IF samples at a rate of 16 Ms/s and employed a 2 MHz front end filter. The antenna and receiver were mounted on a rigid body and carried in the pedestrian’s hand. Under opensky conditions, the subject initially stood for one minute and subsequently traversed a 150 m eastwest path, repeatedly, at a steady walk for a period of four minutes. The antenna was maintained approximately level for the duration of the experiment and, being handheld, the antenna, oscillator, and receiver experienced the typical dynamics of a pedestrian including gross velocity of each traversal and the transient, stepinduced accelerations.
One particular satellite, PRN 17, was observed at azimuth and elevation of approximately and , respectively, and a received of approximately 46 dB Hz. This signal was tracked using a typical tracking configuration, consisting of a secondorder 20 Hz PLL using an arctangent discriminator, defined by (20), and a 0.5 Hz secondorder delaylock loop (DLL). Both tracking loops used a 1 ms update rate. A secondorder noncarrieraided DLL was chosen to ensure that carrierphase tracking errors, induced by the oscillator and pedestrian dynamics, were not propagated to the code tracking loop.
The observed carrier Doppler is presented in Figure 8. The first sixty seconds represent the stationary part of this experiment where only the satelliteinduced Doppler is evident. The remainder of the data represents walking dynamics where both the satellite and pedestrianinduced gross velocities contribute the observed Doppler. In addition to the gross Doppler, the transient accelerations associated with walking have induced quasisinusoidal perturbations to the observed Doppler via the socalled gsensitivity of the oscillator [19]. Typically, a temperature compensated crystal oscillator (TCXO) used for handheld GNSS applications will exhibit a gsensitivity of the order of 1.5 to 2.5 ppb/g, while specialized low gsensitivity oscillators are in the range of 0.35 to 0.5 ppb/g (see, e.g., [15]). The dynamics of a walking stride can be expected to induce acceleration peaks and troughs of approximately 8.0 and −6.0 m/s, respectively [21]. Given these values, the Doppler perturbations visible in Figure 8 appear consistent with what would be expected for a lowpower lowcost device.
To observe the behavior of the carrier tracking loops under weaksignal conditions, the IF data was attenuated prior to reprocessing. This attenuation was achieved by adding white Gaussian noise directly to the IF samples such that the noise power spectral density in the vicinity of the carrier frequency was increased by the required amount. The tracking loops were initialized using the carrier frequency, carrier phase, and code phase estimates gained from the reference, unattenuated trial. An estimate of the tracking performance was then made by comparing the carrier phase estimate of the PLL during the attenuated trial to that of the reference trial. This experiment was then repeated for a selection of signal attenuation values for both the gaincompensated and nongaincompensated PLLs. Specifically, the data was processed for each of 9, 12, 15, and 18 dB of attenuation, which corresponds to average values of 12.1, 9.1, 6.1, and 3.1 dB, respectively. These values of attenuation were chosen to cover an interesting range of discriminator behavior, including the transition from unity gain to progressively reducing gain, including the onset of discriminator variance saturation and including the steepest region of contraction of the linear region. In this way, it is expected that the performance should degrade rapidly with increasing attenuation level and that the application of gaincompensation should improve, to some extent, the performance.
An estimate of the carrier discriminator gain was produced within the gaincompensated PLLs by applying a standard estimator to the correlator values, and [8], and using this estimate in conjunction with the equations provided in the appendix. Details of the accuracy of estimation and the relative sensitivity of the PLL tuning are discussed further in Section 5.3. Figures 9(a) and 9(b) show the measured phase error for the nongaincompensating and the gaincompensating loops. Apart from the obvious observations that cycleslips only occur once the pedestrian has begun to walk (from 60 seconds onwards), and that cycleslips are more frequent in the more highly attenuated trials, there are some more interesting features of these results.
(a) Nongaincompensated PLL
(b) Gaincompensated PLL
Firstly, it is clear that the gaincompensating loop exhibits significantly less slips than the nongaincompensating loop. This is due to the fact that the gain reduction induces a slower response to changes in the received phase, thereby resulting in a failure to adequately track the phase trajectory. This observation is supported by the results presented in Table 2, which shows the number of measured halfcycleslips for five visible satellites. The tabulated data is arranged as follows: each row represents a single satellite, the first and second columns of each row are the PRN and the at which the signal was observed prior to attenuation. The remaining columns represent the total number of halfcycleslips observed during the attenuated trial, with a pair of numbers per attenuation value. The leftmost number represents the number of halfcycleslips observed on the nongain compensated PLL, while the rightmost, italic number represents the cycleslips observed by a PLL implementing live gaincompensation. Secondly, considering the and attenuation numbers from Table 2, it is evident that the benefits of gaincompensation are most pronounced within an range of 5.0 to 12.0 dB, which corresponds to the point at which the linear region of the discriminator begins to contract. In this range, the PLL is most sensitive to large sustained phase errors, resulting from a low discriminator gain, as it drastically increases the probability of a cycleslip. For values below this range, the increased noise present on the phase error estimate, the of the discriminator as significantly reduced, and the contribution of thermal noise error becomes significant. Ultimately, of course, the design bandwidth of the PLL ought to be reduced to effect a more reasonable tradeoff between dynamic and thermal noise errors.

Although this particular experiment only investigates the arctangent discriminator, the general results support the observations made in Section 4.2 and suggest that a similar trend may be observed in the case of other discriminators which exhibit low induced gaindegradation. It should be noted that gain compensation is employed exclusively here; however, in some cases, the problem of gain degradation can be circumvented by simply increasing the coherent integration, thereby increasing and placing the discriminator in its unitygain region. Unfortunately, this approach is not always possible. Firstly the integration period may be limited by data modulation and, secondly, either local oscillator instability and/or excessive user dynamics can induce sufficiently rapid phase variations as to necessitate a high loop update rate to maintain phaselock; that is, net phase dynamics may limit the integration period. Under these circumstances, gain compensation can prove useful.
It is worth commenting on the difference between gaincompensation, as implemented in this experiment, and traditional gainscheduling or adaptive filtering. The process of gain compensation maintains a constant tracking bandwidth, , across a range of conditions. In contrast, gainscheduling prescribes a particular loop filter which is deemed appropriate for the prevailing signal conditions and an adaptive filter will modify its filter parameters in response to features of the received signal (, e.g.) [20]. Gaincompensation does not adapt the PLL bandwidth, it ensures that it remains constant and equal to the design bandwidth, . With this in mind, it is clear that the effective implementation of a gainscheduled or adaptive PLL must consider the effect of on the discriminator and employ appropriate gaincompensation.
5.2. Choosing a Discriminator for Linear Operation
Section 4.3 has shown that the closed loop tracking jitter observed in a PLL can be related directly to the and the PLL bandwidth. By examining the relative values of different discriminators, in Figure 4, in conjunction with their linear regions, in Figure 5, it is possible to choose a discriminator which will minimize for a given loop filter choice.
The first, perhaps obvious, conclusion that can be drawn from these figures is that the better of the two purePLL discriminators always outperforms the better of the two Costas discriminators, in terms of and linear region. Therefore, if the received signal is not data modulated, or if the modulation is known, then one of the purePLL discriminators will always yield the better steadystate tracking performance. A Costas discriminator should only be used when necessitated by the presence of unknown datamodulation. Thus, the choice of discriminator should then be considered for two different discriminator classes, namely, purePLL or Costas.
For the purePLL discriminators, under high conditions (>11 dB), the fourquadrant arctangent discriminator incurs less than a performance degradation, when compared with the quadrature discriminator, yet it exhibits a significantly larger linear region. The fourquadrant arctangent discriminator should, therefore, be used in this region as it provides more robustness than the quadrature discriminator, being capable of absorbing larger phase transients while maintaining linear operation.
In the region −3 dB < < 11 dB, the optimum choice of discriminator may be dependent on the application, the quadrature discriminator significantly outperforms the fourquadrant arctangent discriminator in terms of tracking error but has a notably narrower linear region. For applications where low tracking error is the main priority, the quadrature discriminator should be used whereas, if resilience to signal dynamics is desired, a designer may wish to avail of the larger linear region of the fourquadrant arctangent discriminator. For very low values (−3 dB), the linear regions of both discriminators are similar, yet the quadrature discriminator provides approximately 3 dB less tracking error variance and should, therefore, be used.
Unlike the purePLL discriminators, the choice of discriminator is simpler for the Costas case. At dB, the linear regions of the arctangent discriminator and the decisiondirected discriminator begin to converge. Also, for reducing values around this point, the of the decisiondirected discriminator begins to significantly outperform the arctangent discriminator. Thus, for values above approximately 9 dB, the arctangent discriminator should be used while, for values below this point, the decisiondirected discriminator should be employed.
Although these conclusions have been drawn from inspection of Figures 4 and 5, they can also be inferred from inspection of the MonteCarlo simulation results presented in Section 4.3.
5.3. A Note on Estimation
As the configuration and tuning of the PLL described here are based on the premise that the prevailing is reasonably well known, it is worth briefly commenting on the sensitivity of the PLL tuning to errors in the estimate of . Here, the experiment described in Section 5.1 is taken as an example. In this case, a runtime estimate of was generated using the wellknown estimator described in [8] and further smoothed by a 1 second movingaverage filter. Ideally, the choice of smoothing applied to the estimate should reflect a reasonable tradeoff between noiserejection and the speed of response to changes; however, this empirically derived configuration has proven effective.
Recall that the IF data was digitally attenuated by a precise factor for each trial. A reference measurement of the original unattenuated was taken, and being a very high value of approximately 22 dB, it was considered to be an errorfree estimate. Then, for each attenuated trial, the difference between this reference value and the runtime estimate, minus the applied attenuation, was recorded. This represented the error in the runtime estimate. A plot of the measured standard deviation of the error is shown in Figure 10, along with the CramérRao Lower Bound standard deviation for nondataaided BPSK estimation [22, 23].
In terms of sensitivity to errors in the estimate of , the expressions for discriminator gain, , and tracking bandwidth, , presented in the appendix, can be used to explore the how accurately the loop bandwidth can be restored under low conditions. Assuming a secondorder PLL and the arctangent discriminator, and using (A.3) and (A.5), the envelope of was computed for , using the the CRLB shown in Figure 10 as a reference. These envelopes are shown in Figure 11 along with the tracking bandwidth for the cases of perfect compensation and of no compensation. Interestingly, even for very weaksignal conditions, down to an of 0 dB, the bandwidth can typically be restored to within one Hertz of its design value. For lower values, however, the error becomes noticeable, suggesting that more averaging should be applied in the signaltonoise ratio estimator, in order to provide a less noisy estimate.
5.4. Choosing a Discriminator for Data/Pilot Tracking
This section examines the problem of carrier tracking for a datapilot signal structure, specifically, the Galileo E1 B/C signal is taken as an example. Of particular interest is the scenario when the receiver has synchronized with the secondary code on the pilot component of the signal and is capable of combining both a Costas (E1B) and a purePLL (E1C) estimate of the carrier phase error. The benefits of using the GNR and linear region analysis presented in Section 3, when choosing weights for estimate combining, will be illustrated and some further considerations will be discussed.
A fiveminute IF dataset was collected using a roofmounted antenna during which time the Galileo Prototype Flight Model (PFM) satellite was broadcasting on PRN 11. A complex sample rate of 20 MHz was used and data was collected with a fourteenbit quantizer resolution. The received signal was observed at a of approximately 47 dB Hz.
Reference carrier phase and frequency trajectories were extracted from this dataset for use in the following experiments. This was done by processing the data with a standard pilotonly tracking architecture, comprising a 4 Hz PLL bandwidth operating with an update rate of 20 ms, combined with a 0.1 Hz PLLassisted DLL. As the antenna was stationary and the reference oscillator was particularly stable, the use of a lowbandwidth PLL and long coherent integration period ensured that these reference measurements were of high accuracy. An attenuated copy of this data was then processed using different and pilotonly or data/pilot architectures, using a more typical PLL design. The difference between the estimated carrier phase and frequency for these architectures and that of the reference was used as an indication of relative performance.
When processing a data/pilot signal, a PLL can either produce phase estimates using the pilot signal alone, or combine estimates from both the data and the pilot signal (see, e.g., [10, 14]). When two estimates are combined, they can be weighted such that the tracking jitter is minimized. The combined estimate and the associated tracking error variance for such a combined estimate can be expressed aswhere denotes the weight applied to each estimate and the subscripts and denote data and pilot signals, respectively.
Equating the partial derivative of (42) with respect to to zero and noting that , the (rather intuitive) set of weights which minimize the tracking jitter can be shown to be
It is worth commenting that this result differs from previously reported [10, 14] weighting guidelines which recommend that weights are chosen based upon discriminator variance. Variancebased weighting is inappropriate for some discriminators, such as the arctangent discriminator, as the variance saturates to a moderate value for low while the gain continues to reduce. Variance alone, therefore, does not reflect the true usefulness of the discriminator. based weights, as prescribed by (43), consider both gain and variance and, thus, yield superior performance. Note also that (43) implicitly considers the coherent integration period, which does not need to be equal for both the data and the pilot signals, as the is a direct function of .
Thus, a data/pilot architecture employing (43) should choose the appropriate discriminator for each of the pilot and data signals separately, based on the prevailing and using the guidelines presented in Section 5.2. Subsequently, the combining weights should be calculated based on the values of each of the chosen discriminators. This composite phase error estimate can then be passed to the loop filter.
Note that for very low conditions the receivers estimate of can become noisy and unreliable, as shown in Section 5.3. Thus, in certain cases, it may be beneficial to consider alternative architectures; for example, the pilotonly approach which neglects entirely [10]. Alternatively, the weights can be formed based upon . Noting that, under high conditions, the becomes approximately linearly proportional to , the weights could be computed by replacing in (43) with the coherent integration period of the corresponding signal.
To examine their relative performance, the attenuated IF data was processed with each of a pilotonly PLL and both a timebased and a based combining PLL. A critically damped, 10 Hz loop was employed in all cases, and a coherent integration period of 4 ms was used. The attenuation was timevarying, beginning at 0 dB and increasing at a rate of 0.1 dB/s to a final value of 30 dB at five minutes. The variance of the difference between the carrier phase of the attenuated data and that of the reference was calculated over a 30 second window for each PLL configuration. A plot of the measured tracking error variance versus the average over each 30 second window is shown in Figure 12.
It is clear from the measured results that when appropriate weighting is employed, the based data/pilot PLL outperforms both of the other candidate architectures, specifically in the range dB reaching almost 3 dB. Perhaps more interesting, however, is the relative performance of the pilotonly and the timebased data/pilot schemes. For values higher than 5 dB, the timebased and based architectures perform equally well. This is because the arctangent and fourquadrant arctangent have equal in this region, and the weights in each case are equal. Indeed, the divergence in performance at dB coincides with the divergence in shown in Figure 4, for very low values, the respective performance of the pilotonly and that of the based PLL converge. At this point, the difference in between the purePLL, used for the pilot signal, and Costas, used for the data signal, is so large that is almost zero.
A number of conclusions can be drawn from this experiment. Firstly, it is clear that there may always be an advantage to utilizing the datasignal for carrier phase estimation, provided the estimate can be appropriately weighted.
However, it is evident that the incremental benefit diminishes rapidly for very weaksignals, becoming effectively useless for values below approximately 0 dB. This observation is broadly in line with that of [10], which claim that a pilotonly scheme is optimal under weaksignal conditions.
It is evident, however, that inappropriate weighting can prove detrimental to receiver performance. Specifically, this occurs under low conditions, as evidenced by the performance of the timebased combining architecture which can perform more poorly than a pilotonly PLL. Secondly, it is clear that a reasonably wellperforming suboptimal architecture may be constructed by simply using a timebased combining datapilot PLL for strong and moderate signal strengths and a pilotonly PLL when the signal is weak. An appropriate threshold may, for example, be dB (i.e., dB Hz for ms).
6. Conclusions
Following a thorough analysis of carrier phase discriminators, it is evident that, under weaksignal conditions, traditional performance models fail to fully describe PLL behavior. Both MonteCarlo simulation and live signal tests appear to confirm that induced gain degradation is prevalent in some of these discriminators and that that has a significant impact on overall PLL performance. For the specific case of pedestrian navigation, it appears that the proposed gaincompensation technique can provide substantial performance improvements in terms of dynamic response and cycleslip frequency.
Results pertaining to the closed loop noise performance of the PLL, when operating in its linear region, illustrated that the represents a useful metric which can infer the relative closed loop performance of various discriminators, based on their respective open loop characteristics. Utilizing both this metric and the linear region analysis, experiments have confirmed that the choice of discriminator should consider the prevailing , via the discriminatorspecific function. Moreover, in terms of the design of data/pilot tracking architectures the usefulness of the metric in providing a discriminator weighting scheme appears to provide a corresponding improvement in tracking accuracy.
It is noteworthy that while the analysis presented here considered only four discriminators, the metrics, and the theoretical model employed (, GNR, and LR), can be extended to consider and provide a comparative analysis of a host of carrier phase estimators. Given an expression for these three metrics as a function of , this analysis could be extended to consider any memoryless discriminator.
Appendix
As the integral expressions for the statistics of the fourquadrant arctangent and the arctangent carrier phase discriminators do not appear to yield a closed form, a set of approximate expressions are presented here. The forms of the expressions have been chosen by inspection of numerical evaluations (9) and (10), for each discriminator and the coefficients (, , ) have been optimized to minimize the r.m.s error in the range dB < < 30 dB. Detailed also are the maximum error, of the approximate model, the value of at which this error occurs, and the standard deviation of the percentage error, denoted by , calculated across the entire fit range.
The gain, , and variance, , of the fourquadrant arctangent carrier phase discriminator can be well approximated bywhere and model errors are , and at dB andwhere , , and and model errors are and at dB. Similarly, the gain, , and variance, , of the arctangent carrier phase discriminator can be well approximated bywhere and model errors are and at dB and where and and model errors are and at dB. Approximate expressions for the of the arctangent discriminators can be found by substituting the above expressions into (34). One further interesting result is the solution to (36), given the filter (8) and . This expresses the bandwidth of a secondorder PLL which uses a proportional and integral controller and is given by
Explicit design equations for seconderorder filters can be found in, for example, [5, 6].
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
References
 H. de Bellescize, “La réception synchrone,” L’Onde Électrique, vol. 11, pp. 230–240, 1932. View at: Google Scholar
 J. P. Costas, “Synchronous communications,” Proceedings of the IRE, vol. 44, no. 12, pp. 1713–1718, 1956. View at: Publisher Site  Google Scholar
 J. T. Curran, D. Borio, G. Lachapelle, and C. C. Murphy, “Reducing frontend bandwidth may improve digital GNSS receiver performance,” IEEE Transactions on Signal Processing, vol. 58, no. 4, pp. 2399–2404, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 W. Zhuang, “Performance analysis of GPS carrier phase observable,” IEEE Transactions on Aerospace and Electronic Systems, vol. 32, no. 2, pp. 754–767, 1996. View at: Publisher Site  Google Scholar
 S. A. Stephens and J. B. Thomas, “Controlledroot formulation for digital phaselocked loops,” IEEE Transactions on Aerospace and Electronic Systems, vol. 31, no. 1, pp. 78–95, 1995. View at: Publisher Site  Google Scholar
 J. T. Curran, G. Lachapelle, and C. C. Murphy, “Improving the design of frequency lock loops for GNSS receivers,” IEEE Transactions on Aerospace and Electronic Systems, vol. 48, no. 1, pp. 850–868, 2012. View at: Publisher Site  Google Scholar
 J. G. Proakis, Digital Communications, Electrical Engineering Series, McGraw Hill International Editions, 3rd edition, 1995.
 E. D. Kaplan, Ed., Understanding GPS: Principles and Applications, vol. 1, chapter 5, Artech House, 2006.
 W. Hagmann and J. Habermann, “On the phase error distribution of an open loop phase estimator,” in Proceedings of the IEEE International Conference on Communications (ICC ’88), vol. 2, pp. 1031–1037, IEEE, Philadelphia, Pa, USA, June 1988. View at: Publisher Site  Google Scholar
 O. Julien, “Carrierphase tracking of future data/pilot signals,” in Proceedings of the 18th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS ’05), pp. 113–124, Long Beach, Calif, USA, September 2005. View at: Google Scholar
 A. J. Van Dierendonck, “GPS receivers,” in Global Positioning System: Theory & Applications, vol. 1 of Progress in Astronautics and Aeronautics, chapter 8, pp. 329–408, AIAA (American Institute of Aeronautics & Astronautics), 1996. View at: Google Scholar
 A. Viterbi, “Phaselocked loop dynamics in the presence of noise by fokkerplanck techniques,” Proceedings of the IEEE, vol. 51, no. 12, pp. 1737–1753, 1963. View at: Publisher Site  Google Scholar
 F. M. Gardner, Phaselock Techniques, Wiley, 3rd edition, 2005.
 A. Jovanovic, Y. Tawk, C. Botteron, and P.A. Farine, “Dual channel optimization of tracking schemes for E1 CBOC signal,” in Proceedings of the Vehicular Technology Conference (VTC Fall ’11), pp. 1–5, IEEE, San Francisco, Calif, USA, September 2011. View at: Publisher Site  Google Scholar
 Vectron International, “Temperature Compensated Crystal Oscillator: Low GSensitivity,” July 2010, http://www.vectron.com/products/tcxo/tx508.pdf. View at: Google Scholar
 P. Groves, Principles of GNSS, Inertial, and Multisensor Integrated Navigation Systems, Artech House, Boston, Mass, USA, 2013.
 GPS Creations, “GPS1A,” January 2012, http://www.gpscreations.com/Products GPS1A.html. View at: Google Scholar
 Antcom Corporation, “GPS L1/L2 Antenna P/N2DG1215AMNS4,” http://www.antcom.com/. View at: Google Scholar
 R. L. Filler, “The acceleration sensitivity of quartz crystal oscillators: a review,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 35, no. 3, pp. 297–305, 1987. View at: Google Scholar
 K. J. Åström and B. Wittenmark, Adaptive Control, chapter 9, AddisonWesley, 1989.
 Q. Ladetto, V. Gabaglio, and B. Merminod, “Combining gyroscopes, magnetic compass and GPS for pedestrian navigation,” in Proceedings of the International Symposium on Kinematic Systems in Geodesy, Geomatics and Navigation, pp. 205–212, Banff, Canada, AugustSeptember 2001. View at: Google Scholar
 S. Satyanarayana, D. Borio, and G. Lachapelle, “C/N0 estimation: design criteria and reliability analysis under global navigation satellite system (GNSS) weak signal scenarios,” IET Radar, Sonar and Navigation, vol. 6, no. 2, pp. 81–89, 2012. View at: Publisher Site  Google Scholar
 N. S. Alagha, “CramerRao bounds of SNR estimates for BPSK and QPSK modulated signals,” IEEE Communications Letters, vol. 5, no. 1, pp. 10–12, 2001. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 James T. Curran. 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.