# Derivation of Bessel filters

I was reading Winder's Analog and Digital Filter Design and the section on Bessel filter.

I was hoping to see a complete derivation of the Bessel filter theory, but Winder's book gives only

\begin{align} H(s)&=e^{-sT} \tag{$\scriptstyle\text{pure delay}$}\\ H(s)&=e^{-s}=\frac{1}{\sinh(s)+\cosh(s)} \tag{$\scriptstyle\text{normalized pure delay, T=1}$} \end{align}

and then using series expansion of $\sinh$ and $\cosh$ we're left with

$$H(s)=\frac{a_0}{B_n(s)}$$

where $B_n(s)$ is a Bessel polynomial.

Anyone know where can I get a more exhaustive derivation of the Bessel filters?

• I don't know of any explicit derivation, but the Wikipedia page Example section suggests what is happening: The idea is to derive the group delay $D(\omega)$ and set as many of the Taylor series terms to zero.
– Peter K.
Jun 27, 2016 at 12:13

The optimality criterion for Bessel (low pass) filters is that they have a maximally flat group delay at $$\omega=0$$, as correctly pointed out in a comment by Peter K.. The derivation is slightly tortuous, but you can find it for the continuous-time case in these lecture notes, which are taken from the book Introduction to Circuit Synthesis and Design by Temes and La Patra. I can recommend this book to anybody who wants to gain a deep understanding of electrical network theory.

The interesting fact is that unlike for other standard filters (Butterworth, Chebyshev, Cauer), the discrete-time version of a Bessel filter cannot be derived from its continuous-time counterpart via the bilinear transform. The reason is the frequency warping inherent in the bilinear transform, which doesn't affect the optimality conditions of the other standard filters (maximal flatness of the magnitude, or equiripple behavior of the magnitude, or a combination thereof), but it does not preserve the maximally flat group delay. A discrete-time filter derived from an analog Bessel filter via the bilinear transform will not be optimal in any way, and it will not have a maximally flat group delay. This means that a discrete-time Bessel filter must be derived directly in the discrete-time domain. This has been done in this famous paper by Thiran. Again, the derivation is slightly painful but I'm afraid there is no easier route.

What you can take away:

• Bessel filters have a maximally flat group delay.
• Discrete-time Bessel filters cannot be obtained from continuous-time Bessel filters via the bilinear (or any other) transform.
• Nice answer, Matt! :-)
– Peter K.
Jun 27, 2016 at 16:47
• @PeterK.: Thanks Peter, that stuff reminds me of my student times ... :) Jun 27, 2016 at 16:48
• Personally, I find the Temes-Lapatra method a bit more complicated, but that may be just me. This method (link on ee.se) seems a bit more to the point (but I won't say it's "better", or worse", just another way). Apr 2 at 20:27