Hello! Welcome back to our study of analog filter design.
In our last two lessons, we explored the Butterworth filter, which offers a "maximally flat" passband, and the Chebyshev filter, which sacrifices that flatness for a steeper transition by allowing ripple in either the passband (Type I) or the stopband (Type II). This leads to a natural question: what happens if we allow ripple in both bands?
This is precisely the idea behind the elliptic filter, also known as the Cauer filter after its inventor, the German mathematician Wilhelm Cauer. His work in the 1930s, building on the 19th-century mathematics of Yegor Zolotarev, led to the most efficient of all classical analog filters. For any given filter order, the elliptic filter provides the steepest possible transition from passband to stopband.
Today, our goal is to derive the elliptic (Cauer) filter approximation and analyze its equiripple behavior in both passbands and stopbands. By the end of this lesson, you will understand the mathematical foundation of these optimal filters and how they achieve their remarkable performance.
1. The Elliptic Filter's Magnitude Response
Just as Chebyshev polynomials were the key to Chebyshev filters, the magnitude response of an elliptic filter is defined by a special class of functions called elliptic rational functions.
The magnitude-squared response of an Nth-order low-pass elliptic filter is given by:
Let's break this down:
- is the Nth-order elliptic rational function. This function is the heart of the filter.
- is the ripple factor, a parameter that controls the amount of passband ripple, just like in the Chebyshev filter.
- (xi) is the selectivity factor, defined as the ratio of the stopband edge frequency to the passband edge frequency (). This parameter determines the steepness of the transition band.
The defining characteristic of the elliptic rational function is that it is a ratio of polynomials. This is a crucial difference from the Chebyshev polynomial . Because is a rational function, its transfer function will have both poles and zeros, just like a Chebyshev Type II filter. However, in the elliptic filter, these poles and zeros are placed in a highly specific way to create ripple in both the passband and stopband.
Let's examine the response visually.
The following resource provides a fantastic introduction to the elliptic filter, its history, and its key properties. It also includes helpful graphs of the frequency response.
Please read the introduction and the 'Properties' section of the Wikipedia article on Elliptic filters. Pay close attention to the formula for the gain G_n(ω) and the accompanying graphs. Note how the filter relates to Chebyshev and Butterworth filters in the limiting cases.
As you saw in the reading, the response looks like this:
![]()
The magnitude response of a 4th-order elliptic filter. Note the equal-height ripples in the passband and the equal-depth nulls (zeros) in the stopband, which define the equiripple characteristic in both bands.
Let's analyze this behavior:
- Passband (): The elliptic rational function oscillates in this region, causing the magnitude to ripple between a maximum of 1 and a minimum of . This is very similar to a Chebyshev Type I filter.
- Stopband (): The transfer function has zeros on the -axis, which cause the deep nulls you see in the plot. Between these nulls, the response rises to local peaks. The magic of the Cauer approximation is that the poles and zeros are placed so that these stopband peaks are all of equal height, creating an equiripple stopband.
Test your understanding!
Looking at the magnitude response formula, what happens to the stopband attenuation as the selectivity factor gets larger (i.e., the transition band becomes wider)?
Show answer
A larger means a wider transition band. For a fixed filter order , a wider transition band allows for much greater stopband attenuation. The elliptic rational function grows more rapidly for larger , leading to a much larger value in the denominator of the gain function in the stopband, and thus a smaller overall gain (i.e., more attenuation).
2. Deriving the Elliptic Filter Approximation
The design of an elliptic filter is mathematically intensive and relies on a family of special functions known as Jacobi elliptic functions (e.g., sn, cn) and elliptic integrals. Your background in advanced mathematics will be helpful in appreciating the elegance of this formulation.
The process, similar to other filter types, involves three main steps:
- Determine the minimum filter order to meet specifications.
- Calculate the locations of the transfer function's poles and zeros.
- Construct the final transfer function .
Let's delve into the mathematics.
The Mathematical Machinery
The elliptic rational functions are constructed from Jacobi elliptic functions. While a full derivation is beyond the scope of this lesson, it's important to recognize their role as the fundamental building blocks.
13.2 Analog Elliptic Filter Design
To understand the mathematical machinery behind elliptic filters, this document from PTC Mathcad provides the definitions of the core elliptic functions and integrals.
Please read the section 'Elliptic Function Definitions'. You don't need to memorize the integrals, but recognize that these special functions—like the complete elliptic integral K(k) and the Jacobi elliptic sine sn(u,k)—are the essential tools for this derivation, much like sine and cosine are for Fourier analysis.
Step 1: Determining the Filter Order
Given the desired passband ripple (), stopband attenuation (), passband edge (), and stopband edge (), the first task is to find the minimum filter order that can satisfy these constraints. This calculation directly involves the complete elliptic integral .
13.2 Analog Elliptic Filter Design
Now, let's see how these functions are used in practice. The following section of the Mathcad document demonstrates the calculation of the required filter order N from a set of design specifications.
Read the 'Calculation of Filter Order' section. Follow how the design specifications (ripple, band edges) are translated into intermediate values called moduli (k, k1) and then used in the formula for N. This is the heart of the filter approximation.
As you've just read, the formula for the required order is:
where , is related to the ripple specifications, and denotes the complementary modulus . This formula elegantly captures the relationship between the filter's specifications and the required complexity ().
Step 2: Finding the Poles and Zeros
With the order determined, the next step is to find the locations of the poles and the associated zeros.
- The zeros are located on the -axis, within the stopband.
- The poles are located in the left-half of the s-plane on an ellipse, similar to the Chebyshev filter.
![]()
Pole-zero plot for an 8th-order elliptic lowpass filter. The zeros lie on the imaginary axis, creating stopband nulls. The poles lie on an ellipse in the left-half plane, creating the passband ripple.
The precise locations are calculated using the Jacobi elliptic sine function, sn(u, k).
13.2 Analog Elliptic Filter Design
The calculation of the pole and zero locations is where the Jacobi elliptic sine function, sn(u,k), takes center stage. This next section provides the formulas and shows a resulting pole-zero plot.
Read the section 'Finding the Zeros and Poles of the Transfer Function'. Observe the formulas for the poles (sp) and zeros (sr) and note how they are derived from the elliptic functions. You don't need to perform the calculation, but understand the process.
Step 3: Constructing the Transfer Function
Once the left-half-plane poles and the imaginary-axis zeros are found, the transfer function is simply constructed as a rational function:
Here, is the number of zeros (for a low-pass filter, ), and is a normalization constant chosen to set the gain at DC to the appropriate value (either 1 or , depending on whether is odd or even).
The Mathcad resource you've been using (LINK, section 4) continues from here to construct the final transfer function and plot its frequency response, which provides a satisfying visual confirmation that the intricate mathematics correctly produces the desired equiripple behavior in both bands.
3. The Optimality of Elliptic Filters
The reason elliptic filters are so important is their optimality. For a given order and set of ripple specifications, no other linear time-invariant filter can achieve a faster rate of attenuation between the passband and stopband.
![]()
Comparison of 5th-order Butterworth, Chebyshev, and Elliptic filters. The elliptic filter clearly has the sharpest transition between the passband and stopband, achieved by tolerating ripple in both bands.
This graph perfectly illustrates the trade-offs we've been discussing throughout this module. The Butterworth is smooth but slow. The Chebyshev is faster but has passband ripple. The Elliptic filter is the fastest of all, representing the limit of what is possible for a given filter order.
Conclusion
In this lesson, we have explored the most efficient of the classical analog filter approximations, the elliptic (Cauer) filter. You have seen how it leverages sophisticated mathematics to achieve performance that is, in a very real sense, optimal.
Key Takeaways:
- Elliptic filters achieve the sharpest possible transition for a given order by featuring equiripple behavior in both the passband and the stopband.
- The filter's response is defined by elliptic rational functions, which are constructed from Jacobi elliptic functions and elliptic integrals.
- The transfer function has both poles (arranged on an ellipse in the LHP) and zeros (on the -axis), which are meticulously placed to create the dual-equiripple characteristic.
- The design process involves using elliptic integrals to calculate the required filter order and then using Jacobi elliptic functions to find the pole and zero locations to construct .
Preview of the Next Lesson:
We have now studied the "big three" analog filter approximations: Butterworth, Chebyshev, and Elliptic. Each embodies a different design philosophy and offers a unique set of trade-offs. In our next lesson, we will step back and systematically analyze the trade-offs between these filter types in terms of order, transition bandwidth, and ripple, helping you understand how to select the best filter for a given engineering problem.