Chebychev filter (Type I)

The Chebychev filter of type I is given by the normalized transfer function

$$H(S)=\frac{1}{2^{n-1}\epsilon}\prod_{k=1}^{n}\frac{1}{S-s_k}$$

where n is the order of the filter, ε defines the size of the ripples of the filter's magnitude response in the pass band, and the roots of the denominator sk are

$$s_k=-\sinh(\frac{\mathrm{asinh}(\frac{1}{\epsilon})}{n})\,\sin(\frac{\pi(2k-1)}{2n})$$ $$-j\,\cosh(\frac{\mathrm{asinh}(\frac{1}{\epsilon})}{n})\,\cos(\frac{\pi(2k-1)}{2n})$$

Here sinh and cosh are the hyperbolic sine and cosine functions and asinh and acosh are the inverse hyperbolic sine and cosine functions.

The transfer function above is shown irrespective of the cutoff frequency (i.e., for the cutoff frequency 1) and the various forms of the Chebychev type I filter require substitutions. This transfer function is normalized. Use S = s / ωc to get the Chebychev type I low pass filter, where s = jω, ωc is the cutoff frequency of the filter, and ω is the angular frequency spanning the frequency spectrum between 0 and π. Use S = ωc / s to get the Chebychev type I high pass filter. Use S = (s2 + ωc2) / (B s) to get the Chebychev type I band pass filter, where ωc is the midpoint of the pass band and B is the width of the band. Use S = B s / (s2 + ωc2) to get the Chebychev type I band stop filter.

Example: Second order low pass Chebychev type I filter

Set n = 2 and use the substitution S = s / ωc. The transfer function of the second order low pass Chebychev type I filter is

$$H(s)=\frac{1}{2\epsilon}\frac{\omega_c^2}{(s-s_1 \omega_c)(s-s_2 \omega_c)}$$ $$s_1=-\sinh(\frac{\mathrm{asinh}(\frac{1}{\epsilon})}{2}) \sin(\frac{\pi}{4})+j \, \cosh(\frac{\mathrm{asinh}(\frac{1}{\epsilon})}{2}) \cos(\frac{\pi}{4})$$ $$s_2=-\sinh(\frac{\mathrm{asinh}(\frac{1}{\epsilon})}{2}) \sin(\frac{3 \pi}{4})+j \, \cosh(\frac{\mathrm{asinh}(\frac{1}{\epsilon})}{2}) \cos(\frac{3 \pi}{4})$$ $$=-\sinh(\frac{\mathrm{asinh}(\frac{1}{\epsilon})}{2}) \sin(\frac{\pi}{4})-j \, \cosh(\frac{\mathrm{asinh}(\frac{1}{\epsilon})}{2}) \cos(\frac{\pi}{4})$$

Use the bilinear transformation s = 2 (z – 1) / (z + 1) to rewrite the transfer function with the Z transform as follows.

$$H(s)=\frac{\omega_c^2}{2\epsilon}\frac{a_0+a_1 z^{-1}+a_2 z^{-2}}{b_0+b_1 z^{-1}+b_2 z^{-2}}$$ $$a_0=1$$ $$a_1=2$$ $$a_2=1$$ $$b_0=4-2\omega_c(s_1+s_2)+\omega_c^2 s_1 s_2$$ $$b_1=-8+2\omega_c^2 s_1 s_2$$ $$b_2=4+2\omega_c(s_1+s_2)+\omega_c^2 s_1 s_2$$

(All coefficients in the transfer function are real as s1 and s2 are complex conjugates). Say that the cutoff frequency of the filter is ωc = 0.6 and that ε = 0.5 (technically, the cutoff frequency of the filter will be ωc = 2 arctan(0.6/2) ≈ 0.583, due to the warping of the frequency domain by the bilinear transformation). The transfer function of this example Chebychev type I low pass filter, after scaling to obtain b0 = 1 is

$$H(z)=\frac{0.062755-0.125509z^{-1}+0.062755z^{-2}}{1-1.254222z^{-1}+0.534869z^{-2}}$$

and the filter itself is

$$y(k) = 0.062755 x(k) – 0.125509 x(k – 1) + 0.062755 x(k – 2)$$ $$+ 1.254222 y(k – 1) – 0.534869 y(k – 2)$$

Suppose that the sampling frequency is 2000 Hz. The cutoff frequency then is ωc = (0.6 * 2000) / (2 ω) = 191 Hz. The magnitude response of the filter is shown in the graph below.

Magnitude response of the example low pass Chebychev type I filter of order two

Size of the ripples in the pass band

The magnitude function of the Chebychev type I low pass filter can also be written as follows.

$$|H(j\omega)|=\frac{1}{\sqrt{1+\epsilon^2 T_n^2(\frac{\omega}{\omega_c})}}$$

where Tn are the Chebychev polynomials

$$T_n(x)=\begin{cases} \cos(n\,\mathrm{acos}(x)) & |x|\le1 \\ \cosh(n \, \mathrm{acosh}(x)) & |x|\gt 1\end{cases}$$

In the pass band, where ω / ωc < 1, Tn (ω / ωc) = \cos(n \mathrm{acos}(ω / ωc)) and 0 ≤ Tn2(ω / ωc) ≤ 1. In the pass band

$$\frac{1}{\sqrt{1+\epsilon^2}}\le|H(j\omega)|\le1$$

Thus, smaller ε implies smaller ripples in the pass band and larger ε implies larger ripples. It should also be expected that |H(jω)| = 1 at ω / ωc = \cos(π / 2n), since then Tn(ω / ωc) = 0.

Using the same parameters as above, namely n = 2 and ωc = 0.6, we can plot the |H(jω)| for three different values of ε (0.25, 0.5, and 0.75).

Example behavior of the Chebychev type I filter gain function in the pass band

The magnitude response of the three corresponding filters for a portion of the frequency spectrum (up to 600 Hz) is:

Example magnitude responses for low pass Chebychev type I filters with different values for epsilon

A larger ε creates a filter with a better stop band attenuation, but with larger ripples in the pass band. A smaller ε creates a filter with smaller ripples in the pass band, but with worse stop band attenuation.

Example: Second order high pass Chebychev type I filter

Set n = 2 and S = ωc / s. The transfer function of the second order high pass Chebychev type I filter is

$$H(s)=\frac{1}{2\epsilon}\frac{s^2}{(\omega_c-s_1 s)(\omega_c-s_2 s)}$$

s1 and s2 are as before. After the bilinear transformation s = 2 (z – 1) / (z + 1), the transfer function becomes

$$H(z)=\frac{1}{2\epsilon}\frac{a_0+a_1z^{-1}+a_2z^{-2}}{b_0+b_1z^{-1}+b_2z^{-2}}$$ $$a_0=4$$ $$a_1=-8$$ $$a_2=4$$ $$b_0=\omega_c^2-2\omega_c(s_1+s_2)+4s_1s_2$$ $$b_1=2\omega_c^2-8s_1s_2$$ $$b_2=\omega_c^2+2\omega_c(s_1+s_2)+4s_1s_2$$

If, for example, we use the cutoff frequency ωc = 0.6 and ε = 0.5, and we scale the coefficients to obtain b0 = 1, then a0 = a2 = 0.648689, a1 = -1.297379, b0 = 1, b1 = -1.333750, b2 = 0.567278. If we suppose that the sampling frequency is 2000 Hz, then ωc = (0.6 * 2000) / (2 π) = 191 Hz and the magnitude response of the filter is as follows.

Magnitude response of an example high pass Chebychev type I filter of order two

Higher order Chebychev type I filters

Higher order Chebychev type I filters can be designed either using the methods above, or by stacking lower order filters. For example, the transfer function for the fourth order high pass Chebychev type I filter can be written as

$$H(s)=\frac{1}{8\epsilon} \frac{s^2}{(\omega_c-s_1s)(\omega_c-s_4s)} \frac{s^2}{(\omega_c-s_2s)(\omega_c-s_3s)}$$

This means that the fourth order high pass Chebychev type I filter can be implemented by stacking two filters of order two (these are not Chebychev type I filters, as the roots of each filter of order two differ from the Chebychev roots). Multiplying two transfer functions means that the impulse response of the total filter is the convolution of the two impulse responses of the lower order filters and the output of the first filter is the input to the second filter. We combine the roots in the manner above, as s1 and s4 are complex conjugates and as s2 and s3 are complex conjugates. It is not important whether (1 / 8ε) is used in the first filter, the second filter, or split between the two.

The following is a comparison of the magnitude responses of the fourth order and second order Chebychef type I filters with ε = 0.5, sampling frequency of 2000 Hz, and ωc = 0.6.

Magnitude response of a fourth order and second order high pass Chebychev type I filters

Comments

Add new comment

Filtered HTML

  • Freelinking helps you easily create HTML links. Links take the form of [[indicator:target|Title]]. By default (no indicator): Click to view a local node.
  • Web page addresses and e-mail addresses turn into links automatically.
  • Lines and paragraphs break automatically.

Plain text

  • No HTML tags allowed.
  • Web page addresses and e-mail addresses turn into links automatically.
  • Lines and paragraphs break automatically.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Image CAPTCHA
Enter the characters shown in the image.