3.8  Importance of Linear Phase (or Constant Group Delay)

In some applications such as analog signal transmission and audio amplification, the system should ideally have an output y(t) identical to its input x(t). Given that obtaining y(t) = x(t) is often unfeasible due to propagation delays inherent of communication channels and electronic systems, a more realistic target is to obtain y(t) = x(t t0), a delayed version of the input. From the time-shift Fourier property (see Appendix A.17), the delay by t0 corresponds to a linear phase ej2πft0 in frequency domain. Hence, having a linear phase is an important property of a system to achieve distortionless transmission, i. e., letting a signal pass without distortion.

A linear phase filter with frequency response H(ω) has a phase Θ = ∠H(ω) that corresponds to a line segment in the passband (the phase behavior in the stopband is considered irrelevant). In practice, the line has a negative slope because a positive slope would correspond to a non-causal behavior (the output would be an anticipated version x(t + t0) of the input signal). Hence, it is convenient to define the group delay as:

τg(ω) = .
(3.44)

Similarly, the discrete-time version is

τg(ejΩ) = dΩ
(3.45)

with Θ = ∠H(ejΩ).

Example 3.8. Impact of the system phase on the output signal. To help understanding the practical importance of a system with linear phase, it is adopted a periodic pulse x[n] with N1 = 15 (number of non-zero samples in a period) and period N = 50 (see Example 2.12). The experiment is to obtain Y [k] by multiplying the DTFS X[k] of x[n] by ej2πN0kN, where N0 = 4. Listing 3.10 carries out the operation, including the conversion of Y [k] to y[n].

Listing 3.10: MatlabOctaveCodeSnippets/snip_systems_linear_phase_system.m. [ Python version]
1%Specify: N-period, N1/N-duty cicle, N0-delay, k-frequency 
2N=50; N1=15; N0=4; k=0:N-1; 
3xn=[ones(1,N1) zeros(1,N-N1)]; %x[n] 
4Xk=fft(xn)/N; %calculate the DTFS of x[n] 
5phase = -2*pi/N*N0*k; %define linear phase 
6Yk=Xk.*exp(j*phase); %impose the linear phase 
7yn=ifft(Yk)*N; %recover signal in time domain
  

Figure 3.24 illustrates the result: y[n] = x[n N0] is a perfect delayed version of x[n]. The middle plot shows the phase that was added to the original phase of X[k]. The magnitude of X[k] was left unchanged but that would not be enough to keep the shape of x[n] in case the phase was not linear with frequency.

PIC

Figure 3.24: Effect of adding a linear phase ej2πN0kN (N0 = 4 and N = 50) to the DTFS of x[n] resulting in a delayed version y[n] = x[n 4] .

To be contrasted with the linear phase case, Figure 3.25 provides an example where the phase is 0, 2 or 2 rad, as depicted in the top-most plot. In order to assure that y[n] is real, the phase is an odd function to preserve the Hermitian symmetry.

PIC

Figure 3.25: Effect of adding the specified nonlinear phase (top) to the pulse in Figure 3.24, which leads to a distorted signal (bottom).

Note that in the case of Figure 3.25, the pulses in x[n] are severely distorted due to the effect of a non-linear phase.    

Considering continuous-time (same is valid for discrete-time), a linear phase ej2πft0 avoids distorting an input signal x(t) because the system delays all components of x(t) by the same time interval t0. To obtain that, a linear phase system adds to the input signal component with frequency fi, a phase Θi = 2πfit0 that is proportional to fi (the larger fi, the larger the phase Θi added by the system). The following example concerns this aspect.

Example 3.9. A linear phase allows components with distinct frequencies to be delayed by the same angle. For example, assume components corresponding to ω = 200 and 400 rad/s with t0 = 4 s. The phases Θ = ωt0 are 800 and 1600 rad, respectively. On the other hand, if these two components pass a system that adds a fixed phase Θ = 1200 rad to both, the delays would be 6 and 3 s, respectively, which would potentially distort the delayed version with respect to the original signal x(t).    

For non-linear phase filters, the group delay must be interpreted as an average delay that is frequency-dependent (compare Figure 3.50 and Figure 3.51).11

Example 3.10. Gaussian filter and minimum group delay. In some applications, such as digital communications, minimizing the group delay is important to achieve low latency. In this case, the Gaussian filter is competitive because it has the minimum possible group delay. As a consequence, this filter has no overshoot to an input u(t) while minimizing the rise and fall time intervals. Its non-causal (and ideal, which requires truncation in practice) impulse response is

h(t) = π α e (πt α )2 ,
(3.46)

where α controls the 3-dB bandwidth BW. More specifically, α depends on the product of BW and the symbol period Tsym as follows:

α = 1 BWTsym ln 2 2 .
(3.47)

The frequency response of this filter is H(f) = eα2f2 and decreases with f but is not zero for finite f.