Time and frequency response of a lowpass FIR filter designed with Parks-McClellan algorithm for N=33

Why FIR Filters have Linear Phase

One of the most attractive properties of a Finite Impulse Response (FIR) filter is that a linear phase response is easier to achieve. Not all FIR filters have linear phase though. This is only possible when the coefficients or taps of the filter are symmetric or anti-symmetric around a point.

Today I want to describe the reason behind this kind of phase response in an intuitive manner. We have described Finite Impulse Response (FIR) filters before. Moreover, we have also discussed that the Discrete Fourier Transform (DFT) of a signal is complex in general and therefore both magnitude and phase plots come into the picture.

Background


Based on whether the filter length is even or odd, and the filter taps are symmetric or anti-symmetric, there are 4 types of FIR filters. Our focus is on the simplest case: an odd length filter (so that the taps are symmetric around a single point) with symmetric coefficients.

Shown below in the figure is an odd-length FIR filter with symmetric coefficients. The corresponding magnitude response is also plotted below.

Time and frequency response of a lowpass FIR filter designed with Parks-McClellan algorithm for N=33

To understand the idea, let us consider only the two samples at times $n=+1$ and $n=-1$. Imagine that the remaining samples have disappeared and whatever conclusion we draw for the samples at $n=\pm 1$ can also be applied to them.

An Impulse in Time


What does a time domain impulse correspond to in frequency domain? To answer this question, we refer to the definition of the Discrete-Time Fourier Transform (DFTF).
\[
X(e^{j\omega})=\sum_{n=0}^{N-1} x[n]e^{-j\omega n}
\]

Clearly, when the input $x[n]$ is an impulse at $n=+1$, we have
\[
X(e^{j\omega})=\sum_{n=0}^{N-1} \delta [n-1]e^{-j\omega n} = e^{-j\omega \cdot 1}
\]

The result is a complex sinusoid in frequency domain with period 1! Another way of seeing this ideas is as follows.

Time-Frequency Duality


From the concept of frequency, we know that a complex sinusoid in time domain corresponds to an impulse in frequency domain. Owing to the time-frequency duality, an impulse in time domain corresponds to a complex sinusoid in frequency domain!

But we have two impulses in terms of two filter coefficients at $n=\pm 1$.

From One to Two Impulses


Two impulses at $n=\pm 1$ imply two complex sinusoids rotating in frequency domain in opposite directions. The result is that the frequency domain signal can be written as
\begin{equation}\label{equation-cosine-1}
e^{+j\omega\cdot 1}+e^{-j\omega\cdot 1} = 2\cos \left(\omega \cdot 1\right)
\end{equation}

where we have used the Euler’s identity to simplify the above expression. Observe that the above signal is not a time domain cosine; instead, it is a frequency domain sinusoid.

Symmetric Coefficients


With this background, symmetric coefficients of an FIR filter $h[n]$ come into play. For an $N$-length filter due to symmetry, we have
\[
\begin{aligned}
h[0] &= h[N-1] \\
h[1] &= h[N-2] \\
\vdots &= \vdots\\
\end{aligned}
\]

Since DTFT is given by $H(e^{j\omega})=\sum_{n=0}^{N-1} h[n]e^{-j\omega n}$, the contribution from these two terms can be written as
\[
h[0]e^{-j\omega\cdot 0} + h[N-1]e^{-j\omega(N-1)}
\]

This expression can be simplified using the following two facts.

  • $N$ is odd, so $N-1$ is an even number.
  • The coefficients $h[0]$ and $h[N-1]$ are the same.

Thus, we can write
\[
\begin{aligned}
h[0]e^{-j\omega\cdot 0} + h[0]e^{-j\omega(N-1)} &= e^{-j\omega\cdot \frac{N-1}{2}} h[0]\Big[e^{+j\omega\cdot \frac{N-1}{2}} + e^{-j\omega\cdot \frac{N-1}{2}}\Big]\\
&= e^{-j\omega\cdot \frac{N-1}{2}} \underbrace{ h[0]\cdot 2\cos \left(\omega \frac{N-1}{2}\right)}_{\text{Real frequency domain signal}}
\end{aligned}
\]

where we have used the Euler’s identity for cosine above. The real signal has no imaginary component and hence a zero phase. All that is left in the phase part is the following term.
\[
\text{Phase response} = -\omega\cdot \frac{N-1}{2}
\]

which is clearly a linear expression similar to $y=mx$. The x-axis represents $\omega$ while the slope is given by $-(N-1)/2$.

Conclusion


In the above equation, we have only considered the contributions from impulses at $n=\pm 1$. The impulses or coefficients at other times also add up to real cosines with similar phase contributions as above, i.e., $(N-1)/2$.

This is why it is straightforward to design an FIR filter with linear phase response by simply choosing symmetric or anti-symmetric coefficients.

Leave a Reply

Your email address will not be published. Required fields are marked *