ESC794: Special Topics: Model Predictive Control
|
|
- Isabella Anthony
- 5 years ago
- Views:
Transcription
1 ESC794: Special Topics: Model Predictive Control Discrete-Time Systems Hanz Richter, Professor Mechanical Engineering Department Cleveland State University
2 Discrete-Time vs. Sampled-Data Systems A continuous-time (C-T) system operates on input functions. It evolves according to a set of differential equations where time t is the independent variable, t R. A discrete-time (D-T) system operates on input sequences. It evolves according to a set of difference equations where sequences are indexed with an integer k N 0. A sampled-data system is a D-T system obtained from a C-T system by a sampling process. Inputs, however, are functions defined on a continuous sampling interval. See remarks on G&P, p.21. Sampled-data analysis will not be included in this course. We will describe the dynamics of D-T systems using the difference equation x + = f(x,u) where the transition map f : X U X returns the next state x + given a current state and control value. By repeatedly applying a sequence u(k), we obtain a state sequence x(k). Exercise: find the transition map of a savings account with an interest rate of i and initial balance x 0. Positive u represent deposits and negative u represent withdrawals. 2 / 23
3 Zero-Order Hold Sampling: LTI Systems Clock u(k) u(t) Continuous y(t) y(k) D/A System A/D Sample and hold (ZOH) Sample and hold (ZOH) We may assume that there is no delay between the application of u k and the effect on y k (do not confuse lag with delay). The piecewise-constant signal u(t) is treated as input and y(t) is obtained from the general state-space solution of the continuous-time system Then y(t) is sampled (in sync with u(k)) to give y(k). u(k) maps into y(k) through a discrete-time dynamic system that can be obtained knowing the continuous-time model and the sampling period. 3 / 23
4 Zero-Order Hold Equivalent Calculations Let the LTI continuous-time system be represented by ẋ = Ax+Bu; y = Cx+Du The general solution given initial condition x 0 is well-known: x(t) = e A(t t 0) x(t 0 )+ t t 0 e A(t τ) Bu(τ)dτ A state-space discrete equivalent based on the ZOH is found by substituting a piecewise-constant function u(t k ) and integrating by segments: x(k +1) = Φx(k)+Γu(k); y(k) = Cx(k)+Du(k) where Φ = e AT s Γ = Ts 0 e As Bds 4 / 23
5 Some ways to compute Φ 1. Use series expansion e AT s = I +AT s +A 2 T 2 s/2+... and hope that A n becomes zero for some n (nilpotent matrix). 2. Use the Laplace identity (only adequate for low orders) e At = L 1 (si A) 1 3. Use the Cayley-Hamilton theorem (convenient for analytical evaluations) 4. Numerically in Matlab or symbolically in Maple/Mathematica/Mathcad. 5 / 23
6 Using the Cayley-Hamilton Theorem A square matrix satisfies its own characteristic equation. This can be clearly seen from the eigenvalue equation det(λi A) = 0. If we define a matrix function f(a) (for example e At ), the Cayley-Hamilton theorem can be used to show that f(a) = α 0 A n 1 +α 1 A n α n 1 I It can be shown that the eigenvalues of A must satisfy the above. If we have n distinct eigenvalues, we can determine the α i, and consequently the matrix exponential. If an eigenvalue is repeated m times, we take the first m 1 derivatives of the above to obtain additional conditions. 6 / 23
7 Refresher Exercise : Linear Systems Theory Obtain the ZOH equivalent of the double-integrator G(s) = 1/s 2 for T s = 1 using the four methods. 7 / 23
8 Difference equation representations The state-space discrete-time system presented above can be considered to be the realization of a linear n-th order input-output difference equation: y(k)+a 1 y(k 1)+...a na y(k n a ) = b 0 u(k d)+...+b nb u(k d n b ) Conversely, the difference equation can be obtained from the state-space system by algebraic elimination of states (entirely analogous to continuous-time case). To obtain a state-space realization from a difference equation, proceed as in the continuous-time case. For instance, states may be chosen as y(k) and its delayed versions. To recall state-space realization methods from I/O equations see Dorf, Modern Control Systems. Exercise: for the CT transfer function G(s) = s+1 s 2 +s+1 obtain a ZOH discretization in state-space and find the corresponding difference equation. From it, find a state-space realization that matches the one returned by Matlab s c2d applied to G(s). 8 / 23
9 Exercise - Simulink Use a Simulink model to simulate the C-T double integrator under ZOH and the D-T equivalent in state-space and difference equation forms. Apply some input and check that the sampled-data system and the D-T systems all give the same output sequences. 9 / 23
10 Discretization of nonlinear C-T systems How do we obtain x + = f(x,u) from ẋ = h(x,u) when h is a nonlinear function? We will use two methods in this course: 1. Euler discretization (forward differences) 2. Numerical integration of the C-T system (using ode45 or similar solvers) There are other options, including implicit forms that use a backward difference, and in some cases, exact analytical discretization can be attempted. The numerical integration method is useful to generate online MPC predictions when dealing with C-T system models. We simply apply a constant input over a short period of time (corresponding to the sampling period δ) and use a numerical solver to find the corresponding state trajectories. For this, we write a Matlab function x=dynamics(x0,delta,u) which calls ode45 using an initial point and a constant input to determine the state x + at the end of the simulation. 10 / 23
11 Forward Euler Discretization Given a C-T system in the form we approximate the time derivative with ẋ = h(x,u) ẋ x(t+δ) x(t) δ Substitution gives and using t = kδ gives: x(t+δ) = x(t)+δh(x(t),u(t)) x(k +1) = x(k)+δh(x(k),u(k)) or simply with f(x,u) = x+δh(x,u). x + = f(x,u) The above gives an explicit discretization. If we used the backward difference, we wouldn t be able to solve for x + in the general nonlinear case. In the linear case this is possible. 11 / 23
12 Discretization of LTI and Stability Preservation Suppose ẋ = Ax and a forward Euler discretization is used. Then the discrete equivalent is x + = (I +δa)x The C-T is asymptotically stable if and only if all eigenvalues λ C of A lie in the region defined by Re(λ) < 0. The D-T is asymptotically stable if and only if all eigenvalues λ D of I +δa lie in the unit circle, defined by λ D < 1. What is the relationship between these two regions for a given δ? Does forward Euler preserve stability? Now suppose that a backward Euler discretization is used. What is the discrete equivalent? What can we say about stability preservation for LTI systems? Fact: A+I is invertible is none of the eigenvalues of A is equal to / 23
13 Basic State-Space D-T LTI Design The pole placement and LQR methods follow the same principles as in the C-T case, using the open unit circle as stability region. For low-order designs, pole placement can be useful: 1. Use the 2nd-order C-T response as a prototype: select a damping ratio ζ and natural frequency w n for the dominant poles. 2. Use a pole mapping from C-T to D-T to determine the approximate locations of the desired poles in D-T. Keep in mind that zeroes can appear in the D-T system even when the C-T did not have them. For details see Franklin, Powell & Workman or Åstrom & Wittenmark. 3. The ZOH pole mapping is z = e st s. 13 / 23
14 Pole placement Since the closed-loop system matrix is still A BK, we may use the same formulas and tools of C-T pole placement. 2. For state feedback with setpoint following, check that the (SISO) control law is u = Pr Kx, where the reference gain is P = 1 (C DK)(I (A BK)) 1 B +D 3. Deadbeat response: If we can place all closed-loop poles at zero, the setpoint will be reached in a finite number of steps, which is at most equal to the system order. See FPW. Exercise: Design K and P for the double-integrator plant discretized with ZOH, Ts = 0.1. The settling time should be 1 second and the overshoot 5% or less. Redesign for deadbeat response. Check in Simulink against the D-T and sampled-data systems. 14 / 23
15 D-T LQR As in C-T, an optimal design can be attempted using an infinite-horizon quadratic performance index: J = x T (k)qx(k)+u T (k)ru(k) Assume the following: k=0 1. (A, B) is stabilizable (uncontrollable modes are asymptotically stable) 2. R > 0 and Q 0 3. (Q,A) has no unobservable modes on the unit circle Then the solution is a state feedback controller u = Kx where K = (B T PB +R) 1 B T PA and P is the unique positive-definite solution of the discrete algebraic Riccati equation (DARE): A T PA P A T PB(R+B T PB) 1 B T PA+Q = 0 The Matlab function dlqr calculates K and P. 15 / 23
16 Quadratic Lyapunov Functions and Invariance MPC was developed to seek the best performance for the given constraints. Constraints are usually given as a set of polyhedral regions in the state and control spaces, for example 10 x 1 10, 5 x 2 5, or more generally Fx 1, where F has as many rows as constraint boundaries. When a LTI system is controlled with linear state feedback to give an asymptotically stable closed-loop, ellipsoidal invariant sets appear, which are easy to compute. Invariant ellipsoids inscribed in a polyhedral constraint set are conservative safe operating regions. More importantly, ellipsoids are used as terminal regional constraints to establish stability of MPC, including nonlinear cases. An ellipsoidal set centered at the origin of R n is defined as for some P > 0. E = {x R n x T Px 1} 16 / 23
17 Finding invariant ellipsoids: DLTI Suppose the closed-loop system is x + = A cl x where A cl is asymptotically stable. Consider the quadratic Lyapunov function with arbitrary P = P T > 0 V(x) = x T Px The difference V = V(x + ) V(x) along closed-loop trajectories is V = x T (A T clpa cl P)x Intuitively at this point, we know that V(k) must be monotonically decreasing with k, since x(k) is exponentially converging to zero. This means that V must be negative for any value of x (because any such value could be used as initial condition). Then A T cl PA cl P must be negative definite. We can generate P by selecting Q > 0 and solving the discrete Lyapunov equation A T clpa cl P = Q where Q sets how fast V decreases. In Matlab dlyap can be used, adjusting for a transposed A cl. 17 / 23
18 Finding invariant ellipsoids... The level sets of V define ellipsoidal sets. Therefore, the monotonic decrease of V indicates that E must be positively-invariant, that is: x(0) E = x(k) E for k > 0 We will use the Ellipsoidal Toolbox (Alex Kurzhanskiy, UC Berkeley) for various operations with ellipsoids and hyperplanes. Example: We find a conservative set (E) of initial conditions such that the states of the controlled D-T double-integrator (for r = 0) satisfy unit box constraints. The largest E should be found. Then we plot some trajectories. We minimize the trace of P to obtain the largest volume for E We require P = P T > 0 The Lyapunov inequality must hold We must keep E inside the unit box This can be efficiently programmed using LMIs, but we will use the general-purpose optimizer fmincon and Lagrange multiplier techniques. 18 / 23
19 Example... We first solve the useful sub-problem: given a family of hyperplanes defined by γ = Γx, minimize γ subject to x T Px = 1, where P = P T > 0. Check that the solution is γ = ΓP 1 Γ T and interpret graphically in 2D using the Ellipsoidal Toolbox. The box constraints are defined with Γ 1 = [1 0], Γ 2 = [ 1 0], Γ 3 = [0 1] and Γ 4 = [0 1]. To simplify the optimization problem, we search over Q = Q T > 0, since any such choice will automatically guarantee P = P T > 0 and the Lyapunov inequality. 19 / 23
20 Example... The optimization problem is then formulated as minimize Q trace P subject to Q = Q T > 0 A T clp +PA cl Q = 0 γ 1 1 < 0 γ 3 1 < 0 This is readily solved using fmincon s equality constraints and nonlinear constraints capabilities. Note that the Lyapunov equation is handled by substitution (not programmed as an equality constraint). The algorithm tries a Q, then finds P and checks for constraint satisfaction. The requirement Q = Q T is programmed as an equality constraint based on a vectorized form of Q (using reshape). 20 / 23
21 Example Invariant Ellipsoid and Trajectories: Double Integrator x x 1 The ellipsoidal region is conservative. Initial conditions outside E may still be feasible, we don t know until we do a simulation. 21 / 23
22 Output Admissible Sets The initial state of an unforced linear system is output admissible with respect to a constraint set Y if the resulting output function satisfies the pointwise-in-time condition y(t) Y, t 0. The set of all possible such initial conditions is the maximal output admissible set O. Gilbert E.G. and Tan, K.T., Linear Systems with State and Control Constraints: The Theory and Application of Maximal Output Admissible Sets, IEEE Trans. Automatic Control, 36(9) pp , Gilbert and Tan showed that O preserved the closure, symmetry and convexity of Y and proposed an algorithm to approximate O with a finite intersection of regions (several ellipsoids in the figure below). 22 / 23
23 Semi-Ellipsoidal Sets and Invariant Cylinders For linear systems under state feedback, O Dell proposed that a good tradeoff between the complexity of polyhedrals and the conservatism of ellipsoids could be achieved by intersecting an invariant ellipsoid with the constraints themselves to define a semi-ellipsoidal set. This set was designed by maximizing the ellipsoid volume under applicable constraints. O Dell, B.D. and Misawa, E.A, Semi-Ellipsoidal Controlled Invariant Sets for Constrained Linear Systems, ASME J. Dynamic Systems, Measurements and Control, 124(1) pp , For systems under sliding mode control, Richter followed a similar approach, using infinite cylinders as primary invariant sets. Richter, H., O Dell, B.D. and Misawa, E.A, Robust Positively Invariant Cylinders in Constrained Variable Structure Control, IEEE Trans. Automatic Control, 52 (11) pp , / 23
ESC794: Special Topics: Model Predictive Control
ESC794: Special Topics: Model Predictive Control Nonlinear MPC Analysis : Part 1 Reference: Nonlinear Model Predictive Control (Ch.3), Grüne and Pannek Hanz Richter, Professor Mechanical Engineering Department
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 informationEE C128 / ME C134 Feedback Control Systems
EE C128 / ME C134 Feedback Control Systems Lecture Additional Material Introduction to Model Predictive Control Maximilian Balandat Department of Electrical Engineering & Computer Science University of
More informationIMPROVED MPC DESIGN BASED ON SATURATING CONTROL LAWS
IMPROVED MPC DESIGN BASED ON SATURATING CONTROL LAWS D. Limon, J.M. Gomes da Silva Jr., T. Alamo and E.F. Camacho Dpto. de Ingenieria de Sistemas y Automática. Universidad de Sevilla Camino de los Descubrimientos
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 informationLecture 2 and 3: Controllability of DT-LTI systems
1 Lecture 2 and 3: Controllability of DT-LTI systems Spring 2013 - EE 194, Advanced Control (Prof Khan) January 23 (Wed) and 28 (Mon), 2013 I LTI SYSTEMS Recall that continuous-time LTI systems can be
More information5. Observer-based Controller Design
EE635 - Control System Theory 5. Observer-based Controller Design Jitkomut Songsiri state feedback pole-placement design regulation and tracking state observer feedback observer design LQR and LQG 5-1
More informationMathematics for Control Theory
Mathematics for Control Theory Outline of Dissipativity and Passivity Hanz Richter Mechanical Engineering Department Cleveland State University Reading materials Only as a reference: Charles A. Desoer
More informationControl Systems Lab - SC4070 Control techniques
Control Systems Lab - SC4070 Control techniques Dr. Manuel Mazo Jr. Delft Center for Systems and Control (TU Delft) m.mazo@tudelft.nl Tel.:015-2788131 TU Delft, February 16, 2015 (slides modified from
More informationAutonomous Mobile Robot Design
Autonomous Mobile Robot Design Topic: Guidance and Control Introduction and PID Loops Dr. Kostas Alexis (CSE) Autonomous Robot Challenges How do I control where to go? Autonomous Mobile Robot Design Topic:
More informationIEOR 265 Lecture 14 (Robust) Linear Tube MPC
IEOR 265 Lecture 14 (Robust) Linear Tube MPC 1 LTI System with Uncertainty Suppose we have an LTI system in discrete time with disturbance: x n+1 = Ax n + Bu n + d n, where d n W for a bounded polytope
More informationẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)
EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and
More informationModule 07 Controllability and Controller Design of Dynamical LTI Systems
Module 07 Controllability and Controller Design of Dynamical LTI Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha October
More informationStability lectures. Stability of Linear Systems. Stability of Linear Systems. Stability of Continuous Systems. EECE 571M/491M, Spring 2008 Lecture 5
EECE 571M/491M, Spring 2008 Lecture 5 Stability of Continuous Systems http://courses.ece.ubc.ca/491m moishi@ece.ubc.ca Dr. Meeko Oishi Electrical and Computer Engineering University of British Columbia,
More informationPrashant Mhaskar, Nael H. El-Farra & Panagiotis D. Christofides. Department of Chemical Engineering University of California, Los Angeles
HYBRID PREDICTIVE OUTPUT FEEDBACK STABILIZATION OF CONSTRAINED LINEAR SYSTEMS Prashant Mhaskar, Nael H. El-Farra & Panagiotis D. Christofides Department of Chemical Engineering University of California,
More informationOutline. 1 Linear Quadratic Problem. 2 Constraints. 3 Dynamic Programming Solution. 4 The Infinite Horizon LQ Problem.
Model Predictive Control Short Course Regulation James B. Rawlings Michael J. Risbeck Nishith R. Patel Department of Chemical and Biological Engineering Copyright c 217 by James B. Rawlings Outline 1 Linear
More informationState Regulator. Advanced Control. design of controllers using pole placement and LQ design rules
Advanced Control State Regulator Scope design of controllers using pole placement and LQ design rules Keywords pole placement, optimal control, LQ regulator, weighting matrixes Prerequisites Contact state
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 14: Controllability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 14: Controllability p.1/23 Outline
More informationEE221A Linear System Theory Final Exam
EE221A Linear System Theory Final Exam Professor C. Tomlin Department of Electrical Engineering and Computer Sciences, UC Berkeley Fall 2016 12/16/16, 8-11am Your answers must be supported by analysis,
More informationLecture 2: Discrete-time Linear Quadratic Optimal Control
ME 33, U Berkeley, Spring 04 Xu hen Lecture : Discrete-time Linear Quadratic Optimal ontrol Big picture Example onvergence of finite-time LQ solutions Big picture previously: dynamic programming and finite-horizon
More informationMATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem
MATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem Pulemotov, September 12, 2012 Unit Outline Goal 1: Outline linear
More informationConstrained Biogeography-Based Optimization for Invariant Set Computation
Constrained Biogeography-Based Optimization for Invariant Set Computation Arpit Shah, Dan Simon and Hanz Richter Abstract We discuss the application of biogeography-based optimization (BBO) to invariant
More information4F3 - Predictive Control
4F3 Predictive Control - Lecture 2 p 1/23 4F3 - Predictive Control Lecture 2 - Unconstrained Predictive Control Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 2 p 2/23 References Predictive
More informationTopic # Feedback Control Systems
Topic #17 16.31 Feedback Control Systems Deterministic LQR Optimal control and the Riccati equation Weight Selection Fall 2007 16.31 17 1 Linear Quadratic Regulator (LQR) Have seen the solutions to the
More informationHomework Solution # 3
ECSE 644 Optimal Control Feb, 4 Due: Feb 17, 4 (Tuesday) Homework Solution # 3 1 (5%) Consider the discrete nonlinear control system in Homework # For the optimal control and trajectory that you have found
More information7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system
7 Stability 7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system ẋ(t) = A x(t), x(0) = x 0, A R n n, x 0 R n. (14) The origin x = 0 is a globally asymptotically
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Introduction to Nonlinear Controllability and Observability Hanz Richter Mechanical Engineering Department Cleveland State University Definition of
More informationECE7850 Lecture 7. Discrete Time Optimal Control and Dynamic Programming
ECE7850 Lecture 7 Discrete Time Optimal Control and Dynamic Programming Discrete Time Optimal control Problems Short Introduction to Dynamic Programming Connection to Stabilization Problems 1 DT nonlinear
More informationControl Systems I. Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback. Readings: Emilio Frazzoli
Control Systems I Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 13, 2017 E. Frazzoli (ETH)
More informationSampling of Linear Systems
Sampling of Linear Systems Real-Time Systems, Lecture 6 Karl-Erik Årzén January 26, 217 Lund University, Department of Automatic Control Lecture 6: Sampling of Linear Systems [IFAC PB Ch. 1, Ch. 2, and
More informationTime-Invariant Linear Quadratic Regulators!
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 17 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
More informationContents. 1 State-Space Linear Systems 5. 2 Linearization Causality, Time Invariance, and Linearity 31
Contents Preamble xiii Linear Systems I Basic Concepts 1 I System Representation 3 1 State-Space Linear Systems 5 1.1 State-Space Linear Systems 5 1.2 Block Diagrams 7 1.3 Exercises 11 2 Linearization
More informationControl Systems. Design of State Feedback Control.
Control Systems Design of State Feedback Control chibum@seoultech.ac.kr Outline Design of State feedback control Dominant pole design Symmetric root locus (linear quadratic regulation) 2 Selection of closed-loop
More informationDecentralized and distributed control
Decentralized and distributed control Centralized control for constrained discrete-time systems M. Farina 1 G. Ferrari Trecate 2 1 Dipartimento di Elettronica, Informazione e Bioingegneria (DEIB) Politecnico
More informationTime-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2015
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 15 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
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 informationOPTIMAL CONTROL. Sadegh Bolouki. Lecture slides for ECE 515. University of Illinois, Urbana-Champaign. Fall S. Bolouki (UIUC) 1 / 28
OPTIMAL CONTROL Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 28 (Example from Optimal Control Theory, Kirk) Objective: To get from
More information16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1
16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1 Charles P. Coleman October 31, 2005 1 / 40 : Controllability Tests Observability Tests LEARNING OUTCOMES: Perform controllability tests Perform
More informationFEL3210 Multivariable Feedback Control
FEL3210 Multivariable Feedback Control Lecture 8: Youla parametrization, LMIs, Model Reduction and Summary [Ch. 11-12] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 8: Youla, LMIs, Model Reduction
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationMODERN CONTROL DESIGN
CHAPTER 8 MODERN CONTROL DESIGN The classical design techniques of Chapters 6 and 7 are based on the root-locus and frequency response that utilize only the plant output for feedback with a dynamic controller
More informationExtensions and applications of LQ
Extensions and applications of LQ 1 Discrete time systems 2 Assigning closed loop pole location 3 Frequency shaping LQ Regulator for Discrete Time Systems Consider the discrete time system: x(k + 1) =
More informationOptimal control and estimation
Automatic Control 2 Optimal control and estimation Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011
More informationCDS Solutions to Final Exam
CDS 22 - Solutions to Final Exam Instructor: Danielle C Tarraf Fall 27 Problem (a) We will compute the H 2 norm of G using state-space methods (see Section 26 in DFT) We begin by finding a minimal state-space
More information1 Steady State Error (30 pts)
Professor Fearing EECS C28/ME C34 Problem Set Fall 2 Steady State Error (3 pts) Given the following continuous time (CT) system ] ẋ = A x + B u = x + 2 7 ] u(t), y = ] x () a) Given error e(t) = r(t) y(t)
More informationCDS Solutions to the Midterm Exam
CDS 22 - Solutions to the Midterm Exam Instructor: Danielle C. Tarraf November 6, 27 Problem (a) Recall that the H norm of a transfer function is time-delay invariant. Hence: ( ) Ĝ(s) = s + a = sup /2
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems
CDS 101 1. For each of the following linear systems, determine whether the origin is asymptotically stable and, if so, plot the step response and frequency response for the system. If there are multiple
More informationLinear System Theory
Linear System Theory Wonhee Kim Chapter 6: Controllability & Observability Chapter 7: Minimal Realizations May 2, 217 1 / 31 Recap State space equation Linear Algebra Solutions of LTI and LTV system Stability
More informationJoão P. Hespanha. January 16, 2009
LINEAR SYSTEMS THEORY João P. Hespanha January 16, 2009 Disclaimer: This is a draft and probably contains a few typos. Comments and information about typos are welcome. Please contact the author at hespanha@ece.ucsb.edu.
More informationReview of Controllability Results of Dynamical System
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 13, Issue 4 Ver. II (Jul. Aug. 2017), PP 01-05 www.iosrjournals.org Review of Controllability Results of Dynamical System
More informationAppendix A Solving Linear Matrix Inequality (LMI) Problems
Appendix A Solving Linear Matrix Inequality (LMI) Problems In this section, we present a brief introduction about linear matrix inequalities which have been used extensively to solve the FDI problems described
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 informationAutomatic Control 2. Model reduction. Prof. Alberto Bemporad. University of Trento. Academic year
Lecture: Automatic Control 2 Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011 1 / 17 Lecture:
More informationHybrid Systems - Lecture n. 3 Lyapunov stability
OUTLINE Focus: stability of equilibrium point Hybrid Systems - Lecture n. 3 Lyapunov stability Maria Prandini DEI - Politecnico di Milano E-mail: prandini@elet.polimi.it continuous systems decribed by
More informationEEE582 Homework Problems
EEE582 Homework Problems HW. Write a state-space realization of the linearized model for the cruise control system around speeds v = 4 (Section.3, http://tsakalis.faculty.asu.edu/notes/models.pdf). Use
More informationModel Predictive Control Short Course Regulation
Model Predictive Control Short Course Regulation James B. Rawlings Michael J. Risbeck Nishith R. Patel Department of Chemical and Biological Engineering Copyright c 2017 by James B. Rawlings Milwaukee,
More informationSolution of Linear State-space Systems
Solution of Linear State-space Systems Homogeneous (u=0) LTV systems first Theorem (Peano-Baker series) The unique solution to x(t) = (t, )x 0 where The matrix function is given by is called the state
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationSubject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)
Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must
More informationCourse Outline. Higher Order Poles: Example. Higher Order Poles. Amme 3500 : System Dynamics & Control. State Space Design. 1 G(s) = s(s + 2)(s +10)
Amme 35 : System Dynamics Control State Space Design Course Outline Week Date Content Assignment Notes 1 1 Mar Introduction 2 8 Mar Frequency Domain Modelling 3 15 Mar Transient Performance and the s-plane
More informationModule 03 Linear Systems Theory: Necessary Background
Module 03 Linear Systems Theory: Necessary Background Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha/index.html September
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 informationSYSTEMTEORI - ÖVNING 5: FEEDBACK, POLE ASSIGNMENT AND OBSERVER
SYSTEMTEORI - ÖVNING 5: FEEDBACK, POLE ASSIGNMENT AND OBSERVER Exercise 54 Consider the system: ẍ aẋ bx u where u is the input and x the output signal (a): Determine a state space realization (b): Is the
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 informationFINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES. Danlei Chu, Tongwen Chen, Horacio J. Marquez
FINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES Danlei Chu Tongwen Chen Horacio J Marquez Department of Electrical and Computer Engineering University of Alberta Edmonton
More informationECE7850 Lecture 8. Nonlinear Model Predictive Control: Theoretical Aspects
ECE7850 Lecture 8 Nonlinear Model Predictive Control: Theoretical Aspects Model Predictive control (MPC) is a powerful control design method for constrained dynamical systems. The basic principles and
More informationProblem 2 (Gaussian Elimination, Fundamental Spaces, Least Squares, Minimum Norm) Consider the following linear algebraic system of equations:
EEE58 Exam, Fall 6 AA Rodriguez Rules: Closed notes/books, No calculators permitted, open minds GWC 35, 965-37 Problem (Dynamic Augmentation: State Space Representation) Consider a dynamical system consisting
More informationLinear-Quadratic Optimal Control: Full-State Feedback
Chapter 4 Linear-Quadratic Optimal Control: Full-State Feedback 1 Linear quadratic optimization is a basic method for designing controllers for linear (and often nonlinear) dynamical systems and is actually
More informationStabilization of constrained linear systems via smoothed truncated ellipsoids
Preprints of the 8th IFAC World Congress Milano (Italy) August 28 - September 2, 2 Stabilization of constrained linear systems via smoothed truncated ellipsoids A. Balestrino, E. Crisostomi, S. Grammatico,
More informationLQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin
LQR, Kalman Filter, and LQG Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin May 2015 Linear Quadratic Regulator (LQR) Consider a linear system
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 informationPOLE PLACEMENT. Sadegh Bolouki. Lecture slides for ECE 515. University of Illinois, Urbana-Champaign. Fall S. Bolouki (UIUC) 1 / 19
POLE PLACEMENT Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 19 Outline 1 State Feedback 2 Observer 3 Observer Feedback 4 Reduced Order
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 informationHybrid Systems Course Lyapunov stability
Hybrid Systems Course Lyapunov stability OUTLINE Focus: stability of an equilibrium point continuous systems decribed by ordinary differential equations (brief review) hybrid automata OUTLINE Focus: stability
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationLinear State Feedback Controller Design
Assignment For EE5101 - Linear Systems Sem I AY2010/2011 Linear State Feedback Controller Design Phang Swee King A0033585A Email: king@nus.edu.sg NGS/ECE Dept. Faculty of Engineering National University
More informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 19: Stabilization via LMIs Optimization Optimization can be posed in functional form: min x F objective function : inequality constraints
More informationTheory in Model Predictive Control :" Constraint Satisfaction and Stability!
Theory in Model Predictive Control :" Constraint Satisfaction and Stability Colin Jones, Melanie Zeilinger Automatic Control Laboratory, EPFL Example: Cessna Citation Aircraft Linearized continuous-time
More informationEnlarged terminal sets guaranteeing stability of receding horizon control
Enlarged terminal sets guaranteeing stability of receding horizon control J.A. De Doná a, M.M. Seron a D.Q. Mayne b G.C. Goodwin a a School of Electrical Engineering and Computer Science, The University
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 informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More informationFall 線性系統 Linear Systems. Chapter 08 State Feedback & State Estimators (SISO) Feng-Li Lian. NTU-EE Sep07 Jan08
Fall 2007 線性系統 Linear Systems Chapter 08 State Feedback & State Estimators (SISO) Feng-Li Lian NTU-EE Sep07 Jan08 Materials used in these lecture notes are adopted from Linear System Theory & Design, 3rd.
More informationLinear Systems. Manfred Morari Melanie Zeilinger. Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich
Linear Systems Manfred Morari Melanie Zeilinger Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich Spring Semester 2016 Linear Systems M. Morari, M. Zeilinger - Spring
More informationTopic # Feedback Control
Topic #5 6.3 Feedback Control State-Space Systems Full-state Feedback Control How do we change the poles of the state-space system? Or,evenifwecanchangethepolelocations. Where do we put the poles? Linear
More informationAutomatic Control 2. Nonlinear systems. Prof. Alberto Bemporad. University of Trento. Academic year
Automatic Control 2 Nonlinear systems Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011 1 / 18
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 informationVideo 6.1 Vijay Kumar and Ani Hsieh
Video 6.1 Vijay Kumar and Ani Hsieh Robo3x-1.6 1 In General Disturbance Input + - Input Controller + + System Output Robo3x-1.6 2 Learning Objectives for this Week State Space Notation Modeling in the
More informationImplicit Euler Numerical Simulation of Sliding Mode Systems
Implicit Euler Numerical Simulation of Sliding Mode Systems Vincent Acary and Bernard Brogliato, INRIA Grenoble Rhône-Alpes Conference in honour of E. Hairer s 6th birthday Genève, 17-2 June 29 Objectives
More informationECE504: Lecture 8. D. Richard Brown III. Worcester Polytechnic Institute. 28-Oct-2008
ECE504: Lecture 8 D. Richard Brown III Worcester Polytechnic Institute 28-Oct-2008 Worcester Polytechnic Institute D. Richard Brown III 28-Oct-2008 1 / 30 Lecture 8 Major Topics ECE504: Lecture 8 We are
More informationNOTICE WARNING CONCERNING COPYRIGHT RESTRICTIONS: The copyright law of the United States (title 17, U.S. Code) governs the making of photocopies or
NOTICE WARNING CONCERNING COPYRIGHT RESTRICTIONS: The copyright law of the United States (title 17, U.S. Code) governs the making of photocopies or other reproductions of copyrighted material. Any copying
More informationQuadratic Stability of Dynamical Systems. Raktim Bhattacharya Aerospace Engineering, Texas A&M University
.. Quadratic Stability of Dynamical Systems Raktim Bhattacharya Aerospace Engineering, Texas A&M University Quadratic Lyapunov Functions Quadratic Stability Dynamical system is quadratically stable if
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear
More informationEE451/551: Digital Control. Chapter 8: Properties of State Space Models
EE451/551: Digital Control Chapter 8: Properties of State Space Models Equilibrium State Definition 8.1: An equilibrium point or state is an initial state from which the system nevers departs unless perturbed
More informationSemidefinite Programming Duality and Linear Time-invariant Systems
Semidefinite Programming Duality and Linear Time-invariant Systems Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 2 July 2004 Workshop on Linear Matrix Inequalities in Control LAAS-CNRS,
More informationEML5311 Lyapunov Stability & Robust Control Design
EML5311 Lyapunov Stability & Robust Control Design 1 Lyapunov Stability criterion In Robust control design of nonlinear uncertain systems, stability theory plays an important role in engineering systems.
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 informationarzelier
COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.1 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS STABILITY ANALYSIS Didier HENRION www.laas.fr/ henrion henrion@laas.fr Denis ARZELIER www.laas.fr/
More informationStability of Parameter Adaptation Algorithms. Big picture
ME5895, UConn, Fall 215 Prof. Xu Chen Big picture For ˆθ (k + 1) = ˆθ (k) + [correction term] we haven t talked about whether ˆθ(k) will converge to the true value θ if k. We haven t even talked about
More informationECE504: Lecture 9. D. Richard Brown III. Worcester Polytechnic Institute. 04-Nov-2008
ECE504: Lecture 9 D. Richard Brown III Worcester Polytechnic Institute 04-Nov-2008 Worcester Polytechnic Institute D. Richard Brown III 04-Nov-2008 1 / 38 Lecture 9 Major Topics ECE504: Lecture 9 We are
More informationLearning Model Predictive Control for Iterative Tasks: A Computationally Efficient Approach for Linear System
Learning Model Predictive Control for Iterative Tasks: A Computationally Efficient Approach for Linear System Ugo Rosolia Francesco Borrelli University of California at Berkeley, Berkeley, CA 94701, USA
More information