Modern Optimal Control
|
|
- Allen McLaughlin
- 5 years ago
- Views:
Transcription
1 Modern Optimal Control Matthew M. Peet Arizona State University Lecture 19: Stabilization via LMIs
2 Optimization Optimization can be posed in functional form: min x F objective function : inequality constraints subject to which may have the form min f 0(x) : x F f i (x) 0 subject to i = 1, k Special Cases: Linear Programming If the fi(x) : R n R m are affine (f i(x) = Ax b) EASY: Simplex/Ellipsoid Algorithm Polynomial Programming The fi(x) : R n R m are polynomial. Semidefinite Programming The fi(x) : R n S m are affine. For semidefinite programming, what does f i (x) 0 mean? M. Peet Lecture 19: 2 / 30
3 Positive Matrices Definition 1. A symmetric matrix P S n is Positive Semidefinite, denoted P 0 if x T P x 0 for all x R n Definition 2. A symmetric matrix P S n is Positive Definite, denoted P > 0 if x T P x > 0 for all x 0 P is Negative Semidefinite if P 0 P is Negative Definite if P > 0 A matrix which is neither Positive nor Negative Semidefinite is Indefinite The set of positive or negative matrices is a convex cone. M. Peet Lecture 19: 3 / 30
4 Positive Matrices Lemma 3. P S n is positive definite if and only if all its eigenvalues are positive. Easy Proofs: A Positive Definite matrix is invertible. The inverse of a positive definite matrix is positive definite. If P > 0, then T P T T 0 for any T. If T is invertible, then T P T T > 0. Lemma 4. For any P > 0, there exists a positive square root, P 1 2 P = P 1 2 P 1 2. > 0 such that M. Peet Lecture 19: 4 / 30
5 Convex Optimization Why is Linear Programming easy and polynomial programming hard? min f 0(x) : x F f i (x) 0 subject to i = 1, k Consider the geometric representation of the optimization problem: min f 0(x) : subject to x F x S where S := {x : f i (x) 0, i = 1,, k}. or the more compact pure geometric representation min γ,x F γ : subject to (γ, x) S where S := {(γ, x) : γ f 0 (x) 0, f i (x) 0, i = 1,, k}. M. Peet Lecture 19: 5 / 30
6 Convexity Definition 5. A set is convex if for any x, y Q, {µx + (1 µ)y : µ [0, 1]} Q. The line connecting any two points lies in the set. M. Peet Lecture 19: 6 / 30
7 Convexity Example Convex Optimization: minimize x 1 + x 2 Definitionsubject 6. to 3x 1 + x 2 3 Consider the optimization problem x 2 1 x min γ : subject 1 4 to γ,x F x 1 + 5x 2 20 (γ, x) x S x 2 20 subject to Ax = b Cx d The problem is Convex Optimization if the set S is convex. Convex optimization problems have the property that scaled Newton iteration always converges to the global optimal. The question is, of course, when is the set S convex? For polynomial optimization, a sufficient condition is that all functions f i are convex. The set of convex functions is a convex cone. The level set of a convex function is a convex set. M. Peet Lecture 19: 7 / 30
8 Convexity Newton s Algorithm: Solve f(x) = 0 x k+1 = x k t f(x k) f (x k ) t is the step-size. For non-convex optimization, Newton descent may get stuck at local optima. For Newton descent, constraints are represented by barrier functions. M. Peet Lecture 19: 8 / 30
9 Semidefinite Programming minimize trace CX subject to trace A i X = b i for all i X 0 The variable X is a symmetric matrix X 0 means X is positive semidefinite The feasible set is the intersection of an affine set with the positive semidefinite cone { X S n X 0 } M. Peet Lecture 19: 9 / 30
10 SDPs with Explicit Variables We can also explicitly parameterize the affine set to give minimize c T x subject to F 0 + x 1 F 1 + x 2 F x n F n 0 where F 0, F 1,..., F n are symmetric matrices. The inequality constraint is called a linear matrix inequality (LMI); e.g., x 1 3 x 1 + x 2 1 x 1 + x 2 x x 1 which is equivalent to x x M. Peet Lecture 19: 10 / 30
11 Linear Matrix Inequalities Review Linear Matrix Inequalities are often a Simpler way to solve control problems. Semidefinite Programming min c T x : Y 0 + i x i Y i > 0 Here the Y i are symmetric matrices. Linear Matrix Inequalities Findx : Y 0 + i x i Y i > 0 Commonly Takes the Form FindX : A i XB i + Q > 0 i M. Peet Lecture 19: 11 / 30
12 Linear Matrix Inequalities Common Form: FindX : A i XB i + Q > 0 i The most important Linear Matrix Inequality is the Lyapunov Inequality. There are several very efficient LMI/SDP Solvers for Matlab: SeDuMi Fast, but somewhat unreliable. See LMI Lab (Part of Matlab s Robust Control Toolbox) Universally disliked See YALMIP (a parser for other solvers) See I recommend YALMIP with solver SeDuMi. M. Peet Lecture 19: 12 / 30
13 Semidefinite Programming(SDP): Common Examples in Control Some Simple examples of LMI conditions in control include: Stability Stabilization H 2 Synthesis We will go beyond these examples. A T X + XA 0 X 0 AX + BZ + XA T + Z T B T 0 X 0 min T r(w ) [ ] [ ] X A B2 + [ X Z ] [ ] Z T A T B2 T + B 1 B1 T 0 [ ] X (CX + DZ) T 0 (CX + DZ) W M. Peet Lecture 19: 13 / 30
14 Lyapunov Theory LMIs unite time-domain and frequency-domain analysis Theorem 7 (Lyapunov). Suppose there exists a continuously differentiable function V for which V (0) = 0 and V (x) > 0 for x 0. Furthermore, suppose lim x V (x) = and V (x(t + h)) V (x(t)) lim = d h 0 + h dt V (x(t)) < 0 for any x such that ẋ(t) = f(x(t)). Then for any x(0) R the system of equations ẋ(t) = f(x(t)) has a unique solution which is stable in the sense of Lyapunov. M. Peet Lecture 19: 14 / 30
15 The Lyapunov Inequality Lemma 8. A is Hurwitz if and only if there exists a P > 0 such that A T P + P A < 0 Proof. Suppose there exists a P > 0 such that A T P + P A < 0. Define the Lyapunov function V (x) = x T P x. Then V (x) > 0 for x 0 and V (0) = 0. Furthermore, V (x(t)) = ẋ(t) T P x(t) + x(t) T P ẋ(t) = x(t) T A T P x(t) + x(t) T P Ax(t) = x(t) T ( A T P + P A ) x(t) Hence V (x(t)) < 0 for all x 0. Thus the system is globally stable. Global stability implies A is Hurwitz. M. Peet Lecture 19: 15 / 30
16 The Lyapunov Inequality Proof. For the other direction, if A is Hurwitz, let P = Converges because A is Hurwitz. Furthermore P A = e AT s e As Ads = = e AT s Ae As ds = [ e AT s e As ] 0 = I Thus P A + A T P = I < 0. 0 e AT s e As ds e AT s d ( e As ) ds 0 ds d T 0 ds ea s e As A T e AT s e As = I A T P M. Peet Lecture 19: 16 / 30
17 The Lyapunov Inequality Other Versions: Lemma 9. (A, B) is controllable if and only if there exists a X > 0 such that A T X + XA + BB T 0 Lemma 10. (C, A) is stabilizable if and only if there exists a X > 0 such that AX + XA T + C T C 0 These also yield Lyapunov functions for the system. M. Peet Lecture 19: 17 / 30
18 The Static State-Feedback Problem Lets start with the problem of stabilization. Definition 11. The Static State-Feedback Problem is to find a feedback matrix K such that ẋ(t) = Ax(t) + Bu(t) u(t) = Kx(t) is stable We have already solved the problem earlier using the Controllability Grammian. Find K such that A + BK is Hurwitz. Can also be put in LMI format: Find X > 0, K : X(A + BK) + (A + BK) T X < 0 Problem: Bilinear in K and X. M. Peet Lecture 19: 18 / 30
19 The Static State-Feedback Problem The bilinear problem in K and X is a common paradigm. Bilinear optimization is not convex. To convexify the problem, we use a change of variables. Problem 1: Problem 2: Definition 12. Find X > 0, K : X(A + BK) + (A + BK) T X < 0 Find P > 0, Z : AP + BZ + P A T + ZB T < 0 Two optimization problems are equivalent if a solution to one will provide a solution to the other. Theorem 13. Problem 1 is equivalent to Problem 2. M. Peet Lecture 19: 19 / 30
20 The Dual Lyapunov Equation Problem 1: Find X > 0, : Lemma 14. XA + A T X < 0 Problem 1 is equivalent to problem 2. Problem 2: Find Y > 0, : Y A T + AY < 0 Proof. First we show 1) solves 2). Suppose X > 0 is a solution to Problem 1. Let Y = X 1 > 0. If XA + A T X < 0, then Hence X 1 (XA + A T X)X 1 < 0 X 1 (XA + A T X)X 1 = AX 1 + X 1 A T = AY + Y A T < 0 Therefore, Problem 2 is feasible with solution Y = X 1. M. Peet Lecture 19: 20 / 30
21 The Dual Lyapunov Equation Problem 1: Find X > 0, : XA + A T X < 0 Problem 2: Find Y > 0, : Y A T + AY < 0 Proof. Now we show 2) solves 1) in a similar manner. Suppose Y > 0 is a solution to Problem 1. Let X = Y 1 > 0. Then XA + A T X = X(AX 1 + X 1 A T )X = X(AY + Y A T )X < 0 Conclusion: If V (x) = x T P x proves stability of ẋ = Ax, Then V (x) = x T P 1 x proves stability of ẋ = A T x. M. Peet Lecture 19: 21 / 30
22 The Stabilization Problem Thus we rephrase Problem 1 Problem 1: Find P > 0, K : (A + BK)P + P (A + BK) T < 0 Theorem 15. Problem 1 is equivalent to Problem 2. Problem 2: Find X > 0, Z : AX + BZ + XA T + ZB T < 0 Proof. We will show that 2) Solves 1). Suppose X > 0, Z solves 2). Let P = X > 0 and K = ZP 1. Then Z = KP and (A + BK)P + P (A + BK) T = AP + P A T + BKP + P K T B T = AP + P A T + BZ + Z T B T < 0 Now suppose that P > 0 and K solve 1). Let X = P > 0 and Z = KP. Then AP + P A T + BZ + Z T B T = (A + BK)P + P (A + BK) T < 0 M. Peet Lecture 19: 22 / 30
23 The Stabilization Problem The result can be summarized more succinctly Theorem 16. (A, B) is static-state-feedback stabilizable if and only if there exists some P > 0 and Z such that AP + P A T + BZ + Z T B T < 0 with u(t) = ZP 1 x(t). Standard Format: [ ] [ ] P A B + [ P Z ] [ ] Z T A T B T < 0 There are a number of general-purpose LMI solvers available. e.g. SeDuMi - Free LMI Lab - in Matlab Robust Control Toolbox (in Computer lab) YALMIP - A nice front end for several solvers (Free) M. Peet Lecture 19: 23 / 30
24 The Schur complement Before we get to the main result, recall the Schur complement. Theorem 17 (Schur Complement). For any S S n, Q S m and R R n m, the following are equivalent. [ ] M R 1. R T > 0 Q 2. Q > 0 and M RQ 1 R T > 0 A commonly used property of positive matrices. Also Recall: If X > 0, then X ɛi > 0 for ɛ sufficiently small. M. Peet Lecture 19: 24 / 30
25 The KYP Lemma (AKA: The Bounded Real Lemma) The most important theorem in this class. Lemma 18. Suppose Then the following are equivalent. G H γ. [ A B Ĝ(s) = C D There exists a X > 0 such that [ ] A T X + XA XB B T + 1 [ ] C T [C ] D < 0 X γi γ D T ]. Can be used to calculate the H -norm of a system Originally used to solve LMI s using graphs. (Before Computers) Now used directly instead of graphical methods like Bode. The feasibility constraints are linear Can be combined with other methods. M. Peet Lecture 19: 25 / 30
26 The KYP Lemma Proof. We will only show that ii) implies i). The other direction requires the Hamiltonian, which we have not discussed. We will show that if y = Gu, then y L2 γ u L2. From the 1 x 1 block of the LMI, we know that A T X + XA < 0, which means A is Hurwitz. Because the inequality is strict, there exists some ɛ > 0 such that [ ] A T X + XA XB B T + 1 [ ] C T [C ] D X (γ ɛ)i γ D T [ ] A = T X + XA XB B T + 1 [ ] [ ] C T [C ] 0 0 D + X γi γ D T < 0 0 ɛi Let y = Gu. Then the state-space representation is y(t) = Cx(t) + Du(t) ẋ(t) = Ax(t) + Bu(t) x(0) = 0 M. Peet Lecture 19: 26 / 30
27 The KYP Lemma Proof. Let V (x) = x T Xx. Then the LMI implies [ ] [ T [A ] x(t) T X + XA XB u(t) B T + 1 [ ] C T [C ] ] [ ] x(t) D X (γ ɛ)i γ D T u(t) [ ] T [ ] [ ] x A = T X + XA XB x u B T + 1 [ ] T [ ] x C T [C ] [ ] x D X (γ ɛ)i u γ u D T u [ ] T [ [ ] x A = T X + XA XB x u B T + X (γ ɛ)i] 1 u γ yt y = x T (A T X + XA)x + x T XBu + u T B T Xx (γ ɛ)u T u + 1 γ yt y = (Ax + Bu) T Xx + x T X(Ax + Bu) (γ ɛ)u T u + 1 γ yt y = ẋ(t) T Xx(t) + x(t) T Xẋ(t) (γ ɛ) u(t) γ y(t) 2 = V (x(t)) (γ ɛ) u(t) γ y(t) 2 < 0 M. Peet Lecture 19: 27 / 30
28 The KYP Lemma Proof. Now we have V (x(t)) (γ ɛ) u(t) γ y(t) 2 < 0 Integrating in time, we get T ( V (x(t)) (γ ɛ) u(t) γ y(t) 2) dt 0 = V (x(t )) V (x(0)) (γ ɛ) T 0 u(t) 2 dt + 1 γ Because A is Hurwitz, lim t x(t) = 0. Hence lim t V (x(t)) = 0. Likewise, because x(0) = 0, we have V (x(0)) = 0. T 0 y(t) 2) dt < 0 M. Peet Lecture 19: 28 / 30
29 The KYP Lemma Proof. Since V (x(0)) = V (x( )) = 0, [ lim V (x(t )) V (x(0)) (γ ɛ) T = 0 0 (γ ɛ) 0 u(t) 2 dt + 1 γ = (γ ɛ) u 2 L γ y 2 L 2 dt < 0 T 0 0 u(t) 2 dt + 1 γ y(t) 2 dt T 0 y(t) 2) ] dt Thus y 2 L 2 dt < (γ 2 ɛγ) u 2 L 2 By definition, this means G 2 H (γ 2 ɛγ) < γ 2 or G H < γ M. Peet Lecture 19: 29 / 30
30 The Positive Real Lemma A Passivity Condition A Variation on the KYP lemma is the positive-real lemma Lemma 19. Suppose Then the following are equivalent. [ A B Ĝ(s) = C D G is passive. i.e. ( u, Gu L2 0). There exists a P > 0 such that [ ] A T P + P A P B C T B T P C D T 0 D ]. M. Peet Lecture 19: 30 / 30
LMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 4: LMIs for State-Space Internal Stability Solving the Equations Find the output given the input State-Space:
More informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 22: H 2, LQG and LGR Conclusion To solve the H -optimal state-feedback problem, we solve min γ such that γ,x 1,Y 1,A n,b n,c n,d
More informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 21: Optimal Output Feedback Control connection is called the (lower) star-product of P and Optimal Output Feedback ansformation (LFT).
More informationGrammians. Matthew M. Peet. Lecture 20: Grammians. Illinois Institute of Technology
Grammians Matthew M. Peet Illinois Institute of Technology Lecture 2: Grammians Lyapunov Equations Proposition 1. Suppose A is Hurwitz and Q is a square matrix. Then X = e AT s Qe As ds is the unique solution
More informationLMIs for Observability and Observer Design
LMIs for Observability and Observer Design Matthew M. Peet Arizona State University Lecture 06: LMIs for Observability and Observer Design Observability Consider a system with no input: ẋ(t) = Ax(t), x(0)
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 informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 14: LMIs for Robust Control in the LF Framework ypes of Uncertainty In this Lecture, we will cover Unstructured,
More informationDenis ARZELIER arzelier
COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.2 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS PERFORMANCE ANALYSIS and SYNTHESIS Denis ARZELIER www.laas.fr/ arzelier arzelier@laas.fr 15
More informationThe norms can also be characterized in terms of Riccati inequalities.
9 Analysis of stability and H norms Consider the causal, linear, time-invariant system ẋ(t = Ax(t + Bu(t y(t = Cx(t Denote the transfer function G(s := C (si A 1 B. Theorem 85 The following statements
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 02: Optimization (Convex and Otherwise) What is Optimization? An Optimization Problem has 3 parts. x F f(x) :
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationAnalysis of Systems with State-Dependent Delay
Analysis of Systems with State-Dependent Delay Matthew M. Peet Arizona State University Tempe, AZ USA American Institute of Aeronautics and Astronautics Guidance, Navigation and Control Conference Boston,
More informationLecture 10: Linear Matrix Inequalities Dr.-Ing. Sudchai Boonto
Dr-Ing Sudchai Boonto Department of Control System and Instrumentation Engineering King Mongkuts Unniversity of Technology Thonburi Thailand Linear Matrix Inequalities A linear matrix inequality (LMI)
More informationFast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma
Fast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 8 September 2003 European Union RTN Summer School on Multi-Agent
More informationLinear Matrix Inequalities in Robust Control. Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University MTNS 2002
Linear Matrix Inequalities in Robust Control Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University MTNS 2002 Objective A brief introduction to LMI techniques for Robust Control Emphasis on
More informationSemidefinite Programming Basics and Applications
Semidefinite Programming Basics and Applications Ray Pörn, principal lecturer Åbo Akademi University Novia University of Applied Sciences Content What is semidefinite programming (SDP)? How to represent
More informationEE363 homework 8 solutions
EE363 Prof. S. Boyd EE363 homework 8 solutions 1. Lyapunov condition for passivity. The system described by ẋ = f(x, u), y = g(x), x() =, with u(t), y(t) R m, is said to be passive if t u(τ) T y(τ) dτ
More informationROBUST ANALYSIS WITH LINEAR MATRIX INEQUALITIES AND POLYNOMIAL MATRICES. Didier HENRION henrion
GRADUATE COURSE ON POLYNOMIAL METHODS FOR ROBUST CONTROL PART IV.1 ROBUST ANALYSIS WITH LINEAR MATRIX INEQUALITIES AND POLYNOMIAL MATRICES Didier HENRION www.laas.fr/ henrion henrion@laas.fr Airbus assembly
More informationModern Control Systems
Modern Control Systems Matthew M. Peet Arizona State University Lecture 09: Observability Observability For Static Full-State Feedback, we assume knowledge of the Full-State. In reality, we only have measurements
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 informationLecture Note 5: Semidefinite Programming for Stability Analysis
ECE7850: Hybrid Systems:Theory and Applications Lecture Note 5: Semidefinite Programming for Stability Analysis Wei Zhang Assistant Professor Department of Electrical and Computer Engineering Ohio State
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 informationLinear Matrix Inequalities in Control
Linear Matrix Inequalities in Control Delft Center for Systems and Control (DCSC) Delft University of Technology The Netherlands Department of Electrical Engineering Eindhoven University of Technology
More informationME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ
ME 234, Lyapunov and Riccati Problems. This problem is to recall some facts and formulae you already know. (a) Let A and B be matrices of appropriate dimension. Show that (A, B) is controllable if and
More informationLinear Matrix Inequality (LMI)
Linear Matrix Inequality (LMI) A linear matrix inequality is an expression of the form where F (x) F 0 + x 1 F 1 + + x m F m > 0 (1) x = (x 1,, x m ) R m, F 0,, F m are real symmetric matrices, and the
More informationNonlinear Systems Theory
Nonlinear Systems Theory Matthew M. Peet Arizona State University Lecture 2: Nonlinear Systems Theory Overview Our next goal is to extend LMI s and optimization to nonlinear systems analysis. Today we
More informationIntroduction to linear matrix inequalities Wojciech Paszke
Introduction to linear matrix inequalities Wojciech Paszke Institute of Control and Computation Engineering, University of Zielona Góra, Poland e-mail: W.Paszke@issi.uz.zgora.pl Outline Introduction to
More informationRiccati Equations and Inequalities in Robust Control
Riccati Equations and Inequalities in Robust Control Lianhao Yin Gabriel Ingesson Martin Karlsson Optimal Control LP4 2014 June 10, 2014 Lianhao Yin Gabriel Ingesson Martin Karlsson (LTH) H control problem
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 informationOptimization of Polynomials
Optimization of Polynomials Matthew M. Peet Arizona State University Thanks to S. Lall and P. Parrilo for guidance and supporting material Lecture 03: Optimization of Polynomials Overview In this lecture,
More informationSYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1. P. V. Pakshin, S. G. Soloviev
SYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1 P. V. Pakshin, S. G. Soloviev Nizhny Novgorod State Technical University at Arzamas, 19, Kalinina ul., Arzamas, 607227,
More informationRobust and Optimal Control, Spring 2015
Robust and Optimal Control, Spring 2015 Instructor: Prof. Masayuki Fujita (S5-303B) D. Linear Matrix Inequality D.1 Convex Optimization D.2 Linear Matrix Inequality(LMI) D.3 Control Design and LMI Formulation
More information9 The LQR Problem Revisited
9 he LQR Problem Revisited Problem: Compute a state feedback controller u(t) = Kx(t) that stabilizes the closed loop system and minimizes J := lim E [ z(t) z(t) ] t for the LI system Assumptions: a) D
More informationH 2 Suboptimal Estimation and Control for Nonnegative
Proceedings of the 2007 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 2007 FrC20.3 H 2 Suboptimal Estimation and Control for Nonnegative Dynamical Systems
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 informationMarcus Pantoja da Silva 1 and Celso Pascoli Bottura 2. Abstract: Nonlinear systems with time-varying uncertainties
A NEW PROPOSAL FOR H NORM CHARACTERIZATION AND THE OPTIMAL H CONTROL OF NONLINEAR SSTEMS WITH TIME-VARING UNCERTAINTIES WITH KNOWN NORM BOUND AND EXOGENOUS DISTURBANCES Marcus Pantoja da Silva 1 and Celso
More informationVerification and Control of Safety-Factor Profile in Tokamaks
Verification and Control of Safety-Factor Profile in Tokamaks Aditya Gahlawat, Matthew M. Peet, and Emmanuel Witrant Illinois Institute of Technology Chicago, Illinois International Federation of Automatic
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May
More informationRobust Stability. Robust stability against time-invariant and time-varying uncertainties. Parameter dependent Lyapunov functions
Robust Stability Robust stability against time-invariant and time-varying uncertainties Parameter dependent Lyapunov functions Semi-infinite LMI problems From nominal to robust performance 1/24 Time-Invariant
More informationLecture: Convex Optimization Problems
1/36 Lecture: Convex Optimization Problems http://bicmr.pku.edu.cn/~wenzw/opt-2015-fall.html Acknowledgement: this slides is based on Prof. Lieven Vandenberghe s lecture notes Introduction 2/36 optimization
More informationLinear Matrix Inequalities in Control
C. W. Scherer a S. Weiland b Linear Matrix Inequalities in Control a Mathematical Systems Theory, University of Stuttgart, Pfaffenwaldring 57, 7569 Stuttgart/ Germany Carsten.scherer@mathematik.uni-stuttgart.de
More informationRobust Multivariable Control
Lecture 2 Anders Helmersson anders.helmersson@liu.se ISY/Reglerteknik Linköpings universitet Today s topics Today s topics Norms Today s topics Norms Representation of dynamic systems Today s topics Norms
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 informationModern Control Systems
Modern Control Systems Matthew M. Peet Arizona State University Lecture 1: Introduction MMAE 543: Modern Control Systems This course is on Modern Control Systems Techniques Developed in the Last 50 years
More informationLecture 6: Conic Optimization September 8
IE 598: Big Data Optimization Fall 2016 Lecture 6: Conic Optimization September 8 Lecturer: Niao He Scriber: Juan Xu Overview In this lecture, we finish up our previous discussion on optimality conditions
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 informationLecture 6. Foundations of LMIs in System and Control Theory
Lecture 6. Foundations of LMIs in System and Control Theory Ivan Papusha CDS270 2: Mathematical Methods in Control and System Engineering May 4, 2015 1 / 22 Logistics hw5 due this Wed, May 6 do an easy
More informationRobust Anti-Windup Compensation for PID Controllers
Robust Anti-Windup Compensation for PID Controllers ADDISON RIOS-BOLIVAR Universidad de Los Andes Av. Tulio Febres, Mérida 511 VENEZUELA FRANCKLIN RIVAS-ECHEVERRIA Universidad de Los Andes Av. Tulio Febres,
More informationConvex Optimization. (EE227A: UC Berkeley) Lecture 6. Suvrit Sra. (Conic optimization) 07 Feb, 2013
Convex Optimization (EE227A: UC Berkeley) Lecture 6 (Conic optimization) 07 Feb, 2013 Suvrit Sra Organizational Info Quiz coming up on 19th Feb. Project teams by 19th Feb Good if you can mix your research
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 informationLinear-quadratic control problem with a linear term on semiinfinite interval: theory and applications
Linear-quadratic control problem with a linear term on semiinfinite interval: theory and applications L. Faybusovich T. Mouktonglang Department of Mathematics, University of Notre Dame, Notre Dame, IN
More informationOutline. Linear Matrix Inequalities in Control. Outline. System Interconnection. j _jst. ]Bt Bjj. Generalized plant framework
Outline Linear Matrix Inequalities in Control Carsten Scherer and Siep Weiland 7th Elgersburg School on Mathematical Systems heory Class 3 1 Single-Objective Synthesis Setup State-Feedback Output-Feedback
More information4. Convex optimization problems
Convex Optimization Boyd & Vandenberghe 4. Convex optimization problems optimization problem in standard form convex optimization problems quasiconvex optimization linear optimization quadratic optimization
More informationFrom Convex Optimization to Linear Matrix Inequalities
Dep. of Information Engineering University of Pisa (Italy) From Convex Optimization to Linear Matrix Inequalities eng. Sergio Grammatico grammatico.sergio@gmail.com Class of Identification of Uncertain
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 5. Input-Output Stability DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Input-Output Stability y = Hu H denotes
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 informationModern Control Systems
Modern Control Systems Matthew M. Peet Illinois Institute of Technology Lecture 1: Modern Control Systems MMAE 543: Modern Control Systems This course is on Modern Control Systems Techniques Developed
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 informationLinear Algebra. P R E R E Q U I S I T E S A S S E S S M E N T Ahmad F. Taha August 24, 2015
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 5243 INTRODUCTION TO CYBER-PHYSICAL SYSTEMS P R E R E Q U I S I T E S A S S E S S M E N T Ahmad F. Taha August 24, 2015 The objective of this exercise is to assess
More informationCourse Outline. FRTN10 Multivariable Control, Lecture 13. General idea for Lectures Lecture 13 Outline. Example 1 (Doyle Stein, 1979)
Course Outline FRTN Multivariable Control, Lecture Automatic Control LTH, 6 L-L Specifications, models and loop-shaping by hand L6-L8 Limitations on achievable performance L9-L Controller optimization:
More informationIntroduction to Nonlinear Control Lecture # 4 Passivity
p. 1/6 Introduction to Nonlinear Control Lecture # 4 Passivity È p. 2/6 Memoryless Functions ¹ y È Ý Ù È È È È u (b) µ power inflow = uy Resistor is passive if uy 0 p. 3/6 y y y u u u (a) (b) (c) Passive
More informationMulti-objective Controller Design:
Multi-objective Controller Design: Evolutionary algorithms and Bilinear Matrix Inequalities for a passive suspension A. Molina-Cristobal*, C. Papageorgiou**, G. T. Parks*, M. C. Smith**, P. J. Clarkson*
More informationIntroduction to Linear Matrix Inequalities (LMIs)
ECE 680 Fall 2017 Introduction to Linear Matrix Inequalities (LMIs) by Stanislaw H. Żak and Guisheng Zhai October 18, 2017 Henry Ford offered customers the Model T Ford in any color the customer wants,
More informationOptimality Conditions for Constrained Optimization
72 CHAPTER 7 Optimality Conditions for Constrained Optimization 1. First Order Conditions In this section we consider first order optimality conditions for the constrained problem P : minimize f 0 (x)
More informationInterior Point Methods: Second-Order Cone Programming and Semidefinite Programming
School of Mathematics T H E U N I V E R S I T Y O H F E D I N B U R G Interior Point Methods: Second-Order Cone Programming and Semidefinite Programming Jacek Gondzio Email: J.Gondzio@ed.ac.uk URL: http://www.maths.ed.ac.uk/~gondzio
More informationConvex Optimization 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Convex Optimization 1 Many optimization objectives generated
More informationLecture 6 Verification of Hybrid Systems
Lecture 6 Verification of Hybrid Systems Ufuk Topcu Nok Wongpiromsarn Richard M. Murray AFRL, 25 April 2012 Outline: A hybrid system model Finite-state abstractions and use of model checking Deductive
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
More informationContinuous Optimisation, Chpt 9: Semidefinite Problems
Continuous Optimisation, Chpt 9: Semidefinite Problems Peter J.C. Dickinson DMMP, University of Twente p.j.c.dickinson@utwente.nl http://dickinson.website/teaching/2016co.html version: 21/11/16 Monday
More informationConstrained Nonlinear Optimization Algorithms
Department of Industrial Engineering and Management Sciences Northwestern University waechter@iems.northwestern.edu Institute for Mathematics and its Applications University of Minnesota August 4, 2016
More informationControllability, Observability, Full State Feedback, Observer Based Control
Multivariable Control Lecture 4 Controllability, Observability, Full State Feedback, Observer Based Control John T. Wen September 13, 24 Ref: 3.2-3.4 of Text Controllability ẋ = Ax + Bu; x() = x. At time
More informationEE 227A: Convex Optimization and Applications October 14, 2008
EE 227A: Convex Optimization and Applications October 14, 2008 Lecture 13: SDP Duality Lecturer: Laurent El Ghaoui Reading assignment: Chapter 5 of BV. 13.1 Direct approach 13.1.1 Primal problem Consider
More informationECE 680 Linear Matrix Inequalities
ECE 680 Linear Matrix Inequalities Stan Żak School of Electrical and Computer Engineering Purdue University zak@purdue.edu October 11, 2017 Stan Żak (Purdue University) ECE 680 Linear Matrix Inequalities
More informationFilter Design for Feedback-loop Trade-off of L 1 Adaptive Controller: A Linear Matrix Inequality Approach
AIAA Guidance, Navigation and Control Conference and Exhibit 18-21 August 2008, Honolulu, Hawaii AIAA 2008-6280 Filter Design for Feedback-loop Trade-off of L 1 Adaptive Controller: A Linear Matrix Inequality
More informationIII. Applications in convex optimization
III. Applications in convex optimization nonsymmetric interior-point methods partial separability and decomposition partial separability first order methods interior-point methods Conic linear optimization
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 informationAn Introduction to Linear Matrix Inequalities. Raktim Bhattacharya Aerospace Engineering, Texas A&M University
An Introduction to Linear Matrix Inequalities Raktim Bhattacharya Aerospace Engineering, Texas A&M University Linear Matrix Inequalities What are they? Inequalities involving matrix variables Matrix variables
More informationRobust and Optimal Control, Spring 2015
Robust and Optimal Control, Spring 2015 Instructor: Prof. Masayuki Fujita (S5-303B) G. Sum of Squares (SOS) G.1 SOS Program: SOS/PSD and SDP G.2 Duality, valid ineqalities and Cone G.3 Feasibility/Optimization
More informationLinear Matrix Inequalities for Normalizing Matrices
Linear Matrix Inequalities for Normalizing Matrices Christian Ebenbauer Abstract A real square matrix is normal if it can be diagonalized by an unitary matrix. In this paper novel convex conditions are
More informationEE363 homework 7 solutions
EE363 Prof. S. Boyd EE363 homework 7 solutions 1. Gain margin for a linear quadratic regulator. Let K be the optimal state feedback gain for the LQR problem with system ẋ = Ax + Bu, state cost matrix Q,
More information6-1 The Positivstellensatz P. Parrilo and S. Lall, ECC
6-1 The Positivstellensatz P. Parrilo and S. Lall, ECC 2003 2003.09.02.10 6. The Positivstellensatz Basic semialgebraic sets Semialgebraic sets Tarski-Seidenberg and quantifier elimination Feasibility
More informationThe Q-parametrization (Youla) Lecture 13: Synthesis by Convex Optimization. Lecture 13: Synthesis by Convex Optimization. Example: Spring-mass System
The Q-parametrization (Youla) Lecture 3: Synthesis by Convex Optimization controlled variables z Plant distubances w Example: Spring-mass system measurements y Controller control inputs u Idea for lecture
More informationWritten Examination
Division of Scientific Computing Department of Information Technology Uppsala University Optimization Written Examination 202-2-20 Time: 4:00-9:00 Allowed Tools: Pocket Calculator, one A4 paper with notes
More informationTheorem 1. ẋ = Ax is globally exponentially stable (GES) iff A is Hurwitz (i.e., max(re(σ(a))) < 0).
Linear Systems Notes Lecture Proposition. A M n (R) is positive definite iff all nested minors are greater than or equal to zero. n Proof. ( ): Positive definite iff λ i >. Let det(a) = λj and H = {x D
More informationLecture 17: Primal-dual interior-point methods part II
10-725/36-725: Convex Optimization Spring 2015 Lecture 17: Primal-dual interior-point methods part II Lecturer: Javier Pena Scribes: Pinchao Zhang, Wei Ma Note: LaTeX template courtesy of UC Berkeley EECS
More informationResearch Article An Equivalent LMI Representation of Bounded Real Lemma for Continuous-Time Systems
Hindawi Publishing Corporation Journal of Inequalities and Applications Volume 28, Article ID 67295, 8 pages doi:1.1155/28/67295 Research Article An Equivalent LMI Representation of Bounded Real Lemma
More informationEE363 homework 2 solutions
EE363 Prof. S. Boyd EE363 homework 2 solutions. Derivative of matrix inverse. Suppose that X : R R n n, and that X(t is invertible. Show that ( d d dt X(t = X(t dt X(t X(t. Hint: differentiate X(tX(t =
More informationMultiobjective Optimization Applied to Robust H 2 /H State-feedback Control Synthesis
Multiobjective Optimization Applied to Robust H 2 /H State-feedback Control Synthesis Eduardo N. Gonçalves, Reinaldo M. Palhares, and Ricardo H. C. Takahashi Abstract This paper presents an algorithm for
More informationLecture 7: September 17
10-725: Optimization Fall 2013 Lecture 7: September 17 Lecturer: Ryan Tibshirani Scribes: Serim Park,Yiming Gu 7.1 Recap. The drawbacks of Gradient Methods are: (1) requires f is differentiable; (2) relatively
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 informationRichard Anthony Conway. A dissertation submitted in partial satisfaction of the requirements for the degree of Doctor of Philosophy
Discrete-Time H 2 Guaranteed Cost Control by Richard Anthony Conway A dissertation submitted in partial satisfaction of the requirements for the degree of Doctor of Philosophy in Engineering Mechanical
More informationRobust Multi-Objective Control for Linear Systems
Robust Multi-Objective Control for Linear Systems Elements of theory and ROMULOC toolbox Dimitri PEAUCELLE & Denis ARZELIER LAAS-CNRS, Toulouse, FRANCE Part of the OLOCEP project (includes GloptiPoly)
More informationLinearization problem. The simplest example
Linear Systems Lecture 3 1 problem Consider a non-linear time-invariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and
More informationECEN 605 LINEAR SYSTEMS. Lecture 7 Solution of State Equations 1/77
1/77 ECEN 605 LINEAR SYSTEMS Lecture 7 Solution of State Equations Solution of State Space Equations Recall from the previous Lecture note, for a system: ẋ(t) = A x(t) + B u(t) y(t) = C x(t) + D u(t),
More informationCS675: Convex and Combinatorial Optimization Fall 2016 Convex Optimization Problems. Instructor: Shaddin Dughmi
CS675: Convex and Combinatorial Optimization Fall 2016 Convex Optimization Problems Instructor: Shaddin Dughmi Outline 1 Convex Optimization Basics 2 Common Classes 3 Interlude: Positive Semi-Definite
More informationThe moment-lp and moment-sos approaches
The moment-lp and moment-sos approaches LAAS-CNRS and Institute of Mathematics, Toulouse, France CIRM, November 2013 Semidefinite Programming Why polynomial optimization? LP- and SDP- CERTIFICATES of POSITIVITY
More informationACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016
ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe
More informationControl Systems. LMIs in. Guang-Ren Duan. Analysis, Design and Applications. Hai-Hua Yu. CRC Press. Taylor & Francis Croup
LMIs in Control Systems Analysis, Design and Applications Guang-Ren Duan Hai-Hua Yu CRC Press Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 20: LMI/SOS Tools for the Study of Hybrid Systems Stability Concepts There are several classes of problems for
More informationLearning Control Lyapunov Functions from Counterexamples and Demonstrations
Noname manuscript No. (will be inserted by the editor) Learning Control Lyapunov Functions from Counterexamples and Demonstrations Hadi Ravanbakhsh Sriram Sankaranarayanan Received: date / Accepted: date
More information