SP Slides 5
SP Slides 5
Assuming the input, xa (t) is bandlimited and the sampling frequency is high
enough to avoid aliasing, then overall system behaves like an LTI system:
(
H(ej⌦T ), |⌦| < ⇡/T
He↵ (j⌦) =
0, |⌦| ⇡/T
Passband: Ideal gain: 20 log10 (1) = 0 dB, Max gain: 20 log10 (1.01) = 0.0864 dB
Min gain: 20 log10 (0.99) = 0.0873 dB; Max stopband: 20 log10 (0.001) = 60 dB
Common Analog Filters
Continous-Time Filter Design
The magnitude-squared frequency response for lowpass prototypes are all of the
form:
1
|H(j⌦)|2 =
1 + g(j⌦)
where g is a polynomial or ratio of polynomials.
I Lowpass prototypes can be converted to a di↵erent lowpass filter, or a
highpass/bandpass/bandstop filter by performing a frequency transformation
(either in the analog or digital domain)
I Commonly used analog filters: Butterworth, Chebyshev, Elliptic
Butterworth Lowpass Filters
1
|H(j⌦)|2 =
1 + (j⌦/j⌦c )2N
The roots of the denominator polynomial are the poles of the magnitude-square
response:
sk = ( 1)1/2N (j⌦c ) = ⌦c d(j⇡/2N )(2k+N 1)
⇥ ⇤
I N th -order Chebyshev polynomial: VN (x) = cos N cos 1 (x)
I VN (x) oscillates between ±1 for |x| 1 and increases monotonically for |x| > 1
I |H(j⌦)|2 ripples between 1 and 1/(1 + ✏2 ) for 0 ⌦/⌦c 1
I |H(j⌦)|2 decreases monotonically for ⌦/⌦c > 1
I ✏ specified by allowable passband ripple, ⌦c is desired passband cuto↵
frequency, N is chosen so these specifications are met
Chebyshev Filters (Continued)
I Poles lie on an ellipse in the s-plane
I minor axis: 2a⌦c
I a = 12 (↵1/N ↵ 1/N )
p
I ↵=✏ 1+ 1+✏ 2
I major axis = 2b⌦c
I b = 12 (↵1/N + ↵ 1/N
)
I Identify points on major and minor circles
with symmetry (imaginary axis) and equal
angle spacing of ⇡/N
I Pole: horizontal distance from minor circle
and vertical distance from major circle
I Angle spacing corresponds to poles for
equivalent N th order Butterworth filter
Elliptic Filters
Various methods for converting an analog filter into the digital domain should
possess these desirable properties:
1. The j⌦ axis in the s-plane should map into the unit circle in the z-plane.
There will be a direct relationship between the two frequency variables in the
two domains.
2. The left-half plane (LHP) of the s-plane should map into the inside of the unit
circle in the z-plane. A stable analog filter will convert to a stable digital filter.
We specify the desired filter characteristics for the magnitude and then accept the
phase response that is obtained from the design method.
IIR Filter Design by Approximation of Derivatives
Substitute the backward di↵erence to find the derivative at time t = nT :
Since analog di↵erentiator has system function H(s) = s and the digital system
that produces the output [y[n] y[n 1]] /T has the system function
H(z) = (1 z 1 )/T . Therefore, the frequency domain equivalent relationship is:
1 z 1
s=
T
The kth derivative of y(t) results in this equivalent frequency-domain relationship:
✓ 1
◆k
1 z
sk =
T
Disadvantages of Approximation of Derivatives
Given an analog filter with system function, Ha (s), the system function for a
digital IIR filter using approximation of derivatives by finite di↵erences:
1 z 1 1
s= =) z=
T 1 sT
Using s = j⌦:
1 1 ⌦T
z= = +j
1 j⌦T 1 + ⌦2 T 2 1 + ⌦2 T 2
As ⌦ varies from 1 to 1 the corresponding z-plane is a circle of radius 12 with
center z = 12 . So this is stable, but confined to relatively small frequencies.
IIR Filter Design: Impulse Invariance
Impulse Invariance Design
h[n] = Td hc (nTd )
where Td is the sampling interval. Note that Td is not related to the sampling
period T if the discrete-time system has C/D or D/C conversion.
Relationship between between frequency responses of the discrete-time filter and
the continuous-time filter:
1
X ✓ ◆
j! ! 2⇡
H(e ) = Hc j +j k
Td Td
k= 1
Aliasing in Impulse Invariance Design
If the continuous-time filter is bandlimited:
Hc (j⌦) = 0 |⌦| ⇡/Td
then responses are related by scaling of frequency axis, ! = ⌦Td for |!| ⇡
✓ ◆
j! !
H(e ) = Hc j |!| ⇡
Td
Any practical filter can’t be bandlimited so aliasing will occur:
P
1 ⇣ ⌘
H(ej! ) = Hc j T!d + j 2⇡
Td k
k= 1
Impulse Invariance Design Procedure
Assuming that aliasing is negligible, the we use the frequency scaling to obtain the
continuous-time filter specifications from the desired discrete-time filter:
⌦ = !/Td
XN
Ak
Hc (s) =
s sk
k=1
8
< P A e sk t ,
> N
k t 0
hc (t) = k=1
>
:0, t<0
Impulse Invariance Design Procedure (Continued)
Using the z-transform provides the system function for the discrete-time system:
1 N
! N 1 N
X X X X X T d Ak
sk T d n n sk T d 1 n
H(z) = Td Ak (e ) z = T d Ak (e z ) =
1 e sk T d z 1
n=0 k=1 k=1 n=0 k=1
With the trapezoid rule, the integrator equation is, in discrete time:
T
y[n] = y[n 1] + (x[n] + x[n 1])
2
Going to the z-domain:
1 T 1
Y (z) = z Y (z) +
X(z) + z X(z)
✓2 ◆
Y (z) T 1+z 1
H(z) = =
X(z) 2 1 z 1
Bilinear Transformation (Continued)
✓ 1
◆
1 T 1+z
The continuous time integration is replaced by 1
for discrete time.
s 2 1 z
This mapping from the s-plane to the z-plane is the bilinear transformation (BLT):
✓ 1
◆
2 1 z
s= 1
T 1+z
Given a continuous time system function, Hc (s), the corresponding discrete time
system function is: ✓ ◆
21 z 1
H(z) = Hc
T 1+z 1
Bilinear Transformation Properties
Properties of the Bilinear Transformation: Stability
21 z 1 1 + sT /2
s= 1
! z=
T 1+z 1 sT /2
Using s = + j⌦:
Prewarping:
2 ⇣! ⌘
p
⌦p = tan
T 2
2 ⇣! ⌘
s
⌦s = tan
T 2
Design Procedure Using Bilinear Transformation
2. If you design a continuous filter with ⌦p and ⌦s and apply the BLT, the
resulting critical frequencies will be warped to values lower than desired. So
apply the inverse frequency warping to the desired !p and !s , resulting in
2 ⇣! ⌘
p,s
prewarped continuous filter specs ⌦0p,s = tan .
T 2
3. Design continuous filter Hc (s) with prewarped critical frequencies ⌦0p and ⌦0s .
4. Apply the BLT to get the discrete time filter H(z) = Hc (s) s= T2 1 z 1 .
1+z 1
The critical frequencies will be warped back to the desired values !p and !s .
Design Comparisons:
Discrete-Time Butterworth, Chebyshev, and Elliptic Filters
Example: Lowpass Filter Specifications To Compare
20 log10 ( s ) = 30 =) s = 0.0316
Butterworth Filter: 15th-order
into a bandpass filter with upper and lower band edge frequencies ⌦u and ⌦l .
1
H(s) = s2 +⌦l ⌦u
s 2 + ⌦ l ⌦u s(⌦u ⌦l ) +1
s ! ⌦p
s(⌦u ⌦l ) (⌦u ⌦l )s
=
s2 + (⌦u ⌦l )s + ⌦l ⌦u
This filter has a zero at s = 0 and poles at
q
(⌦u ⌦l ) ± ⌦2u + ⌦2l 6⌦u ⌦l
s=
2