Digital Signal Processing 10EC52

SUBJECT CODE : 10EC52 IA MARKS : 25

NO. OF LECTURE HRS/WEEK : 04 EXAM HOURS : 03
TOTAL NO. OF LECTURE HRS. : 52 EXAM MARKS : 100

PART - A
UNIT - 1
DISCRETE FOURIER TRANSFORMS (DFT): FREQUENCY DOMAIN SAMPLING AND
RECONSTRUCTION OF DISCRETE TIME SIGNALS. DFT AS A LINEAR TRANSFORMATION, ITS
RELATIONSHIP WITH OTHER TRANSFORMS.
7 HOURS
UNIT - 2
PROPERTIES DFT, MULTIPLICATION OF TWO DFTS- THE CIRCULAR CONVOLUTION,
OF
ADDITIONAL DFT PROPERTIES, USE OF DFT IN LINEAR FILTERING, OVERLAP-SAVE AND
OVERLAP-ADD METHOD.
6 HOURS
UNIT - 3
FAST-FOURIER-TRANSFORM (FFT) ALGORITHMS: DIRECT COMPUTATION OF DFT, NEED
FOR EFFICIENT COMPUTATION OF THE DFT (FFT ALGORITHMS).
8 HOURS
UNIT - 4
RADIX-2 FFT ALGORITHM FOR THE COMPUTATION OF DFT IDFT–DECIMATION-IN-
AND
TIME AND DECIMATION-IN-FREQUENCY ALGORITHMS. GOERTZEL ALGORITHM, AND CHIRP-Z
TRANSFORM
6 HOURS
PART - B
UNIT - 5
IIR FILTER DESIGN: CHARACTERISTICS OF COMMONLY USED ANALOG FILTERS –
BUTTERWORTH AND CHEBYSHEVE FILTERS, ANALOG TO ANALOG FREQUENCY
TRANSFORMATIONS.
6 HOURS
UNIT - 6
FIR FILTER DESIGN: INTRODUCTION TO FIR FILTERS, DESIGN OF FIR FILTERS USING -
RECTANGULAR, HAMMING, BARTLET AND KAISER WINDOWS, FIR FILTER DESIGN USING
FREQUENCY SAMPLING TECHNIQUE
6 HOURS

Page 1
Digital Signal Processing 10EC52

UNIT - 7
DESIGN OF IIR (BUTTERWORTH AND CHEBYSHEV) -
FILTERS FROM ANALOG FILTERS
IMPULSE INVARIANCE METHOD. MAPPING OF TRANSFER FUNCTIONS: APPROXIMATION OF
DERIVATIVE (BACKWARD DIFFERENCE AND BILINEAR TRANSFORMATION) METHOD,
MATCHED Z TRANSFORMS, VERIFICATION FOR STABILITY AND LINEARITY DURING MAPPING
7 HOURS
UNIT - 8
IMPLEMENTATION OF DISCRETE-TIME SYSTEMS: STRUCTURES FOR IIR AND FIR SYSTEMS-
DIRECT FORM I AND DIRECT FORM II SYSTEMS, CASCADE, LATTICE AND PARALLEL
REALIZATION.
6 HOURS

TEXT BOOK:
DIGITAL SIGNAL PROCESSING – PRINCIPLES ALGORITHMS & APPLICATIONS, PROAKIS &
MONALAKIS, PEARSON EDUCATION, 4TH EDITION, NEW DELHI, 2007.

REFERENCE BOOKS:
1. DISCRETE TIME SIGNAL PROCESSING, OPPENHEIM & SCHAFFER, PHI, 2003.
2. DIGITAL SIGNAL PROCESSING, S. K. MITRA, TATA MC-GRAW HILL, 2ND EDITION,
2004.
3. DIGITAL SIGNAL PROCESSING, LEE TAN: ELSIVIER PUBLICATIONS, 2007

Page 2
Digital Signal Processing 10EC52

INDEX SHEET
SL.NO TOPIC PAGE NO.
I Unit-1: Discrete Fourier Transforms 1-12
1.1 Frequency Domain Sampling
1.2 Reconstruction of Discrete time signals
1.3 DFT as a Linear Transform
1.4 DFT relationship with other transforms
II UNIT - 2 : Properties of DFT 13-25
2.1 Multiplication of two DFTs –Circular convolution
2.2 Additional DFT properties
2.3 Use of DFT in Linear Filtering
2.4 Overlap save and Overlap Add Method
2.5 Solution to Problems
III UNIT – 3 : Fast Fourier Transform Algorithms 26-31
3.1 Direct computation of DFT
3.2 Need for Efficient computation of DFT
3.2 FFT algorithms
UNIT – 4 : Radix-2 FFT Algorithms for DFT and
IV 32-52
IDFT
4.1 Decimation In Time Algorithms
4.2 Decimation-in-Frequency Algorithms
4.3 Goertzel Algorithms
4.4 Chirp-z-Transforms
V UNIT – 5 : IIR Filter Design 53-70
5.1 Characteristics of commonly used Analog Filters
5.2 Butterworth and Chebyshev Filters
5.3 Analog to Analog Frequency Transforms
5.4 Solution to problems
VI UNIT – 6 : FIR Filter Design 71-108
6.1 Introduction to FIR filters
Design of FIR Filters using Rectangular and Hamming
6.2
window
6.3 Design of FIR Filters using Bartlet and Hamming window
6.4 F IR filter design using frequency sampling technique
VII UNIT – 7 : Design of IIR Filters from Analog Filters 109-128
7.1 Impulse Invariance Method
7.2 Mapping of Transfer Functions
Approximation of derivative (Backward Difference and
7.3
Bilinear Transforms) method
7.4 Matched Z transforms
7.5 Verification for stability and Linearity during mapping
7.6 Solution to problems

Page 3
Digital Signal Processing 10EC52

VIII UNIT – 8 : Implementation of Discrete time systems 129-156

8.1 Structure of IIR and FIR systems
8.2 Direct form I and direct form II Systems
8.3 Cascade Realization
8.4 Lattice Realization
8.5 Parallel Realization
8.6 Solution to problems

Page 4
Digital Signal Processing 10EC52

UNIT 1

DISCRETE FOURIER TRANSFORMS (DFT)

CONTENTS:-

1. FREQUENCY DOMAIN SAMPLING

2. RECONSTRUCTION OF DISCRETE TIME SIGNALS

3. DFT AS A LINEAR TRANSFORMATION

4. DFT RELATIONSHIP WITH OTHER TRANSFORMS.

 DFT RELATIONSHIP WITH FOURIER SERIES

 DFT RELATIONSHIP WITH Z-TRANSFORMS

RECOMMENDED READINGS

1. DIGITAL SIGNAL PROCESSING – PRINCIPLES ALGORITHMS & APPLICATIONS, PROAKIS &

MONALAKIS, PEARSON EDUCATION, 4TH EDITION, NEW DELHI, 2007.

2. DISCRETE TIME SIGNAL PROCESSING, OPPENHEIM & SCHAFFER, PHI, 2003.

3. DIGITAL SIGNAL PROCESSING, S. K. MITRA, TATA MC-GRAW HILL, 2ND EDITION, 2004.

UNIT 1

Page 5
Digital Signal Processing 10EC52

Discrete Fourier Transform

1.1 Introduction:
Before we introduce the DFT we consider the sampling of the Fourier transform of an
aperiodic discrete-time sequence. Thus we establish the relation between the sampled Fourier
transform and the DFT.A discrete time system may be described by the convolution sum, the
Fourier representation and the z transform as seen in the previous chapter. If the signal is
periodic in the time domain DTFS representation can be used, in the frequency domain the
spectrum is discrete and periodic. If the signal is non-periodic or of finite duration the
frequency domain representation is periodic and continuous this is not convenient to
implement on the computer. Exploiting the periodicity property of DTFS representation the
finite duration sequence can also be represented in the frequency domain, which is referred to
as Discrete Fourier Transform DFT.
DFT is an important mathematical tool which can be used for the software
implementation of certain digital signal processing algorithms .DFT gives a method to
transform a given sequence to frequency domain and to represent the spectrum of the sequence
using only k frequency values, where k is an integer that takes N values, K=0, 1, 2,…..N-1.
The advantages of DFT are:
1. It is computationally convenient.
2. The DFT of a finite length sequence makes the frequency domain analysis much
simpler than continuous Fourier transform technique.

1.2 FREQUENCY DOMAIN SAMPLING AND RECONSTRUCTION OF DISCRETE

TIME SIGNALS:

Consider an aperiodic discrete time signal x (n) with Fourier transform, an aperiodic finite
energy signal has continuous spectra. For an aperiodic signal x[n] the spectrum is:

X w x ne jwn
………………………………(1.1)
n

Suppose we sample X[w] periodically in frequency at a sampling of w radians between

successive samples. We know that DTFT is periodic with 2 , therefore only samples in the

Page 6
Digital Signal Processing 10EC52

fundamental frequency range will be necessary. For convenience we take N equidistant

2
samples in the interval (0<=w<2 ). The spacing between samples will be w as shown
N
below in Fig.1.1.
X[w]

w
0 2

Fig 1.1 Frequency Domain Sampling

Let us first consider selection of N, or the number of samples in the frequency domain.
2 k
If we evaluate equation (1) at w
N
2 k
X xne j 2 kn / N
k 0,1,2,......., ( N 1) ………………………. (1.2)
N n

We can divide the summation in (1) into infinite number of summations where each sum
contains N terms.

1 N 1 2N 1
2 k j 2 kn / N j 2 kn / N j 2 kn / N
X ....... xne xne xne
N n N n 0 n N

lN N 1
j 2 kn / N
xne
l n lN

If we then change the index in the summation from n to n-l N and interchange the order of
summations we get:

N 1
2 k
x n lN e j 2 kn / N
for k 0,1,2,......, ( N 1) …….(1.3)
N n 0 l

Page 7
Digital Signal Processing 10EC52

Denote the quantity inside the bracket as xp[n]. This is the signal that is a repeating version of
x[n] every N samples. Since it is a periodic signal it can be represented by the Fourier series.

N 1
xp n ck e j 2 kn / N
n 0,1,2,........, ( N 1)
k 0

With FS coefficients:

N 1
1
ck xp n e j 2 kn / N
k 0,1,2,......., ( N 1) …………… (1.4)
N n 0

Comparing the expressions in equations (1.4) and (1.3) we conclude the following:

1 2
ck X k k 0,1,......., ( N 1) ………………. (1.5)
N N

Therefore it is possible to write the expression x p[n] as below:

N 1
1 2
xp n X k e j2 kn / N
n 0,1,....., ( N 1) ………. (1.6)
N k 0 N

The above formula shows the reconstruction of the periodic signal x p[n] from the samples of
the spectrum X[w]. But it does not say if X[w] or x[n] can be recovered from the samples.

Let us have a look at that:

Since xp[n] is the periodic extension of x[n] it is clear that x[n] can be recovered from x p[n] if
there is no aliasing in the time domain. That is if x[n] is time-limited to less than the period N
of xp[n].This is depicted in Fig. 1.2 below:

Page 8
Digital Signal Processing 10EC52

x[n]

0 L

xp[n] N>=L
No aliasing

0 L N

xp[n] N<L
Aliasing

0 N

Fig. 1.2 Signal Reconstruction

Hence we conclude:
The spectrum of an aperiodic discrete-time signal with finite duration L can be exactly
2 k
recovered from its samples at frequencies wk if N >= L.
N

We compute xp[n] for n=0, 1,....., N-1 using equation (1.6)

Then X[w] can be computed using equation (1.1).

1.3 Discrete Fourier Transform:

The DTFT representation for a finite duration sequence is

∞ -jωn
X (jω) = ∑ x (n) ℮
n= -∞
jωn
X (n) =1/2π ∫X (jω) e dω , Where ω═ 2πk/n

Page 9
Digital Signal Processing 10EC52

2π
Where x(n) is a finite duration sequence, X(jω) is periodic with period 2π.It is
convenient sample X(jω) with a sampling frequency equal an integer multiple of its period =m
that is taking N uniformly spaced samples between 0 and 2π.
Let ωk= 2πk/n, 0≤k≤N-1

∞ -j2πkn/N
Therefore X(jω) = ∑ x(n) ℮
n=−∞
Since X(jω) is sampled for one period and there are N samples X(jω) can be expressed
as
N-1 -j2πkn/N
X(k) = X(jω)│ ω=2πkn/N ═∑ x(n) ℮ 0≤k≤N-1
n=0

1.4 Matrix relation of DFT

The DFT expression can be expressed as

[X] = [x(n)] [WN]

T
Where [X] = [X(0), X(1),……..]

[x] is the transpose of the input sequence. WN is a N x N matrix

WN = 1 1 1 1 ………………1
1 wn1 wn2 wn3……………...wn n-1
1 wn2 wn4 wn6 ……………wn2(n-1)
…………………………………………….
…………………………………………….
1………………………………..wN (N-1)(N-1)

ex;
4 pt DFT of the sequence 0,1,2,3

X(0) 1 1 1 1
X(1) 1 -j -1 j
X(2) = 1 -1 1 -1
X(3) 1 j -1 -j

Solving the matrix X(K) = 6 , -2+2j, -2 , -2-2j

1.5 Relationship of Fourier Transforms with other transforms

Page 10
Digital Signal Processing 10EC52

1.5.1 Relationship of Fourier transform with continuous time signal:

Suppose that xa(t) is a continuous-time periodic signal with fundamental period Tp= 1/F0.The
signal can be expressed in Fourier series as

Where {ck} are the Fourier coefficients. If we sample xa(t) at a uniform rate Fs = N/Tp = 1/T,
we obtain discrete time sequence

Thus {ck‟} is the aliasing version of {ck}

1.5.2 Relationship of Fourier transform with z-transform

Let us consider a sequence x(n) having the z-transform

With ROC that includes unit circle. If X(z) is sampled at the N equally spaced points on the
unit circle Zk = e j2πk/N for K= 0,1,2,………..N-1 we obtain

The above expression is identical to Fourier transform X(ω) evaluated at N equally spaced
frequencies ωk = 2πk/N for K= 0,1,2,………..N-1.

Page 11
Digital Signal Processing 10EC52

If the sequence x(n) has a finite duration of length N or less. The sequence can be recovered
from its N-point DFT. Consequently X(z) can be expressed as a function of DFT as

Fourier transform of a continuous time signal can be obtained from DFT as

Page 12
Digital Signal Processing 10EC52

Recommended Questions with solutions

Question 1
The first five points of the 8-point DFT of a real valued sequence are {0.25, 0.125-j0.318, 0,
0.125-j0.0518, 0}. Determine the remaining three points

Ans: Since x(n) is real, the real part of the DFT is even, imaginary part odd. Thus the
remaining points are {0.125+j0.0518,0,0, 0.125+j0.318}.

Question 2
Compute the eight-point DFT circular convolution for the following sequences.
x2(n) = sin 3πn/8
Ans:

Question 3
Compute the eight-point DFT circular convolution for the following sequence
X3(n) = cos 3πn/8

Question 4

Page 13
Digital Signal Processing 10EC52

Define DFT. Establish a relation between the Fourier series coefficients of a continuous time
signal and DFT

Solution

The DTFT representation for a finite duration sequence is

∞
X (jω) = ∑ x (n) ℮-jωn
n= -∞

X (n) =1/2π ∫X (jω) e jωn dω , Where ω═ 2πk/n

2π
Where x(n) is a finite duration sequence, X(jω) is periodic with period 2π.It is
convenient sample X(jω) with a sampling frequency equal an integer multiple of its period =m
that is taking N uniformly spaced samples between 0 and 2π.
Let ωk= 2πk/n, 0≤k≤N
∞
Therefore X(jω) = ∑ x(n) ℮-j2πkn/N
n=−∞
Since X(jω) is sampled for one period and there are N samples X(jω) can be expressed
as
N-1
X(k) = X(jω)│ ω=2πkn/N ═∑ x(n) ℮-j2πkn/N 0≤k≤N-1
n=0
Question 5

Solution:-

Page 14
Digital Signal Processing 10EC52

Question 6

Find the 4-point DFT of sequence x(n) = 6+ sin(2πn/N), n= 0,1,………N-1

Solution :-

Question 7

Solution

Page 15
Digital Signal Processing 10EC52

Question 8

Solution

Page 16
Digital Signal Processing 10EC52

UNIT 2

PROPERTIES OF DISCRETE FOURIER TRANSFORMS (DFT)

CONTENTS:-

1. MULTIPLICATION OF TWO DFT‟S- THE CIRCULAR CONVOLUTION,

2. ADDITIONAL DFT PROPERTIES

3. USE OF DFT IN LINEAR FILTERING

4. OVERLAP-SAVE AND OVERLAP-ADD METHOD.

RECOMMENDED READINGS

1. DIGITAL SIGNAL PROCESSING – PRINCIPLES ALGORITHMS & APPLICATIONS, PROAKIS &

MONALAKIS, PEARSON EDUCATION, 4TH EDITION, NEW DELHI, 2007.

2. DISCRETE TIME SIGNAL PROCESSING, OPPENHEIM & SCHAFFER, PHI, 2003.

3. DIGITAL SIGNAL PROCESSING, S. K. MITRA, TATA MC-GRAW HILL, 2ND EDITION, 2004.

Page 17
Digital Signal Processing 10EC52

Unit 2
Properties of DFT

2.1 Properties:-

The DFT and IDFT for an N-point sequence x(n) are given as

In this section we discuss about the important properties of the DFT. These properties are
helpful in the application of the DFT to practical problems.

Periodicity:-

2.1.2 Linearity: If

Page 18
Digital Signal Processing 10EC52

Then A x1 (n) + b x2 (n) a X1(k) + b X2(k)

2.1.3 Circular shift:

In linear shift, when a sequence is shifted the sequence gets extended. In circular shift the
number of elements in a sequence remains the same. Given a sequence x (n) the shifted
version x (n-m) indicates a shift of m. With DFTs the sequences are defined for 0 to N-1.

If x (n) = x (0), x (1), x (2), x (3)

X (n-1) = x (3), x (0), x (1).x (2)

X (n-2) = x (2), x (3), x (0), x (1)

2.1.4 Time shift:

If x (n) X (k)
mk
Then x (n-m) WN X (k)

2.1.5 Frequency shift

If x(n) X(k)
+nok
Wn x(n) X(k+no)
N-1 kn
Consider x(k) = x(n) W n
n=0
N-1
(k+ no)n
X(k+no)= \ x(n) WN
n=0
kn non
= x(n) WN WN

non
X(k+no) x(n) WN

2.1.6 Symmetry:

Page 19
Digital Signal Processing 10EC52

For a real sequence, if x(n) X(k)

X(N-K) = X* (k)

For a complex sequence

DFT(x*(n)) = X*(N-K)

If x(n) then X(k)

Real and even real and even

Real and odd imaginary and odd
Odd and imaginary real odd
Even and imaginary imaginary and even

2.2 Convolution theorem;

Circular convolution in time domain corresponds to multiplication of the DFTs

If y(n) = x(n) h(n) then Y(k) = X(k) H(k)

Ex let x(n) = 1,2,2,1 and h(n) = 1,2,2,1

Then y (n) = x(n) h(n)

Y(n) = 9,10,9,8

N pt DFTs of 2 real sequences can be found using a single DFT

If g(n) & h(n) are two sequences then let x(n) = g(n) +j h(n)

G(k) = ½ (X(k) + X*(k))

H(k) = 1/2j (X(K) +X*(k))

2N pt DFT of a real sequence using a single N pt DFT

Let x(n) be a real sequence of length 2N with y(n) and g(n) denoting its N pt DFT

Let y(n) = x(2n) and g(2n+1)

k
X (k) = Y (k) + WN G (k)

Using DFT to find IDFT

The DFT expression can be used to find IDFT

Page 20
Digital Signal Processing 10EC52

X(n) = 1/N [DFT(X*(k)]*

2.3 Digital filtering using DFT

In a LTI system the system response is got by convoluting the input with the impulse
response. In the frequency domain their respective spectra are multiplied. These spectra are
continuous and hence cannot be used for computations. The product of 2 DFT s is equivalent
to the circular convolution of the corresponding time domain sequences. Circular convolution
cannot be used to determine the output of a linear filter to a given input sequence. In this case a
frequency domain methodology equivalent to linear convolution is required. Linear
convolution can be implemented using circular convolution by taking the length of the
convolution as N >= n1+n2-1 where n1 and n2 are the lengths of the 2 sequences.

2.3.1 Overlap and add

In order to convolve a short duration sequence with a long duration sequence x(n) ,x(n)
is split into blocks of length N x(n) and h(n) are zero padded to length L+M-1 . circular
convolution is performed to each block then the results are added. These data blocks may be
represented as

The IDFT yields data blocks of length N that are free of aliasing since the size of the
DFTs and IDFT is N = L+M -1 and the sequences are increased to N-points by appending
zeros to each block. Since each block is terminated with M-1 zeros, the last M-1 points from
each output block must be overlapped and added to the first M-1 points of the succeeding

Page 21
Digital Signal Processing 10EC52

block. Hence this method is called the overlap method. This overlapping and adding yields the
output sequences given below.

2.3.2 Overlap and save method

In this method x (n) is divided into blocks of length N with an overlap of k-1 samples.
The first block is zero padded with k-1 zeros at the beginning. H (n) is also zero padded to
length N. Circular convolution of each block is performed using the N length DFT .The output
signal is obtained after discarding the first k-1 samples the final result is obtained by adding
the intermediate results.

Page 22
Digital Signal Processing 10EC52

In this method the size of the I/P data blocks is N= L+M-1 and the size of the DFts and
IDFTs are of length N. Each data block consists of the last M-1 data points of the previous
data block followed by L new data points to form a data sequence of length N= L+M-1. An N-
point DFT is computed from each data block. The impulse response of the FIR filter is
increased in length by appending L-1 zeros and an N-point DFT of the sequence is computed
once and stored.
The multiplication of two N-point DFTs {H(k)} and {Xm(k)} for the mth block of data yields

Since the data record is of the length N, the first M-1 points of Ym(n) are corrupted by
aliasing and must be discarded. The last L points of Ym(n) are exactly the same as the result
from linear convolution and as a consequence we get

Page 23
Digital Signal Processing 10EC52

Page 24
Digital Signal Processing 10EC52

Recommended Questions with solutions

Question 1

State and Prove the Time shifting Property of DFT

Solution

The DFT and IDFT for an N-point sequence x(n) are given as

Time shift:

If x (n) X (k)
mk
Then x (n-m) WN X (k)

Question 2

State and Prove the: (i) Circular convolution property of DFT; (ii) DFT of Real and even
sequence.

Solution

(i) Convolution theorem

Circular convolution in time domain corresponds to multiplication of the DFTs

If y(n) = x(n) h(n) then Y(k) = X(k) H(k)
Ex let x(n) = 1,2,2,1 and h(n) = 1,2,2,1 Then y (n) = x(n) h(n)

Y(n) = 9,10,9,8
N pt DFTs of 2 real sequences can be found using a single DFT
If g(n) & h(n) are two sequences then let x(n) = g(n) +j h(n)
G(k) = ½ (X(k) + X*(k))
H(k) = 1/2j (X(K) +X*(k))

Page 25
Digital Signal Processing 10EC52

2N pt DFT of a real sequence using a single N pt DFT

Let x(n) be a real sequence of length 2N with y(n) and g(n) denoting its N pt DFT
Let y(n) = x(2n) and g(2n+1)
k
X (k) = Y (k) + WN G (k)
Using DFT to find IDFT
The DFT expression can be used to find IDFT
X(n) = 1/N [DFT(X*(k)]*

(ii)DFT of Real and even sequence.

For a real sequence, if x(n) X(k)
X (N-K) = X* (k)

For a complex sequence

DFT(x*(n)) = X*(N-K)

If x(n) then X(k)

Real and even real and even
Real and odd imaginary and odd
Odd and imaginary real odd
Even and imaginary imaginary and even

Question 3

Distinguish between circular and linear convolution

Solution

1) Circular convolution is used for periodic and finite signals while linear convolution is
used for aperiodic and infinite signals.
2) In linear convolution we convolved one signal with another signal where as in circular
convolution the same convolution is done but in circular pattern depending upon the
samples of the signal
3) Shifts are linear in linear in linear convolution, whereas it is circular in circular
convolution.

Question 4

Page 26
Digital Signal Processing 10EC52

Solution(a)

Solution(b)

Solution(c)

Solution(d)

Page 27
Digital Signal Processing 10EC52

Question 5

Solution

Question 6

Solution

Page 28
Digital Signal Processing 10EC52

Page 29
Digital Signal Processing 10EC52

UNIT 3

FAST FOURIER TRANSFORMS (FFT) ALOGORITHMS

CONTENTS:-

1. FAST-FOURIER-TRANSFORM (FFT) ALGORITHMS

2. DIRECT COMPUTATION OF DFT,

3. NEED FOR EFFICIENT COMPUTATION OF THE DFT (FFT ALGORITHMS).

RECOMMENDED READINGS

1. DIGITAL SIGNAL PROCESSING – PRINCIPLES ALGORITHMS & APPLICATIONS, PROAKIS &

MONALAKIS, PEARSON EDUCATION, 4TH EDITION, NEW DELHI, 2007.

2. DISCRETE TIME SIGNAL PROCESSING, OPPENHEIM & SCHAFFER, PHI, 2003.

3. DIGITAL SIGNAL PROCESSING, S. K. MITRA, TATA MC-GRAW HILL, 2ND EDITION, 2004.

Page 30
Digital Signal Processing 10EC52

UNIT 3
FAST-FOURIER-TRANSFORM (FFT) ALGORITHMS

3.1 Direct Computation of DFT

The problem:
Given signal samples: x[0], . . . , x[N - 1] (some of which may be zero), develop a procedure
to compute

for k = 0, . . . , N - 1 where

We would like the procedure to be fast, simple, and accurate. Fast is the most important, so we will
sacrifice simplicity for speed, hopefully with minimal loss of accuracy

3.2 Need for efficient computation of DFT (FFT Algorithms)

Let us start with the simple way. Assume that has been precompiled and stored in a

table for the N of interest. How big should the table be? is periodic in m with period N,
so we just need to tabulate the N values:

(Possibly even less since Sin is just Cos shifted by a quarter periods, so we could save just Cos
when N is a multiple of 4.)
Why tabulate? To avoid repeated function calls to Cos and sin when computing the DFT. Now
we can compute each X[k] directly form the formula as follows

For each value of k, there are N complex multiplications, and (N-1) complex additions. There
are N values of k, so the total number of complex operations is

Page 31
Digital Signal Processing 10EC52

Complex multiplies require 4 real multiplies and 2 real additions, whereas complex additions
require just 2 real additions N2 complex multiplies are the primary concern.
N2 increases rapidly with N, so how can we reduce the amount of computation? By exploiting
the following properties of W:

The first and third properties hold for even N, i.e., when 2 is one of the prime factors of N.
There are related properties for other prime factors of N.

Divide and conquer approach

We have seen in the preceding sections that the DFT is a very computationally
intensive operation. In 1965, Cooley and Tukey published an algorithm that could be used to
compute the DFT much more efficiently. Various forms of their algorithm, which came to be
known as the Fast Fourier Transform (FFT), had actually been developed much earlier by
other mathematicians (even dating back to Gauss). It was their paper, however, which
stimulated a revolution in the field of signal processing.

It is important to keep in mind at the outset that the FFT is not a new transform. It is
simply a very efficient way to compute an existing transform, namely the DFT. As we saw, a
straight forward implementation of the DFT can be computationally expensive because the
number of multiplies grows as the square of the input length (i.e. N2 for an N point DFT). The
FFT reduces this computation using two simple but important concepts. The first concept,
known as divide-and-conquer, splits the problem into two smaller problems. The second
concept, known as recursion, applies this divide-and-conquer method repeatedly until the
problem is solved.

Page 32
Digital Signal Processing 10EC52

Recommended Questions with solutions

Question1

Solution:-

Question 2

Solution:-

Page 33
Digital Signal Processing 10EC52

Question 3

Solution:-

Question 4

Page 34
Digital Signal Processing 10EC52

Solution:- (a)

(b)

Page 35
Digital Signal Processing 10EC52

UNIT 4

FAST FOURIER TRANSFORMS (FFT) ALOGORITHMS

CONTENTS:-

1. RADIX-2 FFT ALGORITHM FOR THE COMPUTATION OF DFT AND IDFT

2. DECIMATION-IN-TIME AND DECIMATION-IN-FREQUENCY ALGORITHMS.

3. GOERTZEL ALGORITHM,

4. CHIRP-Z TRANSFORM

RECOMMENDED READINGS

1. DIGITAL SIGNAL PROCESSING – PRINCIPLES ALGORITHMS & APPLICATIONS, PROAKIS &

MONALAKIS, PEARSON EDUCATION, 4TH EDITION, NEW DELHI, 2007.

2. DISCRETE TIME SIGNAL PROCESSING, OPPENHEIM & SCHAFFER, PHI, 2003.

3. DIGITAL SIGNAL PROCESSING, S. K. MITRA, TATA MC-GRAW HILL, 2ND EDITION, 2004.

Page 36
Digital Signal Processing 10EC52

UNIT 4
RADIX-2 FFT ALGORITHM FOR THE COMPUTATION OF DFT AND
IDFT
4.1 Introduction:
Standard frequency analysis requires transforming time-domain signal to frequency
domain and studying Spectrum of the signal. This is done through DFT computation. N-point
DFT computation results in N frequency components. We know that DFT computation
through FFT requires N/2 log2N complex multiplications and N log2N additions. In certain
applications not all N frequency components need to be computed (an application will be
discussed). If the desired number of values of the DFT is less than 2 log 2N than direct
computation of the desired values is more efficient that FFT based computation.
4.2 Radix-2 FFT
Useful when N is a power of 2: N = rv for integers r and v. „r‟ is called the radix, which
comes from the Latin word meaning .a root, and has the same origins as the word radish.
When N is a power of r = 2, this is called radix-2, and the natural .divide and conquer
approach. is to split the sequence into two sequences of length N=2. This is a very clever trick
that goes back many years.
4.2.1 Decimation in time

Fig 4.1 First step in Decimation-in-time domain Algorithm

Page 37
Digital Signal Processing 10EC52

Page 38
Digital Signal Processing 10EC52

4.2.2 Decimation-in-frequency Domain

Another important radix-2 FFT algorithm, called decimation-in-frequency algorithm is

obtained by using divide-and-conquer approach with the choice of M=2 and L= N/2.This
choice of data implies a column-wise storage of the input data sequence. To derive the
algorithm, we begin by splitting the DFT formula into two summations, one of which involves
the sum over the first N/2 data points and the second sum involves the last N/2 data points.
Thus we obtain

Now, let us split X(k) into the even and odd-numbered samples. Thus we obtain

Page 39
Digital Signal Processing 10EC52

Fig 4.2 Shuffling of Data and Bit reversal

The computation of the sequences g1 (n) and g2 (n) and subsequent use of these
sequences to compute the N/2-point DFTs depicted in fig we observe that the basic
computation in this figure involves the butterfly operation.

Page 40
Digital Signal Processing 10EC52

The computation procedure can be repeated through decimation of the N/2-point DFTs,
X(2k) and X(2k+1). The entire process involves v = log 2 N of decimation, where each stage
involves N/2 butterflies of the type shown in figure 4.3.

Fig 4.3 First step in Decimation-in-time domain Algorithm

Page 41
Digital Signal Processing 10EC52

Fig 4.4 N=8 point Decimation-in-frequency domain Algorithm

4.2 Example: DTMF – Dual Tone Multi frequency
This is known as touch-tone/speed/electronic dialing, pressing of each button generates a
unique set of two-tone signals, called DTMF signals. These signals are processed at exchange
to identify the number pressed by determining the two associated tone frequencies. Seven
frequencies are used to code the 10 decimal digits and two special characters (4x3 array)

Page 42
Digital Signal Processing 10EC52

In this application frequency analysis requires determination of possible seven (eight)

DTMF fundamental tones and their respective second harmonics .For an 8 kHz sampling freq,
the best value of the DFT length N to detect the eight fundamental DTMF tones has been
found to be 205 .Not all 205 freq components are needed here, instead only those
corresponding to key frequencies are required. FFT algorithm is not effective and efficient in
this application. The direct computation of the DFT which is more effective in this application
is formulated as a linear filtering operation on the input data sequence.

This algorithm is known as Goertzel Algorithm

This algorithm exploits periodicity property of the phase factor. Consider the DFT definition

N 1
X (k ) x(n)WNnk (1)
n 0

WN kN
Since is equal to 1, multiplying both sides of the equation by this results in;
N 1 N 1
X (k ) WN kN x(m)WNmk x(m)WN k ( N m)
(2)
m 0 m 0

y k ( n) x(n) hk (n)
This is in the form of a convolution
N 1
y k ( n) x(m)WN k ( n m)
(3)
m 0

hk (n) WN kn u(n) (4)

Where yk(n) is the out put of a filter which has impulse response of hk(n) and input x(n).

The output of the filter at n = N yields the value of the DFT at the freq ωk = 2πk/N

The filter has frequency response given by

1
H k ( z) (6)
1 WN k z 1

The above form of filter response shows it has a pole on the unit circle at the frequency ωk =
2πk/N.

Entire DFT can be computed by passing the block of input data into a parallel bank of N
single-pole filters (resonators)

Page 43
Digital Signal Processing 10EC52

The above form of filter response shows it has a pole on the unit circle at the frequency ωk =
2πk/N.

Entire DFT can be computed by passing the block of input data into a parallel bank of N
single-pole filters (resonators)

1.3 Difference Equation implementation of filter:

From the frequency response of the filter (eq 6) we can write the following difference
equation relating input and output;
Yk ( z ) 1
H k ( z) k 1
X ( z ) output
The desired 1 W N is z X(k) = yk(n) for k = 0,1,…N-1. The phase factor appearing in the
k
difference equation can be
y k (n) W N y k (n 1) x(n)computed y konce
( 1) and0 stored. (7)

The form shown in eq (7) requires complex multiplications which can be avoided
doing suitable modifications (divide and multiply by 1 WNk z 1 ). Then frequency response of
the filter can be alternatively expressed as
1 WNk z 1
H k ( z) 1 2
(8)
1 2 cos(2 k / N ) z z

This is second –order realization of the filter (observe the denominator now is a second-order
expression). The direct form realization of the above is given by

vk (n) 2 cos(2 k / N )vk (n 1) vk (n 2) x(n) (9)

k
yk (n) vk (n) W v (n 1)
N k vk ( 1) vk ( 2) 0 (10)

The recursive relation in (9) is iterated for n = 0,1,……N, but the equation in (10) is
computed only once at time n =N. Each iteration requires one real multiplication and two
additions. Thus, for a real input sequence x(n) this algorithm requires (N+1) real
multiplications to yield X(k) and X(N-k) (this is due to symmetry). Going through the Goertzel

Page 44
Digital Signal Processing 10EC52

algorithm it is clear that this algorithm is useful only when M out of N DFT values need to be
computed where M≤ 2log2N, Otherwise, the FFT algorithm is more efficient method. The
utility of the algorithm completely depends on the application and number of frequency
components we are looking for.

4.2. Chirp z- Transform

4.2.1 Introduction:

Computation of DFT is equivalent to samples of the z-transform of a finite-length

sequence at equally spaced points around the unit circle. The spacing between the samples is
given by 2π/N. The efficient computation of DFT through FFT requires N to be a highly
composite number which is a constraint. Many a times we may need samples of z-transform
on contours other than unit circle or we my require dense set of frequency samples over a
small region of unit circle. To understand these let us look in to the following situations:

1. Obtain samples of z-transform on a circle of radius „a‟ which is concentric to unit circle
The possible solution is to multiply the input sequence by a -n
2. 128 samples needed between frequencies ω = -π/8 to +π/8 from a 128 point sequence
From the given specifications we see that the spacing between the frequency samples is
π/512 or 2π/1024. In order to achieve this freq resolution we take 1024- point FFT of
the given 128-point seq by appending the sequence with 896 zeros. Since we need
only 128 frequencies out of 1024 there will be big wastage of computations in this
scheme.

For the above two problems Chirp z-transform is the alternative.

Chirp z- transform is defined as:

N 1
X ( zk ) x ( n) z k n k 0,1,...... L 1 (11)
n 0

Page 45
Digital Signal Processing 10EC52

Where zk is a generalized contour. Zk is the set of points in the z-plane falling on an arc which
begins at some point z0 and spirals either in toward the origin or out away from the origin such
that the points {zk}are defined as,

zk r0e j 0 ( R0e j 0 ) k k 0,1,.... L 1 (12)

Note that,

a. if R0< 1 the points fall on a contour that spirals toward the origin

b. If R0 > 1 the contour spirals away from the origin

c. If R0= 1 the contour is a circular arc of radius

d.If r0=1 and R0=1 the contour is an arc of the unit circle.

(Additionally this contour allows one to compute the freq content of the sequence x(n) at
dense set of L frequencies in the range covered by the arc without having to compute a large
DFT (i.e., a DFT of the sequence x(n) padded with many zeros to obtain the desired resolution
in freq.))

e. If r0= R0=1 and θ0=0 Φ0=2π/N and L = N the contour is the entire unit circle similar to the
standard DFT. These conditions are shown in the following diagram.

Page 46
Digital Signal Processing 10EC52

Substituting the value of zk in the expression of X(zk)

N 1 N 1
X ( zk ) x ( n) z k n x(n)(r0e j 0 ) nW nk
(13)
n 0 n 0

where
W R0 e j 0
(14)

4.2.2 Expressing computation of X(zk) as linear filtering operation:

By substitution of
1 2
nk ( n k 2 ( k n) 2 ) (15)
2
we can express X(zk) as

k2 /2
X ( zk ) W y (k ) y(k ) / h(k ) k 0,1,......... .L 1 (16)
Where
2
n2 / 2
h( n) W n /2
g (n) x(n)(r0e j 0 ) n W
N 1
y (k ) g ( n) h( k n) (17)
n 0
both g(n) and h(n) are complex valued sequences

4.2.3 Why it is called Chirp z-transform?

If R0 =1, then sequence h(n) has the form of complex exponential with argument ωn =
n2Φ0/2 = (n Φ0/2) n. The quantity (n Φ0/2) represents the freq of the complex exponential

Page 47
Digital Signal Processing 10EC52

signal, which increases linearly with time. Such signals are used in radar systems are called
chirp signals. Hence the name chirp z-transform.

4.2.4 How to Evaluate linear convolution of eq (17)

1. Can be done efficiently with FFT
2. The two sequences involved are g(n) and h(n). g(n) is finite length seq of length N and
h(n) is of infinite duration, but fortunately only a portion of h(n) is required to compute
L values of X(z), hence FFT could be still be used.
3. Since convolution is via FFT, it is circular convolution of the N-point seq g(n) with an
M- point section of h(n) where M > N

4. The concepts used in overlap –save method can be used

5. While circular convolution is used to compute linear convolution of two sequences we
know the initial N-1 points contain aliasing and the remaining points are identical to
the result that would be obtained from a linear convolution of h(n) and g(n), In view of
this the DFT size selected is M = L+N-1 which would yield L valid points and N-1
points corrupted by aliasing. The section of h(n) considered is for –(N-1) ≤ n≤ (L-1)
yielding total length M as defined
6. The portion of h(n) can be defined in many ways, one such way is,

h1(n) = h(n-N+1) n = 0,1,…..M-1

7. Compute H1(k) and G(k) to obtain

Y1(k) = G(K)H1(k)
8. Application of IDFT will give y1(n), for

Page 48
Digital Signal Processing 10EC52

n =0,1,…M-1. The starting N-1 are discarded and desired values are y1(n) for

N-1 ≤n ≤ M-1 which corresponds to the range 0 ≤n ≤ L-1 i.e.,

y(n)= y1(n+N-1) n=0,1,2,…..L-1

9. Alternatively h2(n) can be defined as
h2 (n) h( n ) 0 n L 1
h(n ( N L 1)) L n M 1
10. Compute Y2(k) = G(K)H2(k), The desired values of y2(n) are in the range

0 ≤n ≤L-1 i.e.,

y(n) = y2(n) n=0,1,….L-1

11. Finally, the complex values X(zk) are computed by dividing y(k) by h(k)
For k =0,1,……L-1

4.3 Computational complexity

In general the computational complexity of CZT is of the order of M log 2M complex
multiplications. This should be compared with N.L which is required for direct evaluation.
If L is small direct evaluation is more efficient otherwise if L is large then CZT is more
efficient.

4.3.1 Advantages of CZT

a. Not necessary to have N =L

b.Neither N or L need to be highly composite

c.The samples of Z transform are taken on a more general contour that includes the unit
circle as a special case.

4.4 Example to understand utility of CZT algorithm in freq analysis

(ref: DSP by Oppenheim Schaffer)

CZT is used in this application to sharpen the resonances by evaluating the z-transform
off the unit circle. Signal to be analyzed is a synthetic speech signal generated by exciting a
five-pole system with a periodic impulse train. The system was simulated to correspond to a
sampling freq. of 10 kHz. The poles are located at center freqs of 270,2290,3010,3500 & 4500
Hz with bandwidth of 30, 50, 60,87 & 140 Hz respectively.

Page 49
Digital Signal Processing 10EC52

Solution: Observe the pole-zero plots and corresponding magnitude frequency response for
different choices of |w|. The following observations are in order:

• The first two spectra correspond to spiral contours outside the unit circle with a resulting
broadening of the resonance peaks
• |w| = 1 corresponds to evaluating z-transform on the unit circle
• The last two choices correspond to spiral contours which spirals inside the unit circle and
close to the pole locations resulting in a sharpening of resonance peaks.

Page 50
Digital Signal Processing 10EC52

4.5 Implementation of CZT in hardware to compute the DFT signals

The block schematic of the CZT hardware is shown in down figure. DFT computation
requires r0 =R0 =1, θ0 = 0 Φ0 = 2π/N and L = N.

The cosine and sine sequences in h(n) needed for pre multiplication and post multiplication are
usually stored in a ROM. If only magnitude of DFT is desired, the post multiplications are
unnecessary,

In this case |X(zk)| = |y(k)| k =0,1,….N-1

Page 51
Digital Signal Processing 10EC52

Recommended Questions with solutions

Question 1

Solution:-

Page 52
Digital Signal Processing 10EC52

Question 2

Solution :- There are 20 real , non trial multiplications

Figure 4.1 DIF Algorithm for N=16

Page 53
Digital Signal Processing 10EC52

Question 3

Solution:-

Question 4

Solution:-

Page 54
Digital Signal Processing 10EC52

Question 5

Solution:-

Question 6

Solution:-

This can be viewed as the convolution of the N-length sequence x(n) with implulse
response of a linear filter

Page 55
Digital Signal Processing 10EC52

Page 56
Digital Signal Processing 10EC52

UNIT 5

IIR FILTER DESIGN

CONTENTS:-

1. IIR FILTER DESIGN:

2. CHARACTERISTICS OF COMMONLY USED ANALOG FILTERS –

3. BUTTERWORTH AND CHEBYSHEVE FILTERS,

4. ANALOG TO ANALOG FREQUENCY TRANSFORMATIONS.

RECOMMENDED READINGS

1. DIGITAL SIGNAL PROCESSING – PRINCIPLES ALGORITHMS & APPLICATIONS, PROAKIS &

MONALAKIS, PEARSON EDUCATION, 4TH EDITION, NEW DELHI, 2007.

2. DISCRETE TIME SIGNAL PROCESSING, OPPENHEIM & SCHAFFER, PHI, 2003.

3. DIGITAL SIGNAL PROCESSING, S. K. MITRA, TATA MC-GRAW HILL, 2ND EDITION, 2004.

Page 57
Digital Signal Processing 10EC52

Unit 5
Design of IIR Filters
5.1 Introduction
A digital filter is a linear shift-invariant discrete-time system that is realized using finite
precision arithmetic. The design of digital filters involves three basic steps:

 The specification of the desired properties of the system.

 The approximation of these specifications using a causal discrete-time system.
 The realization of these specifications using finite precision arithmetic.

These three steps are independent; here we focus our attention on the second step. The
desired digital filter is to be used to filter a digital signal that is derived from an analog signal
by means of periodic sampling. The specifications for both analog and digital filters are often
given in the frequency domain, as for example in the design of low pass, high pass, band pass
and band elimination filters.
Given the sampling rate, it is straight forward to convert from frequency specifications
on an analog filter to frequency specifications on the corresponding digital filter, the analog
frequencies being in terms of Hertz and digital frequencies being in terms of radian frequency
or angle around the unit circle with the point Z=-1 corresponding to half the sampling
frequency. The least confusing point of view toward digital filter design is to consider the filter
as being specified in terms of angle around the unit circle rather than in terms of analog
frequencies.

Page 58
Digital Signal Processing 10EC52

Figure 5.1: Tolerance limits for approximation of ideal low-pass filter

A separate problem is that of determining an appropriate set of specifications on the

digital filter. In the case of a low pass filter, for example, the specifications often take the form
of a tolerance scheme, as shown in Fig. 5.1.

Many of the filters used in practice are specified by such a tolerance scheme, with no
constraints on the phase response other than those imposed by stability and causality
requirements; i.e., the poles of the system function must lie inside the unit circle. Given a set
of specifications in the form of Fig. 5.1, the next step is to and a discrete time linear system
whose frequency response falls within the prescribed tolerances. At this point the filter design
problem becomes a problem in approximation. In the case of infinite impulse response (IIR)
filters, we must approximate the desired frequency response by a rational function, while in the
finite impulse response (FIR) filters case we are concerned with polynomial approximation.

5.1 Design of IIR Filters from Analog Filters:

Page 59
Digital Signal Processing 10EC52

The traditional approach to the design of IIR digital filters involves the transformation
of an analog filter into a digital filter meeting prescribed specifications. This is a reasonable
approach because:
 The art of analog filter design is highly advanced and since useful results can be
achieved, it is advantageous to utilize the design procedures already developed for
analog filters.
 Many useful analog design methods have relatively simple closed-form design
formulas.
Therefore, digital filter design methods based on analog design formulas are rather simple to
implement. An analog system can be described by the differential equation

And the corresponding rational function is

The corresponding description for digital filters has the form

and the rational function

In transforming an analog filter to a digital filter we must therefore obtain either H(z)
or h(n) (inverse Z-transform of H(z) i.e., impulse response) from the analog filter design. In
such transformations, we want the imaginary axis of the S-plane to map into the nit circle of
the Z-plane, a stable analog filter should be transformed to a stable digital filter. That is, if the
analog filter has poles only in the left-half of S-plane, then the digital filter must have poles
only inside the unit circle. These constraints are basic to all the techniques discussed here.

Page 60
Digital Signal Processing 10EC52

5.2 Characteristics of Commonly Used Analog Filters:

From the previous discussion it is clear that, IIT digital filters can be obtained by
beginning with an analog filter. Thus the design of a digital filter is reduced to designing an
appropriate analog filter and then performing the conversion from Ha(s) to H (z). Analog filter
design is a well - developed field, many approximation techniques, viz., Butterworth,
Chebyshev, Elliptic, etc., have been developed for the design of analog low
pass filters. Our discussion is limited to low pass filters, since, frequency transformation can
be applied to transform a designed low pass filter into a desired high pass, band pass and band
stop filters.

5.2.1 Butterworth Filters:

Low pass Butterworth filters are all - pole filters with monotonic frequency response in
both pass band and stop band, characterized by the magnitude - squared frequency response

Where, N is the order of the filter, Ώc is the -3dB frequency, i.e., cutoff frequency, Ώp is the
pass band edge frequency and 1= (1 /1+ε2 ) is the band edge value of │Ha(Ώ)│2. Since the
product Ha(s) Ha(-s) and evaluated at s = jΏ is simply equal to │Ha(Ώ)│2, it follows that

The poles of Ha(s)Ha(-s) occur on a circle of radius Ώc at equally spaced points. From Eq.
(5.29), we find the pole positions as the solution of

And hence, the N poles in the left half of the s-plane are

Page 61
Digital Signal Processing 10EC52

Note that, there are no poles on the imaginary axis of s-plane, and for N odd there will
be a pole on real axis of s-plane, for N even there are no poles even on real axis of s-plane.
Also note that all the poles are having conjugate symmetry. Thus the design methodology to
design a Butterworth low pass filter with δ2 attenuation at a specified frequency Ώs is Find N,

Where by definition, δ2 = 1/√1+δ2. Thus the Butterworth filter is completely

characterized by the parameters N, δ2, ε and the ratio Ώs/Ώp or Ώc.Then, from Eq. (5.31) find
the pole positions Sk; k = 0,1, 2,……..(N-1). Finally the analog filter is given by

5.2.2 Chebyshev Filters:

There are two types of Chebyshev filters. Type I Chebyshev filters are all-pole filters
that exhibit equiripple behavior in the pass band and a monotonic characteristic in the stop
band. On the other hand, type II Chebyshev filters contain both poles and zeros and exhibit a
monotonic behavior in the pass band and an equiripple behavior in the stop band. The zeros of
this class of filters lie on the imaginary axis in the s-plane. The magnitude squared of the
frequency response characteristic of type I Chebyshev filter is given as

Where ε is a parameter of the filter related to the ripple in the pass band as shown in Fig.
(5.7), and TN is the Nth order Chebyshev polynomial defined as

The Chebyshev polynomials can be generated by the recursive equation

Page 62
Digital Signal Processing 10EC52

Where T0(x) = 1 and T1(x) = x.

At the band edge frequency Ώ= Ώp, we have

Figure 5.2: Type I Chebysehev filter characteristic

Or equivalently

Where δ1 is the value of the pass band ripple.

The poles of Type I Chebyshev filter lie on an ellipse in the s-plane with major axis

And minor axis

Where β is related to ε according to the equation

Page 63
Digital Signal Processing 10EC52

The angular positions of the left half s-plane poles are given by

Then the positions of the left half s-plane poles are given by

Where ζk = r2 Cos φk and Ώk = r1 Sinφk. The order of the filter is obtained from

Where, by definition δ2 = 1/√1+δ2.

Finally, the Type I Chebyshev filter is given by

A Type II Chebyshev filter contains zero as well as poles. The magnitude squared response is
given as

Where TN(x) is the N-order Chebyshev polynomial. The zeros are located on the imaginary
axis at the points

and the left-half s-plane poles are given

Page 64
Digital Signal Processing 10EC52

Where

and

Finally, the Type II Chebyshev filter is given by

The other approximation techniques are elliptic (equiripple in both passband and
stopband) and Bessel (monotonic in both passband and stopband).

5.3 Analog to Analog Frequency Transforms

Frequency transforms are used to transform lowpass prototype filter to other filters like
highpass or bandpass or bandstop filters. One possibility is to perform frequency transform in
the analog domain and then convert the analog filter into a corresponding digital filter by a
mapping of the s-plane into z-plane. An alternative approach is to convert the analog lowpass
filter into a lowpass digital filter and then to transform the lowpass digital filter into the
desired digital filter by a digital transformation.

Suppose we have a lowpass filter with pass edge ΩP and if we want convert that into
another lowpass filter with pass band edge Ω‟ P then the transformation used is

To convert low pass filter into highpass filter the transformation used is
Page 65
Digital Signal Processing 10EC52

Thus we obtain

The filter function is

Page 66
 10EC52
Digital Signal Processing
Recommended Questions with answers
Question 1

I Design a digital filter to satisfy the following characteristics.

 -3dB cutoff frequency of 0:5_ rad.
 Magnitude down at least 15dB at 0:75_ rad.
 Monotonic stop band and pass band Using
 Impulse invariant technique
 Approximation of derivatives
 Bilinear transformation technique

Figure 5.8: Frequency response plot of the example

Solution:-
a) Impulse Invariant Technique
From the given digital domain frequency, _nd the corresponding analog domain frequencies.

Where T is the sampling period and 1/T is the sampling frequency and it always corresponds
to 2Π radians in the digital domain. In this problem, let us assume T = 1sec.
Then Ώc = 0:5Π and Ώs = 0:75Π
Let us find the order of the desired filter using

Page 67
Digital Signal Processing 10EC52

Where δ2 is the gain at the stop band edge frequency ωs.

Order of filter N =5.

Then the 5 poles on the Butterworth circle of radius Ώc = 0:5 Π are given by

Then the filter transfer function in the analog domain is

where Ak's are partial fractions coefficients of Ha(s).

Finally, the transfer function of the digital filter is

Page 68
Digital Signal Processing 10EC52

b)

c) For the bilinear transformation technique, we need to pre-warp the digital frequencies
into corresponding analog frequencies.

Then the order of the filter

The pole locations on the Butterworth circle with radius Ώc = 2 are

Then the filter transfer function in the analog domain is

Finally, the transfer function of the digital filter is

Page 69
Digital Signal Processing 10EC52

Question 2
Design a digital filter using impulse invariant technique to satisfy following
characteristics
(i) Equiripple in pass band and monotonic in stop band
(ii) -3dB ripple with pass band edge frequency at 0:5П radians.
(iii) Magnitude down at least 15dB at 0:75 П radians.

Solution: Assuming T=1, Ώ= 0:5 П and s = 0:75 П

The order of desired filter is

Page 70
Digital Signal Processing 10EC52

Page 71
Digital Signal Processing 10EC52

Page 72
Digital Signal Processing 10EC52

Question 3

Solution:-
For the design specifications we have

Page 73
Digital Signal Processing 10EC52

Question 4

Solution:-

Page 74
Digital Signal Processing 10EC52

UNIT 6

FIR FILTER DESIGN

CONTENTS:-

1. INTRODUCTION TO FIR FILTERS,

2. DESIGN OF FIR FILTERS USING

 RECTANGULAR

 HAMMING

 BARTLET

 KAISER WINDOWS,

3. FIR FILTER DESIGN USING FREQUENCY SAMPLING TECHNIQUE

RECOMMENDED READINGS

4. DIGITAL SIGNAL PROCESSING – PRINCIPLES ALGORITHMS & APPLICATIONS, PROAKIS &

MONALAKIS, PEARSON EDUCATION, 4TH EDITION, NEW DELHI, 2007.

5. DISCRETE TIME SIGNAL PROCESSING, OPPENHEIM & SCHAFFER, PHI, 2003.

6. DIGITAL SIGNAL PROCESSING, S. K. MITRA, TATA MC-GRAW HILL, 2ND EDITION, 2004.

Page 75
Digital Signal Processing 10EC52

UNIT 6
Design of FIR Filters

6.1 Introduction:
Two important classes of digital filters based on impulse response type are

Finite Impulse Response (FIR)

Infinite Impulse Response (IIR)

The filter can be expressed in two important forms as:

1 ) System function representation;

M
k
bk z
k 0
H ( z) N
(1)
k
1 ak z
k 1

2) Difference Equation representation;

N M
ak y ( n k ) bk x(n k ) (2)
k 0 k 0

Each of this form allows various methods of implementation. The eq (2) can be viewed
as a computational procedure (an algorithm) for determining the output sequence y(n) of the
system from the input sequence x(n). Different realizations are possible with different
arrangements of eq (2)

The major issues considered while designing a digital filters are :

• Realiability (causal or non causal)
• Stability (filter output will not saturate)
• Sharp Cutoff Characteristics
• Order of the filter need to be minimum (this leads to less delay)
• Generalized procedure (having single procedure for all kinds of filters)
• Linear phase characteristics

Page 76
Digital Signal Processing 10EC52

The factors considered with filter implementation are ,

a. It must be a simple design

b. There must be modularity in the implementation so that any order filter can be obtained with
lower order modules.

c. Designs must be as general as possible. Having different design procedures for different
types of filters( high pass, low pass,…) is cumbersome and complex.

d. Cost of implementation must be as low as possible

e. The choice of Software/Hardware realization

6.2 Features of IIR:

The important features of this class of filters can be listed as:
• Out put is a function of past o/p, present and past i/p‟s
• It is recursive in nature
• It has at least one Pole (in general poles and zeros)
• Sharp cutoff chas. is achievable with minimum order
• Difficult to have linear phase chas over full range of freq.
• Typical design procedure is analog design then conversion from analog to digital

6.3 Features of FIR : The main features of FIR filter are,

• They are inherently stable
• Filters with linear phase characteristics can be designed
• Simple implementation – both recursive and nonrecursive structures possible
• Free of limit cycle oscillations when implemented on a finite-word length digital system

6.3.1 Disadvantages:
• Sharp cutoff at the cost of higher order
• Higher order leading to more delay, more memory and higher cost of implementation

6.4 Importance of Linear Phase:

The group delay is defined as
d ( )
g
d
which is negative differential of phase function.

Nonlinear phase results in different frequencies experiencing different delay and arriving
at different time at the receiver. This creates problems with speech processing and data

Page 77
Digital Signal Processing 10EC52

communication applications. Having linear phase ensures constant group delay for all
frequencies.
The further discussions are focused on FIR filter.
6.5 Examples of simple FIR filtering operations: 1.Unity Gain Filter

y(n)=x(n)

2. Constant gain filter

y(n)=Kx(n)

3. Unit delay filter

y(n)=x(n-1)

4.Two - term Difference filter

y(n) = x(n)-x(n-1)

5. Two-term average filter

y(n) = 0.5(x(n)+x(n-1))

6. Three-term average filter (3-point moving average filter)

y(n) = 1/3[x(n)+x(n-1)+x(n-2)]

7. Central Difference filter

y(n)= 1/2[ x(n) – x(n-2)]

When we say Order of the filter it is the number of previous inputs used to compute the
current output and Filter coefficients are the numbers associated with each of the terms x(n),
x(n-1),.. etc
The table below shows order and filter coefficients of above simple filter types:

Page 78
Digital Signal Processing 10EC52

Ex. order a0 a1 a2

1 0 1 - -

2 0 K - -

3 1 0 1 -

4(HP) 1 1 -1 -

5(LP) 1 1/2 1/2 -

6(LP) 2 1/3 1/3 1/3

7(HP) 2 1/2 0 -1/2

6.6 Design of FIR filters:

The section to follow will discuss on design of FIR filter. Since linear phase can be
achieved with FIR filter we will discuss the conditions required to achieve this.

6.6.1 Symmetric and Antisymmetric FIR filters giving out Linear Phase characteristics:

Symmetry in filter impulse response will ensure linear phase

An FIR filter of length M with i/p x(n) & o/p y(n) is described by the difference equation:

M 1
y(n)= b0 x(n) + b1 x(n-1)+…….+b M-1 x(n-(M-1)) = bk x(n k ) -(1)
k 0

Alternatively. it can be expressed in convolution form

M 1
y ( n) h( k ) x ( n k ) - (2)
k 0

i.e bk= h(k), k=0,1,…..M-1

Filter is also characterized by

Page 79
Digital Signal Processing 10EC52
M 1
H ( z) h( k ) z k
-(3) polynomial of degree M-1 in the variable z-1. The roots of this
k 0
polynomial constitute zeros of the filter.

An FIR filter has linear phase if its unit sample response satisfies the condition
h(n)= ± h(M-1-n) n=0,1,…….M-1 -(4)

Incorporating this symmetry & anti symmetry condition in eq 3 we can show linear phase
chas of FIR filters

1 2 ( M 2) ( M 1)
H ( z) h(0) h(1) z h(2) z .......... . h( M 2) z h( M 1) z

If M is odd
M 1 M 1 M 3
1 M 1 (
2
) M 1 (
2
) M 3 (
2
)
H ( z) h(0) h(1) z .......... h( )z h( )z h( )z .......... .
2 2 2
( M 2) ( M 1)
h( M 2) z h( M 1) z

M 1 M 1 M 3 M 1
(
2
) (
2
) (
2
) M 1 M 1 1 M 3 2
(
2
)
z h(0) z h(1) z .......... .. h( ) h( )z h( )z .....h( M 1) z
2 2 2
Applying symmetry conditions for M odd

h(0) h( M 1)
h(1) h( M 2)
.
.
M 1 M 1
h( ) h( )
2 2
M 1 M 3
h( ) h( )
2 2
.
.
h( M 1) h(0)

Page 80
Digital Signal Processing 10EC52
M 3
M 1 2
( ) M 1
H ( z) z 2
h( ) h(n){z ( M 1 2n) / 2
z ( M 1 2n) / 2
}
2 n 0

similarly for M even

M
M 1 1
( ) 2
H ( z) z 2
h(n){z ( M 1 2n) / 2
z ( M 1 2n) / 2
}
n 0

6.6.2 Frequency response:

If the system impulse response has symmetry property (i.e.,h(n)=h(M-1-n)) and M is odd
H (e j ) e j ( ) | H r (e j ) | where

M 3
2
M 1 M 1
H r (e j ) h( ) 2 h(n) cos ( n)
2 n 0 2

M 1
( ) ( ) if | H r (e j ) | 0
2
M 1
( ) if | H r (e j ) | 0
2
In case of M even the phase response remains the same with magnitude response expressed as

M
1
2
M 1
H r (e j ) 2 h(n) cos ( n)
n 0 2

If the impulse response satisfies anti symmetry property (i.e., h(n)=-h(M-1-n))then for
M odd we will have
M 1 M 1 M 1
h( ) h( ) i.e., h( ) 0
2 2 2
M 3
2
M 1
H r (e j ) 2 h(n) sin ( n)
n 0 2

If M is even then,

Page 81
Digital Signal Processing 10EC52
M
1
2
M 1
H r (e j ) 2 h(n) sin ( n)
n 0 2

In both cases the phase response is given by

M 1
( ) ( ) /2 if | H r (e j ) | 0
2
M 1
( ) 3 / 2 if | H r (e j ) | 0
2

Which clearly shows presence of Linear Phase characteristics.

6.6.3 Comments on filter coefficients:

• The number of filter coefficients that specify the frequency response is (M+1)/2 when is M
odd and M/2 when M is even in case of symmetric conditions
• In case of impulse response antisymmetric we have h(M-1/2)=0 so that there are (M-1/2)
filter coefficients when M is odd and M/2 coefficients when M is even

6.6.5 Choice of Symmetric and antisymmetric unit sample response

When we have a choice between different symmetric properties, the particular one is
picked up based on application for which the filter is used. The following points give an
insight to this issue.
• If h(n)=-h(M-1-n) and M is odd, Hr(w) implies that Hr(0)=0 & Hr(π)=0, consequently not
suited for lowpass and highpass filter. This condition is suited in Band Pass filter design.
• Similarly if M is even Hr(0)=0 hence not used for low pass filter
• Symmetry condition h(n)=h(M-1-n) yields a linear-phase FIR filter with non zero response
at w = 0 if desired.
Looking at these points, antisymmetric properties are not generally preferred.

Page 82
Digital Signal Processing 10EC52

6.6.6 Zeros of Linear Phase FIR Filters:

Consider the filter system function

M 1
n
H ( z) h( n) z
n o

Expanding this equation

H ( z ) h(0) h(1) z 1 h(2) z 2  h( M 2) z ( M 2)

h( M 1) z ( M 1)

sin ce for Linear phase we need

h(n) h( M 1 n) i.e.,
h(0) h( M 1); h(1) h( M 2);...... h( M 1) h(0);
then
1 ( M 2) ( M 1)
H ( z) h( M 1) h( M 2) z ........ h(1) z h(0) z
( M 1)
H ( z) z [h( M 1) z ( M 1)
h( M 2) z ( M 2)
..... h(1) z h(0)]
M 1
( M 1)
H ( z) z [ h(n)( z 1 ) n ] z ( M 1)
H (z 1)
n 0

This shows that if z = z1 is a zero then z=z1 -1 is also a zero

The different possibilities:

1. If z1 = 1 then z1 = z1-1 =1 is also a zero implying it is one zero
2. If the zero is real and |z|<1 then we have pair of zeros
3. If zero is complex and |z|=1then and we again have pair of complex zeros.
4. If zero is complex and |z|≠1 then and we have two pairs of complex zeros

Page 83
Digital Signal Processing 10EC52

The plot above shows distribution of zeros for a Linear – phase FIR filter. As it can be seen
there is pattern in distribution of these zeros.

6.7 Methods of designing FIR filters:

The standard methods of designing FIR filter can be listed as:

1. Fourier series based method

2. Window based method
3. Frequency sampling method

6.7.1 Design of Linear Phase FIR filter based on Fourier Series method:

Motivation: Since the desired freq response Hd(ejω) is a periodic function in ω with
period 2π, it can be expressed as Fourier series expansion

H d (e j ) hd (n)e j n

where hd (n) are fourier series coefficients

1
hd (n) H d (e j )e j n d
2

This expansion results in impulse response coefficients which are infinite in duration and non
causal. It can be made finite duration by truncating the infinite length. The linear phase can be
obtained by introducing symmetric property in the filter impulse response, i.e., h(n) = h(-n). It
can be made causal by introducing sufficient delay (depends on filter length)

6.7.2 Stepwise procedure:

1. From the desired freq response using inverse FT relation obtain hd(n)
2. Truncate the infinite length of the impulse response to finite length with ( assuming
M odd)
h(n) hd (n) for ( M 1) / 2 n ( M 1) / 2
0 otherwise

3. Introduce h(n) = h(-n) for linear phase characteristics

4. Write the expression for H(z); this is non-causal realization
5. To obtain causal realization H‟(z) = z -(M-1)/2 H(z)

Page 84
Digital Signal Processing 10EC52

Exercise Problems
Problem 1 : Design an ideal bandpass filter with a frequency response:

3
H d (e j ) 1 for
4 4
0 otherwise
Find the values of h(n) for M = 11 and plot the frequency response.

1
hd (n) H d (e j )e j n d
2
/4 3 /4
1
e j nd e j nd
2 3 /4 /4

1 3
sin n sin n n
n 4 4
truncating to 11 samples we have h(n) hd (n) for | n | 5
0 otherwise
For n = 0 the value of h(n) is separately evaluated from the basic integration

h(0) = 0.5

Other values of h(n) are evaluated from h(n) expression

h(1)=h(-1)=0

h(2)=h(-2)=-0.3183

h(3)=h(-3)=0

h(4)=h(-4)=0

h(5)=h(-5)=0

The transfer function of the filter is

Page 85
Digital Signal Processing 10EC52

( N 1) / 2
H ( z) h ( 0) [h(n){z n z n }]
n 1

0.5 0.3183( z 2 z 2 )
the transfer function of the realizable filter is
H ' ( z) z 5 [0.5 0.3183( z 2 z 2 )]
0.3183z 3 0.5 z 5
0.3183z 7

the filter coeff are

h ' (0) h' (10) h' (1) h' (9) h ' ( 2) h' (8) h ' ( 4) h ' ( 6) 0
h' (3) h' (7) 0.3183
h' (5) 0.5

The magnitude response can be expressed as

( N 1) / 2
| H (e j ) | a(n) cos n
n 1

comparing this exp with

5
| H (e j ) | | z 5 [h(0) 2 h(n) cos n] |
n 1

We have

a(0)=h(0)

a(1)=2h(1)=0

a(2)=2h(2)=-0.6366

a(3)=2h(3)=0

a(4)=2h(4)=0

a(5)=2h(5)=0

The magnitude response function is

|H(e jω)| = 0.5 – 0.6366 cos 2ω which can plotted for various values of ω

ω in degrees =[0 20 30 45 60 75 90 105 120 135 150 160 180];

Page 86
Digital Signal Processing 10EC52

|H(e jω)| in dBs= [-17.3 -38.17 -14.8 -6.02 -1.74 0.4346 1.11 0.4346 -1.74 -6.02 -14.8 -38.17 -
17.3];

Problem 2: Design an ideal lowpass filter with a freq response

H d (e j ) 1 for
2 2
0 for
2

Find the values of h(n) for N =11. Find H(z). Plot the magnitude response

From the freq response we can determine hd(n),

n
/2 sin
1 2
hd (n) e j nd n and n 0
2 /2
n
Truncating hd(n) to 11 samples

h(0) = 1/2
h(1)=h(-1)=0.3183
h(2)=h(-2)=0
h(3)=h(-3)=-0.106

Page 87
Digital Signal Processing 10EC52

h(4)=h(-4)=0
h(5)=h(-5)=0.06366

The realizable filter can be obtained by shifting h(n) by 5 samples to right h‟(n)=h(n-5)

h‟(n)= [0.06366, 0, -0.106, 0, 0.3183, 0.5, 0.3183, 0, -0.106, 0, 0.06366];

2 4 5 6 8 10
H ' ( z) 0.06366 0.106z 0.3183z 0.5 z 0.3183z 0.106z 0.06366z

Using the result of magnitude response for M odd and symmetry

M 3
2
M 1 M 1
H r (e j ) [ h( ) h(n) cos ( n)]
2 n 0 2
| H r (e j ) | | [0.5 0.6366cos w 0.212 cos 3w 0.127 cos 5w] |

Problem 3 :

Design an ideal band reject filter with a frequency response:

2
H d (e j ) 1 for and
3 3
0 otherwise

Find the values of h(n) for M = 11 and plot the frequency response

Ans:h(n)= [0 -0.1378 0 0.2757 0 0.667 0 0.2757 0 -0.1378 0];

6.8 Window based Linear Phase FIR filter design

The other important method of designing FIR filter is by making use of windows. The
arbitrary truncation of impulse response obtained through inverse Fourier relation can lead to
distortions in the final frequency response.The arbitrary truncation is equivalent to multiplying
infinite length function with finite length rectangular window, i.e.,
h(n) = hd(n) w(n) where w(n) = 1 for n = ±(M-1)/2
The above multiplication in time domain corresponds to convolution in freq domain, i.e.,

Page 88
Digital Signal Processing 10EC52

H( e jω ) = Hd(e jω) * W(e jω ) where W(e jω ) is the FT of window function w(n).

The FT of w(n) is given by

sin( M / 2)
W (e j )
sin( / 2)

The whole process of multiplying h(n) by a window function and its effect in freq domain are
shown in below set of figures.

Page 89
Digital Signal Processing 10EC52

Suppose the filter to be designed is Low pass filter then the convolution of ideal filter freq
response and window function freq response results in distortion in the resultant filter freq
response. The ideal sharp cutoff chars are lost and presence of ringing effect is seen at the band
edges which is referred to Gibbs Phenomena. This is due to main lobe width and side lobes of
the window function freq response.The main lobe width introduces transition band and side
lobes results in rippling characters in pass band and stop band. Smaller the main lobe width
smaller will be the transition band. The ripples will be of low amplitude if the peak of the first
side lobe is far below the main lobe peak.

6.8.1 How to reduce the distortions?

1. Increase length of the window

- as M increases the main lob width becomes narrower, hence the transition band width is
decreased

-With increase in length the side lobe width is decreased but height of each side lobe
increases in such a manner that the area under each sidelobe remains invariant to changes in
M. Thus ripples and ringing effect in pass-band and stop-band are not changed.
2. Choose windows which tapers off slowly rather than ending abruptly - Slow tapering
reduces ringing and ripples but generally increases transition width since main lobe width
of these kind of windows are larger.

Page 90
Digital Signal Processing 10EC52

6.8.2 What is ideal window characteristics?

Window having very small main lobe width with most of the energy contained with it
(i.e.,ideal window freq response must be impulsive).Window design is a mathematical
problem, more complex the window lesser are the distortions. Rectangular window is one of
the simplest window in terms of computational complexity. Windows better than rectangular
window are, Hamming, Hanning, Blackman, Bartlett, Traingular,Kaiser. The different
window functions are discussed in the following sention.

6.8.3 Rectangular window: The mathematical description is given by,

wr (n) 1 for 0 n M 1

6.8.4 Hanning windows:

It is defined mathematically by,

2 n
whan (n) 0.5(1 cos ) for 0 n M 1
M 1

Page 91
Digital Signal Processing 10EC52

6.8.5 Hamming windows:

This window function is given by,

2 n
wham (n) 0.54 0.46 cos for 0 n M 1
M 1

6.8.6 Blackman windows:

This window function is given by,
2 n 4 n
wblk (n) 0.42 0.5 cos 0.08 cos for 0 n M 1
M 1 M 1

Page 92
Digital Signal Processing 10EC52

6.8.7 Bartlett (Triangular) windows:

The mathematical description is given by,

M 1
2|n |
wbart (n) 1 2 for 0 n M 1
M 1

6.8.8 Kaiser windows: The mathematical description is given by,

2
2
M 1 M 1
I0 n
2 2
wk (n) for 0 n M 1
M 1
I0
2

Page 93
Digital Signal Processing 10EC52

Type of window Appr. Transition Peak

width of the main lobe sidelobe (dB)

Rectangular 4π/M -13

Bartlett 8π/M -27

Hanning 8π/M -32

Hamming 8π/M -43

Blackman 12π/M -58

Looking at the above table we observe filters which are mathematically simple do not
offer best characteristics. Among the window functions discussed Kaiser is the most complex
one in terms of functional description whereas it is the one which offers maximum flexibility
in the design.

6.8.9 Procedure for designing linear-phase FIR filters using windows:

1. Obtain hd(n) from the desired freq response using inverse FT relation
2. Truncate the infinite length of the impulse response to finite length with

Page 94
Digital Signal Processing 10EC52

( assuming M odd) choosing proper window

h(n) hd (n) w(n) where

w(n) is the window function defined for ( M 1) / 2 n ( M 1) / 2

3. Introduce h(n) = h(-n) for linear phase characteristics

4. Write the expression for H(z); this is non-causal realization

5. To obtain causal realization H‟(z) = z -(M-1)/2 H(z)

Exercise Problems
Prob 1: Design an ideal highpass filter with a frequency response:

H d (e j ) 1 for
4
0 | |
4

using a hanning window with M = 11 and plot the frequency response.

/4
1
hd (n) [ e j nd e j nd ]
2 /4

Page 95
Digital Signal Processing 10EC52

1 n
hd (n) [sin n sin ] for n and n 0
n 4
/4
1 3
hd (0) [ d d ] 0.75
2 /4
4

hd(1) = hd(-1)=-0.225
hd(2) = hd(-2)= -0.159
hd(3) = hd(-3)= -0.075
hd(4) = hd(-4)= 0
hd(5) = hd(-5) = 0.045
The hamming window function is given by

2 n M 1 M 1
whn (n) 0.5 0.5 cos ( ) n ( )
M 1 2 2
0 otherwise

for N 11
n
whn (n) 0.5 0.5 cos 5 n 5
5

whn(0) = 1
whn(1) = whn(-1)=0.9045
whn(2)= whn(-2)=0.655
whn(3)= whn(-3)= 0.345
whn(4)= whn(-4)=0.0945
whn(5)= whn(-5)=0

h(n)= whn(n)hd(n)

h(n)=[0 0 -0.026 -0.104 -0.204 0.75 -0.204 -0.104 -0.026 0 0]

h' ( n) h(n 5)
2 3 4 5 6 7 8
H ' ( z) 0.026z 0.104z 0.204z 0.75z 0.204z 0.104z 0.026z

Using the equation

Page 96
Digital Signal Processing 10EC52
M 3
2
M 1 M 1
H r (e jw ) [h( ) 2 h(n) cos ( n)
2 n 0 2
4
H r (e jw ) 0.75) 2 h(n) cos (5 n)
n 0

The magnitude response is given by,

|Hr(e jω)| = |0.75 - 0.408cosω - 0.208 cos2ω - 0.052cos3ω|

ω in degrees = [0 15 30 45 60 75 90 105 120 135 150 165 180]

|H(e jω)| in dBs = [-21.72 -17.14 -10.67 -6.05 -3.07 -1.297 -0.3726

-0.0087 0.052 0.015 0 0 0.017]

Page 97
Digital Signal Processing 10EC52

Prob 2 : Design a filter with a frequency response:

H d (e j ) e j3
for
4 4
0 | |
4

using a Hanning window with M = 7

Soln:
The freq resp is having a term e –jω(M-1)/2 which gives h(n) symmetrical about
n = M-1/2 = 3 i.e we get a causal sequence.

/4
1 j3
hd (n) e e j nd
2 /4

sin (n 3)
4
(n 3)
this gives hd (0) hd (6) 0.075
hd (1) hd (5) 0.159
hd (2) hd (4) 0.22
hd (3) 0.25

The Hanning window function values are given by

whn(0) = whn(6) =0

whn(1)= whn(5) =0.25

whn(2)= whn(4) =0.75

whn(3)=1

h(n)=hd(n) whn(n)

h(n)=[0 0.03975 0.165 0.25 0.165 0.3975 0]

Page 98
 Digital Signal Processing 10EC52

6.9 Design of Linear Phase FIR filters using Frequency Sampling method:

6.9.1 Motivation: We know that DFT of a finite duration DT sequence is obtained by sampling
FT of the sequence then DFT samples can be used in reconstructing original time domain
samples if frequency domain sampling was done correctly. The samples of FT of h(n) i.e., H(k)
are sufficient to recover h(n).

Since the designed filter has to be realizable then h(n) has to be real, hence even
symmetry properties for mag response |H(k)| and odd symmetry properties for phase response
can be applied. Also, symmetry for h(n) is applied to obtain linear phase chas.

Fro DFT relationship we have

N 1
1
h( n) H ( k )e j 2 kn / N
for n 0,1,...... N 1
N k 0
N 1
j 2 kn / N
H (k ) h ( n )e for k 0,1,......... N 1
n 0

Also we know H(k) = H(z)|z=e j2πkn/N

The system function H(z) is given by

N 1
n
H ( z) h( n ) z
n 0

Substituting for h(n) from IDFT relationship

N N 1
1 z H (k )
H ( z)
N k 0 1 e j2 kn / N
z 1

Page 99
Digital Signal Processing 10EC52

Since H(k) is obtained by sampling H(ejω) hence the method is called Frequency Sampling
Technique.

Since the impulse response samples or coefficients of the filter has to be real for filter to be
realizable with simple arithmetic operations, properties of DFT of real sequence can be used.
The following properties of DFT for real sequences are useful:

H*(k) = H(N-k)

|H(k)|=|H(N-k)| - magnitude response is even

θ(k) = - θ(N-k) – Phase response is odd

N 1
1
h( n) H ( k )e j 2 kn / N
can be rewritten as (for N odd)
N k 0

N 1
1
h( n) H ( 0) H ( k )e j 2 kn / N

N k 1
N 1/ 2 N 1
1
h( n) H ( 0) H ( k )e j 2 kn / N
H ( k )e j 2 kn / N

N k 1 k N 1/ 2

Using substitution k = N – r or r = N- k in the second substitution

with r going from now (N- 1)/2 to 1 as k goes from 1 to (N-1)/2

( N 1) / 2 ( N 1) / 2
1 j 2 kn / N j 2 kn / N
h( n) H (0) H ( k )e H (N k )e
N k 1 k 1
( N 1) / 2 ( N 1) / 2
1
h( n) H (0) H ( k )e j 2 kn / N
H * ( k )e j 2 kn / N

N k 1 k 1
( N 1) / 2 ( N 1) / 2
1
h( n) H (0) H ( k )e j 2 kn / N
( H ( k )e j 2 kn / N
)*
N k 1 k 1
( N 1) / 2
1
h( n) H (0) ( H ( k )e j 2 kn / N
( H ( k )e j 2 kn / N
)*
N k 1
( N 1) / 2
1
h( n) H (0) 2 Re( H (k )e j 2 kn / N

N k 1

Similarly for N even we have

( N 1) / 2
1
h( n) H (0) 2 Re( H (k )e j 2 kn / N

N k 1

Page 100
Digital Signal Processing 10EC52

Using the symmetry property h(n)= h (N-1-n) we can obtain Linear phase FIR filters using the
frequency sampling technique.

Exercise problems
Prob 1 : Design a LP FIR filter using Freq sampling technique having cutoff freq of π/2
rad/sample. The filter should have linear phase and length of 17.

The desired response can be expressed as

M 1
j ( )
H d (e j ) e 2
for | | c
0 otherwise
with M 17 and c /2

H d (e j ) e j 8
for 0 /2
0 for /2

2 k 2 k
Selecting k for k 0,1,...... 16
M 17

H (k ) H d (e j ) | 2 k
17

2 k
j
17
8 2 k
H (k ) e for 0
17 2
2 k
0 for /2
17
16 k
j 17
H (k ) e 17 for 0 k
4
17 17
0 for k
4 2

The range for “k” can be adjusted to be an integer such as

0 k 4
and 5 k 8

Page 101
Digital Signal Processing 10EC52

The freq response is given by

2 k
j 8
17
H (k ) e for 0 k 4
0 for 5 k 8

Using these value of H(k) we obtain h(n) from the equation

( M 1) / 2
1
h( n) ( H (0) 2 Re( H (k )e j 2 kn / M ))
M k 1
4
1
i.e., h(n) (1 2 Re(e j16 k / 17 e j 2 kn / 17 ))
17 k 1
4
1 2 k (8 n)
h( n) ( H (0) 2 cos( ) for n 0,1,........ 16
17 k 1 17

Even though k varies from 0 to 16 since we considered ω varying between 0 and π/2
only k values from 0 to 8 are considered
While finding h(n) we observe symmetry in h(n) such that n varying 0 to 7 and 9 to 16
have same set of h(n)

6.10 Design of FIR Differentiator

Differentiators are widely used in Digital and Analog systems whenever a derivative
of the signal is needed. Ideal differentiator has pure linear magnitude response in the freq
range –π to +π. The typical frequency response characteristics is as shown in the below figure.

Page 102
Digital Signal Processing 10EC52

Problem 2: Design an Ideal Differentiator using a) rectangular window and b)Hamming

window with length of the system = 7.

Solution:
As seen from differentiator frequency chars. It is defined as

H(ejω) = jω between –π to +π

1 cos n
hd (n) j ej n d n and n 0
2 n
The hd(n) is an add function with hd(n)=-hd(-n) and hd(0)=0

a) rectangular window

h(n)=hd(n)wr(n)

h(1)=-h(-1)=hd(1)=-1
h(2)=-h(-2)=hd(2)=0.5
h(3)=-h(-3)=hd(3)=-0.33

h‟(n)=h(n-3) for causal system

thus,

1 2 4 5 6
H ' ( z) 0.33 0.5 z z z 0.5 z 0.33z

Also from the equation

( M 3) / 2
M 1
H r (e j ) 2 h(n) sin ( n)
n 0 2

For M=7 and h‟(n) as found above we obtain this as

H r (e j ) 0.66 sin 3 sin 2 2 sin

H (e j ) jH r (e j ) j (0.66 sin 3 sin 2 2 sin )

b) Hamming window
h(n)=hd(n)wh(n)

where wh(n) is given by

Page 103
Digital Signal Processing 10EC52

2 n
wh ( n ) 0.54 0.46 cos ( M 1) / 2 n ( M 1) / 2
( M 1)
0 otherwise

For the present problem

n
wh (n) 0.54 0.46 cos 3 n 3
3
The window function coefficients are given by for n=-3 to +3
Wh(n)= [0.08 0.31 0.77 1 0.77 0.31 0.08]

Thus h‟(n) = h(n-5) = [0.0267, -0.155, 0.77, 0, -0.77, 0.155, -0.0267]

Similar to the earlier case of rectangular window we can write the freq response of
differentiator as
H (e j ) jH r (e j ) j (0.0534sin 3 0.31sin 2 1.54 sin )

We observe
With rectangular window, the effect of ripple is more and transition band width is
small compared with hamming window
With hamming window, effect of ripple is less whereas transition band is more

6.11 Design of FIR Hilbert transformer:

Hilbert transformers are used to obtain phase shift of 90 degree. They are also called j
operators. They are typically required in quadrature signal processing. The Hilbert transformer

Page 104
Digital Signal Processing 10EC52

is very useful when out of phase component (or imaginary part) need to be generated from
available real component of the signal.

Problem 3: Design an ideal Hilbert transformer using a) rectangular window and b)

Blackman Window with M = 11

Solution:
As seen from freq chars it is defined as

H d (e j ) j 0
j 0

The impulse response is given by

0
1 (1 cos n)
hd (n) [ je j n d je j n d ] n except n 0
2 0
n
At n = 0 it is hd(0) = 0 and hd(n) is an odd function

a) Rectangular window
h(n) = hd(n) wr(n) = hd(n) for -5 ≥n ≥5

h‟(n)=h(n-5)

h(n)= [-0.127, 0, -0.212, 0, -0.636, 0, 0.636, 0, 0.212, 0, 0.127]

Page 105
Digital Signal Processing 10EC52
4
H r (e j ) 2 h(n) sin (5 n)
n 0

H (e j ) j | H r (e j ) | j{0.254sin 5 0.424sin 3 1.272sin }

b) Blackman Window
window function is defined as
n 2 n
wb (n) 0.42 0.5 cos 0.08 cos 5 n 5
5 5
0 otherwise

Wb(n) = [0, 0.04, 0.2, 0.509,0.849,1,0.849, 0.509, 0.2, 0.04,0] for -5≥n≥5

h‟(n) = h(n-5) = [0, 0, -0.0424, 0, -0.5405, 0, 0.5405, 0, 0.0424, 0, 0]

H (e j ) j[0.0848sin 3 1.0810sin ]

Page 106
Digital Signal Processing 10EC52

Recommended questions with solution

Question1

Solution:-

(b) Magnitude plot

Page 107
Digital Signal Processing 10EC52

Phase plot

(c) Hamming window

(d) Bartlett window

Page 108
Digital Signal Processing 10EC52

Question 2

Solution:-

Page 109
Digital Signal Processing 10EC52

Page 110
Digital Signal Processing 10EC52

Question 3

Solution:-

Magnitude and phase response

Page 111
Digital Signal Processing 10EC52

Question 4

Solu

tion:-

Page 112
Digital Signal Processing 10EC52

UNIT 7

Design of IIR Filters from Analog Filters

CONTENTS:-

1. DESIGN OF IIR FILTERS FROM ANALOG FILTERS (BUTTERWORTH AND CHEBYSHEV)

2. IMPULSE INVARIANCE METHOD

3. MAPPING OF TRANSFER FUNCTIONS

4. APPROXIMATION OF DERIVATIVE (BACKWARD DIFFERENCE AND BILINEAR

TRANSFORMATION) METHOD, MATCHED Z TRANSFORMS

5. VERIFICATION FOR STABILITY AND LINEARITY DURING MAPPING

RECOMMENDED READINGS:-

1. DIGITAL SIGNAL PROCESSING – PRINCIPLES ALGORITHMS & APPLICATIONS, PROAKIS &

MONALAKIS, PEARSON EDUCATION, 4TH EDITION, NEW DELHI, 2007.

2. DISCRETE TIME SIGNAL PROCESSING, OPPENHEIM & SCHAFFER, PHI, 2003.

3. DIGITAL SIGNAL PROCESSING, S. K. MITRA, TATA MC-GRAW HILL, 2ND EDITION, 2004.

Page 113
Digital Signal Processing 10EC52

UNIT - 7
DESIGN OF IIR FILTERS FROM ANALOG FILTERS
(BUTTERWORTH AND CHEBYSHEV)

7.1 Introduction
A digital filter is a linear shift-invariant discrete-time system that is realized using finite
precision arithmetic. The design of digital filters involves three basic steps:

 The specification of the desired properties of the system.

 The approximation of these specifications using a causal discrete-time system.
 The realization of these specifications using _nite precision arithmetic.

These three steps are independent; here we focus our attention on the second step.
The desired digital filter is to be used to filter a digital signal that is derived from an analog
signal by means of periodic sampling. The speci_cations for both analog and digital filters are
often given in the frequency domain, as for example in the design of low
pass, high pass, band pass and band elimination filters. Given the sampling rate, it is straight
forward to convert from frequency specifications on an analog _lter to frequency
speci_cations on the corresponding digital filter, the analog frequencies being in terms of Hertz
and digital frequencies being in terms of radian frequency or angle around the unit circle with
the point Z=-1 corresponding to half the sampling frequency. The least confusing point of
view toward digital filter design is to consider the filter as being specified in terms of angle
around the unit circle rather than in terms of analog frequencies.

Figure 7.1: Tolerance limits for approximation of ideal low-pass filter

Page 114
Digital Signal Processing 10EC52

A separate problem is that of determining an appropriate set of specifications on the digital

filter. In the case of a low pass filter, for example, the specifications often take the
form of a tolerance scheme, as shown in Fig. 4.1

Many of the filters used in practice are specified by such a tolerance scheme, with no
constraints on the phase response other than those imposed by stability and causality
requirements; i.e., the poles of the system function must lie inside the unit circle. Given a set
of specifications in the form of Fig. 7.1, the next step is to and a discrete time linear system
whose frequency response falls within the prescribed tolerances. At this point the filter design
problem becomes a problem in approximation. In the case of infinite impulse response (IIR)
filters, we must approximate the desired frequency response by a rational function, while in the
finite impulse response (FIR) filters case we are concerned with polynomial approximation.

7.2 Design of IIR Filters from Analog Filters:

The traditional approach to the design of IIR digital filters involves the transformation of an
analog filter into a digital filter meeting prescribed specifications. This is a reasonable
approach because:

 The art of analog filter design is highly advanced and since useful results can be
achieved, it is advantageous to utilize the design procedures already developed for
analog filters.
 Many useful analog design methods have relatively simple closed-form design
formulas.

Therefore, digital filter design methods based on analog design formulas are rather simple to
implement.
An analog system can be described by the differential equation

------------------------------------------------------------7.1

And the corresponding rational function is

Page 115
Digital Signal Processing 10EC52

---------------------------------------------------------7.2

The corresponding description for digital filters has the form

--------------------------------------------------7.3
and the rational function

--------------------------------------------------------7.4

In transforming an analog filter to a digital filter we must therefore obtain either H(z)or h(n)
(inverse Z-transform of H(z) i.e., impulse response) from the analog filter design. In such
transformations, we want the imaginary axis of the S-plane to map into the finite circle of the
Z-plane, a stable analog filter should be transformed to a stable digital filter. That is, if the
analog filter has poles only in the left-half of S-plane, then the digital filter must have poles
only inside the unit circle. These constraints are basic to all the techniques discussed

7.3 IIR Filter Design by Impulse Invariance:

This technique of transforming an analog filter design to a digital filter design corresponds to
choosing the unit-sample response of the digital filter as equally spaced samples of the impulse
response of the analog filter. That is,

-------------------------------------------------------------------------7.5
Where T is the sampling period. Because of uniform sampling, we have

---------------------------------------------7.6

Or

---------------------------------------------7.7
Page 116
Digital Signal Processing 10EC52

Figure 7.2: Mapping of s-plane into z-plane

Where s = jω and Ω=ω/T, is the frequency in analog domain and ω is the frequency in digital
domain.
From the relationship Z = e ST it is seen that strips of width 2π/T in the S-plane map into the
entire Z-plane as shown in Fig. 7.2. The left half of each S-plane strip maps into interior of the
unit circle, the right half of each S-plane strip maps into the exterior of the unit circle, and the
imaginary axis of length 2π/T of S-plane maps on to once round the unit circle of Z-plane.
Each horizontal strip of the S-plane is overlaid onto the Z-plane to form the digital filter
function from analog filter function. The frequency response of the digital filter is related to
the frequency response of the

Figure 7.3: Illustration of the effects of aliasing in the impulse invariance technique

Page 117
Digital Signal Processing 10EC52

analog filter as

------------------------------------------------7.8

From the discussion of the sampling theorem it is clear that if and only if

Then

Unfortunately, any practical analog filter will not be band limited, and consequently there is
interference between successive terms in Eq. (7.8) as illustrated in Fig. 7.3. Because of the
aliasing that occurs in the sampling process, the frequency response of the resulting digital
filter will not be identical to the original analog frequency response. To get the filter design
procedure, let us consider the system function of the analog filter expressed in terms of a
partial-fraction expansion

-----------------------------------------------------------------------7.9

The corresponding impulse response is

--------------------------------------------------------------- 7.10

And the unit-sample response of the digital filter is then

--------------7.11

The system function of the digital filter H(z) is given by

------------------------------------------------------------7.12

Page 118
Digital Signal Processing 10EC52

In comparing Eqs. (7.9) and (7.12) we observe that a pole at s=sk in the S-plane transforms to
a pole at expskT in the Z-plane. It is important to recognize that the impulse invariant design
procedure does not correspond to a mapping of the S-plane to the Z-plane.

7.4 IIR Filter Design By Approximation Of Derivatives:

A second approach to design of a digital filter is to approximate the derivatives in Eq. (4.1) by
finite differences. If the samples are closer together, the approximation to the derivative would
be increasingly accurate. For example, suppose that the first derivative is approximated by the
first backward difference

--------------------------7.13

Where y(n)=y(nT). Approximation to higher-order derivatives are obtained by repeated

application of Eq. (7.13); i.e.,

-------------------------- 7.14
For convenience we define

-------------------------------------------------------------------7.15

Applying Eqs. (7.13), (7.14) and (7.15) to (7.1), we obtain

---------------------------------------------7.16

Where y(n) = ya(nT) and x(n) = xa(nT). We note that the operation ∆(1)[ ] is a linear shift-
invariant operator and that ∆(k)[ ] can be viewed as a cascade of (k) operators ∆(1)[ ]. In
particular

And

Page 119
Digital Signal Processing 10EC52

Thus taking the Z-transform of each side in Eq. (7.16), we obtain

------------------------------------------------------------7.17

Comparing Eq. (7.17) to (7.2), we observe that the digital transfer function can be obtained
directly from the analog transfer function by means of a substitution of variables

---------------------------------------------------------------------------------7.18

So that, this technique does indeed truly correspond to a mapping of the S-plane to the Z-
plane, according to Eq. (7.18). To investigate the properties of this mapping, we must express
z as a function of s, obtaining

Substituting s = jΩ, i.e., imaginary axis in S-plane

------------------------------------------------------7.19

Which corresponds to a circle whose center is at z =1/2 and radius is 1/2, as shown in Fig. 7.4.
It is easily verified that the left half of the S-plane maps into the inside of the small circle and
the right half of the S-plane maps onto the outside of the small circle. Therefore, although the

Page 120
Digital Signal Processing 10EC52

requirement of mapping the jΩ-axis to the unit circle is not satisfied, this mapping does satisfy
the stability condition.

Figure 4.4: Mapping of s-plane to z-plane corresponding to first backward-difference

approximation to the derivative

In contrast to the impulse invariance technique, decreasing the sampling period T, theoretically
produces a better filter since the spectrum tends to be concentrated in a very small region of
the unit circle. These two procedures are highly unsatisfactory for anything but low pass
filters. An alternative approximation to the derivative is a forward difference and it provides a
mapping into the unstable digital filters.

7.5 IIR Filter Design By The Bilinear Transformation:

In the previous section a digital filter was derived by approximating derivatives by differences.
An alternative procedure is based on integrating the differential equation and then using a
numerical approximation to the integral. Consider the first - order equation

-----------------------------------------------------------7.20

Where y‟a(t) is the first derivative of ya(t). The corresponding analog system function is

We can write ya(t) as an integral of y‟a(t), as in

Page 121
Digital Signal Processing 10EC52

In particular, if t = nT and t 0 = (n - 1)T,

If this integral is approximated by a trapezoidal rule, we can write

----------------------7.21

However, from Eq. (7.20),

Substituting into Eq. (4.21) we obtain

Where y(n) = y(nT) and x(n) = x(nT). Taking the Z-transform and solving for H(z) gives

--------------------------------------------7.22

From Eq. (7.22) it is clear that H(z) is obtained from Ha(s) by the substitution

-------------------------------------------------------------------7.23

That is,

--------------------------------------------------------------7.24

This can be shown to hold in general since an Nth - order differential equation of the form of
Eq. (7.1) can be written as a set of N first-order equations of the form of Eq. (7.20). Solving
Eq. (7.23) for z gives

Page 122
Digital Signal Processing 10EC52

----------------------------------------------------------------------------7.25
The invertible transformation of Eq. (7.23) is recognized as a bilinear transformation. To see
that this mapping has the property that the imaginary axis in the s-plane maps onto the unit
circle in the z-plane, consider z = ejω, then from Eq. (7.23), s is given by

Figure 7.5: Mapping of analog frequency axis onto the unit circle using the bilinear
Transformation

Thus for z on the unit circle, ζ = 0 and Ω and ω are related by

T Ω/2 = tan (ω/2)

or
ω = 2 tan -1(T Ω/2)

Page 123
Digital Signal Processing 10EC52

This relationship is plotted in Fig. (7.5), and it is referred as frequency warping. From the
_gure it is clear that the positive and negative imaginary axis of the s-plane are mapped,
respectively, into the upper and lower halves of the unit circle in the z-plane. In addition to the
fact that the imaginary axis in the s-plane maps into the unit circle in the z-plane, the left half
of the s-plane maps to the inside of the unit circle and the right half of the s-plane maps to the
outside of the unit circle, as shown in Fig. (7.6). Thus we see that the use of the bilinear
transformation yields stable digital filter from analog filter. Also this transformation avoids the
problem of aliasing encountered with the use of impulse invariance, because it maps the entire
imaginary axis in the s-plane onto the unit circle in the z-plane. The price paid for this,
however, is the introduction of a distortion in the frequency axis.

Figure 4.6: Mapping of the s-plane into the z-plane using the bilinear transformation

7.6 The Matched-Z Transform:

Another method for converting an analog filter into an equivalent digital filter is to map
the poles and zeros of Ha(s) directly into poles and zeros in the z-plane. For analog filter

-----------------------------------------------------------------7.26

the corresponding digital filter is

---------------------------------------------------------7.27

Page 124
Digital Signal Processing 10EC52

Where T is the sampling interval. Thus each factor of the form (s-a) in Ha(s) is mapped
into the factor (1- eaT z-1).
Recommended questions with solution
Question 1

Page 125
Digital Signal Processing 10EC52

Question 2

Question 3

Page 126
Digital Signal Processing 10EC52

Question 4

Question 5

Page 127
Digital Signal Processing 10EC52

Page 128
Digital Signal Processing 10EC52

Question 6

Page 129
Digital Signal Processing 10EC52

Question 7

Page 130
Digital Signal Processing 10EC52

Page 131
Digital Signal Processing 10EC52

Page 132
Digital Signal Processing 10EC52

UNIT 8
Implementation of Discrete time systems

CONTENTS:-

1. IMPLEMENTATION OF DISCRETE-TIME SYSTEMS

2. STRUCTURES FOR IIR AND FIR SYSTEMS

3. DIRECT FORM I AND DIRECT FORM II SYSTEMS,

4. CASCADE, LATTICE AND PARALLEL REALIZATION.

RECOMMENDED READINGS:-

1. DIGITAL SIGNAL PROCESSING – PRINCIPLES ALGORITHMS & APPLICATIONS, PROAKIS &

MONALAKIS, PEARSON EDUCATION, 4TH EDITION, NEW DELHI, 2007.

2. DISCRETE TIME SIGNAL PROCESSING, OPPENHEIM & SCHAFFER, PHI, 2003.

3. DIGITAL SIGNAL PROCESSING, S. K. MITRA, TATA MC-GRAW HILL, 2ND EDITION, 2004.

Page 133
Digital Signal Processing 10EC52

UNIT 8
Implementation of Discrete-Time Systems

8.1 Introduction
The two important forms of expressing system leading to different realizations of FIR & IIR
filters are
a) Difference equation form
N M
y ( n) a k y (n k ) bk x(n k )
k 1 k 1
b) Ration of polynomials
M
k
bk Z
k 0
H (Z ) N
k
1 ak Z
k 1

The following factors influence choice of a specific realization,

Computational complexity
Memory requirements
Finite-word-length
Pipeline / parallel processing

8.1.1 Computation Complexity

This is do with number of arithmetic operations i.e. multiplication, addition & divisions. If
the realization can have less of these then it will be less complex computationally.
In the recent processors the fetch time from memory & number of times a comparison between
two numbers is performed per output sample is also considered and found to be important
from the point of view of computational complexity.

8.1.2 Memory requirements

This is basically number of memory locations required to store the system parameters,
past inputs, past outputs, and any intermediate computed values. Any realization requiring less
of these is preferred.

8.1.3 Finite-word-length effects

These effects refer to the quantization effects that are inherent in any digital
implementation of the system, either in hardware or in software. No computing system has
infinite precision. With finite precision there is bound to be errors. These effects are basically
to do with truncation & rounding-off of samples. The extent of this effect varies with type of
arithmetic used(fixed or floating). The serious issue is that the effects have influence on
system characteristics. A structure which is less sensitive to this effect need to be chosen.

Page 134
Digital Signal Processing 10EC52

8.1.4 Pipeline / Parallel Processing

This is to do with suitability of the structure for pipelining & parallel processing. The
parallel processing can be in software or hardware. Longer pipelining make the system more
efficient.

8.2 Structure for FIR Systems:

FIR system is described by,
M 1
y ( n) bk x(n k )
k 0
Or equivalently, the system function
M 1
k
H (Z ) bk Z
k 0

bn 0 n n 1
Where we can identify h(n)
0 otherwise
Different FIR Structures used in practice are,
1. Direct form
2. Cascade form
3. Frequency-sampling realization
4. Lattice realization

8.2.1 Direct – Form Structure

Convolution formula is used to express FIR system given by,

M 1
y ( n) h( k ) x ( n k )
k 0

It is Non recursive in structure

As can be seen from the above implementation it requires M-1 memory locations for
storing the M-1 previous inputs
It requires computationally M multiplications and M-1 additions per output point
It is more popularly referred to as tapped delay line or transversal system
Efficient structure with linear phase characteristics are possible where
h(n) h(M 1 n)

Page 135
Digital Signal Processing 10EC52

Prob:
Realize the following system function using minimum number of multiplication
1 1 1 2 1 3 1 4
(1) H ( Z ) 1 Z Z Z Z Z 5
3 4 4 3
1 1 1 1
We recognize h(n) 1, , , , , 1
3 4 4 3
M is even = 6, and we observe h(n) = h(M-1-n) h(n) = h(5-n)
i.e h(0) = h(5) h(1) = h(4) h(2) = h(3)
Direct form structure for Linear phase FIR can be realized

Exercise: Realize the following using system function using minimum number of
multiplication.
1 1 1 2 1 3 1 5 1 6 1 7
H (Z ) 1 Z Z Z Z Z Z Z 8
4 3 2 2 3 4
1 1 1 1 1 1
m=9 h(n) 1, , , , , , , 1
4 3 2 2 3 4
odd symmetry
h(n) = -h(M-1-n); h(n) = -h(8-n); h(m-1/2) = h(4) = 0
h(0) = -h(8); h(1) = -h(7); h(2) = -h(6); h(3) = -h(5)

Page 136
Digital Signal Processing 10EC52

8.2.2 Cascade – Form Structure

The system function H(Z) is factored into product of second – order FIR system
K
H (Z ) H k (Z )
k 1

Where H k (Z ) bk 0 bk1 Z 1 bk 2 Z 2 k = 1, 2, ….. K

and K = integer part of (M+1) / 2
The filter parameter b0 may be equally distributed among the K filter section, such that b0
= b10 b20 …. bk0 or it may be assigned to a single filter section. The zeros of H(z) are grouped
in pairs to produce the second – order FIR system. Pairs of complex-conjugate roots are
formed so that the coefficients {bki} are real valued.

Page 137
Digital Signal Processing 10EC52

In case of linear –phase FIR filter, the symmetry in h(n) implies that the zeros of H(z)
also exhibit a form of symmetry. If zk and zk* are pair of complex – conjugate zeros then
1/zk and 1/zk* are also a pair complex –conjugate zeros. Thus simplified fourth order
sections are formed. This is shown below,

H k ( z) C k 0 (1 z k z 1 )(1 z k * z 1 )(1 z 1 / z k )(1 z 1 / z k * )

1 2 3 4
Ck 0 C k1 z Ck 2 z C k1 z z

Problem: Realize the difference equation

y(n) x(n) 0.25x(n 1) 0.5x(n 2) 0.75x(n 3) x(n 4)
in cascade form.

1 2 3
Y ( z) X ( z ){1 0.25z 0.5 z 0.75z z 4)
1 2
H ( z ) 1 0.25z 0.5 z 0.75z _ 3 z 4
Soln: 1
H ( z) (1 1.1219z 1.2181z 2 )(1 1.3719z 1
0.821z 2 )
H ( z) H1 ( z)H 2 ( z)

8.3 Frequency sampling realization:

We can express system function H(z) in terms of DFT samples H(k) which is given by
1 N 1 H (k )
H ( z ) (1 z N )
N k 0 1 WN k z 1

Page 138
Digital Signal Processing 10EC52

This form can be realized with cascade of FIR and IIR structures. The term (1-z-N) is realized
1 N 1 H (k )
as FIR and the term as IIR structure.
N k 0 1 WN k z 1

The realization of the above freq sampling form shows necessity of complex arithmetic.
Incorporating symmetry in h(n) and symmetry properties of DFT of real sequences the
realization can be modified to have only real coefficients.

8.4 Lattice structures

Lattice structures offer many interesting features:

1. Upgrading filter orders is simple. Only additional stages need to be added instead of
redesigning the whole filter and recalculating the filter coefficients.
2. These filters are computationally very efficient than other filter structures in a filter
bank applications (eg. Wavelet Transform)
3. Lattice filters are less sensitive to finite word length effects.

Consider

m
Y ( z)
H ( z) 1 a m (i) z i
X ( z) i 1
m is the order of the FIR filter and am(0)=1

when m = 1 Y(z)/ X(z) = 1+ a1(1) z-1

Page 139
Digital Signal Processing 10EC52

y(n)= x(n)+ a1(1)x(n-1)

f1(n) is known as upper channel output and r1(n)as lower channel output.

f0(n)= r0(n)=x(n)

The outputs are

f 1 ( n) f 0 (n) k1 r0 (n 1) 1a
r1 (n) k1 f 0 (n) r0 (n 1) 1b
if k1 a1 (1), then f 1 ( n) y ( n)

If m=2

Y ( z)
1 a 2 (1) z 1 a 2 (2) z 2
X ( z)
y (n) x(n) a 2 (1) x(n 1) a 2 (2) x(n 2)
y ( n) f1 (n) k 2 r1 (n 1) (2)

Substituting 1a and 1b in (2)

y ( n) f 0 ( n) k1r0 (n 1) k 2 [k1 f 0 (n 1) r0 (n 2)]

f 0 ( n) k1r0 (n 1) k 2 k1 f 0 (n 1) k 2 r0 (n 2)]
sin ce f 0 ( n) r0 (n) x ( n)
y ( n) x ( n) k1 x(n 1) k 2 k1 x(n 1) k 2 x(n 2)]
x ( n) (k1 k1k 2 ) x(n 1) k 2 x(n 2)

We recognize

Page 140
Digital Signal Processing 10EC52

a 2 (1) k1 k1 k 2
a 2 (1) k2

Solving the above equation we get

a 2 (1)
k1 and k 2 a 2 (2) (4)
1 a 2 (2)

Equation (3) means that, the lattice structure for a second-order filter is simply a cascade of
two first-order filters with k1 and k2 as defined in eq (4)

Similar to above, an Mth order FIR filter can be implemented by lattice structures with
M – stages

8.4.1 Direct Form –I to lattice structure

For m = M, M-1, ………..2, 1 do

km a m ( m)
a m (i ) a m (m)a m (m i )
a m 1 (i ) 1 i m 1
1 k m2

Page 141
Digital Signal Processing 10EC52

The above expression fails if k m=1. This is an indication that there isa zero on the unit
circle. If km=1, factor out this root from A(z) and the recursive formula can be applied
for reduced order system.

for m 2 and m 1
k2 a 2 ( 2) & k 1 a1 (1)

for m 2&i 1
a 2 (1) a 2 (2)a 2 (1) a 2 (1)[1 a 2 (2)] a 2 (1)
a1 (1)
1 k 22 1 a 22 (2) 1 a 2 (2)
a 2 (1)
Thus k1
1 a 2 ( 2)

8.4.2 Lattice to direct form –I

For m = 1,2,…….M-1

a m (0) 1
a m ( m) km
a m (i ) a m 1 (i ) a m (m)a m 1 (m i ) 1 i m 1

Problem:
Given FIR filter H (Z ) 1 2Z 1 1
3 Z 2
obtain lattice structure for the same
Given a1 (1) 2 , a2 (2) 13
Using the recursive equation for
m = M, M-1, ……, 2, 1
here M=2 therefore m = 2, 1
if m=2 k 2 a2 (2) 13
if m=1 k1 a1 (1)
also, when m=2 and i=1
a 2 (1) 2 3
a1 (1)
1 a 2 (2) 1 3 2 1

Hence k1 a1 (1) 3 2

Page 142
Digital Signal Processing 10EC52

Recommended questions with solution

Problem:1
1 1 1
Consider an FIR lattice filter with co-efficients k1 , k2 , k3 . Determine the FIR
2 3 4
filter co-efficient for the direct form structure
( H (Z ) a3 (0) a3 (1)Z 1 a3 (2)Z 2 a3 (3)Z 3 )
1 1
a 3 ( 0) 1 a3 (3) k3 a2 (2) k2
4 3
1
a1 (1) k1
2

for m=2, i=1

a2 (1) a1 (1) a2 (2)a1 (1)

1 1
= a1 (1)[1 a 2 (2)] 1
2 3

4 2
=
6 3

for m=3, i=1

a3 (1) a 2 (1) a3 (3)a 2 (2)

2 1 1
= .
3 4 3
2 1 8 1
= =
3 12 12

9 3
=
12 4

for m=3 & i=2

a3 (2) a2 (2) a3 (3)a2 (1)

1 1 2
= .
3 4 3
1 1 2 1
=
3 6 6

Page 143
Digital Signal Processing 10EC52

3 1
=
6 2

3 1 1
a 3 ( 0) 1 , a3 (1) , a3 (2) , a3 (3)
4 2 4

8.5 Structures for IIR Filters

The IIR filters are represented by system function;
M
k
bk z
k 0
H(Z) = N
k
1 ak z
k 1
and corresponding difference equation given by,
N N
y ( n) a k y (n k ) bk x(n k )
k 1 k 0

Different realizations for IIR filters are,

1. Direct form-I
2. Direct form-II
3. Cascade form
4. Parallel form
5. Lattice form

8.5.1 Direct form-I

This is a straight forward implementation of difference equation which is very simple.
Typical Direct form – I realization is shown below . The upper branch is forward path and
lower branch is feedback path. The number of delays depends on presence of most previous
input and output samples in the difference equation.

Page 144
Digital Signal Processing 10EC52

8.5.2 Direct form-II

The given transfer function H(z) can be expressed as,

Y ( z) V ( z) Y ( z)
H ( z) .
X ( z) X ( z) V ( z)
where V(z) is an intermediate term. We identify,

V ( z) 1
N
-------------------all poles
X ( z) k
1 ak z
k 1

M
Y ( z)
1 bk z k -------------------all zeros
V ( z) k 1

The corresponding difference equations are,

N
v ( n) x ( n) a k v(n k )
k 1
M
y ( n) v ( n) bk v(n 1)
k 1

Page 145
Digital Signal Processing 10EC52

This realization requires M+N+! multiplications, M+N addition and the maximum of
{M, N} memory location

8.5.3 Cascade Form

The transfer function of a system can be expressed as,

Page 146
Digital Signal Processing 10EC52

H ( z) H 1 ( z ) H 2 ( z )....H k ( z )

Where H k (Z ) could be first order or second order section realized in Direct form – II form
i.e.,
bk 0 bk1 Z 1 bk 2 Z 2
H k (Z )
1 ak1 Z 1 ak 2 Z 2
where K is the integer part of (N+1)/2
Similar to FIR cascade realization, the parameter b0 can be distributed equally among the
k filter section B0 that b0 = b10b20…..bk0. The second order sections are required to realize
section which has complex-conjugate poles with real co-efficients. Pairing the two complex-
conjugate poles with a pair of complex-conjugate zeros or real-valued zeros to form a
subsystem of the type shown above is done arbitrarily. There is no specific rule used in the
combination. Although all cascade realizations are equivalent for infinite precision arithmetic,
the various realizations may differ significantly when implemented with finite precision
arithmetic.

8.5.4 Parallel form structure

In the expression of transfer function, if N M we can express system function

N N
Ak
H (Z ) C C H k (Z )
k 11 pk Z 1 k 1

Where {pk} are the poles, {Ak} are the coefficients in the partial fraction expansion, and the
constant C is defined as C bN a N , The system realization of above form is shown below.

bk 0 bk1 Z 1
Where H k ( Z )
1 a k1 Z 1 a k 2 Z 2

Page 147
Digital Signal Processing 10EC52

Once again choice of using first- order or second-order sections depends on poles of the
denominator polynomial. If there are complex set of poles which are conjugative in nature then
a second order section is a must to have real coefficients.

Problem 2
Determine the
(i)Direct form-I (ii) Direct form-II (iii) Cascade &
(iv)Parallel form realization of the system function

1 1 2 1 1
10 1 2 Z 1 3 Z 1 2Z
H (Z ) 3 1 1 1 1 1 1
1 4 Z 1 8 Z 1 2 j 12 Z 1 1
2 j 12 Z

7 1 1 2 1
10 1 6 Z 3 Z 1 2Z
7 1 3 2 1 1 2
1 8 Z 32 Z 1 Z 2 Z

5 1 2 2 3
10 1 6 Z 2Z 3 Z
H (Z ) 15 1 47 2 17 3 3 4
1 8 Z 32 Z 32 Z 64 Z

( 14.75 12.90z 1 ) (24.50 26.82 z 1 )

H ( z)
7 1 3 2 1 2
(1 z z ) (1 z 1 z )
8 32 2

Page 148
Digital Signal Processing 10EC52

Cascade Form
H(z) = H1(z) H2(z)
Where
7 1 1 2
1 z z
H1 ( z) 6 3
7 1 3 2
1 z z
8 32

10(1 2 z 1 )
H1 ( z)
1 2
1 z 1 z
2

Page 149
Digital Signal Processing 10EC52

Parallel Form
H(z) = H1(z) + H2(z)

( 14.75 12.90z 1 ) (24.50 26.82 z 1 )

H ( z)
7 1 3 2 1 2
(1 z z ) (1 z 1 z )
8 32 2

Problem: 3
Obtain the direct form – I, direct form-II
Cascade and parallel form realization for the following system,
y(n)= -0.1 y(n-1)+0.2y(n-2)+3x(n)+3.6 x(n-1)+0.6 x(n-2)

Solution:
The Direct form realization is done directly from the given i/p – o/p equation, show in below
diagram

Direct form –II realization

Taking ZT on both sides and finding H(z)

1 2
Y ( z) 3 3.6 z 0.6 z
H ( z) 1 2
X ( z) 1 0.1z 0. 2 z

Page 150
Digital Signal Processing 10EC52

Cascade form realization

The transformer function can be expressed as:
(3 0.6 z 1 )(1 z 1 )
H ( z)
(1 0.5 z 1 )(1 0.4 z 1 )

which can be re written as

1
3 0.6 z 1 z 1
where H 1 ( z ) 1
and H 2 ( z )
1 0.5 z 1 0.4 z 1

Parallel Form realization

The transfer function can be expressed as

H(z) = C + H1(z) + H2(z) where H1(z) & H2(z) is given by,

7 1
H ( z) 3 1 1
1 0.4 z 1 0.5 z

Page 151
Digital Signal Processing 10EC52

8.6 Lattice Structure for IIR System:

Consider an All-pole system with system function.

1 1
H (Z ) N
k AN ( Z )
1 a N (k ) Z
k 1

The corresponding difference equation for this IIR system is,

N
y ( n) a N (k ) y (n k ) x ( n)
k 1
OR
N
x ( n) y ( n) a N (k ) y (n k )
k 1

For N=1
x(n) y(n) a1 (1) y(n 1)
Which can realized as,

We observe
x ( n) f 1 ( n)
y ( n) f 0 ( n) f1 (n) k1 g 0 (n 1)
x(n) k1 y(n 1) k1 a1 (1)
g 1 ( n) k1 f 0 (n) g 0 (n 1) k1 y (n) y (n 1)
For N=2, then
y(n) x(n) a2 (1) y(n 1) a2 (2) y(n 2)

Page 152
Digital Signal Processing 10EC52

This output can be obtained from a two-stage lattice filter as shown in below fig

f 2 (n) x(n)
f1 (n) f 2 (n) k 2 g1 (n 1)
g 2 (n) k 2 f1 (n) g1 (n 1)
f 0 ( n ) f 1 ( n ) k1 g 0 ( n 1)
g 1 ( n ) k1 f 0 ( n ) g 0 ( n 1)

y ( n) f 0 ( n) g 0 ( n) f1 (n) k1 g 0 (n 1)
f 2 (n) k 2 g1 (n 1) k1 g 0 (n 1)
f 2 (n) k 2 k1 f 0 (n 1) g 0 (n 2) k1 g 0 (n 1)
x(n) k 2 k1 y(n 1) y(n 2) k1 y(n 1)
x(n) k1 (1 k 2 ) y(n 1) k 2 y(n 2)
Similarly
g 2 (n) k 2 y(n) k1 (1 k 2 ) y(n 1) y(n 2)
We observe
a2 (0) 1; a2 (1) k1 (1 k 2 ); a2 (2) k2
N-stage IIR filter realized in lattice structure is,

f N ( n) x ( n)
f m 1 ( n) f m (n) k m g m 1 (n 1) m=N, N-1,---1
g m ( n) k m f m 1 ( n) g m 1 (n 1) m=N, N-1,---1

Page 153
Digital Signal Processing 10EC52

y ( n) f 0 ( n) g 0 ( n)

8.6.1 Conversion from lattice structure to direct form:

a m ( m) km ; a m (0) 1
a m (k ) a m 1 ( k ) a m ( m) a m 1 ( m k )

Conversion from direct form to lattice structure

a m 1 ( 0) 1 km a m (m)

a m ( k ) a m ( m) a m ( m k )
a m 1 (k )
1 a m2 (m)

8.6.2 Lattice – Ladder Structure:

A general IIR filter containing both poles and zeros can be realized using an all pole
lattice as the basic building block.

If,
M
k
bM ( k ) Z
BM ( Z ) k 0
H (Z ) N
AN ( Z ) k
1 a N (k ) Z
k 1
Where N M
A lattice structure can be constructed by first realizing an all-pole lattice co-efficients
k m , 1 m N for the denominator AN(Z), and then adding a ladder part for M=N. The
output of the ladder part can be expressed as a weighted linear combination of {g m(n)}.
Now the output is given by
M
y ( n) C m g m ( n)
m 0
Where {Cm} are called the ladder co-efficient and can be obtained using the recursive relation,
M
Cm bm Ci ai (i m); m=M, M-1, ….0
i m 1

Page 154
Digital Signal Processing 10EC52

Problem:4
Convert the following pole-zero IIR filter into a lattice ladder structure,
1 2Z 1 2Z 2 Z 3
H (Z ) 1 2 3
1 1324 Z
5
8 Z
1
3Z
Solution:
Given bM ( Z ) 1 2 Z 1 2 Z 2 Z 3
And AN (Z ) 1 13
24Z 1 85 Z 2
Z 3
1
3
13 5 1
a3 (0) 1; a3 (1) 24; a3 (2) 8; a3 (3) 3

k3 a3 (3) 13
Using the equation
a m ( k ) a m ( m) a m ( m k )
a m 1 (k )
1 a 2 m( m)
for m=3, k=1
13 1 5
a3 (1) a3 (3)a3 (2) 24 .
3 8 3
a 2 (1) 8
1 a32 (3) 1 1 2
3

for m=3, & k=2

a3 (2) a3 (3)a3 (1)
a 2 (2) k2
1 a32 (3)
5 1 13 45 13
8 .
3 24 72 1
1 8 2
1 9 9
for m=2, & k=1
a 2 (1) a 2 (2)a 2 (1)
a1 (1) k1
1 a 22 (2)

Page 155
Digital Signal Processing 10EC52
3 1 3 3 3
8 2 8 . 8 16 1
1 2 1 4
1 ( ) 2
1 4

1 1 1
for lattice structure k1 4 , k2 2 , k3 3
For ladder structure
M
Cm bm C1 .a1 (1 m) m=M, M-1,1,0
i m 1

C3 b3 1; C 2 b2 C3 a3 (1)
M=3 13
2 1.( ) 1.4583
24
3
C1 b1 c1 a1 (i m) m=1
i 2

b1 c2 a2 (1) c3 a3( 2)
2 1.4583 ( 83 ) 5
8 0.8281
3
c0 b0 c1a1 (i m)
i 1

b0 c1a1 (1) c2 a 2 (2) c3 a3 (3)

1 08281( 14 ) 1.4583( 12 ) 1
3 02695

To convert a lattice- ladder form into a direct form, we find an equation to obtain
a N (k ) from k m (m=1,2,………N) then equation for c m is recursively used to compute bm
(m=0,1,2,………M).

Problem 5

Page 156
Digital Signal Processing 10EC52

Page 157
Digital Signal Processing 10EC52

Question 6

Consider a FIR filter with system function:

H(z) = 1+2.82 Z-1 +3.4048z-2 +1.74z- 3. Sketch the direct form and lattice
realizations of the filter.

Page 158
Digital Signal Processing 10EC52

Page 159
Digital Signal Processing 10EC52

Page 160