Finite Impulse Response (FIR)
Digital Filters (I)
Types of linear phase FIR filters
Yogananda Isukapalli
1
Key characteristic features of FIR filters
1. The basic FIR filter is characterized by the following two
equations:
N -1 N -1
y ( n) = å h( k ) x ( n - k ) (1) H ( z) = å
k =0
h( k ) z - k (2)
k =0
where h(k), k=0,1,…,N-1, are the impulse response
coefficients of the filter, H(z) is the transfer function
and N the length of the filter.
2. FIR filters can have an exactly linear phase response.
3. FIR filters are simple to implement with all DSP processors
available having architectures that are suited to FIR filtering.
2
Linear phase response
• Consider a signal that consists of several frequency
components passing through a filter.
1. The phase delay (T p) of the filter is the amount of
time delay each frequency component of the signal
suffers in going through the filter.
2. The group delay (T g ) is the average time delay the
composite signal suffers at each frequency.
3. Mathematically, T p = -q (w ) / w (3)
T g = -dq (w ) / dw (4)
where q(w) is the phase angle.
3
• A filter is said to have a linear phase response if,
J (w ) = -aw (5)
J (w ) = b - aw (6)
where a and b are constants.
Example) Ideal Lowpass filter
k
Magnitude
0 wc p w
Phase
0 a wc p w
4
k e-jwa passband
H(ejw) = 0 otherwise
Magnitude response = |H(ejw)| = k
Phase response (q(w)) = < H(ejw) = -wa
Follows: y[n] = kx[n-a] : Linear phase implies that the
output is a replica of x[n] {LPF} with a time shift of a
-p -wu -wl 0 wl wu p w
-p -wu -wl 0 wl wu p w
5
Linear phase FIR filters
• Symmetric impulse response will yield linear phase FIR filters.
1. Positive symmetry of impulse response:
n = 0,1,…, (N-1)/2 (N odd)
h(n) = h(N-n-1),
n = 0,1,…, (N/2)-1 (N even)
a = (N-1)/2 in eqn (5)
1. Negative symmetry of impulse response:
n = 0,1,…, (N-1)/2 (N odd)
h(n) = -h(N-n-1),
n = 0,1,…, (N/2)-1 (N even)
a = (N-1)/2
b = p/2 in eqn (6)
6
Example) For positive symmetry and N = 11 odd length
h[n] = h[N-1-n] n = 0,1,….(N-1)/2
h[0] = h[10]
h[1] = h[9]
h[2] = h[8]
h[3] = h[7]
h[4] = h[6]
h[5] = h[5]
7
N = 11
0 1 2 3 4 5 6 7 8 9 10
N = 10
0 1 2 3 4 5 6 7 8 9
h[n] = h[10-1-n] = h[9-n] h[0] = h[9]
h[1] = h[8]
h[2] = h[7]
h[3] = h[6]
h[4] = h[5] 8
Example : Proof of Linear Phase :
4
H ( z) = å h[n]z
n =0
-n
H ( z ) = h[0] + h[1] z -1 + h[ 2] z -2 + h[3] z -3 + h[ 4] z -4
For N = 5; Symmetric impulse Response Implies :
h[n] = h[5-1-n] = h[4-n]
h[0] = h[4]
h[1] = h[3]
h[2] = h[2]
Now : H ( z ) = h[ 2] z -2 + [ h[1] z + h[3] z -1 ] z -2 +
+ [ h[0] z 2 + h[ 4] z -2 ] z -2
9
Consider Frequency Response :
H (e jw ) = H ( z ) z = e jw
( T = 1)
H (e jw ) = h[ 2]e - j 2 w + [ h[1]e jw + h[3]e - jw ]e - j 2 w +
[h[0]e j 2 w + h[4]e - j 2 w ]e - j 2 w
= e - j 2 w [[h[2] + 2h[1] cos w + 2h[0] cos 2 w]
é 1
ù
=e - j 2w
ê
ë
[ h[ 2] + å
n =0
2 h[ n ] cos( w( n - 2)) ú
û
= e - j 2 w | H (e jqw ) | ( Linear Phase form)
Phase = -2w
10
Group Delay :
-d
Tg = ( phase)
dw
H =2
passband
0 wp p w
Tg
2w
0 wp p w
Group delay is constant over the passband for linear phase
filters.
11
Types of FIR linear phase systems
1. Type I FIR linear phase system
The impulse response is positive symmetric and N an
odd integer
h[n] = h[ N - 1 - n], 0 £ n £ ( N - 1) / 2
The frequency response is
( N -1) / 2
H (e jw
)= ån =0
h[n]e - jwn
( N -1) / 2
H (e jw
)=e - jw( N -1) / 2
å a[n] cos(wn)
n =0
where
a[0] = h[( N - 1) / 2],
a[n] = 2h[(( N - 1) / 2) - n],
n = 1,2,...( N - 1) / 2.
12
2. Type II FIR linear phase system
The impulse response is positive symmetric and
N is an even integer
h[n] = h[ N - 1 - n], 0 £ n £ ( N / 2) - 1
The frequency response is
( N / 2 ) -1
H (e jw
)= ån =0
h[n]e - jwn
ìN / 2
- jw( N -1) / 2 ï 1 üï
H (e jw
)=e å
í b[n] cos[ w(n - )]ý,
ïî n =1 2 ïþ
where
b[n] = 2h[ N / 2 - n],
n = 1,2,...N / 2.
13
3. Type III FIR linear phase system
The impulse response is negative-symmetric and
N an odd integer.
h[n] = -h[ N - 1 - n], 0 £ n £ ( N - 1) / 2
The frequency response is
( N -1) / 2
H (e jw
)= å
n =0
h[n]e - jwn
ì( N -1) / 2
- jw( N -1) / 2 ï
üï
H (e jw
) = je í å
ïî n =1
a[n] sin( wn)ý,
ïþ
where
a[n] = 2h[(( N - 1) / 2) - n],
n = 1,2,...( N - 1) / 2.
14
4. Type IV FIR linear phase system
The impulse response is negative-symmetric and N an
even integer.
h[n] = -h[ N - 1 - n], 0 £ n £ ( N / 2) - 1
The frequency response is
( N / 2 ) -1
H (e jw
)= ån =0
h[n]e - jwn
ìN / 2
- jw( N -1) / 2 ï 1 üï
H (e jw
) = je í å
ïî n =1
b[ n ] sin[ w( n -
2
)]ý,
ïþ
where
b[n] = 2h[ N / 2 - n],
n = 1,2,..., N / 2.
15
Fig: A summary of four types of linear phase FIR filters1
16
Fig: A comparison of the impulse of the four types of linear phase FIR filters1
17
Notes:
• The frequency response of a Type 2 filter is always zero at f=0.5
( half the sampling frequency as all the frequencies are normalized
to the sampling frequency) and thus is unsuitable as a highpass
filter.
• The frequency response of a Type 3 filter is always zero at f=0
and 0.5, while that of Type 4 filter is zero at f=0. Thus, Type 3
filter cannot be used as either a lowpass or high pass filter
whereas Type 4 cannot be used as a lowpass filter.
• Both Type 3 and 4 filters introduce a 90o phase shift and are
often used to design differentiators and Hilbert transformers.
• Type 1 is the most versatile of the four.
18
• The phase delay (for type 1 and 2 filters) or group delay (for
all four types) is expressible in terms of the number of
coefficients of the filter.
• Thus they can be corrected to give a zero phase or group delay
response.
• For example,
æ N -1ö Phase delay for
Tp = ç ÷T
è 2 ø types 1 and 2
æ N -1- p ö Group delay for
Tg = ç ÷T
è 2 ø types 3 and 4
where T is the sampling period.
19
References
1. “Digital Signal Processing – A Practical Approach” -
Emmanuel C. Ifeachor and Barrie W. Jervis
Second Edition
20