<< Chapter < Page Chapter >> Page >

This is the procedure used in the design Program 8 in the appendix.

When designing a bandpass elliptic-function filter, four frequencies must be specified: the lower stopband edge, the lowerpassband edge, the upper passband edge, and the upper stopband edge. All four must be prewarped to the equivalent analog values.A problem occurs when the two transition bands of the bandpass filter are converted into the single transition band of thelowpass prototype filter. In general they will be inconsistant; therefore, the narrower of the two transition bands should beused to specify the lowpass filter. The same problem occurs in designing a bandreject elliptic-function filter. Program 8 inthe appendix should be studied to understand how this is carried out.

An alternative to the process of converting a lowpass analog into a bandpass analog filter which is then converted into a digitalfilter, is to first convert the prototype lowpass analog filter into a lowpass digital filter and then make the conversion into abandpass filter. If the prototype digital filter transfer function is H p ( z ) and the frequency transformation is f ( z ) , the desired transformed digital filter is described by

H ( z ) = H p ( f ( z ) )

Since the frequency response of both H ( z ) and H p ( z ) is obtained by evaluating them on the unit circle in the- z plane, f ( z ) should map the unit circle onto the unit circle ( | z | = 1 = > | f ( z ) | = 1 ) . Both H ( z ) and H p ( z ) should be stable; therefore, f ( z ) should map the interior of the unit circle into the interior of the unit circle ( | z | < 1 = > | f ( z ) | < 1 ) . If f ( z ) were viewed as a filter, it would be an “all-pass" filter with a unity magnitude frequency response of the form

f ( z ) = p ( z ) z n p ( 1 / z ) = a n z n + a n - 1 z n - 1 + + a 0 a 0 z n + a 1 z n - 1 + + a n

The prototype digital lowpass filter is usually designed with bandedges at ± π / 2 . Determining the frequency transformation then becomes the problem of solving the n + 1 equations

f ( e j ω i ) = e ± j π / 2 = ( - 1 ) i j

for the unknown a k where i = 0 , 1 , 2 , n and the ω i are the bandedges of the desired transformed frequency response put in ascending order. The resulting simultaneous equationshave a special structure that allow a recursive solution. Details of this approach can be found in [link] .

This is an extremely general approach that allows multiple passbands of arbitrary width. If elliptic-function approximationsare used, only one of the transition bandwidths can be independently specified. If more than one passband or rejectband is desired, f(z)will be higher order than second order and, therefore, the transformed transfer function H ( f ( z ) ) will have to be factored using a root finder.

To illustrate the results of using transform methods to design filters, three examples are given which are designed withProgram 8 from the appendix.

Design of an chebyshev highpass filter

The specifications are given for a highpass Chebyshev frequency response with a passband edge at f p = 0 . 3 Hertz with a sampling rate of one sample per second. The order is set at N = 5 and the passband ripple at 0.91515 dB. The transfer function is

H ( z ) = ( z - 1 ) ( z 2 - 2 z + 1 ) ( z 2 - 2 z + 1 ) ( z + 0 . 64334 ) ( z 2 + 0 . 97495 z + 0 . 55567 ) ( z 2 + 0 . 57327 z + 0 . 83827 )

The frequency response plot is given in [link] .

Figure three is titled chebyshev highpass filter. Its horizontal axis is labeled Normalized frequency, and ranges in value from 0 to 1 in increments of 0.2. Its vertical axis is labeled Magnitude Response and ranges in value from 0 to 1 in increments of 0.2. There are two blue curves on the graph. The first is a simple horizontal line extending from (0, 0) to (1, 0). The second stretches out from this horizontal axis at approximately (0.3, 0) and begins increasing at an increasing rate all the way up to nearly the top of the page at approximately (0.6, 1). At this point the curve peaks and continues downward sharply, then back up to another peak more shallowly that occurs at approximately (0.75, 1), then back down and up in an even more shallow fashion to its termination point at (1, 1).
Fifth Order Digital Chebyshev Highpass Filter

Design of an elliptic-function bandpass filter

This filter requires a bandpass frequency response with an elliptic-function approximation. The maximum passpand ripple is one dB,the minimum stopband attenuation is 30 dB, the lower stopband edge f 1 = 0 . 19 , the lower passband edge f 2 = 0 . 2 , the upper passband edge f 3 = 0 . 3 , and the upper stopband edge f 4 = 0 . 31 Hertz with a sampling rate of one sample per second. The design program calculated a requiredprototype order of N = 6 and, therefore, a total order of 10. The frequency response plot is shown in [link] .

Figure four is titled tenth-order elliptic function filter. The horizontal axis is labeled normalized frequency and ranges in value from 0 to 1 in increments of 0.2. The vertical axis is labeled magnitude response and ranges in value from 0 to 1 in increments of 0.2. There are two blue curves on this graph. The first is a simple horizontal line from (0, 0) across the figure to (1, 0). The second is a more complex, but symmetric, curve that begins from the left with a very shallow positive slope. It decreases back down to the horizontal axis and wavers at (0.4, 0), where it then increases sharply to (0.4, 1). At this point, the graph peaks and follows a sinusoidal chape of an amplitude of approximately 0.05 with an increasing wavelength. There are two troughs and a wide peak in the middle of the graph. At this peak, the graph moves to the right in a completely symmetric fashion, mirroring the aforementioned description.
Tenth Order Digital Elliptic Bandpass Filter

Design of an inverse-chebyshev bandreject filter

The specifications are given for a bandreject Inverse- Chebyshev frequency response with bandedges at f s = 0 . 1 and 0 . 2 Hertz with a sampling rate of one sample per second. The order is set at N = 11 and the minimum stopband attenuation at 30 dB. The frequency response plot is given in [link] .

Figure five is a graph titled Inverse chebyshev band reject filter. The horizontal axis is labeled Normalized frequency, and ranges in value from 0 to 1 in increments of 0.2. The vertical axis is labeled Magnitude response and ranges in value from 0 to 1 in increments of 1. There are two blue curves on the graph. The first is a simple horizontal line from (0, 0) to (1, 0). The second begins at (0, 1) and extends to the right horizontally until just before (0.2, 1) where it sharply moves downward to (0.2, 0). At this point, the curve contains 10 peaks of a sinusoidal wave-shape of an extremely small amplitude. The wavelength of this section of the curve starts out small, then by the fifth  peak is its longest, then by the tenth peak is back to the wavelength of the first trough. This sinusoidal segment occurs horizontally from 0.2 to 0.4. After (0.4, 0), the curve returns sharply to the top of the figure at (0.4, 1), where it then continues horizontally to the right edge of the figure and terminates at (1, 1).
Twenty second Order Digital Inverse Chebyshev Band Reject Filter

Summary

This section has described the two most popular and useful methods for transforming a prototype analog filter into a digitalfilter. The analog frequency variable is used because a literature on analog filter design exists, but more importantly,many approximation theories are more straightforward in terms of the Laplace-transform variable than the z-transform variable. Theimpulse-invariant method is particularly valuable when time- domain characteristics are important. The bilinear-transformmethod is the most common when frequency-domain performance is the main interest. Use of the BLT warps the frequency scale and,therefore, the digital band edges must be prewarped to obtain the necessary band edges for the analog filter design. Formulas thattransform the analog prototype filters into the desired digital filters and for prewarping specifications were derived.

The use of frequency transformations to convert lowpass filters into highpass, bandpass, and bandreject filters wasdiscussed as a particularly useful combination with the bilinear transformation. These are implemented in Program 8 and designexamples from this program were shown.

There are cases where no analytic results are possible or where the desired frequency response is not piecewise constant andtransformation methods are not appropriate. Direct methods for these cases are developed in other sections.

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Digital signal processing and digital filter design (draft). OpenStax CNX. Nov 17, 2012 Download for free at http://cnx.org/content/col10598/1.6
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Digital signal processing and digital filter design (draft)' conversation and receive update notifications?

Ask