0% found this document useful (0 votes)
46 views11 pagesFIR Filter
Uploaded by
Ehtisham BhattiCopyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF or read online on Scribd
0% found this document useful (0 votes)
46 views11 pagesFIR Filter
Uploaded by
Ehtisham BhattiCopyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF or read online on Scribd
/ 11
“EEE305", “EEESO] Part A": Digital Signal Processing Chapter 4: Design of FIR Filters
Chapter 4
Design of FIR Filters
4.1 Introduction
Digital FIR filters have many favourable properties, which is why they are extremely popular in digital signal
processing. One of these properties is that they may exhibit linear phase, which means that signals in the passband will
suffer no dispersion, Dispersion ocours when different frequency components of a signal have a different delay through
asystem.
The simplest design method for FIR filters is impulse response mumeation (IRT), but unfortunately it has undesirable
frequency-domain characteristics, owing to the Giéb’s phenomenon. The second design method for a FIR filter that we
shall cover in this Chapter is the windowing technique. The windowing method can be used to mitigate the adverse
effects of impulse response truncation,
4.2 Fourier Transform Relationship
‘The frequency response of a generalised FIR filter is defined by Equation 4.1 below
ay
‘To find the frequency response. can be replaced bye, as follows:
ef seh ae a2
Equation 4.1 can be reswriten as Equation 4.3
HQ)= Yh) 43)
The inverse Fourier transform of H (2) gives an expression for the impulse response of the FIR filter:
Hml=L fr1()-2*% 0 aa)
ms
Equations (4.3) and (4.4) form a Fourier transform pair. There is a similarity between this pair and the pair for a
periodic signal x(t) , as shown in equations (4.5):
x= SxW 2
(ita)
fue ee at as)
a
Both of the functions x(t) and H(Q) are periodic. The Fourier transform of x(¢) gives the discrete function of fy
— in other words the discrete harmonies of the spectrum. The inverse Fourier transform of H(Q) gives the diserete
funtion of in other words the discrete samples ofthe impulse response of the filter:
42.1 Low-pass Filter design
Consider the ideal low-pass filter frequency response. as illustrated in Figure 4.1 below, with a normalised angular eut—
off frequency ©. Usually the subscript D is used to distinguish between the ideal and actual, impulse and frequency
responses of a filter.
University of Neweastle upon Tyne Page 4.1 “EEE305", “EEESO] Part A": Digital Signal Processing Chapter 4: Design of FIR Filters
HO)
De a a, 2 a
Figure 4.1: Ideal low-pass filter frequency response.
The impulse response of an ideal low-pass filter ho[n] is found by substituting Hp(@) = 1 in Equation 4.5 and
integrating between the limits of the cut-off frequencies [-Q., QJ,
hot = fred {|| =x
al ae | |, 72a
oe 1 2jsin@aQ.)
jn
jn
‘Multiplying both the numerator and denominator by Q.. the Equation above becomes:
holr) =o we
‘The impulse response of an ideal low-pass filler can also be re-written by replacing Q. ~ 2RF. in the Equation above, to
obtain it in terms of the normalised cut-off frequency F.. This is illustrated in Equation 4.6 below,
nQ,
Agfa) =2F, Sine) ao
2.)
In Chapter 2, we found that the Fourier transform of a rectangular window is a sine funetion, which is same for the
Jmpolse response ofa lowpass filter, a6 illustrated in Figure 4.2 below
fe 4.2: Impulse response of an ideal low-pass filter.
4.3. FIR Filter Design by Impulse Response Truncation (IRT)
With reference to Figure 4.2, although fila] decays to either side of n = 0 it theoretically continues for ever in both
directions. This reflects a general antithesis between band limitation and time limitation; since we have chosen a
frequency response with a sharp cut-off (or brick wall response), then the time-domain response continues forever. To
realise such a filter the impulse response is sruncated in some way or other. One approach is to ignore the small sample
values atthe ends and shift [a] to begin at 1 = 0, giving a causal flter as depicted in Figure 4.3.
In general, the more samples we include of hifi] the closer we get to the desired form of Hp(Q), but the less economic
the filter becomes due to the relative number of computations. In practice we must settle for an approximation to the
ideal frequency response. It is usually customary to truncate the impulse response to N= (2M + 1) terms. In Figure 4.3,
the impulse response of the ideal low-pass filter is truncated to M= 9 samples and is delayed by M samples
University of Neweastle upon Tyne Page 4. “EEE305", “EEESO] Part A": Digital Signal Processing Chapter 4: Design of FIR Filters
Figure 4.3: Truncated impulse response to 2M=1 samples, with a delay of M samples.
‘The 2-transfer function of the filter now becomes:
Yano an
The ideal impulse responses for a low-pass, high-pass, band-pass and band-stop filters are depicted in Table 4.1 below.
Filter type hotel. a#0 hl,
sin@oQ,)
ow-pass ap, nee.) 2
Lowe a F,
sin(nd,)
igh-pass 1-2p, Sneee.) 2
High a 1
sin(vQ,) sin(nQ,)
088 op, WO) op i 2F,-2.
Band-ps = 712, Fy OF,
Band-stoy ror, Sl) 9p sine) | _ be oR]
° 0 22. poet
‘Table 4.1: deal impulse responses for various FIR filter types.
4.4 Summary of FIR Filter Design Using The IRT Method
© Choose the ideal frequency response Hy(Q), depending on the type of filter (e.g. low-pass, high-pass ete), rom
Table 4.1
* Caleulate the impulse response of the ideal filter hofn), using the inverse Fourier transform formula of Equation
45.
Finally, truncate the ideal impulse response to N= (2M + 1) terms,
4.5 Optimality of the IRT Method
‘A measure of the error between the frequency response of the ideal filter H((2) and that of the designed filter H(Q) is
the integral of the squared error, defined by Equation 4.8.
E(Q)
+ ftom -Heafan as)
ae}
Filters designed by the IRT method are considered optimal in the sense of minimising the integral of the squared error.
However, the IRT approach is not considered practical because of the oscillatory nature of the frequency response near
the discontinuity points of the desired response Hy(Q). It is the case that as NV’ increases, the magnitude of the
‘oscillations remains the same, but they become more localised to the discontinuity, with the result that E(Q) decreases
‘with increasing NV. This property is known as the Gibb's phenomenon and is described in detail in Section 4.6
University of Neweastle upon Tyne Page 43 “EEE305", “EEESO] Part A": Digital Signal Processing Chapter 4: Design of FIR Filters
4.6 Gibb’s Phenomenon
‘Truncating the impulse response introduces undesirable ripples and overshoots in the Srequency response. This effect is
known as the Gibb’s phenomenon and is illustrated in Figure 4.4. As an example, the impulse response for a low-pass
filter is truncated with M = 9, 25 and an infinite number of samples. Notice how the undesirable ripples and avershoots
are clearly visible in the frequency response.
pop
ELA
‘o. a” a, a
Figure 4.4: Effects on the frequency response of truncating the ideal impulse response.
In the time-domain truncation is achieved by multiplying the impulse response with a window fimetion w(w)
Conversely, in the frequency-domain truncation is achieved by convolving the frequency response H(Q) with W(Q),
where W(Q) is the Fourier transform of the window function. In this example, the truncation was achieved by
multiplying the impulse response with a rectangular window function, The height of the ripples is nor dependent on the
value of M, it only affects the width. The heights of the ripples are however controlled by the type of window function
cused for the time-domain multiplication. Table 4.2 illustrates the mathematical definitions, and Figure 4.5 shows the
‘characteristic shape and amplitude responses of various window functions.
Name of window finetion wD Mathematical efinition
Rectangular 1
2a
Hanning 05-0. seo{ 2]
=
Hamas 0.54-0s6co =
N-1
2
Blackman 042-0. seo aa I 0.08c0:
Where, I,(x)=
-1,(B) “
Kaiser
Table 4.2: Mathematical definitions for various window functions.
University of Neweastle upon Tyne Page 44 “EEE305", “EEESO] Part A”: Digital Signal Processing
Figure 4,5: Window functions and the
Chapter 4: Design of FIR Filters
t
respective amplitude responses in dB,
University of Neweastle upon Tyne
Page 45 “EEE305", “EEESO] Part A": Digital Signal Processing Chapter 4: Design of FIR Filters
4.7 Filter Specification Requirements
Figure 4.6 provides a graphical description of the specifications of a normalised low-pass filter. The shaded areas in the
pass band and in the stop band indicate the forbidden magnitude values in these bands, In the transition bend there are
no forbidden values but itis usually required that the magnitude decreases monotonically in this band,
oh,
ie
fated Tented Sips
Figure 4.6: A normalised low-pass filter specification.
‘The parameter 6, is the tolerance of the magnitude response in the stop-band and the desired magnitude response is
always lose to zero, The quantity 8, is known as the stap-Pand attenuation and is sometimes expressed in terms of GBs
using Equation 4.9
‘Minimum stop-band attenuation, A,
2010g(6) as)
‘The parameter dj is the tolerance of the magnitude response in the pass-band and the desired magnitude response is,
