Chapter 3 Data Acquisition and Manipulation

Size: px
Start display at page:

Download "Chapter 3 Data Acquisition and Manipulation"

Transcription

1 1 Chapter 3 Data Acquisition and Manipulation In this chapter we introduce z transf orm, or the discrete Laplace Transform, to solve linear recursions. Section 3.1 z-transform Given a data stream x {x 0, x 1, x 2,...}, let X Z (x x k z k x X is called a z transform Example: x {a, a, a,...}, X az (z 1 1 Example: x { 1, a, a 2, a 3,... }, X z (z a 1 Data x i could be a complex number. For instance, the Sine wave x (t sin ωt is sampled every T seconds to yield the signals x k sin (ωt k Im ( e kωt i, k 0, 1, 2,...

2 2 So by the Euler formula X sin (ωt k z k ( e ωt i k 1 e kωt i e kωt i 2iz k ( e ωt i k 1 2i z 2i (1 eωt i 1 1 z (1 e ωt i 1 Linearity Property: 1 2i z 2i 1 z 2i z e 1 z ωt i 2i z e ωt i z sin ωt z 2 2z cos ωt + 1 z Delay of x : x ( 1 {0, x 0, x 1,...} Z (ax + bx az (x + bz (x X 1 Z ( x ( 1 k1 x k 1 z k X z

3 3 Delay of x ( 1 : x ( 2 {0, 0, x 0, x 1,...} ( x ( 1 ( 1 In general Forward of x : x (+1 {x 1, x 2,...} X +1 Z ( x (+1 X 2 Z ( x ( 2 1 z Z ( x ( 1 X z 2 x k+1 z k z Z ( x ( j X z j (1 x k+1 z z k+1 x k z k x 0z z (X x 0 zx x 0 z In general for kth forward of x : x (+k, its Z transform X +k Z ( x (+k X +k z k X x 0 z k x 1 z k 1 x 2 z k 2... x k 1 z For any sequence x {x k }, for convenience we may write Z (x k Z (x.

4 4 For a shift {x k 1 } {x 1, x 0, x 1, x 2,...} of x with given initial value x 1,we write X k 1 Z (x k 1 Z ({x k 1 } x 1 + x 0 z + x 1 z x ( x 0 + x 1 z z +... x 1 + Z (x z For a pth shift {x k p } {x p, x p+1,...x 1, x 0, x 1, x 2,...}with given initial values x p, x p+1,...x 1,we have accordingly X k p Z (x k p Z ( {x k p } x p + x p+1 z x p + x p+1 z x p + x p+1 z x 1 z + x 0 p 1 z + x p zp x 1 z + 1 ( x p 1 z p 0 + x 1 z x 1 z + 1 p 1 zpz (x (2

5 5 For instance, If y k ax k 1 + bx k 2, and Y Z (y k Z (ax k 1 + bx k 2 az (x k 1 + bz (x k 2 ( a x 1 + Z (x ( k + b x 2 + x 1 z z + Z (x k z 2 ax 1 + bx 2 + bx ( 1 a + z z + b Z (x z 2 (3 Section 3.2 Linear Recursions Consider the equation for all k n, n + 1,... x k a 1 x k 1 + a 2 x k a n x k n where a k are fixed constant. x 1, x 2,..., x n are given initially, and are called initial data. We can use z transform to solve x by using the above calculation: (i Apply z transform to both sides. (ii note that Z (x k p is given by (2 (iii Solve X Z (x (iv Recover x k Example (page 34 Solve x k x k 1 + 2x k 2, x 1 1/2, x 2 1/4

6 6 Solution: By (3 with a 1, b 2, we have Z (x k ax 1 + bx 2 + bx 1 z + ( a z + b Z (x z 2 k Set w 1/z, then 1 z + ( 1 z + 2 z 2 Z (x k ( 1 1 z 2 z 2 1 Z (x k 1 z 1 3 (1 2w 1 3 (1 + w 1 (2w k 1 ( w k k ( 1 k 3 1 z k w 1 w 2w w 2 (1 2w (1 + w [ 2 k ( 1 k] w k Ans: x k 2k ( 1 k 3 for k > 0

7 Problem in page 35: x k ax k 1 + b, for k 1, 2,..., x 0 0.7%. Find x k. There are two ways to solve it. (1 Set y k x k+1. So y 1 x 0 0.7% is given, and y k ay k 1 +b for k 0, 1,... Use the relation we derived earlier, Y ay 1 + Z (b a ( z 1 Y + y 1 + ( 1 az 1 Y ay 1 + bz z 1 bz z 1 ay 1 Y 1 az + bz 1 (z 1 (1 az 1 (b + ay 1 z 2 ay 1 z (z 1 (z a (b + ax 0 + A (z 1 + B (z a (2 From x 0 ax 1 + b, we solve x 1 (x 0 b /z. Then, we proceed to solve X (b + ax 1 z 2 ax 1 z (z 1 (z a x 0 + C (z 1 + D (z a x 0z 2 ax 1 z (z 1 (z a 7

8 Convolution: For x {x k }, h {h k }, the discrete convolution product y h x is defined as follows k k y k h 0 x k + h 1 x k h k u 0 h j x k j x j h k j In the frequency domain, Z (h u Z (h Z (u, or Y HU. This is because the above expression indicate power series product. In fact ( ( a k z k b k z k j0 ( a 0 + a 1 z + a 2 z ( b 0 + b 1 z + b 2 z j0 a 0 b 0 + (a 0 b 1 + a 1 b 0 z + (a 0 b 2 + a 1 b 1 + a 2 b 0 z Examples: (1 for δ (0 {1, 0, 0, 0,...}, u δ (0 u, (2 δ (1 {0, 1, 0, 0, 0,...} δ ( 1, u δ (1 u ( 1, (3 u δ (j u ( j Section 3.3 Filters A filter F is a device or algorithm that turns one stream of signals to another more useful one. We assume following three properties (i F is linear: F (au + bv af (u + bf (v (ii F is causal, i.e., the output depends on past and current inputs but not future inputs. In other 8

9 9 words, if y F (u, then u k 0 for all k < k 0 implies y k 0 for all k < k 0. (iii F is time invariant: if y F (u, then the same is true for delay shift, i.e., F { { } u ( k0 } y( k0 In other words, if F : u y, then F : {u k, u k+1,...} {y k, y k+1,...} Let δ {1, 0, 0, 0,...} be the unit impulse, and h F (δ {h 0, h 1,...}.Then for any u y F (u h u where stands for discrete convolution product defined as follows y k h 0 u k + h 1 u k h k u 0 k h j u k j j0 k u j h k j j0 Proof: Let δ (j δ ( j {0,..., 0, 1, 0,...} be a j delay, where 1 is in j th position. Then by

10 property #3, F ( δ (j h ( j {0,..., 0, h 0, h 1,...}Now 10 u u j δ ( j j0. Note that so for any k y k y F (u u j F ( δ ( j u j h ( j j0 ( h( j k 0 if k < j, ( h( j ( u j h( j k u j h k j j0 j0 j0 k h k j if k j k u j h k j (h u k H (z Z (h is called transfer function of F with impulse response h Examples: (1 u δ (j u ( j, (2 Z ( u ( j Z (u δ (j ( UZ δ (j Uz j j0

11 11 Example (page 38 Consider the filter with the finite impulse response (FIR h {1/2, 1/2, 0, 0,...} 1 2 δ( δ(1 For any signal u {u 0, u 1,...}, y u h 1 2 u δ( u δ(1 1 2 u h ( 1. This can also be seen directly y k k j0 u j h k j u k h 0 + u k 1 h 1 u k + u k 1, y 0 u 0 h 0 u So y u + u ( 1, Y U + 1 2z U ( 1 + z 1 U 2 The transfer function is H 1 + z 1 2 Problem (page 38: Suppose that in a room with various sounds. One want to eliminate frequency 60 Hz sound wave. In other words, the component A sin (120πt needs to be eliminated. In theory,

12 to do so, we need to add A sin (120πt A sin (120πt π A sin ( 120π(t The black curve is sine wave of frquency 60 Hz with period 1/60 the red has frequency 720 Hz with period 1/720. If we sample using 720 Hz, i.e., pick up a signal every T 1/720 second, it will meet the max of the black wave at t 1/120 6T, 3/120 18T, 5/120 30T, 42T,..., i.e., every 12T. 12 y x

13 1 scale 1 : 120 6T Let u k signal picked up at time kt. So u 6, u 18, u 30,... will be 60 Hz signal. To cancel it, we use the min point: at the second max point u 18, we use the previous min point u 12 to cancel it. Thus the filter y k u k + u k 6, or Y ( 1 + z 6 U Downside: Green wave has frequency 180 Hz. This filter could also cancel this frequency sound. Section 3.4 Stability We say a filter with impulse response h is stable if bounded inputs u yield bounded outputs y. Theorem: A filter is stable iff the transfer function H (z absolutely converges on the unit circle z 1, i.e., h j < If H has a pole at z z 0 of order k, i.e., H (z G (z j0 (z z 0 k, G (z is bounded near z 0 then it is stable only if z 0 < 1. This is because for all z > 1, z z 0 z z 0 > 1 z 0 > 0. 13

14 14 Thus Example in page 40. H (z G (z G (z is bounded for all z. k k z z 0 (1 z 0 Section 3.5 Polar and Bode Plots We shall introduce two graphing methods to exam efficacy of filters Theorem. Let F be a stable filter with real impulse response h {h 0, h 1,...} and transfer function H (z Z (h.then after transients have died away, the response to the sinusoidal signal u k sin ωkt is also a sinusoid y k r sin (ωkt + φ of the same frequency but different amplitude and phase angle r H ( e iωt, φ arg H ( e iωt. Proof: Recall that y F (u h u. In particular y k k h j u k j j0

15 15 For u k ξ k, k k y k h j ξ k j ξ k h j ξ j j0 ξ k h j ξ j ξ k j0 j0 jk+1 h j ξ j ξ k H (ξ + o (1, where o (1 0 as k (by assumption of after dying away Set then, the response to u k ξ k e iωkt is H (ξ re iφ, ξ e iωt y k ξ k H (ξ + o (1 re iφ ξ k + o (1 re i(ωkt +φ + o (1 r [cos (ωkt + φ + i sin (ωkt + φ] + o (1 Since h k are real, the response to the imaginary parts of e iωkt, i.e., Im (u k sin ωkt, should be the imaginary parts of the response to u k Im (y k r sin (ωkt + φ + o (1. The graph (r (ω, φ (ω in polar coordinate is called Polar plots. T smapling period, 1/T sampling

16 16 rate. Two separate graphs, r r (ω, and φ φ (ω are called Bode plots (Bode-ee. Bode plot for r is often scaled in decibels (db of powers 20 log 10 r (see section 3.8 Example in page 41: h {1/2, 1/2, 0, 0, 0,...}, y k (u k + u k 1 /2, H (z ( 1 + z 1 2. So So Polar plot H ( e iωt 1 + e iωt 2 r 1 + cos ωt, φ arctan cos ωt i sin ωt 2 sin ωt 1 + cos ωt

17 17 y x Bode plots 1.0 r fi (degree w w Section 3.7 Closing the loop

18 18 See chart in page 47. We add a feedback filter H to improve stability For any plant (U HY P Y Y P 1 1 az z 1 z a P 1 + HP U that has a pole at z a. So for a > 1, it is unstable. If we add a constant filter H k. Then that has a pole at Y z P 1 + kp U z (1 + k z a U a 1 + k < 1 if k > a 1 Example (at bottom of page 48: For the planr P ( 1 2z 1 1, a 2. We add filter k 3/2 > 1 will make it stable. Homework: 3.1, 3.2, 3.7, 3.8, 3.18, 3.20, 3.21 Find explicit formula for Fibonacci sequence: F 0 F 1 1, F k+2 F k + F k+1 for k 0, 1, 2,... Then

19 19 write a Matlab routine to varify your answer. Would answer be different if we start with F 0 1, F 1 2, F k+2 F k + F k+1 for k 0, 1, 2,...? Or F 0 2, F 1 3, F k+2 F k + F k+1 for k 0, 1, 2,...? Project: 3.10 (optional

This homework will not be collected or graded. It is intended to help you practice for the final exam. Solutions will be posted.

This homework will not be collected or graded. It is intended to help you practice for the final exam. Solutions will be posted. 6.003 Homework #14 This homework will not be collected or graded. It is intended to help you practice for the final exam. Solutions will be posted. Problems 1. Neural signals The following figure illustrates

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

21.4. Engineering Applications of z-transforms. Introduction. Prerequisites. Learning Outcomes

21.4. Engineering Applications of z-transforms. Introduction. Prerequisites. Learning Outcomes Engineering Applications of z-transforms 21.4 Introduction In this Section we shall apply the basic theory of z-transforms to help us to obtain the response or output sequence for a discrete system. This

More information

GATE EE Topic wise Questions SIGNALS & SYSTEMS

GATE EE Topic wise Questions SIGNALS & SYSTEMS www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)

More information

Problem Set 2: Solution Due on Wed. 25th Sept. Fall 2013

Problem Set 2: Solution Due on Wed. 25th Sept. Fall 2013 EE 561: Digital Control Systems Problem Set 2: Solution Due on Wed 25th Sept Fall 2013 Problem 1 Check the following for (internal) stability [Hint: Analyze the characteristic equation] (a) u k = 05u k

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

Dynamic Response. Assoc. Prof. Enver Tatlicioglu. Department of Electrical & Electronics Engineering Izmir Institute of Technology.

Dynamic Response. Assoc. Prof. Enver Tatlicioglu. Department of Electrical & Electronics Engineering Izmir Institute of Technology. Dynamic Response Assoc. Prof. Enver Tatlicioglu Department of Electrical & Electronics Engineering Izmir Institute of Technology Chapter 3 Assoc. Prof. Enver Tatlicioglu (EEE@IYTE) EE362 Feedback Control

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

Lecture 9 Infinite Impulse Response Filters

Lecture 9 Infinite Impulse Response Filters Lecture 9 Infinite Impulse Response Filters Outline 9 Infinite Impulse Response Filters 9 First-Order Low-Pass Filter 93 IIR Filter Design 5 93 CT Butterworth filter design 5 93 Bilinear transform 7 9

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 05 IIR Design 14/03/04 http://www.ee.unlv.edu/~b1morris/ee482/

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

CMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation

CMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation CMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University January 23, 2006 1 Exponentials The exponential is

More information

David Weenink. First semester 2007

David Weenink. First semester 2007 Institute of Phonetic Sciences University of Amsterdam First semester 2007 Digital s What is a digital filter? An algorithm that calculates with sample values Formant /machine H 1 (z) that: Given input

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

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

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

A system that is both linear and time-invariant is called linear time-invariant (LTI).

A system that is both linear and time-invariant is called linear time-invariant (LTI). The Cooper Union Department of Electrical Engineering ECE111 Signal Processing & Systems Analysis Lecture Notes: Time, Frequency & Transform Domains February 28, 2012 Signals & Systems Signals are mapped

More information

Discrete Systems. Step response and pole locations. Mark Cannon. Hilary Term Lecture

Discrete Systems. Step response and pole locations. Mark Cannon. Hilary Term Lecture Discrete Systems Mark Cannon Hilary Term 22 - Lecture 4 Step response and pole locations 4 - Review Definition of -transform: U() = Z{u k } = u k k k= Discrete transfer function: Y () U() = G() = Z{g k},

More information

Analog vs. discrete signals

Analog vs. discrete signals Analog vs. discrete signals Continuous-time signals are also known as analog signals because their amplitude is analogous (i.e., proportional) to the physical quantity they represent. Discrete-time signals

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

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

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

MAE 143B - Homework 9

MAE 143B - Homework 9 MAE 143B - Homework 9 7.1 a) We have stable first-order poles at p 1 = 1 and p 2 = 1. For small values of ω, we recover the DC gain K = lim ω G(jω) = 1 1 = 2dB. Having this finite limit, our straight-line

More information

GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM. COURSE: ECE 3084A (Prof. Michaels)

GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM. COURSE: ECE 3084A (Prof. Michaels) GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM DATE: 30-Apr-14 COURSE: ECE 3084A (Prof. Michaels) NAME: STUDENT #: LAST, FIRST Write your name on the front page

More information

Generalizing the DTFT!

Generalizing the DTFT! The Transform Generaliing the DTFT! The forward DTFT is defined by X e jω ( ) = x n e jωn in which n= Ω is discrete-time radian frequency, a real variable. The quantity e jωn is then a complex sinusoid

More information

Automatique. A. Hably 1. Commande d un robot mobile. Automatique. A.Hably. Digital implementation

Automatique. A. Hably 1. Commande d un robot mobile. Automatique. A.Hably. Digital implementation A. Hably 1 1 Gipsa-lab, Grenoble-INP ahmad.hably@grenoble-inp.fr Commande d un robot mobile (Gipsa-lab (DA)) ASI 1 / 25 Outline 1 2 (Gipsa-lab (DA)) ASI 2 / 25 of controllers Signals must be sampled and

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 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

Signals, Instruments, and Systems W5. Introduction to Signal Processing Sampling, Reconstruction, and Filters

Signals, Instruments, and Systems W5. Introduction to Signal Processing Sampling, Reconstruction, and Filters Signals, Instruments, and Systems W5 Introduction to Signal Processing Sampling, Reconstruction, and Filters Acknowledgments Recapitulation of Key Concepts from the Last Lecture Dirac delta function (

More information

Review of Fundamentals of Digital Signal Processing

Review of Fundamentals of Digital Signal Processing Solution Manual for Theory and Applications of Digital Speech Processing by Lawrence Rabiner and Ronald Schafer Click here to Purchase full Solution Manual at http://solutionmanuals.info Link download

More information

UNIT - III PART A. 2. Mention any two techniques for digitizing the transfer function of an analog filter?

UNIT - III PART A. 2. Mention any two techniques for digitizing the transfer function of an analog filter? UNIT - III PART A. Mention the important features of the IIR filters? i) The physically realizable IIR filters does not have linear phase. ii) The IIR filter specification includes the desired characteristics

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

Control Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich

Control Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich Control Systems I Lecture 6: Poles and Zeros Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 27, 2017 E. Frazzoli (ETH) Lecture 6: Control Systems I 27/10/2017

More information

Question Paper Code : AEC11T02

Question Paper Code : AEC11T02 Hall Ticket No Question Paper Code : AEC11T02 VARDHAMAN COLLEGE OF ENGINEERING (AUTONOMOUS) Affiliated to JNTUH, Hyderabad Four Year B. Tech III Semester Tutorial Question Bank 2013-14 (Regulations: VCE-R11)

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

Lecture 8 Finite Impulse Response Filters

Lecture 8 Finite Impulse Response Filters Lecture 8 Finite Impulse Response Filters Outline 8. Finite Impulse Response Filters.......................... 8. oving Average Filter............................... 8.. Phase response...............................

More information

Frequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability

Frequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability Lecture 6. Loop analysis of feedback systems 1. Motivation 2. Graphical representation of frequency response: Bode and Nyquist curves 3. Nyquist stability theorem 4. Stability margins Frequency methods

More information

It is common to think and write in time domain. creating the mathematical description of the. Continuous systems- using Laplace or s-

It is common to think and write in time domain. creating the mathematical description of the. Continuous systems- using Laplace or s- It is common to think and write in time domain quantities, but this is not the best thing to do in creating the mathematical description of the system we are dealing with. Continuous systems- using Laplace

More information

Digital Signal Processing

Digital Signal Processing COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #21 Friday, October 24, 2003 Types of causal FIR (generalized) linear-phase filters: Type I: Symmetric impulse response: with order M an even

More information

Review of Fundamentals of Digital Signal Processing

Review of Fundamentals of Digital Signal Processing Chapter 2 Review of Fundamentals of Digital Signal Processing 2.1 (a) This system is not linear (the constant term makes it non linear) but is shift-invariant (b) This system is linear but not shift-invariant

More information

EC CONTROL SYSTEM UNIT I- CONTROL SYSTEM MODELING

EC CONTROL SYSTEM UNIT I- CONTROL SYSTEM MODELING EC 2255 - CONTROL SYSTEM UNIT I- CONTROL SYSTEM MODELING 1. What is meant by a system? It is an arrangement of physical components related in such a manner as to form an entire unit. 2. List the two types

More information

Step Response Analysis. Frequency Response, Relation Between Model Descriptions

Step Response Analysis. Frequency Response, Relation Between Model Descriptions Step Response Analysis. Frequency Response, Relation Between Model Descriptions Automatic Control, Basic Course, Lecture 3 November 9, 27 Lund University, Department of Automatic Control Content. Step

More information

GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM. COURSE: ECE 3084A (Prof. Michaels)

GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM. COURSE: ECE 3084A (Prof. Michaels) GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM DATE: 09-Dec-13 COURSE: ECE 3084A (Prof. Michaels) NAME: STUDENT #: LAST, FIRST Write your name on the front page

More information

Dr. Ian R. Manchester

Dr. Ian R. Manchester Dr Ian R. Manchester Week Content Notes 1 Introduction 2 Frequency Domain Modelling 3 Transient Performance and the s-plane 4 Block Diagrams 5 Feedback System Characteristics Assign 1 Due 6 Root Locus

More information

The loop shaping paradigm. Lecture 7. Loop analysis of feedback systems (2) Essential specifications (2)

The loop shaping paradigm. Lecture 7. Loop analysis of feedback systems (2) Essential specifications (2) Lecture 7. Loop analysis of feedback systems (2). Loop shaping 2. Performance limitations The loop shaping paradigm. Estimate performance and robustness of the feedback system from the loop transfer L(jω)

More information

Dynamic circuits: Frequency domain analysis

Dynamic circuits: Frequency domain analysis Electronic Circuits 1 Dynamic circuits: Contents Free oscillation and natural frequency Transfer functions Frequency response Bode plots 1 System behaviour: overview 2 System behaviour : review solution

More information

Review of Frequency Domain Fourier Series: Continuous periodic frequency components

Review of Frequency Domain Fourier Series: Continuous periodic frequency components Today we will review: Review of Frequency Domain Fourier series why we use it trig form & exponential form how to get coefficients for each form Eigenfunctions what they are how they relate to LTI systems

More information

Solutions to Problems in Chapter 4

Solutions to Problems in Chapter 4 Solutions to Problems in Chapter 4 Problems with Solutions Problem 4. Fourier Series of the Output Voltage of an Ideal Full-Wave Diode Bridge Rectifier he nonlinear circuit in Figure 4. is a full-wave

More information

Lecture 16: Filter Design: Impulse Invariance and Bilinear Transform

Lecture 16: Filter Design: Impulse Invariance and Bilinear Transform EE58 Digital Signal Processing University of Washington Autumn 2 Dept. of Electrical Engineering Lecture 6: Filter Design: Impulse Invariance and Bilinear Transform Nov 26, 2 Prof: J. Bilmes

More information

EE 3054: Signals, Systems, and Transforms Spring A causal discrete-time LTI system is described by the equation. y(n) = 1 4.

EE 3054: Signals, Systems, and Transforms Spring A causal discrete-time LTI system is described by the equation. y(n) = 1 4. EE : Signals, Systems, and Transforms Spring 7. A causal discrete-time LTI system is described by the equation Test y(n) = X x(n k) k= No notes, closed book. Show your work. Simplify your answers.. A discrete-time

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

Chapter 7. Digital Control Systems

Chapter 7. Digital Control Systems Chapter 7 Digital Control Systems 1 1 Introduction In this chapter, we introduce analysis and design of stability, steady-state error, and transient response for computer-controlled systems. Transfer functions,

More information

Problem Weight Score Total 100

Problem Weight Score Total 100 EE 350 EXAM IV 15 December 2010 Last Name (Print): First Name (Print): ID number (Last 4 digits): Section: DO NOT TURN THIS PAGE UNTIL YOU ARE TOLD TO DO SO Problem Weight Score 1 25 2 25 3 25 4 25 Total

More information

Andrea Zanchettin Automatic Control AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Linear systems (frequency domain)

Andrea Zanchettin Automatic Control AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Linear systems (frequency domain) 1 AUTOMATIC CONTROL Andrea M. Zanchettin, PhD Spring Semester, 2018 Linear systems (frequency domain) 2 Motivations Consider an LTI system Thanks to the Lagrange s formula we can compute the motion of

More information

EE Experiment 11 The Laplace Transform and Control System Characteristics

EE Experiment 11 The Laplace Transform and Control System Characteristics EE216:11 1 EE 216 - Experiment 11 The Laplace Transform and Control System Characteristics Objectives: To illustrate computer usage in determining inverse Laplace transforms. Also to determine useful signal

More information

LTI Systems (Continuous & Discrete) - Basics

LTI Systems (Continuous & Discrete) - Basics LTI Systems (Continuous & Discrete) - Basics 1. A system with an input x(t) and output y(t) is described by the relation: y(t) = t. x(t). This system is (a) linear and time-invariant (b) linear and time-varying

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

Transform analysis of LTI systems Oppenheim and Schafer, Second edition pp For LTI systems we can write

Transform analysis of LTI systems Oppenheim and Schafer, Second edition pp For LTI systems we can write Transform analysis of LTI systems Oppenheim and Schafer, Second edition pp. 4 9. For LTI systems we can write yœn D xœn hœn D X kd xœkhœn Alternatively, this relationship can be expressed in the z-transform

More information

Time Series Analysis

Time Series Analysis Time Series Analysis Informatics and Mathematical Modelling Technical University of Denmark DK-2800 Kgs. Lyngby 1 Outline of todays lecture Descriptions of (deterministic) linear systems. Chapter 4: Linear

More information

Sinusoids. Amplitude and Magnitude. Phase and Period. CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation

Sinusoids. Amplitude and Magnitude. Phase and Period. CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation Sinusoids CMPT 889: Lecture Sinusoids, Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University September 6, 005 Sinusoids are

More information

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 9. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -

More information

CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation

CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation CMPT 889: Lecture 2 Sinusoids, Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University September 26, 2005 1 Sinusoids Sinusoids

More information

Control for. Maarten Steinbuch Dept. Mechanical Engineering Control Systems Technology Group TU/e

Control for. Maarten Steinbuch Dept. Mechanical Engineering Control Systems Technology Group TU/e Control for Maarten Steinbuch Dept. Mechanical Engineering Control Systems Technology Group TU/e Motion Systems m F Introduction Timedomain tuning Frequency domain & stability Filters Feedforward Servo-oriented

More information

Methods for analysis and control of. Lecture 6: Introduction to digital control

Methods for analysis and control of. Lecture 6: Introduction to digital control Methods for analysis and of Lecture 6: to digital O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.inpg.fr www.lag.ensieg.inpg.fr/sename 6th May 2009 Outline Some interesting books:

More information

Digital Filters Ying Sun

Digital Filters Ying Sun Digital Filters Ying Sun Digital filters Finite impulse response (FIR filter: h[n] has a finite numbers of terms. Infinite impulse response (IIR filter: h[n] has infinite numbers of terms. Causal filter:

More information

Contents. PART I METHODS AND CONCEPTS 2. Transfer Function Approach Frequency Domain Representations... 42

Contents. PART I METHODS AND CONCEPTS 2. Transfer Function Approach Frequency Domain Representations... 42 Contents Preface.............................................. xiii 1. Introduction......................................... 1 1.1 Continuous and Discrete Control Systems................. 4 1.2 Open-Loop

More information

Dr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Review

Dr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Review Week Date Content Notes 1 6 Mar Introduction 2 13 Mar Frequency Domain Modelling 3 20 Mar Transient Performance and the s-plane 4 27 Mar Block Diagrams Assign 1 Due 5 3 Apr Feedback System Characteristics

More information

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!

Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 3.. 24 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -

More information

Time Response Analysis (Part II)

Time Response Analysis (Part II) Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary

More information

Frequency response. Pavel Máša - XE31EO2. XE31EO2 Lecture11. Pavel Máša - XE31EO2 - Frequency response

Frequency response. Pavel Máša - XE31EO2. XE31EO2 Lecture11. Pavel Máša - XE31EO2 - Frequency response Frequency response XE3EO2 Lecture Pavel Máša - Frequency response INTRODUCTION Frequency response describe frequency dependence of output to input voltage magnitude ratio and its phase shift as a function

More information

EE Control Systems LECTURE 9

EE Control Systems LECTURE 9 Updated: Sunday, February, 999 EE - Control Systems LECTURE 9 Copyright FL Lewis 998 All rights reserved STABILITY OF LINEAR SYSTEMS We discuss the stability of input/output systems and of state-space

More information

Homework 7 - Solutions

Homework 7 - Solutions Homework 7 - Solutions Note: This homework is worth a total of 48 points. 1. Compensators (9 points) For a unity feedback system given below, with G(s) = K s(s + 5)(s + 11) do the following: (c) Find the

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

EET 3212 Control Systems. Control Systems Engineering, 6th Edition, Norman S. Nise December 2010, A. Goykadosh and M.

EET 3212 Control Systems. Control Systems Engineering, 6th Edition, Norman S. Nise December 2010, A. Goykadosh and M. NEW YORK CITY COLLEGE OF TECHNOLOGY The City University of New York 300 Jay Street Brooklyn, NY 11201-2983 Department of Electrical and Telecommunications Engineering Technology TEL (718) 260-5300 - FAX:

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

Some solutions of the written exam of January 27th, 2014

Some solutions of the written exam of January 27th, 2014 TEORIA DEI SISTEMI Systems Theory) Prof. C. Manes, Prof. A. Germani Some solutions of the written exam of January 7th, 0 Problem. Consider a feedback control system with unit feedback gain, with the following

More information

Table of Laplacetransform

Table of Laplacetransform Appendix Table of Laplacetransform pairs 1(t) f(s) oct), unit impulse at t = 0 a, a constant or step of magnitude a at t = 0 a s t, a ramp function e- at, an exponential function s + a sin wt, a sine fun

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

Module 4 : Laplace and Z Transform Problem Set 4

Module 4 : Laplace and Z Transform Problem Set 4 Module 4 : Laplace and Z Transform Problem Set 4 Problem 1 The input x(t) and output y(t) of a causal LTI system are related to the block diagram representation shown in the figure. (a) Determine a differential

More information

Discrete-Time David Johns and Ken Martin University of Toronto

Discrete-Time David Johns and Ken Martin University of Toronto Discrete-Time David Johns and Ken Martin University of Toronto (johns@eecg.toronto.edu) (martin@eecg.toronto.edu) University of Toronto 1 of 40 Overview of Some Signal Spectra x c () t st () x s () t xn

More information

CM2202: Scientific Computing and Multimedia Applications General Maths: 3. Complex Numbers

CM2202: Scientific Computing and Multimedia Applications General Maths: 3. Complex Numbers CM2202: Scientific Computing and Multimedia Applications General Maths: 3. Complex Numbers Prof. David Marshall School of Computer Science & Informatics A problem when solving some equations There are

More information

Chap. 3 Laplace Transforms and Applications

Chap. 3 Laplace Transforms and Applications Chap 3 Laplace Transforms and Applications LS 1 Basic Concepts Bilateral Laplace Transform: where is a complex variable Region of Convergence (ROC): The region of s for which the integral converges Transform

More information

Chapter 1 Fundamental Concepts

Chapter 1 Fundamental Concepts Chapter 1 Fundamental Concepts 1 Signals A signal is a pattern of variation of a physical quantity, often as a function of time (but also space, distance, position, etc). These quantities are usually the

More information

ELEC2400 Signals & Systems

ELEC2400 Signals & Systems ELEC2400 Signals & Systems Chapter 7. Z-Transforms Brett Ninnes brett@newcastle.edu.au. School of Electrical Engineering and Computer Science The University of Newcastle Slides by Juan I. Yu (jiyue@ee.newcastle.edu.au

More information

Review of Linear Time-Invariant Network Analysis

Review of Linear Time-Invariant Network Analysis D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x

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

Digital Filter Structures. Basic IIR Digital Filter Structures. of an LTI digital filter is given by the convolution sum or, by the linear constant

Digital Filter Structures. Basic IIR Digital Filter Structures. of an LTI digital filter is given by the convolution sum or, by the linear constant Digital Filter Chapter 8 Digital Filter Block Diagram Representation Equivalent Basic FIR Digital Filter Basic IIR Digital Filter. Block Diagram Representation In the time domain, the input-output relations

More information

MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions

MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions Question 1. SIGNALS: Design of a noise-cancelling headphone system. 1a. Based on the low-pass filter given, design a high-pass filter,

More information

MAE 143B - Homework 7

MAE 143B - Homework 7 MAE 143B - Homework 7 6.7 Multiplying the first ODE by m u and subtracting the product of the second ODE with m s, we get m s m u (ẍ s ẍ i ) + m u b s (ẋ s ẋ u ) + m u k s (x s x u ) + m s b s (ẋ s ẋ u

More information

Cosc 3451 Signals and Systems. What is a system? Systems Terminology and Properties of Systems

Cosc 3451 Signals and Systems. What is a system? Systems Terminology and Properties of Systems Cosc 3451 Signals and Systems Systems Terminology and Properties of Systems What is a system? an entity that manipulates one or more signals to yield new signals (often to accomplish a function) can be

More information

Shift Property of z-transform. Lecture 16. More z-transform (Lathi 5.2, ) More Properties of z-transform. Convolution property of z-transform

Shift Property of z-transform. Lecture 16. More z-transform (Lathi 5.2, ) More Properties of z-transform. Convolution property of z-transform Shift Property of -Transform If Lecture 6 More -Transform (Lathi 5.2,5.4-5.5) then which is delay causal signal by sample period. If we delay x[n] first: Peter Cheung Department of Electrical & Electronic

More information

UNIVERSITY OF CALIFORNIA AT BERKELEY College of Engineering Department of Electrical Engineering and Computer Sciences. EE105 Lab Experiments

UNIVERSITY OF CALIFORNIA AT BERKELEY College of Engineering Department of Electrical Engineering and Computer Sciences. EE105 Lab Experiments UNIVERSITY OF CALIFORNIA AT BERKELEY College of Engineering Department of Electrical Engineering and Computer Sciences EE15 Lab Experiments Bode Plot Tutorial Contents 1 Introduction 1 2 Bode Plots Basics

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

Digital Signal Processing

Digital Signal Processing COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #24 Tuesday, November 4, 2003 6.8 IIR Filter Design Properties of IIR Filters: IIR filters may be unstable Causal IIR filters with rational system

More information

SAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015

SAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015 FACULTY OF ENGINEERING AND SCIENCE SAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015 Lecturer: Michael Ruderman Problem 1: Frequency-domain analysis and control design (15 pt) Given is a

More information

Fourier Methods in Digital Signal Processing Final Exam ME 579, Spring 2015 NAME

Fourier Methods in Digital Signal Processing Final Exam ME 579, Spring 2015 NAME Fourier Methods in Digital Signal Processing Final Exam ME 579, Instructions for this CLOSED BOOK EXAM 2 hours long. Monday, May 8th, 8-10am in ME1051 Answer FIVE Questions, at LEAST ONE from each section.

More information

LECTURE 12 Sections Introduction to the Fourier series of periodic signals

LECTURE 12 Sections Introduction to the Fourier series of periodic signals Signals and Systems I Wednesday, February 11, 29 LECURE 12 Sections 3.1-3.3 Introduction to the Fourier series of periodic signals Chapter 3: Fourier Series of periodic signals 3. Introduction 3.1 Historical

More information

EEE 188: Digital Control Systems

EEE 188: Digital Control Systems EEE 88: Digital Control Systems Lecture summary # the controlled variable. Example: cruise control. In feedback control, sensors and measurements play an important role. In discrete time systems, the control

More information

( ) John A. Quinn Lecture. ESE 531: Digital Signal Processing. Lecture Outline. Frequency Response of LTI System. Example: Zero on Real Axis

( ) John A. Quinn Lecture. ESE 531: Digital Signal Processing. Lecture Outline. Frequency Response of LTI System. Example: Zero on Real Axis John A. Quinn Lecture ESE 531: Digital Signal Processing Lec 15: March 21, 2017 Review, Generalized Linear Phase Systems Penn ESE 531 Spring 2017 Khanna Lecture Outline!!! 2 Frequency Response of LTI System

More information