DSP Algorithms for RF Systems

Trending

Buy the Book!

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

Single Pole IIR Filter: Substitute for Moving Average Filter
August 10, 2022

Table of Contents

Introduction

The single pole IIR filter is incredibly efficient and can be a great stand-in for a moving average filter! In this blog I describe how to derive the 3 dB bandwidth of the IIR filter to be equivalent to a moving average filter by designing the \alpha weight. The design process builds on the IIR analysis from this DSP Stack Exchange question and my previous blog post.

More blog posts on DSP:

Single Pole IIR

The impulse response of the single pole IIR filter is

(1)   \begin{equation*}y[n] = \alpha x[n] + (1-\alpha) y[n-1]\end{equation*}

which can also be written as

(2)   \begin{equation*}y[n] = \alpha x[n] + \beta y[n-1]\end{equation*}

where 

(3)   \begin{equation*}\beta = 1-\alpha.\end{equation*}

The magnitude of the frequency response of a single pole IIR filter is given by

(4)   \begin{equation*}\left|H\left(e^{j\omega}\right)\right| = \frac{\alpha}{\sqrt{\left(1 - \beta e^{-j\omega}\right)\left(1 - \beta e^{j\omega}\right)}},\end{equation*}

where the magnitude-squared is

(5)   \begin{equation*}\left|H\left(e^{j\omega}\right)\right|^2 = \frac{\alpha^2}{\left(1 - \beta e^{-j\omega}\right)\left(1 - \beta e^{j\omega}\right)}.\end{equation*}

Setting up the Design Problem

A previous blog post demonstrated how the 3 dB bandwidth of a moving average filter of length N can be approximated by 

(6)   \begin{equation*}\omega_{3dB} = \pi/N.\end{equation*}

The goal is to find the value \alpha from (1) such that the 3 dB cutoff frequency of (5) is \pi/N.

IIR Filter Design Equation

The 3 dB cutoff corresponds to a magnitude-squared of 1/2. The design equation is therefore

(7)   \begin{equation*}\frac{\alpha^2}{\left(1 - \beta e^{-j\omega}\right)\left(1 - \beta e^{j\omega}\right)} = \frac{1}{2}.\end{equation*}

Multiplying the complex exponential terms,

(8)   \begin{equation*}\frac{\alpha^2}{1 + \beta^2 - \beta e^{j\omega} - \beta e^{-j\omega}} = \frac{1}{2},\end{equation*}

which is simplified as

(9)   \begin{equation*}\frac{\alpha^2}{1 + \beta^2 - \beta \left( e^{j\omega} + e^{-j\omega}\right)} = \frac{1}{2}.\end{equation*}

The complex exponentials can be written in terms of cosine,

(10)   \begin{equation*}2\cos\left( \omega \right) = e^{j\omega} + e^{-j\omega}.\end{equation*}

Substituting (10) into (9),

(11)   \begin{equation*}\frac{\alpha^2}{1 + \beta^2 - 2\beta \cos\left( \omega \right)} = \frac{1}{2}.\end{equation*}

Multiplying the terms of (11) results in

(12)   \begin{equation*}\begin{split}2\alpha^2 & = 1 + \beta^2 - 2 \beta \cos(\omega)\\0 & = - 2\alpha^2+ 1 + \beta^2 - 2 \beta \cos(\omega)\\\end{split}\end{equation*}

Substituting \beta = 1-\alpha,

(13)   \begin{equation*}\begin{split}0 & = -2\alpha^2 +1 + (1-\alpha)^2 - 2 (1-\alpha)\cos\left( \omega \right) \\0 & = -2\alpha^2 +1 + (1+\alpha^2-2\alpha) - 2 (1-\alpha)\cos\left( \omega \right) \\0 & = -\alpha^2 + 2\alpha\left(\cos\left( \omega \right)-1\right) + 2\left(1-\cos\left( \omega \right) \right) \\0 & = \alpha^2 + 2\alpha\left( 1 - \cos\left( \omega \right)\right) -2\left(1-\cos\left(\omega\right)\right).\end{split}\end{equation*}

Solving with the Quadratic Equation

The design equation (13) can be solved for \alpha using the quadratic equation,

(14)   \begin{equation*}\alpha = \frac{- b \pm \sqrt{b^2-4ac}}{2a}\end{equation*}

where

(15)   \begin{equation*}a = 1,\end{equation*}

(16)   \begin{equation*}b = 2\left( 1-\cos\left( \omega \right)\right),\end{equation*}

(17)   \begin{equation*}c = -2\left(1-\cos\left( \omega \right)\right).\end{equation*}

Substituting (15), (16) and (17) into (14),

(18)   \begin{equation*}\alpha = \frac{- 2\left( 1-\cos\left( \omega \right)\right) \pm \sqrt{\left(2\left( 1-\cos\left( \omega \right)\right)\right)^2-4\left(-2\left(1-\cos\left( \omega \right)\right)\right)}}{2}.\end{equation*}

Simplifying (18),

(19)   \begin{equation*}\begin{split}\alpha & = \frac{- 2\left( 1-\cos\left( \omega \right)\right) \pm \sqrt{4\left( 1-\cos\left( \omega \right)\right)^2 +8\left(1-\cos\left( \omega \right)\right)}}{2} \\& = \frac{- 2\left( 1-\cos\left( \omega \right)\right) \pm 2\sqrt{\left( 1-\cos\left( \omega \right)\right)^2 +2\left(1-\cos\left( \omega \right)\right)}}{2} \\& = -\left( 1-\cos\left( \omega \right)\right) \pm \sqrt{\left( 1-\cos\left( \omega \right)\right)^2 +2\left(1-\cos\left( \omega \right)\right)} \\& = \cos\left( \omega \right)-1 \pm \sqrt{1-2\cos\left( \omega \right)+\cos^2\left(\omega\right) +2-2\cos\left( \omega \right)} \\& = \cos\left( \omega \right)-1 \pm \sqrt{3-4\cos\left( \omega \right)+\cos^2\left(\omega\right)} \\\end{split}\end{equation*}

Picking the Alpha Value

The design equation in (19) has two solutions due to the \pm term. Figure 1 shows the two solutions for 0 \le \omega \le \pi.

Figure 1: The two solutions of the quadratic equation for all values of omega.
Figure 1: The two solutions of the quadratic equation for all values of omega.

The + from the \pm is chosen from (19) to select the positive coefficient, therefore the final design equation is

(20)   \begin{equation*}\alpha = \cos\left( \omega \right)-1 + \sqrt{3-4\cos\left( \omega \right)+\cos^2\left(\omega\right)}.\end{equation*}

To design a high pass filter (HPF), choose the – term in the \pm from (19).

Approximating the Moving Average Filter

The 3 dB cutoff for a length N moving average filter from (6) is substituted into the design equation (20) and the magnitude of the frequency response is plotted in Figure 2 for multiple values of N.

Figure 2: The frequency response of single pole IIR filters with different cutoff frequencies. The magnitudes are normalized to a maximum gain of 0 dB.
Figure 2: The frequency response of single pole IIR filters with different cutoff frequencies. The magnitudes are normalized to a maximum gain of 0 dB.

Comparison to Moving Average Filter

Figure 2 shows that increasing the filter cutoff frequency also reduces the sidelobe attenuation. Although the single pole IIR filter uses less computation than a moving average filter, it also has lower attenuation of out of band frequency content. Figure 3 shows that the single pole IIR filter approaches the response  of the moving average filter for N=8 but with about 4 dB worse attenuation.

Figure 3: Comparing a moving average filter N=8 with the single pole IIR filter approximation.
Figure 3: Comparing a moving average filter N=8 with the single pole IIR filter approximation.

Conclusion

The single pole IIR filter is a great substitute to the moving average filter! The substitution trades slightly worse sidelobe attenuation for a decrease in computation needed to implement the filter. This blog post described how to design the single pole IIR filter weight to approximate the frequency response of a moving average filter.

More blog posts on DSP:

2 Responses

  1. There is an excellent paper “Low-pass filters for signal averaging” by
    Edward Voigtman and James D. Winefordner. I provide the citation link below.

    In it, they explain how the ideal filter for this application is the integrator, one embodiment of which is the moving average together with a proper antialiasing filter. The single-pole IIR filter unfortunately does a terrible job at the proper metric for this application which is the product of settling time x noise bandwidth. In other words, tracking of sidelobes doesn’t matter when you’re averaging; what matters is how fast it settles to the right value and how much noise it lets through.

    If simplicity over performance is the target while still retaining low memory usage, you will do much much better by cascading three to six identical IIR single-lowpass of the type described in this article and adjusting the corner to meet your needs. It’s called a synchronous filter (hard to find in filter literature since it’s not a good filter for anything but averaging time series) and will approximate a Bessel filter for this data-averaging application.

    Citation: 57, 957 (1986); doi: 10.1063/1.1138645
    View online: http://dx.doi.org/10.1063/1.1138645
    View Table of Contents: http://aip.scitation.org/toc/rsi/57/5 Published by the American Institute of Physics

Leave a Reply

For everything there is a season, and a time for every matter under heaven. A time to cast away stones, and a time to gather stones together. A time to embrace, and a time to refrain from embracing. Ecclesiastes 3:1,5
The earth was without form and void, and darkness was over the face of the deep. And the Spirit of God was hovering over the face of the waters. Genesis 1:2
Behold, I am toward God as you are; I too was pinched off from a piece of clay. Job 33:6
Enter His gates with thanksgiving, and His courts with praise! Give thanks to Him; bless His name! Psalm 100:4
Lift up your hands to the holy place and bless the Lord! Psalm 134:2
Blessed is the man who trusts in the Lord, whose trust is the Lord. He is like a tree planted by water, that sends out its roots by the stream, and does not fear when heat comes, for its leaves remain green, and is not anxious in the year of drought, for it does not cease to bear fruit. Jeremiah 17:7-8
He said to him, “You shall love the Lord your God with all your heart and with all your soul and with all your mind. This is the great and first commandment. And a second is like it: You shall love your neighbor as yourself. On these two commandments depend all the Law and the Prophets.” Matthew 22:37-39
Then He said to me, “Prophesy over these bones, and say to them, O dry bones, hear the word of the Lord. Thus says the Lord God to these bones: Behold, I will cause breath to enter you, and you shall live." Ezekiel 37:4-5
Riches do not profit in the day of wrath, but righteousness delivers from death. Proverbs 11:4
The angel of the Lord appeared to him in a flame of fire out of the midst of a bush. He looked, and behold, the bush was burning, yet it was not consumed. And Moses said, “I will turn aside to see this great sight, why the bush is not burned.” When the Lord saw that he turned aside to see, God called to him out of the bush, “Moses, Moses!” And he said, “Here I am.” Exodus 3:2-3
Daniel answered and said: “Blessed be the name of God forever and ever, to whom belong wisdom and might. He changes times and seasons; He removes kings and sets up kings; He gives wisdom to the wise and knowledge to those who have understanding." Daniel 2:20-21
Now the Lord is the Spirit, and where the Spirit of the Lord is, there is freedom. 2 Corinthians 3:17
Previous slide
Next slide

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

© 2021-2024 Wave Walker DSP