View source for Elliptic filter

{{Linear analog electronic filter|filter2=hide|filter3=hide}}

An '''elliptic filter''' (also known as a '''Cauer filter''', named after [[Wilhelm Cauer]]) is a [[filter (signal processing)|signal processing filter]] with equalized [[ripple (filters)|ripple]] ([[equiripple]]) behavior in both the [[passband]] and the [[stopband]]. The amount of ripple in each band is independently adjustable, and no other filter of equal order can have a faster transition in [[gain]] between the [[passband]] and the [[stopband]], for the given values of ripple (whether the ripple is equalized or not). Alternatively, one may give up the ability to independently adjust the passband and stopband ripple, and instead design a filter which is maximally insensitive to component variations.

As the ripple in the stopband approaches zero, the filter becomes a type I [[Chebyshev filter]]. As the ripple in the passband approaches zero, the filter becomes a type II [[Chebyshev filter]] and finally, as both ripple values approach zero, the filter becomes a [[Butterworth filter]].

The gain of a [[lowpass]] elliptic filter as a function of angular frequency ω is given by:

:<math>
G_n(\omega) = {1 \over \sqrt{1 + \epsilon^2 R_n^2(\xi,\omega/\omega_0)}}
</math>

where R<sub>n</sub> is the ''n''th-order [[Elliptic rational functions|elliptic rational function]] (sometimes known as a Chebyshev rational function) and

:<math>\omega_0</math>  is the cutoff frequency
:<math>\epsilon</math> is the ripple factor
:<math>\xi</math> is the selectivity factor

The value of the ripple factor specifies the passband ripple, while the combination of the ripple factor and the selectivity factor specify the stopband ripple.

== Properties ==

[[File:Elliptic Filter Response (4th Order).svg|thumb|right|350px|The frequency response of a fourth-order elliptic low-pass filter with &epsilon;=0.5 and &xi;=1.05. Also shown are the minimum gain in the passband and the maximum gain in the stopband, and the transition region between normalized frequency 1 and &xi;]]
[[File:Elliptic Filter Response (4th Order, Detail).svg|thumb|right|340px|A closeup of the transition region of the above plot.]]

* In the passband, the elliptic rational function varies between zero and unity. The passband of the gain therefore will vary between 1 and <math>1/\sqrt{1+\epsilon^2}</math>.

* In the stopband, the elliptic rational function varies between infinity and the discrimination factor <math>L_n</math> which is defined as:

:<math>L_n=R_n(\xi,\xi)\,</math>

:The gain of the stopband therefore will vary between 0 and <math>1/\sqrt{1+\epsilon^2L_n^2}</math>.

* In the limit of <math>\xi \rightarrow \infty</math> the elliptic rational function becomes a [[Chebyshev polynomial]], and therefore the filter becomes a [[Chebyshev filter|Chebyshev type I filter]], with ripple factor ε

* Since the Butterworth filter is a limiting form of the Chebyshev filter, it follows that in the limit of <math>\xi \rightarrow \infty</math>, <math>\omega_0 \rightarrow 0</math> and <math>\epsilon \rightarrow 0</math> such that <math>\epsilon\,R_n(\xi,1/\omega_0)=1</math> the filter becomes a [[Butterworth filter]]

* In the limit of <math>\xi \rightarrow \infty</math>, <math>\epsilon \rightarrow 0</math> and  <math>\omega_0\rightarrow 0</math> such that <math>\xi\omega_0=1</math> and <math>\epsilon L_n=\alpha</math>, the filter becomes a [[Chebyshev filter|Chebyshev type II filter]] with gain

::<math>G(\omega)=\frac{1}{\sqrt{1+\frac{1}{\alpha^2 T^2_n(1/\omega)}}}</math>
<br style="clear:both;" />

== Poles and zeroes ==

[[File:Elliptic Filter s-Plane Response (8th Order).svg|right|thumb|300px|Log of the absolute value of the gain of an 8th order elliptic filter in [[complex frequency space]] (s=&sigma;+j&omega;) with &epsilon;=0.5, &xi;=1.05 and <math>\omega_0=1</math>. The white spots are poles and the black spots are zeroes. There are a total of 16 poles and 8 double zeroes. What appears to be a single pole and zero near the transition region is actually four poles and two double zeroes as shown in the expanded view below. In this image, black corresponds to a gain of 0.0001 or less and white corresponds to a gain of 10 or more.]]
[[File:Elliptic Filter s-Plane Response (8th Order, Detail).svg|right|thumb|300px|An expanded view in the transition region of the above image, resolving the four poles and two double zeroes.]]

The zeroes of the gain of an elliptic filter will coincide with the poles of the elliptic rational function, which are derived in the article on [[elliptic rational functions]].

The poles of the gain of an elliptic filter may be derived in a manner very similar to the derivation of the poles of the gain of a type I [[Chebyshev filter]]. For simplicity, assume that the cutoff frequency is equal to unity. The poles <math>(\omega_{pm})</math> of the gain of the elliptical filter will be the zeroes of the denominator of the gain. Using the complex frequency <math>s=\sigma+j\omega</math> this means that:

:<math>1+\epsilon^2R_n^2(-js,\xi)=0\,</math>

Defining <math>-js=\mathrm{cd}(w,1/\xi)</math> where cd() is the [[Jacobi elliptic functions|Jacobi elliptic cosine function]] and using the definition of the elliptic rational functions yields:

:<math>1+\epsilon^2\mathrm{cd}^2\left(\frac{nwK_n}{K},\frac{1}{L_n}\right)=0\,</math>

where <math>K=K(1/\xi)</math> and <math>K_n=K(1/L_n)</math>. Solving for ''w''

:<math>w=\frac{K}{nK_n}\mathrm{cd}^{-1}\left(\frac{\pm j}{\epsilon},\frac{1}{L_n}\right)+\frac{mK}{n}</math>

where the multiple values of the inverse cd() function are made explicit using the integer index ''m''.

The poles of the elliptic gain function are then:

:<math>s_{pm}=i\,\mathrm{cd}(w,1/\xi)\,</math>

As is the case for the Chebyshev polynomials, this may be expressed in explicitly complex form {{harv|Lutovac|et al.|2001|loc=&sect; 12.8}}

:<math>s_{pm}=\frac{a+jb}{c}</math>

:<math>a=-\zeta_n\sqrt{1-\zeta_n^2}\sqrt{1-x_m^2}\sqrt{1-x_m^2/\xi^2}</math>

:<math>b=x_m\sqrt{1-\zeta_n^2(1-1/\xi^2)}</math>

:<math>c=1-\zeta_n^2+x_i^2\zeta_n^2/\xi^2</math>

where <math>\zeta_n</math> is a function of <math>n,\,\epsilon</math> and <math>\xi</math> and <math>x_m</math> are the zeroes of the elliptic rational function. <math>\zeta_n</math> is expressible for all ''n'' in terms of Jacobi elliptic functions, or algebraically for some orders, especially orders 1,2, and 3. For orders 1 and 2 we have

:<math>\zeta_1=\frac{1}{\sqrt{1+\epsilon^2}}</math>

:<math>\zeta_2=\frac{2}{(1+t)\sqrt{1+\epsilon^2}+\sqrt{(1-t)^2+\epsilon^2(1+t)^2}}</math>

where 

:<math>t=\sqrt{1-1/\xi^2}</math>

The algebraic expression for <math>\zeta_3</math> is rather involved (See {{harvtxt|Lutovac|et al.|2001|loc=&sect; 12.8.1}}).

The nesting property of the [[elliptic rational functions]] can be used to build up higher order expressions for <math>\zeta_n</math>:

:<math>\zeta_{m\cdot n}(\xi,\epsilon)=
\zeta_m\left(\xi,\sqrt{\frac{1}{\zeta_n^2(L_m,\epsilon)}-1}\right)</math>

where <math>L_m=R_m(\xi,\xi)</math>.

== Minimum Q-factor elliptic filters ==
[[Image:Elliptic8 Qfactor.png|300px|thumb|right|
The normalized Q-factors of the poles of an 8-th order elliptic filter with ξ=1.1 as a function of ripple factor ε. Each curve represents four poles, since complex conjugate pole pairs and positive-negative pole pairs have the same Q-factor. (The blue and cyan curves nearly coincide). The Q-factor of all poles are simultaneously minimized at ε<sub>Qmin</sub>=1/√L<sub>n</sub>=0.02323...]]

See {{harvtxt|Lutovac|et al.|2001|loc=&sect; 12.11, 13.14}}.

Elliptic filters are generally specified by requiring a particular value for the passband ripple, stopband ripple and the sharpness of the cutoff. This will generally specify a minimum value of the filter order which must be used. Another design consideration is the sensitivity of the gain function to the values of the electronic components used to build the filter. This sensitivity is inversely proportional to the quality factor ([[Q-factor]]) of the poles of the transfer function of the filter. The Q-factor of a pole is defined as:

:<math>Q =-\frac{|s_{pm}|}{2\mathrm{Re} (s_{pm})} = -\frac{1}{2\cos(\arg(s_{pm}))}</math>

and is a measure of the influence of the pole on the gain function. For an elliptic filter, it happens that, for a given order, there exists a relationship between the ripple factor and selectivity factor which simultaneously minimizes the Q-factor of all poles in the transfer function:

:<math>\epsilon_{Qmin}=\frac{1}{\sqrt{L_n(\xi)}}</math>

This results in a filter which is maximally insensitive to component variations, but the ability to independently specify the passband and stopband ripples will be lost. For such filters, as the order increases, the ripple in both bands will decrease and the rate of cutoff will increase. If one decides to use a minimum-Q elliptic filter in order to achieve a particular minimum ripple in the filter bands along with a particular rate of cutoff, the order needed will generally be greater than the order one would otherwise need without the minimum-Q restriction. An image of the absolute value of the gain will look very much like the image in the previous section, except that the poles are arranged in a circle rather than an ellipse. They will not be evenly spaced and there will be zeroes on the ω axis, unlike the [[Butterworth filter]], whose poles are also arranged in a circle.

==Comparison with other linear filters==

Here is an image showing the elliptic filter next to other common kind of filters obtained with the same number of coefficients:

[[Image:Electronic linear filters.svg|500px|center]]

As is clear from the image, elliptic filters are sharper than all the others, but they show ripples on the whole bandwidth.

== References ==

* {{cite book |last=Daniels |first=Richard W. |authorlink= |coauthors= |title=Approximation Methods for Electronic Filter Design |year=1974 |publisher=McGraw-Hill |location=New York |isbn=0-07-015308-6 }}

* {{cite book |last=Lutovac |first=Miroslav D. |coauthors= Tosic, Dejan V., Evans, Brian L. |title=Filter Design for Signal Processing using MATLAB&copy; and Mathematica&copy; |year=2001 |publisher=Prentice Hall |location=New Jersey, USA |isbn=0-201-36130-2|ref=CITEREFLutovac_et_al.}}

{{DEFAULTSORT:Elliptic Filter}}
[[Category:Linear filters]]
[[Category:Network synthesis filters]]
[[Category:Electronic design]]

[[ca:Filtre el·líptic]]
[[de:Cauer-Filter]]
[[es:Filtro de Cauer]]
[[eo:Elipsa filtrilo]]
[[fr:Filtre elliptique]]
[[it:Filtro ellittico]]
[[ja:楕円フィルタ]]
[[pt:Filtro elíptico]]
[[ru:Эллиптический фильтр]]
[[sv:Cauerfilter]]
[[zh:椭圆函数滤波器]]
[[Category:2Fix]]