Chapter 3 Data Acquisition and Manipulation
|
|
- Fay Rogers
- 5 years ago
- Views:
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.
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 informationSignals 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 information21.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 informationGATE 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 informationProblem 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 informationRecursive, 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 informationDynamic 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 informationEEL3135: 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 informationLecture 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 informationEE482: 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 informationEE482: 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 informationCMPT 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 informationDavid 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 informationLet 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 informationDiscrete 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 informationLecture 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 informationA 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 informationDiscrete 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 informationAnalog 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 informationE : 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 informationDiscrete-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 informationReview 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 informationMAE 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 informationGEORGIA 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 informationGeneralizing 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 informationAutomatique. 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 informationCh. 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 informationLecture 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 informationSignals, 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 informationReview 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 informationUNIT - 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 informationResponses 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 informationControl 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 informationQuestion 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 informationz-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 informationLecture 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 informationFrequency 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 informationIt 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 informationDigital 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 informationReview 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 informationEC 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 informationStep 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 informationGEORGIA 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 informationDr. 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 informationThe 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 informationDynamic 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 informationReview 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 informationSolutions 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 informationLecture 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 informationEE 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 informationCMPT 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 informationDigital 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 informationChapter 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 informationProblem 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 informationAndrea 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 informationEE 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 informationLTI 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 informationUse: 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 informationTransform 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 informationTime 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 informationSinusoids. 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!
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 informationCMPT 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 informationControl 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 informationMethods 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 informationDigital 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 informationContents. 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 informationDr 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!
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 informationTime 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 informationFrequency 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 informationEE 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 informationHomework 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 informationVU 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 informationEET 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 informationDigital 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 informationSome 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 informationTable 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 information2.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 informationModule 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 informationDiscrete-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 informationCM2202: 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 informationChap. 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 informationChapter 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 informationELEC2400 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 informationReview 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 informationZ 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 informationDigital 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 informationMAE143a: 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 informationMAE 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 informationCosc 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 informationShift 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 informationUNIVERSITY 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 informationHow 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 informationDigital 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 informationSAMPLE 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 informationFourier 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 informationLECTURE 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 informationEEE 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 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