Group Delay and Phase Delay in GNSS Systems
Group Delay and Phase Delay in GNSS Systems
To cite this article: Reuben Berry, Philip G. Mattos & Izzet Kale (2013) Group delay and
phase delay in GNSS systems, Geo-spatial Information Science, 16:3, 210-219, DOI:
10.1080/10095020.2013.834111
GNSS signals have previously been modulated using binary phase shift keying but this modulation scheme is being
replaced by binary offset carrier (BOC) modulation. Research has considered how the BOC signals might be affected
differently when passed through a surface acoustic wave (SAW) filter. The concern has been that because of the split
spectrum nature of the BOC signals, the upper and lower side-lobes will be delayed significantly differently. This was
suggested because SAW filters have nonlinear phase characteristics and therefore different frequencies are delayed
differently. It was suggested that this difference in delay will result in greater distortion of the correlation triangle.
A delay magnification effect was also mentioned when analyzing the delay of a BOC signal. It was not understood why
the theoretical delay calculations did not match up with the actual results in both hardware and simulation. This paper
clarifies some of the confusion and explains why the “delay magnification” applies to phase delay but not group delay.
This paper also takes a look at how the code phase delay can vary with frequency and correlator spacing as a result of
the SAW filter properties.
Keywords: group delay; phase delay; SAW filter; non-linear phase; GNSS delay compensation
1. A simple sinusoid model carrier and information signal sinusoid are clearly visible
on the plots.
To explain the delay magnification idea proposed in Ref.
Figure 2 shows the two sinusoids f1 and f2 . It clearly
(1), we start with a model of two sinusoids being passed
shows that the phases of the two sinusoids start aligned
through a filter. These two sinusoids are generated by
but gradually separate due to the difference in frequency.
multiplying a carrier sinusoid with a signal sinusoid and
If the sinusoid of frequency f2 was delayed by a
are the input to a filter:
small amount, we would see that the envelope of the sig-
nal would be delayed much more. Let the filter have no
xðtÞ ¼ cosð2pfc tÞ cosð2pfs t Þ ð1Þ delay for frequency f 1 but have a phase delay of 16 ns
at f2 . Figure 3 shows the output of the filter yðtÞ. This
where x(t) is the input, fs is the frequency of the infor- shows that the envelope of the signal has been delayed
mation signal which for our purposes is assumed to be a by approximately 250 ns even though the phase delay of
simple sinusoid, and fc is the frequency of the carrier. the filter for f2 was only 16 ns. This demonstrates the
The input to the filter is equivalent to two sinusoids idea of delay magnification, as proposed in Ref. (1),
added together: where a small delay results in a much larger delay.
We have already used the term “phase delay” in the
xðtÞ ¼ cosð2pf1 t Þ þ cosð2pf2 tÞ ð2Þ previous paragraph so we had better also define it here.
The delay of a filter can be measured in terms of the
where f1 ¼ fc fs and f2 ¼ fc þ fs . “phase delay” but more commonly the “group delay” is
In the original analysis in Ref. (1), when these two used instead. Both the phase and the group delays are
sinusoids (f1 and f2 ) were passed through a SAW filter, calculated using the phase of the filter. The equations
the sinusoids were delayed according to the group delay describing the different delay definitions were taken from
of the SAW filter at the two frequencies. Because of this, Ref. (2) and are defined here as:
any slight difference in group delay between the two fre-
quencies results in a much greater delay in the informa- h
sp ¼ ð3Þ
tion signal sinusoid. This effect is demonstrated below. 2pf
Figure 1 shows the input to the filter xðt Þ. In this
example, fs ¼ 1 MHz and fc ¼ 30 MHz. This results in dh
f1 ¼ 29 MHz and f2 ¼ 31 MHz. Although these are not sg ¼ ð4Þ
2pdf
GPS frequencies, they were chosen so that both the
Figure 1. Filter input signal xðtÞ. Figure 3. Output of the filter yðtÞ
y0 ð t Þ ¼ x0 t sp ð7Þ
y0 ðt Þ ¼ cosð2pf0 ðt sp ÞÞ ð8Þ
approximation to the actual response of the filter if the signals but it defines the delay of the envelope of a
group delay varies only very slightly over the frequency signal. Equation (10) also shows that the phase delay
range of ye ðt Þ. If the group delay of the filter varies defines the steady-state response of a filter.
greatly with frequency then the bandwidth of the signal
xe ðt Þ can be reduced until the changes in group delay are
small enough to be neglected. 2. Modeling using spread-spectrum signals
In the case of a linear phase filter, the values for the Up to this point, the signals used in the analysis have
phase delay sp are identical to the group delay sg values been useful for the initial explanation of delay but they
for all frequencies. In this case, even though the phase are inadequate for the analysis of how the real signals
and group delay definitions are not equivalent, the group and the correlation peaks will be affected by a SAW
delay values sg could be applied in the same way as the filter in the front end of a receiver. The signals used
phase delay values sp to get the same results. However, previously have consisted of pure sinusoids meaning that
one must be cautious because, in the case of a nonlinear they are continuous signals with no transients and there-
phase filter, the group and phase delay values are not the fore their spectra contained only two frequencies. Also
same and therefore, if applied incorrectly the calculations the signals previously used carried no real information.
will not match the measured or simulated results. GPS signals have a spread spectrum meaning that,
We now return to our original model of two their spectra do not contain just a few frequencies but,
sinusoids being passed through a filter as described by they use a range of frequencies to represent the signals.
Equations (11) and (12). The output of the filter yðtÞ is: In order to be able to model such signals, the modeling
was done in MATLAB. This allows the effects that a
SAW filter has on the signals and their correlation peaks
yðtÞ ¼ cosð2pf1 t þ h1 Þ þ cosð2pf2 t þ h2 Þ ð11Þ to be observed. For the analysis, both BPSK and BOC
(1, 1) signals were simulated separately by being passed
where h1 is the phase of the filter at f1 and h2 is the through a model of a SAW filter. The first SAW filter
phase of the filter at f2 . Therefore, in terms of the modeled was a B3521 GPS filter from EPCOS.
original carrier and information signal, the output of This shows whether or not the BOC(1, 1) signals will
the filter is: be delayed differently or will have a different distortion
of the correlation triangles when compared with the
BPSK signals.
h1 þ h2 h2 h1
yðtÞ ¼ cos 2pfc t þ cos 2pfs t þ The model of the signals and SAW filter calculates
2 2
the response of the filter for 1 ms of signal, this being
ð12Þ the length of a complete chipping sequence for the C/A
code. The signals were first generated in the time domain
This equation can be compared with Equation (10) and then passed through a model of the SAW filter by
where the sinusoid in (10) is the carrier signal in (12) multiplying the spectrum of the signals in the frequency
and the envelope function in (10) is the information domain with the frequency response of the SAW filter.
signal in (12). It should be noted that the carrier signal is
The data provided by the manufacturer for the SAW
the cosine with the fc term and the information signal is
filter contained the magnitude and phase response of the
the cosine with the fs term.
filter in 312.5 kHz steps. Because the simulation length
In Equation (10), the group delay sg of the filter
is 1 ms, the points in the frequency domain are spaced
delays the envelope signal envðtÞ. In Equation (12), the
every 1 kHz meaning that the SAW filter data had to be
information signal is equivalent to the envelope signal in
Equation (10). Therefore, the group delay sg of the filter
can be found from the phase of the information signal in
equation (12) and is defined as:
h2 h1
sg ¼ ð13Þ
2p2fs
up-sampled and interpolated using a first-order hold to domain using an IFFT. After this the signals are demod-
have data points every 1 kHz. ulated by multiplying them with the carrier signal. And
Once the signals had been passed through the SAW finally, the signals are correlated. The correlation is
filter by means of multiplication in the frequency achieved by multiplication in the frequency domain as
domain, the signals are converted back to the time this requires less processing than direct correlation in the
time domain.
Figure 5. Frequency spectrum of BOC(1, 1) signal. Figure 8. Close-up of correlation triangle for BPSK.
Figure 6. Correlation response of BPSK signal. Figure 9. Correlation response of BOC(1, 1) signal.
Figure 7. Correlation triangle of BPSK signal. Figure 10. Correlation triangle of BOC(1, 1) signal.
214 R. BERRY et al.
one from another to find the error. This error was found the signals at each stage in their modulation, filtering
to be approximately 8 105 relative to the actual value and demodulation. The example given in this section
at any given point. After zooming in on the errors, it uses low frequencies to aid explanation, but the process
was found that the errors were oscillatory and were at a is still the same for the actual GPS frequencies.
frequency twice the carrier frequency. This shows that Figure 14 shows the spectrum of a 10 Hz carrier
the first way of simulating the system that modulates the signal. Figure 15 shows the spectrum of the baseband
signal, filters, and then demodulates it, does not com- signal being transmitted. For both of these signals, they
pletely filter out the 2fc frequencies. This is because no have a purely positive and real spectrum and so only the
special filtering was used after demodulation. However, magnitude is shown.
when the demodulated signal is convolved with the Figure 16 shows the transmitted signal that is
flipped version of the original signal, this filters out generated by multiplying the baseband signal with the
the high frequencies. Because of this, it was assumed carrier. In the frequency domain, the two signals are
that the demodulated signal did not need to be low-pass convolved resulting in the spectra in Figure 16. Again
filtered first. the spectrum of the transmitted signal is purely positive
and real.
Figure 17 shows the spectrum of the filter’s response.
4. Modeling the system at baseband This time both real and imaginary parts are shown on
This section goes about to explain why it is possible to the same figure. When the transmitted signal is passed
model the system only at baseband frequencies rather through the filter it results in the signal shown in
than at RF frequencies. Modeling at baseband reduces Figure 18.
the sampling rates required and hence reduces computa- The spectrum of the modulated and filtered signal
tion time. To explain, we first look at what happens to shown in Figure 18 is next demodulated. This results in
Figure 14. Spectrum of carrier signal. Figure 16. Spectrum of transmitted signal.
Figure 15. Spectrum of baseband signal. Figure 17. Spectrum of filter response.
216 R. BERRY et al.
Figure 18. Spectrum of filtered signal. Figure 20. Spectrum of low-pass filtered demodulated signal.
Notes on contributors
Reuben Berry was awarded Master of Engineering (MEng)
Figure 26. Signal delay against correlator spacing. with Distinction in Electronic Engineering at the University of
Geo-spatial Information Science 219
Westminster in 2012. He is currently working for the Applied power DSP algorithms/architectures and sigma–delta modulator
DSP and VLSI Research Group (ADVRG) at the University of structures for use in the communications and biomedical
Westminster contributing to the design of a research enabling industries and hybrid Global Navigation Satellite Systems GPS/
GNSS receiver. Galileo/Glonass/Compass receiver structures and systems.