ELEG 305: Digital Signal Processing

Size: px
Start display at page:

Download "ELEG 305: Digital Signal Processing"

Transcription

1 ELEG 305: Digital Signal Processing Lecture 19: Lattice Filters Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 2008 K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19 Outline 1 Review of Previous Lecture 2 Lecture Objectives 3 Lattice to Direct Form Translation K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19

2 Review of Previous Lecture Review of Previous Lecture Efficient FFT computation for real sequences For x 1 (n) and x 2 (n) real, set x(n) x 1 (n) + jx 2 (n) then X 1 (k) 1 2 [X(k) + X (N k)]; X 2 (k) 1 2j [X(k) X (N k)] Similar tricks hold for length 2N real sequences Linear filtering computation of the DFT the Goertzel algorithm y k (n) W k N y k(n 1) + x(n) gives X(k) y k (n) nn Implementations of FIR discrete time systems Direct, Cascade, and Lattice filter structures K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19 Lecture Objectives Lecture Objectives Objective Derive lattice structures for FIR filters; Develop a method for converting lattice filter coefficients to direct form (FIR) coefficients Reading Chapters 9 (9.2); Next lecture, complete lattice filters, structures for IIR systems (Chapter ); start filter design (Chapter ) K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19

3 Objective: Developed a lattice structure realization of FIR filters Approach: Suppose we have a sequence of FIR filters H m (z) A m (z) m 0, 1,..., M 1 where by definition A m (z) 1 + α m (k)z k m 1 k1 Suppose m 1. Then the output of H 1 (z), for input x(n), is y(n) x(n) + α 1 (1)x(n 1) ( ) To build a cascade of 1 st order stages, generalize the notation. Let f 0 (n) x(n) g 0 (n) x(n) f 1 (n) f 0 (n) + K 1 g 0 (n 1) g 1 (n) K 1 f 0 (n) + g 0 (n 1) Question: For what K 1 value is ( ) realized, i.e., f 1 (n) y(n)? K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19 Single Stage Lattice Filter To realize ( ), let K 1 α 1 (1). Then f 1 (n) f 0 (n) + K 1 g 0 (n 1) x(n) + α 1 (1)x(n 1) y(n) [order m 1 filter output] Similarly, g 1 (n) α 1 (1)x(n) + x(n 1) Note: The stage 1 governing equations are f 1 (n) f 0 (n) + K 1 g 0 (n 1) g 1 (n) K 1 f 0 (n) + g 0 (n 1) where we set the reflection coefficient as K 1 α 1 (1) to realize ( ) K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19

4 Objective: Add a second stage and equate to the m 2 filter Note that [m 2 filter] H 2 (z) A 2 (z) 1 + α 2 (1)z 1 + α 2 (2)z 2 y(n) x(n) + α 2 (1)x(n 1) + α 2 (2)x(n 2) ( ) The two stage lattice filter is given by The governing equations are f 1 (n) f 0 (n) + K 1 g 0 (n 1) f 2 (n) f 1 (n) + K 2 g 1 (n 1) g 1 (n) K 1 f 0 (n) + g 0 (n 1) g 2 (n) K 2 f 1 (n) + g 1 (n 1) Procedure: Set K 1 and K 2 such that ( ) is realized Note: f 0 (n) g 0 (n) x(n) [pipeline input] and f 2 (n) y(n) [output] K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19 [lattice f 1 (n) f 0 (n) + K 1 g 0 (n 1) f 2 (n) f 1 (n) + K 2 g 1 (n 1) equations] g 1 (n) K 1 f 0 (n) + g 0 (n 1) g 2 (n) K 2 f 1 (n) + g 1 (n 1) By substitution f 2 (n) f 1 (n) + K 2 g 1 (n 1) [f 0 (n) + K 1 g 0 (n 1)] + K 2 [K 1 f 0 (n 1) + g 0 (n 2)] [x(n) + K 1 x(n 1)] + K 2 [K 1 x(n 1) + x(n 2)] x(n) + K 1 (1 + K 2 )x(n 1) + K 2 x(n 2) ( ) Recalling ( ) [direct FIR form] y(n) x(n) + α 2 (1)x(n 1) + α 2 (2)x(n 2) and equating with ( ) [lattice form] yields or, equivalently, α 2 (2) K 2 and α 2 (1) K 1 (1 K 2 ) K 2 α 2 (2) and K 1 α 2(1) 1 + α 2 (2) K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19 ( )

5 Process: This is repeated for M 1 stages, with general recursion f 0 (n) g 0 (n) x(n) f m (n) f m 1 (n) + K m g m 1 (n 1), m 1, 2,..., M 1 g m (n) K m f m 1 (n) + g m 1 (n 1), m 1, 2,..., M 1 At the final stage y(n) f m 1 (n) m 1 α m 1 (k)x(n k) K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19 Observation: At each stage there are two outputs, f m (n) and g m (n) Consider g 2 (n). First & second stage equations: f 1 (n) f 0 (n) + K 1 g 0 (n 1) f 2 (n) f 1 (n) + K 2 g 1 (n 1) g 1 (n) K 1 f 0 (n) + g 0 (n 1) g 2 (n) K 2 f 1 (n) + g 1 (n 1) Apply substitution similarly to the previous case g 2 (n) K 2 f 1 (n) + g 1 (n 1) K 2 [f 0 (n) + K 1 g 0 (n 1)] + [K 1 f 0 (n 1) + g 0 (n 2)] K 2 x(n) + K 1 (1 + K 2 )x(n 1) + x(n 2) using the prior result K 2 α 2 (2) and K 1 (1 K 2 ) α 2 (1), g 2 (n) α 2 (2)x(n) + α 2 (1)x(n 1) + x(n 2) Result: The coefficients for g m (n) are reverse order those for f m (n) K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19

6 Result: The two outputs at arbitrary stage m are expressed as f m (n) α m (k)x(n k) and g m (n) β m (k)x(n k) where β m (k) α m (m k), k 0, 1,..., m Note: f m (n) is the forward prediction; g m (n) the backward prediction Taking the z transform Similarly F m (z) A m (z)x(z) or A m (z) F m(z) X(z) F m(z) F 0 (z) G m (z) B m (z)x(z) or B m (z) G m(z) X(z) G m(z) G 0 (z) K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19 Recall A m (z) α m (k)z k Question: How does B m (z) relate to A m (z)? B m (z) β m (k)z k [substitute β m (k) α m (m k)] α m (m k)z k [let l m k] α m (l)z l m l0 m z m α m (l)z l l0 z m A m (z 1 ) Result: B(z) has reciprocal zeros of A(z); B(z) is the reciprocal, or reverse, polynomial of A(z) K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19

7 Lattice Filter Representation Summary Case 1: FIR filter representations f m (n) α m (k)x(n k) and g m (n) β m (k)x(n k) Case 2: z domain representations F m (z) A m (z)x(z) or A m (z) F m(z) X(z) G m (z) B m (z)x(z) or B m (z) G m(z) X(z) Also, B m (z) z m A m (z 1 ), m 1, 2,..., M 1 Case 3: Recursion lattice representations f 0 (n) g 0 (n) x(n) f m (n) f m 1 (n) + K m g m 1 (n 1), m 1, 2,..., M 1 g m (n) K m f m 1 (n) + g m 1 (n 1), m 1, 2,..., M 1 K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19 Taking z transform of the recursion lattice representations F 0 (z) G 0 (z) X(z) F m (z) F m 1 (z) + K m z 1 G m 1 (z), m 1, 2,..., M 1 G m (z) K m F m 1 (z) + z 1 G m 1 (z), m 1, 2,..., M 1 Dividing by X(z) and using A m (z) Fm(z) X(z) and B m(z) Gm(z) X(z) Case 4: z domain recursion lattice representations A 0 (z) B 0 (z) 1 A m (z) A m 1 (z) + K m z 1 B m 1 (z), m 1, 2,..., M 1 B m (z) K m A m 1 (z) + z 1 B m 1 (z), m 1, 2,..., M 1 or in matrix notation [ ] [ ] [ ] Am (z) 1 Km A m 1 (z) B m (z) K m 1 z 1 B m 1 (z) K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19

8 Lattice to Direct Form Translation Design Questions: How do we convert lattice coefficients (K terms) to direct-form filter coefficients (α terms)? How do we convert direct form coefficients to lattice coefficients? Case 1: Lattice to Direct Form Translation Approach: Utilize relations A 0 (z) B 0 (z) 1 A m (z) A m 1 (z) + K m z 1 B m 1 (z), m 1, 2,..., M 1 B m (z) z m A m (z 1 ), m 1, 2,..., M 1 Solve for α terms recursively, starting with m 1 Example Determine the direct form realization of a three stage (M 4) lattice filter with coefficients K 1 1/2, K 2 1/2, K 3 1/4. K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19 Lattice to Direct Form Translation Start recursion set m 1 A 1 (z) A 0 (z) + K 1 z 1 B 0 (z) z 1 ( ) Set ( ) equal to α 1 polynomial A 1 (z) 1 + α 1 (1)z 1 α 1 (1) 1 2 Note, B m (z) is the reverse polynomial of A m (z) B 1 (z) z 1 A 1 (z 1 ) z 1 ( z ) z 1 K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19

9 Lattice to Direct Form Translation For m 2, use recursion A m (z) A m 1 (z) + K m z 1 B m 1 (z) A 2 (z) A 1 (z) + K 2 z 1 B 1 (z) [ ] [ ( )] 1 1 z z z 1 equating α terms yields z z 2 α 2 (1) 3 4 and α 2 (2) 1 2 Reversing the polynomial of A m (z) to get B m (z) B 2 (z) z 2 A 2 (z 1 ) z 1 + z 2 K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19 For m 3, Thus A 3 (z) A 2 (z) + K 3 z 1 B 2 (z) [ z ] z z z z 3 Lattice to Direct Form Translation [ ( z )] 4 z 1 + z 2 α 3 (1) 7 8, α 3(2) 11 16, α 3(3) 1 4 Result: The system in direct-form: (where α 3 (0) 1) y(n) 3 α 3 (k)x(n k) x(n) x(n 1) x(n 2) + 1 x(n 3) 4 K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19

10 Lecture Summary Lecture Summary Lattice Filters general recursion f 0 (n) g 0 (n) x(n) f m (n) f m 1 (n) + K m g m 1 (n 1) g m (n) K m f m 1 (n) + g m 1 (n 1), Lattice to Direct Form Translation utilize recursion on A 0 (z) B 0 (z) 1 A m (z) A m 1 (z) + K m z 1 B m 1 (z), m 1, 2,..., M 1 B m (z) z m A m (z 1 ), m 1, 2,..., M 1 Next lecture Complete lattice filters, structures for IIR systems (Chapter ); start filter design (Chapter ) K. E. Barner (Univ. of Delaware) ELEG 305: Digital Signal Processing Fall / 19

ELEG 305: Digital Signal Processing

ELEG 305: Digital Signal Processing ELEG 305: Digital Signal Processing Lecture 18: Applications of FFT Algorithms & Linear Filtering DFT Computation; Implementation of Discrete Time Systems Kenneth E. Barner Department of Electrical and

More information

ELEG 305: Digital Signal Processing

ELEG 305: Digital Signal Processing ELEG 5: Digital Signal Processing Lecture 6: The Fast Fourier Transform; Radix Decimatation in Time Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 8 K.

More information

Lecture 19 IIR Filters

Lecture 19 IIR Filters Lecture 19 IIR Filters Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/5/10 1 General IIR Difference Equation IIR system: infinite-impulse response system The most general class

More information

ELEG 305: Digital Signal Processing

ELEG 305: Digital Signal Processing ELEG 305: Digital Signal Processing Lecture : Design of Digital IIR Filters (Part I) Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 008 K. E. Barner (Univ.

More information

ELEG 305: Digital Signal Processing

ELEG 305: Digital Signal Processing ELEG 305: Digital Signal Processing Lecture 4: Inverse z Transforms & z Domain Analysis Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 008 K. E. Barner

More information

ELEG 305: Digital Signal Processing

ELEG 305: Digital Signal Processing ELEG 305: Digital Signal Processing Lecture 1: Course Overview; Discrete-Time Signals & Systems Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 2008 K. E.

More information

R13 SET - 1

R13 SET - 1 R13 SET - 1 III B. Tech II Semester Regular Examinations, April - 2016 DIGITAL SIGNAL PROCESSING (Electronics and Communication Engineering) Time: 3 hours Maximum Marks: 70 Note: 1. Question Paper consists

More information

Discrete-time signals and systems

Discrete-time signals and systems Discrete-time signals and systems 1 DISCRETE-TIME DYNAMICAL SYSTEMS x(t) G y(t) Linear system: Output y(n) is a linear function of the inputs sequence: y(n) = k= h(k)x(n k) h(k): impulse response of the

More information

DHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A

DHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A DHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A Classification of systems : Continuous and Discrete

More information

APPLIED SIGNAL PROCESSING

APPLIED SIGNAL PROCESSING APPLIED SIGNAL PROCESSING DIGITAL FILTERS Digital filters are discrete-time linear systems { x[n] } G { y[n] } Impulse response: y[n] = h[0]x[n] + h[1]x[n 1] + 2 DIGITAL FILTER TYPES FIR (Finite Impulse

More information

Introduction to Digital Signal Processing

Introduction to Digital Signal Processing Introduction to Digital Signal Processing 1.1 What is DSP? DSP is a technique of performing the mathematical operations on the signals in digital domain. As real time signals are analog in nature we need

More information

Lecture 7 Discrete Systems

Lecture 7 Discrete Systems Lecture 7 Discrete Systems EE 52: Instrumentation and Measurements Lecture Notes Update on November, 29 Aly El-Osery, Electrical Engineering Dept., New Mexico Tech 7. Contents The z-transform 2 Linear

More information

Digital Signal Processing Lecture 4

Digital Signal Processing Lecture 4 Remote Sensing Laboratory Dept. of Information Engineering and Computer Science University of Trento Via Sommarive, 14, I-38123 Povo, Trento, Italy Digital Signal Processing Lecture 4 Begüm Demir E-mail:

More information

Signals and Systems. Problem Set: The z-transform and DT Fourier Transform

Signals and Systems. Problem Set: The z-transform and DT Fourier Transform Signals and Systems Problem Set: The z-transform and DT Fourier Transform Updated: October 9, 7 Problem Set Problem - Transfer functions in MATLAB A discrete-time, causal LTI system is described by the

More information

EEL3135: Homework #4

EEL3135: Homework #4 EEL335: Homework #4 Problem : For each of the systems below, determine whether or not the system is () linear, () time-invariant, and (3) causal: (a) (b) (c) xn [ ] cos( 04πn) (d) xn [ ] xn [ ] xn [ 5]

More information

ECE4270 Fundamentals of DSP Lecture 20. Fixed-Point Arithmetic in FIR and IIR Filters (part I) Overview of Lecture. Overflow. FIR Digital Filter

ECE4270 Fundamentals of DSP Lecture 20. Fixed-Point Arithmetic in FIR and IIR Filters (part I) Overview of Lecture. Overflow. FIR Digital Filter ECE4270 Fundamentals of DSP Lecture 20 Fixed-Point Arithmetic in FIR and IIR Filters (part I) School of ECE Center for Signal and Information Processing Georgia Institute of Technology Overview of Lecture

More information

Analog LTI system Digital LTI system

Analog LTI system Digital LTI system Sampling Decimation Seismometer Amplifier AAA filter DAA filter Analog LTI system Digital LTI system Filtering (Digital Systems) input output filter xn [ ] X ~ [ k] Convolution of Sequences hn [ ] yn [

More information

Linear Convolution Using FFT

Linear Convolution Using FFT Linear Convolution Using FFT Another useful property is that we can perform circular convolution and see how many points remain the same as those of linear convolution. When P < L and an L-point circular

More information

Z - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.

Z - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals. Z - Transform The z-transform is a very important tool in describing and analyzing digital systems. It offers the techniques for digital filter design and frequency analysis of digital signals. Definition

More information

z-transforms Definition of the z-transform Chapter

z-transforms Definition of the z-transform Chapter z-transforms Chapter 7 In the study of discrete-time signal and systems, we have thus far considered the time-domain and the frequency domain. The z- domain gives us a third representation. All three domains

More information

SIDDHARTH GROUP OF INSTITUTIONS:: PUTTUR Siddharth Nagar, Narayanavanam Road QUESTION BANK (DESCRIPTIVE)

SIDDHARTH GROUP OF INSTITUTIONS:: PUTTUR Siddharth Nagar, Narayanavanam Road QUESTION BANK (DESCRIPTIVE) SIDDHARTH GROUP OF INSTITUTIONS:: PUTTUR Siddharth Nagar, Narayanavanam Road 517583 QUESTION BANK (DESCRIPTIVE) Subject with Code : Digital Signal Processing(16EC422) Year & Sem: III-B.Tech & II-Sem Course

More information

LECTURE NOTES DIGITAL SIGNAL PROCESSING III B.TECH II SEMESTER (JNTUK R 13)

LECTURE NOTES DIGITAL SIGNAL PROCESSING III B.TECH II SEMESTER (JNTUK R 13) LECTURE NOTES ON DIGITAL SIGNAL PROCESSING III B.TECH II SEMESTER (JNTUK R 13) FACULTY : B.V.S.RENUKA DEVI (Asst.Prof) / Dr. K. SRINIVASA RAO (Assoc. Prof) DEPARTMENT OF ELECTRONICS AND COMMUNICATIONS

More information

Review of Discrete-Time System

Review of Discrete-Time System Review of Discrete-Time System Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides were based on class notes developed by Profs. K.J. Ray Liu and Min Wu.

More information

EE 521: Instrumentation and Measurements

EE 521: Instrumentation and Measurements Aly El-Osery Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA November 1, 2009 1 / 27 1 The z-transform 2 Linear Time-Invariant System 3 Filter Design IIR Filters FIR Filters

More information

DSP Configurations. responded with: thus the system function for this filter would be

DSP Configurations. responded with: thus the system function for this filter would be DSP Configurations In this lecture we discuss the different physical (or software) configurations that can be used to actually realize or implement DSP functions. Recall that the general form of a DSP

More information

Let H(z) = P(z)/Q(z) be the system function of a rational form. Let us represent both P(z) and Q(z) as polynomials of z (not z -1 )

Let H(z) = P(z)/Q(z) be the system function of a rational form. Let us represent both P(z) and Q(z) as polynomials of z (not z -1 ) Review: Poles and Zeros of Fractional Form Let H() = P()/Q() be the system function of a rational form. Let us represent both P() and Q() as polynomials of (not - ) Then Poles: the roots of Q()=0 Zeros:

More information

EE 123: Digital Signal Processing Spring Lecture 23 April 10

EE 123: Digital Signal Processing Spring Lecture 23 April 10 EE 123: Digital Signal Processing Spring 2007 Lecture 23 April 10 Lecturer: Anant Sahai Scribe: Chuo J. Liu 23.1 Outline In this lecture, we talked about implementation issue for rational transfer function.

More information

/ (2π) X(e jω ) dω. 4. An 8 point sequence is given by x(n) = {2,2,2,2,1,1,1,1}. Compute 8 point DFT of x(n) by

/ (2π) X(e jω ) dω. 4. An 8 point sequence is given by x(n) = {2,2,2,2,1,1,1,1}. Compute 8 point DFT of x(n) by Code No: RR320402 Set No. 1 III B.Tech II Semester Regular Examinations, Apr/May 2006 DIGITAL SIGNAL PROCESSING ( Common to Electronics & Communication Engineering, Electronics & Instrumentation Engineering,

More information

VU Signal and Image Processing

VU Signal and Image Processing 052600 VU Signal and Image Processing Torsten Möller + Hrvoje Bogunović + Raphael Sahann torsten.moeller@univie.ac.at hrvoje.bogunovic@meduniwien.ac.at raphael.sahann@univie.ac.at vda.cs.univie.ac.at/teaching/sip/18s/

More information

DEPARTMENT OF EI DIGITAL SIGNAL PROCESSING ASSIGNMENT 1

DEPARTMENT OF EI DIGITAL SIGNAL PROCESSING ASSIGNMENT 1 This PDF is Created by Simpo PDF Creator unregistered version - http://wwwsimpopdfcom Study hard, for the well is deep, and our brains are shallow DEPARTMENT OF EI DIGITAL SIGNAL PROCESSING ASSIGNMENT

More information

ECE503: Digital Signal Processing Lecture 6

ECE503: Digital Signal Processing Lecture 6 ECE503: Digital Signal Processing Lecture 6 D. Richard Brown III WPI 20-February-2012 WPI D. Richard Brown III 20-February-2012 1 / 28 Lecture 6 Topics 1. Filter structures overview 2. FIR filter structures

More information

Lecture 8 - IIR Filters (II)

Lecture 8 - IIR Filters (II) Lecture 8 - IIR Filters (II) James Barnes (James.Barnes@colostate.edu) Spring 2009 Colorado State University Dept of Electrical and Computer Engineering ECE423 1 / 27 Lecture 8 Outline Introduction Digital

More information

2. Typical Discrete-Time Systems All-Pass Systems (5.5) 2.2. Minimum-Phase Systems (5.6) 2.3. Generalized Linear-Phase Systems (5.

2. Typical Discrete-Time Systems All-Pass Systems (5.5) 2.2. Minimum-Phase Systems (5.6) 2.3. Generalized Linear-Phase Systems (5. . Typical Discrete-Time Systems.1. All-Pass Systems (5.5).. Minimum-Phase Systems (5.6).3. Generalized Linear-Phase Systems (5.7) .1. All-Pass Systems An all-pass system is defined as a system which has

More information

VALLIAMMAI ENGINEERING COLLEGE. SRM Nagar, Kattankulathur DEPARTMENT OF INFORMATION TECHNOLOGY. Academic Year

VALLIAMMAI ENGINEERING COLLEGE. SRM Nagar, Kattankulathur DEPARTMENT OF INFORMATION TECHNOLOGY. Academic Year VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur- 603 203 DEPARTMENT OF INFORMATION TECHNOLOGY Academic Year 2016-2017 QUESTION BANK-ODD SEMESTER NAME OF THE SUBJECT SUBJECT CODE SEMESTER YEAR

More information

UNIT V FINITE WORD LENGTH EFFECTS IN DIGITAL FILTERS PART A 1. Define 1 s complement form? In 1,s complement form the positive number is represented as in the sign magnitude form. To obtain the negative

More information

DISCRETE-TIME SIGNAL PROCESSING

DISCRETE-TIME SIGNAL PROCESSING THIRD EDITION DISCRETE-TIME SIGNAL PROCESSING ALAN V. OPPENHEIM MASSACHUSETTS INSTITUTE OF TECHNOLOGY RONALD W. SCHÄFER HEWLETT-PACKARD LABORATORIES Upper Saddle River Boston Columbus San Francisco New

More information

Discrete Time Systems

Discrete Time Systems Discrete Time Systems Valentina Hubeika, Jan Černocký DCGM FIT BUT Brno, {ihubeika,cernocky}@fit.vutbr.cz 1 LTI systems In this course, we work only with linear and time-invariant systems. We talked about

More information

Digital Signal Processing Lecture 5

Digital Signal Processing Lecture 5 Remote Sensing Laboratory Dept. of Information Engineering and Computer Science University of Trento Via Sommarive, 14, I-38123 Povo, Trento, Italy Digital Signal Processing Lecture 5 Begüm Demir E-mail:

More information

E : Lecture 1 Introduction

E : Lecture 1 Introduction E85.2607: Lecture 1 Introduction 1 Administrivia 2 DSP review 3 Fun with Matlab E85.2607: Lecture 1 Introduction 2010-01-21 1 / 24 Course overview Advanced Digital Signal Theory Design, analysis, and implementation

More information

ECSE 512 Digital Signal Processing I Fall 2010 FINAL EXAMINATION

ECSE 512 Digital Signal Processing I Fall 2010 FINAL EXAMINATION FINAL EXAMINATION 9:00 am 12:00 pm, December 20, 2010 Duration: 180 minutes Examiner: Prof. M. Vu Assoc. Examiner: Prof. B. Champagne There are 6 questions for a total of 120 points. This is a closed book

More information

Ch. 7: Z-transform Reading

Ch. 7: Z-transform Reading c J. Fessler, June 9, 3, 6:3 (student version) 7. Ch. 7: Z-transform Definition Properties linearity / superposition time shift convolution: y[n] =h[n] x[n] Y (z) =H(z) X(z) Inverse z-transform by coefficient

More information

Lecture 11 FIR Filters

Lecture 11 FIR Filters Lecture 11 FIR Filters Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/4/12 1 The Unit Impulse Sequence Any sequence can be represented in this way. The equation is true if k ranges

More information

Tunable IIR Digital Filters

Tunable IIR Digital Filters Tunable IIR Digital Filters We have described earlier two st-order and two nd-order IIR digital transfer functions with tunable frequency response characteristics We shall show now that these transfer

More information

Multi-rate Signal Processing 3. The Polyphase Representation

Multi-rate Signal Processing 3. The Polyphase Representation Multi-rate Signal Processing 3. The Polyphase Representation Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides were based on class notes developed by

More information

An Fir-Filter Example: Hanning Filter

An Fir-Filter Example: Hanning Filter An Fir-Filter Example: Hanning Filter Josef Goette Bern University of Applied Sciences, Biel Institute of Human Centered Engineering - microlab Josef.Goette@bfh.ch February 7, 2018 Contents 1 Mathematical

More information

Advanced Digital Signal Processing -Introduction

Advanced Digital Signal Processing -Introduction Advanced Digital Signal Processing -Introduction LECTURE-2 1 AP9211- ADVANCED DIGITAL SIGNAL PROCESSING UNIT I DISCRETE RANDOM SIGNAL PROCESSING Discrete Random Processes- Ensemble Averages, Stationary

More information

Digital Signal Processing Lecture 10 - Discrete Fourier Transform

Digital Signal Processing Lecture 10 - Discrete Fourier Transform Digital Signal Processing - Discrete Fourier Transform Electrical Engineering and Computer Science University of Tennessee, Knoxville November 12, 2015 Overview 1 2 3 4 Review - 1 Introduction Discrete-time

More information

Use: Analysis of systems, simple convolution, shorthand for e jw, stability. Motivation easier to write. Or X(z) = Z {x(n)}

Use: Analysis of systems, simple convolution, shorthand for e jw, stability. Motivation easier to write. Or X(z) = Z {x(n)} 1 VI. Z Transform Ch 24 Use: Analysis of systems, simple convolution, shorthand for e jw, stability. A. Definition: X(z) = x(n) z z - transforms Motivation easier to write Or Note if X(z) = Z {x(n)} z

More information

Recursive, Infinite Impulse Response (IIR) Digital Filters:

Recursive, Infinite Impulse Response (IIR) Digital Filters: Recursive, Infinite Impulse Response (IIR) Digital Filters: Filters defined by Laplace Domain transfer functions (analog devices) can be easily converted to Z domain transfer functions (digital, sampled

More information

EE 225D LECTURE ON DIGITAL FILTERS. University of California Berkeley

EE 225D LECTURE ON DIGITAL FILTERS. University of California Berkeley University of California Berkeley College of Engineering Department of Electrical Engineering and Computer Sciences Professors : N.Morgan / B.Gold EE225D Digital Filters Spring,1999 Lecture 7 N.MORGAN

More information

INF3440/INF4440. Design of digital filters

INF3440/INF4440. Design of digital filters Last week lecture Today s lecture: Chapter 8.1-8.3, 8.4.2, 8.5.3 INF3440/INF4440. Design of digital filters October 2004 Last week lecture Today s lecture: Chapter 8.1-8.3, 8.4.2, 8.5.3 Last lectures:

More information

Z-Transform. x (n) Sampler

Z-Transform. x (n) Sampler Chapter Two A- Discrete Time Signals: The discrete time signal x(n) is obtained by taking samples of the analog signal xa (t) every Ts seconds as shown in Figure below. Analog signal Discrete time signal

More information

Digital Signal Processing: Signal Transforms

Digital Signal Processing: Signal Transforms Digital Signal Processing: Signal Transforms Aishy Amer, Mohammed Ghazal January 19, 1 Instructions: 1. This tutorial introduces frequency analysis in Matlab using the Fourier and z transforms.. More Matlab

More information

Question Bank. UNIT 1 Part-A

Question Bank. UNIT 1 Part-A FATIMA MICHAEL COLLEGE OF ENGINEERING & TECHNOLOGY Senkottai Village, Madurai Sivagangai Main Road, Madurai -625 020 An ISO 9001:2008 Certified Institution Question Bank DEPARTMENT OF ELECTRONICS AND COMMUNICATION

More information

Sistemas de Aquisição de Dados. Mestrado Integrado em Eng. Física Tecnológica 2015/16 Aula 6-26 de Outubro

Sistemas de Aquisição de Dados. Mestrado Integrado em Eng. Física Tecnológica 2015/16 Aula 6-26 de Outubro Sistemas de Aquisição de Dados Mestrado Integrado em Eng. Física Tecnológica 2015/16 Aula 6-26 de Outubro Flash Decoder Thermometer code Wired NOR based decoder 2 Successive Approximation ADC (SAR) CONVERT

More information

Very useful for designing and analyzing signal processing systems

Very useful for designing and analyzing signal processing systems z-transform z-transform The z-transform generalizes the Discrete-Time Fourier Transform (DTFT) for analyzing infinite-length signals and systems Very useful for designing and analyzing signal processing

More information

1. Z-transform: Initial value theorem for causal signal. = u(0) + u(1)z 1 + u(2)z 2 +

1. Z-transform: Initial value theorem for causal signal. = u(0) + u(1)z 1 + u(2)z 2 + 1. Z-transform: Initial value theorem for causal signal u(0) lim U(z) if the limit exists z U(z) u(k)z k u(k)z k k lim U(z) u(0) z k0 u(0) + u(1)z 1 + u(2)z 2 + CL 692 Digital Control, IIT Bombay 1 c Kannan

More information

INFINITE-IMPULSE RESPONSE DIGITAL FILTERS Classical analog filters and their conversion to digital filters 4. THE BUTTERWORTH ANALOG FILTER

INFINITE-IMPULSE RESPONSE DIGITAL FILTERS Classical analog filters and their conversion to digital filters 4. THE BUTTERWORTH ANALOG FILTER INFINITE-IMPULSE RESPONSE DIGITAL FILTERS Classical analog filters and their conversion to digital filters. INTRODUCTION 2. IIR FILTER DESIGN 3. ANALOG FILTERS 4. THE BUTTERWORTH ANALOG FILTER 5. THE CHEBYSHEV-I

More information

DSP Algorithm Original PowerPoint slides prepared by S. K. Mitra

DSP Algorithm Original PowerPoint slides prepared by S. K. Mitra Chapter 11 DSP Algorithm Implementations 清大電機系林嘉文 cwlin@ee.nthu.edu.tw Original PowerPoint slides prepared by S. K. Mitra 03-5731152 11-1 Matrix Representation of Digital Consider Filter Structures This

More information

EE482: Digital Signal Processing Applications

EE482: Digital Signal Processing Applications Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 02 DSP Fundamentals 14/01/21 http://www.ee.unlv.edu/~b1morris/ee482/

More information

ESE 531: Digital Signal Processing

ESE 531: Digital Signal Processing ESE 531: Digital Signal Processing Lec 22: April 10, 2018 Adaptive Filters Penn ESE 531 Spring 2018 Khanna Lecture Outline! Circular convolution as linear convolution with aliasing! Adaptive Filters Penn

More information

Lecture 2 Discrete-Time LTI Systems: Introduction

Lecture 2 Discrete-Time LTI Systems: Introduction Lecture 2 Discrete-Time LTI Systems: Introduction Outline 2.1 Classification of Systems.............................. 1 2.1.1 Memoryless................................. 1 2.1.2 Causal....................................

More information

(Refer Slide Time: )

(Refer Slide Time: ) Digital Signal Processing Prof. S. C. Dutta Roy Department of Electrical Engineering Indian Institute of Technology, Delhi FIR Lattice Synthesis Lecture - 32 This is the 32nd lecture and our topic for

More information

Lecture 18: Stability

Lecture 18: Stability Lecture 18: Stability ECE 401: Signal and Image Analysis University of Illinois 4/18/2017 1 Stability 2 Impulse Response 3 Z Transform Outline 1 Stability 2 Impulse Response 3 Z Transform BIBO Stability

More information

Examples. 2-input, 1-output discrete-time systems: 1-input, 1-output discrete-time systems:

Examples. 2-input, 1-output discrete-time systems: 1-input, 1-output discrete-time systems: Discrete-Time s - I Time-Domain Representation CHAPTER 4 These lecture slides are based on "Digital Signal Processing: A Computer-Based Approach, 4th ed." textbook by S.K. Mitra and its instructor materials.

More information

Z Transform (Part - II)

Z Transform (Part - II) Z Transform (Part - II). The Z Transform of the following real exponential sequence x(nt) = a n, nt 0 = 0, nt < 0, a > 0 (a) ; z > (c) for all z z (b) ; z (d) ; z < a > a az az Soln. The given sequence

More information

III. Time Domain Analysis of systems

III. Time Domain Analysis of systems 1 III. Time Domain Analysis of systems Here, we adapt properties of continuous time systems to discrete time systems Section 2.2-2.5, pp 17-39 System Notation y(n) = T[ x(n) ] A. Types of Systems Memoryless

More information

Digital Filter Structures

Digital Filter Structures Chapter 8 Digital Filter Structures 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 8-1 Block Diagram Representation The convolution sum description of an LTI discrete-time system can, in principle, be used to

More information

2.161 Signal Processing: Continuous and Discrete Fall 2008

2.161 Signal Processing: Continuous and Discrete Fall 2008 MI OpenCourseWare http://ocw.mit.edu.6 Signal Processing: Continuous and Discrete Fall 008 For information about citing these materials or our erms of Use, visit: http://ocw.mit.edu/terms. Massachusetts

More information

UNIVERSITY OF OSLO. Please make sure that your copy of the problem set is complete before you attempt to answer anything.

UNIVERSITY OF OSLO. Please make sure that your copy of the problem set is complete before you attempt to answer anything. UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Examination in INF3470/4470 Digital signal processing Day of examination: December 9th, 011 Examination hours: 14.30 18.30 This problem set

More information

Lecture 3 Matlab Simulink Minimum Phase, Maximum Phase and Linear Phase Systems

Lecture 3 Matlab Simulink Minimum Phase, Maximum Phase and Linear Phase Systems Lecture 3 Matlab Simulink Minimum Phase, Maximum Phase and Linear Phase Systems Lester Liu October 31, 2012 Minimum Phase, Maximum Phase and Linear Phase LTI Systems In this section, we will explore the

More information

Lecture 8 - IIR Filters (II)

Lecture 8 - IIR Filters (II) Lecture 8 - IIR Filters (II) James Barnes (James.Barnes@colostate.edu) Spring 24 Colorado State University Dept of Electrical and Computer Engineering ECE423 1 / 29 Lecture 8 Outline Introduction Digital

More information

CMPT 889: Lecture 5 Filters

CMPT 889: Lecture 5 Filters CMPT 889: Lecture 5 Filters Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University October 7, 2009 1 Digital Filters Any medium through which a signal passes may be regarded

More information

Digital Filters. Linearity and Time Invariance. Linear Time-Invariant (LTI) Filters: CMPT 889: Lecture 5 Filters

Digital Filters. Linearity and Time Invariance. Linear Time-Invariant (LTI) Filters: CMPT 889: Lecture 5 Filters Digital Filters CMPT 889: Lecture 5 Filters Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University October 7, 29 Any medium through which a signal passes may be regarded as

More information

ESE 531: Digital Signal Processing

ESE 531: Digital Signal Processing ESE 531: Digital Signal Processing Lec 6: January 30, 2018 Inverse z-transform Lecture Outline! z-transform " Tie up loose ends " Regions of convergence properties! Inverse z-transform " Inspection " Partial

More information

Lecture 7 - IIR Filters

Lecture 7 - IIR Filters Lecture 7 - IIR Filters James Barnes (James.Barnes@colostate.edu) Spring 204 Colorado State University Dept of Electrical and Computer Engineering ECE423 / 2 Outline. IIR Filter Representations Difference

More information

Homework #1 Solution

Homework #1 Solution February 7, 4 Department of Electrical and Computer Engineering University of Wisconsin Madison ECE 734 VLSI Array Structures for Digital Signal Processing Homework # Solution Due: February 6, 4 in class.

More information

Analysis of Finite Wordlength Effects

Analysis of Finite Wordlength Effects Analysis of Finite Wordlength Effects Ideally, the system parameters along with the signal variables have infinite precision taing any value between and In practice, they can tae only discrete values within

More information

Responses of Digital Filters Chapter Intended Learning Outcomes:

Responses of Digital Filters Chapter Intended Learning Outcomes: Responses of Digital Filters Chapter Intended Learning Outcomes: (i) Understanding the relationships between impulse response, frequency response, difference equation and transfer function in characterizing

More information

We have shown earlier that the 1st-order lowpass transfer function

We have shown earlier that the 1st-order lowpass transfer function Tunable IIR Digital Filters We have described earlier two st-order and two nd-order IIR digital transfer functions with tunable frequency response characteristics We shall show now that these transfer

More information

Detailed Solutions to Exercises

Detailed Solutions to Exercises Detailed Solutions to Exercises Digital Signal Processing Mikael Swartling Nedelko Grbic rev. 205 Department of Electrical and Information Technology Lund University Detailed solution to problem E3.4 A

More information

DFT & Fast Fourier Transform PART-A. 7. Calculate the number of multiplications needed in the calculation of DFT and FFT with 64 point sequence.

DFT & Fast Fourier Transform PART-A. 7. Calculate the number of multiplications needed in the calculation of DFT and FFT with 64 point sequence. SHRI ANGALAMMAN COLLEGE OF ENGINEERING & TECHNOLOGY (An ISO 9001:2008 Certified Institution) SIRUGANOOR,TRICHY-621105. DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING UNIT I DFT & Fast Fourier

More information

! Circular Convolution. " Linear convolution with circular convolution. ! Discrete Fourier Transform. " Linear convolution through circular

! Circular Convolution.  Linear convolution with circular convolution. ! Discrete Fourier Transform.  Linear convolution through circular Previously ESE 531: Digital Signal Processing Lec 22: April 18, 2017 Fast Fourier Transform (con t)! Circular Convolution " Linear convolution with circular convolution! Discrete Fourier Transform " Linear

More information

How to manipulate Frequencies in Discrete-time Domain? Two Main Approaches

How to manipulate Frequencies in Discrete-time Domain? Two Main Approaches How to manipulate Frequencies in Discrete-time Domain? Two Main Approaches Difference Equations (an LTI system) x[n]: input, y[n]: output That is, building a system that maes use of the current and previous

More information

IT6502 DIGITAL SIGNAL PROCESSING Unit I - SIGNALS AND SYSTEMS Basic elements of DSP concepts of frequency in Analo and Diital Sinals samplin theorem Discrete time sinals, systems Analysis of discrete time

More information

Digital Signal Processing

Digital Signal Processing DIGITAL SIGNAL PROCESSING SUBJECT CODE : NO. OF LECTURE HRS/WEEK : 04 TOTAL NO. OF LECTURE HRS. : 52 IA MARKS : 25 EXAM HOURS : 03 EXAMMARKS : 100 UNIT - 1 DISCRETE FOURIER TRANSFORMS (DFT): FREQUENCY

More information

Your solutions for time-domain waveforms should all be expressed as real-valued functions.

Your solutions for time-domain waveforms should all be expressed as real-valued functions. ECE-486 Test 2, Feb 23, 2017 2 Hours; Closed book; Allowed calculator models: (a) Casio fx-115 models (b) HP33s and HP 35s (c) TI-30X and TI-36X models. Calculators not included in this list are not permitted.

More information

Multidimensional digital signal processing

Multidimensional digital signal processing PSfrag replacements Two-dimensional discrete signals N 1 A 2-D discrete signal (also N called a sequence or array) is a function 2 defined over thex(n set 1 of, n 2 ordered ) pairs of integers: y(nx 1,

More information

14 EE 2402 Engineering Mathematics III Solutions to Tutorial 3 1. For n =0; 1; 2; 3; 4; 5 verify that P n (x) is a solution of Legendre's equation wit

14 EE 2402 Engineering Mathematics III Solutions to Tutorial 3 1. For n =0; 1; 2; 3; 4; 5 verify that P n (x) is a solution of Legendre's equation wit EE 0 Engineering Mathematics III Solutions to Tutorial. For n =0; ; ; ; ; verify that P n (x) is a solution of Legendre's equation with ff = n. Solution: Recall the Legendre's equation from your text or

More information

Lecture 2 OKAN UNIVERSITY FACULTY OF ENGINEERING AND ARCHITECTURE

Lecture 2 OKAN UNIVERSITY FACULTY OF ENGINEERING AND ARCHITECTURE OKAN UNIVERSITY FACULTY OF ENGINEERING AND ARCHITECTURE EEE 43 DIGITAL SIGNAL PROCESSING (DSP) 2 DIFFERENCE EQUATIONS AND THE Z- TRANSFORM FALL 22 Yrd. Doç. Dr. Didem Kivanc Tureli didemk@ieee.org didem.kivanc@okan.edu.tr

More information

Voiced Speech. Unvoiced Speech

Voiced Speech. Unvoiced Speech Digital Speech Processing Lecture 2 Homomorphic Speech Processing General Discrete-Time Model of Speech Production p [ n] = p[ n] h [ n] Voiced Speech L h [ n] = A g[ n] v[ n] r[ n] V V V p [ n ] = u [

More information

Lecture 04: Discrete Frequency Domain Analysis (z-transform)

Lecture 04: Discrete Frequency Domain Analysis (z-transform) Lecture 04: Discrete Frequency Domain Analysis (z-transform) John Chiverton School of Information Technology Mae Fah Luang University 1st Semester 2009/ 2552 Outline Overview Lecture Contents Introduction

More information

Digital Signal Processing, Homework 1, Spring 2013, Prof. C.D. Chung

Digital Signal Processing, Homework 1, Spring 2013, Prof. C.D. Chung Digital Signal Processing, Homework, Spring 203, Prof. C.D. Chung. (0.5%) Page 99, Problem 2.2 (a) The impulse response h [n] of an LTI system is known to be zero, except in the interval N 0 n N. The input

More information

ESE 531: Digital Signal Processing

ESE 531: Digital Signal Processing ESE 531: Digital Signal Processing Lec 6: January 31, 2017 Inverse z-transform Lecture Outline! z-transform " Tie up loose ends " Regions of convergence properties! Inverse z-transform " Inspection " Partial

More information

Test 2 Electrical Engineering Bachelor Module 8 Signal Processing and Communications

Test 2 Electrical Engineering Bachelor Module 8 Signal Processing and Communications Test 2 Electrical Engineering Bachelor Module 8 Signal Processing and Communications (201400432) Tuesday May 26, 2015, 14:00-17:00h This test consists of three parts, corresponding to the three courses

More information

Fourier Series Representation of

Fourier Series Representation of Fourier Series Representation of Periodic Signals Rui Wang, Assistant professor Dept. of Information and Communication Tongji University it Email: ruiwang@tongji.edu.cn Outline The response of LIT system

More information

Discrete-time first-order systems

Discrete-time first-order systems Discrete-time first-order systems 1 Start with the continuous-time system ẏ(t) =ay(t)+bu(t), y(0) Zero-order hold input u(t) =u(nt ), nt apple t

More information

Digital Signal Processing:

Digital Signal Processing: Digital Signal Processing: Mathematical and algorithmic manipulation of discretized and quantized or naturally digital signals in order to extract the most relevant and pertinent information that is carried

More information

Implementation of Discrete-Time Systems

Implementation of Discrete-Time Systems EEE443 Digital Signal Processing Implementation of Discrete-Time Systems Dr. Shahrel A. Suandi PPKEE, Engineering Campus, USM Introduction A linear-time invariant system (LTI) is described by linear constant

More information