05.Filter design and approximation
.pdf
94 FILTER DESIGN AND APPROXIMATION
5.5THE ELLIPTIC FILTER
The low-pass filter can be characterized as having a pass band from ω D 0 to ω D ωp with an attenuation no greater than H0 plus a small ripple. In addition it is characterized as having a stop band from ωs to 1 with an insertion loss no less than some high value, ˛min (Fig. 5.5). In the Chebyshev filter the pass-band ripple is fixed to a certain maximum, but small, value while the attenuation in the stop band increases monotonically with ω. The inverse Chebyshev filter produces an equal ripple in the stop band and a monotonically decreasing insertion loss for ω going from ωp toward ω D 0. The elliptic function filter equal ripple response in both the pass band and in the stop band. This design provides a way of not throwing away excess stop-band attenuation at high frequencies by allowing redistribution of the attenuation over the whole stop band. As a consequence the rate of cutoff may be increased by putting some of the transmission zeros near the pass band. The cost for having equal ripple response in both the pass band and in the stop band is a slightly more complicated circuit topology for the elliptic filter (Fig. 5.5).
There is no simple recursion formula for the design of elliptic function filters. Typically tables of values are derived numerically [4,5] and are used for the lowpass prototype filter. These tabulated values have been incorporated in a program called ELLIPTIC. In this program the desired maximum attenuation level in the pass band, minimum attenuation in the stop band, the frequencies where the pass band ends and the stop band begins, and finally the number of poles, n, are balanced against each other to provide an elliptic filter design. If so desired, the program will produce a SPICE net list that can be used to analyze the design. In the PSPICE version of SPICE, the voltage is plotted using V(21) or VDB(21) to display the insertion loss on a linear or log scale, respectively.
5.6MATCHING BETWEEN UNEQUAL RESISTANCES
For a low-pass filter, perfect match cannot in principle be achieved when
impedance |
matching is |
used. In |
the preceding Butterworth and Chebyshev |
|||||||||||||||||||
|
|
|
|
|
|
|
L2 |
|
L4 |
L6 |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C2 |
|
|
C4 |
|
|
C6 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
50 |
|
|
C1 |
|
C3 |
|
|
C5 |
|
|
C7 |
|
|
|
|
|
50 |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
FIGURE 5.5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
A |
|
seven-pole low-pass |
elliptic filter |
topology. |
When |
fp D 0.8 GHz |
||||||||||||||||
and fs D 1 GHz, C1 D 3.285 pF, C2 D 0.547 pF, |
L2 D 12.653 nH, |
C3 D 5.459 pF, |
||||||||||||||||||||
C4 |
D 2.682 pF, |
L4 D 9.947 nH, C5 D 4.846 pF, C6 D 2.040 pF, |
L6 D 8.963 nH, and |
|||||||||||||||||||
C7 |
D 2.231 pF. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
MATCHING BETWEEN UNEQUAL RESISTANCES |
95 |
functions, the constant, H0, is 1, since a passive filter cannot produce gain greater than 1. When the input and output resistance levels are equal, then H0 is 1. The ratio of the load to generator resistances introduces a constraint on H0. For Butterworth filters this constraint is
RL |
D |
1 |
p |
|
|
|
š1 |
||
1 |
H0 |
||||||||
|
1 C p |
|
|
5.41 |
|||||
RG |
|||||||||
1 |
H0 |
||||||||
|
|
|
|
|
|
|
|
|
|
while for Chebyshev filters this constraint is given as follows [1]:
RL |
|
1 |
|
|
p |
|
|
|
|
š1 |
|
|
|
|
|
|||
|
1 |
|
H0 |
|
|
|
|
|
||||||||||
|
D |
1 C p |
|
|
|
|
, |
n odd |
5.42 |
|||||||||
RG |
||||||||||||||||||
1 |
H0 |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
D |
p |
|
|
|
|
|
|
|
|
š1 |
|
||||||
RL |
1 ε2 |
C |
|
|
1 |
ε2 |
H0 |
|
||||||||||
|
|
C |
1 C |
|
H0 |
, n even |
5.43 |
|||||||||||
RG |
p1 ε2 |
|
ε2 |
|||||||||||||||
|
|
|
|
|
C |
|
|
|
C |
|
|
|
|
|
||||
One might wonder if the generator and load consisted of complex impedances, what technique might be used for matching. Without getting too involved with that issue, it is known that such matching is not always possible. The impedances must be “compatible” for matching to occur. One thing a designer can do, though, is try to incorporate the reactive part of the load into the filter as much as possible.
5.6.1The Darlington Procedure
A doubly terminated filter can be designed for any physically realizable transfer function. A variety of different circuit realizations may be possible, but only one will be described. However, this particular realization method is widely used and provides practical filter design. Approximation theory determines the transfer function jH jω j2 that comes closest to the ideal filter characteristics. In a lossless, low-pass circuit with possibly unequal termination resistances (Fig. 5.6), the reflected power j jω j2 can be found:
1 jH jω j2 D j jω j2 D jω jω |
5.44 |
L1 |
L3 |
|
RG |
C2 |
RL |
FIGURE 5.6 Butterworth low-pass filter with unequal resistance terminations. When RG D 20 and RL D 80 , L1 D 127.75 H, C2 D 0.01804 F, and L3 D 43.400 H.
96 FILTER DESIGN AND APPROXIMATION
The final expression results from the magnitude being the product of the reflection coefficient and its complex conjugate. This can be generalized by replacing jω with the complex frequency s:
s s D 1 jH sj2 |
5.45 |
The right-hand side is a known function that is given in the form of a ratio of polynomials in s. A requirement for realizability of an impedance or reflection coefficient is that it be positive real. All the poles of the function must lie in the left half side of the complex plane in order to avoid unrealizable growing exponentials. Half of the poles of j sj2 lie in the left half side and half in the right half side of the complex frequency plane. The function s can be extracted from s s by choosing only those poles in the left half side. The zeros of the function need some further consideration. The choice of which zeros to choose is more arbitrary, since there is not the same realizability restrictions on the zeros. If the choice is made to use only the left half-plane zeros, the resulting reflection coefficient and the corresponding driving point impedance is the minimum phase function. The jω-axis zeros are even multiples of complex conjugate pairs and are divided equally between s and s.
The problem of actually finding the poles and zeros requires finding the roots of the denominator and numerator polynomials. While these roots can be found analytically for the Butterworth and Chebyshev filters, the roots for other functions such as the Thompson-Bessel filter function must be found numerically. The transfer function takes the form
jH sj2 D H0
F s
where F0 D 1.
The dc transfer function is
j 0 j2 D 1 H0
0 D š 1 H0
5.46
5.47
5.48
Since at dc in a low-pass circuit the series reactive elements are short circuits and the shunt reactive elements are open circuits, the reflection coefficient is
0 |
|
RL RG |
|
5.49 |
|
D RL C RG |
|||||
|
|
||||
Consequently |
4RL/RG |
|
|||
H0 D |
5.50 |
||||
RL/RG C 1 2 |
|||||
which of course equals 1 when both sides of the filter have equal terminations.
MATCHING BETWEEN UNEQUAL RESISTANCES |
97 |
Once the reflection coefficient is determined, the Darlington synthesis procedure is used to obtain the circuit elements. The input impedance to the filter at any frequency is given in terms of the reflection coefficient:
Z |
R |
|
1 C s |
5.51 |
|
G 1 s |
|||||
in D |
|
|
|||
The Cauer extraction technique for a ladder network can now be used. The polynomials in the numerator and denominator are arranged in descending powers of s. It will always be the case for a lossless transfer function that the highest power of the numerator and denominator will differ by at least 1. If the numerator is the higher-order polynomial, then an impedance pole at s D 1 (i.e., a series inductor) can be extracted from the impedance function. This is done by synthetic division. The fractional remainder is now inverted, and synthetic division again carried out to extract an admittance pole at s D 1 (i.e., a shunt capacitor). The process continues until only the load resistance or conductance remains.
As an example, consider a three-pole Butterworth filter with a 3 dB cutoff frequency at 1 rad/s. The input resistance RG D 20 and the output resistance is RL D 80 . The Butterworth transfer function is therefore
|
|
|
jH ω j D |
|
|
H0 |
|
|
5.52 |
|||||
|
|
|
1 C ω6 |
|
||||||||||
where from Eq. (5.50), |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
H0 D |
16 |
|
|
|
5.53 |
|||||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
25 |
|
|
|
|||||||
and |
|
|
|
|
|
|
|
|
|
|
1 C ω6 H0 |
|
||
j |
ω 2 |
D |
1 |
j |
H ω 2 |
D |
5.54 |
|||||||
|
||||||||||||||
j |
|
|
j |
|
1 C ω6 |
|
||||||||
Now replace ω with js, factor the denominator into the six roots of 1, and recombine into two cubic factors where one factor contains the left half-plane roots and the other the right half-plane roots. This is the standard Butterworth polynomial:
s |
|
s |
|
9/25 s6 |
|
5.55 |
|
D s3 C 2s2 2s C 1 s3 C 2s2 |
|
||||
|
|
C 2s C 1 |
||||
In this case the denominator is readily factored analytically, but the roots of the numerator when H0 D6 1 must be found numerically. The program POLY can provide the complex roots of a polynomial with complex coefficients. In this example all values are calculated using double precision arithmetic, though for clarity only three or four significant figures are shown.
98 FILTER DESIGN AND APPROXIMATION
The reflection coefficient containing only left half-plane poles and zeros is
s |
D |
s3 C 1.687s2 C 1.423s C 0.599 |
|
5.56 |
|
s3 C 2s2 C 2s C 1 |
|||||
|
|
||||
The input impedance is found from Eq. (5.51):
Z |
20 |
2s3 C 3.687s2 C 3.423s C 1.599 |
5.57 |
|
0s3 C 0.313s2 C 0.577s C 0.400s |
||||
in D |
|
|
Extraction of the impedance pole at s D 1 is done by synthetic division:
6.387s
0.313s2 C 0.577s C 0.400 2s3 C 3.687s2 C 3.423s C 1.599
2s3 C 3.687s2 C 2.555s
0s3 C 0s2 C 0.868s C 1.599
The remainder is inverted, and an admittance pole at s D 1 is extracted:
0.361s
0.868s C 1.599 0.313s2 C 0.577s C 0.400 0.313s2 C 0.577s
0s2 C 0s C 0.868s C 0.400
By inversion again and performing synthetic division, once more another impedance pole at s D 1 is removed:
2.170s
0.400 0.868s C 1.599
The final remainder, 1.599/0.400 D 4.000 represents the normalized load resistance, which is the expected value. Hence L10 D 6.387 H, C02 D 0.361 F, and L30 D 2.170 H. The impedance level of the circuit is now adjusted from 1 to RG D 20 by multiplying all the inductances and dividing all capacitances by 20 . Thus L10 becomes L1 D 127.75 H, C02 becomes C2 D 0.01804 F, and L30 becomes L3 D 43.400 H. The final circuit is shown in Fig. 5.6. Verification of this circuit is shown by a SPICE analysis found in Fig. 5.7. Near zero frequency the insertion loss is 0.8 or 1.938 dB and at 1 rad/s (0.159 Hz); the loss has increased by 3 dB.
Easier analytical methods are available for the Chebyshev filter, and these are in fact used in the Chebyshev impedance transforming circuit described in
MATCHING BETWEEN UNEQUAL RESISTANCES |
99 |
Transmission
1.00
0.90
0.80
0.70
0.60
0.50
0.40
0.30
0.20
0.10
0.00
0.0 |
0.1 |
0.2 |
0.3 |
0.4 |
0.5 |
|
|
Frequency, Hz |
|
|
|
FIGURE 5.7 Frequency response of the Butterworth low-pass filter.
Section 5.6.3. The Darlington method shown here can be used where a closed form solution for the roots is not available.
5.6.2Filter Type Transformation
Filter design is based on the design of a low-pass prototype circuit whose impedance level is 1 and whose low-pass cutoff frequency is ωc D 1 rad/s. If the desired impedance level is to be changed from 1 to RL, then all inductors and resistors should be multiplied by RL and all capacitors should be divided by RL. If the circuit elements of the low-pass prototype are denoted by a “p” subscript, then the new adjusted values can be found:
L D RLLp |
5.58 |
||
C D |
Cp |
5.59 |
|
RL |
|
||
R D RLRp |
5.60 |
||
To adjust the cutoff frequency from 1 rad/s to ωc, the low-pass circuit elements are further modified in the following way:
L0 |
L |
5.61 |
D ωc |
100 FILTER DESIGN AND APPROXIMATION
C0 |
D |
C |
5.62 |
ωc |
|||
R0 |
D R |
5.63 |
|
Transformation of the low-pass filter to a high-pass filter can be accomplished by another frequency transformation. The normalized complex frequency variable for the low-pass prototype circuit is sn. On the jω-axis the pass band of the low-pass filter occurs between ω D 1 and C1. If the cutoff frequency for the high-pass filter is ωc, then the high-pass frequency variable is
ωc |
5.64 |
s D sn |
Applying this transformation will transform the pass-band frequencies of the lowpass filter to the pass band of the high-pass filter. This is illustrated in Fig. 5.8. The reactance of an inductor, L, in the low-pass filter becomes a capacitance, Ch, in the high-pass filter:
Lsn D |
Lωc |
1 |
5.65 |
|||
|
|
D |
|
|||
s |
Chs |
|||||
or |
|
|
1 |
|
|
|
Ch D |
|
|
|
5.66 |
||
|
|
|||||
Lωc |
||||||
Similarly application of the frequency transformation Eq. (5.64) will convert a capacitor in the low-pass filter to an inductor in the high-pass filter:
1 |
5.67 |
||||||||
|
|
|
Lh D |
|
|||||
Cωc |
|||||||||
jω |
|
|
|
+jω |
|
+jω o |
|||
|
|
|
|||||||
+1j |
|
|
|
|
|
|
|||
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
σ |
|
|
|
σ |
|||
–1j |
|
|
|
|
–j ω |
|
–jω o |
||
|
|
|
|||||||
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
||||
Low pass |
|
High pass |
|||||||
FIGURE 5.8 Low-pass to high-pass transformation.
MATCHING BETWEEN UNEQUAL RESISTANCES |
101 |
A band-pass filter is specified to have a pass band from ω1 to ω2. The “center” of the pass band is the geometric mean of the band edge frequencies, ω0 D pω1ω2. The fractional bandwidth is w D ω2 ω1 /ω0. A band-pass circuit can be formed from the low-pass prototype by using a frequency transformation that will map the pass band of the low-pass filter to the pass band of the band-pass filter. The desired frequency transformation is
sn D w |
ω0 |
C |
s0 |
5.68 |
|
1 |
|
s |
|
ω |
|
where s is the frequency variable for the bandpass circuit. To verify this expression for the jω-axis, Eq. (5.58) is rewritten as
ωn D w |
ω0 |
ω0 |
5.69 |
|
1 |
|
ω |
ω |
|
A short table of specific values for the normalized low-pass prototype circuit and the corresponding band-pass frequencies are shown in Table 5.1.
A graphic illustration of the frequency transformation is shown in Fig. 5.9. A consequence of this transformation is that an inductor L in the low-pass prototype filter becomes a series LC circuit in the band-pass circuit:
Lsn D |
Ls |
C |
Lω0 |
5.70 |
wω0 |
ws |
Similarly a capacitance in the low-pass filter is transformed to a parallel LC circuit:
Csn D |
Cs |
C |
Cω0 |
5.71 |
wω0 |
ws |
Finally, the low-pass to band-stop filter frequency transformation is the reciprocal of Eq. (5.68):
sn D w |
s |
|
ω |
|
1 |
|
C |
0 |
5.72 |
||
ω0 |
s |
All these transformations from the low-pass prototype filter are summarized in Fig. 5.10.
TABLE 5.1 Low-Pass to Bandpass Mapping
Bandpass ω |
Low-pass ωn |
ω2 |
C1 |
ω0 |
0 |
ω1 |
1 |
ω1 |
C1 |
ω0 |
0 |
ω2 |
1 |
102 FILTER DESIGN AND APPROXIMATION
jω |
+jω |
ω2 |
ω1
+1j
σ |
σ |
–1j 
–ω 1
–jω
–ω 2
Low pass |
Band pass |
FIGURE 5.9 Low-pass to band-pass transformation.
Low pass |
High pass |
Band pass |
Band stop |
||
|
|
|
|
Lw |
|
|
|
|
|
ωo |
|
L |
1 |
L |
w |
|
|
|
|
|
|||
|
ω c L |
ωoW |
Lωo |
1 |
|
|
|
|
|
|
|
|
|
|
|
Lωow |
|
|
|
|
C |
|
|
C |
1 |
ωoW |
1 |
Cw |
|
|
|
||||
ωcC |
|
|
|||
|
|
w |
Cωow |
ωo |
|
|
|
|
|||
|
|
|
|
|
|
|
|
ωoC |
|
|
|
FIGURE 5.10 Filter conversion chart.
5.6.3Chebyshev Bandpass Filter Example
The analytical design technique for a Chebyshev filter with two unequal resistances has been implemented in the program called CHEBY. As an example of its use, we will consider the design of a Chebyshev filter that matches a 15 to a 50 load resistance. It will have n D 3 poles, center frequency of 1.9 GHz, a fractional bandwidth w D f2 f1 /f0 D 20%. The program CHEBY is used
MATCHING BETWEEN UNEQUAL RESISTANCES |
103 |
to find the filter circuit elements. The program could have used the Darlington procedure described in Section 5.6.1, but instead it used the simpler analytical formulas [1]. The following is a sample run of the program:
Generator AND Load resistances |
15.,50. |
|
|||||
Passband ripple |
(dB) |
0.2 |
|
|
|
|
|
Bandpass Filter? |
hY/Ni |
Y |
|
|
hA/Ni |
|
|
Specify stopband attenuation OR n, |
N |
||||||
Number of transmission poles n = |
3 |
|
|
||||
L(1) = |
.62405E + 02 |
C(2) = |
.25125E — 01 |
L(3) = |
.36000E + 02 |
||
Number of poles = 3 Ripple = |
.20000E + 00 dB |
|
|||||
Center Frequency, Fo (Hz), AND Fractional Bandwidth, |
|||||||
w |
1.9E9,.2 |
|
|
|
|
|
|
Through series LC. L1( 1) = |
.261370E — 07 C1( 1) |
||||||
=.268458E — 12
L 1 |
C 1 |
|
L 3 |
C 3 |
R G |
C 2 |
L 2 |
|
R L |
FIGURE 5.11 A 15 : 50 Ohm Chebyshev band-pass filter, where L1 D 26.14 nH, C1 D 0.2685 pF, L2 D 0.6668 nH, C1 D 10.52 pF, L3 D 15.08 nH, and C3 D 0.4654 pF.
Amplitude, dB
0.00
–5.00
–10.00
–15.00
–20.00
–25.00
–30.00 |
1.5 |
1.6 |
1.7 |
1.8 |
1.9 |
2.0 |
2.1 |
2.2 |
2.3 |
2.4 |
1.4 |
Frequency, GHz
FIGURE 5.12 SPICE analysis of a Chebyshev filter.
