Wave Walker DSP

DSP Algorithms for RF Systems

Most Popular

Buy the Book!

DSP for Beginners: Simple Explanations for Complex Numbers! The second edition includes a new chapter on complex sinusoids.

Deriving the Ideal Low Pass Filter (LPF)
October 1, 2022

Table of Contents

Introduction

This blog derives the ideal low pass filter (LPF) step by step using simple math.

More DSP blogs!

Frequency Domain Ideal Low Pass Filter

An ideal low pass filter (LPF) allows signals below the cutoff frequency \omega_c to pass through unmodified with a linear gain of 1 and no phase change, while simultaneously completely rejecting all frequencies above the cutoff frequency with a linear gain of 0.

Mathematically, the ideal LPF is defined in the frequency domain by [oppenheim1999, p.43]

(1)   \begin{equation*}H\left(e^{j\omega}\right) = \begin{cases}1, & -\omega_c \le \omega \le \omega_c, \\0, & \text{otherwise}.\end{cases}\end{equation*}

The frequency response for the ideal LPF in (1) for \omega_c = \pi/8 is given in Figure 1.

Figure 1: The frequency response of the ideal low pass filter (LPF) for omega c = 1/8.
Figure 1: The frequency response of the ideal low pass filter (LPF) for omega c = 1/8.

Time Domain Ideal Low Pass Filter

The ideal LPF in the time domain is derived by taking the inverse discrete-time Fourier transform (IDTFT) of (1),

(2)   \begin{equation*}h[n] = \mathcal{F}^{-1} \left\{ H\left( e^{j\omega} \right) \right\}.\end{equation*}

The IDTFT is defined as

(3)   \begin{equation*}x[n] = \frac{1}{2\pi} \int_{-\pi}^{\pi} X\left( e^{j\omega} \right) e^{j\omega n}d\omega.\end{equation*}

Substituting (1) into (3),

(4)   \begin{equation*}h[n] = \frac{1}{2\pi} \int_{-\omega_c}^{\omega_c} e^{j\omega n} d\omega.\end{equation*}

The integral in (4) is simplified via Euler’s formula as

(5)   \begin{equation*}\begin{split}h[n] & = \frac{1}{j2\pi n} e^{j \omega n} \Big|_{-\omega_c}^{\omega_c} \\& = \frac{1}{j2\pi n} \left( e^{j\omega_c n} - e^{-j\omega_c n} \right) \\& = \frac{1}{\pi n} \sin \left( \omega_c n \right). \\\end{split}\end{equation*}

The cutoff frequency \omega_c is the frequency in radians, defined as

(6)   \begin{equation*}\omega_c = 2\pi f_{n},\end{equation*}

where f_{n} is the frequency f_c in Hz normalized by the sampling rate f_s,

(7)   \begin{equation*}f_{n} = \frac{f_c}{f_s}.\end{equation*}

The impulse response h[n] can therefore be written as

(8)   \begin{equation*}h[n] = \frac{1}{\pi n} \sin \left( 2\pi f_{n} n \right).\end{equation*}

Equation (8) can then be transformed into a sinc function where

(9)   \begin{equation*}\text{sinc} \left( x \right) = \frac{ \sin \left( \pi x \right) }{\pi x}.\end{equation*}

The impulse response h[n] can therefore be written as

(10)   \begin{equation*}\begin{split}h[n] & = \frac{1}{\pi n} \sin \left( 2 \pi f_n n \right) \\& = \frac{2 f_n}{2 f_n \pi n} \left( 2 \pi f_n n \right) \\& = 2 f_n \text{sinc} \left( 2 f_n n \right) \\\end{split}\end{equation*}

The impulse response from (10) is given in Figure 2.

Figure 2: The impulse response for the ideal low pass filter (LPF) for fn = 1/8.
Figure 2: The impulse response for the ideal low pass filter (LPF) for fn = 1/8.

Conclusion

The ideal low pass filter is a finite-length rectangle in the frequency domain but an infinitely long sinc in the time domain. Selecting the cutoff frequency determines which frequencies will be passed by the filter and which will be attenuated.

More DSP blogs!

God, the Lord, is my strength; he makes my feet like the deer’s; he makes me tread on my high places. Habakkuk 3:19

Leave a Reply

This website participates in the Amazon Associates program. As an Amazon Associate I earn from qualifying purchases.