Ece V Digital Signal Processing (10ec52) Notes

Digital Signal Processing
10EC52
SUBJECT CODE
: 10EC52
IA M ARKS EXAM HOURS EXAM M ARKS
: 25 : 03 : 100
NO. OF LECTURE HRS/WEEK : 04 TOTAL NO. OF LECTURE HRS. : 52
PART - A UNIT - 1 DISCRETE FREQUENCY DOMAIN SAMPLING AND RECONSTRUCTION OF DISCRETE TIME SIGNALS. DFT AS A LINEAR TRANSFORMATION, ITS RELATIONSHIP WITH OTHER TRANSFORMS. 7 HOURS UNIT - 2 PROPERTIES
OF
FOURIER
TRANSFORMS
(DFT):
EN
ALGORITHMS:
DFT, MULTIPLICATION OF TWO DFTS- THE CIRCULAR CONVOLUTION, ADDITIONAL DFT PROPERTIES, USE OF DFT IN LINEAR FILTERING, OVERLAP -SAVE AND OVERLAP -ADD METHOD. 6 HOURS
COMPUTATION OF
UNIT - 3
FAST-FOURIER-TRANSFORM (FFT)
TS .IN
DIRECT DFT
AND TO ANALOG
DFT,
NEED
FOR EFFICIENT COMPUTATION OF THE DFT (FFT ALGORITHMS).
TU D
8 HOURS
UNIT - 4 RADIX-2 FFT

TRANSFORM
ALGORITHM FOR THE COMPUTATION OF
TIME AND DECIMATION-IN-FREQUENCY
IDFTDECIMATION-INALGORITHMS. GOERTZEL ALGORITHM, AND CHIRP -Z 6 HOURS
TS
PART - B
USED ANALOG FILTERS
UNIT - 5
CI
IIR FILTER DESIGN: CHARACTERISTICS OF COMMONLY B UTTERWORTH AND CHEBYSHEVE FILTERS, ANALOG TRANSFORMATIONS. UNIT - 6
FREQUENCY
6 HOURS FIR FILTER DESIGN: INTRODUCTION TO FIR FILTERS, DESIGN OF FIR FILTERS USING RECTANGULAR, HAMMING, B ARTLET AND KAISER WINDOWS, FIR FILTER DESIGN USING
FREQUENCY SAMPLING TECHNIQUE
6 HOURS
CITSTUDENTS.IN
Page 1

UNIT - 7 DESIGN
OF
10EC52
IIR
FILTERS FROM ANALOG FILTERS
(BUTTERWORTH
IMPULSE INVARIANCE METHOD. DERIVATIVE
MAPPING
OF TRANSFER FUNCTIONS:
(BACKWARD DIFFERENCE AND BILINEAR M ATCHED Z TRANSFORMS, VERIFICATION FOR STABILITY AND LINEARITY DURING MAPPING 7 HOURS UNIT - 8
CHEBYSHEV) APPROXIMATION OF TRANSFORMATION) METHOD,

AND
TEXT BOOK:
REFERENCE BOOKS:
1. DISCRETE TIME SIGNAL PROCESSING, OPPENHEIM & SCHAFFER, PHI, 2003.
3. DIGITAL SIGNAL PROCESSING, LEE TAN: ELSIVIER PUBLICATIONS, 2007
CI
TS
TU D
2. DIGITAL SIGNAL PROCESSING, S. K. M ITRA, TATA M C-GRAW HILL, 2ND EDITION, 2004.
EN
DIGITAL SIGNAL PROCESSING PRINCIPLES ALGORITHMS & APPLICATIONS, PROAKIS & M ONALAKIS, PEARSON EDUCATION, 4TH EDITION, NEW DELHI, 2007.
CITSTUDENTS.IN
TS .IN
IMPLEMENTATION OF DISCRETE-TIME SYSTEMS: STRUCTURES FOR IIR AND FIR SYSTEMSDIRECT FORM I AND DIRECT FORM II SYSTEMS, CASCADE, LATTICE AND PARALLEL REALIZATION. 6 HOURS
Page 2
10EC52
INDEX SHEET
SL.NO I 1.1 1.2 1.3 1.4 II 2.1 2.2 2.3 2.4 2.5 III 3.1 3.2 3.2 IV 4.1 4.2 4.3 4.4 V 5.1 5.2 5.3 5.4 VI 6.1 6.2 TOPIC Unit-1: Discrete Fourier Transforms Frequency Domain Sampling Reconstruction of Discrete time signals DFT as a Linear Transform DFT relationship with other transforms UNIT - 2 : Properties of DFT Multiplication of two DFTs Circular convolution Additional DFT properties Use of DFT in Linear Filtering Overlap save and Overlap Add Method Solution to Problems UNIT 3 : Fast Fourier Transform Algorithms Direct computation of DFT Need for Efficient computation of DFT FFT algorithms UNIT 4 : Radix-2 FFT Algorithms for DFT and IDFT Decimation In Time Algorithms Decimation-in-Frequency Algorithms Goertzel Algorithms Chirp-z-Transforms UNIT 5 : IIR Filter Design Characteristics of commonly used Analog Filters Butterworth and Chebyshev Filters Analog to Analog Frequency Transforms Solution to problems UNIT 6 : FIR Filter Design Introduction to FIR filters Design of FIR Filters using Rectangular and Hamming window Design of FIR Filters using Bartlet and Hamming window F IR filter design using frequency sampling technique UNIT 7 : Design of IIR Filters from Analog Filters Impulse Invariance Method Mapping of Transfer Functions Approximation of derivative (Backward Difference and Bilinear Transforms) method Matched Z transforms Verification for stability and Linearity during mapping Solution to problems PAGE NO. 1-12
EN
TU D
TS
CI
6.3 6.4 VII 7.1 7.2 7.3 7.4 7.5 7.6
CITSTUDENTS.IN
TS .IN
13-25 26-31 32-52 53-70 71-108 109-128
Page 3

VIII 8.1 8.2 8.3 8.4 8.5 8.6 UNIT 8 : Implementation of Discrete time systems Structure of IIR and FIR systems Direct form I and direct form II Systems Cascade Realization Lattice Realization Parallel Realization Solution to problems
10EC52
129-156
CI
TS
TU D
CITSTUDENTS.IN
Page 4
EN
TS .IN

UNIT 1
10EC52
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.
CI
UNIT 1
CITSTUDENTS.IN
Page 5
TS
TU D
EN
TS .IN
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
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.
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
CI
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:
TS
FREQUENCY DOMAIN SAMPLING AND RECONSTRUCTION OF DISCRETE TIME SIGNALS :
TU D
X w
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
EN
xn e
jw n
DFT is an important mathematical tool which can be used for the so ftware
CITSTUDENTS.IN
TS .IN
periodic in the time domain DTFS representation can be used, in the frequency domain the
(1.1)
Page 6
10EC52
fundamental frequency range will be necessary. For convenience we take N equidist ant samples in the interval (0<=w<2 ). The spacing between samples will be w below in Fig.1.1. X[w] 2 as shown N
TS .IN
2 1) . (1.2)
2N 1 j 2 kn / N n N
Fig 1.1 Frequency Domain Sampling
If we evaluate equation (1) at w X 2 k N xne

n j 2 kn / N
2 k N
We can divide the summation in (1) into infinite number of summations where each sum contains N terms.
TS
2 k N
TU D
1 N 1
.......
xne
j 2 kn / N n 0
EN
k 0,1,2,......., ( N xne
j 2 kn / N
Let us first consider selection of N, or the number of samples in the frequency domain.
xne
j 2 kn / N
lN N 1 n lN
CI
xne
j 2 kn / N
If we then change the index in the summation from n to n-l N and interchange the order of
summations we get:
2 k N
N 1
x n lN e
n 0 l
for k
0,1,2,......, ( N 1) .(1.3)
CITSTUDENTS.IN
Page 7
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
k 0
ck e j 2
kn / N
n 0,1,2,........, ( N 1)
With FS coefficients: 1 N
N 1
ck
xp n e
n 0
j 2 kn / N
Comparing the expressions in equations (1.4) and (1.3) we conclude the following: ck 1 2 X k N N k
Therefore it is possible to write the expression xp[n] as below:

N 1
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
CI
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:
TS
TU D
X
k 0
xp n
1 N
2 k e j 2 kn / N N
EN
CITSTUDENTS.IN
TS .IN
0,1,2,......., ( N 1) (1.4)
0,1,......., ( N 1) . (1.5)
n 0,1,....., ( N 1) . (1.6)
Page 8

x[n]
10EC52
n 0 L xp[n] N>=L No aliasing
0 xp[n]
EN
N Fig. 1.2 Signal Reconstruction
CI
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
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 2k/n
TS
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
Discrete Fourier Transform:
TU D
CITSTUDENTS.IN
TS .IN
n N N<L Aliasing n
Page 9
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= 2k/n, 0kN-1 -j2kn/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 X(k) = X(j) =2kn/N n=0 N-1 -j2kn/N x(n) 0kN-1
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)
CI
X(0) X(1) X(2) X(3) = 1 1 1 1
TS
ex; 4 pt DFT of the sequence 0,1,2,3 1 -j -1 j 1 -1 1 -1 1 j -1 -j
Solving the matrix X(K) = 6 , -2+2j, -2 , -2-2j
1.5 Relationship of Fourier Transforms with other transforms

CITSTUDENTS.IN
Page 10
TU D
EN
TS .IN

1.5.1 Relationship of Fourier transform with continuous time signal:
10EC52
Suppose that xa(t) is a continuous-time periodic signal with fundamental period T p= 1/F0.The signal can be expressed in Fourier series as
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 j2k/N for K= 0,1,2,..N-1 we obtain
CI
The above expression is identical to Fourier transform X() evaluated at N equally spaced frequencies k = 2k/N for K= 0,1,2,..N-1.
TS
TU D
EN
Page 11
CITSTUDENTS.IN
TS .IN
Where {ck} are the Fourier coefficients. If we sample xa(t) at a uniform rate Fs = N/Tp = 1/T,
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
CI
TS
TU D
CITSTUDENTS.IN
Page 12
EN
TS .IN
Digital Signal Processing Recommended Questions with solutions
10EC52
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 3n/8 Ans:
CI
Question 4
TS
Question 3 Compute the eight-point DFT circular convolution for the following sequence X3(n) = cos 3n/8
TU D
EN
CITSTUDENTS.IN
Page 13
TS .IN
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 jn d , Where 2k/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= 2k/n, 0kN Therefore X(j) = x(n) -j2kn/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) =2kn/N x(n) -j2kn/N 0kN-1 n=0 Question 5
Solution:-
CI
TS
TU D
EN
CITSTUDENTS.IN
TS .IN
Page 14

Question 6 Find the 4-point DFT of sequence x(n) = 6+ sin(2 n/N), n= 0,1,N-1 Solution :-
10EC52
Question 7
CI
Solution
TS
TU D
CITSTUDENTS.IN
Page 15
EN
TS .IN
10EC52
Question 8
Solution
CI
TS
TU D
CITSTUDENTS.IN
Page 16
EN
TS .IN
10EC52
UNIT 2
PROPERTIES OF DISCRETE FOURIER TRANSFORMS (DFT) CONTENTS:1. 2. 3. 4.
MULTIPLICATION OF TWO DFTS- THE CIRCULAR CONVOLUTION, ADDITIONAL DFT PROPERTIES USE OF DFT IN LINEAR FILTERING
OVERLAP-SAVE AND OVERLAP-ADD METHOD.
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.
CI
TS
TU D
1. DIGITAL SIGNAL PROCESSING PRINCIPLES ALGORITHMS & APPLICATIONS, PROAKIS &
EN
CITSTUDENTS.IN
TS .IN
Page 17
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:-
CI
2.1.2 Linearity: If
TS
TU D
Page 18
CITSTUDENTS.IN
EN
TS .IN

Then A x1 (n) + b x2 (n) a X1(k) + b X2(k)
10EC52
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) +nok Wn Consider X(k) x(n)
TS
x(n) WN
X(k+no) N-1 kn x(k) = x(n) W n n=0 N-1
CI
X(k+no)
X(k+no)= \ x(n) WN n=0 kn non = x(n) WN WN non
2.1.6 Symmetry:
CITSTUDENTS.IN
Page 19
TU D
(k+ no)n
EN
TS .IN

For a real sequence, if x(n) X(N-K) = X* (k) For a complex sequence DFT(x*(n)) = X*(N-K) X(k)
10EC52
If x(n)
then
X(k)
2.2 Convolution theorem;
If y(n) = x(n)
h(n) then Y(k) = X(k) H(k)
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))
CI
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
TS
TU D
Ex let x(n) = 1,2,2,1 and h(n) = 1,2,2,1 Then y (n) = x(n) h(n)
EN
Circular convolution in time domain corresponds to multiplication of the DFTs
CITSTUDENTS.IN
TS .IN
Page 20
Real and even Real and odd Odd and imaginary Even and imaginary
real and even imaginary and odd real odd imaginary and even

X(n) = 1/N [DFT(X*(k)]*
10EC52
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.
In order to convolve a short duration sequence wit h 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
represented as
CI
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 fro m each output block must be overlapped and added to the first M-1 points of the succeeding
TS
TU D
convolution is performed to each block then the results are added. These data blocks may be
EN
2.3.1 Overlap and add
CITSTUDENTS.IN
TS .IN
Page 21
10EC52
block. Hence this method is called the overlap method. This overlapping and adding yields the output sequences given below.
CI
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
CITSTUDENTS.IN
TS
TU D
EN
TS .IN
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 Npoint 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
CI
TS
TU D
CITSTUDENTS.IN
Page 23
EN
from linear convolution and as a consequence we get
TS .IN

Input signal
10EC52
f----- L ---lf---L--1-o----L-----l
Output signal
Discard
M-1 points
/
Discard M-1
points M-1
/
Discard points
CI
TS
TU D
CITSTUDENTS.IN
Page 24
EN
TS .IN
f1pn: 5.10 Linear FIR filtering by the overlap-save method.

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
10EC52
Time shift: If x (n) X (k) mk WN X (k)
Then x (n-m) Question 2
Solution
CI
(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))
TS
State and Prove the: (i) Circular convolution property of DFT; (ii) DFT of Real and even sequence.
TU D
EN
Page 25
CITSTUDENTS.IN
TS .IN

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 Real and even Real and odd Odd and imaginary Even and imaginary X(k) real and even imaginary and odd real odd imaginary and even
10EC52
Distinguish between circular and linear convolution Solution
CI
Question 4
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.
TS
TU D
Question 3
EN
Page 26
CITSTUDENTS.IN
TS .IN
10EC52
Solution(a)
Solution(b)
Solution(c)
CI
Solution(d)
TS
TU D
CITSTUDENTS.IN
Page 27
EN
TS .IN
10EC52
Question 5
Solution
Question 6
CI
Solution
TS
TU D
CITSTUDENTS.IN
Page 28
EN
TS .IN
10EC52
CI
TS
TU D
CITSTUDENTS.IN
Page 29
EN
TS .IN
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).
CI
TS
TU D
EN
CITSTUDENTS.IN
TS .IN
Page 30
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
3.2 Need for efficient computation of DFT (FFT Algorithms)
TU D
Let us start with the simple way. Assume that
table for the N of interest. How big should the table be? so we just need to tabulate the N values:
EN
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
(Possibly even less since Sin is just Cos shifted by a quarter periods, so we could save just Cos
CI
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
TS
CITSTUDENTS.IN
TS .IN
has been precompiled and stored in a is periodic in m with period N,
Page 31
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 computationall y 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,
CI
problem is solved.
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
TS
TU D
EN
CITSTUDENTS.IN
TS .IN
Page 32

Question1
10EC52
Solution:-
Solution:-
CI
TS
TU D
Page 33
Question 2
CITSTUDENTS.IN
EN
TS .IN
10EC52
Question 3
Solution:-
CI
Question 4
TS
TU D
CITSTUDENTS.IN
Page 34
EN
TS .IN

Solution:- (a)
10EC52
(b)
z-plane
TU D
cir'Cle of radius 0.8
X(k)
EN
X( z)l z:;:z.
TS
Ez(n) [o.se'[lfk+tl]
t :O
f(n) =
zo(n)O.se-if"
CI
CITSTUDENTS.IN
TS .IN
n
-n
Page 35
Digital Signal Processing UNIT 4
10EC52
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
CI
TS
TU D
EN
CITSTUDENTS.IN
TS .IN
Page 36
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
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
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
that goes back many years.
4.2.1 Decimation in time
CI
Fig 4.1 First step in Decimation-in-time domain Algorithm
TS
TU D
approach. is to split the sequence into two sequences of length N=2. This is a very clever trick
EN
Useful when N is a power of 2: N = rv for integers r and v. r is called the radix, which
CITSTUDENTS.IN
TS .IN
DFT computation results in N frequency components. We know that DFT computation
Page 37
10EC52
N = 8-point decimation-in-ti me FFf algorithm

Stage 1 Stage 2 Stage 3 X[O]
X[l ]
EN TU D TS CI
Each clot repre ents a complex addition. Each arrow repre ents a complex multiplication.
CITSTUDENTS.IN
TS .IN
X[2]
X[3]
X[4]
X[5]
X[6]
X[?]
Page 38

4.2.2 Decimation-in-frequency Domain
10EC52
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.
Now, let us split X(k) into the even and odd-numbered samples. Thus we obtain
CI
TS
TU D
CITSTUDENTS.IN
Page 39
EN
TS .IN
Thus we obtain
10EC52
CI
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.
TS
TU D
CITSTUDENTS.IN
Page 40
EN
TS .IN
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 = log2 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
CI
TS
TU D
EN
CITSTUDENTS.IN
Page 41
TS .IN
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
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)
CI
TS
TU D
unique set of two-tone signals, called DTMF signals. These signals are processed at exchange
EN
CITSTUDENTS.IN
TS .IN
Page 42
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 fundament al 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 )
n 0
nk x(n)W N
(1)
N 1
N 1 mk x( m)WN m 0
X ( k ) WN
kN m 0
x( m)WN k ( N
TU D
x(m)WN k ( n
m)
This is in the form of a convolution

N 1
yk (n)
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 = 2k/N
CI
The filter has frequency response given by H k ( z) 1 1 WN k z

1
TS
(6)
The above form of filter response shows it has a pole on the unit circle at the frequency k = 2k/N. Entire DFT can be computed by passing the block of input data into a parallel bank of N single-pole filters (resonators)
EN
m)
Since
WN kN
is equal to 1, multiplying both sides of the equation by this results in; (2)
yk ( n )
(3)
CITSTUDENTS.IN
TS .IN
x ( n) hk ( n )
Page 43
10EC52
The above form of filter response shows it has a pole on the unit circle at the frequency k = 2k/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;
The form shown in eq (7) requires complex multiplications which can be avoided k z 1 ). Then frequency response of doing suitable modifications (divide and multiply by 1 WN the filter can be alternatively expressed as
1
vk (n) yk (n)
2 cos(2 k / N )vk (n 1) vk (n 2) vk (n) W v (n 1)

k N k
TU D
vk ( 1)
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 x(n) (9) (10)
CI
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
TS
EN
2
H k ( z)
1 WNk z 1 1 2 cos(2 k / N ) z
(8)
vk ( 2) 0
CITSTUDENTS.IN
TS .IN
Yk ( z ) 1 H k ( z) k 1 X (z 1utp W z X(k) N is The d es ir) ed o ut = yk(n) for k = 0,1,N-1. The phase factor appearing in the k diffy ere c e e qua t ion can be c o m puted nc1 e) and WN y k (n 1) x( n) yko( 0 stored. (7) k ( n)
Page 44
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 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
CI
N 1
For the above two problems Chirp z-transform is the alternative.
Chirp z- transform is defined as: X ( zk )

n 0
TS
scheme. x(n) zk n
only 128 frequencies out of 1024 there will be big wastage of computations in this
TU D
k 0,1,...... L 1
The possible solution is to multiply the input sequence by a -n
EN
(11)
CITSTUDENTS.IN
TS .IN
Page 45

Where zk is a generalized contour.
Zk is
10EC52
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, Note that, zk r0e j 0 ( R0e j 0 ) k k 0,1,.... L 1 (12)
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
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.
CI
TS
TU D
CITSTUDENTS.IN
Page 46
EN
in freq.))
TS .IN
10EC52
Substituting the value of zk in the expression of X(zk)

N 1 N 1
X ( zk )
n 0
x(n) z k n
n 0
x(n)(r0e j 0 ) nW
nk
where
R 0e j 0
(14)
4.2.2 Expressing computation of X(zk) as linear filtering operation: By substitution of nk
X ( zk ) W Where
TS
k2 /2
1 2 2 (n k (k 2 we can express X(zk) as y (k )
TU D
n) 2 ) (15) y(k ) / h(k ) g (n) x(n)(r0e j 0 ) n W (17)
EN
k 0,1,......... .L 1 (16)
n2 / 2
CI
h( n) W n y( k )
2
/2
N 1
g (n)h(k n)
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
n 0
CITSTUDENTS.IN
TS .IN
(13)
Page 47
10EC52
signal, which increases linearly with time. Such signals are used in radar systems are called chirp signals. Hence the name chirp z-transform.
CI
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
TS
TU D
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
EN
CITSTUDENTS.IN
TS .IN
Page 48
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
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)
CI
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.
TS
CZT is used in this application to sharpen the resonances by evaluating the z-transform
TU D
4.3.1 Advantages of CZT
EN
In general the computational complexity of CZT is of the order of M log 2 M 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.
CITSTUDENTS.IN
TS .IN
Page 49
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.
CI
TS
TU D
CITSTUDENTS.IN
Page 50
EN
TS .IN

4.5 Implementation of CZT in hardware to compute the DFT signals
10EC52
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
CI
TS
TU D
CITSTUDENTS.IN
Page 51
EN
TS .IN

Question 1
Compute the 16-point OFT of the sequence x(n) =cos
7f
10EC52
"2'1
0 5 n 5 15
using the radix-4 decimation-in-time algorithm.
Solution:1
1
-1
1
1
A=
A !.1 {
!_2 {
[ 1
-J ]
-1 -1 -)
-]
)T
z(l) z (S) z(9) z(13) (
[ z(2)
z ( 6) z(lO) z (14) ]T
!A( z(3) %(7) z(ll) z(l5) ] T
[ ml
2) F(6 F(lO)
F(l2)
CI
[ )) F(7 3
TS
F( (5 l) ] F(9) . [ F(13)
l [ l [l
= dl3 =
0
0
TU D
= &.1
=[
= &3 =
[0 0 0
l l
-0 4 ]
EN
Page 52
z(O)
z(4) z(8) z(l2)
CITSTUDENTS.IN
TS .IN

F(l4) F(ll) [ F(l5)
10EC52
= &:.. =
o
0
CI
TS
TU D
CITSTUDENTS.IN
Page 53
EN
TS .IN
10EC52
Solution :- There are 20 real , non trial multiplications
CI
Figure 4.1 DIF Algorithm for N=16
TS
TU D
CITSTUDENTS.IN
Page 54
EN
TS .IN
Question 2

Question 3
10EC52
Solution:-
Question 4
Solution:-
CI
TS
TU D
CITSTUDENTS.IN
Page 55
EN
TS .IN

Question 5
10EC52
Solution:-
Question 6
Solution:-
CI
This can be viewed as the convolution of the N-length sequence x(n) with implulse
response of a linear filter
TS
TU D
Page 56
CITSTUDENTS.IN
EN
TS .IN
10EC52
CI
TS
TU D
CITSTUDENTS.IN
Page 57
EN
TS .IN
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.
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.
CI
TS
TU D
EN
CITSTUDENTS.IN
TS .IN
Page 58
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
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
frequencies.
CI
TS
as being specified in terms of angle around the unit circle rather than in terms of analog
TU D
EN
by means of periodic sampling. The specifications for both analog and digital filters are often
CITSTUDENTS.IN
TS .IN
Page 59
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
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 causalit y 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
CI
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:
TS
TU D
of a tolerance scheme, as shown in Fig. 5.1.
EN
CITSTUDENTS.IN
TS .IN
Page 60
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
CI
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.
TS
In transforming an analog filter to a digital filter we must therefore obtain either H(z)
TU D
EN
Page 61
CITSTUDENTS.IN
TS .IN
Digital Signal Processing 5.2 Characteristics of Commonly Used Analog Filters:
10EC52
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:
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)H a(-s) occur on a circle of radius c at equally spaced points. From Eq.
CI
(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
TS
TU D
EN
Low pass Butterworth filters are all - pole filters with monotonic frequency response in
CITSTUDENTS.IN
TS .IN
Page 62
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,
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
CI
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
TS
TU D
EN
CITSTUDENTS.IN
TS .IN
Where by definition, 2 = 1/1+2. Thus the Butterworth filter is completely
Page 63
10EC52
Where T0(x) = 1 and T1(x) = x. At the band edge frequency = p, we have
Or equivalently
Where 1 is the value of the pass band ripple.
CI
And minor axis
The poles of Type I Chebyshev filter lie on an ellipse in the s-plane with major axis
Where is related to according to the equation
TS
TU D
Figure 5.2: Type I Chebysehev filter characteristic
EN
CITSTUDENTS.IN
Page 64
TS .IN
10EC52
The angular positions of the left half s-plane poles are given by
Where k = r2 Cos k and k = r1 Sink. The order of the filter is obtained from
Where, by definition 2 = 1/1+2.
Finally, the Type I Chebyshev filter is given by
CI
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
TS
A Type II Chebyshev filter contains zero as well as poles. The magnitude squared response is given as
TU D
EN
Page 65
CITSTUDENTS.IN
TS .IN
Then the positions of the left half s-plane poles are given by

Where
10EC52
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
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
CI
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
TS
TU D
Frequency transforms are used to transform lowpass prototype filter to other filters like
EN
CITSTUDENTS.IN
TS .IN
Page 66
10EC52
Thus we obtain
The filter function is
CI
TS
TU D
Page 67
CITSTUDENTS.IN
EN
TS .IN
10EC52
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 Approximation of derivatives Bilinear transformation technique Impulse invariant technique
CI
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
TS
TU D
Figure 5.8: Frequency response plot of the example
EN
CITSTUDENTS.IN
Page 68
TS .IN
10EC52
Where 2 is the gain at the stop band edge frequency s.
Then the filter transfer function in the analog domain is
CI
where Ak's are partial fractions coefficients of Ha(s). Finally, the transfer function of the digital filter is
TS
TU D
CITSTUDENTS.IN
Page 69
EN
Order of filter N =5. Then the 5 poles on the Butterworth circle of radius c = 0:5 are given by
TS .IN

b)
10EC52
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
CI
Finally, the transfer function of the digital filter is
TS
TU D
EN
Page 70
CITSTUDENTS.IN
TS .IN
10EC52
Question 2 Design a digital filter using impulse invariant technique to satisfy following characteristics (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
CI
TS
TU D
CITSTUDENTS.IN
Page 71
EN
TS .IN
(i) Equiripple in pass band and monotonic in stop band

c52 = 10-{).'7S = 0.177
Hence
10EC52
[( j 1 - (0.177 ) 2 + J1 - (0.177 ) 2 (1
+ 0.9952)) /0.9976 X
1]
0.177
j
]o!T[os . . 0 0.:>.
+ v/(o.. '75)2 _ 0.5..
2.4 "' 3
The order of fiJter l K = 3.
The 3 poles on the ellipse are det ermined by
EN
1
0. 9976
N
[
(3 =
V1
+ E2 + 1
E
[V1+ 0.99762
TU D TS
r2 -
n,
{32
2/3
X -'-----'-----
0 .01f
1.639
n,
(32- 1
2/3
X
0 .01f
CI
0.469
The angles,
. 1r
k =
"2 +
(2k + 1)n , k JV 2
The poles ar e at
CITSTUDENTS.IN
TS .IN
+ l. t
j3
=
1.342
+1
(1.341)2 + 1 2 X 1.341
(1.341)2 - 1 2 X 1.341
0,1 ,2
Page 72
10EC52
s0
0.469cos(
) +jl.639sin ( 11' )
- 0.2345 + j1.419
St
0.469 cos(n ) + jl.639sin(n) -0.469 + jO.O
0. 469 cos(
!!_ )+ j1.6939 sin( _!!_) 6 6
-0.2345 - jl.419 The analog filter transfer function is given by
where Ak 's are t he par t ial fraction coefficients. Finally. the digital filter t ransfer function is given by
CI
TS
TU D
Ak k=t (s - s k )
s + 0.2345 - j1.419)(s + 0.469) + (s + 0.2345 + jl.419)
EN
1
CITSTUDENTS.IN
TS .IN
Page 73

Question 3
10EC52
Solution:For the design specifications we have
CI
TS
TU D
CITSTUDENTS.IN
Page 74
EN
TS .IN
10EC52
Question 4
Solution:-
CI
TS
TU D
Page 75
CITSTUDENTS.IN
EN
TS .IN
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
5. DISCRETE TIME SIGNAL PROCESSING, OPPENHEIM & SCHAFFER, PHI, 2003.
CI
6. DIGITAL SIGNAL PROCESSING, S. K. MITRA, TATA MC-GRAW HILL, 2ND EDITION, 2004.
TS
MONALAKIS, PEARSON EDUCATION, 4TH EDITION, NEW DELHI, 2007.
TU D
EN
CITSTUDENTS.IN
TS .IN
Page 76
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
H ( z) 1
k 0 N
(1) ak z
k
k 1
2) Difference Equation representation;

N M
ak y(n k )
k 0
CI

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
TS
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)
TU D
bk x(n k ) (2)
k 0
EN
Page 77
bk z
CITSTUDENTS.IN
TS .IN

The factors considered with filter implementation are , a. It must be a simple design
10EC52
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:
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
CI
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
TS
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
TU D
EN
Out put is a function of past o/p, present and past i/ps 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
CITSTUDENTS.IN
TS .IN
Page 78
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)
7. Central Difference filter
y(n)= 1/2[ x(n) x(n-2)]
CI
TS
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:
TU D
y(n) = 1/3[x(n)+x(n-1)+x(n-2)]
EN
Page 79
CITSTUDENTS.IN
TS .IN

Ex. 1 2 3 order 0 0 1 a0 1 K 0 1 1/2 1/3 1/2 a1 1 -1 1/2 1/3 0 a2 -
10EC52
5(LP) 1 6(LP) 2 7(HP) 2
1/3 -1/2
6.6 Design of FIR filters:
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)) =

k 0
TS
TU D
- (2)
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.
EN
CI
Alternatively. it can be expressed in convolution form

M 1
y(n)
k 0
h(k ) x(n k )
i.e bk= h(k), k=0,1,..M-1 Filter is also characterized by
CITSTUDENTS.IN
TS .IN
bk x(n k ) -(1)
Page 80
4(HP) 1

M 1
10EC52
H ( z)
k 0
h( k ) z
-(3) polynomial of degree M-1 in the variable z -1. The roots of this
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
H ( z) h (0) h(1) z
1
h(2) z
.......... . h( M
2) z
If M is odd H ( z) h( M h(0) h(1) z 2) z

( M 2) 1
..........
h(
M 1 )z 2
M 1 2
TS .IN
( M 2)
h( M
1) z
( M 1)
h(
)z
M 1 2
h(
)z
M 3 2
.......... .
h( M 1) z
( M 1)
M 1 ) 2
h(0) z
M 1 ) 2
h(1) z
M 3 ) 2
EN
M 1 2 ) h( M 2
.......... .. h(
)z
h(
M 2
)z
.....h( M
1) z
M 1 ) 2
Applying symmetry conditions for M odd h(0) h (1) . . h( h( M M h( M h( M 1)
TS
1 2 2 ) h( M 2 1 ) h( M 2 1) h(0)
CI
. . h( M
TU D
2) 1 ) ) 3
Page 81
CITSTUDENTS.IN

( M 1 ) 2
10EC52
1 2n ) / 2 ( M 1 2n ) / 2
H ( z)
h(
M 1 ) 2
M 3 2 n 0
h(n){z ( M
similarly for M even

M (
H ( z)
M 1 ) 2
h(n){z ( M
1 2n ) / 2
( 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 1 M 1 ) 2 h(n) cos ( n) 2 2 n 0
M 3 2
H r (e j )
h(
H r (e j )
M 1 2
CI
H r (e j )
M 3 2 n 0
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 2 h( n ) sin ( M 2 1 n)
If M is even then,
TS
TU D
h( n) cos ( M 1
n 0
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 ( ) (
EN
n)
CITSTUDENTS.IN
TS .IN
Page 82

M j 2 1
10EC52
H r (e )
2
n 0
h(n) sin (
M 2
n)
In both cases the phase response is given by ( ) ( ( M 2 1 ) /2 if | H r (e j ) | 0
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
CI
TS
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.
TU D
EN
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 filter coefficients when M is odd and M/2 coefficients when M is even
CITSTUDENTS.IN
TS .IN
M 1 ) 2
3 / 2 if | H r (e j ) | 0
(M-1/2)
Page 83

6.6.6 Zeros of Linear Phase FIR Filters: Consider the filter system function
M 1
10EC52
H ( z)
n o
h( n ) z
Expanding this equation H ( z) h(0) h(1) z

1
h(2) z
h( M
2) z
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) then H ( z) H ( z) H ( z) h( M 1) h( M z z
( M 1) M 1 ( M 1)
h(0);
( M 2) 2)
2) z
1)
........ h(1) z h( M z 2) z ( M
[h( M 1) z ( M [
n 0
h(n)( z 1 ) n ]
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
CI
TS
TU D
This shows that if z = z1 is a zero then z=z1 -1 is also a zero
EN
( M 1)
H ( z 1)
CITSTUDENTS.IN
TS .IN
( M 2)
h( M 1) z
( M 1)
h(0) z
( M 1)
..... h(1) z h(0)]
Page 84
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
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 )
n
hd (n)e
j n
where hd (n) are hd (n) 1 2
fourier series coefficients
TS
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 M odd) h(n) hd (n) for ( M 1) / 2 0 otherwise n ( M 1) / 2
TU D
This expansion results in impulse response coefficients which are infinite in duration and non causal. It can be made finite duration by truncat ing 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)
EN
H d (e j )e j n d
CI
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)
CITSTUDENTS.IN
TS .IN
6.7.1 Design of Linear Phase FIR filter based on Fourier Series method:
( assuming
Page 85
Digital Signal Processing Exercise Problems

Problem 1 : Design an ideal bandpass filter with a frequency response: H d (e j ) 1 for 3 4
10EC52
4 0 otherwise Find the values of h(n) for M = 11 and plot the frequency response.
1 2 1
/4
3 /4
e j nd
3 /4
e j nd
/4
n n sin n 4 n 4 truncating to 11 samples we have h(n) hd (n) for | n | 5 sin 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
CI
h(3)=h(-3)=0 h(4)=h(-4)=0 h(5)=h(-5)=0
h(1)=h(-1)=0
h(2)=h(-2)=-0.3183
The transfer function of the filter is
TS
TU D
3
EN
Page 86
hd (n)
1 2
H d (e j )e j n d
CITSTUDENTS.IN
TS .IN

( N 1) / 2
10EC52
H ( z)
h(0)
n 1
[h(n ){z n
z n }]
0.5 0.3183( z 2 H ' ( z)
z 2) z 2 )] 0.3183z h' (2)

7
the transfer function of the realizable filter is z 5 [0.5 0.3183( z 2 0.3183z h ' (0) h' (3) h' (5) h' (10) h' (7) 0.5
3
0.5 z
the filter coeff are
The magnitude response can be expressed as

( N 1) / 2
| H (e j ) |
n 1
a(n) cos n
comparing this exp with

5
| H (e j ) | | z 5 [h(0) 2
n 1
h(n) cos n] |
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
CI
The magnitude response function is |H(e j)| = 0.5 0.6366 cos 2 which can plotted for various values o f
in degrees =[0 20 30 45 60 75 90 105 120 135 150 160 180];
TS
TU D
Page 87
CITSTUDENTS.IN
EN
TS .IN
h' (1) h' (9) 0.3183
h' (8)
h' (4)
h' (6)
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 0 for 2 2
Find the values of h(n) for N =11. Find H(z). Plot the magnitude response
CI
h(0) = 1/2 h(1)=h(-1)=0.3183 h(2)=h(-2)=0 h(3)=h(-3)=-0.106
From the freq response we can determine hd(n), n sin /2 1 2 n and e j nd hd ( n ) n 2 /2 Truncating hd(n) to 11 samples
TS
TU D
for 2
EN
n 0
CITSTUDENTS.IN
TS .IN
Page 88

h(4)=h(-4)=0 h(5)=h(-5)=0.06366
10EC52
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];
H ' (z) 0.06366 0.106z
2
0.3183z
0.5 z
0.3183z
0.106z
0.06366z
10
Using the result of magnitude response for M odd and symmetry

M 3 2 M 1 M 1 h(n) cos ( H r (e j ) [h( ) n)] 2 2 n 0 j | H r (e ) | | [0.5 0.6366cos w 0.212cos3 w 0.127cos5 w] |
Problem 3 :
Design an ideal band reject filter with a frequency response: H d (e j ) 1 0 for and 2 3
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
TU D
3 otherwise
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.,
CI
TS
EN
CITSTUDENTS.IN
TS .IN
-0.1378 0];
Page 89
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 W (e j ) sin( M / 2) 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.
CI
TS
TU D
CITSTUDENTS.IN
Page 90
EN
TS .IN
10EC52
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.
CI
TS
TU D
EN
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.
CITSTUDENTS.IN
TS .IN
Page 91

6.8.2 What is ideal window characteristics?
10EC52
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
whan (n)
TS
0.5(1 cos
6.8.4 Hanning windows: It is defined mathematically by,
CI
TU D
2 n ) for 0 n M 1 M 1
Page 92
CITSTUDENTS.IN
EN
TS .IN
10EC52
6.8.5 Hamming windows: This window function is given by,
6.8.6 Blackman windows: This window function is given by, 2 n wblk (n) 0.42 0.5 cos M 1
CI
TS
TU D
0.08cos 4n for 0 n M 1 M 1
EN
Page 93
wham(n)
0.54 0.46cos
2n for 0 n M 1
CITSTUDENTS.IN
TS .IN
M 1
10EC52
6.8.7 Bartlett (Triangular) windows: The mathematical description is given by, M 1 | 2 M 1
2|n wbart (n ) 1
for 0 n
6.8.8 Kaiser windows: The mathematical description is given by,
TS
I0 wk (n)
TU D
M 1 2 I0
2
EN
2
n M 1 2
M 1 2 for 0 n M 1
CI
CITSTUDENTS.IN
TS .IN
M 1
Page 94
10EC52
Type of window
Appr. Transition width of the main lobe
Peak
sidelobe (dB)
Bartlett
8/M 8/M 8/M
Hanning Hamming Blackman
TU D
12/M
CI
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 flexibilit y 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
TS
EN
-27 -32 -43 -58
Page 95
Rectangular
4/M
-13
CITSTUDENTS.IN
TS .IN

( assuming M odd) choosing proper window
10EC52
h(n) hd (n) w( n) where w( n) is the window function defined for ( M 1) / 2 n ( M 1) / 2
3. 4. 5.
Introduce h(n) = h(-n) for linear phase characteristics Write the expression for H(z); this is non-causal realization 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: Hd (e j ) 1 0 for | |
using a hanning window with M = 11 and plot the frequency response.
CI
1 [ 2
/4
TS
hd ( n) e j nd e j nd ]
/4
TU D
EN
4 4
Page 96
CITSTUDENTS.IN
TS .IN

hd (n) hd (0) 1 [sin n sin
/4
10EC52
n 0.75 and n 0
n 1 [ 2
n ] 4 d ]
for 3 4
d
/4
whn (n)
0.5 0.5 cos 0
2 n M 1 otherwise
M 1 ) 2
M 1 ) 2
for N whn (n)
11 0.5 0.5 cos n 5 5 n 5
h(n)= whn(n)hd(n)
CI
h' (n) h(n 5)
2
h(n)=[0 0 -0.026 -0.104 -0.204 0.75 -0.204 -0.104 -0.026 0 0]
TS
TU D
0.104z
3
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
EN
4
H ' ( z)
0.026z
0.204z
0.75z
Using the equation
CITSTUDENTS.IN
TS .IN
5
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
0.204z
0.104z
0.026z
Page 97

M 3
10EC52
H r (e jw ) H r (e jw )
2 M 1 M 1 [ h( ) 2 h(n) cos ( n) 2 2 n 0 4
0.75) 2
n 0
h(n) cos (5 n)
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]
CI
TS
TU D
CITSTUDENTS.IN
Page 98
EN
TS .IN

Prob 2 : Design a filter with a frequency response:
H d (e j ) e 0
j3
10EC52
for | |
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. hd (n) 1 2
/4
e
/4
j3
e j nd
hd (1) hd (5) hd (2) hd (4) hd (3) 0.25
0.159 0.22
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
CI
h(n)=hd(n) whn(n) h(n)=[0 0.03975 0.165 0.25 0.165 0.3975 0]
whn(3)=1
TS
TU D
EN
Page 99
(n 3) 4 (n 3) this gives hd (0) hd (6) 0.075
sin
CITSTUDENTS.IN
TS .IN
10EC52
6.9 Design of Linear Phase FIR filters using Frequency Sampling method:
Fro DFT relationship we have h(n) H (k ) 1 N

N 1
H ( k )e j 2
k 0
N 1
TS
h(n)e
n 0 N 1
Also we know H(k) = H(z)|z=e j2kn/N
CI
H ( z) h(n) z
n n 0
The system function H(z) is given by
Substituting for h(n) from IDFT relationship H ( z) 1 z N

N N 1
k 0
TU D
kn / 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.
for n
j 2 kn / N
for k
0,1,......... N 1
H (k ) 1 e j2
kn / N
EN
0,1,...... N 1
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).
CITSTUDENTS.IN
TS .IN
Page 100
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 h ( n) 1 N
N 1
H ( k )e j 2
k 0 N 1
kn / N
can be rewritten as (for N odd)
h( n) h( n)
1 H (0) N 1 H (0) N
H (k )e j 2 kn / N
k 1 N 1/ 2 N 1
H (k )e j 2 kn / N
k 1
h( n) h( n) h( n) h( n)
1 H (0) N 1 H (0) N 1 H (0) N 1 H (0) N
TU D
( N 1) / 2 ( 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 H ( k )e j 2 H ( k )e j 2 H ( k )e j 2
kn / N
k 1
( N 1) / 2
kn / N
k 1
TS
( N 1) / 2 k 1 ( N 1) / 2 k 1 ( N 1) / 2
kn / N k 1 kn / N
CI
h( n) h ( n) 1 H (0) N
( H ( k )e j 2
1 H (0) 2 N
Re( H (k )e j 2 kn / N
k 1
Similarly for N even we have

( N 1) / 2
2
k 1
Re( H (k )e j 2
EN
H ( k )e j 2
kn / N k N 1/ 2
H (N
k 1
( N 1) / 2
H * ( k )e
k 1
( N 1) / 2
( H ( k )e j 2
( H (k )e j 2 kn / N ) *
kn / N
CITSTUDENTS.IN
TS .IN
k )e
j 2 kn / N j 2 kn / N kn / N *
Page 101
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.
H d (e ) 0 with M
j (
M 1 ) 2
for | c
| /2
otherwise 17 and
for
/2
Selecting H (k )
H d (e j ) |
2 k 8 17
CI
0 k 4 and 5 k 8
The range for k can be adjusted to be an integer such as
TS
H (k )
2 k 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 e

j
TU D
for k
2 k 17
2 k M
2 k 17
for 0
EN
0,1,...... 16
Page 102
H d (e j )
j 8
for 0
/2
CITSTUDENTS.IN
TS .IN
The desired response can be expressed as

The freq response is given by
j
10EC52
H (k )
e 0
2 k 8 17
for 0 5 k
k 8
for
Using these value of H(k) we obtain h(n) from the equation

( M 1) / 2 1 h( n) Re( H (k )e j 2 kn / M )) ( H (0) 2 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) ) for n h(n) ( H (0) 2 cos( 17 17 k 1
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.
CI
TS
TU D
EN
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)
CITSTUDENTS.IN
TS .IN
0,1,........ 16
Page 103
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 hd (n) between to + 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,
H ' ( z) 0.33 0.5 z
1
Also from the equation

( M 3) / 2
H r (e )
2
n 0
For M=7 and h(n) as found above we obtain this as H r (e j ) 0.66sin 3 sin 2 2 sin sin 2 2 sin )
CI
H (e j )
b) Hamming window h(n)=hd(n)wh(n) where wh(n) is given by
TS
jH r (e j )
TU D
h( n) sin ( M 1 2 n) j (0.66sin 3
0.5 z
EN
5
0.33z
CITSTUDENTS.IN
TS .IN
6
cos n 1 j ej n d and n n 2 n The hd(n) is an add function with hd(n)=-hd(-n) and hd(0)=0
Page 104

wh (n) 0.54 0.46 cos 0 otherwise For the present problem n 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] wh ( n) 0.54 0.46cos 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 differentiator as H (e j ) jH r (e j ) j (0.0534sin 3 0.31sin 2 1.54sin ) 2 n ( M 1) ( M 1) / 2 n ( M 1) / 2
10EC52
CI
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
TS
TU D
CITSTUDENTS.IN
Page 105
EN
TS .IN
write the freq response of
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 j 0
CI
(1 cos n) 1 0 jn je j n d ] [ je d n 2 0 At n = 0 it is hd(0) = 0 and hd(n) is an odd function hd (n)
TS
The impulse response is given by
TU D
0 n except n 0
Page 106
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]
CITSTUDENTS.IN
EN
TS .IN

H r (e j ) H (e j )
4
10EC52
2
n 0
h(n) sin (5 n) 0.424sin 3 1.272sin }
j | H r (e j ) | j{0.254sin 5
b) Blackman Window window function is defined as n 2 n wb (n) 0.42 0.5 cos 0.08cos 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 -5n5 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 ]
CI
TS
TU D
CITSTUDENTS.IN
Page 107
EN
TS .IN
10EC52
Recommended questions with solution

Question1
Solution:-
(b) Magnitude plot
CI
TS
TU D
CITSTUDENTS.IN
Page 108
EN
TS .IN

Phase plot
10EC52
(c) Hamming window
(d) Bartlett window
CI
TS
TU D
CITSTUDENTS.IN
Page 109
EN
TS .IN
10EC52
Question 2
Solution:-
CI
TS
TU D
Page 110
CITSTUDENTS.IN
EN
TS .IN
10EC52
CI
TS
TU D
CITSTUDENTS.IN
Page 111
EN
TS .IN

Question 3
10EC52
Solution:-
Magnitude and phase response
CI
TS
TU D
Page 112
CITSTUDENTS.IN
EN
TS .IN

Question 4
10EC52
tion:-
CI
TS
TU D
EN
CITSTUDENTS.IN
Page 113
TS .IN
Solu
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
TRANSFORMATION) METHOD,
MATCHED Z TRANSFORMS
5. VERIFICATION FOR STABILITY AND LINEARITY DURING MAPPING
RECOMMENDED READINGS:-
CI
TS
TU D
EN
CITSTUDENTS.IN
TS .IN
DIFFERENCE AND
BILINEAR
Page 114
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:
CI
Figure 7.1: Tolerance limits for approximation of ideal low-pass filter
TS
TU D
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.
EN
CITSTUDENTS.IN
TS .IN
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.
Page 115
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
7.2 Design of IIR Filters from Analog Filters:
CI
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
TS
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.
And the corresponding rational function is
TU D
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:
------------------------------------------------------------7.1
EN
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 causalit y 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.
CITSTUDENTS.IN
TS .IN
Page 116
10EC52
---------------------------------------------------------7.2
The corresponding description for digital filters has the form
--------------------------------------------------7.3 and the rational function
--------------------------------------------------------7.4
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,
CI
Or
------------------------------------------------------------------------- 7.5 Where T is the sampling period. Because of uniform sampling, we have
TS
TU D
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
EN
---------------------------------------------7.6
---------------------------------------------7.7
CITSTUDENTS.IN
TS .IN
Page 117
10EC52
Figure 7.2: Mapping of s-plane into z-plane
CI
Figure 7.3: Illustration of the effects of aliasing in the impulse invariance technique
TS
TU D
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
EN
CITSTUDENTS.IN
TS .IN
Page 118

analog filter as
10EC52
------------------------------------------------7.8 From the discussion of the sampling theorem it is clear that if and only if
Then
The corresponding impulse response is
And the unit-sample response of the digital filter is then
CI
TS
TU D
------------------------------------------------------------------ -----7.9
--------------------------------------------------------------- 7.10
EN
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
The system function of the digital filter H(z) is given by
------------------------------------------------------------7.12
CITSTUDENTS.IN
TS .IN
--------------7.11
Page 119
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 exp skT 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
Where y(n)=y(nT). Approximation to higher-order derivatives are obtained by repeated application of Eq. (7.13); i.e.,
EN
Applying Eqs. (7.13), (7.14) and (7.15) to (7.1), we obtain
TS CI
And
Where y(n) = ya(nT) and x(n) = xa(nT). We note that the operation (1)[ ] is a linear shiftinvariant operator and that (k)[ ] can be viewed as a cascade of (k) operators (1)[ ]. In particular
TU D
For convenience we define
-------------------------------------------------------------------7.15
---------------------------------------------7.16
CITSTUDENTS.IN
TS .IN
--------------------------7.13
-------------------------- 7.14
Page 120
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 Zplane, 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
CI
------------------------------------------------------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
TS
TU D
EN
CITSTUDENTS.IN
TS .IN
Page 121
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
TS
CI
Where ya(t) is the first derivative of ya(t). The corresponding analog system function is
We can write ya(t) as an integral of ya(t), as in
TU D
-----------------------------------------------------------7.20
EN
CITSTUDENTS.IN
TS .IN
Page 122

In particular, if t = nT and t 0 = (n - 1)T,
10EC52
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
TU D
From Eq. (7.22) it is clear that H(z) is obtained from Ha(s) by the substitution
TS CI
That is,
-------------------------------------------------------------------7.23
--------------------------------------------------------------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 123
CITSTUDENTS.IN
EN
--------------------------------------------7.22
TS .IN
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
CI
Thus for z on the unit circle, = 0 and and are related by T /2 = tan (/2) or = 2 tan -1(T /2)
TS
Figure 7.5: Mapping of analog frequency axis onto the unit circle using the bilinear Transformation
TU D
CITSTUDENTS.IN
Page 124
EN
TS .IN
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
CI
the corresponding digital filter is
TS
TU D
-----------------------------------------------------------------7.26
EN
---------------------------------------------------------7.27
CITSTUDENTS.IN
TS .IN
Page 125
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).

Question 1
CI
TS
TU D
CITSTUDENTS.IN
Page 126
EN
TS .IN
10EC52
Question 2
CI
TS
Question 3
TU D
Page 127
CITSTUDENTS.IN
EN
TS .IN

Question 4
10EC52
Question 5
CI
TS
TU D
CITSTUDENTS.IN
Page 128
EN
TS .IN

H(z)
1-e
10EC52
S;7 -1
Here
C,
H(z)
c;
1 1- t'
-T -1
62 0.577e- 12
t'
= = =
+1
-c
-O..'iT J0.866T
1-e-T z-1
+ -------:-::-::::----------:::---:----1- 2t..osr cos(0.866T)z-1 +e-T :-2
2(0.577) cos(-2.62)- 2(0.577)e-osr ;:
TS .IN
"
1 l-e-r:1
0.577{' J - ,_ -+ _e(-05+,0&.6)7z-' t----1 1-L' cu.s ,o N>>f -!

1- e
0.577e-12 ' 62
,t,
0. 77e'2 62 +-1 -05T J0.866T -1
cos(262- 0. 66T)
EN
z_ + -:
T 2 - l!Se"" : 2 - 2e-o.5I
_.,
Multiplying the numerator and denominator of'first term on RHS by z and by t for second term on RHS, above equation becomes,
5
r cos(+ 0.8667 z
l' -
H(z)
z-c
cos(0.866T): +
In terms of sampling interval T = 1, transfer function is,

H(z)
Where b 0
TU D
=
1
b0 z
-1
-1
+b 1 z-2
-a,z
-n 1z
- n 3:'
5
-2e-{l.5T cos(0.866T) + e-T + 1154e-{).<;T cos(
+ 0.866T ) = 1.0773
TS
11 : = 1
5 b: + 0.866T J0.12S4 1 = e r + 1154e tsr cos(
e- T + 2r-0.5T cos( 0.8661') = ..I538

15
CI
a2
113
-l'- T - 2et'-ZT
cos(0.866T) =- 0.657
= 0.1353
CITSTUDENTS.IN
Page 129

Question 6
10EC52
Des1gn a lmear phase high pass filter using the Hamming window for tlte following desired Jreq11.ency response. -i 3w n I c h
wliC
Hd (m) =
0
l ool<
2 N -1 )
oo (n) = 0.54 - 0.46 cos ( n:n ), where N is the lengtlz of the Hamming winaow.
Sol. :
_!_
27t -It
-n/n
1t
H J (ro)e jom d ro
= =
1 2rt -n
1
J -}
1'
EN
n(n- 3)
TU D
''d(3)
Also,
2(5: 5:):
+
Let
N = 7 Impulse response of FIR filter is :

h(n) = hd(n)w(n)
CI
TS
CITSTUDENTS.IN
TS .IN
3w c Jiwn t Iro +
1 J c-1 3we 21t rc 6 6
n:
<llud"'' \.U
[siniTI(11 - 3)] - sin[1t (n -
3)]] n 3
Page 130
10EC52
Question 7
CI
TS
TU D
CITSTUDENTS.IN
Page 131
EN
TS .IN

Given : T = 2 sec.
Q' p =
10EC52
n; =
N =
tan(ro; J tan( ta{i}tan( 0 ")=07265

log[(lO-kpflO -lJ/(lo-ks/10 2log (Q'p
0 1t ) = 0.3249
-1)]
/0.;)
= log(0.3335J30.6228) = . 2 8083 3 2/og(0.3249, 07265) Q'

= (10-Apf l0 -l) li2N
P
The required prewarped analog filter is obtained by applying lowpass to lowpass transforma bon.
TU D
ll 11(s)
111(s)ls--.
2
(s + s + l)(s +1) -+s

2 1
TS
CI
= =
EN
1
2
Referring to normalized Iowpass butterworth filter tables 1 H3(s) = ---(s2 + s + l)(s + 1)
(': + i +1)(+1 )
1
(:+';2)(';2)
1
2
( 2s
+84s+8 )(';2)
CITSTUDENTS.IN
TS .IN
=
= 1.7688 Page 132
10EC52
CI
TS
TU D
CITSTUDENTS.IN
Page 133
EN
TS .IN
10EC52
UNIT 8 Implementation of Discrete time systems

CONTENTS:1. IMPLEMENTATION OF DISCRETE-TIME SYSTEMS
3. 4.
DIRECT FORM I AND DIRECT FORM
II SYSTEMS,
CASCADE, LATTICE AND PARALLEL REALIZATION.
CI
TS
TU D
EN
RECOMMENDED READINGS:-
CITSTUDENTS.IN
TS .IN
Page 134
2. STRUCTURES FOR IIR AND FIR SYSTEMS
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)
k 1
k 1
b) Ration of polynomials
M
bk Z H (Z ) 1
k 1 k 0 N
ak Z
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.
CI
8.1.3 Finite-word-length effects

These effects refer to the quantization effects that are inherent in any digital implementation of the system, eit her 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.
TS
TU D
EN
The following factors influence choice of a specific realization, Computational complexity Memory requirements Finite-word-length Pipeline / parallel processing
CITSTUDENTS.IN
TS .IN
Page 135
a k y (n k )
bk x(n k )

8.1.4 Pipeline / Parallel Processing
10EC52
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)
k 0
bk x(n k )
M 1
Or equivalently, the system function H (Z )

k 0
bk Z
Where we can identify h( n)
bn
n 1
8.2.1 Direct Form Structure Convolution formula is used to express FIR system given by,
M 1
y( n)
k 0
h(k ) x(n k )
It is Non recursive in structure
CI
TS
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)
TU D
Page 136
CITSTUDENTS.IN
EN
0 otherwise Different FIR Structures used in practice are, 1. Direct form 2. Cascade form 3. Frequency-sampling realization 4. Lattice realization
TS .IN
10EC52
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 4 3 2 2 3 1 1 1 1 1 1 h(n) 1, , , , m=9 , , , 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)
CI
TS
TU D
CITSTUDENTS.IN
Page 137
EN
TS .IN
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 3 4 4 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
10EC52
8.2.2 Cascade Form Structure

K
The system function H(Z) is factored into product of second order FIR system H (Z )
k 1
H k (Z )
CI
TS
TU D
Where H k ( Z ) bk 0 bk 1 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.
EN
CITSTUDENTS.IN
TS .IN
Page 138
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 co mplex 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 * ) Ck 0 Ck 1 z
1
Ck 2 z
Ck 1 z
Problem: Realize the difference equation y(n) x(n) 0.25x(n 1) 0.5x(n 2) 0.75x(n in cascade form.
TU D
1
Y ( z) Soln: H ( z) H ( z)
X ( z ){1 0.25z (1 1.1219z

1
0.5 z
H ( z ) 1 0.25z
0.5 z
0.75z _ 3
1.2181z 2 )(1 1.3719z
H 1 ( z) H 2 ( z)
CI
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
TS
EN
3)
2
0.75z
4 1
CITSTUDENTS.IN
TS .IN
x(n 4) z 4) 0.821z 2 )
Page 139
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.
CI
Consider H ( z)
m Y ( z) 1 a m (i ) z i X ( z) i 1 m is the order of the FIR filter and a m(0)=1
when m = 1 Y(z)/ X(z) = 1+ a1(1) z-1
TS
TU D
EN
CITSTUDENTS.IN
Page 140
TS .IN

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)
10EC52
The outputs are f 1 ( n) f 0 ( n) k1 r0 (n 1) 1a 1b y ( n)
If m=2
CI
y( n) f 0 (n )
Substituting 1a and 1b in (2)
TS
k1r0 (n 1) k1r0 ( n r0 ( n) k1 x( n 1) f 0 ( n) f 0 (n ) x( n ) x( n) ( k1
Y ( z) 1 a 2 (2) z 2 1 a 2 (1) z X ( z) y (n) x(n) a 2 (1) x(n 1) a 2 (2) x(n 2) (2) y (n) f 1 (n) k 2 r1 (n 1)
TU D
k 2 [ k1 f 0 (n k 2 k1 f 0 (n k 2 k1 x( n 1) 1) k 2 x (n 1) 1) r0 (n k 2 r0 (n 2)] 2)] 2)] 1) x(n ) k 2 x( n 2) k1k 2 ) x(n
Page 141
r1 (n) k1 f 0 (n) r0 (n 1) if k1 a1 (1), then f 1 (n)
sin ce y( n)
We recognize
CITSTUDENTS.IN
EN
TS .IN

a 2 (1) a 2 (1) k1 k2 k 1k 2
10EC52
Solving the above equation we get k1 a 2 (1) 1 a 2 (2) and k 2 a2 (2) (4)
Similar to above, an Mth order FIR filter can be implemented by lattice structures with M stages
CI
km a m (m) a m 1 (i )
8.4.1 Direct Form I to lattice structure
For m = M, M-1, ..2, 1 do
TS
a m (i ) a m (m)a m (m i ) 2 1 km 1 i m 1
TU D
EN
CITSTUDENTS.IN
Page 142
TS .IN
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)
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 k2 2 and m 1 a 2 (2) & k1 a1 (1) 2&i 1 a 2 (1) a 2 (2)a 2 (1) 1 k a 2 (1) 1 a 2 (2)
2 2
for m a1 (1)
a 2 (1)[1 a 2 (2)] 1 a (2)

2 2
Thus k1
8.4.2 Lattice to direct form I For m = 1,2,.M-1 a m (0) 1 a m (m) k m a m (i )
a m 1 (i ) a m (m)a m 1 (m i ) 1
CI
Problem: Given FIR filter H (Z ) 1 2Z 1 Given a1 (1) 2 , a2 (2) 1 3 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) 1 3 if m=1 k1 a1 (1) also, when m=2 and i=1 a 2 (1) 2 3 a1 (1) 1 1 a 2 (2) 1 3 2 Hence k1 a1 (1) 3 2
TS
TU D
1 3
EN
i m 1 Z
2
obtain lattice structure for the same
CITSTUDENTS.IN
TS .IN
Page 143
a 2 (1) 1 a 2 (2)
10EC52

Problem:1 Consider an FIR lattice filter with co-efficients k1 filter ( H (Z )
a3 (0) 1
co-efficient a3 (0) a3 (1) Z 1 a3 (3) k3
for a3 (2) Z 2 1 4
the a3 (3)Z 3 )
1 1 , k3 , k2 3 2 direct
1 . Determine the FIR 4 form structure 1 3 1 2
for m=2, i=1 a2 (1) a1 (1) a2 (2)a1 (1) 1 1 1 2 3
= for m=3, i=1
4 6
2 3
a3 (1)
TS
CI
for m=3 & i=2
TU D
a 2 (1) a3 (3)a 2 (2)
2 1 1 . 3 4 3 2 1 8 1 = = 3 12 12 = 9 12 3 4
a3 (2) a2 (2) a3 (3)a 2 (1) 1 1 2 = . 3 4 3 1 1 2 1 = 3 6 6
EN
Page 144
= a1 (1)[1 a 2 (2)]
CITSTUDENTS.IN
TS .IN
a1 (1) k1
a2 (2)
k2

= 3 6 1 2
10EC52
a3 (0)
1,
a3 (1)
3 , 4
a3 (2)
1 , 2
a3 (3)
1 4
8.5 Structures for IIR Filters

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

N
y ( n)
k 1
a k y (n k )
Different realizations for IIR filters are,
CI
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.
TS
1. Direct form-I 2. Direct form-II 3. Cascade form 4. Parallel form 5. Lattice form
TU D
bk x ( n k )
k 0
EN
Page 145
CITSTUDENTS.IN
TS .IN
10EC52
8.5.2 Direct form-II The given transfer function H(z) can be expressed as, Y ( z ) V ( z) Y ( z) . X ( z) X ( z) V ( z) where V(z) is an intermediate term. We identify, H ( z) V ( z) X ( z) 1
N
-------------------all poles
k
1
k 1
ak z
v(n) y (n)
x(n)
k 1 M
v(n)
k 1
CI
TS
TU D
a k v( n k ) bk v(n 1)
M Y ( z) -------------------all zeros 1 bk z k V ( z) k 1 The corresponding difference equations are,
EN
Page 146
CITSTUDENTS.IN
TS .IN
10EC52
CI
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,
TS
TU D
CITSTUDENTS.IN
Page 147
EN
TS .IN

H (z) H 1 ( z ) H 2 ( z )....H k ( z )
10EC52
8.5.4 Parallel form structure
CI
Where H k ( Z ) bk 0 bk 1 Z 1 1 a k1 Z 1 a k 2 Z
2
TS
TU D
In the expression of transfer function, if N M we can express system function N N Ak H (Z ) C C H k (Z ) pk Z 1 k 11 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.
EN
CITSTUDENTS.IN
TS .IN
Where H k ( Z ) could be first order or second order section realized in Direct form II form i.e., bk 0 bk 1Z 1 bk 2 Z 2 H k (Z ) 1 a k1 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 complexconjugate 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 precisio n arithmetic.
Page 148
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
H (Z )
3 4
1
7 6 1
1 8
Z
1
1 Z
2
1 2
10 1 1
7 8
Z
3 32
1 3 2
1 2Z
1
1 Z Z
47 32 1
H (Z )
10 1 1
15 8
5 6 1
2Z
2
H ( z)
CI
TS
TU D
CITSTUDENTS.IN
( 14.75 12.90 z 1 ) 7 1 3 2 (1 z z ) 8 32
EN
(24.50 26.82 z 1 ) 1 2 (1 z 1 z ) 2
TS .IN
j1 2 Z
1 1
10 1
1 2
2 3
1 2Z 1
1 2
j1 2 Z
1 2
17 32
2 3 3
3 64
Page 149
10EC52
H1 ( z)
CI
TS
10(1 2 z 1 ) 1 2 1 z 1 z 2
TU D
Page 150
Cascade Form H(z) = H1(z) H2(z) Where 7 1 1 z 6 H1 ( z) 7 1 z 1 8
1 2 z 3 3 2 z 32
CITSTUDENTS.IN
EN
TS .IN

Parallel Form H(z) = H1(z) + H2(z)
10EC52
H(z) = (-14.75-12.90z-
1 )
+ (24.50+26.82z2
1 )
(1 +7 -z-1 +3 -z-2) 8 32
(1 -z -1 +1 -z-2)
Solution: The Direct form realization is done directly from the given i/p - o/p equation, show in below diagram
CI
Direct form -II realization Taking ZT on both sides and fmding H(z)
H (z) = Y(z) = 3 + 3.6z- + 0.6zX(z) 1 + 0.1z-1 - 0.2z-2
1 2
TS
TU D
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)
EN
CITSTUDENTS.IN
TS .IN
Page 151
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 where H 1 ( z ) 3 0.6 z 1 0.5 z
1 1
and H 2 ( z )
Parallel Form realization
CI
H ( z) 3
The transfer function can be expressed as H(z) = C + H1(z) + H2(z) where H1(z) & H2(z) is given by, 7 1 0.4 z 1 1 0.5 z
TS
1
TU D
1
EN
Page 152
1 z 1 1 0.4 z 1
CITSTUDENTS.IN
TS .IN
10EC52
8.6 Lattice Structure for IIR System:

Consider an All-pole system with system function. H (Z ) 1
k 1
1
N
a N (k ) Z
The corresponding difference equation for this IIR system is,

N
y (n)
k 1
a N (k ) y (n k )
N
OR x ( n) For N=1 x( n) Which can realized as, y (n)
CI
We observe
TS
x(n) y (n )
g1 (n )
TU D
f 1 (n) f 0 (n)
f1 ( n )
y( n) a1 (1) y(n 1)
x(n) k1 y(n 1) k1 f 0 ( n) g 0 (n 1) k1 y (n )
EN
x(n) a N (k ) y (n k )
k 1
k1 g 0 (n 1)
For N=2, then y ( n) x( n) a2 (1) y(n 1) a2 (2) y(n 2)
CITSTUDENTS.IN
TS .IN
1 AN ( Z ) k1
y (n 1)
a1 (1)
Page 153

This output can be obtained from a two-stage lattice filter as shown in below fig
10EC52
f 2 (n) x( n) f 1 (n) f 2 (n) k 2 g1 (n 1) g 2 ( n) k 2 f 1 (n) g1 (n 1) f 0 (n ) f1 ( n) k1 g 0 (n 1)

g1 (n ) k1 f 0 (n ) g 0 (n 1)
y( n)
f 0 ( n)
g 0 (n)
Similarly We observe
g 2 ( n)
TU D
x ( n) k m g m 1 (n 1) g m 1 (n 1) f m (n )
f 2 ( n) f 2 (n ) x( n) x( n)
k 2 g1 (n 1) k1 g 0 (n 1) k 2 k1 f 0 (n 1) g 0 (n 2) k1 g 0 (n 1) k 2 k1 y(n 1) y(n 2) k1 y(n 1) k1 (1 k 2 ) y(n 1) k 2 y(n 2) y(n 2) k2
k 2 y( n) k1 (1 k 2 ) y(n 1)
CI
f N (n) f m 1 (n) g m ( n) k m f m 1 (n)
TS
a2 (0) 1; a 2 (1) k1 (1 k 2 ); a 2 (2) N-stage IIR filter realized in lattice structure is,
EN
f 1 (n) k1 g 0 (n 1)
CITSTUDENTS.IN
TS .IN
m=N, N-1,---1 m=N, N-1,---1
Page 154

y (n ) f 0 (n) g 0 (n)
10EC52
8.6.1 Conversion from lattice structure to direct form:

a m ( m) am (k ) km ; a m 1 (k ) a m (0) 1 k)
a m (m)a m 1 (m
Conversion from direct form to lattice structure
a m 1 (k )
a m ( k ) a m ( m) a m ( m k ) 2 ( m) 1 am
8.6.2 Lattice Ladder Structure:
H (Z )
BM ( Z ) AN ( Z )
CI
y(n)
m 0 M
Where {C m} are called the ladder co-efficient and can be obtained using the recursive relation, Cm bm Ci ai (i
i m 1
TS
M
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 C m g m (n)
TU D
M
If,
bM (k ) Z
k 0 N
a N (k ) Z
k 1
m);
EN
m=M, M-1, .0
A general IIR filter containing both poles and zeros can be realized using an all pole lattice as the basic building block.
CITSTUDENTS.IN
TS .IN
a m 1 (0) 1
km
a m (m )
Page 155
10EC52
k 3 a3 (3) Using the equation
for m=3, k=1
TS
a 2 (1) a 2 (2)
5 8 1 13 3 24 1 9
a m 1( k )
TU D
1 3
Problem:4 Convert the following pole-zero IIR filter into a lattice ladder structure, 1 2Z 1 2Z 2 Z 3 H (Z ) 5 1 Z 1 8 Z 2 3 Z 3 1 13 24 Solution: 1 2Z 2 Z 3 Given b M ( Z ) 1 2 Z 1 And AN ( Z ) 1 13 Z 1 5 Z 2 3 Z 3 8 24 13 ; a3 (2) 5 ; a3 (3) 1 a3 (0) 1; a3 (1) 24 8 3
a m ( k ) a m ( m) a m ( m k ) 1 a 2 m( m)
13 24 1 5 3 8 1 2 3
CI
for m=3, & k=2 1 for m=2, & k=1 a1 (1) k1
a 3 (1) a 3 (3)a 3 (2) 2 1 a3 (3) k2
a3 (2) a3 (3)a 3 (1) 2 1 a3 (3)

45 13 72 8 9 1 2
a 2 (1) a 2 (2)a 2 (1) 2 1 a2 (2)

Page 156
CITSTUDENTS.IN
EN
. 1
3 8
TS .IN

3 8 1 3 2 8 1 2 2
10EC52
3 8
1 ( )
1
1 4
3 16 1 4
1 4
for lattice structure k1 For ladder structure

M
, k2
1 2
, k3
1 3
Cm M=3 C1 b1
2
bm b3
13 24
C1 .a1 (1 m)
i m 1
m=M, M-1,1,0 C 3 a3 (1)
C3
3
1; C 2
b2
b1
i 2
c1 a1 (i
m)
m=1
c2 a 2 (1)
3
c3 a3( 2)
5 8
3 1.4583 ( 8 )
0.8281
c0
b0
i 1
c1a1 (i m) b0 1 c1 a1 (1) c2 a 2 (2) c3 a3 (3)

1 1 08281 (4 ) 1.4583( 2 ) 1 3
CI
Problem 5
TS
TU D
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).
EN
02695
CITSTUDENTS.IN
TS .IN
Page 157
2 1.( ) 1.4583
10EC52
CI
TS
TU D
CITSTUDENTS.IN
Page 158
EN
TS .IN
10EC52
Question 6
CI
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.
TS
TU D
CITSTUDENTS.IN
Page 159
EN
TS .IN
10EC52
CI
TS
TU D
CITSTUDENTS.IN
Page 160
EN
TS .IN

Lattice realization :
10EC52
CI
TS
TU D
CITSTUDENTS.IN
Page 161
EN
TS .IN
Fig. 6

Ece V Digital Signal Processing (10ec52) Notes

Uploaded by

Document Informationclick to expand document informationNotes on Digital signal processing

Document Informationclick to expand document information

Copyright:

Available Formats

Ece V Digital Signal Processing (10ec52) Notes

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Ece V Digital Signal Processing (10ec52) Notes

Uploaded by

Copyright:

Available Formats

Digital Signal Processing

IA M ARKS EXAM HOURS EXAM M ARKS

NO. OF LECTURE HRS/WEEK : 04 TOTAL NO. OF LECTURE HRS. : 52

FOR EFFICIENT COMPUTATION OF THE DFT (FFT ALGORITHMS).

UNIT - 4 RADIX-2 FFT

ALGORITHM FOR THE COMPUTATION OF

TIME AND DECIMATION-IN-FREQUENCY

IDFTDECIMATION-INALGORITHMS. GOERTZEL ALGORITHM, AND CHIRP -Z 6 HOURS

Digital Signal Processing

FILTERS FROM ANALOG FILTERS

IMPULSE INVARIANCE METHOD. DERIVATIVE

CHEBYSHEV) APPROXIMATION OF TRANSFORMATION) METHOD,

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

3. DIGITAL SIGNAL PROCESSING, LEE TAN: ELSIVIER PUBLICATIONS, 2007

Digital Signal Processing

6.3 6.4 VII 7.1 7.2 7.3 7.4 7.5 7.6

Digital Signal Processing

Digital Signal Processing

DISCRETE FOURIER TRANSFORMS (DFT) CONTENTS:1. FREQUENCY DOMAIN SAMPLING 2.

3. DFT AS A LINEAR TRANSFORMATION

4. DFT RELATIONSHIP WITH OTHER TRANSFORMS.

DFT RELATIONSHIP WITH FOURIER SERIES DFT RELATIONSHIP WITH Z- TRANSFORMS

Digital Signal Processing

Discrete Fourier Transform

FREQUENCY DOMAIN SAMPLING AND RECONSTRUCTION OF DISCRETE TIME SIGNALS :

Digital Signal Processing

Fig 1.1 Frequency Domain Sampling

If we evaluate equation (1) at w X 2 k N xne

Digital Signal Processing

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

Let us have a look at that:

Digital Signal Processing

n 0 L xp[n] N>=L No aliasing

Discrete Fourier Transform:

Digital Signal Processing

Matrix relation of DFT

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

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

1.5 Relationship of Fourier Transforms with other transforms

Digital Signal Processing

we obtain discrete time sequence

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

Digital Signal Processing

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

Digital Signal Processing Recommended Questions with solutions

Digital Signal Processing

Digital Signal Processing

Digital Signal Processing

Digital Signal Processing

OVERLAP-SAVE AND OVERLAP-ADD METHOD.

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

Digital Signal Processing

Unit 2 Properties of DFT

Digital Signal Processing

2.1.3 Circular shift:

2.1.4 Time shift:

2.1.5 Frequency shift