Lecture9-IIR Filter
Lecture9-IIR Filter
Response Filters
Presenteed By
Dr M.Murugappan
School of Mechatronic Engineering
Universiti Malaysia Perlis
Introduction
A digital filter is a linear time invariant (LTI) discrete time system.
The FIR and IIR filters are of type of non-recursive and recursive
type, respectively.
In FIR filter design, the present output sample depends on the present
and previous input samples.
In IIR filter design, the present output sample depends on the present,
past and output samples.
The Impulse response for realizable filter and The stability condition
must satisfy.
Disadvantage:
Quantization error arises due to finite length of the representation of signals and
parameters.
Analog Lowpass Filter Design
General form analog filter transfer function is:
N M must satisfied and H(s) must lie in left half of the s-plane.
Analog lowpass Butterworth filter
Magnitude function of Butterworth lowpass filter is given by
Round N = 7
s
Lowpass with cutoff 0 c
( 2 1 )s
Bandstop with cutoffs 1 and 2
s 2 1 2
Design of IIR filters from analog filters
The conversion technique should be effective it should posses following
desirable properties.
The jΩ –axis in the s-plane should map into the unit circle in the z-
plane. Thus, have direct relationship between two frequency variable in two
domain.
The left–half plane of the s-plane should map into the inside of the unit
circle in the z-plane. Thus, we can convert stable analog to stable digital
filter.
4 most widely use Methods for digitizing Analog filter to digital filter
Approximation of derivatives.
Impulse invariant transformation.
Bilinear transformation.
matched z-transformation technique.
Design of IIR Filter using Impulse Invariance Technique
IIR filter is design such that unit impulse response h(n) of digital filter is the
sampled version of the impulse response of analog filter. The z-transform of
infinite impulse response given by
Let us consider the mapping points from the s-plane to the z-plane by the
relation z=esT. Substitute s=σ+jΩ and express the complex variable z in polar
form: z=rejω
rejω = e(σ+jΩ)T , we r = eσT, ω = ΩT.
Therefore, analog is mapped to a place in the z plane of magnitude e σT and
angle ΩT
Real part of analog pole =radius z-
plane,
Imaginary part=angle of digital pole,
Consider any pole on jΩ -axis, where
σ=0. Poles maps at the z-plane at a
radius r=e0.T=1. Therefore, the impulse
invariance had map poles from the s-
plane’s jΩ -axis to z-plane’s unit
circle.
2nd case
Consider pole on left–half s-plane
where σ < 0.Therefore, all s-plane
poles with negative real parts map to
z-plane poles inside the unit circle –
stable analog poles are mapped to
stable digital poles. Because r= e σT<1
for <0.
Unstable pole mapping occur when all poles at right half of the s-plane map
to the digital poles outside the unit circle.
Third case
many point in s-plane are mapped in one point in z-plane .
Easiest way to explain is to consider two poles in the s=plane with identical
real parts.
S 1 = , S 2=
Let Ha(s) is the system function of an analog filter and {ck} are the coefficients and
{pk} are the poles of analog filter.
For high sampling rates (small T), the digital gain is high, we can use
Step to design a digital filter using impulse
invariance method
For given specifications, find Ha(s), transfer function of analog filter.
Select sampling rate of the digital filter, T second per sample.
Express analog transfer function as sum of single-pole filters.
Which an be written
From differential eq
Which implies
The system function of the digital filter is
Solution:
In bilinear transformation
Ass T= 1 sec.
Then,
Realization of Digital Filters
There are two type of realization of digital filter transfer
function.
Structure call
Direct form 1
Realize the second order digital filter
y(n) =
Direct form II realization
Consider the difference equation
Which gives
The realization Eq.(5.114) and Eq.(5.115) shown in Fig.(5.35) ,(5.36)
Realize the second order system y(n)
2 1
H(jΩ ) N 1, 2, .......
Ω
1 ε 2 C N2
ΩP
- - - - - (1)
1
Where ε is a parameter of the filter related to C N (x) cos(Ncos x), | x | 1 (Passband)
the ripple in the passband
CN(x) is the Nth order Chebyshev 1
C N (x) cosh(Ncosh x) ,| x | 1 (Stopband)
polynominal
α P 10 log (1 ε 2 ) C N (1) 1 Pole locations for Chebyshev Filter
0.1α p
ε (10 1)0.5 μ ε 1 1 ε 2
2 Ωs
α s 10 log 1 ε C N
2
The poles of a Chebyshev filter
Ω P μ 1/N μ 1/N
a ΩP
2
1 10 0.1 α s
1
cosh 0.1 α p μ 1/N μ 1/N
10 1 b ΩP
N 2
-1 Ω s π (2k 1) π
cosh φk k 1, 2, ..., N
ΩP 2 2N
s k a cos φ k jbsin φ k
Comparison between Butterworth and Chebyshev Filter
The magnitude response of Butterworth filter decreases monotonically as the frequency Ω
increases from 0 to ∞, whereas the magnitude response of the Chebyshev filter exhibits
ripples in the passband or stopband according to the type.
The transition band is more in Butterworth filter when compared to Chebyshev filter.
The poles of the Butterworth filter lie on a circle, whereas the poles of the Chebyshev filter
lie on the ellipse.
s k a cos φ k jbsin φ k
5. Find the denominator polynomial of the transfer function using above
poles.
6. The numerator of the transfer function depends on the value of N.
(a) For N odd substitute s = 0 in the denominator polynomial and find
the value. This value is equal to the numerator of the transfer
function.
(b) For N even substitute s = 0 in the numerator polynomial and divide
the result by √1+ε2. This value is equal to the numerator.
Determine the order and the poles of a type I lowpass Chebyshev filter that
has a 1 dB ripple in the passband and passband frequency Ωp = 1000π, a
stopband frequency of 2000π and an attenuation of 40dB or more.
10 0.1 α s 1
1
cosh 0.1 α
10 p 1
N 4.536
-1 Ω s π (2k 1) π
cosh φk k 1, 2, ..., 5
N= 5 ΩP 2 2N
φ 1 108 ; φ 2 144 ; φ 3 180 ;
ε (10 0.1 α P
1) 0.5
0.508 φ 4 216
; φ 5 252
μ ε 1 1 ε 2 4.17
s 1 a cos φ 1 jbsin φ 1 89.5 π j989 π
s 2 a cos φ 2 jbsin φ 2 234.2 π j612 π
The poles of a Chebyshev filter
s 3 a cos φ 3 jbsin φ 3 289.5 π
μ 1/N μ 1/N
a Ω P 289.5 π s 4 a cos φ 4 jbsin φ 4 234.2 π j612 π
2
s 5 a cos φ 5 jbsin φ 5 89.5 π j989 π
μ 1/N
μ 1/N
b Ω P 1041 π
2
Given the specifications αp = 3dB ; αs = 16 dB ; fp = 1kHz, and fs = 2kHz, Determine the order of the
filter using Chebyshev approximation. Find H(s).
Step 1: cosh 1
0.1 α
10 p 1
10 0.1 α s 1
Find N N
-1 Ω s
1.91
cosh
ΩP
Step 3: The values of minor axis and major axis can be found as below
=(1414.38)2π2
Design a Chebyshev low pass filter with the specifications α p = 1 dB ripple in the
passband 0 ≤ ω ≤ 0.2π, αs = 15 dB ripple in the stopband 0.3π ≤ ω ≤ π, using (a), bilinear
transformation, (b). Impulse invariance.
2 ω 0.3 π
Ωs tan s 2tan 1.02
T 2 2
π (2k 1) π
10 1
0.1 α p
φk k 1, 2,3,4
cosh 1 2 2N
10 0.1 α s 1 φ1 112.5 ; φ 2 157.5 ; φ 3 202.5 ;
N 3.01
Ω
cosh -1 s φ 4 247.5 ;
ΩP
Let us take N 4
s 1 a cos φ 1 jbsin φ 1 0.0907 j0.639
s 2 a cos φ 2 jbsin φ 2 0.2189 j0.2647
ε (10 0.1 α P
1) 0.5
0.508
s 3 a cos φ 3 jbsin φ 3 0.2189 j 0.2647
1 2
μ ε 1ε 4.17 s 4 a cos φ 4 jbsin φ 4 0.0907 j0.639
The poles of a Chebyshev filter The denominator of H(s) =[(s+0.0907)2 +(0.639)2] [(s+0.2189)2
μ 1/N μ 1/N +(0.2647)2]
a ΩP 0.237 =(s2+0.1814s+0.4165) (s2+0.4378s+0.118)
2
As N is even, the numerator of H(s) =(0.4165) (0.118)/√1+ε2
μ 1/N μ 1/N
b ΩP 0.6918
2 =0.04381
0.001836(1 z 1 ) 4
(1 1.499z 1 0.8482z 2 ) ( 1 1.5548z 1
0.6493z 2 )
0.1 α
Impulse Invariance Method: cosh
10 p 1
1
μ ε 1 1 ε 2 4.17 π (2k 1) π
φk k 1, 2,3,4
2 2N
φ 1 112.5 ; φ 2 157.5 ; φ 3 202.5 ;
The poles of a Chebyshev filter
φ 4 247.5 ;
μ 1/N
μ 1/N
a ΩP 0.229
2
s 1 a cos φ 1 jbsin φ 1 0.0876 j0.619
μ 1/N μ 1/N
b ΩP 0.67 s 2 a cos φ 2 jbsin φ 2 0.2115 j0.2564
2
s 3 a cos φ 3 jbsin φ 3 0.2115 j0.2564
s 4 a cos φ 4 jbsin φ 4 0.0876 j0.619
The denominator of H(s) =[(s+0.0876)2 +(0.619)2] [(s+0.2115)2 +(0.2564)2]
=(s2+0.175s+0.391) (s2+0.423s+0.11)
As N is even, the numerator of H(s) =(0.391) (0.11)/√1+ε2
=0.03834