Group delay and phase delay


In signal processing, group delay and phase delay are functions that describe in different ways the delay times experienced by a signal’s various sinusoidal [|frequency components] as they pass through a linear time-invariant (LTI) system.
These delays are sometimes frequency dependent, which means that different sinusoid frequency components experience different time delays. As a result, the signal's waveform experiences distortion as it passes through the system. This distortion can cause problems such as poor fidelity in analog video and analog audio, or a high bit-error rate in a digital bit stream.

Background

Frequency components of a signal

Fourier analysis reveals how signals in time can alternatively be expressed as the sum of sinusoidal frequency components, each based on the trigonometric function with a fixed amplitude and phase and no beginning and no end.
Linear time-invariant systems process each sinusoidal component independently; the property of linearity means they satisfy the superposition principle.

Introduction

The group delay and phase delay properties of a linear time-invariant system are functions of frequency, giving the time from when a frequency component of a time varying physical quantity—for example a voltage signal—appears at the LTI system input, to the time when a copy of that same frequency component—perhaps of a different physical phenomenon—appears at the LTI system output.
A varying phase response as a function of frequency, from which group delay and phase delay can be calculated, typically occurs in devices such as microphones, amplifiers, loudspeakers, magnetic recorders, headphones, coaxial cables, and antialiasing filters. All frequency components of a signal are delayed when passed through such devices, or when propagating through space or a medium, such as air or water.
While a phase response describes phase shift in angular units, the phase delay is in units of time and equals the negative of the phase shift at each frequency divided by the value of that frequency. Group delay is the negative derivative of phase shift with respect to frequency.

Phase delay

A linear time-invariant system or device has a phase response property and a phase delay property, where one can be calculated exactly from the other. Phase delay directly measures the device or system time delay of individual sinusoidal frequency components. If the phase delay function at any given frequency—within a frequency range of interest—has the same constant of proportionality between the phase at a selected frequency and the selected frequency itself, the system/device will have the ideal of a flat phase delay property, a.k.a. linear phase. Since phase delay is a function of frequency giving time delay, a departure from the flatness of its function graph can reveal time delay differences among the signal’s various sinusoidal frequency components, in which case those differences will contribute to signal distortion, which is manifested as the output signal waveform shape being different from that of the input signal.
The phase delay property in general does not give useful information if the device input is a modulated signal. For that, group delay must be used.

Group delay

The group delay is a convenient measure of the linearity of the phase with respect to frequency in a modulation system. For a modulation signal, the information carried by the signal is carried exclusively in the wave envelope. Group delay therefore operates only with the frequency components derived from the envelope.

Basic modulation system

A device's group delay can be exactly calculated from the device's phase response, but not the other way around. The simplest use case for group delay is illustrated in Figure 1 which shows a conceptual modulation system, which is itself an LTI system with a baseband output that is ideally an accurate copy of the baseband signal input. This system as a whole is referred to here as the outer LTI system/device, which contains an inner LTI system/device. As is often the case for a radio system, the inner red LTI system in Fig 1 can represent two LTI systems in cascade, for example an amplifier driving a transmitting antenna at the sending end and the other an antenna and amplifier at the receiving end.

Amplitude Modulation

Amplitude modulation creates the passband signal by shifting the baseband frequency components to a much higher frequency range. Although the frequencies are different, the passband signal carries the same information as the baseband signal. The demodulator does the inverse, shifting the passband frequencies back down to the original baseband frequency range. Ideally, the output signal is a time delayed version of the input signal where the waveform shape of the output is identical to that of the input.
In Figure 1, the outer system phase delay is the meaningful performance metric. For amplitude modulation, the inner red LTI device group delay becomes the outer LTI device phase delay. If the inner red device group delay is completely flat in the frequency range of interest, the outer device will have the ideal of a phase delay that is also completely flat, where the contribution of distortion due to the outer LTI device's phase response—determined entirely by the inner device's possibly different phase response—is eliminated. In that case, the group delay of the inner red device and the phase delay of the outer device give the same time delay figure for the signal as a whole, from the baseband input to the baseband output. It is significant to note that it is possible for the inner device to have a very non-flat phase delay, while the outer device has the ideal of a perfectly flat phase delay. This is fortunate because in LTI device design, a flat group delay is easier to achieve than a flat phase delay.

Angle Modulation

In an angle-modulation system—such as with frequency modulation or phase modulation —the passband signal applied to an LTI system input can be analyzed as two separate passband signals, an in-phase amplitude modulation AM passband signal and a quadrature-phase amplitude modulation AM passband signal, where their sum exactly reconstructs the original angle-modulation passband signal. While the passband signal is not amplitude modulation, and therefore has no apparent outer envelope, the I and Q passband signals do indeed each have their own separate amplitude modulation envelopes. So, for each of the I and Q passband signals, a flat group delay ensures that neither the I pass band envelope nor the Q passband envelope will have wave shape distortion, so when the I passband signal and the Q passband signal are added back together, the sum is the original FM/PM passband signal, which will also be unaltered.

Theory

According to LTI system theory, the output signal of an LTI system can be determined by convolving the time-domain impulse response of the LTI system with the input signal. expresses this relationship as:
where denotes the convolution operation, and are the Laplace transforms of the input and impulse response, respectively, is the complex frequency, and is the inverse Laplace transform. is called the transfer function of the LTI system and, like the impulse response, fully defines the input-output characteristics of the LTI system. This convolution can be evaluated by using the integral expression in the time domain, or by using multiplication in the Laplace domain and then applying the inverse transform to return to time domain.

LTI system response to wave packet

Suppose that such a system is driven by a wave packet formed by a sinusoid multiplied by an amplitude envelope, so the input can be expressed in the following form:
Also suppose that the envelope is slowly changing relative to the sinusoid's frequency. This condition can be expressed mathematically as:
Applying the earlier convolution equation would reveal that the output of such an LTI system is very well approximated as:
Here is the group delay and is the phase delay, and they are given by the expressions below. The phase of the sinusoid, as indicated by the positions of the zero crossings, is delayed in time by an amount equal to the phase delay,. The envelope of the sinusoid is delayed in time by the group delay,.

Mathematical definition of group delay and phase delay

The group delay,, and phase delay,, are frequency-dependent and can be computed from the unwrapped phase shift. The phase delay at each frequency equals the negative of the phase shift at that frequency divided by the value of that frequency:
The group delay at each frequency equals the negative of the slope of the phase at that frequency:
In a linear phase system, both and are constant and equal, and their common value equals the overall delay of the system; and the unwrapped phase shift of the system is negative, with magnitude increasing linearly with frequency.

LTI system response to complex sinusoid

More generally, it can be shown that for an LTI system with transfer function driven by a complex sinusoid of unit amplitude,
the output is
where the phase shift is

1st order low- or high-pass RC filter example

The phase of a 1st-order low-pass filter formed by a RC circuit with cutoff frequency is:
Similarly, the phase for a 1st-order RC high-pass filter is:
Taking the negative derivative with respect to for either this low-pass or high-pass filter yields the same group delay of:
For frequencies significantly lower than the cutoff frequency, the phase response is approximately linear, so the group delay simplifies to a constant value of:
Similarly, right at the cutoff frequency,
As frequencies get even larger, the group delay decreases with the inverse square of the frequency and approaches zero as frequency approaches infinity.

Negative group delay

Filters will have negative group delay over frequency ranges where its phase response is positively-sloped. If a signal is band-limited within some maximum frequency B, then it is predictable to a small degree. A filter whose group delay is negative over that signal's entire frequency range is able to use the signal's predictability to provide an illusion of a non-causal time advance. However, if the signal contains an unpredictable event, then the illusion breaks down. Circuits with negative group delay are possible, though causality is not violated.
Negative group delay filters can be made in both digital and analog domains. Applications include compensating for the inherent delay of low-pass filters, to create zero phase filters, which can be used to quickly detect changes in the trends of sensor data or stock prices.

Group delay in audio

Group delay has some importance in the audio field and especially in the sound reproduction field.

Threshold of audibility

Many components of an audio reproduction chain, notably loudspeakers and multiway loudspeaker crossover networks, introduce group delay in the audio signal. It is therefore important to know the threshold of audibility of group delay with respect to frequency, especially if the audio chain is supposed to provide high fidelity reproduction. The best thresholds of audibility table has been provided by Blauert and Laws.
Frequency
Threshold
Periods
0.53.21.6
122
212
41.56
8216

Flanagan, Moore and Stone conclude that at 1, 2 and 4 kHz, a group delay of about 1.6 ms is audible with headphones in a non-reverberant condition. Other experimental results suggest that when the group delay in the frequency range from 300 Hz to 1 kHz is below 1.0 ms, it is inaudible.
It is possible to use digital signal processing techniques to correct the group delay distortion that arises due to the use of crossover networks in multi-way loudspeaker systems. This involves considerable computational modeling of loudspeaker systems in order to successfully apply delay equalization, using the Parks-McClellan FIR equiripple filter design algorithm.

Calculating minimum amplifier bandwidths for a specified phase linearity

As a criterion for calculating amplifier bandwidths for a specified phase linearity, Leach introduced the concept of differential time-delay distortion, defined as:
which is the difference between the phase delay and the group delay, where the phase delay refers to the time delay of a single relatively higher frequency sinusoid, and the group delay refers to the much lower frequencies derived from the envelope of the amplitude modulation applied to that sinusoid, resulting in what is perceptually a narrow band musical note. However, exhibiting zero or negligible differential time-delay distortion is a necessary but not sufficient property of an ideal system which must have flat phase delay. That's because when two narrow band amplitude modulated notes, one much higher up in the human audio range than the other so their frequency spectrums do not overlap, the results being presented separately are perceptionally ideal, because it's not likely that a human can detect that the begin time of one note took longer from device input to device output than that of the other note. However, in the same system, if both notes are presented together at the device input, then, at the device output, one note will take longer to begin than the other, and can alter the perception of the resulting combination note if the begin times at the output are separated too much in time.

Group delay in optics

Group delay is important in physics, and in particular in optics.
In an optical fiber, group delay is the transit time required for optical power, traveling at a given mode's group velocity, to travel a given distance. For optical fiber dispersion measurement purposes, the quantity of interest is group delay per unit length, which is the reciprocal of the group velocity of a particular mode. The measured group delay of a signal through an optical fiber exhibits a wavelength dependence due to the various dispersion mechanisms present in the fiber.
It is often desirable for the group delay to be constant across all frequencies; otherwise there is temporal smearing of the signal. Because group delay is, it therefore follows that a constant group delay can be achieved if the transfer function of the device or medium has a linear phase response. The degree of nonlinearity of the phase indicates the deviation of the group delay from a constant value.
The differential group delay is the difference in propagation time between the two eigenmodes X and Y polarizations. Consider two eigenmodes that are the 0° and 90° linear polarization states. If the state of polarization of the input signal is the linear state at 45° between the two eigenmodes, the input signal is divided equally into the two eigenmodes. The power of the transmitted signal ET,total is the combination of the transmitted signals of both x and y modes.
The differential group delay Dt is defined as the difference in propagation time between the eigenmodes: Dt = |tt,''xt''t,''y''|.

True time delay

A transmitting apparatus is said to have true time delay if the time delay is independent of the frequency of the electrical signal. TTD allows for a wide instantaneous signal bandwidth with virtually no signal distortion such as pulse broadening during pulsed operation.
TTD is an important characteristic of lossless and low-loss, dispersion free, transmission lines. reveals that signals propagate through them at a speed of for a distributed inductance and capacitance. Hence, any signal's propagation delay through the line simply equals the length of the line divided by this speed.

Group delay from transfer function polynomials

If a transfer function or Sij of a scattering parameter, is in a polynomial Laplace transform form, then the mathematical definition for group delay above may be solved analytically in closed form. A polynomial transfer function may be taken along the axis and defined as. may be determined from, and then the group delay may be determined by solving for.
to determine from, use the definition of. Given that is always real, and is always imaginary, may be redefined as where even and odd refer to the polynomials that contain only the even or odd order coefficients respectively. The in the numerator merely converts the imaginary numerator to a real value, since by itself is purely imaginary.
The above expressions contain four terms to calculate:
The equations above may be used to determine the group delay of polynomial in closed form, shown below after the equations have been reduced to a simplified form.

Polynomial ratio

A polynomial ratio of the form, such as that typically found in the definition of filter designs, may have the group delay determined by taking advantage of the phase relation,.

Simple filter example

A four pole Legendre filter transfer function used in the Legendre filter example is shown below.
The numerator group delay by inspection is zero, so only the denominator group delay need be determined.
Evaluating at = 1 rad/sec:
The group delay calculation procedure and results may be confirmed to be correct by comparing them to the results derived from the digital derivative of the phase angle,, using a small delta of +/-1.e-04 rad/sec.
Since the group delay calculated by the digital derivative using a small delta is within 7 digits of accuracy when compared to the precise analytical calculation, the group delay calculation procedure and results are confirmed to be correct.

Deviation from Linear Phase

Deviation from Linear Phase,, sometimes referred to as just, "phase deviation", is the difference between the phase response,, and the linear portion of the phase response, and is a useful measurement to determine the linearity of.
A convenient means to measure is to take the simple linear regression of sampled over a frequency range of interest, and subtract it from the actual. The of an ideal linear phase response would be expected to have a value of 0 across the frequency range of interest, while the of a real-world approximately linear phase response may deviate from 0 by a small finite amount across the frequency range of interest.

Advantage over group delay

An advantage of measuring or calculating over measuring or calculating group delay,, is always converges to 0 as the phase becomes linear, whereas converges on a finite quantity that may not be known ahead of time. Given this, a linear phase optimizing function may more easily be executed with a goal than with a goal when the value for is not necessarily already known.