Digital Signal Processing

( DSP )
Digital Signal Processing
COURSE FILE
(A.Y.2021-2022)
(Academic Regulations -2018)

DEPARTMENT OF Electronics and Communication Engineering

Courses file contents:

1. Cover Page
2. Vision of the Institute
3. Mission of the Institute
4. Vision of the Department
5. Mission of the Department
6. PEOs, POs and PSOs
7. Syllabus copy (scanned copy from the syllabus book)
8. Course objectives and Outcomes
9. Brief note on the course & how it fits in to the curriculum
10. Prerequisites, if any.
11. Instructional Learning Outcomes
12. Course mapping with PEOs, POs and PSOs.
13. Lecture plan with methodology being used/adopted.
14. Assignment questions (Unit wise)
15. Tutorial problems (Unit wise)
16. a) Unit wise short and long answer question bank
b) Unit wise Quiz Questions
17. Detailed notes##
18. Additional topics, if any.
19. Known gaps, if any.
20. Discussion topics, if any.
21. University/Autonomous Question papers of previous years.
22. References, Journals, websites and E-links, if required.
23. Quality Control Sheets. (to be submitted at the end of the semester)
a. Course end survey
b. Feedback on Teaching Learning Process(TLP)
c. CO- attainment
24. Student List (can be submitted later)
25. Group-Wise students list for discussion topics (can be submitted later)

DEPARTMENT OF Electronics and Communication Engineering

GEETHANJALICOLLEGE OF ENGINEERING AND TECHNOLOGY
DEPARTMENT OF Electronics and Communication Engineering
(Name of the Subject ) : Digital Signal Processing
(GCET CODE –18EC3201) Programme : UG

Branch: ECE
Year: III ECE Document No. GCET/ECE/----/-----
Semester: II No. of pages :232

Classification status (Unrestricted / Restricted ) : Unrestricted

Distribution List : Dept. Library, Dept Office, Concerned Faculty

Prepared by Updated by

1) Name : Mr.A.Srinivas 1) Name : Mr.A.Srinivas

2) Sign : 2) Sign :

3) Design : Asst. Professor 3) Design : Asst. Professor

4) Date : 20/12/2012 4) Date : 4/1/2022

Verified by : * For Q.C Only.

1) Name : Prof.P.Srihari 1) Name :

2) Sign : 2) Sign :

3) Design :Professor 3) Design :

4) Date : 4) Date :

Approved by : (HOD ) 1) Name : Dr. S.Suryanarayana

2) Sign :

3) Date :

DEPARTMENT OF Electronics and Communication Engineering

2. VISION OF THE INSTITUTE
“Geethanjali visualizes dissemination of knowledge and skills to students, who would
eventually contribute to the well being of the people of the nation and global community.”

3. MISSION OF THE INSTITUTE

 To Impart adequate fundamental knowledge in all basic sciences and engineering,
technical and inter-personal skills to students.
 To bring out creativity in students that would promote innovation, research and
entrepreneurship.
 To preserve and promote cultural heritage, humanistic and spiritual values
promoting peace and harmony in society.

4. VISION OF THE DEPARTMENT

To impart quality technical education in Electronics and Communication Engineering
emphasizing analysis, design/synthesis and evaluation of hardware/embedded software using
various Electronic Design Automation (EDA) tools with accent on creativity, innovation and
research thereby producing competent engineers who can meet global challenges with societal
commitment.

5. MISSION OF THE DEPARTMENT

i. To impart quality education in fundamentals of basic sciences, mathematics, electronics and

communication engineering through innovative teaching-learning processes.
ii. To facilitate Graduates define, design, and solve engineering problems in the field of
Electronics and Communication Engineering using various Electronic Design Automation
(EDA) tools.
iii. To encourage research culture among faculty and students thereby facilitating them to be
creative and innovative through constant interaction with R & D organizations and Industry.
iv. To inculcate teamwork, imbibe leadership qualities, professional ethics and social
responsibilities in students and faculty

6. PEOs, POs and PSOs

DEPARTMENT OF Electronics and Communication Engineering

6.1 Program Educational Objectives of B. Tech (ECE) Program:

I. To prepare students with excellent comprehension of basic sciences, mathematics

and engineering subjects facilitating them to gain employment or pursue
postgraduate studies with an appreciation for lifelong learning.

II. To train students with problem solving capabilities such as analysis and design with
adequate practical skills wherein they demonstrate creativity and innovation that
would enable them to develop state of the art equipment and technologies of
multidisciplinary nature for societal development.

III. To inculcate positive attitude, professional ethics, effective communication and

interpersonal skills which would facilitate them to succeed in the chosen profession
exhibiting creativity and innovation through research and development both as team
member and as well as leader.

6.2 Program Outcomes (POs)

1. Engineering knowledge: Apply the knowledge of mathematics, science, engineering
fundamentals, and an engineering specialization to the solution of complex engineering
problems.
2. Problem analysis: Identify, formulate, research literature, and analyze complex engineering
problems reaching substantiated conclusions using first principles of mathematics, natural
sciences, and engineering sciences.
3. Design/development of solutions: Design solutions for complex engineering problems and
design system components or processes that meet the specified needs with appropriate
consideration for the public health and safety, and the cultural, societal, and environmental
considerations.
4. Conduct investigations of complex problems: The problems
 That cannot be solved by straightforward application of knowledge, theories and
techniques applicable to the engineering discipline.
 That may not have a unique solution. For example, a design problem can be solved in
many ways and lead to multiple possible solutions.
 That requires consideration of appropriate constraints/requirements not explicitly given
in the problem statement. (like: cost, power requirement, durability, product life, etc.).
 Which need to be defined (modeled) within appropriate mathematical framework?
 That often requires use of modern computational concepts and tools.
5. Modern tool usage: Create, select, and apply appropriate techniques, resources, and
modern engineering and IT tools including prediction and modeling to complex engineering
activities with an understanding of the limitations.

DEPARTMENT OF Electronics and Communication Engineering

 6. The engineer and society: Apply reasoning informed by the contextual knowledge to
assess societal, health, safety, legal and cultural issues and the consequent responsibilities
relevant to the professional engineering practice.
7. Environment and sustainability: Understand the impact of the professional engineering
solutions in societal and environmental contexts, and demonstrate the knowledge of, and
need for sustainable development.
8. Ethics: Apply ethical principles and commit to professional ethics and responsibilities and
norms of the engineering practice.
9. Individual and team work: Function effectively as an individual, and as a member or
leader in diverse teams, and in multidisciplinary settings.
10. Communication: Communicate effectively on complex engineering activities with the
engineering community and with society at large, such as, being able to comprehend and
write effective reports and design documentation, make effective presentations, and give
and receive clear instructions.
11. Project management and finance: Demonstrate knowledge and understanding of the
engineering and management principles and apply these to one’s own work, as a member
and leader in a team, to manage projects and in multidisciplinary environments.
12. Life-long learning: Recognize the need for, and have the preparation and ability to engage
in independent and life-long learning in the broadest context of technological change.

6.3 PROGRAM SPECIFIC OUTCOMES (PSOs):

1. An ability to design an Electronics and Communication Engineering system, component, or
process and conduct experiments, analyze, interpret data and prepare a report with
conclusions to meet desired needs within the realistic constraints such as economic,
environmental, social, political, ethical, health and safety, manufacturability and
sustainability.
2. An ability to use modern Electronic Design Automation (EDA) tools, software and electronic
equipment to analyze, synthesize and evaluate Electronics and Communication Engineering
systems for multidisciplinary tasks.

7. Syllabus:

DEPARTMENT OF Electronics and Communication Engineering

DEPARTMENT OF Electronics and Communication Engineering
8. Course Objectives and Outcomes:
Objectives: Develop ability to

DEPARTMENT OF Electronics and Communication Engineering

 1. Understand fundamental concepts involved in the analysis and processing of discrete signals.
2. Distinguish between various discrete -time signals and Systems.
3. Understand frequency domain analysis of discrete signals and systems using DTFT, DFT and
FFT tools.
4. Understand the design of Infinite Impulse Response (IIR) and Finite Impulse Response (FIR)
filters for a given specifications.
5. Understand Multi-rate signal processing Techniques and finite word length effects.
Course Outcomes:
CO1. Perform analysis on discrete time signals and systems in the frequency domain using DFS,
DTFT and Z transform
CO2. Compute the DFT of a given discrete time sequence and plot the spectrum respectively.
CO3. Compute radix-2 FFT for a given sequence.
CO4. Design IIR and FIR filters for given specifications
CO5. Convert from one sampling rate to another. Analyze finite word length effects in digital
filters.

9. Brief note on the course & how it fits into the curriculum
Digital Signal Processing (DSP) is concerned with the representation, transformation and
manipulation of signals on a computer. After half a century advances, DSP has become an
important field, and has penetrated a wide range of application systems, such as consumer
electronics, digital communications, medical imaging and so on. With the dramatic increase of
the processing capability of signal processing microprocessors, it is the expectation that the
importance and role of DSP is to accelerate and expand.
Discrete-Time Signal Processing is a general term including DSP as a special case. This course
will introduce the basic concepts and techniques for processing discrete-time signal. By the end
of this course, the students should be able to understand the most important principles in DSP.
The course emphasizes understanding and implementations of theoretical concepts, methods and
algorithms.
The following answers will give an idea of how it fits in to the curriculum.

i. What role does this course play within the Program?

This course strengthens analysis and design of digital filter capabilities of the students.

ii. How is the course unique or different from other courses of the Program?

It makes the students to learn the concepts of Digital Filters and Multirate Systems.

iii. What essential knowledge or skills should they gain from this experience?

DEPARTMENT OF Electronics and Communication Engineering

 Students acquire design, analysis and simulation capabilities with this course

iv. What knowledge or skills from this course will students need to have mastered to
perform
well in future classes or later (Higher Education / Jobs)?

Thorough knowledge on fundamentals as well as analytical and design skills

v. Why is this course important for students to take?

In order to design, simulate and develop the LSI(Linear Shift Invariant)system, this
course is essential.

vi. What is/are the prerequisite(s) for this course?

1. Theory of Signals and Systems.

vii. When students complete this course, what do they need know or be able to do?

Able to design, analyze, simulate, compare and evaluate the digital circuits.

viii. Is there specific knowledge that the students will need to know in the future?

In addition to the concepts of z-Transforms, DFT, FFT and the concept of finite word
length effects are needed for future courses.

ix. Are there certain practical or professional skills that students will need to apply in the
future?

YES

x. Five years from now, what do you hope students will remember from this course?

The concepts of different types of Digital filters (FIR and IIR) and importance of
Sampling rate conversion.

xi. What is it about this course that makes it unique or special?

After completion of this course, the student can able design any Digital Filters as per the
specifications.

xii. Why does the program offer this course?

DEPARTMENT OF Electronics and Communication Engineering

 This course is a prerequisite for DSPA and DIP. And also it cultivates analytical and
design skills.

xiii. Why can’t this course be “covered” as a sub-section of another course?

It is not possible as it covers widely about different types of digital filters.(Frequency

Approach)

xiv. What unique contributions to students’ learning experience does this course make?

It helps in executing mini and major projects having digital circuits and DSP processors.

xv. What is the value of taking this course? How exactly does it enrich the program? The
“Course Purpose” describes how the course fits into the student's educational experience
and curriculum in the program and how it helps in his/her professional career.

This course plays a vital role in design and development of Electronic and digital
Communication system useful to the society and this course also helps for the student’s
professional career growth.
10. Prerequisite, if any

18PH1202- Signals and Systems

11. Instructional learning outcomes

UNIT-I :Introduction to Digital Signal Processing

1) Students can understand the concept of discrete time signals & sequences.
2) Analyze and implement digital signal processing systems in time domain.
3) They can solve linear constant coefficient difference equations.
4) They can understand Frequency domain representation of discrete time signals and
systems.
5) They can understand the practical purpose of stability and causality.
6) To determine stability, causality for a given impulse response.
7) Understand how analog signals are represented by their discrete-time samples, and in
what ways digital filtering is equivalent to analog filtering.
8) The basics of Z-transforms and its applications are studied.
9) Ability to understand discrete time domain and frequency domain representation of
signals and systems using DFS and DTFT.
10) Calculate the response of applying a given input signal to a system described by a linear
constant coefficient differential equation.

UNIT-II: Discrete Fourier Transform (DFT) and FFT

1) Ability to understand discrete time domain and frequency domain representation of

DEPARTMENT OF Electronics and Communication Engineering

 signals and systems.
2) Compute convolution and the discrete Fourier transform (DFT) of discrete-time signals.
3) Analyze and implement digital systems using the DFT.
4) Ability to understand Discrete Fourier Series and Transforms and comparison with other
transforms like Z transforms.
5) Ability to represent discrete-time signals in the frequency domain.
6) Calculate exponential Fourier series coefficients using properties of Fourier series
7) Graphically portray the magnitude and phase of the Fourier series coefficients versus ω.
Fast Fourier Transforms
8) Ability to develop Fast Fourier Transform algorithms for faster realization of signals and
systems.
9) Analyze and implement digital systems using the FFT.
10) Describe how and why Fourier Transforms and Fourier series are related.

UNIT-III :Design of DIGITAL FILTERS (IIR)- Structures of IIR systems

1) The student is able to solve basic digital signal processing algorithms.
2) Assess signal acquisition, processing, and reconstruction.
3) Digital filters are realized using difference equations.
4) Ability to understand the concepts of Digital filters IIR like Chebychev, Butter worth
filters.
5) The student is able to solve the impulse and frequency response of IIR filters given as
difference equations, transfer functions, or realization diagrams, and can present analyses
of the aliasing and imaging effects based on the responses of the filters.
6) Learn the basic forms of IIR filters, and how to design filters with desired frequency
responses.
UNIT-IV :Design of DIGITAL FILTERS (FIR) – Structure of FIR Systems

1) Ability to understand the characteristics of linear-phase finite impulse response (FIR)

filters
2) Ability to understand Digital Filters with special emphasis on realization of FIR and IIR
filters.
3) Ability to design linear-phase FIR filters according to predefined specifications using the
window and frequency sampling methods
4) Ability to understand the concepts of Digital FIR filters.
5) The student is able to specify and design respective frequency selective FIR and IIR
filters using the most common methods.
6) The student is able to solve for the impulse and frequency responses of FIR and IIR
filters given as difference equations, transfer functions, or realization diagrams, and can
present analyses of the aliasing and imaging effects based on the responses of the filters.
7) Measure the effectiveness of FIR filters.

UNIT-V: Introduction to Multi rate Digital Signal Processing

1) Ability to understand the concepts of sampling rate conversions, Decimation and
Interpolation as part of Signal Processing techniques.

DEPARTMENT OF Electronics and Communication Engineering

 2) Able to explain how the multi rate implementation of ADC and DAC converters works.
3) Able to describe basic sampling rate conversion algorithms.
4) Able to draw and describe different kinds of interpolator and decimator.
5) Able to analyze how the interpolated FIR filter works.
6) Able to do sampling rate conversion.
Finite word length effects in Digital filters

1) Ability to analyze the concepts of Signal scaling.

2) Ability to analyze the concepts Quantization errors in FFT algorithms.
3) Ability to analyze Number representation, Floating point numbers, quantization
Noise, Overflow limit cycle oscillations, Signal scaling, Finite Word Length effects
in FIR, Quantization errors in FFT algorithms.
4) The student is able to explain the impact of finite word length in filter design.
5) The student is able to design simple digital filters.
6) Concept of quantization Noise is analyzed and errors are identified.

2. Course Outcomes mapping with POs, PSOs and PEOs

Relationship of the course to program outcomes:

Course Outcomes- Program Outcomes Mapping Matrix

Note: Enter Correlation Levels as mentioned below

1: Slight (Low) 2: Moderate (Medium) 3: Substantial (High) if there is no correlation put “__”

CO# PO1 PO2 PO3 PO PO PO PO PO PO PO1 PO1 PO1 PSO

4 5 6 7 8 9 0 1 2 PSO1 2
18EC3201. 3 3 2 - - - - - - - 1 1
1 3 1
18EC3201. 3 3 2 - - - - - - - 1 1
2 3 2
18EC3201. 3 3 2 - - - - - - - 1 1
3 3 1
18EC3201. 3 3 3 - - - - - - - 1 1
4 3 1
18EC3201. 3 3 2 - - - - - - - 1 1
5 3 2
18EC3201. 3 3 2 - - - - - - - 1 1
6 3 2
3 3 2.166 - - - - - - - 1 1 3 1.5

Course Outcomes- Program Specific Outcomes Mapping Matrix

DEPARTMENT OF Electronics and Communication Engineering

 S.No Course Code Course Semeste PEO PEO PEO 3
componen r 1 2
t
1 Core 18EC3201 DSP III-II √ √ -

Note: Enter Correlation Levels as mentioned below

1: Slight (Low) 2: Moderate (Medium) 3: Substantial (High) if there is no correlation put “__”

Relationship of the course to the Program Educational Objectives:

I. To prepare students with excellent comprehension of basic sciences, mathematics
and engineering subjects facilitating them to gain employment or pursue
postgraduate studies with an appreciation for lifelong learning.

II. To train students with problem solving capabilities such as analysis and design
with adequate practical skills wherein they demonstrate creativity and innovation
that would enable them to develop state of the art equipment and technologies of 
multidisciplinary nature for societal development.

III. To inculcate positive attitude, professional ethics, effective communication and

interpersonal skills which would facilitate them to succeed in the chosen profession
exhibiting creativity and innovation through research and development both as team
member and as well as leader

13. Lecture schedule with methodology being used/adopted.

Section A: Name of the Faculty: Mr.K. Victor
Designation: Assistant Professor
Section B: Name of the Faculty: Mr. U Nageswar Rao
Designation: Assistant Professor
Section C: Name of the Faculty: Mr. M. Jagan Mohan Rao
Designation: Assistant Professor
Section D: Name of the Faculty: Mr. A. Srinivas
Designation: Assistant Professor
Section E: Name of the Faculty: Mr. R. Santosh
Designation: Assistant Professor

13.1 COURSE SCHEDULE (Expected)

DEPARTMENT OF Electronics and Communication Engineering

Academic Year: 2021-2022 Year: III B. Tech Semester: II

Branch: ECE Section: A,B,C,D,E

Subject Name : Digital signal Processing Subject Code: 18EC3201

Faculty : Mr. A. Srinivas Total No of Classes: 65 (Estimated)

Contact No. : 9603296965 Commencement of Instruction:

E-mail :
Completion of Instruction:

Duration
S.N (Dates) Total
Unit No Title
o No. Periods
From To
1. I Introduction to Digital Signal Processing 15
2. II Discrete Fourier Transform 14
3. III Design of Digital IIR Filter 13
4. IV Design of FIR Digital Filters 12
5. V Introduuction to multirate signal Processing 11
Total 65

13.2 LESSON PLAN

DEPARTMENT OF Electronics and Communication Engineering

 Academic Year: 2020-21 Year: III B. Tech Semester:II

Branch: ECE Section: A,B,C,D,E

Subject Name : Digital Signal Processing Subject Code : 18EC3201
Faculty : Mr. A. Srinivas
Total No of Classes: 65 (Estimated)
Contact No. 9948876998 Commencement of Instruction: 17-01-2022
E-mail : [email protected]
Completion of Instruction :

Time Table of the Subject:

Wednesda
Day Monday Tuesday Thursday Friday Saturday
y
Period
2 1 1 1
s

S.N Total Topic to be covered in Regular/ Teaching Remarks

o number of One Additional/ Aids used
Unit Date lecture Missing LCD
No. periods /OHP/BB
1 Introduction to Digital Regular BB
Signal
18/01/2022 processing
2 18/01/2022 Review of Discrete time Regular BB
signals and
sequences
3 19/01/2022 Analog to Digital Regular BB
Conversion Process

4 19/01/2022 Sampling of LP and BP Regular BB

Signals
difference equations
5 20/01/2022 Analysis of Discrete time Regular BB/OHP
Invariant Systems, Causal
Linear shift invariant
U
systems,
N
6 15 25/01/2022 representation of discrete Regular BB/OHP
I time signals and
systemsstability
7 T 27/01/2022 causality, problemsLinear Regular
constant coefficient
- difference Equations BB
8 27/01/2022 solution of difference Regular BB

DEPARTMENT OF Electronics and Communication Engineering

 I equations
9 28/01/2022 Regular BB

10 01/02/2022 Introduction to Regular BB

DFS, Properties of
DFS

11 01/02/2022 Discrete Time Regular BB

Fourier Transform
(DTFT)

12 02/02/2022 Relation between Z- Regular BB

Transform and DFS

13 02/02/2022 Problems Regular BB

14 03/02/2022 Assignment-1 Regular BB

15 08/02/2022 Tutorial class-1

16 09/02/2022 Discrete Fourier Regular BB

Transform: Its Properties
and Applications
17 09/02/2022 Relation of DFT to other Regular
(DTFS and Z-Transforms)
BB
transforms.
and DFS.
18 10/02/2022 Inverse DFT(IDFT) Regular BB

19 10/02/2022 Linear convolution of Regular BB/OHP/

sequences using DFT, LCD
20 15/02/2022 Computation of DFT and Regular BB/OHP/
IDFT LCD
21 U 15/02/2022 Fast Fourier transformer: Regular BB/OHP/
Efficient computation of LCD
N DFT,
22 14 16/02/2022 FFT Algorithms, Direct Regular BB/OHP/
I computation of DFT. LCD
23 16/02/2022 Radix-2, decimation in Regular BB/OHP/
T
time decimation in LCD
frequency FFT

DEPARTMENT OF Electronics and Communication Engineering

 - algorithms
24 17/02/2022 Divide and conquer Regular BB/OHP/
II approach to computation of LCD
DFT(Radix-N algorithms)
25 22/02/2022 Circular convolution of Regular BB/OHP/
sequences using DFT LCD
26 22/02/2022 Overlap and Overlap save Regular BB
method
27 23/02/2022 PRACTICE PROBLEMS Regular BB

28 23/02/2022 Assignment-2 Regular BB

29 24/02/2022 Tutorial class-2 Regular BB

30 Design of IIR DIGITAL Regular BB/OHP/

FILTERS: Direct Form I LCD
25/02/2022
and II, and
31 01/03/2022 Cascade form Regular BB/OHP

32 01/03/2022 Parallel form structures, Regular BB/OHP

33 02/03/2022 Design of IIR Filters from Regular BB/OHP

analog filters:
34 02/03/2022 Characteristics of Regular BB/OHP
commonly used analog
U filters
35 03/03/2022 Analog filter Regular BB/OHP
N approximations-
13 Butterworth and
36 I 08/03/2022 Chebyshev, Regular BB/OHP

37 T 08/03/2022 IIR filter design by Impulse Regular BB/OHP

invariance
-
38 09/03/2022 Bilinear Transformation Regular BB
III
method,
39 09/03/2022 Frequency transformations. Regular OHP

40 10/03/2022 Assignment-3 Regular BB/OHP

41 15/03/2022 Problems practice Regular BB/OHP

42 15/03/2022 Tutorial-3 Regular BB/OHP

43 16/03/2022 Untroduction Regular BB/OHP

44 16/03/2022 Design of digital filters Regular OHP

(FIR):
U Structure of FIR Systems:

DEPARTMENT OF Electronics and Communication Engineering

 N Direct form,
45 17/03/2022 Cascade realization,and Regular BB/OHP
I Linear phase realization;
45 12 22/03/2022 Characterstics of linear Regular BB/OHP
T
phase FIR filter and its
frequency response;
47 - 23/03/2022 Regular BB/OHP
Comparison of IIR and FIR
IV filters;
48 24/03/2022 Design of linear phase FIR Regular BB/OHP
filters using windows
method
49 25/03/2022 Rectangular Regular BB/OHP
window,Hanning window,
Hamming window,
50 25/03/2022 Bartlett window,Kaiser Regular BB/OHP
window
51 05/04/2022 frequency-sampling Regular BB/OHP
method.
52 05/04/2022 Problem from FIR Filters Regular BB/OHP

53 06/04/2022 Assignment-4 Regular BB/OHP

54 06/04/2022 Tutorial-4 Regular BB/OHP

57 07/04/2022 Introduction to Multirate Regular BB/OHP

Digital signal Processing
58 07/04/2022 Decimation by a factor D Regular BB/OHP

59 08/04/2022 Interpolation by a factor I Regular BB/OHP

53 U 12/04/2022 Sampling rate conversion Regular BB

by a rational factor I/D
54 N 13/04/2021 Multistage implementation Regular BB/OHP
of sampling rate conversion
55 I 14/04/2021 Aplication of multirate Regular BB/OHP
T signal processing.
56 15/04/2021 Introduction to Finite word Regular BB/OHP
- length effects in fixed point
DSP System.
57 V 15/04/2021 Assignment-5 Regular BB

58 19/04/2021 Tutorial-5 Regular BB

59 11 20/04/2021 Revision Regular BB

60 21/04/2021 Revision Regular BB

DEPARTMENT OF Electronics and Communication Engineering

61 22/04/2021 Revision Regular BB

62 26/04/2021 Solution of University Additional BB

Question Paper and
Revision

Books / Material

Text Books (TB)

TB1: Digital signal Processing: Principles, Algorithms and Applications-John G.Proakis,
D.G.Manolakis, 4thEdition,Pearson/PHI,2009.
TB2: Digital Signal Processing, S K Mitra,3/e, TMH, 2006. References 1. Discrete time signal
Processing- 2.
Reference Books (R)
R1: A.V.Oppenheim and R.W.Schaffer,PHI,2009.
R2: Digital signal Processing-A Practical Approach-Emmanuel C.Ifeacher, Barrie.W.Jervis, 2nd
Edition, Pearson Education,2009.
R3: Fundamentals of Digital signal Processing using MAT Lab-Robert J.Schilling, Sandra
L.Harris,Thomson,2007.

GUIDELINES

Distribution of periods:
No. of classes required to cover GCET syllabus : 60
No of classes required to solve University papers :5
Total classes required : 65

14.Assignment Questions

DEPARTMENT OF Electronics and Communication Engineering

 Unit 1
1 ) For each of the following systems, determine whether the system is stable and time invariant
or not :

i. T(x[n]) =

ii. T(x[n]) = ex(n)

ii. 2. y(n) = Cos(w0n)x(n)

2) Compute the convolution of the following signals

x(n)=u(n+1)-u(n-4)-u(n-5); h(n)=[u(n+2)- u(n-3)]

3) Implement direct form II realization of the following system:

4) Determine the zero -input response of the system described by the second order difference
equation x(n)-3y(n-1)-4y(n-2)=0

5) Compute the convolution of the following signals x(n)=a n u(n), h(n)=b n u(n) when
a ≠b and a=b

6 a) Determine the particular solution of the difference equation

y(n)=5/6 y(n-1)-1/6 y(n-2)+x(n); When forced by x(n)=2nu(n)

|x(n) |
b) Determine whether the given system y(n) = ne has time invariance and stability using
appropriate test

7) Obtain the Direct form-I realization for the system described by the difference equation
y(n)=0.5y(n-1)-0.25y(n-2)+x(n)+0.4x(n-1)

8. a) Prove the Initial value theorem and final value theorem of Z Transforms.

b) List and prove any two properties of the DFS.

9. Determine the Fourier series spectrum of signals.

i) x(n) = cos πn/3 ii) x(n) = cos

iii) x(n) is periodic with period N=4 and x(n) = [1 1 0 0 ]

DEPARTMENT OF Electronics and Communication Engineering

 10..a)Define an LTI System and show that the output of an LTI system is given by the convolution of

Input sequence and impulse response.

b) Prove that the system defined by the following difference

onequati is an LTI system y(n) = x(n+1)-3x(n)+x(n-1) ; n≥0.

11.aWrite short notes on classification of systems.

b) Derive BIBO stability criteria to achieve stability of a system.

12 .a)Discuss various discrete time sequences.

b) Give the Basic block diagram of Digital Signal Processor.

13.a)Define Linearity, Time Invariant, Stability and Causality.

b)The discrete time system is represented by the following difference equations in which x(n) is

input and y(n) is output. Y(n) = 3y2(n-1)-nx(n)+4x(n-1)-2x(n-1).

14.(a) Discuss impulse invariance method of deriving IIR digital filter from corre-sponding analog

filter.

(b) Use the Bilinear transformation to convert the analog filter with system func- tion

H (S) = S + 0.1/(S + 0.1)2 + 9 into a digital IIR filters. Select T = 0.1 and compare

the location of the zeros in H(Z) with the locations of the zeros obtained by applying

the impulse invariance method in the conversion of H(S).

15. a)Define Z-Transform and List out its properties.

b)Discuss Direct form, Cascade and Linear phase realization structures of FIR filters.

16. a)Discuss transposed form structures with an example.

b)Discuss Direct form, Cascade realization structures of FIR filters.

17. Discuss and draw various IIR realization structures like Direct form – I, Direct form-II,

DEPARTMENT OF Electronics and Communication Engineering

 Parallel and cascade forms for the difference equation given by

y(n) = - 3/8 Y(n-1) + 3/32 y(n-2) + 1/64 y(n-3) + x(n) + 3 x(n-1)

18. What are the various basic building blocks in realization of Digital Systems and hence discuss

transposed form realization structures.

(a) Implement the decimation in time FFT algorithm for N=16.

(b) In the above Question how many non - trivial multiplications are Required.

UNIT-02

1) List the properties of the DFT and prove the following properties:

i) Circular frequency shift ii) Multiplication of two sequences.

iii) Convolution Property

2) Find the 4-point DFT of the following sequences

a) x(n) = {1 -1 0 0} b) x(n) = {1 1 -2 -2}

3) Compute IDFT of sequence X(K)={7 7-j70 7 –j 0.7 -j0.707 0.707+j0.707 1 707+j.707}

4) Compute 8-point for the following sequences using DIT- FFT algorithm.

a) x1(n)= 1 for -3≤n≤3 b) x2(n) = 1 for 0≤n≤6

0 otherwise 0 otherwise

5) Find the IDFT of sequence X(k)={4 1-j2.414 0 1-j.414 0 1+j.414 0 1+j2.414} using DIF algorithm.

6) Draw the signal flow graph for 16-point DFT using a) DIT algorithm b) DIF algorithm.

7) Compute 8-point DFT of the sequence x(n)={1/2,1/2,1/2,1/2} using the in- place radix-2

Decimation – In-Time ( DIT-FFT) algorithm.

8) a).Define DFS. State any Four properties of DFS.

b)Find the IDFT of the given sequence x(K) = {2, 2-3j, 2+3j, -2}.

9) Define Convolution. Compare Linear and Circular Convolution techniques. b) Find the

Linear convolution of the given two sequences x(n)={1,2} and

DEPARTMENT OF Electronics and Communication Engineering

 h(n) ={1,2,3} using DFT and IDFT.

10) a)Design a high pass filter using hamming window with a cut-off frequency 1.2 radians/second and

N=9

(b) Compare FIR and IIR filters.

11 ) a)Define DFS. State any Four properties of DFS.

b) Find the IDFT of the given sequence x(K) = {2, 2-3j, 2+3j, -2}.

12) a) Define DFT and IDFT. State any Four properties of DFT.

b)Find 8-Point DFT of the given time domain sequence x(n) = {1, 2, 3, 4}.

13) a)Find IFFT of the given X(K) = { 1,2,3,4,4,3,2,1}using DIF algorithm

b)Bring out the relationship between DFT and Z-transform.

14) a)Develop DIT-FFT algorithm and draw signal flow graphs for decomposing the

DFT for N=6 by considering the factors for N = 6 = 2.3.

b)Bring out the relationship between DFT and Z-transform.

15). a) For each of the following systems, determine whether or not the system is

i. stable ii. Causal iii. linear iv. shift-invariant.

A. T [x(n)] = x(n − n0 )

B. T [x(n)] = ex(n)

C. T[x(n)] = a x(n) + b. Justify your answer.

(b) A system is described by the difference equation y(n)-y(n-1)-y(n-2) = x(n-1). Assuming that

the system is initially relaxed, determine its unit sample response h(n).

16) .a)Derive the expressions for computing the FFT using DIT algorithm and hence draw the

standard butterfly structure.

DEPARTMENT OF Electronics and Communication Engineering

 b)Compare the computational complexity of FFT and DFT.

17 .a)Find X(K) of the given sequence x(n) = { 1,2,3,4,4,3,2,1}using DIT- FFTalgorithm.

b)Compare the computational complexity of FFT and DFT.

UNIT-03

1).Obtain the Direct form-I, Direct form-II,Cascade Form and Parallel Form realization for the

system described by the difference equation y(n)=0.5y(n-1)-0.25y(n-2)+x(n)+0.4x(n-1)

2) Determine the order of the filter for the given specifications

αp =1dB αs = 30 dB ; Ωp = 200 rad/sec ; Ωs = 200 rad/sec.

3) List out the steps involved in designing an analog Butterworth low pass filter.

4) Compute the poles of an analog Chebyshev filter transfer function that

Satisfies the constraints

0.707 ≤ |H(jΩ) | ≤ 1 ; 0≤Ω≤2

|H(ejw) | ≤ 0.1 ; Ω≥4

And determine Ha (S) and hence obtain H(Z) using Bilinear Transformation method.

5). Compare Impulse Invariant Method with Bilinear Transformation method of

Designing IIR filters.

6) A 4-th order Butterworth filter has cut off frequency c  200  rad/sec.

a) Determine the zeros and poles of the transfer function.

b) What would be its pass band and stopband frequencies if we want

1dB ripple in the pass band and 40dB attenuation in the stopband?

c) If we apply a Bilinear Transformation with sampling frequency Fs  1 kHz,

Determine the zeros and poles in the z-plane.

7). Given the specification αp=3dB, αs=16dB, fp=1KHz, fs=2KHz. Determine the

DEPARTMENT OF Electronics and Communication Engineering

 Order of the filter Using chebyhev approximation. Find H(s).

8). An analog filter has a transfer function H(s)=2/(s+1)(s+2).Design a digital filter

Equivalent to this using impulse invariant method for T=1Sec.

9). Determine the order of low pass Butterworth filter that has a 3 dB at 500Hz and an

Attenuation of 40 dB at 1000 Hz.

10) Compare an analog filter with that of a digital filter.

UNIT-04

1) Determine the coefficients of a linear phase FIR filter of length M=15 has a asymmetric unit sample

response and a frequency response that satisfies the conditions

H(2πk/15) =1, k =0,1,2,3 and

=0 k = 4,5,6,7

2) Explain in detail the comparison of different windows on FIR filter design.

3) Explain how to design FIR filter using frequency sampling technique.

4) Design an ideal low pass filter with frequency response

Hd(ejw) = 1 for –π/2 ≤ w ≤ π/2

0 for π/2 ≤ | w | ≤ π

5).Find the value of h(n) for N=11, find H(Z) . Plot the magnitude response.

6) Design an ideal Band reject filter with a desired frequency response

Hd(ejw) = 1 for | w| ≤ π/3 & | w| ≥ 2π/3

0 otherwise

Find the value of h(n) for N = 11. Find H(Z), plot the magnitude response.
7) Design an ideal differentiate H(ejw) = j ω -π≤ ω≤π

Using a) rectangular window b)Hamming window with N=8.plot frequency response in both cases.

DEPARTMENT OF Electronics and Communication Engineering

8) using frequency sampling method design a bandpass filter with following specifications

Sampling frequency F=8000Hz Cut off frequency fc1=1000Hz fc2=3000Hz Determine the filter

coefficients for N=7.

9) a)Derive the conditions to achieve Linear Phase characteristics of FIR filters

b) Design an FIR Digital Low pass filter using Hanning window whose cut 1ff freq is

2 rad/s and length of window N=9.

10) .a)Compare FIR and IIR filters

b)Design an FIR Digital High pass filter using Hamming window whose cut off

freq is 1.2 rad/s and length of window N=9.

11 a)Compare various windowing functions

b)Design an FIR Digital Band pass filter using rectangular window whose

upper and lower cut off freq.’s are 1 & 2 rad/s and length of window N = 9.

12).a)Compare various windowing functions.

b)Design an FIR Digital Low pass filter using rectangular window whose cut off

freq is 2 rad/s and length of window N=9.

UNIT-05

1) Give the frequency domain analysis of interpolation by a factor of ‘ I ’ .

2) For the sequence x(n)={ 5,6,8,4, 2,1,3,12,10,7,11} find the output sequence y(z) which is

Down sampled version of x(n) by 2.

3.Derive the expression for decimation by factor D.

4. Derive the expression for interpolation by factor I.

5. Write notes on sampling rate conversion by a rational factor I/D.

6. a) Discuss the implementation of Polyphase filters for Interpolators with an example

DEPARTMENT OF Electronics and Communication Engineering

 b) Discuss the sampling rate conversion by a factor I/D with the help of a Neat block Diagram.

7. a) Define Interpolation and Decimation.

b) Discuss the sampling rate conversion by a factor I/D with the help of aNeat block Diagram.

8.a) Define Interpolation and Decimation. List out the advantages of Sampling

rate conversion.

b)Discuss the sampling rate conversion by a factor I with the help of a Neat

block Diagram.

9. a) Define Multirate systems and Sampling rate conversion

b) Discuss the process of n Decimation by a factor D and explain how the aliasing

effect can be eliminated

10 .(a) An LTI system is described by the equation

y(n)=x(n)+0.81x(n-1)-0.81x(n-2)-0.45y(n-2).

Determine the transfer function of the system. Sketch the poles and zeroes on the Z-

plane.

(b) Define stable and unstable systems. Test the condition for stability

of the first-order IIR filter governed by the equation y(n)=x(n)+bx(n-).

11 (a) Compute Discrete Fourier transform of the following finite length sequence

considered to be of length N.

i. x(n) = δ(n + n0 ) where 0 < n0 < N

ii. x(n) = an where 0 < a < 1.

DEPARTMENT OF Electronics and Communication Engineering

 (b) If x(n) denotes a finite length sequence of length N, show

that x((−n))N=x((N − n))N

15. Tutorial Problems

Tutorial -1

1) Determine the response y(n), n>0 of the system described by the second order
difference equation y(n) -5y(n-1) +6y(n-2) = x(n), for x(n) = n

2) Determine the cascade-form realization for each of the following system:

3) Determine the particular solution of the difference equation

y(n)=(5/6) y(n-1) - (1/6) y(n-2) + x(n) when the forcing function is

x(n)=2nu(n).

4). Obtain the i) Direct forms ii) cascade iii) parallel form realizations for the

Following systems y (n) = 3/4(n-1) – 1/8 y(n-2) + x(n) +1/3 x(n-1) .

5). Use the one-sided Z-transform to determine y(n) n ≥ 0 in the following cases.

(a) y(n)−1.5y(n−1) +0.5y(n−2) = 0; y(−1) = 1; y(−2) = 0

i) Compute the 10 first samples of its impulse response.

ii) Find the input-output relation.

iii) Understand the input x(n) = {1 1. . . .} and compute the first 10 samples of the output.

iv) Compute the first 10 samples of the output for the input given in part (c) by

Using convolution. is the system causal? Is it stable?

Tutorial -2
1. Calculate the 4 – point IDFT of X (K) = [1, -1, 2, -2] using DIT FFT algorithm. Compare

DEPARTMENT OF Electronics and Communication Engineering

 the number of calculations required to find DFT of a sequence using Radix -2 FFT
Algorithms and using DFT formula.
2. The DTFT of a real signal x(n) is X(F). How is the DTFT of the following?
Signals related to X(F). (a) y(n)=x(-n) (b) r(n)=x(n/4) (c) h(n) =jnx(n)
Find the IDFT of sequence X(k)={4 1-j2.414 0 1-j.414 0 1+j.414 0 1+j2.414} using DIF
Algorithm.
3.Draw the signal flow graph for 16-point DFT using a) DIT algorithm b) DIF algorithm.
Given x(n)=2n and N=8 find X(k) using DIT-FFT algorithm.
4. Given x(n)=2n and N=8 find X(k) using DIT-FFT algorithm.
5. Find the DFT of a sequence x(n)={0 1 2 3 4 5 6 7}.
Tutorial -3
1) Convert Analog filter with a transfer function H(s) = S+0.1/ ((S+0.1)2+9) into digital IIR f
filter equivalent using impulse invariant method .
2). Design a digital Chebyshev filter to meet the following specifications:
0.8 ≤ |H(ejw) | ; 0 ≤ w ≤ 0.2π
jw
|H(e ) | ≤ 0.2 ; 0 .6 ≤ w ≤ π
Using Bilinear Transformation method.
3).Design an ideal high pass filter with frequency response
Hd(ejw) = 1 for π/4 ≤ w ≤ 3π/4
0 other wise
Find the value of h(n) for N = 11. Find H(Z), plot the magnitude response.
4). given the specification αp=3dB, αs=16dB, fp=1KHz, fs=2KHz. Determine the order of the
Filter Using chebyshev approximation. find H(s).
5). Design an analog Butterworth filter that as a -2dB pass band attenuation at a
Frequency of 20rad/sec and at least -10dB stop band attenuation at 30rad/sec.
6). Determine the order and the poles of type 1 low pass chebyshev filter that has a 1 dB ripple
in the pass band and pass band frequency Ωp =1000π and a stop band of frequency of
2000π and an attenuation of 40dB or more.

Tutorial -4
1. Design an ideal low pass filter with frequency response
Hd(ejw) = 1 for –π/2 ≤ w ≤ π/2
0 for π/2 ≤ | w | ≤ π

DEPARTMENT OF Electronics and Communication Engineering

2. Find the value of h(n) for N=11, find H(Z) . Plot the magnitude response.
3. Design an ideal Hilbert transformer having frequency response H (jώ) = j
-π ≤ ω ≤ 0 -j 0 ≤ ω ≤ π for N=11 using rectangular window
4. using a rectangular window technique design a low pass filter with passband gain of unity,
Cutoff frequency of 100Hz and working at a sampling frequency of 5KHz.The length of the
Impulse response should be7.
5. Design a HPF of length 7 with cut off frequency of 2 rad/sec using Hamming window.
Plot the magnitude and phase response.

Tutorial -5
1. Understand the decimation process with a neat block diagram.
2. Consider a signal x(n)=sin(∏n)U(n). Obtain a signal with an interpolation factor of ‘2’.
3. Why multirate digital signal processing is needed?
4. Design a two state decimator for the following specifications. Decimation factor = 50
Pass band = 0<f<50 Transitive band = 50≤f≤ 55 Input sampling = 10 KHz
Ripple = δ1=0.1, δ2=0.001.
5. Design a linear pahse FIR filter that satisfies the following specifications based on a single-
stage and two-stage multirate structure.
Input sampling rate: 10K Hz
Passband : 0 ≤ F ≤ 60
Transition band : 60 ≤ F ≤ 65
Ripple : δ1 =10-1,δ2 = 10-3

16.1: Unit-Wise Short and Long Answer Question Bank

DEPARTMENT OF Electronics and Communication Engineering

 Unit - 1

1. Discuss the condition for causality and stability?

2. State the Sampling theorem?
3. Discuss the advantages of DSP?
4. What is the need for a sample and hold circuit?
5. Define Parseval’s theorem?
6. State commutative and associative properties of linear convolution?
7. Calculate the energy and power of the signal Aejwnu(n)?
8. What is aliasing effect? How can it be avoided?
9. Define Z-transform and region of convergence (ROC)?
10. What are the properties of Z-transform?
11. What are the properties of ROC?
12. With reference to Z-transform, what is initial value theorem and final value theorem?
13.Distinguish between DFS and DTFT?
UNIT- 2

1. Distinguish between DFT and DTFT?

2. State and prove time shifting property of DFT?
3. Explain the relation between DFT and Z-transform?
4. Distinguish between linear and circular convolution of two sequences?
5. Differences between overlap-add and overlap-save method?
6. Why FFT is needed?
7. What are the differences between DIT-FFT and DIF-FFT algorithms?
8. How can we calculate IDFT using the FFT algorithm?
9. How is linear convolution calculated using the DFT?
10. Discuss about power density spectrum of periodic signals?
UNIT-3

1. Give any two properties of Butterworth filter.

2. What are the properties of Chebyshev filter.
3. What is meant by frequency warping? What is the cause of this effect?
4. Distinguish between Butterworth and Chebyshev filter.
5. Distinguish between IIR and FIR filter?
6. What is bilinear transformation?
7. What is impulse invariant method of designing IIR filter?
8. What are the advantages of Chebyshev filter over Butterworth filter?
9. Write a note on pre-warping?
10. Distinguish between analog and digital filters?
UNIT -4

1. What is meant by FIR filter? What are the advantages of FIR filter?
2. List the design techniques for FIR filter design?
3. What is Gibbs phenomenon?

DEPARTMENT OF Electronics and Communication Engineering

 4. Under what condition, the finite duration sequence h(n) will yield constant group delay in
the frequency response characteristics and not the phase delay?
5. What are the desirable characteristics of windows?
6. Compare Hamming window with Kaiser window.
7. Draw impulse response of ideal low pass filter?
8. What is the principle of designing FIR filter using frequency sampling method?
9. What is the necessary and sufficient condition for linear phase characteristics in FIR
filter?
10. Explain the procedure for designing FIR filters using windows?

UNIT-5

1. Explain the term up-sampling and down-sampling ?

2. What are the advantages of multirate DSP?
3. What is decimation by a factor D?
4. What is interpolation by a factor I?
5. What happens if I violate Nyquist criterion in decimation or down sampling?
6. What are the finite word length effects in digital filters?
7. Find the spectrum of exponential signal decimated by a factor of 2?
8. Find the spectrum of exponential signal interpolated by a factor of 2?
9. List the errors which arise due to quantization process.
10. What are the advantages of floating point arithmetic.

16.2: Unit wise Quiz Questions

UNIT1
1.Discrete time signals are _________ in time and _______________ in amplitude.
2.The representation of a signal by mathe matical expression is known as____________
3.An LTI system is one which satisfies te properties of ____________ and ___________.
4.The response of the system due to input alone ,when the initial conditions are neglected is called
the ____________ of the systems.forced response
5.The Z-Transform converts _____________ equaton into ____________ equation.
(Difference,Algebraic)
Unit 2
1.The DTFT is a periodic ____________ function of ωwith a period of(Continuous ,2П)
2.The DFI X(K) of a discrete time sequence x(n) is defined as____________
3.____________ is known as tfe twiddle factor.

DEPARTMENT OF Electronics and Communication Engineering

4.The basic FFT algorithms are _________ and ____________.
5.The computation of 32-point DFT by radix -2 DIT FFT involves ___________ stages of
computation.(five)
Unit3
1.Filters designed by considering _______________ samples of the impulse reponce are called IIR
filters.
2.In______________transformation the impulse response of digital filter is the sampled version
of the impulse response of analog version.(impulse invariant)
3.The phenomenon of high frequency components acquiring the identity of low frequency
components is called_______________.(aliasing)
4.The type-2 chebyshev response is also called ____________ response.(inverse chebyshev)
5.The two popular techniques used to approximate the ideal frequency response are __________ and
_____________approximatoiopns.(Butterworth,chebyshev).
Multiple choice Questions:
1.IIR filters are designed by considering all the
a).Infinite samples of frequency response.
b).finite samples of frequency response
c).Infinite samples of impulse respoinse
d).none of the above.
2.For the analog and Digital IIR filters o be causal, the number of zeros should be
a).>=number of poles b).<=numbers of poles c).=number of poles d).Zero.
3.The zeros of the butterworth filters exists at
a)left half of s-plane b).origin c).infinity d).Right half of s-plane
4.The poles of Butter worth transfer unction lie,
a).Symmetrically on a circle in S-plane
b). Symmetrically on an ellipse in S-plane
c).Antisymmetrically on a circle in S-plane
d).Antisymmetrically on a ellipse in S-plane
5.The relation between analog and digital frequency is nonlinear in case o
a).Impulse invariant transformation b).Bilinear transformation c).Frequency sampling d).All the
above.

DEPARTMENT OF Electronics and Communication Engineering

Unit 4
1.The ___________is due to nonlinear phase characteristics of the filter.
2.In FIR filters ______________function is a linear function of ω
3.In rectangular window the width of main-lobe is equal to___________
4.The _________ window spectrum has the highest attenuation for side lobes.
5.The width of the main-lobe in window spectrum can be reduced by increasing the length
of___________.
Multiple choice Questions:
1.1.The frequency response of a digital filter is periodic in the range
a).0<ω<2π b).-π<ω<π c).0<ω<π d). 0<ω<2π or -π<ω<π
2.If ωc is the cutoff frequency of highpass filter,then the response lies only in the range of,
a)-ωc<ω<ωc b). a)-ωc<ω<π a)-π<ω<ωc a)-ωc<ω<π
3.Raised cosine windows also called generalized
a).Hamming window b).Hanning Window c).Rectangular window
d).Blackman window
4.Symmetric impulse response having odd numbe of samples,N=7 with centre of symmetry α is
equal to
a).2 b).5 c).3.5 d).3
5.The Symmetric impulse response having even number of samples cannotbe used to design,
a).Lowpass filter b).Band pass filter c).High pass filter d).Band stop filter.

Unit 5
1.The processing of a signal at different sampling rates is called_________________.
2.The ________ is the process of increasing the sampling rate.
3. The ________ is the process of decreasing the sampling rate.
4.The process of dividing a filter into a number of sub-filters is called_____________
5.The ___________banks are filter banks with complementary frequency resonse.

Multiple choice Questions:

1.If x(n) and y(n ) are input and output of a decimator with sampling rate conversion factor A,then,

DEPARTMENT OF Electronics and Communication Engineering

a).y(n)=x(n-A) b).y(n)=x(n/A) c).y(n)=x(n+A) d).y(n)=x(An)
2.To avoid aliasing at output during decimation by D,the input signal of a decimator should be
bandlimited to,
a).∏/2D a).2∏/D a).∏/D a).∏/D 2

3.If x(n) and y(n ) are input and output of a interpolator with sampling rate conversion factor B,then,
a).y(n)=x(Bn) b).y(n)=x(n/B) c).y(n)=x(n)/B d).y(n)=Bx(n)
4.To eliminate multiple images at the output, during interpolation by I,the output is filtered to have
a bandwidth of,
a).∏I b). ∏I c).I/∏ d).∏/I2
5.If A and B are integer sampling rate conversion factor fr decimation and interpolation
respectively,then sampling rate convertion factor for conversion by rational factor is,
a)A/B b)B/A c)A2/B d).B/A2

17. Detailed notes (some material is presented here and remaining material is
provided in pdf format)

UNIT-I

Introduction to Digital Signal Processing

What is DSP?

DSP, or Digital Signal Processing, as the term suggests, is the processing of signals by digital
means. A signal in this context can mean a number of different things. Historically the origins of
signal processing are in electrical engineering, and a signal here means an electrical signal
carried by a wire or telephone line, or perhaps by a radio wave. More generally, however, a
signal is a stream of information representing anything from stock prices to data from a remote-
sensing satellite.

DEPARTMENT OF Electronics and Communication Engineering

Analog and digital signals

In many cases, the signal is initially in the form of an analog electrical voltage or current,
produced for example by a microphone or some other type of transducer. In some situations the
data is already in digital form - such as the output from the readout system of a CD (compact
disc) player. An analog signal must be converted into digital (i.e. numerical) form before DSP
techniques can be applied. An analog electrical voltage signal, for example, can be digitized
using an integrated electronic circuit (IC) device called an analog-to-digital converter or ADC.
This generates a digital output in the form of a binary number whose value represents the
electrical voltage input to the device.

Signal processing

Signals commonly need to be processed in a variety of ways. For example, the output signal
from a transducer may well be contaminated with unwanted electrical "noise". The electrodes
attached to a patient's chest when an ECG is taken measure tiny electrical voltage changes due to
the activity of the heart and other muscles. The signal is often strongly affected by "mains
pickup" due to electrical interference from the mains supply. Processing the signal using a filter
circuit can remove or at least reduce the unwanted part of the signal. Increasingly nowadays the
filtering of signals to improve signal quality or to extract important information is done by DSP
techniques rather than by analog electronics.

Development of DSP

The development of digital signal processing dates from the 1960's with the use of mainframe
digital computers for number-crunching applications such as the Fast Fourier Transform (FFT),
which allows the frequency spectrum of a signal to be computed rapidly. These techniques were
not widely used at that time, because suitable computing equipment was available only in
universities and other scientific research institutions.

Digital Signal Processors (DSPs)

The introduction of the microprocessor in the late 1970's and early 1980's made it possible for
DSP techniques to be used in a much wider range of applications. However, general-purpose
microprocessors such as the Intel x86 family are not ideally suited to the numerically-intensive
requirements of DSP, and during the 1980's the increasing importance of DSP led several major
electronics manufacturers (such as Texas Instruments, Analog Devices and Motorola) to develop
Digital Signal Processor chips - specialized microprocessors with architectures designed
specifically for the types of operations required in digital signal processing. (Note that the
acronym DSP can variously mean Digital Signal Processing, the term used for a wide range of
techniques for processing signals digitally, or Digital Signal Processor, a specialized type of
microprocessor chip). Like a general-purpose microprocessor, a DSP is a programmable device,
with its own native instruction code. DSP chips are capable of carrying out millions of floating
point operations per second, and like their better-known general-purpose cousins, faster and
more powerful versions are continually being introduced.

DEPARTMENT OF Electronics and Communication Engineering

Applications of DSP

DSP technology is nowadays commonplace in such devices as mobile phones, multimedia

computers, video recorders, CD players, hard disc drive controllers and modems, and will soon
replace analog circuitry in TV sets and telephones. An important application of DSP is in signal
compression and decompression. In CD systems, for example, the music recorded on the CD is
in a compressed form (to increase storage capacity) and must be decompressed for the recorded
signal to be reproduced. Signal compression is used in digital cellular phones to allow a greater
number of calls to be handled simultaneously within each local "cell". DSP signal compression
technology allows people not only to talk to one another by telephone but also to see one another
on the screens of their PCs, using small video cameras mounted on the computer monitors, with
only a conventional telephone line linking them together.

Although the mathematical theory underlying DSP techniques such as Fast Fourier and Hilbert
Transforms, digital filter design and signal compression can be fairly complex, the numerical
operations required to implement these techniques are in fact very simple, consisting mainly of
operations that could be done on a cheap four-function calculator. The architecture of a DSP chip
is designed to carry out such operations incredibly fast, processing up to tens of millions of
samples per second, to provide real-time performance: that is, the ability to process a signal
"live" as it is sampled and then output the processed signal, for example to a loudspeaker or
video display. All of the practical examples of DSP applications mentioned earlier, such as hard
disc drives and mobile phones, demand real-time operation.

The major electronics manufacturers have invested heavily in DSP technology. Because they
now find application in mass-market products, DSP chips account for a substantial proportion of
the world market for electronic devices. Sales amount to billions of dollars annually, and seem
likely to continue to increase rapidly.

DEPARTMENT OF Electronics and Communication Engineering

Introduction to DSP

Most DSP applications deal with analog signals.

 the analog signal has to be converted to digital form

The analog signal - a continuous variable defined with infinite precision - is converted to a
discrete sequence of measured values which are represented digitally.

Information is lost in converting from analogue to digital, due to:

 inaccuracies in the measurement

 uncertainty in timing
 limits on the duration of the measurement

These effects are called quantization errors.

The continuous analog signal has to be held before it can be sampled.

Otherwise, the signal would be changing during the measurement.

Only after it has been held can the signal be measured, and the measurement converted to a
digital value.

DEPARTMENT OF Electronics and Communication Engineering

The sampling results in a discrete set of digital numbers that represent measurements of the
signal - usually taken at equal intervals of time.

Note that the sampling takes place after the hold. This means that we can sometimes use a slower
Analogue to Digital Converter (ADC) than might seem required at first sight. The hold circuit
must act fast - fast enough that the signal is not changing during the time the circuit is acquiring
the signal value - but the ADC has all the time that the signal is held to make its conversion.

We don't know what we don't measure.

In the process of measuring the signal, some information is lost.

Sometimes we may have some a priori knowledge of the signal, or be able to make some
assumptions that will let us reconstruct the lost information.

DEPARTMENT OF Electronics and Communication Engineering

UNIT-II
Discrete Fourier Transform (DFT):

1. Introduction

1. x(t) --- Continuous-time signal

X(f) --- Fourier Transform, frequency characteristics

Can we find
Z (r)
if we don’t have a mathematical equation for x(t) ? No!

2. What can we do?

(1) Sample x(t) =>

x0, x1, … , xN-1 over T (for example 1000 seconds)

Sampling period (interval) Z(r)

N (samples) over T => Z ( r )

Can we have infinite T and N? Impossible!

(2) Discrete Fourier Transform (DFT):

=>
Z (r)
for the line spectrum at frequency
Z (r)
DEPARTMENT OF Electronics and Communication Engineering
3. Limited N and T =>

limited frequency resolution

Z(r)
limited frequency band (from Z ( r ) to Z(r) in Fourier transform to):

Z (r)

4.
Z (r) ---- periodic function (period N)

x(t) --- general function

 sampling and inverse transform

xn --- periodic function

5.
Z ( r ) line spectrum)
Z (r) period function (period N)

2 Error Sources in the DFT

1. Preparations

(1) Ideal sampling waveform

Z (r) :

DEPARTMENT OF Electronics and Communication Engineering

 Z (r) Center of the
window

in DFT Z (r)

(2) Rectangular Pulse (window)

Z ( r ) Z (r )
when t0 = 0
Z (r)
(3) Z ( r )

(4) Z ( r )

2. Illustration of Error Sources

Example 10-1: Continuous – time signal: two-sided exponential signal

Z (r)
Its Fourier transform

DEPARTMENT OF Electronics and Communication Engineering

 Z (r)
(1) If we sample x(t) in
Z (r) with sampling frequency fs :

sampled signal
Z (r)
its Fourier transform:

Z (r)
sampling => (1) possible overlapping if Z ( r ) is not held.
(2) periodic function, introduce frequencies beyond fs .

(2) Limited T (over which x(t) is sampled to collect data for DFT)

window
Z(r)

Z (r)
Fourier transform given by sampled data in limited window (T)

DEPARTMENT OF Electronics and Communication Engineering

  Z ( r ) is a worse estimate of X(f) than Xs(f) due to the introduction of

( Tsinc( Tf ) ) for convolution!

Effect of limited T

(3) Dose DFT give Z ( r ) for every f ?

No! only discrete frequencies.

DFT as an estimate for X(f): even worse than Z ( r ) due to the limited frequency
resolution.

DEPARTMENT OF Electronics and Communication Engineering

3. Effect of sampling frequency (or number of points) on accuracy when T is given: Example

use
Z(r) for
Z(r)
4. Effect of T (window size)

Compare
Z(r) and
Z(r) for
Z(r)

DEPARTMENT OF Electronics and Communication Engineering

5. DFT Errors superposition of the analog signal spectrum X(f) a
(1) Aliasing

Z (r)
Caused by sampling

Overlapping of X(f) and its translates: aliasing (sampling effect)

(2) Leakage Effect

limited window size T (

Z(r) )
Z ( r )

DEPARTMENT OF Electronics and Communication Engineering

 worse than Xs(f) as approximation of X(f).

Z ( r ) : contribution of Z ( r ) to Z ( r ) : determined by

weight Z ( r )

frequency energy “leaks” from one frequency to another!

(3) Picket – Fence Effect:

As an estimation of X(f), does Z ( r ) have picket fence effect? No!
DFT: discrete frequencies (not blocked by the fence).

6. Minimization of DFT Error Effects.

 Major ways: increase T and fs

Problem: DFT for large N.

3 Examples Illustrating the computation of the DFT

(Preparation for Mathematical Derivation of FFT)

1. DFT Algorithm

DEPARTMENT OF Electronics and Communication Engineering

 Z (r )
Denote Z (r) , then

Z (r)
Properties of Z (r) :

(1) Z ( r )
(2) Z ( r )

Z (r)
(3) Z ( r )
Z ( r )
Z ( r )
2. Examples
Example 10-3: Two-Point DFT

x(0), x(1):
Z (r)
Z (r)

DEPARTMENT OF Electronics and Communication Engineering

 Z (r)
Example -4: Generalization of derivation in example 10-3 to a four-point DFT

x(0), x(1), x(2), x(3)

Z (r)
Z ( r )
Z (r)
Z (r)
DEPARTMENT OF Electronics and Communication Engineering
 Z (r)

Z (r)

Two – point DFT

If we denote z(0) = x(0), z(1) = x(2) => Z(0) = z(0) + z(1) = x(0) + x(2)

Z(1) = z(0) - z(1) = x(0) - x(2)

v(0) = x(1), v(1) = x(3) => V(0) = v(0) + v(1) = x(1) + x(3)

V(1) = v(0) - v(1) = x(1) - x(3)

Four – point DFT Two-point DFT

DEPARTMENT OF Electronics and Communication Engineering

  X(0) = Z(0) + V(0)

X(1) = Z(1) + (-j)V(1)

X(2) = Z(0) - V(0)

X(3) = Z(1) + jV(1)

 One Four – point DFT Two Two – point DFT

Fast Fourier Transform (FFT):

4 Mathematical Derivation of the FFT

4A Decimation-in-Time FFT Algorithm

x(0), x(1), … , x(N-1) Z (r)

DEPARTMENT OF Electronics and Communication Engineering

 =>
Z (r)
Z (r)
Z (r)
( G(k): N/2 point DFT output (even indexed), H(k) : N/2 point DFT output (odd
indexed))

Z (r)
DEPARTMENT OF Electronics and Communication Engineering
Question: X(k) needs G(k), H(k), k=… N-1

How do we obtain G(k), H(k), for k > N/2-1 ?

G(k) = G(N/2+k) k <= N/2-1

H(k) = H(N/2+k) k <= N/2-1

Future Decimation
g(0), g(1), …, g(N/2-1) G(k)

h(0), h(1), …, h(N/2-1) H(k)

DEPARTMENT OF Electronics and Communication Engineering

 Z(r) Z(r)
even indexed g odd indexed g

(N/4 point) (N/4 point)

Z (r)

Z (r)
=> Z (r )

Similarly,

Z ( r )
even indexed odd indexed

h (N/4 point) h (N/4 point)

For 8 – point

DEPARTMENT OF Electronics and Communication Engineering

Z (r) Z(r)

DEPARTMENT OF Electronics and Communication Engineering

4B Decimation-in-Frequency FFT Algorithm

x(0), x(1), … , x(N-1) Z (r)

DEPARTMENT OF Electronics and Communication Engineering

 Z (r)
let m = n-N/2 (n = N/2+m) n = N/2 => m = N/2-N/2 = 0

n = N-1 => m = N-1-N/2 = N/2-1

Z (r)

Z (r)
Z (r)
N/2 point DFT

(
H ( z )=
0.98−3.43z−1+4.96z−2−3.43z−3+0.98z−4
1−3.47z−1+4.96z−2−3.39z−3+0.96z−4 ) Z (r )

DEPARTMENT OF Electronics and Communication Engineering

 ( )
−1 −2 −3 −4
0.98−3.43 z +4.96 z −3.43 z +0.98 z
H (z)= −1 −2 −3 −4
1−3.47 z +4.96 z −3.39 z +0.96 z
H ( z )=
( 0 . 98− 3. 43 z−1 + 4 . 96 z−2− 3. 43 z−3 + 0 .98 z− 4
1−3. 47 z −1 + 4 . 96 z−2 −3 .39 z−3 + 0. 96 z−4 )
X(k) : N-point DFT of x(0), …, x(N)  two N/2 point DFT

DEPARTMENT OF Electronics and Communication Engineering

One N/2 point DFT => two N/4 point DFT

… two point DFTs

Consider N/2 point DFT

y(0), y(1), …, y(N/2-1)

(
H ( z )=
0.98−3.43z−1+4.96z−2−3.43z−3+0.98z−4
−1 −2 −3
1−3.47z +4.96z −3.39z +0.96z −4 )
DEPARTMENT OF Electronics and Communication Engineering
 ( )
−1 −2 −3 −4
0.98−3.43z +4.96z −3.43z +0.98z
H ( z )= −1 −2 −3 −4
1−3.47z +4.96z −3.39z +0.96z

( )
−1 −2 −3 −4
0.98−3.43z +4.96z −3.43z +0.98z
H ( z )= −1 −2 −3 −4
1−3.47z +4.96z −3.39z +0.96z

DEPARTMENT OF Electronics and Communication Engineering

4 C Computation

N – point DFT : 4N(N-1) real multiplications

4N(N-1) real additions

N – point FFT : 2Nlog2N real multiplications

(N = 2m) 3Nlog2N real additions

Computation ration

DEPARTMENT OF Electronics and Communication Engineering

 ( )
−1 −2 −3 −4
0.98−3.43z +4.96z −3.43z +0.98z
H ( z )=
1−3.47z−1+4.96z−2−3.39z−3+0.96z−4

5 Properties of the DFT

Assumptions

(1)

( 0.98−3.43z−1+4.96z−2−3.43z−3+0.98z−4
H ( z)= −1 −2 −3 −4
1−3.47z +4.96z −3.39z +0.96z ) (0.98−3.43 z−1+4.96 z−2−3.43 z−3+0.98 z−4
H (z)= −1 −2 −3 −4
1−3.47 z +4.96 z −3.39 z +0.96 z )
(2) A, B: arbitrary constants

(3) Subscript e:

Subscript o:

DEPARTMENT OF Electronics and Communication Engineering

 x (n) : even about
e
(
0. 98−3. 43z−1+4. 96z−2−3. 43z−3+0.98z−4
H ( z)= −1 −2 −3 −4
1−3.47z +4. 96z −3.39z +0.96z )
x (n) : odd about
o
(
0. 98−3. 43z−1+4. 96z−2−3. 43z−3+0.98z−4
H ( z)= −1 −2 −3 −4
1−3.47z +4. 96z −3.39z +0.96z )
N = 10, xe(n)

( )
−1 −2 −3 −4
0.98−3.43z +4.96z −3.43z +0.98z
H ( z)= −1 −2 −3 −4
1−3.47z +4.96z −3.39z +0.96z
N = 9, xe(n)

( )
−1 −2 −3 −4
0.98−3.43 z +4.96 z −3.43 z +0.98 z
H (z)= −1 −2 −3 −4
1−3.47 z +4.96 z −3.39 z +0.96 z
(4) Any real sequence can be expressed in terms of its even and odd parts according to

DEPARTMENT OF Electronics and Communication Engineering

 H ( z )=
(
0.98−3.43 z−1+4.96 z−2−3.43 z−3+0.98 z−4
1−3.47 z−1+4.96 z−2−3.39 z−3+0.96 z−4 )
even odd

Question 1: x(n) = 1/2[ ] +1/2 [ ] ? Yes!

Question 2: x(n) + x(N-n) even ?

x(n) - x(N-n) odd ?

Example: N = 9 => N/2 = 4.5

Consider n = 2

x(2) + x(9-2) = x(2) + x(7)

is x(2) + x(7) = x(4.5+(4.5-2)) +x(4.5-(7-4.5))?

4.5 + (4.5-2) = 9-2 = 7

4.5 - (7-4.5) = 9-7 = 2

x(2) + x(7) = x(7) + x(2) ?

Yes! => x(n) + x(N-n) even

Is x(2) - x(7) = - [x(7) + x(2)] ?

Yes! => x(n) - x(N-n) odd

(5) subscript r : xr(n) a real sequence

subscript i : xi(n) Imaginary part of a complex sequence

(6) (
H ( z )=
0.98−3.43 z−1+4.96 z−2−3.43 z−3+0.98 z−4
1−3.47 z−1+4.96 z−2−3.39 z−3+0.96 z−4 )
left right side:

side: DFT

DEPARTMENT OF Electronics and Communication Engineering

 sequence

(7) sequences are assumed periodically repeated if necessary

Properties

1. Linearity :
H ( z )= ( 0 . 98− 3. 43 z−1 + 4 . 96 z−2− 3. 43 z−3 + 0 .98 z− 4
1−3. 47 z −1 + 4 . 96 z−2 −3 .39 z−3 + 0. 96 z−4 )

2. Time Shift:
H ( z )= ( 0 . 98− 3. 43 z−1 + 4 . 96 z−2− 3. 43 z−3 + 0 .98 z− 4
1−3. 47 z −1 + 4 . 96 z−2 −3 .39 z−3 + 0. 96 z−4 )

3. Frequency Shift:

(
H ( z )=
0 . 98−3. 43 z−1 +4 . 96 z−2−3. 43 z−3 +0 .98 z− 4
1−3. 47 z −1 +4 . 96 z−2 −3 .39 z−3 +0. 96 z−4 )
4. Duality : (
H ( z )=
0.98−3.43 z−1+4.96 z−2−3.43 z−3+0.98 z−4
1−3.47 z −1+4.96 z−2−3.39 z−3+0.96 z−4 )

why?
H ( z )=
(
0.98−3.43z−1+4.96z−2−3.43z−3+0.98z−4
1−3.47z−1+4.96z−2−3.39z−3+0.96z−4 )
(
H ( z )=
0.98−3.43 z−1+4.96 z−2−3.43 z−3+0.98 z−4
1−3.47 z−1+4.96 z−2−3.39 z−3+0.96 z−4 )
DFT of x(m)

DEPARTMENT OF Electronics and Communication Engineering

 ( )
−1 −2 −3 −4
0.98−3.43 z +4.96 z −3.43 z +0.98 z
H (z)= −1 −2 −3 −4
1−3.47 z +4.96 z −3.39 z +0.96 z
H ( z )=
(
0 . 98−3. 43 z−1 +4 . 96 z−2−3. 43 z−3 +0 .98 z− 4
1−3. 47 z −1 +4 . 96 z−2 −3 .39 z−3 +0. 96 z−4 )
5. Circular convolution

(
H ( z )=
0. 98−3. 43 z−1 +4 . 96 z−2−3. 43 z−3 +0 .98 z− 4
1−3.47 z −1 +4. 96 z−2 −3 .39 z−3 +0.96 z−4 ) circular convolution

6. Multiplication

H ( z )=
(
0.98−3.43 z−1 +4.96 z−2−3.43 z−3 +0.98 z− 4
1−3.47 z −1 +4.96 z−2−3.39 z−3 +0.96 z−4 )
7. Parseval’s Theorem

H ( z )=
( 0.98−3.43 z−1+4.96 z−2−3.43 z−3+0.98 z−4
1−3.47 z−1+4.96 z−2−3.39 z−3+0.96 z−4 )
8. Transforms of even real functions:
(
H ( z )=
0.98−3.43 z−1+4.96 z−2−3.43 z−3+0.98 z−4
1−3.47 z−1+4.96 z−2−3.39 z−3+0.96 z−4 )
(the DFT of an even real sequence is even and real )

DEPARTMENT OF Electronics and Communication Engineering

9. Transform of odd real functions:

H ( z )= ( 0.98−3.43 z−1+4.96 z−2−3.43 z−3+0.98 z−4

1−3.47 z −1+4.96 z−2−3.39 z−3+0.96 z−4 )
(the DFT of an odd real sequence is odd and imaginary )

10. z(n) = x(n) + jy(n)

z(n)  Z(k) = X(k) + jY(k)

Example -7

H ( z )= ( 0 . 98−3. 43 z−1 +4 . 96 z−2−3. 43 z−3 +0 .98 z− 4

1−3. 47 z −1 +4 . 96 z−2 −3 .39 z−3 +0. 96 z−4 )
Four – point DFT for x(0), x(1), x(2), x(3):

X(0) = [x(0) + x(2)] + [x(1) + x(3)]

X(1) = [x(0) - x(2)] + (-j)[x(1) - x(3)]

X(2) = [x(0) + x(2)] - [x(1) + x(3)]

X(3) = [x(0) - x(2)] + j[x(1) - x(3)]

For (
H ( z )=
0.98−3.43 z−1 +4.96 z−2−3.43 z−3 +0.98 z− 4
1−3.47 z −1 +4.96 z−2−3.39 z−3 +0.96 z−4 ) =>

H ( z )=
(
0.98−3.43 z−1+4.96 z−2−3.43 z−3+0.98 z−4
1−3.47 z−1+4.96 z−2−3.39 z−3+0.96 z−4 )
DEPARTMENT OF Electronics and Communication Engineering
 For (
H ( z )=
0.98−3.43 z−1 +4.96 z−2−3.43 z−3 +0.98 z− 4
1−3.47 z −1 +4.96 z−2−3.39 z−3 +0.96 z−4 ) =>

H ( z )=
(
0.98−3.43 z−1+4.96 z−2−3.43 z−3+0.98 z−4
1−3.47 z−1+4.96 z−2−3.39 z−3+0.96 z−4 )


( )
0.98−3.43z−1+4.96z−2−3.43z−3+0.98z−4
H(z)= −1 −2 −3 −4
1−3.47z +4.96z −3.39z +0.96z

Example -8

DFT of (
H ( z )=
0.98−3.43z−1+4.96z−2−3.43z−3+0.98z−4
1−3.47z−1+4.96z−2−3.39z−3+0.96z−4 ):

H ( z )=
( 0.98−3.43 z−1 +4.96 z−2−3.43 z−3 +0.98 z− 4
1−3.47 z −1 +4.96 z−2 −3.39 z−3 +0.96 z−4 )
Time-shift property

H ( z )=
(
0.98−3.43z−1+4.96z−2−3.43z−3+0.98z−4
1−3.47z−1+4.96z−2−3.39z−3+0.96z−4 )
Example -9: Circular Convolution

H ( z )=
( 0 . 98− 3. 43 z−1 + 4 . 96 z−2− 3. 43 z−3 + 0 .98 z− 4
1−3. 47 z −1 + 4 . 96 z−2 −3 .39 z−3 + 0. 96 z−4 )

Define
H ( z )=
(
0.98−3.43 z−1 +4.96 z−2−3.43 z−3 +0.98 z− 4
1−3.47 z −1 +4.96 z−2 −3.39 z−3 +0.96 z−4 )

DEPARTMENT OF Electronics and Communication Engineering

 ( )
−1 −2 −3 −4
0.98−3.43z +4.96z −3.43z +0.98z
H ( z)= −1 −2 −3 −4
1−3.47z +4.96z −3.39z +0.96z

( )
−1 −2 −3 −4
0.98−3.43 z +4.96 z −3.43 z +0.98 z
H (z)= −1 −2 −3 −4
1−3.47 z +4.96 z −3.39 z +0.96 z
6 Applications of FFT

1. Filtering
x(0), …, x(N-1) FFT (DFT) =>

X(0), … , X(1), … , X(N-1)

X(k): Line spectrum at

H ( z )= ( 0 . 98−3. 43 z−1 +4 . 96 z−2−3. 43 z−3 +0 .98 z− 4
1−3. 47 z −1 +4 . 96 z−2 −3 .39 z−3 +0. 96 z−4 )
(Over T: x(0), …, x(N-1) are sampled.)

Inverse DFT:

( 0.98−3.43z−1+4.96z−2−3.43z−3+0.98z−4
H ( z)= −1 −2 −3 −4
1−3.47z +4.96z −3.39z +0.96z )
DEPARTMENT OF Electronics and Communication Engineering
 0.98−3.43z−1+4.96z−2−3.43z−3+0.98z−4
−1 −2 −3 −4 (
H ( z )=
)
Frequencies with 1−3.47z +4.96z −3.39z +0.96z have been filtered!

Example -10

H ( z )=
(
0.98−3.43 z−1+4.96 z−2−3.43 z−3+0.98 z−4
1−3.47 z −1+4.96 z−2−3.39 z−3+0.96 z−4 )
x(0), x(1), …, x(7)

H ( z )=
(
0.98−3.43 z−1+4.96 z−2−3.43 z−3+0.98 z−4
1−3.47 z−1+4.96 z−2−3.39 z−3+0.96 z−4 )
H(z)= −1 −2 −3 −4
How to filter frequency higher than 1−3.47z +4.96z −3. 9z +0.96z ?
()
0.98−3.4 z−1+4.96z−2−3.4 z−3+0.98z−4

2. Spectrum Analyzers
Analog oscilloscopes => time-domain display

Spectrum Analyzers: Data Storage, FFT

3. Energy Spectral Density

x(0), …, x(N-1): its energy definition

( 0.98−3.43z−1+4.96z−2−3.43z−3+0.98z−4
H ( z )= −1 −2 −3 −4
1−3.47z +4.96z −3.39z +0.96z )
Parseval’s Theorem

( 0.98−3.43z−1+4.96z−2−3.43z−3+0.98z−4
H ( z )= −1 −2 −3 −4
1−3.47z +4.96z −3.39z +0.96z )
DEPARTMENT OF Electronics and Communication Engineering
 UNIT-III
Design of DIGITAL FILTERS (IIR)- Structures of IIR systems:

3.1. Introduction:

A general causal digital filter has the difference equation:

N M

y[n] = åå a i x[n-i] - åå b k y[n-k]

i=0 k=1

which is of order max{ N,M }, and is recursive if any of the b j coefficients are non-zero. A
second order recursive digital filter therefore has the difference equation:

y[n] = a 0 x[n] + a 1 x[n-1] + a 2 x[n-2] - b 1 y[n-1] - b 2 y[n-2]

A digital filter with a recursive linear difference equation can have an infinite impulse-response.
Remember that the frequency-response of a digital filter with impulse-response {h[n]} is:
¥¥

H(e j  ) = åå h[n]e - j  n

n=-¥¥

DEPARTMENT OF Electronics and Communication Engineering

3.2. The z-transform:

Consider the response of a causal stable LTI digital filter to the special sequence {z n } where z is
a complex. If {h[n]} is the impulse-response, by discrete time convolution, the output is a
sequence {y[n]} where

¥¥ ¥¥

y[n] = åå h[k] z n - k = z n åå h[k] z - k

k=-¥¥ k=-¥¥

¥¥

= zn H(z) with H(z) = åå h[k] z - k

k=-¥¥

The expression obtained for H(z) is the ‘z-transform’ of the impulse-response. H(z) is a complex
number when evaluated for a given complex value of z.

It may be shown that for a causal stable system, H(z) must be finite when evaluated for a
complex number z with modulus greater than or equal to one.
H ( z )=
( 0 . 98− 3. 43 z−1 + 4 . 96 z−2− 3. 43 z−3 + 0 .98 z− 4
1−3. 47 z −1 + 4 . 96 z−2 −3 .39 z−3 + 0. 96 z−4 )
it is clear that replacing z by e j  in H(z) gives H(ejWW) .

DEPARTMENT OF Electronics and Communication Engineering

DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
3.3. The ‘z-plane’:

It is useful to represent complex numbers on an ‘Argand diagram’ as illustrated below. The

main reason for doing this is that the modulus of the difference between two complex numbers
a+jb and c+jd say i.e. | (a+jb) - (c+jd) |is represented graphically as the length of the line
between the two complex numbers as plotted on the Argand diagram.

Imaginary part

Rej

-3+3j
R


Real part

1-2j

If one of these complex numbers, c +jd say is zero i.e. 0+j0, then the modulus of the other
number |a+jb| is the distance of a+jb from the origin 0+j0 on the Argand diagram.

Of course, any complex number, a+jb say, can be converted to polar form Re j where R= |a+jb|
and = tan-1(b/a). Plotting a complex number expressed as Re j on an Argand diagram is also
illustrated above. We draw an arrow of length R starting from the origin and set at an angle 
from the ‘real part’ axis (measured anti-clockwise). Re j is then at the tip of the arrow. In the
illustration above,  is about /4 or 45 degrees. If R=1, Re j = ej and on the Argand diagram
would be a point at a distance 1 from the origin. Plotting e j for values  in the range 0 to 2
(360O) produces points all of which lie on a ‘unit circle’ , i.e. a circle of radius 1, with centre at
the origin.

Where the complex numbers plotted on an Argand diagram are values of z for which we are
interested in H(z), the diagram is referred to as ‘the z-plane’. Points with z = e jWW lie on a unit
circle, as shown in Fig 5.1. Remember that |e j  | = |cos() +jsin()| = [cos2() + sin2()] =
1. Therefore evaluating the frequency-response H(e i) for  in the range 0 to  is equivalent to

DEPARTMENT OF Electronics and Communication Engineering

evaluating H(z) for z =e j  which goes round the upper part of the unit circle as  goes from 0 to
.

3.4. Relating H(z), signal-flow graphs and difference equations

This is surprisingly straightforward. Consider non-recursive and recursive difference equations

separately.

Example 3.1: Find H(z) for the difference equation: y[n] = x[n] + x[n-1]

Solution: The impulse response is: {h[n]} = { ... , 0, 1, 1, 0, ... }

(
H ( z )=
0.98−3.43 z−1 +4.96 z−2−3.43 z−3 +0.98 z− 4
1−3.47 z −1 +4.96 z−2 −3.39 z−3 +0.96 z−4 )
Example 3.2:

Find H(z) for the recursive difference equation: y[n] = a 0 x[n] + a 1 x[n-1] - b 1 y[n-1]

Solution:

The method used in Example 5.1 is not so easy because the impulse-response can now be
infinite. Fortunately there is another way. Remember that if x[n] = z n then y[n] = H(z) z n ,
y[n-1] = H(z) z n - 1 etc. Substitute into the difference equation to obtain:

H(z) z n = a 0 z n + a 1 z n - 1 - b 1 H(z) z n - 1

DEPARTMENT OF Electronics and Communication Engineering

 -1
a0 + a1 z

Therefore, H(z) = ¾¾¾¾¾¾¾¾¾¾

1 + b1 z-1

except when z = - b 1 . When z = -b 1 , H(z) = ¥¥.

By the same method, H(z) for a general digital filter whose difference-equation was given
earlier is:

a 0 + a 1 z - 1 + a 2 z - 2 + ... + a N z - N

H(z) = ¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾
(with b 0 = 1)

b 0 + b 1 z - 1 + b 2 z - 2 + ... + b M z - M

Given H(z) in this form, we can easily go back to its difference-equation and hence its signal-
flow graph, as illustrated by the following example.

Example 3.3: Give a signal flow graph for the second order digital filter with:

a 0 + a 1 z -1 + a 2 z - 2

H(z) = ¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾

1 + b1 z -1 + b2 z-2

Solution: The difference-equation is:

y[n] = a 0 x[n] + a 1 x[n-1] + a 2 x[n-2] - b 1 y[n-1] - b 2 y[n-2]

DEPARTMENT OF Electronics and Communication Engineering

The signal-flow graph in Fig 5.2 is readily deduced from this difference-equation. It is referred
to as a second order or ‘bi-quadratic’ IIR section in ‘direct form 1’.

( )
−1 −2 −3 −4
0.98−3.43 z +4.96 z −3.43 z +0.98 z
H (z)= −1 −2 −3 −4
1−3.47 z +4.96 z −3.39 z +0.96 z
Fig 3.2: “ Direct Form I ” Biquadratic section

Alternative signal flow graphs can generally be found for a given difference-equation.
Considering again the ‘Direct Form I’ bi-quadratic section in Fig 5.2, re-ordering the two halves
in this signal flow graph gives Fig 5.3 which, by Problem 5.9, will have the same impulse-
response as the signal-flow graph in fig 5.2. Now observe that Fig 5.3 may be simplified to the
signal-flow graph in Fig 5.4 which is known as a ‘Direct Form II’ implementation of a bi-
quadratic section. It has the minimum possible number of delay boxes and is said to be
‘canonical’. Its system function is identical to that of the ‘Direct Form I’ signal-flow graph, and
therefore it can implement any second order bi-quadratic system function.

( )
−1 −2 −3 −4
0. 98−3. 43z +4. 96z −3. 43z +0.98z
H ( z)= −1 −2 −3 −4
1−3. 47z +4. 96z −3.39z +0. 96z
DEPARTMENT OF Electronics and Communication Engineering
Fig. 3.3: “ Direct Form I ” rearranged

( )
−1 −2 −3 −4
0.98−3.43z +4.96z −3.43z +0.98z
H (z)= −1 −2 −3 −4
1−3.47z +4.96z −3.39z +0.96z
Fig 3.4: “ Direct Form II ” Biquadratic Section

Example 3.4:

Given values for a 1 , a 2 ,a 0 , b 1 and b 2 , write a program to implement Direct Form II.

Solution: W1 = 0; W2 = 0; {Assign vars W1,W2 to delay outputs}

while 1 {Infinite 'while' loop}

X=input('X = '); {Assign X to receive a single input sample}

W = X-b1*W1-b2*W2; {Recursive part of filter}

Y =W*a0+W1*a1+W2*a2; {Non-recursive part}

W2 =W1; W1 =W; {Set up delay outputs for next time}

disp([' Y=' num2str(Y)]); {Assign var Y to output};

end; {Go back for next sample}

DEPARTMENT OF Electronics and Communication Engineering

  1
3.5. System function: The expression obtained for H(z) is a ratio of polynomials in z .
H(z) is the ‘system function’. When ½½z½½ < 1, H(z) need not be finite.

3.6. Poles and zeros of H(z):

The expression above for H(z) for a general digital filter may be re-expressed as:

(a0z N + a 1 z N -1 + ... + a N )

H(z) = zM-N ¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾

(z M + b 1 z N-1 + ... + b M )

The denominator and numerator polynomials may now be expressed in factorised form to obtain:

(z - z 1 )(z - z 2 )(z - z 3 )...( z - z N )

H(z) = a0zM-N ¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾

(z - p 1 )( z - p 2 )(z - p 3 )...(z - p M )

The roots of the numerator: z 1 , z 2 ,..., z N , are called the ‘zeros’ of H(z).

The roots of the denominator: p 1 ,p 2 ,..., p M , are called the ‘poles’ of H(z)

H(z) will be infinite when evaluated with z equal to a pole, and will become zero with z equal
to a zero except in the special case where the zero coincides exactly with one of the poles.
For a causal stable system, H(z) must be finite for ½½ z ½½ ³³ 1. Therefore there cannot be a
pole whose modulus is greater than or equal to 1. All poles must satisfy ½½ z ½½ < 1, and when
plotted on the Argand diagram, this means that they must lie inside the unit circle. There is no
restriction on the positions of zeros.

DEPARTMENT OF Electronics and Communication Engineering

3.7. Distance rule: Omitted in 2009-10

3.8. Estimation of gain response for Example 5.5 from poles and zeros: Omitted in 2009-10

Example 3.6: Omitted in 2001-10

3.9. Design of a notch filter by MATLAB: Modified in 2009-10

Assume we wish to design a 4th order 'notch' digital filter to eliminate an unwanted sinusoid
at 800 Hz without severely affecting rest of signal. The sampling rate is FS = 10 kHz.

One simple way is to use the MATLAB function ‘butter’ as follows:

FS=10000;

FL = 800 – 25 ; FU = 800+25;

[a b] = butter(2, [FL FU]/(FS/2),’stop’);

a = [0.98 -3.43 4.96 -3.43 0.98]

b= [ 1 -3.47 4.96 -3.39 0.96]

freqz(a, b);

freqz(a, b, 512, FS); % Better graph

axis([0 FS/2 -50 5]); % Scales axes

The frequency-responses (gain and phase) produced by the final two MATLAB statements are as
follows:

DEPARTMENT OF Electronics and Communication Engineering

 0
Magnitude (dB)

-20

-40

0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000
Frequency (Hz)

0
Phase (degrees)

-100

-200

-300

-400
0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000
Frequency (Hz)

Since the Butterworth band-stop filter will have -3dB gain at the two cut-off frequencies

FL = 800-25 and FU=800+25, the notch has ‘-3 dB frequency bandwidth’: 25 + 25 = 50 Hz.

Now consider how to implement the 4th order digital filter. The MATLAB function gave us:

a = [0.98 -3.43 4.96 -3.43 0.98]

b= [ 1 -3.47 4.96 -3.39 0.96]

The transfer (System) Function is, therefore:

H ( z )=
( 0 . 98−3. 43 z−1 +4 . 96 z−2−3. 43 z−3 +0 .98 z− 4
1−3. 47 z −1 +4 . 96 z−2 −3 .39 z−3 +0. 96 z−4 )

DEPARTMENT OF Electronics and Communication Engineering

A ‘Direct Form II’ implementation of the 4 th order notch filter would have the signal-flow graph
below:

DEPARTMENT OF Electronics and Communication Engineering

DEPARTMENT OF Electronics and Communication Engineering
This implementation works fine in MATLAB. But ‘direct form’ IIR implementations of order
greater than two are rarely used. Sensitivity to round-off error in coefficient values will be high.

Also the range of ‘intermediate’ signals in the z-1 boxes will be high.

High word-length floating point arithmetic hides this problem, but in fixed point arithmetic,
great difficulty occurs. Instead we use ‘cascaded bi-quad sections’

Given a 4th order transfer function H(z). Instead of the direct form realization below:

x[n] H(z) y[n]

we prefer to arrange two bi-quad sections, with a single leading multiplier G, as follows:

x[n] G H1(z) H2(z) y[n]

To convert the 4th order transfer function H(z) to this new form is definitely a job for MATLAB.
th
Do it as follows after getting a & b for the 4 order transfer function, H(z), as before:

[a b] = butter(2, [FL FU]/(FS/2),’stop’);

[SOS G] = tf2sos(a,b)

• MATLAB responds with:

SOS = 1 -1.753 1 1 -1.722 0.9776

1 -1.753 1 1 -1.744 0.9785

G = 0.978

DEPARTMENT OF Electronics and Communication Engineering

In MATLAB, ‘SOS’ stands for ‘second order section’ (i.e. bi-quad) and the function ‘tf2SOS’
converts the coefficients in arrays ‘a’ and ‘b’ to the new set of coefficients stored in array ‘SOS’
and the constant G. The array SOS has two rows: one row for the first bi-quad section and one
row for the second bi-quad section. In each row, the first three terms specify the non-recursive
part and the second three terms specify the recursive part. Therefore H(z) may now be realized
as follows:

DEPARTMENT OF Electronics and Communication Engineering

DEPARTMENT OF Electronics and Communication Engineering
This is now a practical and realizable IIR digital ‘notch’ filter, though we sometimes implement
the single multiplier G =0.918 by two multipliers, one for each bi-quad section. More about this
later.

3.10 Calculation of gain-response of notch filter:

How good is a notch filter? We can start to answer this question by specifying the filter's 3dB
bandwidth i.e. the difference between the frequencies where the gain crosses 0.707 (-3dB ). We
should also ask what is the gain at the notch frequency (800 Hz in previous example); i.e. what is
the ‘depth’ of the notch. If it is not deep enough either (i) increase the -3 dB bandwidth or (ii)
increase the order. Do both if necessary. To ‘sharpen’ the notch, decrease the -3dB bandwidth,
but this will make the notch less deep; so it may be necessary to increase the order to maintain a
deep enough notch. This is an ‘ad-hoc’ approach – we can surely develop some theory later. It
modifies the more formal approach, based on poles and zeroes, adopted last year.

DEPARTMENT OF Electronics and Communication Engineering

 00.71 Gain
800- 800 800 + -30 dB
5kHz f

Example: A digital filter with a sampling rate of 200 Hz is required to eliminate an

unwanted 50 Hz sinusoidal component of an input signal without affecting the magnitudes of
other components too severely. Design a 4th order "notch" filter for this purpose whose 3dB
bandwidth is not greater than 3.2 Hz. (MATLAB calls this 2nd order.) How deep is the notch?

Solution:

FS=200; FL=50-1.6; FU=50+1.6;

[a b]=butter(2,[FL,FU]/(FS/2), ‘stop’);

[SOS G] = tf2sos(a,b)

3.11. IIR digital filter design by bilinear transformation

Many design techniques for IIR discrete time filters have adopted ideas and terminology
developed for analogue filters, and are implemented by transforming the system function, H a(s),
of an analogue ‘prototype’ filter into the system function H(z) of a digital filter with similar, but
not identical, characteristics.

For analogue filters, there is a wide variety of techniques for deriving H a(s) to have a specified
type of gain-response. For example, it is possible to deriving H a(s) for an n th order analogue
Butterworth low-pass filter, with gain response:

DEPARTMENT OF Electronics and Communication Engineering

 1
Ga ( ) 
1  ( / C ) 2 n

It is then possible to transform H a(s) into H(z) for an equivalent digital filter. There are many
ways of doing this, the most famous being the ‘bilinear transformation’. It is not the only
possible transformation, but a very useful and reliable one.

The bilinear transformation involves replacing s by (2/T) (z-1)/(z+1)], but fortunately, MATLAB
takes care of all the detail and we can design a Butterworth low pass filter simply by executing
the MATLAB statement:

[a b] = butter(N, fc)

N is the required order and fc is the required ‘3 dB’ cut-off frequency normalised (as usual with
MATLAB) to fS/2. Analogue Butterworth filters have a gain which is zero in the pass-band and
falls to -3 dB at the cut-off frequency. These two properties are preserved by the bilinear
transformation, though the traditional Butterworth shape is changed. The shape change is
caused by a process referred to as ‘frequency warping’. Although the gain-response of the digital
filter is consequently rather different from that of the analogue Butterworth gain response it is
derived from, the term ‘Butterworth filter’ is still applied to the digital filter. The order of H(z)
is equal to the order of Ha(s)

Frequency warping:

It may be shown that the new gain-response G() = Ga() where  = 2 tan(/2). The graph of
 against  below, shows how  in the range - to  is mapped to  in the range - to . The
mapping is reasonably linear for  in the range -2 to 2 (giving  in the range -/2 to /2), but as
 increases beyond this range, a given increase in  produces smaller and smaller increases in
. The effect of frequency warping is well illustrated by considering the analogue gain-response
shown in fig 5.17(a). If this were transformed to the digital filter gain response shown in fig
5.17(b), the latter would become more and more compressed as    .

DEPARTMENT OF Electronics and Communication Engineering

 Fig 3.16 Frequency Warping

‘Prototype’ analogue transfer function: Although the shape changes, we would like G(W) at
its cut off WC to the same as Ga(w) at its cut-off frequency. If Ga(w) is Butterworth, it is -3dB at
its cut-off frequency. So we would like G(W) to be -3 dB at its cut-off WC.

Achieved if the analogue

Fig. 3.17(a): Analogueprototype is designed to
Gain Response have
Fig. its cut-off
3.17(b): frequency
Effect of Bilinear at wC = 2
Transformation
tan(WC/2).

wC is then called the ‘pre-warped’ cut-off frequency.

Designing the analogue prototype with cut-off frequency 2 tan(WC/2) guarantees that the digital
filter will have its cut-off at WC.

Design of a 2nd order IIR low-pass digital filter by the bilinear transform method (‘by
hand’)

DEPARTMENT OF Electronics and Communication Engineering

Let the required cut-off frequency WC = p/4 radians/sample. We need a prototype transfer
function Ha(s) for a 2nd order analogue Butterworth low-pass filter with 3 dB cut-off at C =
2tan(WC/2) = 2 tan(p/8) radians/second. Therefore, wC = 2 tan(p/8) = 0.828. It is well known
by analogue filter designers that the transfer function for a 2nd order Butterworth low-pass filter
with cut-off frequency =1 radian/second is:

1
H a (s )=
1+( √2 )s+s 2

When the cut-off frequency is  = C rather than  = 1, the second order expression for H(s)
becomes:
1
H a (s )=
1+ √ 2( s/ω C )+( s /ωC )2

Replacing s by jw and taking the modulus of this expression gives G() = 1/Ö[1+(w/wC)2n] with
n=2. This is the 2nd order Butterworth low-pass gain-response approximation. Deriving the
above expression for Ha(s), and corresponding expressions for higher orders, is not part of our
syllabus. It will not be necessary since MATLAB will take care of it.

Setting wC = 0.828 in this formula, then replacing s by 2(z-1)/(z+1) gives us H(z) for the
required IIR digital filter. You can check this ‘by hand’, but fortunately MATLAB does all this
for us.

Example 3.7

Using MATLAB, design a second order Butterworth-type IIR low-pass filter with  c =  / 4.

Solution:

[a b] = butter(2, 0.25)

a = [0.098 0.196 0.098]

b = [1 -0.94 0.33]

DEPARTMENT OF Electronics and Communication Engineering

The required expression for H(z) is

0.098 + 0.196 z-1 + 0.098 z-2

H(z) = ¾¾¾¾¾¾¾¾¾¾¾

1 - 0.94 z-1 + 0.33z-2

(
H ( z )=
0.00093+0.0037z−1+.0056z−2+.0037z−3+0.00093z−4
1−2.977z−1+3.422z−2−1.786 z−3+0.3556z−4 )
which may be realised by the signal flow graph in fig 5.18. Note the saving of two multipliers
by using a multipler to scale the input by 0.098.

3.12: Higher order IIR digital filters:

Recursive filters of order greater than two are highly sensitive to quantisation error and
overflow. It is normal, therefore, to design higher order IIR filters as cascades
Fig. 5.18 of bi-quadratic
sections. MATLAB does not do this directly as demonstrated by Example 5.8.

Example 3.8: Design a 4th order Butterworth-type IIR low-pass digital filter is needed with 3dB
cut-off at one sixteenth of the sampling frequency fS.

Solution: Relative cut-off frequency is p/8. The MATLAB command below produces the
arrays a and b with the numerator and denominator coefficients for the 4 th order system function
H(z).

DEPARTMENT OF Electronics and Communication Engineering

 [a b] = butter(4, 0.125)

Output produced by MATLAB is:

a = 0.00093 0.0037 0.0056 0.0037 0.00093

b = 1 -2.9768 3.4223 -1.7861 0.3556

The system function is therefore as follows:

H ( z )=
( 0 . 00093+0 . 0037 z−1 +. 0056 z−2 +. 0037 z−3 +0 . 00093 z−4
1−2. 977 z−1 +3 . 422 z−2−1 . 786 z−3 +0 . 3556 z −4 )

This corresponds to the ‘4th order ‘direct form’ signal flow graph shown below.
+
+ 2.977
-3.422
1.79
-0.356 z-1
z-1
z-1 0.000939
0.0009
+ +

Figure 3.19: A 4th order ‘direct form II’ realisation (not commonly used)

DEPARTMENT OF Electronics and Communication Engineering

Higher order IIR digital filters are generally not implemented like this. Instead, they are
implemented as cascaded biquad or second order sections (SOS). Fortunately MATLAB can
transform the ‘direct form’ coefficients to second order section (SOS) coefficients using a
‘Signal Processing Toolbox’ function ‘tf2sos’ as follows:

[a b] = butter(4, 0.125)

[sos G] = tf2sos(a,b)

Executing these statements gives the following response:

[a b] = butter(4, 0.125)

a = [0.0009 0.0037 0.0056 0.0037 0.0009]

b = [1 -2.9768 3.4223 -1.7861 0.3556 ]

[sos G] = tf2sos(a,b)

sos = [1 2 1 1 -1.365 0.478

1 2 1 1 -1.612 0.745 ]

G = 0.00093

This produces a 2-dimensional array ‘sos’ containing two sets of biquad coefficients and a ‘gain’
constant G. A mathematically correct system function based on this data is as follows:

∞ M
H (e jΩ
) = ∑ − jΩ n
h[ n ] e = ∑ an e− j Ω n
n=−∞ n=0

In practice, especially in fixed point arithmetic, the effect of G is often distributed among the two
sections. Noting that 0.033 x 0.028  0.00093, and noting also that the two sections can be in
either order, an alternative expression for H(z) is as follows:

DEPARTMENT OF Electronics and Communication Engineering

 ∞ M
H (e jΩ
) = ∑ h[ n ] e − jΩ n
= ∑ an e− j Ω n
n=−∞ n=0

This alternative expression for H(z) may be realised in the form of cascaded bi-quadratic
sections as shown in fig 5.20.

∞ M
H(e ) =jΩ
∑ h[n]e −jΩn
= ∑ an e − jΩn

n=−∞ n=0
Fig. 3.20 Fourth order IIR Butterworth LP filter with cut-off fs/16

0.707=1/2 at ‘3dB point’

Fig. 3.21(a) Analogue 4th order Butterworth LP gain response

DEPARTMENT OF Electronics and Communication Engineering

 1/2

0.4

Fig. 3.21(a) shows the 4th order Butterworth low-pass gain response:

∞ M
jΩ
H(e ) = ∑ h[n]e−jΩn
= ∑ an e−jΩn
n=−∞ n=0

(with cut-off frequency normalised to 1) as used by MATLAB as a prototype. Fig 5.21(b) shows
the gain-response of the derived digital filter which, like the analogue filter, is 1 at zero
frequency and 0.707 (-3dB) at the cut-off frequency (/8 0.39 radians/sample). Note however
that the analogue gain approaches 0 as    whereas the gain of the digital filter becomes
exactly zero at  = . The shape of the Butterworth gain response is ‘warped’ by the bilinear
transformation. However, the 3dB point occurs exactly at  c for the digital filter, and the cut-off
rate becomes sharper and sharper as    because of the compression as   .

3.13: IIR digital high-pass band-pass and band-stop filter design:

The bilinear transformation may be applied to analogue system functions which are high-pass,
band-pass or band-stop to obtain digital filter equivalents. For example a ‘high-pass’ digital filter
may be designed as illustrated below:

Example 3.9 Design a 4th order high-pass IIR filter with cut-off frequency fs/16.

Solution: Execute the following MATLAB commands and proceed as for low-pass

[a b] = butter(4,0.125,’high’);

freqz(a,b);

[sos G] = tf2sos(a,b)

DEPARTMENT OF Electronics and Communication Engineering

Wide-band band-pass and band-stop filters (fU >> 2fL) may be designed by cascading low-pass
and high-pass sections, but 'narrow band' band-pass/stop filters (fU not >> 2fL) will not be very
accurate if this cascading approach is used. The MATLAB band-pass approach always works,
i.e. for narrowband and wideband. A possible source of confusion is that specifying an order ‘2’
produces what many people (including me, Barry) would call a 4th order IIR digital filter. The
design process carried out by ‘butter’ involves the design of a low-pass prototype and then
applying a low-pass to band-pass transformation which doubles the complexity. The order
specified is the order of the prototype. So if we specify 2nd order for band-pass we get a 4th
order system function which can be re-expressed (using tf2sos) as TWO biquad sections.

Example 3.10: Design a 2nd (4th)order bandpass filter with F L = p/4 , Fu = p/2.

Solution: Execute the following MATLAB statements:

[a b] = butter(2,[0.25 0.5])

freqz(a,b);

[sos G] = tf2sos(a,b)

MATLAB output:is:

a = 0.098 0 -0.195 0 0.098

b = 1 -1.219 1.333 -0.667 0.33

sos = 1 2 1 1 -0.1665 0.5348

1 -2 1 1 -1.0524 0.6232

G = 0.098

DEPARTMENT OF Electronics and Communication Engineering

 Magnitude (dB) -10

-20

-30

-40

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9

Normalized Frequency ( rad/sample)

200
Phase (degrees)

100

-100

-200
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Normalized Frequency (  rad/sample)

Example 3.11: Design a 4th (8th)order bandpass filter with  L = p/4 ,  u = p/2.

Solution: Execute the following MATLAB statements

[a b] = butter(4,[0.25 0.5])

freqz(a,b); axis([0 1 -40 0]);

[sos G] = tf2sos(a,b)

to obtain the MATLAB output:

a = 0.01 0 -0.041 0 0.061 0 -0.041 0 0.01

b = 1 -2.472 4.309 -4.886 4.477 -2.914 1.519 -0.5 0.12

sos =1 2 1 1 -0.351 0.428

1 -2. 1 1 -0.832 0.49

1 2. 1 1 -0.046 0.724

DEPARTMENT OF Electronics and Communication Engineering

 1 -2 1 1 -1.244 0.793

G = 0.01

4th (8th) order IIR Band-pass (Fs/8 - Fs/4)

0

-10
Magnitude (dB)

-20

-30

-40
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Normalized Frequency ( rad/sample)

0
Phase (degrees)

-200

-400

-600

-800
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Normalized Frequency (  rad/sample)

Example 3.12: Design a 4th (8th)order band-stop filter with F L = p/4 , F u = p/2.

Solution: Execute the following MATLAB statements

[a b] = butter(4,[0.25 0.5], ‘stop’)

freqz(a,b); axis([0 1 -40 0]);

[sos G] = tf2sos(a,b)

to obtain the MATLAB output:

a = 0.347 -1.149 2.815 -4.237 5.1 -4.237 2.815 -1.149 0.347

b = 1 -2.472 4.309 -4.886 4.477 -2.914 1.519 -0.5 0.12

sos = 1 -0.828 1 1 -0.351 0.428

DEPARTMENT OF Electronics and Communication Engineering

 1 -0.828 1 1 -0.832 0.49

1 -0.828 1 1 -0.046 0.724

1 -0.828 1 1 -1.244 0.793

G = 0.347

4th (8th) order IIR Band-stop (Fs/8 - Fs/4)

0
Magnitude (dB)

-10

-20

-30

-40
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Normalized Frequency ( rad/sample)

0
Phase (degrees)

-200

-400

-600

-800
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Normalized Frequency ( rad/sample)

3.14: Comparison of IIR and FIR digital filters:

IIR type digital filters have the advantage of being economical in their use of delays, multipliers
and adders. They have the disadvantage of being sensitive to coefficient round-off inaccuracies
and the effects of overflow in fixed point arithmetic. These effects can lead to instability or
serious distortion. Also, an IIR filter cannot be exactly linear phase.

FIR filters may be realised by non-recursive structures which are simpler and more convenient
for programming especially on devices specifically designed for digital signal processing. These
structures are always stable, and because there is no recursion, round-off and overflow errors are

DEPARTMENT OF Electronics and Communication Engineering

easily controlled. A FIR filter can be exactly linear phase. The main disadvantage of FIR filters
is that large orders can be required to perform fairly simple filtering tasks.

Problems:

3.1 Find H(z) for the following difference equations

(a) y[n] = 2x[n] - 3x[n-1] + 6x[n-4]

(b) y[n] = x[n-1] - y[n-1] - 0.5y[n-2]

3.2 Show that passing any sequence {x[n]} through a system with H(z) = z - 1 produces

{x[n-1]} i.e. all samples are delayed by 1 sampling interval.

3.3. Calculate the impulse-response of the digital filter with

H(z) = ¾¾¾¾¾¾¾¾¾¾

1 - 2 z -1

3.4 Draw the signal-flow graph for example 5.3, and plot its poles and zeros.

3.5 If discrete time LTI systems L1 and L2, with impulse responses {h 1 [n] } and {h 2 [n] }

respectively, are serially cascaded as shown below, calculate the overall impulse

response. Show that this will not be affected by interchanging L 1 & L 2 .

∞ M
H(e jΩ ) = ∑ h[n]e−jΩ n = ∑n
a e− jΩ n

n=−∞ n=0

DEPARTMENT OF Electronics and Communication Engineering

3.6. Design a 4th order band-pass IIR digital filter with lower & upper cut-off
frequencies at 300 Hz & 3400 Hz when fS = 8 kHz.

3.7. Design a 4th order band-pass IIR digital filter with lower & upper cut-off
frequencies at 2000 Hz & 3000 Hz when fS = 8 kHz.

3.8. What limits how good a notch filter we can implement on a fixed point DSP processor?

In theory we can make notch sharper & sharper by reducing the -3dB bandwidth and/or
increasing the order. What limits us in practice.
How sharp a notch can we get in 16-bit fixed pt arithmetic?

3.9. What order of FIR filter would be required to implement a p/4 notch approximately as good

as a 2nd order IIR p/4 notch with 3 dB bandwidth 0.2 radians/sample?

3.10. What order of FIR low-pass filter would be required to be approx as good as the 2nd order

IIR low-pass filter (p/4 cut-off) designed in these notes?

UNIT-IV
Design of DIGITAL FILTERS (FIR) – Structure of FIR Systems:

4.1. Introduction

An FIR digital filter of order M may be implemented by programming the signal-flow-graph

shown below. Its difference equation is:

y[n] = a0x[n] + a1x[n-1] + a2x[n-2] + ... + aMx[n-M]

x[
n]
a
0
z
-
a z
1 -
. z-
1
aM z-
-1 1
a
M
y[n
]
1 1 .
.

Fig. 4.1

DEPARTMENT OF Electronics and Communication Engineering

Its impulse-response is {..., 0, ..., a0, a1, a2,..., aM, 0, ...} and its frequency-response is the DTFT
of the impulse-response, i.e.

∞ M
H (e jΩ ) = ∑ h[ n ]e− jΩ n = ∑ an e− jΩ n
n=−∞ n=0

Now consider the problem of choosing the multiplier coefficients. a 0, a1,..., aM such that H( ej
) is close to some desired or target frequency-response H(ej) say. The inverse DTFT of
H’(ejW) gives the required impulse-response :

The methodology is to use the inverse DTFT to get an impulse-response {h¢[n]} & then realise
some approximation to it Note that the DTFT formula is an integral, it has complex numbers
and the range of integration is from -p to p, so it involves negative frequencies.

Reminders about integration

dx
(1) If x (t )=e at then =aeat
dt

[ ]
π
π π 1 at 1 aπ −aπ
∴ ∫−π x(t )dt = ∫− π e at
dt = e = [ e −e ]
a −π a
b
(2) For any x(t ), ∫a x(t )dt is area under curve

a b t
(Have +ve & -ve areas)

Reminder about complex numbers:

Let x = a + j.b, j = Ö[-1]

Modulus: |x| = Ö[a2 + b2] Arg(x)=tan-1(b/a) + {p.sign(b) if a < 0}

= tan2(b,a) = angle(a + j.b) : range - p to pPolar: x = Rejf

where R = |x| & f = Arg(x)

DEPARTMENT OF Electronics and Communication Engineering

De Moivre: ejf = cos(f) + j.sin(f) , e-jf = cos(f) - j.sin(f)

\ ejf + e-jf = 2cos(f) & ejf - e-jf = 2 j.sin(f)

Complex conjugate: x* = a - j.b = Re-jf.

What about the negative frequencies?

Examine the DTFT formula for H(ejW).

∞ ∞
H (e )= ∑ h[ n]e
jΩ − jΩ n
∴ H (e − jΩ
)= ∑ h[ n]e j Ωn
n=−∞ n=−∞

If h[n] real then h[n]ejW is complex-conjugate of h[n]e-jW. Adding up terms gives H(e-jW ) as
complex conj of H(ejW).

G(W) = G(-W) since G(W) = |H(ejW)| & G(-W) = H(e-jW)|

(Mod of a complex no. is Mod of its complex conj.)

Because of the range of integration (- to ) of the DTFT formula, it is common to plot graphs of
G() and () over the frequency range - to  rather than 0 to . As G(W) = G(-W) for a real
filter the gain-response will always be symmetric about W=0
4.2. Design of an FIR low-pass digital filter
Assume we require a low-pass filter whose gain-response approximates the ideal 'brick-wall'
- in Figure 4.2. 1/3
gain-response 0G() /3

Fig. 4.2
If we take the phase-response f() to be zero for all , the required frequency-response is:-

and by the inverse DTFT,

DEPARTMENT OF Electronics and Communication Engineering

 = (1/3)sinc(n/3) for all n. where

A graph Mainagainst
-4 of sinc(x)
-3 ‘lobe’-2x is shown-1below:
1 sinc(x)
1
‘Zero-crossings’ at 2x =1,‘Ripples’
2, 3,3 etc. x

Fig 4.3a

The ideal impulse-response {h[n]} with each sample h[n] = (1/3)sinc(n/3) is therefore as
follows:

-12 -9 -6 -3 1/3 h[n] 3 6 Fig.

9 4.3b n
Ideal impulse response for low-pass filter cut-off 

Reading from the graph, or evaluating the formula, we get:

{h¢[n]} = { ..., -0.055, -0.07, 0, 0.14, 0.28, 0.33, 0.28, 0.14, 0, -0.07, -0.055, ... }

A digital filter with this impulse-response would have exactly the ideal frequency-response we
applied to the inverse-DTFT i.e. a ‘brick-wall’ low-pass gain response & phase = 0 for all W.
But {h¢[n]} has non-zero samples extending from n = -¥ to ¥, It is not a finite impulse-
response. It is also not causal since h[n] is not zero for all n < 0. It is therefore not realisable
in practice.

DEPARTMENT OF Electronics and Communication Engineering

To produce a realisable impulse-response of even order M:

{
' −M M
h [n] : ≤ n ≤
(1) Set h[ n ] = 2 2
0 : otherwise
(2) Delay resulting sequence by M/2 samples to ensure that the first non-zero sample occurs at
n = 0.

The resulting causal impulse response may be realised by setting a n = h[n] for n=0,1,2,...,M.
Taking M=4, for example, the finite impulse response obtained for the p/3 cut-off low-pass
specification is :

{..,0,..,0, 0.14, 0.28, 0.33 , 0.28 , 0.14 , 0 ,..,0,..}

The resulting FIR filter is as shown in Figure 4.1 with a 0=0.14, a1=0.28, a2=0.33, a3=0.28,
a4=0.14. ( Note: a 4th order FIR filter has 4 delays & 5 multiplier coefficients ).

The gain & phase responses of this FIR filter are sketched below.
dB
-10
0
-30
-20 /3
G(
-6 dB    
  

Fig. 4.4

Clearly, the effect of the truncation of {h[n]} to ±M/2 and the M/2 samples delay is to produce
gain and phase responses which are different from those originally specified.Considering the
gain-response first, the cut-off rate is by no means sharp, and two ‘ripples’ appear in the stop-
band, the peak of the first one being at about -21dB.

The phase-response is not zero for all values of  as was originally specified, but is linear phase
( i.e. a straight line graph through the origin ) in the pass-band of the low-pass filter ( - p/3 to p/3
) with slope arctan( M/2 ) with M = 4 in this case. This means that f(  ) =  ( M/2 ) for |  |
£ p/3; i.e. we get a linear phase-response ( for |  | £ p/3 ) with a phase-delay of M/2
samples.It may be shown that the phase-response is linear phase because the truncation was done
symmetrically about n=0.Now let’s try to improve the low-pass filter by increasing the order to
ten. Taking 11 terms of { (1 / 3) sinc (n / 3) } we get, after delaying by 5 samples:

DEPARTMENT OF Electronics and Communication Engineering

 {...0,-0.055,-.069, 0,.138,.276,.333,.276,.138,0,-.069,-.055,0,...}.The signal-flow
graph of the resulting 10th order FIR filter is shown below:
z-1 + -.07
x[n] -.055 z-1 + 0 z-1 + 0.14 z-1 + .28z-1 + .33z-1 + .28z-1 + 0.14
z-1 + 0 z-1 + -.07
z-1 + -.055
y[n]

Notice that the coefficients are again symmetric about the centre one (of value 0.33) and this
again ensures that the FIR filter is linear phase. The gain and phase responses of this tenth order
FIR filter are produced by the MATLAB statement:

freqz( [-0.055, -0.069, 0, 0.138, 0.276, 0.333, 0.276, 0.138, 0, -0.069, -0.055] );

In may be seen in the gain-response, as reproduced below, that the cut-off rate for the 10 th order
FIR filter is sharper than for the 4th order case, there are more stop-band ripples and, rather
disappointingly, the gain at the peak of the first ripple after the cut-off remains at about -21 dB.
This effect is due to a well known property of Fourier series approximations, known as Gibb's
phenomenon. The phase-response is linear phase in the passband ( -p/3 to p/3 ) with a phase
delay of 5 samples. As seen in fig 4.6, going to 20th order produces even faster cut-off rates and
more stop-band ripples, but the main stop-band ripple remains at about -21dB. This trend
continues with 40th and higher orders as may be easily verified. To improve matters we need to
discuss ‘windowing’.

DEPARTMENT OF Electronics and Communication Engineering

Fig 4.5: Gain response of tenth order lowpass FIR filter with C = /3

Fig 4.6: Gain response of 20th order lowpass FIR filter with C = /3

4.3. Windowing:

To design the FIR filters considered up to now we effectively multiplied {h[n]} as calculated by
the inverse DTFT by a rectangular window sequence {rM+1[n]} where

This causes a sudden transition to zero at the window edges and it is these transitions that
produce the stop-band ripples. To understand why, we need to know that the DTFT of {r M+1[n]}
is as follows:

DEPARTMENT OF Electronics and Communication Engineering

 {
sin(( M +1 )Ω/2 )
jΩ : Ω≠0 ,±2 π , .. .
∴ R M+1 ( e )= sin(Ω/2 )
M +1 : otherwise
Note that {rM+1[n]} has non-zero values at M+1 values of n which include n = 0. A graph of
RM+1(ej) against  for M=20 is shown below. It is purely real and this is because because
rM+1(n) is symmetric about n=0.
Dirichelet K of order 20

20

15

10
R

-5
-3 -2 -1 0 1 2 3
radians/sample

It looks rather like a ‘sinc’ function in that it has a ‘main lobe’ and ‘ripples’. It is a little
different from a sinc in that the ripples do not die away quite as fast. It is these ripples that cause
the ripples in the stop-bands of FIR digital filters. It may be shown that:
(a) the height of RM+1(ej) is M+1
(b) the area under the main lobe remains approximately 1 for all values of M.
Frequency-domain convolution:
Multiplying two sequences {x[n]} & {[y[n]} to produce {x[n].y[n]} is called ‘time-domain
multiplication’.

It may be shown that if {x[n]} has DTFT X(e jW) & [y[n]} has DTFT Y(ejW) then the DTFT of
{x[n].y[n]} is:

π
 ( 1/2 π )∫−π X ( e jθ )  R ( e j (Ω−θ ) )  dθ
This is the ‘frequency-domain convolution’ of X(ejW) with R(ejW).
Therefore: ‘time-domain multiplication’ is equivalent to ‘frequency-domain convolution’.

DEPARTMENT OF Electronics and Communication Engineering

Windowing is time-domain multiplication.
Now apply the above formula to explain the effect of rectangular windowing on an ideal gain-
response. Multiplying an ideal impulse-response {h¢[n]} by the rectangular window {rM+1[n]}

causes H¢(ejW) to be ‘convolved’ with RM+1(ejW):

π
H (e )= ( 1/2π )∫− π H ' ( e jθ )  R M+1 ( e j ( Ω−θ ) )  dθ
jΩ

If H¢(ejW) has an ideal ‘brick-wall’ gain-response, & RM+1(ejW) is as shown in previous graph,

convolving H¢(ejW) with RM+1(ejW) reduces the cut-off rate and introduces stop-band ripples.

The sharper the main-lobe of RM+1(ejW) and the lower the ripples, the better.

To illustrate for an ideal low-pass frequency response with cut-off /3 :

H (e j Ω )=
'
{10 : −π: /3≤Ω≤π /3
otherwise

{
w M+1[n]= 0.54+0.46cos(πn/M) :-M/2 ≤ n ≤M/2
0 : |n| > M/2
For any given frequency , this is 1/(2) times the area under the curve as from /3 to
+/3. Consider what happens as  increases from 0 towards ? When  = 0, the whole area
of the main lobe of RM+1(ej ) is included in the integral and the gain G() is close to 1 (= 0 dB).
When  = /3 the gain drops to approximately 0.5 (= -6 dB) since only half the area of the main
lobe is included. When  is further increased above /3, the area of the main lobe becomes less
and less included. So the gain becomes very small as only the ripples are being integreated.
Because of the ripples, for some values of  in the stop-band, some the positive and negative
areas will cancel out, giving a gain of zero ( dB). For other values of  in the stop-band, there
will be slightly more positive area than negative area, or vice–versa. Thus the stop-band ripples
are created.
So we conclude that ripples in the gain-responses of FIR digital filters designed with rectangular
windows arise from the frequency-domain convolution between the ideal (target) frequency
response and the DTFT RM+1(ejW) of the rectangular window. The gain at the cut-off
frequency will be approximately –6 dB less than the pass-band gain (normally 0 dB) because
only half the main lobe lies within the ideal filter’s pass-band. Ripples also occur in the pass-
band but we can hardly notice them.

DEPARTMENT OF Electronics and Communication Engineering

4.4. Non-rectangular windows
The levels of the ripples may be reduced if {r M+1[n]} is replaced by a non-rectangular window
sequence, denoted { wM+1[n] } say, which produces a more gradual transition at the window
edges. The simplest non-rectangular window sequence in common use is the Hann or von-Hann
window which is essentially a ‘raised cosine’:

Many other types of window sequence exist (e.g. Hamming, Kaiser ). Also, slightly different
formulae exist for the Hann window. The formula used by MATLAB is different so stick to
mine for the moment. Perhaps the most well known window is the Hamming window whose
formula is very similar to that of the Hann window. It is also a 'raised cosine' window and has a
similar effect except that it is generally considered to be slightly better. The (M+1) th order
Hamming window formula is:

{
w M +1 [ n ]= 0 . 54+0 . 46 cos( πn/ M ) :-M/2 ≤ n ≤M/2
0 : |n| > M/2

Note that this has non-zero values at M+1 values of n including n=0.

Effect of non-rectangular windows:

For any of these windows, multiplying {h[n]} by {wM+1[n]} instead of {rM+1[n]} gradually tapers
the impulse-response towards zero at the window edges. To understand why this reduces stop-
band ripples, compare the DTFT of {w M+1[n]} with the DTFT of {rM+1[n]}. Take a Hann window
with M=20.

DEPARTMENT OF Electronics and Communication Engineering

 DTFT of Hann & Rec t windows , order 20

20 Hann
Rec t

15

10
R

-5
-3 -2 -1 0 1 2 3
radians /s am ple

To draw this graph it was necessary to derive the DTFT of {w M[n]} for a Hann window. Don’t
worry about this for now. Just look at the graph.

Comparing the DTFT of a rectangular and a Hann window of the same order:
(a) The height of WM+1(ej) is (M+1)/2 which is half that of RM+1(ej)

(b) The main lobe of WM+1(ejW) is about twice the width of the main lobe of RM+1(ejW).

(c) It may be shown that the area of the main lobe of WM+1(ej) remains about 1.

(d) The ripples in WM+1(ejW) are greatly reduced in comparison to R M+1(ejW). This is
good!

The main lobe of WM+1(ejW), being less sharp than that of RM+1(ejW), will reduce the
sharpness of the cut-off of the filter designed with the Hann window. This is the price to be paid
for reducing stop-band ripples. Note that the ripples have not been eliminated entirely.

Applying a 4th order Hann window

Consider again the low-pass filter example with cut-off frequency p/3 radians per sampling
interval. The ideal impulse-response was found to be:
{h[n]} = {... , 0.14, 0.28, 0.33, .28, .14, ...}

When M = 4, the Hann window {wM+1[n]} = {..,0,..,0, 0.25, 0.75, 1, .75, .25, 0,..,0,..}

DEPARTMENT OF Electronics and Communication Engineering

Multiplying term by term and delaying by M / 2 = 2 samples we obtain the finite impulse-
response:

{..,0,..,0, .04, .21,.33,.21,.04, 0,..,0,..}

The resulting "Hann-windowed" FIR filter of order 4 is as shown in Figure 4.1 with a 0=0.04,
a1=0.21, etc. Its gain-response is as shown in Figure 4.9.

Figure 4.9: 4th order FIR filter (Hann)

The Hann window gradually reduces the amplitude of the ideal impulse-response towards zero at
the edges of the window rather than truncating it as does the rectangular window. As seen in the
graph above, the effect on the gain-response of the FIR filter obtained is:
i) to greatly reduce stop-band ripples ( good ).

ii) to reduce the cut-off rate ( bad ).

The phase-response is not affected in the passband. We can improve the cut-off rate by going to
higher orders. The graphs below are for 10th and 20th order ( Hann windowed ):

DEPARTMENT OF Electronics and Communication Engineering

Fig 4.10: Tenth order FIR filter with C = /3 ( Hann window )

MATLAB program to design & graph 10th order FIR lowpass filter with Hann window:

M=10; Fs=8000;

for n=-M/2 : M/2

w=0.5*(1+cos(pi*n/(1+M/2)));

a(1+M/2+n) = 0.3333*sinc(n/3)*w;

end;

freqz(a,1,1000,Fs);

DEPARTMENT OF Electronics and Communication Engineering

Fig 4.11: 20th order FIR filter with C = /3 (Hann window)

Use of MATLAB function ‘FIR’

‘fir1’ uses ‘windowing method’ we have just seen. By default it uses a Hamming window which
is very similar to the Hann window.

It also scales the coefficients to make the gain exactly 0 dB at 0 Hz

To design 10th order FIR lowpass filter with Hamming window & cut-off frequency p/3 (º
Fs/6):

c=fir1(10, 0.33);

Effect of Hann window on first stop-band ripple

It may be observed that the gain at the peak of the first stop-band ripple has been reduced from
about
-21 dB to about –44 dB. At the frequency of the first stop-band ripple since 20log 10(1/10) =-20
dB, the amplitude of the signal being filtered is reduced by factor »10 when the FIR filter is
designed with a rectangular window.

Since 20 log10(1/100) = -40, the same amplitude will be reduced by factor >100 when the FIR
filter is designed with a Hann window.

DEPARTMENT OF Electronics and Communication Engineering

 x(t) = 5 sin(wt) becomes y(t) = 0.5 sin(wt) with rectangular or y(t) = 0.05 sin(wt) with
Hann.

If this is not good enough we must use a different window (e.g. Hamming, Kaiser)

Hamming window is similar to Hann but slightly better.

4.5. Kaiser window

The Kaiser window offers a range of options from rectangular, through Hamming towards even
lower stop-band ripples.

Ripple reduction is at the expense of a less sharp cut-off rate.

The formula is too complicated to quote, but MATLAB command:

KW = kaiser(M+1,beta)

produces a Kaiser window array of length M+1 for any value of beta > 0.

When beta (b) = 0, this is a rectangular window, and when beta = 5.4414 we get a Hamming
window.

Increasing beta further gives further reduced stop-band ripples with a reduced cut-off sharpness.

A slight complication is that MATLAB arrays must start at index 1, rather than –M/2. So we
have to add 1+M/2 to each value of n to find the right entry in array KW. It may be shown that:

(a) the area of the main-lobe remains about 1

(b) the height of the main lobe depends on M and beta.

See Slides 84-90 for gain responses of FIR low-pass filters designed using Kaiser windows with
different values of beta. These are compared with filters designed with Hamming windows.

4.6. High-pass, band-pass & band-stop linear phase FIR filters

DEPARTMENT OF Electronics and Communication Engineering

These can be designed almost as easily as low-pass as long as you remember to define the
required gain-response G() correctly from - to +; i.e. make G(-) = G(). For example, a
band-pass filter with pass-band from FS/8 to FS/4 would have the following gain response
ideally:-

Fig 4.12: Gain response of ideal band-pass filter

Taking () = 0 for all  initially as before, we obtain:

 0 : ||   /4
j 
H (e )  1 :  / 4  |  |   / 2
 0 :  /2  ||  


and applying the inverse DTFT gives the following formula (not forgetting negative values of :

1  1  /4 1  / 2  jn
h[n] 
2 

H (e j )e  jn d 
2 
 /2
1 e  jn d 
2 
 /4
1 e d

Evaluating the integrals in this expression will produce a nice formula for h[n].

Truncating symmetrically, windowing and delaying this sequence, as for the low-pass case
earlier, produces the causal impulse-response of an FIR filter of the required order. This will be
linear phase for the same reasons as in the low-pass case. Its impulse-response will be
symmetric about n=M/2 when M is the order. MATLAB is able to design high-pass, band-pass
and band-stop linear phase FIR digital filters by this method.

FIR filter design in MATLAB:

c = fir1(10,0.33,'high') designs a 10th order high-pass filter.

c = fir1(10,[0.2 0.4],'bandpass') designs a 10th order band-pass filter with cut-off frequencies
0.2p and 0.4p.

c = fir1(20,[0.2 0.4],'stop') designs a 20th order band-stop filter with cut-off frequencies 0.2p
and 0.4p.

By default ‘fir1’ uses a ‘Hamming’ window (very similar to a Hann) and scales the pass-band
gain to 0 dB. Linear phase response is obtained.

DEPARTMENT OF Electronics and Communication Engineering

4.7. Summary of ‘windowing’ design technique for FIR filters

To design an FIR digital filter of even order M, with gain response G() and linear phase, by
the windowing method,

1) Set H(ej) = G() the required gain-response. This assumes f() = 0.

2) Inverse DTFT to produce the ideal impulse-response {h[n]}.

3) Window to ±M/2 using chosen window.

4) Delay windowed impulse-response by M/2 samples.

5) Realise by setting multipliers of FIR filter.

Instead of obtaining H( ej ) = G(  ), we get e-jM/2G() with G() a distorted version of
G() the distortion being due to windowing.

The phase-response is therefore () = -M/2 which is a linear phase-response with phase-
delay M/2 samples at all frequencies  in the range 0 to . This is because -() /  = M/2 for
all .

Notice that the filter coefficients, and hence the impulse-response of each of the digital filters we
have designed so far are symmetric in that h[n] = h[M-n] for all n in the range 0 to M where M is
the order. If M is even, this means that h[M/2 - n] = h[M/2 + n] for all n in the range 0 to M/2.
The impulse response is then said to be 'symmetric' about sample M/2. The following example
illustrates this for an example where M=6 and there are seven non-zero samples within {h[n]}:
{… 0, …, 0, 2, -3, 5, 7, 5, -3, 2, 0, …,0, … }
The most usual case is where M is even, but, for completeness, we should briefly consider the
case where M is odd. In this case, we can still say that {h[n]} is 'symmetric about M/2' even
though sample M/2 does not exist. The following example illustrates the point for an example
where M=5 and {h[n]} therefore has six non-zero sample:
(…, 0,…, 0, 1, 3, 5, 5, 3, 1, 0, …, 0, …}
When M is odd, h[(M-1)/2 - n] = h[(M+1)/2 + n] for n = 0, 1, …, (M-1)/2.
It may be shown that FIR digital filters whose impulse-responses are symmetric in this way are
linear phase. We can easily illustrate this for either of the two examples just given. Take the
second. Its frequency-response is the DTFT of {h[n]} i.e.

It is also possible to design FIR filters which are not linear phase. The technique described in
this section is known as the ‘windowing’ technique or the ‘Fourier series approximation
technique’.

DEPARTMENT OF Electronics and Communication Engineering

4.8. Designing & implementing FIR digital filters using the Mat lab Signal Processing
Toolbox

11) The design of linear phase FIR digital filters by the windowing technique discussed
above is readily carried out by the command 'FIR1' provided by the 'signal processing toolbox' in
MATLAB. The filter may be applied to a segment of sampled sound stored in an array by the
command 'filter'. To illustrate the use of these commands we now design and implement a
128th order FIR band-pass digital filter with lower and upper cut-off frequencies 300 Hz and 3.4
k Hz respectively and apply it to a wav file containing mono music sampled at 11.025 kHz.
(This makes the music sound like it has been transmitted over a wired telephone line.)
12) Notes:
13) (1) The FIR cut-off frequencies must be specified relative to fS/2 rather than in
radians/sample. This is achieved by dividing each cut-off frequency in radians/sample by .
14) (2) By default FIR1 uses a Hamming window. Other available windows such as Hann
can be specified with an optional trailing argument.
15) (3) By default, the filter is scaled so the center of the first pass-band has magnitude
exactly one after windowing.
16)
clear all;

[x, Fs, nbits] = wavread('caprice.wav');

17) wlow = 2 * pi * 300 / Fs ; % radians/sample rel to fs

18) wup = 2 * pi * 3400 / Fs ; % radians/sample rel to fs
a = fir1(128, [wlow wup] / pi ) ;

freqz ( a , 1,1000,Fs) ; % plot gain & phase

19) y = filter(a, 1, x );
wavwrite(x,fs,nbits,'capnew.wav');

4.9. Remez Exchange Algorithm method:

An FIR digital filter design technique which is better than the windowing technique, but more
complicated, is known as the ‘Remez exchange algorithm’. It was developed by McClelland and
Parks and is available in MATLAB. The following MATLAB program designs a 40 th order FIR
low-pass filter whose gain is specified to be unity ( i.e. 0 dB ) in the range 0 to 0.3 
radians/sample and zero in the range 0.4 to . The gain in the “ transition band ” between 0.3
and 0.4 is not specified. The 41 coefficients will be found in array ‘a’. Notice that, in contrast
to the gain-responses obtained from the 'windowing' technique, the Remez exchange algorithm
produces 'equi-ripple' gain-responses (fig 4.14) where the peaks of the stop-band ripples are

DEPARTMENT OF Electronics and Communication Engineering

equal rather than decreasing with increasing frequency. The highest peak in the stop-band will
be lower than that of an FIR filter of the same order designed by the windowing technique to
have the same cut-off rate. Although they are a little difficult to see, there are 'equi-ripple' pass-
band ripples.

a = remez (40, [0, 0.3, 0.4,1],[1, 1, 0, 0] );

freqz (a,1,1000,Fs);

Fig 4.14: Gain response of 40th order FIR lowpass filter designed by “ Remez ”

4.10. Fixed point implementation of FIR digital filters

MATLAB programs were presented in Section 3 to illustrate how an FIR digital filter as
designed above may be implemented on a computer or microprocessor. These programs made
full use of the powerful and highly accurate floating point arithmetic operations available in
MATLAB. However FIR digital filters are often implemented in mobile battery powered
equipment such as a mobile phone where a floating point processor would be impractical as it
would consume too much power. Low power fixed point DSP processors are the norm for such
equipment, typically with a basic 16-bit word-length. Such processors must be programmed
essentially using only integer arithmetic, and it is interesting now to consider how an FIR filter
could be programmed in this way.

DEPARTMENT OF Electronics and Communication Engineering

Take as a simple example the 4th order FIR low-pass digital filter designed earlier with impulse
response: {..,0,..,0, .04, .21,.33,.21,.04, 0,..,0,..}. Rounding each coefficient to the nearest
integer would clearly be a mistake because they would all become trivially zero. The solution is
simple: multiply all coefficients by a constant that is large enough such that when we round to
integers, the resulting error is not too large. So we could make the FIR filter coefficients A 0= 4,
A1=21, A2=33, A3=21 and A4=4. To compensate for this action, we must divide the output
produced by the same constant, in this case 100. In practice instead of 100 or maybe 1000, we
choose a power of two for the constant, since dividing by a power of two (e.g. 1024) is very
simple, requiring only an arithmetic 'shift'. Dividing by 1024 requires only 10 bit-shifts right.

The larger the constant, the less the effect or rounding and the more accurate the coefficients.
However we must be careful not to choose too large a constant because the number of bits per
word is limited often to 16. If the integers produced during the calculation get too large, we risk
overflow and very serious non-linear distortion. Where the result of an addition of positive
integers is too large for the available 16-bits, it may become, by 'wrap-around', a large negative
number which may cause very serious distortion. Similarly the addition of two negative
numbers could wrap around to a large positive number if overflow occurs. So here we have an
difficult balancing act to perform between coefficient inaccuracy and overflow. Modern fixed
point DSP processors offer useful facilities and extra bits to help with the design of fixed point
DSP programs, but this remains a difficult area in general.

Fortunately, FIR digital filters are particularly easy to program in fixed point arithmetic and this
is one of their main advantages. For one thing, unlike IIR digital filters, they can never become
unstable as there is no feedback. The effects of rounding and overflow can be a nightmare when
they are fed back recusively as with IIR filters.

In some cases, overflows (with wrap-around) can be allowed to occur repeatedly with an FIR
filter with the sure knowledge that they will eventually be cancelled out by corresponding wrap-
around overflows in the opposite sense. So we can generally risk overflow more readily with
FIR digital filters than with IIR digital filters, and thus have greater coefficient accuracy.

Fixed point DSP programmers often describe the coefficient scaling proceess used above, i.e.
multiplying by constants which are powers of two, in a slightly different way. Scaling by 1024,
is adopting a 'Q-format' of ten and the programmer is effectively assuming a decimal point (or
binary point) to exist ten bit positions from the right within the 16-bit word. The decimal point
is 'invisible' and often nothing is stored, apart from judicious comments. The programmer
him/herself must keep track of the Q-formats used in different parts of complex programs and
make sure the correct compensation (e.g. dividing by 1024) is applied when needed to produce a
correct output. The effects of limited word-length fixed point arithmetic may be studied in
MATLAB by restricting programs to use integers and integer operations only. A MATLAB

DEPARTMENT OF Electronics and Communication Engineering

implementation of the 4th order low-pass filter mentioned above using integer arithmetic only is
now given.

A = [4 21 33 21 4 ] ;

x = [0 0 0 0 0 ] ;

while 1

x(1) = input( 'X = ');

Y = A(1)*x(1);

for k = 5 : -1: 2

Y = Y + A(k)*x(k);

x(k) = x(k-1);

end;

Y = round( Y / 100) ;

disp([' Y = ' num2str(Y)]);

end;

4.11. Advantages of FIR filters compared with IIR

FIR filters are easy to program in fixed point arithmetic. They never become unstable as there is
no feedback. They can be exactly linear phase.

In some cases, overflows can be allowed to occur since if gain is never greater than 1, you know
that a positive overflow will eventually be cancelled out by a negative overflow or vice versa.•
This works if you do not use ‘saturation mode’ arithmetic which avoids ‘wrap-round’. Can
risk overflow more readily with FIR digital filters than with IIR digital filters, & thus have
greater coefficient accuracy.

Disadvantage: FIR need higher orders than IIR (later).

DEPARTMENT OF Electronics and Communication Engineering

DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
PROBLEMS:

4.1 Design a 10th order FIR low-pass digital filter with cut-off at Fs/4 frequency with & without
a Hann window. Use MATLAB to compare the gain-responses obtained.
4.2 Design a tenth order FIR bandpass digital filter with lower and upper cut-off frequencies at
p/8 and p/3 respectively.

4.3 Write a MATLAB program for one of these filters using integer arithmetic only.

4.4 Design a 4th order FIR high-pass filter with cut-off frequency at p/3.

4.5 Do all FIR filters have exactly linear phase responses?

4.6 Given that for a linear phase filter, its impulse-response must be symmetric about n=N for
some N, (i.e. h[N + n] = h[N-n] for all n), why cannot an IIR filter be linear phase?

4.7. Design a sixth order linear phase FIR filter whose gain-response approximates that
shown in fig 4.13. Plot its gain-response using MATLAB.

4.8. Show that a linear phase lead of ff( WW ) = -kWW corresponds to a delay of k samples.

Answer to 4.6: For linear phase, impulse-response must be symmetric about some value of n, say
n=M. If it is an IIR it goes on for ever as n ® ¥. So it must go on for ever backwards as n ® -
¥. Would have to be non-zero for values on n<0; i.e. non-causal

If h[n] symmetric & IIR it must

h[n]
be non-causal
M n

DEPARTMENT OF Electronics and Communication Engineering

DEPARTMENT OF Electronics and Communication Engineering
UNIT V – Introduction to Multirate Digital Signal Processing:
1.1 What does multirate mean?
Multirate simply means “multiple sampling rates”. A multirate DSP system uses multiple sampling rates
within the system. Whenever a signal at one rate has to be used by a system that expects a different rate, the
rate has to be increased or decreased, and some processing is required to do so. Therefore “Multirate DSP”
really refers to the art or science of changingsampling rates.

1.2 Why should I do multirate DSP?

The most immediate reason is when you need to pass data between two systems which use incompatible
sampling rates. For example, professional audio systems use 48 kHz rate, but consumer CD players use 44.1
kHz; when audio professionals transfer their recorded music to CDs, they need to do a rate conversion.

But the most common reason is that multirate DSP can greatly increase processing efficiency (even by
orders of magnitude!), which reduces DSP system cost. This makes the subject of multirate DSP vital to all
professional DSP practitioners.

1.3 What are the categories of multirate?

Multirate consists of:

1. Resampling:To combine decimation and interpolation in order to change the sampling rate by a
fractional value that can be expressed as a ratio. For example, to resample by a factor of 1.5, you just
interpolate by a factor of 3 then decimate by a factor of 2 (to change the sampling rate by a factor of
3/2=1.5.)

2. Decimation: To decrease the sampling rate,

3. Interpolation: To increase the sampling rate,

Decimation

2.1 Basics
2.1.1 What are “decimation” and “downsampling”?
Loosely speaking, “decimation” is the process of reducing the sampling rate. In practice, this usually implies
lowpass-filtering a signal, then throwing away some of its samples.

“Downsampling” is a more specific term which refers to just the process of throwing away samples, without
the lowpass filtering operation. Throughout this FAQ, though, we’ll just use the term “decimation” loosely,
sometimes to mean “downsampling”.

2.1.2 What is the “decimation factor”?

The decimation factor is simply the ratio of the input rate to the output rate. It is usually symbolized by “M”,
so input rate / output rate=M.

Tip: You can remember that “M” is the symbol for decimation factor by thinking of “deci-M-ation”.
(Exercise for the student: which letter is used as the symbol for interpo-L-ation factor?)

2.1.3 Why decimate?

The most immediate reason to decimate is simply to reduce the sampling rate at the output of one system so
a system operating at a lower sampling rate can input the signal. But a much more common motivation for

DEPARTMENT OF Electronics and Communication Engineering

decimation is to reduce the cost of processing: the calculation and/or memory required to implement a DSP
system generally is proportional to the sampling rate, so the use of a lower sampling rate usually results in a
cheaper implementation.

To that, Jim Thomas adds:

Almost anything you do to/with the signal can be done with fewer operations at a lower sample rate, and the
workload is almost always reduced by more than a factor of M.
For example, if you double the sample rate, an equivalent filter will require four times as many operations to
implement. This is because both amount of data (per second) and the length of the filter increase by two, so
convolution goes up by four. Thus, if you can halve the sample rate, you can decrease the work load by a
factor of four. I guess you could say that if you reduce the sample rate by M, the workload for a filter goes
down to (1/M)^2.

2.1.4 Is there a restriction on decimation factors I can use?

Yes. Decimation involves throwing away samples, so you can only decimate by integer factors; you cannot
decimate by fractional factors. (However, you can do interpolation prior to decimation to achieve an overall
rational factor, for example, “4/5”; see Part 4: Resampling.)

2.1.5 Which signals can be downsampled?

A signal can be downsampled (without doing any filtering) whenever it is “oversampled”, that is, when a
sampling rate was used that was greater than the Nyquist criteria required. Specifically, the signal’s highest
frequency must be less than half the post-decimation sampling rate. (This just boils down to applying the
Nyquist criteria to the input signal, relative to the new sampling rate.)

In most cases, though, you’ll end up lowpass-filtering your signal prior to downsampling, in order to
enforce the Nyquist criteria at the post-decimation rate. For example, suppose you have a signal sampled at
a rate of 30 kHz, whose highest frequency component is 10 kHz (which is less than the Nyquist frequency
of 15 kHz). If you wish to reduce the sampling rate by a factor of three to 10 kHz, you must ensure that you
have no components greater than 5 kHz, which is the Nyquist frequency for the reduced rate. However,
since the original signal has components up to 10 kHz, you must lowpass-filter the signal prior to
downsampling to remove all components above 5 kHz so that no aliasing will occur when downsampling.

This combined operation of filtering and downsampling is called decimation.

2.1.6 What happens if I violate the Nyquist criteria in downsampling or decimating?

You get aliasing–just as with other cases of violating the Nyquist criteria. (Aliasing is a type of distortion
which cannot be corrected once it occurs.)

2.2 Multistage
2.2.1 Can I decimate in multiple stages?
Yes, so long as the decimation factor, M, is not a prime number. For example, to decimate by a factor of 15,
you could decimate by 5, then decimate by 3. The more prime factors M has, the more choices you have.
For example you could decimate by a factor of 24 using:

 one stage: 24
 two stages: 6 and 4, or 8 and 3
 three stages: 4, 3, and 2
 four stages: 3, 2, 2, and 2

2.2.2 Cool. But why bother with all that?

DEPARTMENT OF Electronics and Communication Engineering

If you are simply downsampling (that is, throwing away samples without filtering), there’s no benefit. But
in the more common case of decimating (combining filtering and downsampling), the computational and
memory requirements of the filters can usually be reduced by using multiple stages.

2.2.3 OK, so how do I figure out the optimum number of stages, and the decimation factor at each stage?
That’s a tough one. There isn’t a simple answer to this one: the answer varies depending on many things, so
if you really want to find the optimum, you have to evaluate the resource requirements of each possibility.

However, here are a couple of rules of thumb which may help narrow down the choices:

 Using two or three stages is usually optimal or near-optimal.

 Decimate in order from the largest to smallest factor. In other words, use the largest factor at the
highest sampling rate. For example, when decimating by a factor of 60 in three stages, decimate by 5,
then by 4, then by 3.
The multirate book references give additional, more specific guidance.

2.3 Implementation
2.3.1 How do I implement decimation?
Decimation consists of the processes of lowpass filtering, followed by downsampling.

To implement the filtering part, you can use either FIR or IIR filters.

To implement the downsampling part (by a downsampling factor of “M”) simply keep every Mth sample,
and throw away the M-1samples in between. For example, to decimate by 4, keep every fourth sample, and
throw three out of every four samples away.

2.3.2 That almost sounds too easy…

Beauty, eh?

2.3.3 If I’m going to throw away most of the lowpass filter’s outputs, why bother to calculate them in the
first place?
You may be onto something. In the case of FIR filters, any output is a function only of the past inputs
(because there is no feedback). Therefore, you only have to calculate outputs which will be used.

For IIR filters, you still have to do part or all of the filter calculation for each input, even when the
corresponding output won’t be used. (Depending on the filter topology used, certain feed-forward parts of
the calculation can be omitted.),. The reason is that outputs you do use are affected by the feedback from the
outputs you don’t use.

The fact that only the outputs which will be used have to be calculated explains why decimating filters are
almost always implemented using FIR filters!

2.4 FIR Decimators

2.4.1 What computational savings do I gain by using a FIR decimator?
Since you compute only one of every M outputs, you save M-1 operations per output, or an overall
“savings” of (M-1)/M. Therefore, the larger the decimation factor is, the larger the savings, percentage-wise.

A simple way to think of the amount of computation required to implement a FIR decimator is that it is
equal to the computation required for a non-decimating N-tap filter operating at the output rate.

2.4.2 How much memory savings do I gain by using a FIR decimator?

DEPARTMENT OF Electronics and Communication Engineering

None. You still have to store every input sample in the FIR’s delay line, so the memory requirement is the
same size as for a non-decimated FIR having the same number of taps.

2.4.3 How do I design a FIR decimator?

Just use your favorite FIR design method. The design criteria are:

1. The passband lower frequency is zero; the passband upper frequency is whatever information
bandwidth you want to preserve after decimating. The passband ripple is whatever your application
can tolerate.
2. The stopband lower frequency is half the output rate minus the passband upper frequency. The
stopband attenuation is set according to whatever aliasing your application can stand. (Note that there
will always be aliasing in a decimator, but you just reduce it to a negligible value with the decimating
filter.)
3. As with any FIR, the number of taps is whatever is required to meet the passband and stopband
specifications.

2.4.4 How do I implement a FIR decimator?

A decimating FIR is actually the same as a regular FIR, except that you shift M samples into the delay line
for each output you calculate. More specifically:

1. Store M samples in the delay line.

2. Calculate the decimated output as the sum-of-products of the delay line values and the filter
coefficients.
3. Shift the delay line by M places to make room for the inputs of the next decimation.
Also, just as with ordinary FIRs, circular buffers can be used to eliminate the requirement to literally shift
the data in the delay line.

2.4.5 Where can I get source code to implement a FIR decimator in C?

Iowegian’s ScopeFIR comes with a free set of multirate algorithms, including FIR decimation functions in
C. Just download and install the ScopeFIR distribution file.

2.4.6 Where can I get assembly code to implement a FIR decimator?

The major DSP vendors provide examples of FIR decimators in their data books and application notes;
check their web sites.

2.4.7 How do I test a FIR decimator?

You can test a decimating FIR in most of the ways you might test an ordinary FIR:

1. A special case of a decimator is an “ordinary” FIR. When given a value of “1” for M, a decimator
should act exactly like an ordinary FIR. You can then do impulse, step, and sine tests on it just like you
can on an ordinary FIR.
2. If you put in a sine whose frequency is within the decimator’s passband, the output should be
distortion-free (once the filter reaches steady-state), and the frequency of the output should be the same
as the frequency of the input, in terms of absolute Hz.
3. You also can extend the “impulse response” test used for ordinary FIRs by using a “fat impulse”,
consisting of M consecutive “1” samples followed by a series of “0” samples. In that case, if the
decimator has been implemented correctly, the output will not be the literal FIR filter coefficients, but
will be the sum of every subset of M coefficients.
4. You can use a step response test. Given a unity-valued step input, the output should be the sum of the
FIR coefficients once the filter has reached steady state.

DEPARTMENT OF Electronics and Communication Engineering

 Interpolation

3.1 Basics
3.1.1 What are “upsampling” and “interpolation”?
“Upsampling” is the process of inserting zero-valued samples between original samples to increase the
sampling rate. (This is called “zero-stuffing”.) Upsampling adds to the original signal undesired spectral
images which are centered on multiples of the original sampling rate.

“Interpolation”, in the DSP sense, is the process of upsampling followed by filtering. (The filtering removes
the undesired spectral images.) As a linear process, the DSP sense of interpolation is somewhat different
from the “math” sense of interpolation, but the result is conceptually similar: to create “in-between” samples
from the original samples. The result is as if you had just originally sampled your signal at the higher rate.

3.1.2 Why interpolate?

The primary reason to interpolate is simply to increase the sampling rate at the output of one system so that
another system operating at a higher sampling rate can input the signal.

3.1.3 What is the “interpolation factor”?

The interpolation factor is simply the ratio of the output rate to the input rate. It is usually symbolized by
“L”, so output rate / input rate=L.

Tip: You can remember that “L” is the symbol for interpolation factor by thinking of “interpo-L-ation”.

3.1.4 Is there a restriction on interpolation factors I can use?

Yes. Since interpolation relies on zero-stuffing you can only interpolate by integer factors; you cannot
interpolate by fractional factors. (However, you can combine interpolation and decimation to achieve an
overall rational factor, for example, 4/5; see Part 4: Resampling.)

3.1.4 Which signals can be interpolated?

All. There is no restriction.

3.1.5 When interpolating, do I always need to do filtering?

Yes. Otherwise, you’re doing upsampling.

3.1.6 OK, you know what I mean…do I always need to do interpolation (upsampling followed by filtering)
or can I get by with doing just upsampling?
Upsampling adds undesired spectral images to the signal at multiples of the original sampling rate, so unless
you remove those by filtering, the upsampled signal is not the same as the original: it’s distorted.

Some applications may be able to tolerate that, for example, if the images get removed later by an analog
filter, but in most applications you will have to remove the undesired images via digital filtering. Therefore,
interpolation is far more common that upsampling alone.

3.2 Multistage
3.2.1 Can I interpolate in multiple stages?
Yes, so long as the interpolation ratio, L, is not a prime number. For example, to interpolate by a factor of
15, you could interpolate by 3 then interpolate by 5. The more factors L has, the more choices you have. For
example you could interpolate by 16 in:

DEPARTMENT OF Electronics and Communication Engineering

  one stage: 16
 two stages: 4 and 4
 three stages: 2, 2, and 4
 four stages: 2, 2, 2, and 2

3.2.2 Cool. But why bother with all that?

Just as with decimation, the computational and memory requirements of interpolation filtering can often be
reduced by using multiple stages.

3.2.3 OK, so how do I figure out the optimum number of stages, and the interpolation ratio at each stage?
There isn’t a simple answer to this one: the answer varies depending on many things. However, here are a
couple of rules of thumb:

 Using two or three stages is usually optimal or near-optimal.

 Interpolate in order of the smallest to largest factors. For example, when interpolating by a factor of 60
in three stages, interpolate by 3, then by 4, then by 5. (Use the largest ratio on the highest rate.)
The multirate book references give additional, more specific guidance.

3.3 Implementation
3.3.1 How do I implement interpolation?
Interpolation always consists of two processes:

1. Inserting L-1 zero-valued samples between each pair of input samples. This operation is called “zero
stuffing”.
2. Lowpass-filtering the result.
The result (assuming an ideal interpolation filter) is a signal at L times the original sampling rate which has
the same spectrum over the input Nyquist (0 to Fs/2) range, and with zero spectral content above the
original Fs/2.

3.3.2 How does that work?

1. The zero-stuffing creates a higher-rate signal whose spectrum is the same as the original over the
original bandwidth, but has images of the original spectrum centered on multiples of the original
sampling rate.
2. The lowpass filtering eliminates the images.

3.3.3 Why do interpolation by zero-stuffing? Doesn’t it make more sense to create the additional samples by
just copying the original samples?
This idea is appealing because, intuitively, this “stairstep” output seems more similar to the original than the
zero-stuffed version. But in this case, intuition leads us down the garden path. This process causes a “zero-
order hold” distortion in the original passband, and still creates undesired images (see below).

Although these effects could be un-done by filtering, it turns out that zero-stuffing approach is not only
more “correct”, it actually reduces the amount of computation required to implement a FIR interpolation
filter. Therefore, interpolation is always done via zero-stuffing.

3.4 FIR Interpolators

3.4.1 How does zero-stuffing reduce computation of the interpolation filter?
The output of a FIR filter is the sum each coefficient multiplied by each corresponding input sample. In the
case of a FIR interpolation filter, some of the input samples are stuffed zeros. Each stuffed zero gets

DEPARTMENT OF Electronics and Communication Engineering

multiplied by a coefficient and summed with the others. However, this adding-and-summing processing has
no effect when the data sample is zero–which we know in advance will be the case for L-1 out of each L
input samples of a FIR interpolation filter. So why bother to calculate these taps?

The net result is that to interpolate by a factor of L, you calculate L outputs for each input using L different
“sub-filters” derived from your original filter.

3.4.2 Can you give me an example of a FIR interpolator?

Here’s an example of a 12-tap FIR filter that implements interpolation by a factor of four. The coefficients
are h0-h11, and three data samples, x0-x2 (with the newest, x2, on the left) have made their way into the
filter’s delay line:

h h h h h h h h h h h1 h1
Result
0 1 2 3 4 5 6 7 8 9 0 1

x2 0 0 0 x1 0 0 0 x0 0 0 0 x2·h0+x1·h4+x0·h8

0 x2 0 0 0 x1 0 0 0 x0 0 0 x2·h1+x1·h5+x0·h9

x2·h2+x1·h6+x0·h1
0 0 x2 0 0 0 x1 0 0 0 x0 0
0

x2·h3+x1·h7+x0·h1
0 0 0 x2 0 0 0 x1 0 0 0 x0
1

3.4.3 What can I generalize from that example?

The table suggests the following general observations about FIR interpolators:

 Since the interpolation ratio is four (L=4), there are four “sub-filters” (whose coefficient sets are
marked here with matching colors.) These sub-filters are officially called “polyphase filters”.
 For each input, we calculate L outputs by doing L basic FIR calculations, each using a different set of
coefficients.
 The number of taps per polyphase filter is 3, or, expressed as a formula: Npoly=Ntotal / L.
 The coefficients of each polyphase filter can be determined by skipping every Lth coefficient, starting
at coefficients 0 through L-1, to calculate corresponding outputs 0 through L-1.
 Alternatively, if you rearranged your coefficients in advance in “scrambled” order like this:
h0, h4, h8, h1, h5, h9, h2, h6, h10, h3, h7, h11
then you could just step through them in order.
 We have hinted here at the fact that N should be a multiple of L. This isn’t absolutely necessary, but if
N isn’t a multiple of L, the added complication of using a non-multiple of L often isn’t worth it. So if
the minimum number of taps that your filter specification requires doesn’t happen to be a multiple of
L, your best bet is usually to just increase N to the next multiple of L. You can do this either by adding
some zero-valued coefficients onto the end of the filter, or by re-designing the filter using the larger N
value.

3.4.4 What computational savings do I gain by using a FIR interpolator?

Since each output is calculated using only N/L coefficients (rather than N coefficients), you get an overall
computational “savings” of (N – N/L) per output .

DEPARTMENT OF Electronics and Communication Engineering

A simple way to think of the amount of computation required to implement a FIR interpolator is that it is
equal to the computation required for a non-interpolating N-tap filter operating at the input rate. In effect,
you have to calculate L filters using N/L taps each, so that’s N total taps calculated per input.

3.4.5 How much memory savings do I gain by using a FIR interpolator?

Compared to the straight-forward implementation of interpolation by upsampling the signal by stuffing it
with L-1 zeros , then filtering it, you save memory by a factor of (L-1)/L. In other words, you don’t have to
store L-1 zero-stuffed “upsamples” per actual input sample.

3.4.6 How do I design a FIR interpolator?

Just use your favorite FIR design method. The design criteria are:

1. TBD

3.5 Implementation
3.5.1 How do I implement a FIR interpolator?
An interpolating FIR is actually the same as a regular FIR, except that, for each input, you calculate L
outputs per input using L polyphase filters, each having N/L taps. More specifically:

1. Store a sample in the delay line. (The size of the delay line is N/L.)
2. For each of L polyphase coefficient sets, calculate an output as the sum-of-products of the delay line
values and the filter coefficients.
3. Shift the delay line by one to make room for the next input.
Also, just as with ordinary FIRs, circular buffers can be used to eliminate the requirement to literally shift
the data in the delay line.

3.5.2 Where can I get source code to implement a FIR interpolator in C?

Iowegian’s ScopeFIR comes with a free set of multirate algorithms, including FIR interpolation functions in
C. Just download and install the ScopeFIR distribution file.

3.5.3 Where can I get assembly code to implement a FIR interpolator?

The major DSP vendors provide examples of FIR interpolators in their data books and application notes, so
check their web sites.

3.5.4 How do I test a FIR interpolator?

You can test an interpolating FIR in most of the ways you might test an ordinary FIR:

1. A special case of an interpolator is an ordinary FIR. When given a value of 1 for L, an interpolator
should act exactly like an ordinary FIR. You can then do impulse, step, and sine tests on it just like you
can on an ordinary FIR.
2. If you put in a sine whose frequency is within the interpolator’s passband, the output should be
distortion-free (once the filter reaches steady state), and the frequency of the output should be the same
as the frequency of the input, in terms of absolute Hz.
3. You can use a step response test. Given a unity-valued step input, every group of L outputs should be
the same as the sums of the coefficients of the L individual polyphase filters, once the filter has
reached steady state.

DEPARTMENT OF Electronics and Communication Engineering

DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
 A

DEPARTMENT OF Electronics and Communication Engineering

20. ADDITIONAL TOPICS IF ANY
Additional/missing topics

 Correlation
 Geortzel algorithm
 FIR Least square design
methods
 Multi stage implementation of sampling rate
conversion

Correlation

The concept of correlation can best be presented with an example. Figure 7-13 shows the key
elements of a radar system. A specially designed antenna transmits a short burst of radio wave
energy in a selected direction. If the propagating wave strikes an object, such as the helicopter in
this illustration, a small fraction of the energy is reflected back toward a radio receiver located
near the transmitter. The transmitted pulse is a specific shape that we have selected, such as the
triangle shown in this example. The received signal will consist of two parts: (1) a shifted and
scaled version of the transmitted pulse, and (2) random noise, resulting from interfering radio
waves, thermal noise in the electronics, etc. Since radio signals travel at a known rate, the speed
of light, the shift between the transmitted and received pulse is a direct measure of the distance to
the object being detected. This is the problem: given a signal of some known shape, what is the
best way to determine where (or if) the signal occurs in another signal. Correlation is the answer.

Correlation is a mathematical operation that is very similar to convolution. Just as with

convolution, correlation uses two signals to produce a third signal. This third signal is called the
cross-correlation of the two input signals. If a signal is correlated with itself, the resulting signal
is instead called the autocorrelation. The convolution machine was presented in the last chapter
to show how convolution is performed. Figure 7-14 is a similar illustration of a correlation
machine. The received signal, x[n], and the cross-correlation signal, y[n], are fixed on the page.
The waveform we are looking for, t[n], commonly called the target signal, is contained within
the correlation machine. Each sample in y[n] is calculated by moving the correlation machine
left or right until it points to the sample being worked on. Next, the indicated samples from
the received signal fall into the correlation machine, and are multiplied by the corresponding
points in the target signal. The sum of these products then moves into the proper sample in the
cross- correlation signal.

DEPARTMENT OF Electronics and Communication Engineering

The amplitude of each sample in the cross-correlation signal is a measure of how much the
received signal resembles the target signal, at that location. This means that a peak will occur in
the cross-correlation signal for every target signal that is present in the received signal. In other
words, the value of the cross-correlation is maximized when the target signal is aligned with the
same features in the received signal.

What if the target signal contains samples with a negative value? Nothing changes. Imagine that
the correlation machine is positioned such that the target signal is perfectly aligned with the
matching waveform in the received signal. As samples from the received signal fall into the
correlation machine, they are multiplied by their matching samples in the target signal.
Neglecting noise, a positive sample will be multiplied by itself, resulting in a positive number.
Likewise, a negative sample will be multiplied by itself, also resulting in a positive number.

DEPARTMENT OF Electronics and Communication Engineering

Even if the target signal is completely negative, the peak in the cross -correlation will still be
positive.

If there is noise on the received signal, there will also be noise on the cross-correlation signal. It
is an unavoidable fact that random noise looks a certain amount like any target signal you can
choose. The noise on the cross-correlation signal is simply measuring this similarity. Except for
this noise, the peak generated in the cross-correlation signal is symmetrical between its left
and right. This is true even if the target signal isn't symmetrical. In addition, the width of the
peak is twice the width of the target signal. Remember, the cross-correlation is trying to detect
the target signal, not recreate it. There is no reason to expect that the peak will even look like
the target signal.

Correlation is the optimal technique for detecting a known waveform in random noise. That is,
the peak is higher above the noise using correlation than can be produced by any other linear
system. (To be perfectly correct, it is only optimal for random white noise). Using correlation to
detect a known waveform is frequently called matched filtering.

The correlation machine and convolution machine are identical, except for one small difference.
As discussed in the last chapter, the signal inside of the convolution machine is flipped left-for-
right. This means that samples numbers: 1, 2, 3 … run from the right to the left. In the
correlation machine this flip doesn't take place, and the samples run in the normal direction.

DEPARTMENT OF Electronics and Communication Engineering

DEPARTMENT OF Electronics and Communication Engineering
Geortzel algorithm

The Goertzel algorithm is a digital signal processing (DSP) technique for identifying
frequency components of a signal, published by Gerald Goertzel in 1958. While the general Fast
Fourier transform (FFT) algorithm computes evenly across the bandwidth of the incoming
signal, the Goertzel algorithm looks at specific, predetermined frequencies.

A practical application of this algorithm is recognition of the DTMF tones produced by

the buttons pushed on a telephone keypad

It can also be used "in reverse" as a sinusoid synthesis function, which requires only
1 multiplication and 1 subtraction per sample.

Explanation of algorithm

The Goertzel algorithm computes a sequence, s(n), given an input sequence, x(n):

s(n) = x(n) + 2cos(2πω)s(n − 1) − s(n − 2)

where s( − 2) = s( − 1) = 0 and ω is some frequency of interest, in cycles per sample, which

should be less than 1/2. This effectively implements a second-order IIR filter with poles at
e
+ 2πiω − 2πiω
and e , and requires only one multiplication (assuming 2cos(2πω) is pre-
computed), one addition and one subtraction per input sample. For real inputs, these
operations are real.

The Z transform of this process is

Applying an additional, FIR, transform of the form

will give an overall transform of

The time-domain equivalent of this overall transform is

DEPARTMENT OF Electronics and Communication Engineering

 ,

DEPARTMENT OF Electronics and Communication Engineering

 which becomes, assuming x(k) = 0 for all k < 0

or, the equation for the (n + 1)-sample DFT of x, evaluated for ω and
+ 2πiωn
multiplied by the scale factor e .

Note that applying the additional transform Y(z)/S(z) only requires the
last two samples of the s sequence. Consequently, upon processing N
samples x(0)...x(N − 1), the last two samples from the s sequence can be
used to compute the value of a DFT bin, which corresponds to the
chosen frequency ω as
− 2πiω(N − 1) − 2πiω − 2πiω(N − 1)
X(ω) = y(N − 1)e = (s(N − 1) − e s(N − 2))e

For the special case often found when computing DFT bins, where
ωN = k for some integer, k, this simplifies to
− 2πiω + 2πiω + 2πiω
X(ω) = (s(N − 1) − e s(N − 2))e =e s(N − 1) − s(N − 2)

In either case, the corresponding power can be computed using

the same cosine term required to compute s as

2 2
X(ω)X'(ω) = s(N − 2) + s(N − 1) − 2cos(2πω)s(N − 2)s(N − 1

Least-Squares Linear-Phase FIR Filter

Design

Let the FIR filter length be L+1 samples, with even, and suppose we'll initially design it
to be centered about the time origin. Then the frequency response is given on our frequency

grid by

Enforcing even symmetry in the impulse response, i.e., , gives a zero phase

DEPARTMENT OF Electronics and Communication Engineering

FIR filter which we can later right-shift samples to make a causal, linear phase
filter. In this case, the frequency response reduces to a sum of cosines:

or in matrix
form:

(Note that Remez exchange algorithms are also based on this formulation
internally.)

Matrix Formulation: Optimal Design,

Cont'd

In matrix notation, our filter design problem can be

stated

where

DEPARTMENT OF Electronics and Communication Engineering

and is the desired frequency response at the specified frequencies.

Least Squares Optimization

Hence we can minimize

Expanding this, we have:

This is quadratic in , hence it has a global minimum which we can find by taking
the derivative, setting it to zero, and solving for . Doing this yields:

These are the famous normal equations whose solution is given by:

The matrix

is known as the (Moore-Penrose) pseudo-inverse of the matrix .

Geometrical Interpretation of Least Squares

Typically, the number of frequency constraints is much greater than the number
of design variables (filter taps). In these cases, we have an overdetermined
system of equations (more equations than unknowns). Therefore, we cannot

DEPARTMENT OF Electronics and Communication Engineering

generally satisfy all the equations, and we are left with minimizing some error
criterion to find the ``optimal compromise'' solution.

In the case of least-squares approximation, we are minimizing the Euclidean

distance, which suggests the geometrical interpretation shown in Fig.4.28.

Thus, the desired vector is the vector sum of its best least-squares approximation
plus an orthogonal error :

In practice, the least-squares solution can be found by minimizing the sum of squared
errors:

Figure 4.28 suggests that the error vector is orthogonal to the column space of
the matrix , hence it must be orthogonal to each column in :

This is how the orthogonality principle can be used to derive the fact that the best least
squares solution is given by

DEPARTMENT OF Electronics and Communication Engineering

Note that the pseudo-inverse projects the vector onto the column space of .

(Note: To obtain the best numerical algorithms for least-squares solution in Matlab, it is

usually better to use ``x = A b'' rather than explicitly computing the pseudo-inverse as in
``x = pinv(A) * b''.)

Sampling Rate Conversion by Stages

The decimator and interpolator discussed so far are of a single-stage structure. When

large changes in sampling rate are required, multiple stages of sample rate conversion

are found

to be more computationally efficient. Most practical systems employ a multi-stage

structure, resulting in a considerable relaxation in the specifications of anti-aliasing

(decimation)

or anti-imaging (interpolation) filters in each stage compared to a single stage

realization. The decimation in Figure 3.23 can be realized in two stages if the

decimation factor

D can be expressed as a product of two integers, D1 and D2. Referring to Figure

3.24, in the first stage, the signal x(n) is decimated by a factor of D1. The output, v(p)

is further decimated by D2 in the second stage resulting in an overall decimation of

x(n) by

DEPARTMENT OF Electronics and Communication Engineering

 Figure 3.23: Decimation in a Single Stage.

Figure 3.24: Decimation in Two Stages.

D = (D1D2). The filters H1(z) and H2(z) are so designed that the aliasing in the band of

interest is belowa prescribed level and that the overall passband and stopband tolerances

are met. This multi-stage sampling rate conversion system offers less computation and

more flexibility in filter design. An example is given below to illustrate the idea of multi -

stage sampling rate conversion.

Example: Multi-Stage Sampling Rate Conversion

We have a discrete time signal with a sampling rate of 90 kHz. The signal has the

desired information in the frequency band from 0 to 450 Hz (passband), and the band from

450 to 500 Hz is the transition band. The signal is to be decimated by a factor of

ninety. The required tolerances are a passband ripple of 0.002 and a stopband ripple of

0.001. Decimation in a Single Stage

First we consider a single-stage design as shown in Figure 3.25(a). The specifications

of the required LPF are shown in Figure 3.25(b).

According to the formula by Kaiser, the approximate length of an FIR filter is

given by

DEPARTMENT OF Electronics and Communication Engineering

where peak passband ripple (linear) δp = 0.002, peak stopband ripple (linear) δs = 0.001,

normalized transition bandwidth passband edge frequency fp = 450 Hz,

stopband edge frequency fs = 500 Hz, and sampling frequency Fs = 90 kHz.

From Equation 3.37, the lowpass FIR filter H(z) has a length of N ≈ 5424. Therefore,

the number of multiplications per second, Msec, needed for this single-stage decimator is

Since only one out of ninety samples is actually used, the computation rate is based on

the decimated signal rate.

Decimation in Two Stages

Let us now consider the two-stage implementation of the decimation process as shown in

Figure 3.26.

DEPARTMENT OF Electronics and Communication Engineering

Figure 3.25: (a) Block Diagram for Single-Stage Decimation, (b) The Filter Specification.

Figure 3.26: Block Diagram for Multi-Stage Decimation.

Due to the cascade decomposition, each of the two filters, H1(z) and H2(z), must have

a linear passband ripple specification half of that specified for the single-stage filter,

H(z). The stopband ripple specifications for these two filters can be the same as that o f

H(z) since the cascade connection will only reduce the stopband ripple.

Stage One

The first stage will decimate the input signal x(n) by a factor of forty-five. The filter

DEPARTMENT OF Electronics and Communication Engineering

specifications for the first-stage LPF H1(z) are

Figure 3.27: Decimation Filter Design for Stage One.

These specifications are shown in Figure 3.27.

The reason for choosing this value of the stopband edge is that, after decim ation by a

factor of forty-five, the residual energy of the signal in the band from 1000 to 2000 Hz

will be aliased back to the band from 0 to 1000 Hz. Due to the attenuation in the stopband,

the energy of the signal in the band from 1800 to 2000 Hz is very small compared to that

in

1000 to 1800 Hz. So the amount of aliasing in the desired band of interest (0 to 450 Hz)

will also be small, resulting in very little signal distortion.

According to Equation 3.37, the approximate length of the FIR filter, H1(z) is N1 =

276. The number of multiplications per second for the first stage is

DEPARTMENT OF Electronics and Communication Engineering

Stage Two

The specifications for the second-stage filter, H2(z), are

Figure 3.28 shows the characteristics of H2(z). This stage will perform a decimation

of factor two on the output signal of the first stage. So, the total decimation of x(n) is by

a factor of ninety as required.

Figure 3.28: Decimation Filter Design for Stage Two.

For the second stage, the length of the filter, as calculated from Equation 3.37, is N2 =

129. The number of multiplications required for this stage is

The total number of multiplications per second required for the two-stage implementation

DEPARTMENT OF Electronics and Communication Engineering

of the decimator is

So, the two-stage implementation requires only of the operation required of the

single-stage implementation

Decimation in Three Stages

To further illustrate the concept of multi-stage implementation of decimator and

interpolator, we will now consider the three-stage implementation as shown in Figure 3.29.

Stage One

In this stage, decimation by fifteen is performed on the input signal x(n). The

characteristics of the LPF, H1(z), are shown in Figure 3.30. The filter specifications

are

DEPARTMENT OF Electronics and Communication Engineering

DEPARTMENT OF Electronics and Communication Engineering
 Figure 3.30: Decimation Filter Design for Stage One.

As in the two-stage case, the choice of stopband edge frequency can be extended to

the point for which negligible aliasing occurs in the passband (band of interest).

The approximate length of the filter as given by Equation 3.37 is N1 = 60. The

number of multiplications per second for this stage is calculated as

Stage Two

In this stage, a decimation by a factor of three is done. The specifications of the LPF in

this stage, H2(z), are

As before, the stopband edge frequency can be stretched out to 1800 Hz. The filter

DEPARTMENT OF Electronics and Communication Engineering

characteristics are shown in Figure 3.31.

The length of filter required for this stage is N2 = 20 and the number of multiplications per

second is

Stage Three

The third stage performs a decimation of factor two on the output of the second stage. The

specifications of the LPF, H3(z), in this stage are

Figure 3.31: Decimation Filter Design for Stage Two.

Figure 3.32 shows the specifications for H3(z). As before, the

approximate length of

the filter, as calculated from Equation 3.37, is N3 = 134. The

number of multiplications required per second in the third stage is

DEPARTMENT OF Electronics and Communication Engineering

 Figure 3.32: Decimation Filter Design for Stage Three.

The total number of multiplications in the three stages of implementation is

Compared to the single-stage implementation, the number of multiplications per second

are reduced by a factor of by using three stages.

From this example, we can see that a significant saving in computation as well as in storage can

be achieved by a multi-stage decimator and interpolator design. These savings depend on the

optimum design of the number of stages and the choice of decimation factor for the individual

stages.

The examples illustrate the many different combinations and ordering possible. One approach is

to determine the sets of I and D factors that satisfy the filtering requirements and then estimate

the storage and computational costs for each set. The lowest cost solution is then selected.

Walsh transform:
1D walsh transform :

DEPARTMENT OF Electronics and Communication Engineering

Inverse 1D Walsh transform:

2D Walsh transform:

2D Inverse Walsh transform:

Hadmard Transform:

Discrete Cosine transform:

DCT of 2D image of size (NXM) is defined by the following equation

where

DEPARTMENT OF Electronics and Communication Engineering

18. Additional Topics

19. Known Gaps, if any.

Known Gaps:

1. Comparison of received signal with the reference signal (not just by cross correlation but
by using time shift parameter).

2. Filters to alter time domain characteristics.

3. To differentiate between the signals based on geographical position of the signal source.

4. To study the systems those are non-linear and time- variant.

20. Discussion Topics, if any.

. Find the IDFT of Y (k) = (1, 0, 1, 0)

2. Compute the Fourier transform of the signal x(n) = u(n) – u(n-1).
3. Determine the Discrete Fourier transform x (n) = (1, 1, 1, 1)
4. Derive and draw the 8 point FFT-DIT butterfly structure.
5. Discuss the properties of DFT.
6. Obtain the cascade and parallel form realizations for the following systems (16)
Y (n) = -0.1(n-1) + 0.2 y (n-2) + 3x (n) +3.6 x (n-1) +0.6 x (n-2)
7. Design a HPF of length 7 with cut off frequency of 2 rad/sec using Hamming window. Plot
the magnitude and phase response.
8.Convert the analog filter H(s) = 0.5 (s+4) / (s+1)(s+2) using impulse invariant
transformation T=0.31416s.

8. Obtain the cascade form realizations of FIR systems

H (z) = 1+5/2 z-1+ 2z-2 +2 z-3

9. Determine the system function H(z) of the lowest order Chebyshev digital filter that meet the
following specifications:

a. 3-dB ripple in the passband 0≤│w│≤0.24 .

DEPARTMENT OF Electronics and Communication Engineering

 b. At least 50-dB attenuation in the stop band 0.35 ≤│w│≤ by using Bilinear transformatio
technique.

10. Design a Butterworth analog low-pass filter is required to meet the following specifications:

Pass band ripple: < 1 dB

Pass band edge frequency: 4 kHz
Stop band attenuation: > 40 dB
Stop band edge frequency: 6kHz
Sample rate: 24 kHz

DEPARTMENT OF Electronics and Communication Engineering

21. University/Autonomous Question papers of previous years.

DEPARTMENT OF Electronics and Communication Engineering

DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
DEPARTMENT OF Electronics and Communication Engineering
 III B.TECH - II SEMESTER EXAMINATIONS, APRIL/MAY, 2011
DIGITAL SIGNAL PROCESSING
(COMMON TO EEE, ECE, EIE, ETM, ICE)
Time: 3hours Max. Marks: 80
Answer any FIVE questions
All Questions Carry Equal Marks

1.a) Define an LTI System and show that the output of an LTI system is
given by the convolution of Input sequence and impulse response.
b) Prove that the system defined by the following difference equation
is an LTI
system y(n) = x(n+1)-3x(n)+x(n-1) ; n≥0.

[8+8]

2.a) Define DFT and IDFT. State any Four properties of DFT.
b) Find 8-Point DFT of the given time domain sequence x(n) = {1, 2, 3,
4}. [8+8]

3.a) Derive the expressions for computing the FFT using DIT algorithm
and hence draw the standard butterfly structure.
b) Compare the computational complexity of FFT and DFT. [8+8]

4. Discuss and draw various IIR realization structures like Direct form
– I, Direct form-II, Parallel and cascade forms for the difference
equation given y(n) = - 3/8 Y(n-1) + 3/32 y(n-2) + 1/64 y(n-3) + x(n) +
3 x(n-1) + 2 x(n-2).
5.a) Compare Butterworth and Chebyshev approximation techniques.
b) Design a Digital Butterworth LPF using Bilinear transformation
technique for the following specifications
0.707 ≤ | H(w) | ≤ 1 ; 0
≤ w ≤ 0.2π
| H(w) | ≤ 0.08 ; 0.4 π ≤ w ≤ [ 8+8]

6.a) Compare FIR and IIR filters

b) Design an FIR Digital High pass filter using Hamming window whose
cut off freq is 1.2 rad/s and length of window N=9. [8+8]

7.a) Define Multirate systems and Sampling rate conversion

b) Discuss the process of n Decimation by a factor D and explain how
the aliasing effect can be eliminated. [8+8]

8. Discuss various Modified Bus structures of Programmable DSP

Processors.[16]

DEPARTMENT OF Electronics and Communication Engineering

 III B.TECH - II SEMESTER EXAMINATIONS, APRIL/MAY, 2011
DIGITAL SIGNAL PROCESSING
(COMMON TO EEE, ECE, EIE, ETM, ICE)
Time: 3hours Max. Marks: 80
Answer any FIVE questions
All Questions Carry Equal Marks

1.a) Write short notes on classification of systems.

b) Derive BIBO stability criteria to achieve stability of a system. [8+8]

2.a) Define DFS. State any Four properties of DFS.

b) Find the IDFT of the given sequence x(K) = {2, 2-3j, 2+3j, -2}. [8+8]

3.a) Find X(K) of the given sequence x(n) = { 1,2,3,4,4,3,2,1}using

DIT- FFT
algorithm.
b) Compare the computational complexity of FFT and DFT. [8+8]

4. What are the various basic building blocks in realization of Digital

Systems and hence discuss transposed form realization structures.

5.a) Compare Impulse Invariant and Bilinear transformation techniques.

b) Compute the poles of an Analog Chebyshev filter TF that satisfies the
Constraints
0.707 ≤ | H(jΩ)| ≤ 1 ; 0 ≤ Ω≤ 2
| H(jΩ)| ≤ 0.1 ; Ω ≥ 4
and determine Ha(s) and hence obtain H(z) using Bilinear transformation. [16]

6.a) Derive the conditions to achieve Linear Phase characteristics of FIR filters
b) Design an FIR Digital Low pass filter using Hanning window whose
cut off freq is 2 rad/s and length of window N=9. [8+8]

7.a) Discuss the implementation of Polyphase filters for Interpolators with

an example b) Discuss the sampling rate conversion by a factor I/D
with the help of a Neat block Diagram. [8+8]

8. Write short notes on:

a) VLIW Architecture of Programmable Digital Signal Processors
b) Multiplier and Multiplier Accumulator [8+8]

DEPARTMENT OF Electronics and Communication Engineering

DEPARTMENT OF Electronics and Communication Engineering
40
41
42
 III B.TECH - II SEMESTER EXAMINATIONS,
APRIL/MAY, 2011
DIGITAL SIGNAL PROCESSING
(COMMON TO EEE, ECE, EIE, ETM, ICE)
Time: 3hours Max. Marks: 80
Answer any FIVE questions
All Questions Carry Equal Marks

1.a) Discuss various discrete time sequences.

b) Give the Basic block diagram of Digital Signal Processor. [8+8]

2.a) Define DFS. State any Four properties of DFS.

b) Find the IDFT of the given sequence x(K) = {2, 2-3j, 2+3j, -2}. [8+8]

3.a) Find IFFT of the given X(K) = { 1,2,3,4,4,3,2,1}using DIF algorithm

b) Bring out the relationship between DFT and Z-transform. [8+8]

4.a) Define Z-Transform and List out its properties.

b) Discuss Direct form, Cascade and Linear phase realization
structures of FIR
filters.

[8+8]

5.a) Discuss digital and analog frequency transformation techniques.

b) Discuss IIR filter design using Bilinear transformation and hence discuss
frequency warping effect. [8+8]

6.a) Compare various windowing functions.

b) Design an FIR Digital Low pass filter using rectangular window
whose cut off freq is 2 rad/s and length of window N=9. [8+8]

7.a) Define Interpolation and Decimation. List out the advantages of

Sampling rate conversion.
b) Discuss the sampling rate conversion by a factor I with the help of a
Neat block
Diagram.[8+8]

8.a) Discuss Various Addressing modes of Programmable

Digital Signal Processors. b) Give the Internal Architecture of
TMS320C5X 16 bit fixed point processor.[ 8+8]

43
III B.TECH - II SEMESTER EXAMINATIONS, APRIL/MAY, 2011
DIGITAL SIGNAL PROCESSING
(COMMON TO EEE, ECE, EIE, ETM, ICE)
Time: 3hours Max. Marks: 80
Answer any FIVE questions
All Questions Carry Equal Marks

1.a) Define Linearity, Time Invariant, Stability and Causality.

b) The discrete time system is represented by the following difference
2
equations in which x(n) is input and y(n) is output. Y(n) = 3y (n-1)-
nx(n)+4x(n-1)-2x(n-1).
[8+8]

2.a) Define Convolution. Compare Linear and Circular

Convolution techniques. b) Find the Linear convolution of the
given two sequences x(n)={1,2} and
h(n) ={1,2,3} using DFT and IDFT. [8+8]

3.a) Develop DIT-FFT algorithm and draw signal flow graphs for
decomposing the
DFT for N=6 by considering the factors for N = 6 = 2.3.
b) Bring out the relationship between DFT and Z-transform. [8+8]

4.a) Discuss transposed form structures with an example.

b) Discuss Direct form, Cascade realization structures of FIR filters. [8+8]

5.a) Discuss digital and analog frequency transformation techniques.

b) Discuss IIR filter design using Impulse Invariant transformation
and list out its advantages and Limitations. [8+8]

6.a) Compare various windowing functions

b) Design an FIR Digital Band pass filter using rectangular window
whose upper and lower cut off freq.’s are 1 & 2 rad/s and length of
window N = 9. [8+8]

7.a) Define Interpolation and Decimation.

b) Discuss the sampling rate conversion by a factor I/D with the help of
a Neat block
Diagram. [8+8]

8.a) Write a short notes on On-Chip peripherals of Programmable DSP’s.

b)Give the Internal Architecture of TMS320C5X 16 bit fixed point
processor. [8+8]

44
 III B.TECH - II SEMESTER EXAMINATIONS, APRIL/MAY, 2011
DIGITAL SIGNAL PROCESSING
(COMMON TO EEE, ECE, EIE, ETM, ICE)
Time: 3hours Max. Marks: 80
Answer any FIVE questions
All Questions Carry Equal Marks
1. (a) Discuss impulse invariance method of deriving IIR digital filter
from corre- sponding analog filter.
(b) Use the Bilinear transformation to convert the analog filter with
system func- tion H (S) = S + 0.1/(S + 0.1)2 + 9 into a digital
IIR filters. Select T = 0.1 and compare the location of the zeros
in H(Z) with the locations of the zeros obtained by applying the
impulse invariance method in the conversion of H(S). [8+8]

2. (a) Design a high pass filter using hamming window with a cut-off frequency of
1.2 radians/second and N=9
(b) Compare FIR and IIR filters. [10+6]

3. (a) For each of the following systems, determine whether or not the
system is i. stable
ii. causal
iii. linear
iv. shift-invariant.

A. T [x(n)] = x(n − n0 ) B. T
[x(n)] = ex (n)
C. T[x(n)] = a x(n) + b.
Justify your answer.
(b) A system is described by the difference equation y(n)-y(n-1)-y(n-2) = x(n-
1). Assuming that the system is initially relaxed, determine its unit sample
response h(n). [8+8]

4. (a) Implement the decimation in time FFT algorithm for N=16.

n - trivial multiplications are

5. (a) Discuss the frequency-domain representation of discrete-time

systems and sig- nals by deriving the necessary relation.
(b) Draw the frequency response of LSI system with impulse response
h(n) = an u(−n) (|a| < 1)

6. (a) Describe how targets can be decided using RADAR

(b) Give an expression for the following parameters relative to

45
 RADAR
i. Beam width
ii. Maximum unambiguous range
(c) Discuss signal processing in a RADAR system. [6+6+4]

7. (a) An LTI system is described by the equation y(n)=x(n)+0.81x(n-

1)-0.81x(n-
2)-0.45y(n-2). Determine the transfer function of the system.
Sketch the poles and zeroes on the Z-plane.
(b) Define stable and unstable systems. Test the condition for
stability of the first-order IIR filter governed by the equation
y(n)=x(n)+bx(n-1). [8+8]

8. (a) Compute Discrete Fourier transform of the following finite

length sequence considered to be of length N.
i. x(n) = δ(n + n0 ) where 0 < n0 < N
ii. x(n) = an where 0 < a < 1.
(b) If x(n) denotes a finite length sequence of length N, show that
x((−n))N =x((N − n))N . [8+8]

46
III B.TECH - II SEMESTER EXAMINATIONS, APRIL/MAY, 2011
DIGITAL SIGNAL PROCESSING
(COMMON TO EEE, ECE, EIE, ETM, ICE)
Time: 3hours Max. Marks: 80
Answer any FIVE questions
All Questions Carry Equal Marks

1. (a) An LTI system is described by the equation y(n)=x(n)+0.81x(n-1)-0.81x(n-

2)-0.45y(n-2). Determine the transfer function of the system.
Sketch the poles and zeroes on the Z-plane.
(b) Define stable and unstable systems. Test the condition for
stability of the first-order IIR filter governed by the equation
y(n)=x(n)+bx(n-1). [8+8]

2. (a) Discuss impulse invariance method of deriving IIR digital filter

from corre- sponding analog filter.
(b) Use the Bilinear transformation to convert the analog filter with
system func- tion H (S) = S + 0.1/(S + 0.1)2 + 9 into a digital
IIR filters. Select T = 0.1 and compare the location of the zeros
in H(Z) with the locations of the zeros obtained by applying the
impulse invariance method in the conversion of H(S). [8+8]
3. (a) Describe how targets can be decided using RADAR
(b) Give an expression for the following parameters relative to
RADAR
i. Beam width
ii. Maximum unambiguous range
(c) Discuss signal processing in a RADAR system. [6+6+4]

4. (a) Discuss the frequency-domain representation of discrete-time

systems and sig- nals by deriving the necessary relation.
(b) Draw the frequency response of LSI system with impulse response
h(n) = an u(−n) (|a| < 1) [8+8]

5. (a) For each of the following systems, determine whether or not

the system is i. stable
ii. causal
iii. linear
iv. shift-invariant.

A. T [x(n)]

47
 = x(n −
n0 ) B. T
[x(n)] =
ex (n)
C. T[x(n)] =
a x(n) + b.
Justify
your
answer.
(b) A system is described by the difference equation y(n)-y(n-1)-y(n-
2) = x(n-
1). Assuming that the system is initially relaxed, determine its
unit sample response h(n). [8+8]

6. (a) Implement the decimation in time FFT algorithm for N=16.

n - trivial multiplications are

7. (a) Design a high pass filter using hamming window with a cut-off
frequency of
1.2 radians/second and N=9
(b) Compare FIR and IIR filters. [10+6]

8. (a) Compute Discrete Fourier transform of the following finite

length sequence considered to be of length N.
i. x(n) = δ(n + n0 ) where 0 < n0 < N
ii. x(n) = an where 0 < a < 1.
(b) If x(n) denotes a finite length sequence of length N, show that
x((−n))N =x((N − n))N .[8+8]

48
 III B.TECH - II SEMESTER EXAMINATIONS, APRIL/MAY, 2011
DIGITAL SIGNAL PROCESSING
(COMMON TO EEE, ECE, EIE, ETM, ICE)
Time: 3hours Max. Marks: 80
Answer any FIVE questions
All Questions Carry Equal Marks

1. (a) Design a high pass filter using hamming window with a cut-off
frequency of
1.2 radians/second and N=9
(b) Compare FIR and IIR filters. [10+6]

2. (a) Describe how targets can be decided using RADAR

(b) Give an expression for the following parameters relative to RADAR
i. Beam width
ii. Maximum unambiguous range
(c) Discuss signal processing in a RADAR system. [6+6+4]

3. (a) An LTI system is described by the equation y(n)=x(n)+0.81x(n-

1)-0.81x(n-
2)-0.45y(n-2). Determine the transfer function of the system. Sketch the poles and zeroes
on the Z-plane.
(b) Define stable and unstable systems. Test the condition for
stability of the first-order IIR filter governed by the equation y(n)=x(n)+bx(n-1). [8+8]

4. (a) Compute Discrete Fourier transform of the following finite length sequence
considered to be of length N.
i. x(n) = δ(n + n0 ) where 0 < n0 < N
ii. x(n) = an where 0 < a < 1.
(b) If x(n) denotes a finite length sequence of length N, show that x((−n))N =
x((N − n))N . [8+8]

5. (a) For each of the following systems, determine whether or not the system is i. stable
ii. causal
iii. linear
iv. shift-invariant.

A. T [x(n)] = x(n − n0 ) B. T [x(n)] =

49
 ex (n)
C. T[x(n)] = a x(n) + b.
Justify your answer
(b) A system is described by the difference equation y(n)-y(n-1)-y(n-2) = x(n-
1). Assuming that the system is initially relaxed, determine its unit sample response h(n).
[8+8]

6. (a) Discuss the frequency-domain representation of discrete-time systems and sig- nals by
deriving the necessary relation.
(b) Draw the frequency response of LSI system with impulse response
h(n) = an u(−n) (|a| < 1) [8+8]

7. (a) Implement the decimation in time FFT algorithm for N=16.

(b) In the above Question how many non - trivial multiplications are R e q u i r e d .
8. (a) Discuss impulse invariance method of deriving IIR digital filter from corre- sponding
analog filter.
(b) Use the Bilinear transformation to convert the analog filter with system
func- tion H (S) = S + 0.1/(S + 0.1)2 + 9 into a digital IIR filters. Select
T = 0.1 and compare the location of the zeros in H(Z) with the locations of
the zeros obtained by applying the impulse invariance method in the
conversion of H(S). [8+8

50
 III B.TECH - II SEMESTER EXAMINATIONS, APRIL/MAY, 2011
DIGITAL SIGNAL PROCESSING
(COMMON TO EEE, ECE, EIE, ETM, ICE)
Time: 3hours Max. Marks: 80
Answer any FIVE questions
All Questions Carry Equal Marks

1. (a) Compute Discrete Fourier transform of the following finite

length sequence considered to be of length N.
i. x(n) = δ(n + n0 ) where 0 < n0 < N
ii. x(n) = an where 0 < a < 1.
(b) If x(n) denotes a finite length sequence of length N, show that x((−n))N =
x((N − n))N . [8+8]
2. (a) For each of the following systems, determine whether or
not the system is i. stable
ii. causal
iii. linear
iv. shift-invariant.

A. T [x(n)] = x(n − n0 ) B. T [x(n)] = ex (n)

C. T[x(n)] = a x(n) + b.
Justify your answer.
(b) A system is described by the difference equation y(n)-y(n-1)-y(n-2) = x(n-
1). Assuming that the system is initially relaxed, determine its unit sample
response h(n). [8+8]

[6+6+4]
4. (a) Design a high pass filter using hamming window with a cut-off frequency of
1.2 radians/second and N=9
(b) Compare FIR and IIR filters. [10+6]
5. (a) An LTI system is described by the equation y(n)=x(n)+0.81x(n-1)-0.81x(n-
2)-0.45y(n-2). Determine the transfer function of the system. Sketch the
poles and zeroes on the Z-plane

51
(b) Define stable and unstable systems. Test the condition for
stability of the first-order IIR filter governed by the equation y(n)=x(n)+bx(n-1).
[8+8]

6. (a) Discuss the frequency-domain representation of discrete-time systems

and sig- nals by deriving the necessary relation.
(b) Draw the frequency response of LSI system with impulse response
h(n) = an u(−n) (|a| < 1) [8+8]

7. (a) Discuss impulse invariance method of deriving IIR digital filter

from corre- sponding analog filter.
(b) Use the Bilinear transformation to convert the analog filter with system function H
(S) = S + 0.1/(S + 0.1)2 + 9 into a digital IIR filters. Select T = 0.1 and compare the
location of the zeros in H(Z)with the locations of the zeros obtained by applying the
impulse invariance method in the conversion of H(S). [8+8]

8. (a) Implement the decimation in time FFT algorithm forN=16.

(b) In the above Question how many non - trivial
Multiplications are required.

2)Design an FIR digital high pass filter using hamming window whose cutoff frequency is 1.2
rad/s and length of window N=5.Compare the same using rectangular window .Draw the
frequency response curve for both the cases. (5M)

3a) Determine the order and poles of a type- I low pass Chebyshev filter that has a 1-dB ripple
in the passband, a cutoff frequency Ωp=1000π, a stopband frequency of 2000π, and an
attenuation of 40 dB or more for Ω ≥ Ωs .
(3 M)

3b) Compare frequency sampling and windowing method of filter design (2 M)

( OR)

4a) Explain how aliasing effect can be avoided while performing decimation process of a
signal by a factor of D. (3 M)

4bWhat is meant by overflow error and how it can be avoided? (2 M)

52
53
54
55
56
22. References, Journals, Websites and E-Links, if required.
References

1. Discrete time signal Processing-A.V.Oppenheim and R.W.Schaffer,PHI,2009.

2. Digital signal Processing-Fundamentals and Applications-LiTan,Elsevier,2008
3. Fundamentals of Digital signal Processing using MAT Lab-Robert J.Schilling,Sandra
L.Harris,Thomson,2007.
4. Digital signal Processing-S.Salivahanan,A.Vallavaraj,C.Gnanapriya,TMH,2009.

Text Books

1. Digital signal Processing: Principles,Algorithams and Applications-John

G.Proakis,D.G.Manolakis,4th Edition,Pearson/PHI,2009.
2. Digital signal Processing-A Practical Approach-Emmanuel C.Ifeacher,Barrie.W.Jervis,2 nd
Edition,Pearson Education,2009.

WEBSITES

1.www.google.com
2.www.dspguru.com
3.www.nptel.ac.in
4.www.nptelonlinecourses.iitm.ac.in

23. Quality Control Sheets.

a) Course end survey
b) Feedback on Teaching Learning Process.
c) CO – attainment.

Hard copy will be attach at the end of the semester.

24. Student List

To be attached

25. Group-wise students list for discussion topics

To be attached

57
58