Methods for analysis and control of. Lecture 6: Introduction to digital control
|
|
- Augustine Hodges
- 6 years ago
- Views:
Transcription
1 Methods for analysis and of Lecture 6: to digital O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.inpg.fr 6th May 2009
2 Outline
3 Some interesting books: K.J. Astrom and B. Wittenmark, Computer-Controlled Systems, Information and sciences series. Prentice Hall, New Jersey, 3rd edition, R.C. Dorf and R.H. Bishop, Modern Control Systems, Prentice Hall, USA, G.C. Goodwin, S.F. Graebe, and M.E. Salgado, Control System Design, Prentice Hall, New Jersey, G. Franklin, J. Powell, A. Emami-Naeini, Feedback Dynamic Systems, Prentice Hall, 2005
4 Toward digital Digital Usually lers are implemented in a digital computer as: This requires the use of the discrete theory. (Sampling theory + )
5
6 Definitions Mathematical definition Because the output of the ideal sampler, x (t), is a series of impulses with values x(kt e ), we have: x (t) = k=0 by using the Laplace transform, L [x (t)] = x(kt e )δ(t kt e ) k=0 x(kt e )e kst e Noting z = e st e, we can derive the so called X(z) = Z [x(k)] = k=0 x(k)z k
7 Properties Definition Properties X(z) = Z [x(k)] = k=0 x(k)z k Z [αx(k) + βy(k)] = αx(z) + βy (z) Z [x(k n)] = z n Z [x(k)] Z [kx(k)] = z d dz Z [x(k)] Z [x(k) y(k)] = X(z).Y (z) lim x(k) = lim (z 1)X(z) k 1 z 1 The z 1 can be interpreted as a pure delay operator.
8 Exercise Determine the of the step function (1) and of the ramp function (2) x step (k) = 1 if k 0 x ramp (k) = k if k 0 = 0 if k < 0 = 0 if k < 0
9 Exercise Determine the of the step function (1) and of the ramp function (2) x step (k) = 1 if k 0 x ramp (k) = k if k 0 = 0 if k < 0 = 0 if k < 0 Solution 1) Step X step (z) = 1 + z 1 + z 2 + = 2) Ramp (note that x ramp (k) = kx step (k)) 1 1 z 1 = z z 1 X ramp (z) = z d ( z ) dz z 1 z = (z 1) 2
10 Zero order holder Sampler and Zero order holder A sampler is a switch that close every T e seconds. A Zero order holder holds the signal x for T e seconds to get h as: h(t + kt e ) = x(kt e ), 0 t < T e
11 Zero order holder (cont d) Model of the Zero order holder The transfer function of the zero-order holder is given by: G BOZ (s) = 1 s e st e s = 1 e st e s Influence of the D/A and A/D Note that the precision is also limited by the available precision of the converters (either A/D or D/A). This error is also called the amplitude quantization error.
12 Representation of the discrete linear The discrete output of a system can be expressed as: y(k) = n=0 h(k n)u(n) hence, applying the Z-transform leads to Y (z) = Z [h(k)]u(z) = H(z)U(z) H(z) = b 0 + b 1 z + + b m z m a 0 + a 1 z + + a n z n where n ( m) is the order of the system Corresponding difference equation: = Y U y(k) = 1 a n [ b0 u(k n) + b 1 u(k n + 1) + + b m u(k n + m) a 0 y(k n) a 2 y(k n + 1) a n 1 y(k 1) ]
13 Some useful transformations x(t) X(s) X(z) δ(t) 1 1 δ(t kt e ) e kst e z k u(t) 1 s z 1 t s 2 e at s+a 1 1 e at 1 sin(ωt) cos(ωt) z 1 zt e (z 1) 2 z z e at e z(1 e at e ) s(s+a) (z 1)(z e at e ) ω zsin(ωt e ) s 2 +ω 2 z 2 2zcos(ωT e )+1 s z(z cos(ωt e )) s 2 +ω 2 z 2 2zcos(ωT e )+1 Exercise Discretize (sampling time T e ) the system described by the Laplace function (using a Zero order holder): H(s) = Y (s) U(s) = 1 s(s + 1)
14 Exercise Discretize the system described by the Laplace function (using a Zero order holder): H(s) = Y (s) U(s) = 1 s(s + 1) Adding the Zero order holder leads to: hence G BOZ (s)h(s) = 1 e st e s = 1 e st e s 2 (s + 1) 1 s(s + 1) = (1 e st e ) ( 1 s 2 1 s + 1 s + 1 Z [G BOZ (s)h(s)] = (1 z 1 )Z [ 1 s 2 1 s + 1 ] s + 1 )
15 Exercise (cont d) Z [G BOZ (s)h(s)] = (1 z 1 )Z [ 1 s 2 1 s + 1 ] s + 1 = (1 z 1 ) [ zt e (z 1) 2 z z 1 + z ] z e T e if T e = 1, we have = (ze T e z + zt e ) + (1 e T e T e e T e ) (z 1)(z e T e ) Z [G BOZ (s)h(s)] = (ze T e z + zt e ) + (1 e T e T e e T e ) (z 1)(z e T e ) = ze e 1 (z 1)(z e 1 ) b 1 z + b 0 = z 2 + a 1 z + a 0
16 Exercise (cont d) Let us return back to sampled-time domain Y (z) U(z) = Y (z) = b 1 z+b 0 z 2 +a 1 z+a 0 b 1 z+b 0 z 2 +a 1 z+a 0 U(z) Y (z)(z 2 + a 1 z + a 0 ) = (b 1 z + b 0 )U(z) y(n + 2) + a 1 y(n + 1) + a 0 y(n) = b 1 u(n + 1) + b 0 u(n) With an unit feedback, the closed loop function is given by: F cl (z) = G(z) 1 + G(z)
17
18 Equivalence {s} {z} {s} {z} The equivalence between the Laplace domain and the Z domain is obtained by the following transformation: z = e st e Two poles with a imaginary part witch differs of 2π/T e give the same pole in Z. domain
19 Approximations Forward difference (Rectangle inferior) s = z 1 T e Backward difference (Rectangle superior) s = z 1 zt e
20 Approximations (cont d) Trapezoidal difference (Tustin) s = 2 T e z 1 z + 1
21 Systems definition A state space system is as follows: { x((k + 1)h) = Ad x(kh) + B d u(kh), x(0) = x 0 y(kh) = C d x(kh) + D d u(kh) where h is the sampling period. Matlab : ss(a d,b d,c d,d d,h) creates a SS object SYS representing a state-space model (1)
22 Relation with transfer function For, { x((k + 1)h) = Ad x(kh) + B d u(kh), x(0) = x 0 y(kh) = C d x(kh) + D d u(kh) the discrete transfer function is given by (2) G(z) = C d (zi n A d ) 1 B d + D d (3) where z is the shift operator, i.e. zx(kh) = x((k + 1)h)
23 Recall Laplace & Z-transform From Transfer Function to State Space H(s) to state space X U = den(s) Y X = num(s) Ẋ = AX + BU Y = CX + DU H(z) to state space X U = den(z) Y X = num(z) X k+1 = FX k + GU k Y k = CX k + DU k Y (s) = [ C[sI A] 1 B + D ] U(s) Y (z) = [ C[zI F] 1 G + D ] U(z) }{{}}{{} H(s) H(z)
24 About sampling period and time response Influence of the sampling period on the time response Impose a maximal time response to a discrete system is equivalent to place the poles inside a circle defined by the upper bound of the bound given by this time response. The more the poles are close to zero, the more the system is fast.
25 Frequency analysis As in the continuous time, the Bode diagram can also be used. Example with sampling Time T e = 1s w e = 2π): Note that, in our case, the Bode is cut at the pulse w = π. see SYSD = c2d(sysc,ts,method) in MATLAB.
26 Frequency analysis As in the continuous time, the Bode diagram can also be used. Example with sampling Time T e = 1s w e = 2π): Note that, in our case, the Bode is cut at the pulse w = π. see SYSD = c2d(sysc,ts,method) in MATLAB. Sampling Limitations Recall the Shannon theorem that impose the sampling frequency at least 2 times higher that the system maximum frequency. Related to the anti-aliasing filter...
27 About sampling period and robustness Influence of the sampling period on the poles In theory, smaller the sampling period T e is, closer the discrete system is from the continuous one. But reducing the sampling time modify poles location... Poles and zeros become closer to the limit of the unit circle can introduce instability (decrease robustness). Sampling influences stability and robustness Over sampling increase noise sensitivity
28 Zeros Influence of the sampling period on the poles A discrete system with one or few zeros at the origin is faster than one with no origin zeros. In the time domain a zero at the origin induces a sample advance.
29 Recall A linear continuous feedback system is stable if all poles of the closed-loop transfer function T (s) lie in the left half s-plane. The Z-plane is related to the S-plane by z = e st e = e (σ+jω)t e. Hence z = e σt e and z = ωt e
30 (cont d) Jury criteria The denominator polynomial (den(z) = a 0 z n + a 1 z n a n = 0) has all its roots inside the unit circle if all the first coefficients of the odd row are positive. a n b 1 a 0 a 1 a 2... a n k... a 0 = a 0 a n n a 0 2 a n a n 1 a n 2... a k... a 0 a n b 3 b 0 b 1 b 2... b 1 = a 1 a n 1 n 1 2 b n 1 b n 2 b n 3... b 0.. 2n + 1 s 0 b k = a k a n k a n a 0 a 0 c k = b k b n 1 k b n 1 b 0
31 Example Find the stability region of D(z) = z 2 + a 1 z + a 2
32 Example Find the stability region of D(z) = z 2 + a 1 z + a 2 Solution 1 1 > 0 a 1 a 2 2 a 2 a a2 2 > 0? a 1 a 1 a 2 4 a 1 a 1 a 2 1 a2 2 > 0? 5 (1 a2 2 )2 (a 2 1 (1 a 2) 2 ) 1 a 2 2 hence, 1 a 2 2 > 0 (1 + a 2 ) 2 > a 2 1
33 How to get a discrete ler First way Obtain a plant model (by discretization) Design a ler Derive the difference equation Second way Design a continuous-time ler Converse the continuous-time ler to discrete time (c2d) Derive the difference equation Now the question is how to implement the computed ler on a real-time (embedded) system, and what are the precautions to take before?
34
35 Anti-aliasing & Sampling Anti-aliasing Practically it is smart to use a constant high sampling frequency with an analog filter matching this frequency. Then, after the A/D converter, the signal is down-sampled to the frequency used by the ler. Remember that the pre-filter introduce phase shift. Sampling frequency choice The sampling time for are based on the desired speed of the closed loop system. A rule of thumb is that one should sample 4 10 times per rise time T r of the closed loop system. N sample = T r T e 4 10 where T e is the sampling period, and N sample the number of samples.
36 Delay Problematic Sampled theory assume presence of clock that synchronizes all measurements and signal. Hence in a computer based there always is delays ( delay, computational delay, I/O latency). Origins There are several reasons for delay apparition Execution time (code) Preemption from higher order process Interrupt Communication delay Data dependencies Hence the delay is not constant. The delay introduce a phase shift Instability!
37 Delay (cont d) Admissible delay (Bode) Measure the phase margin: PM = ϕ w0 [], where ϕ w0 is the phase at the crossover frequency w 0, i.e. G(jw 0 ) = 1 Then the delay margin is DM = PMπ 180w 0 [s]
38 : some key issues Analysis of allows to account for implementation constraints (network-induced delays, CPU resource availability, coding,...) Digital vs Continuous : how to choose the ad hoc method? Design method can be systematic whatever the framework is.
39
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 informationPole placement control: state space and polynomial approaches Lecture 1
: state space and polynomial approaches Lecture 1 dynamical O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.fr www.gipsa-lab.fr/ o.sename November 7, 2017 Outline dynamical dynamical
More informationModelling, analysis and control of linear systems using state space representations
, analysis and of linear systems using state space representations Grenoble INP / GIPSA-lab Pole placement February 2018 -based optimal digital of dynamical systems as state space representations of the
More informationMethods for analysis and control of. Lecture 4: The root locus design method
Methods for analysis and control of Lecture 4: The root locus design method O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.inpg.fr www.lag.ensieg.inpg.fr/sename Lead Lag 17th March
More informationMethods for analysis and control of dynamical systems Lecture 4: The root locus design method
Methods for analysis and control of Lecture 4: The root locus design method O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.inpg.fr www.gipsa-lab.fr/ o.sename 5th February 2015 Outline
More informationModelling, analysis and control of linear systems using state space representations
Modelling, analysis and of linear systems using state space O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.fr www.gipsa-lab.fr/ o.sename 12th January 2015 feedback optimal digital
More informationControl System Design
ELEC4410 Control System Design Lecture 19: Feedback from Estimated States and Discrete-Time Control Design Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science
More informationIdentification Methods for Structural Systems
Prof. Dr. Eleni Chatzi System Stability Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can be defined from
More informationOutline. Classical Control. Lecture 1
Outline Outline Outline 1 Introduction 2 Prerequisites Block diagram for system modeling Modeling Mechanical Electrical Outline Introduction Background Basic Systems Models/Transfers functions 1 Introduction
More informationProblem Set 3: Solution Due on Mon. 7 th Oct. in class. Fall 2013
EE 56: Digital Control Systems Problem Set 3: Solution Due on Mon 7 th Oct in class Fall 23 Problem For the causal LTI system described by the difference equation y k + 2 y k = x k, () (a) By first finding
More informationDiscrete-time Controllers
Schweizerische Gesellschaft für Automatik Association Suisse pour l Automatique Associazione Svizzera di Controllo Automatico Swiss Society for Automatic Control Advanced Control Discrete-time Controllers
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 informationPole placement control: state space and polynomial approaches Lecture 2
: state space and polynomial approaches Lecture 2 : a state O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.fr www.gipsa-lab.fr/ o.sename -based November 21, 2017 Outline : a state
More informationCONTROL OF DIGITAL SYSTEMS
AUTOMATIC CONTROL AND SYSTEM THEORY CONTROL OF DIGITAL SYSTEMS Gianluca Palli Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI) Università di Bologna Email: gianluca.palli@unibo.it
More informationChapter 13 Digital Control
Chapter 13 Digital Control Chapter 12 was concerned with building models for systems acting under digital control. We next turn to the question of control itself. Topics to be covered include: why one
More informationClassify a transfer function to see which order or ramp it can follow and with which expected error.
Dr. J. Tani, Prof. Dr. E. Frazzoli 5-059-00 Control Systems I (Autumn 208) Exercise Set 0 Topic: Specifications for Feedback Systems Discussion: 30.. 208 Learning objectives: The student can grizzi@ethz.ch,
More informationEECS C128/ ME C134 Final Thu. May 14, pm. Closed book. One page, 2 sides of formula sheets. No calculators.
Name: SID: EECS C28/ ME C34 Final Thu. May 4, 25 5-8 pm Closed book. One page, 2 sides of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 4 2 4 3 6 4 8 5 3
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014
Prof. Dr. Eleni Chatzi System Stability - 26 March, 24 Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can
More informationDigital Control Systems
Digital Control Systems Lecture Summary #4 This summary discussed some graphical methods their use to determine the stability the stability margins of closed loop systems. A. Nyquist criterion Nyquist
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 informationImplementation Issues for the Virtual Spring
Implementation Issues for the Virtual Spring J. S. Freudenberg EECS 461 Embedded Control Systems 1 Introduction One of the tasks in Lab 4 is to attach the haptic wheel to a virtual reference position with
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 informationUniversity of Alberta ENGM 541: Modeling and Simulation of Engineering Systems Laboratory #7. M.G. Lipsett & M. Mashkournia 2011
ENG M 54 Laboratory #7 University of Alberta ENGM 54: Modeling and Simulation of Engineering Systems Laboratory #7 M.G. Lipsett & M. Mashkournia 2 Mixed Systems Modeling with MATLAB & SIMULINK Mixed systems
More informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries
. AERO 632: of Advance Flight Control System. Preliminaries Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Preliminaries Signals & Systems Laplace
More informationDIGITAL CONTROL OF POWER CONVERTERS. 2 Digital controller design
DIGITAL CONTROL OF POWER CONVERTERS 2 Digital controller design Outline Review of frequency domain control design Performance limitations Discrete time system analysis and modeling Digital controller design
More informationLinear dynamical systems with inputs & outputs
EE263 Autumn 215 S. Boyd and S. Lall Linear dynamical systems with inputs & outputs inputs & outputs: interpretations transfer function impulse and step responses examples 1 Inputs & outputs recall continuous-time
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 informationRELAY CONTROL WITH PARALLEL COMPENSATOR FOR NONMINIMUM PHASE PLANTS. Ryszard Gessing
RELAY CONTROL WITH PARALLEL COMPENSATOR FOR NONMINIMUM PHASE PLANTS Ryszard Gessing Politechnika Śl aska Instytut Automatyki, ul. Akademicka 16, 44-101 Gliwice, Poland, fax: +4832 372127, email: gessing@ia.gliwice.edu.pl
More informationDepartment of Electronics and Instrumentation Engineering M. E- CONTROL AND INSTRUMENTATION ENGINEERING CL7101 CONTROL SYSTEM DESIGN Unit I- BASICS AND ROOT-LOCUS DESIGN PART-A (2 marks) 1. What are the
More informationOptimal Polynomial Control for Discrete-Time Systems
1 Optimal Polynomial Control for Discrete-Time Systems Prof Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning this paper should
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 information1. 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 informationDEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2010
[E2.5] IMPERIAL COLLEGE LONDON DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2010 EEE/ISE PART II MEng. BEng and ACGI SIGNALS AND LINEAR SYSTEMS Time allowed: 2:00 hours There are FOUR
More information4F3 - Predictive Control
4F3 Predictive Control - Discrete-time systems p. 1/30 4F3 - Predictive Control Discrete-time State Space Control Theory For reference only Jan Maciejowski jmm@eng.cam.ac.uk 4F3 Predictive Control - Discrete-time
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 21: Stability Margins and Closing the Loop Overview In this Lecture, you will learn: Closing the Loop Effect on Bode Plot Effect
More informationDistributed Real-Time Control Systems
Distributed Real-Time Control Systems Chapter 9 Discrete PID Control 1 Computer Control 2 Approximation of Continuous Time Controllers Design Strategy: Design a continuous time controller C c (s) and then
More information1 x(k +1)=(Φ LH) x(k) = T 1 x 2 (k) x1 (0) 1 T x 2(0) T x 1 (0) x 2 (0) x(1) = x(2) = x(3) =
567 This is often referred to as Þnite settling time or deadbeat design because the dynamics will settle in a Þnite number of sample periods. This estimator always drives the error to zero in time 2T or
More informationDESIGN OF CMOS ANALOG INTEGRATED CIRCUITS
DESIGN OF CMOS ANALOG INEGRAED CIRCUIS Franco Maloberti Integrated Microsistems Laboratory University of Pavia Discrete ime Signal Processing F. Maloberti: Design of CMOS Analog Integrated Circuits Discrete
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 informationLecture 11. Frequency Response in Discrete Time Control Systems
EE42 - Discrete Time Systems Spring 28 Lecturer: Asst. Prof. M. Mert Ankarali Lecture.. Frequency Response in Discrete Time Control Systems Let s assume u[k], y[k], and G(z) represents the input, output,
More informationProfessor Fearing EE C128 / ME C134 Problem Set 7 Solution Fall 2010 Jansen Sheng and Wenjie Chen, UC Berkeley
Professor Fearing EE C8 / ME C34 Problem Set 7 Solution Fall Jansen Sheng and Wenjie Chen, UC Berkeley. 35 pts Lag compensation. For open loop plant Gs ss+5s+8 a Find compensator gain Ds k such that the
More information1. Cancellation of Left Half Plane Poles with Zeros
1. Cancellation of Left Half Plane Poles with Zeros 1. Cancellation of Left Half Plane Poles with Zeros Control of DC Motor (Astrom & Wittenmark) Velocity of motor shaft - x 1, position - x 2. 1. Cancellation
More information# FIR. [ ] = b k. # [ ]x[ n " k] [ ] = h k. x[ n] = Ae j" e j# ˆ n Complex exponential input. [ ]Ae j" e j ˆ. ˆ )Ae j# e j ˆ. y n. y n.
[ ] = h k M [ ] = b k x[ n " k] FIR k= M [ ]x[ n " k] convolution k= x[ n] = Ae j" e j ˆ n Complex exponential input [ ] = h k M % k= [ ]Ae j" e j ˆ % M = ' h[ k]e " j ˆ & k= k = H (" ˆ )Ae j e j ˆ ( )
More informationIngegneria dell Automazione - Sistemi in Tempo Reale p.1/28
Ingegneria dell Automazione - Sistemi in Tempo Reale Selected topics on discrete-time and sampled-data systems Luigi Palopoli palopoli@sssup.it - Tel. 050/883444 Ingegneria dell Automazione - Sistemi in
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 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 informationESC794: Special Topics: Model Predictive Control
ESC794: Special Topics: Model Predictive Control Discrete-Time Systems Hanz Richter, Professor Mechanical Engineering Department Cleveland State University Discrete-Time vs. Sampled-Data Systems A continuous-time
More informationTheory of Linear Systems Exercises. Luigi Palopoli and Daniele Fontanelli
Theory of Linear Systems Exercises Luigi Palopoli and Daniele Fontanelli Dipartimento di Ingegneria e Scienza dell Informazione Università di Trento Contents Chapter. Exercises on the Laplace Transform
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 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 informationState Feedback and State Estimators Linear System Theory and Design, Chapter 8.
1 Linear System Theory and Design, http://zitompul.wordpress.com 2 0 1 4 2 Homework 7: State Estimators (a) For the same system as discussed in previous slides, design another closed-loop state estimator,
More informationAN EXTENSION OF GENERALIZED BILINEAR TRANSFORMATION FOR DIGITAL REDESIGN. Received October 2010; revised March 2011
International Journal of Innovative Computing, Information and Control ICIC International c 2012 ISSN 1349-4198 Volume 8, Number 6, June 2012 pp. 4071 4081 AN EXTENSION OF GENERALIZED BILINEAR TRANSFORMATION
More informationAndrea Zanchettin Automatic Control 1 AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Discrete time linear systems
Andrea Zanchettin Automatic Control 1 AUTOMATIC CONTROL Andrea M. Zanchettin, PhD Spring Semester, 2018 Discrete time linear systems Andrea Zanchettin Automatic Control 2 Discrete time linear systems Modern
More informationSTABILITY ANALYSIS TECHNIQUES
ECE4540/5540: Digital Control Systems 4 1 STABILITY ANALYSIS TECHNIQUES 41: Bilinear transformation Three main aspects to control-system design: 1 Stability, 2 Steady-state response, 3 Transient response
More informationFrom Continuous-Time Domain to Microcontroller Code
APPLICAION NOE UnitedSiC_AN0019 October 018 From Continuous-ime Domain to Microcontroller Code By Jonathan Dodge, P.E. Introduction Control theory is one of the many aspects of electronic theory required
More informationCITY UNIVERSITY LONDON. MSc in Information Engineering DIGITAL SIGNAL PROCESSING EPM746
No: CITY UNIVERSITY LONDON MSc in Information Engineering DIGITAL SIGNAL PROCESSING EPM746 Date: 19 May 2004 Time: 09:00-11:00 Attempt Three out of FIVE questions, at least One question from PART B PART
More informationLABORATORY INSTRUCTION MANUAL CONTROL SYSTEM II LAB EE 693
LABORATORY INSTRUCTION MANUAL CONTROL SYSTEM II LAB EE 693 ELECTRICAL ENGINEERING DEPARTMENT JIS COLLEGE OF ENGINEERING (AN AUTONOMOUS INSTITUTE) KALYANI, NADIA EXPERIMENT NO : CS II/ TITLE : FAMILIARIZATION
More informationDESIGN USING TRANSFORMATION TECHNIQUE CLASSICAL METHOD
206 Spring Semester ELEC733 Digital Control System LECTURE 7: DESIGN USING TRANSFORMATION TECHNIQUE CLASSICAL METHOD For a unit ramp input Tz Ez ( ) 2 ( z ) D( z) G( z) Tz e( ) lim( z) z 2 ( z ) D( z)
More informationExercise 5: Digital Control
Gioele Zardini Control Systems II FS 017 Exercise 5: Digital Control 5.1 Digital Control 5.1.1 Signals and Systems A whole course is dedicated to this topic (see Signals and Systems of professor D Andrea).
More informationE2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2)
E.5 Signals & Linear Systems Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & ) 1. Sketch each of the following continuous-time signals, specify if the signal is periodic/non-periodic,
More informationV. IIR Digital Filters
Digital Signal Processing 5 March 5, V. IIR Digital Filters (Deleted in 7 Syllabus). (dded in 7 Syllabus). 7 Syllabus: nalog filter approximations Butterworth and Chebyshev, Design of IIR digital filters
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 2.161 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. I Reading:
More informationEE451/551: Digital Control. Chapter 3: Modeling of Digital Control Systems
EE451/551: Digital Control Chapter 3: Modeling of Digital Control Systems Common Digital Control Configurations AsnotedinCh1 commondigitalcontrolconfigurations As noted in Ch 1, common digital control
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 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 informationIntroduction to Digital Control. Week Date Lecture Title
http://elec3004.com Introduction to Digital Control 2016 School of Information Technology and Electrical Engineering at The University of Queensland Lecture Schedule: Week Date Lecture Title 1 29-Feb Introduction
More informationLike bilateral Laplace transforms, ROC must be used to determine a unique inverse z-transform.
Inversion of the z-transform Focus on rational z-transform of z 1. Apply partial fraction expansion. Like bilateral Laplace transforms, ROC must be used to determine a unique inverse z-transform. Let X(z)
More informationChapter 3 Data Acquisition and Manipulation
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
More informationLecture Discrete dynamic systems
Chapter 3 Low-level io Lecture 3.4 Discrete dynamic systems Lecture 3.4 Discrete dynamic systems Suppose that we wish to implement an embedded computer system that behaves analogously to a continuous linear
More informationExam. 135 minutes, 15 minutes reading time
Exam August 6, 208 Control Systems II (5-0590-00) Dr. Jacopo Tani Exam Exam Duration: 35 minutes, 5 minutes reading time Number of Problems: 35 Number of Points: 47 Permitted aids: 0 pages (5 sheets) A4.
More informationAnalysis of Discrete-Time Systems
TU Berlin Discrete-Time Control Systems 1 Analysis of Discrete-Time Systems Overview Stability Sensitivity and Robustness Controllability, Reachability, Observability, and Detectabiliy TU Berlin Discrete-Time
More informationEE480.3 Digital Control Systems. Part 7. Controller Design I. - Pole Assignment Method - State Estimation
EE480.3 Digital Control Systems Part 7. Controller Design I. - Pole Assignment Method - State Estimation Kunio Takaya Electrical and Computer Engineering University of Saskatchewan February 10, 2010 **
More informationEE 3CL4: Introduction to Control Systems Lab 4: Lead Compensation
EE 3CL4: Introduction to Control Systems Lab 4: Lead Compensation Tim Davidson Ext. 27352 davidson@mcmaster.ca Objective To use the root locus technique to design a lead compensator for a marginally-stable
More informationTime Response of Systems
Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =
More informationELEG 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 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 informationINTRODUCTION TO DIGITAL CONTROL
ECE4540/5540: Digital Control Systems INTRODUCTION TO DIGITAL CONTROL.: Introduction In ECE450/ECE550 Feedback Control Systems, welearnedhow to make an analog controller D(s) to control a linear-time-invariant
More informationExercise 1 (A Non-minimum Phase System)
Prof. Dr. E. Frazzoli 5-59- Control Systems I (HS 25) Solution Exercise Set Loop Shaping Noele Norris, 9th December 26 Exercise (A Non-minimum Phase System) To increase the rise time of the system, we
More informationDigital Control & Digital Filters. Lectures 1 & 2
Digital Controls & Digital Filters Lectures 1 & 2, Professor Department of Electrical and Computer Engineering Colorado State University Spring 2017 Digital versus Analog Control Systems Block diagrams
More information9.5 The Transfer Function
Lecture Notes on Control Systems/D. Ghose/2012 0 9.5 The Transfer Function Consider the n-th order linear, time-invariant dynamical system. dy a 0 y + a 1 dt + a d 2 y 2 dt + + a d n y 2 n dt b du 0u +
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 informationControl Systems I Lecture 10: System Specifications
Control Systems I Lecture 10: System Specifications Readings: Guzzella, Chapter 10 Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 24, 2017 E. Frazzoli (ETH) Lecture
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 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 informationCDS 101/110a: Lecture 8-1 Frequency Domain Design
CDS 11/11a: Lecture 8-1 Frequency Domain Design Richard M. Murray 17 November 28 Goals: Describe canonical control design problem and standard performance measures Show how to use loop shaping to achieve
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 informationRobust and Optimal Control, Spring A: SISO Feedback Control A.1 Internal Stability and Youla Parameterization
Robust and Optimal Control, Spring 2015 Instructor: Prof. Masayuki Fujita (S5-303B) A: SISO Feedback Control A.1 Internal Stability and Youla Parameterization A.2 Sensitivity and Feedback Performance A.3
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 18: State Feedback Tracking and State Estimation Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 18:
More informationA FEEDBACK STRUCTURE WITH HIGHER ORDER DERIVATIVES IN REGULATOR. Ryszard Gessing
A FEEDBACK STRUCTURE WITH HIGHER ORDER DERIVATIVES IN REGULATOR Ryszard Gessing Politechnika Śl aska Instytut Automatyki, ul. Akademicka 16, 44-101 Gliwice, Poland, fax: +4832 372127, email: gessing@ia.gliwice.edu.pl
More informationExercise 1 (A Non-minimum Phase System)
Prof. Dr. E. Frazzoli 5-59- Control Systems I (Autumn 27) Solution Exercise Set 2 Loop Shaping clruch@ethz.ch, 8th December 27 Exercise (A Non-minimum Phase System) To decrease the rise time of the system,
More informationLecture 5 - Assembly Programming(II), Intro to Digital Filters
GoBack Lecture 5 - Assembly Programming(II), Intro to Digital Filters James Barnes (James.Barnes@colostate.edu) Spring 2009 Colorado State University Dept of Electrical and Computer Engineering ECE423
More informationLecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.
ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition
More informationLABORATORY OF AUTOMATION SYSTEMS Analytical design of digital controllers
LABORATORY OF AUTOMATION SYSTEMS Analytical design of digital controllers Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI) Università di Bologna email: claudio.melchiorri@unibo.it
More informationChapter 2 z-transform. X(z) = Z[x(t)] = Z[x(kT)] = Z[x(k)] x(kt)z k = x(kt)z k = x(k)z k. X(z)z k 1 dz 2πj c = 1
One-sided -transform Two-sided -transform Chapter 2 -TRANSFORM X() Z[x(t)] Z[x(kT)] Z[x(k)] x(kt) k x(k) k X() k0 k x(kt) k k0 k x(k) k Note that X() x(0) + x(t ) + x(2t ) 2 + + x(kt) k + Inverse -transform
More information6.1 Sketch the z-domain root locus and find the critical gain for the following systems K., the closed-loop characteristic equation is K + z 0.
6. Sketch the z-domain root locus and find the critical gain for the following systems K (i) Gz () z 4. (ii) Gz K () ( z+ 9. )( z 9. ) (iii) Gz () Kz ( z. )( z ) (iv) Gz () Kz ( + 9. ) ( z. )( z 8. ) (i)
More informationControl Systems I. Lecture 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 7: Feedback and the Root Locus method Readings: Jacopo Tani Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 2, 2018 J. Tani, E. Frazzoli (ETH) Lecture 7:
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 informationSystems and Control Theory Lecture Notes. Laura Giarré
Systems and Control Theory Lecture Notes Laura Giarré L. Giarré 2017-2018 Lesson 7: Response of LTI systems in the transform domain Laplace Transform Transform-domain response (CT) Transfer function Zeta
More informationINTRODUCTION TO TRANSFER FUNCTIONS
INTRODUCTION TO TRANSFER FUNCTIONS The transfer function is the ratio of the output Laplace Transform to the input Laplace Transform assuming zero initial conditions. Many important characteristics of
More informationLTI system response. Daniele Carnevale. Dipartimento di Ing. Civile ed Ing. Informatica (DICII), University of Rome Tor Vergata
LTI system response Daniele Carnevale Dipartimento di Ing. Civile ed Ing. Informatica (DICII), University of Rome Tor Vergata Fondamenti di Automatica e Controlli Automatici A.A. 2014-2015 1 / 15 Laplace
More information