Computation of the Constrained Infinite Time Linear Quadratic Optimal Control Problem
|
|
- Dinah Harmon
- 5 years ago
- Views:
Transcription
1 Computation of the Constrained Infinite Time Linear Quadratic Optimal Control Problem July 5, Introduction Abstract Problem Statement and Properties In this paper we will consider discrete-time linear time-invariant systems x(t + ) = Ax(t) + Bu(t) (a) with n states and m inputs, i.e. A R n n and B R m n.. Finite-Time Unconstrained LQR (FTLQR) Consider the unconstrained finite-time optimal control problem N JN (x()) = min { U x(n) P x(n) + u(i) Ru(i) + x(i) Qx(i)} () for system () with Q, P and R. The solution to () yields the state-feedback optimal control law i= u (i) = K i x(i) i =,..., N (3) where the gain matrices K i are given by the equation K i = (B P i+ B + R) B P i+ A, (4) and where the symmetric positive semidefinite matrices P i are given recursively by the algorithm P N = Q (5) P i = A (P i+ P i+ B(B P i+ B + R) BP i+ )A + Q. (6) The optimal cost is given by J N (x()) = x() P x(). (7). Infinite-Time Unconstrained LQR (LQR) If in () we set N = + and we assume that (A, B) stabilizable and (A, C) detectable (Q = C C), then we obtain the standard infinite time linear quadratic regulator (LQR) problem J (x()) = min U { u(i) Ru(i) + x(i) Qx(i)} (8) i=
2 whose solution can be expressed as the state feedback control law where the gain matrix K is given by u (i) = Kx(i), i =,..., + (9) K = (B P B + R) B P A () and P is the unique solution of the algebraic matrix equation P = A (P P B(B P B + R) BP )A + Q () within the class of positive semidefinite symmetric matrices..3 Finite-Time Constrained LQR (FTCLQR) Assume now that the states and the inputs of system () are subject to the following constraints x H (a) u D (b) where H and D are polyhedral sets, and consider the finite-time constrained optimal control problem (FTCLQR) N JN (x()) = min { U x(n) P x(n) + u(i) Ru(i) + x(i) Qx(i)} (3a) subj. to x(i) H i [,..., N] (3b) i= u(i) D i [,..., N ] (3c) The solution to problem (3) has been studied by Bemporad et.al. in [?]. We will briefly summarize the main results. By substituting x(i) = A i x() + i j= Aj Bu(i j), problem (3) can be reformulated as J (x()) = x() Y x() + min U { U HU + x() F U} (4a) s.t. GU W + Ex() (4b) where the column vector U [u(),..., u(n ) ] R s, is the optimization vector, H = H, and H, F, Y, G, W, E are easily obtained from Q, R, P, () and () (see [?] for details). We denote with X f R n the set of initial states x for which the optimal control problem (3) is feasible, i.e. X f = {x R n U R s, GU W + Ex} and with U (x) the optimizer of (4) for x() = x. Before going further, we will introduce the following definitions Definition We define a polyhedral partition as follows... Definition The set of active constraints A(x) at point x of problem (3) is defined as A(x) = {i I G i U (x) W i E i x = }, I = {,,..., m G } (5) whereby G i, W i and E i denote the i-th row of the matrices G, W and E, respectively. Problem (4) depends on x(), therefore the implementation of FTCLQR can be performed either by: solving the quadratic program (4) for a given x(), or as shown in [?,?], by solving problem (4) for all x() within a given range of values, i.e., by considering (4) as a multi-parametric Quadratic Program (mp-qp). Analysis of the mp-qp solution [?] leads to the following results on the solution of (3)
3 Theorem Consider the FTCLQR (3). Then, the set of feasible parameters X f is convex, the optimizer U : X f R s is continuous and piecewise affine, i.e. U (x()) = F k x() + G k if x() P k = {x R n H k x K k }, k =,..., N r (6) and the optimal solution J : X f R is continuous, convex and piecewise quadratic. Furthermore, {P k } N r k= is a polyhedral partition of X f and the set of active constraints in each P k is constant, i.e. A(x) = A k, x P k, A k I, k =,..., N r Once a certain set of active constraints A k is identified, the authors in [?] use KKT conditions to compute the polyhedron P k and the feedback law F k, G k as described next. First, new matrices G Ak, W Ak and S Ak are formed by extracting the rows indexed by A k from G, W and S, where S = E + GH F. Then, the feedback law is given by F k = H G A k (G Ak H G A k ) S Ak H F (7) G k = H G A k (G Ak H G A k ) W Ak (8) (9) and the polyhedron P k is computed as [ ] G(Fk + H H k = F ) S (G Ak H G A k ) S Ak [ ] W GG K k = k (G Ak H G A k ) W Ak () () Various methods to identify the sets of active constraints A k, without enumerating them all, are proposed in the literature....4 Infinite-Time Constrained LQR (CLQR) If in (4) we set N = + and we assume that (A, B) stabilizable and (A, C) detectable (Q = C C), then we obtain the infinite-time constrained LQR (CLQR) problem: J (x()) = x() Y x() + min U { U HU + x() F U} (a) s.t. GU W + Ex() (b) where the column vector U [u(), u(),...] R s, is the optimization vector and all other matrices are computed as in (4). The following theorem is needed to establish the solutions properties of the CLQR: Theorem Given system () with (A, B) stabilizable and (C, A) detectable, there exists a finite horizon ˆN(x()) such that the finite-time solution to () and (3) matches the infinite-time solution of (8) and () for a give x(). Proof cite Kalman(?) or some other guy... In view of the results of the previous section and Theorem, the implementation of CLQR can be performed either by: solving the quadratic program (4) for a given x() with N = ˆN(x()), or we can solve problem () in a given compact set S of the initial conditions by solving the mp-qp (4) for N = ˆN S max x() S ˆN(x()). Theorem 3 Given a compact set S, the horizon ˆN S is finite and by theorem the state feedback solution of problem () in S is PWA. 3
4 Various methods have been proposed in literature which allow an estimation or computation of ˆNS and ˆN(x()). Chmielewski and Manousiouthakis presented an approach, which allows for a conservative estimate of the finite horizon N est ˆN S for a compact set S. They solve a single, finite dimensional, convex program of known size to obtain N est. Rawlings et.al. presented an algorithm which iteratively attempts to identify ˆN(x()). The key definitions and theorems will be reformulated here for completeness. Definition 3 X K R n denotes the set of states x t for which the control law given in (9) satisfies the constraints in () for all time. If the constraints in () are symmetric, X K is a positive invariant set containing an open neighborhood of the origin (cite 6 in rawlings). Theorem 4 The control laws and the associated cost for any given initial state described by (4) are equal to those for the infinite horizon case in () if the terminal state x(n) lies in the unconstrained positive invariant set X K for (4), i.e. if x(n) X K, then J N (x()) = J (x()). Proof This holds because the terminal cost x(n) P x(n) computed from () is equal to the cost for the horizon segment N [N +, ], if and only if x N X K : [ x(n) P x(n) = min U i=n ] x(i) Qx(i) + u(i) Ru(i) x N X K (3) Rawlings method solves the QP for a set horizon N. If the final state is not in X K m N is increased according to a predefined iteration law. One approach would be to continually increment N by, which would yield the minimal horizon N such that x(n) X K. An alternative is to increase N by a factor of at each iteration,i.e. N = N, which results in fewer QPs to be solved at the cost of a larger ˆN(x()). 3 Comparison of Available Techniques for solving the CLQR Based on two examples, some of the drawbacks of the approaches of Rawlings (cite) and Manousiouthakis (cite) will be illustrated. Example We will consider the following system presented in moraribemporad y(t) = s u(t) (4) + 3s + sampled with T =.s to obtain the state-space representation: ( ) ( ) x(t + ) = x(t) + u(t) (5) y(t) = [.44] x(t) The task is to regulate the system to the origin while fulfilling the input constraint u Example We will consider a standard double-integrator system y(t) = u(t) (6) s 4
5 and its equivalent discrete-time state-space representation: ( ) ( ) x(t + ) = x(t) + u(t) (7) y(t) = [ ] x(t) The task is to regulate the system to the origin while fulfilling the input constraint u In order to be able to apply the Manousiouthakis approach, we will first define the compact set S as S M = {x R 5 x i 5, i =, }. Using the method in (Manousiouthakis) we obtain N est = for Example and N est = 743 for Example. Obviously, solving quadratic programs with these horizons will not yield any results in reasonable time. Moreover, the true, non-conservative minimal horizon is ˆN S = 5 for Example and ˆN S = for Example. Rawlings approach is well suited for small compact sets. However, if the set S is larger, e.g. S R = {x R x i, i =, }, the necessary minimal horizon grows quickly, drastically increasing run-times. For random values of x() S R the run-times of the algorithm are presented for an increment policy of N = N + and N = N in the table below: W orst Case Average Case N = N + sec 5 sec N = N 8 sec 6 sec While May be insert comparison with N of manu.. and N of rawlings..., and tehrefore motivate our approach: that is: off-line computation of the PWA but with small N. Chmielewski and Manousiouthakis presented an approach, which allows for a conservative estimate of the finite horizon N. 4 Computation of the Constrained Infinite Time Linear Quadratic Optimal Control Problem In this section, the main contribution of this paper is presented. Our procedure is initialized by computing the positive-invariant unconstrained set X K such that: for x(t) X K x(t + ) X K if x(t + ) = Ax(t) + Bu(t) u(t) = K LQR x(t) The polyhedron X K can be computed by continually increasing the prediction horizon from N = until an additional increase has no effect on X K, i.e. until X K converges (reference Blanchini). X K is depicted in Figure (a). We can describe the polyhderon X K through the matrix-inequalities H X x K X. Note that this approach only works for the unconstrained region. The exploration algorithm finds a point x by stepping over a facet of X K with step-size ɛ (ɛ ). The solution of (??) with horizon N = for x = x when pb?? is fixed fix x to have the same name provides the active constraints that define the neighboring critical region CR (see Figure (b)). By Theorem 4, the finite time optimal solution computed above equals the infinite time optimal solution if x(n) X K, therefore we extract from CR the set of points which will enter X K at the next time-step, provided that 5
6 Feasible Region Fragments 3 Feasible Region Fragments 3 Feasible Region Fragments.5.5 Region for Horizon N=.5 Extract Reachability SubsetITCR Step over Favet θ θ θ epsilon step over facet Pos. Inv. Unconstrained Region Pos. Inv. Unconstrained Region θ θ θ (a) Compute pos.-inv. region X K and step over facet with step-size ɛ (b) Solve QP for new point with horizon N to create the CR (c) Compute reachability subset ITCR Figure : Constrained Infinite-Time Region Exploration the optimal control-law associated with CR (u (t) = F x(t) + G ) is applied. According to Lemma?? the open loop prediction x(t + N) = A N x(t) + [A N B,..., AB, B](F x(t) + G) (8) equals the closed loop response, and the Infinite-Time Critical Region (ITCR) is defined by the intersection of the following polyhedrons and half-spaces: H X x(t + N) K X (9a) Gz W + Sx (9b) λ (9c) Whereby (9b) and (9c) are identical to (??) and (??) which were derived in Section.3. The first identified region will be referred to as ITCR (see Figure (c)) according to the following convention. Definition 4 The region ITCR N i is the i-th Infinite-Time-Critical-Region, computed by optimizing over the horizon N and satisfying the N-steps reachability condition 9a. In the following, we will assume that ITCR N i is defined as the non-redundant set of inequalities Hi N x KN i. This polyhedron has two types of facets: Type I: The facet originated from constraint restrictions in (9b) and (9c). Type II: The facet originated from reachability restrictions in (9a). The exploration procedure continues for all the facets of Type I according to the same scheme: Subsequently this procedure is repeated for all facets, which originate from constraint restrictions (Type I). This also applies to all facets of X K. Type II facets are not considered, since, by definition, those constraints cut off the points that will not reach X K in N steps. Once all facets have been explored for the finite horizon N, the horizon is increased to N + and the entire procedure is restarted. For all regions which have been previously computed, only facets of Type II are considered while for all newly identified regions, only facets of Type I are examined. Theorem 5 The intersection of the interior of ITCR N i and ITCR M j is nonempty, if and only if i = j and M = N. Proof 3 If trivial. Only If, from Theorem 4 and (9), we can conclude that the ITCR region partition is identical to the infinite-time region partition. Therefore each region has a distinct set of active constraints and from Alberto we can conclude that if two IT CRs have a non-empty intersection then they are identical. 6
7 x Trajectory for initial state [.8;.8] Horizon = Horizon = Horizon = 3 Horizon = 4 Horizon = 5 Horizon = 6 Horizon = 7 Horizon = 8 Horizon = 9 Horizon = Horizon = Horizon = Horizon = 3 Horizon = 4 Horizon = x Figure : State trajectory for x = [.8;.8]. We summarize the ideas from Section and this Section in the following algorithm: Algorithm : function R = computec-lqr. ITCR = X K, N =, R = {};. R = R {ITCR, K LQR }; 3. forall ITCR M i R, M < N, explore(itcr M i,type II ); 4. forall ITCR M i R, M = N, explore(itcr M i,type I ); 5. if N < N max, N = N + goto 3, else, return(r). Algorithm : function explore(itcr,type);. forall f facets(itcr),. if f Type, goto 9; 3. Step over f and get x; 4. if ITCR N i R s.t. x ITCR N i, i = i +, goto 9; 5. solve (??) for x = x; 6. compute ITCR new according to (9); 7. compute F and G according to (??); 8. R = R {ITCR new, F, G}; 9. end forall. If the system is subject to symmetric constraints, the computation time can be almost halved, since the region partition is also symmetric to the origin (cite tondels friend). 7
8 R_ Reachability N R_ R_3 R_4 Reachability N- Figure 3: Region R and R are connected to R 4 and region R 3 is adjacent to R 4. Theorem 6 If all regions with a horizon smaller than N have been computed using algorithm, the set of those regions is positive invariant. Proof 4 Follows from Theorem??. The state will always move to a region with horizon N at the next time step until the unconstrained region X K is reached. Theorem 7 In the ITCOC region partition, each region with a reachability of N is adjacent or connected to a region with a reachability of N. We define two regions as connected, if they re not directly adjacent but reachable by moving through adjacent regions with same horizon as the originating region (see Figure 3). Proof 5 Assume that the entire region partition of set S for an arbitrary horizon of N is computed. Subsequently we compute the reachability subset for a horizon of N over the entire set S. Now we recompute the partition and the reachability set for a horizon of N. Obviously the reachability set for horizon N is a subset for the reachability set for horizon N. Therefore all infinite horizon regions with reachability N are adjacent are connected to a region with reachability N. By computing all regions with a reachability horizon of N = and by applying Theorem 7, we can show that the presented algorithm provides the minimal solution to the problem at hand. The solution is minimal in the sense that the smallest possible horizon is always used to compute the QPs, i.e. if a given point can reach X K in N steps, the QP solved for that point will never be computed with a horizon larger than N. Therefore this algorithm is also an efficient method to compute ˆN S for any set sets without having to rely on conservative estimates. 5 Results Figure 4 shows the regions which were computed by using Algorithm on the double integrator example presented in morari et.al.. The times needed to compute the regions using various algorithms with a.ghz System are given in the table below. Inf-MPQP denotes algorithm presented in the previous section; MPQP using N min signifies that the MPQP presented in Section 4 is solved for the lowest possible horizon setting; MPQP using N estimate signifies the the horizon estimates from Chmielewski et.al. were adopted. 8
9 θ Feasible Region Fragments Horizon = Horizon = Horizon = 3 Horizon = 4 Horizon = 5 Horizon = 6 Horizon = 7 Horizon = 8 Horizon = 9 Horizon = Horizon = Horizon = Horizon = 3 Horizon = 4 Horizon = θ Figure 4: Regions computed by Algorithm θ Feasible Region Fragments Horizon = Horizon = Horizon = 3 Horizon = 4 Horizon = 5 Horizon = 6 Horizon = 7 Horizon = 8 Horizon = 9 Horizon = Horizon = Horizon = Horizon = 3 Horizon = 4 Horizon = θ Figure 5: Computation of the invariant set which covers the black 5x5 region. 9
10 explored state-space Inf-MPQP MPQP using N min MPQP using N estimate 5x5.5 sec sec 8 sec x sec sec sec tracking not possible could be used in combination with rawlings point to numerical problems? 6 Conclusions
An Algorithm for Multi-Parametric Quadratic Programming and Explicit MPC Solutions
An Algorithm for Multi-Parametric Quadratic Programming and Explicit MPC Solutions P. Tøndel 1, T.A. Johansen 1, A. Bemporad 2 TuP11-4 Abstract Explicit solutions to constrained linear MPC problems can
More informationRELATIVELY OPTIMAL CONTROL: THE STATIC SOLUTION
RELATIVELY OPTIMAL CONTROL: THE STATIC SOLUTION Franco Blanchini,1 Felice Andrea Pellegrino Dipartimento di Matematica e Informatica Università di Udine via delle Scienze, 208 33100, Udine, Italy blanchini@uniud.it,
More informationLecture 2 September 3
EE 381V: Large Scale Optimization Fall 2012 Lecture 2 September 3 Lecturer: Caramanis & Sanghavi Scribe: Hongbo Si, Qiaoyang Ye 2.1 Overview of the last Lecture The focus of the last lecture was to give
More informationOutline. Robust MPC and multiparametric convex programming. A. Bemporad C. Filippi. Motivation: Robust MPC. Multiparametric convex programming
Robust MPC and multiparametric convex programming A. Bemporad C. Filippi D. Muñoz de la Peña CC Meeting Siena /4 September 003 Outline Motivation: Robust MPC Multiparametric convex programming Kothares
More informationEfficient On-Line Computation of Constrained Optimal Control
Efficient On-Line Computation of Constrained Optimal Control Mato Baotić, Francesco Borrelli, Alberto Bemporad, Manfred Morari Automatic Control Laboratory, ETH Zentrum - ETL, Physikstrasse 3, CH-8092
More informationThe Facet-to-Facet Property of Solutions to Convex Parametric Quadratic Programs and a new Exploration Strategy
The Facet-to-Facet Property of Solutions to Convex Parametric Quadratic Programs and a new Exploration Strategy Jørgen Spjøtvold, Eric C. Kerrigan, Colin N. Jones, Petter Tøndel and Tor A. Johansen Abstract
More informationInfinite Time Optimal Control of Hybrid Systems with a Linear Performance Index
Infinite Time Optimal Control of Hybrid Systems with a Linear Performance Index Mato Baotić, Frank J. Christophersen, and Manfred Morari Automatic Control Laboratory, ETH Zentrum, ETL K 1, CH 9 Zürich,
More information4 LINEAR PROGRAMMING (LP) E. Amaldi Fondamenti di R.O. Politecnico di Milano 1
4 LINEAR PROGRAMMING (LP) E. Amaldi Fondamenti di R.O. Politecnico di Milano 1 Mathematical programming (optimization) problem: min f (x) s.t. x X R n set of feasible solutions with linear objective function
More informationOn the Complexity of Explicit MPC Laws
On the Complexity of Explicit MPC Laws Francesco Borrelli, Mato Baotić, Jaroslav Pekar and Greg Stewart Abstract Finite-time optimal control problems with quadratic performance index for linear systems
More informationOn the facet-to-facet property of solutions to convex parametric quadratic programs
Automatica 42 (2006) 2209 2214 www.elsevier.com/locate/automatica Technical communique On the facet-to-facet property of solutions to convex parametric quadratic programs JZrgen SpjZtvold a,, Eric C. Kerrigan
More informationAdvanced Operations Research Techniques IE316. Quiz 1 Review. Dr. Ted Ralphs
Advanced Operations Research Techniques IE316 Quiz 1 Review Dr. Ted Ralphs IE316 Quiz 1 Review 1 Reading for The Quiz Material covered in detail in lecture. 1.1, 1.4, 2.1-2.6, 3.1-3.3, 3.5 Background material
More informationInteger Programming Theory
Integer Programming Theory Laura Galli October 24, 2016 In the following we assume all functions are linear, hence we often drop the term linear. In discrete optimization, we seek to find a solution x
More informationComplexity Reduction of Explicit Model Predictive Control via Combining Separator Function and Binary Search Trees
American Journal of Computer Science and Technology 2018; 1(1): 19-23 http://www.sciencepublishinggroup.com/j/ajcst doi: 10.11648/j.ajcst.20180101.13 Complexity Reduction of Explicit Model Predictive Control
More informationFurther Results on Multiparametric Quadratic Programming
Further Results on Multiparametric Quadratic Programming P. Tøndel, T.. Johansen,. Bemporad bstract In this paper we extend recent results on strictly convex multiparametric quadratic programming (mpqp)
More informationLesson 17. Geometry and Algebra of Corner Points
SA305 Linear Programming Spring 2016 Asst. Prof. Nelson Uhan 0 Warm up Lesson 17. Geometry and Algebra of Corner Points Example 1. Consider the system of equations 3 + 7x 3 = 17 + 5 = 1 2 + 11x 3 = 24
More informationControlling Hybrid Systems
Controlling Hybrid Systems From Theory to Application Manfred Morari M. Baotic, F. Christophersen, T. Geyer, P. Grieder, M. Kvasnica, G. Papafotiou Lessons learned from a decade of Hybrid System Research
More informationLecture Notes 2: The Simplex Algorithm
Algorithmic Methods 25/10/2010 Lecture Notes 2: The Simplex Algorithm Professor: Yossi Azar Scribe:Kiril Solovey 1 Introduction In this lecture we will present the Simplex algorithm, finish some unresolved
More informationComputation of the Constrained Infinite Time Linear Quadratic Regulator
Computation of the Constrained Infinite Time Linear Quadratic Regulator P. Grieder' '. F. Borrelli'. F. Torrisi' and M. Morari' Abstract This paper presents an efficient algorithin for coiiiputing the
More informationCS675: Convex and Combinatorial Optimization Spring 2018 Convex Sets. Instructor: Shaddin Dughmi
CS675: Convex and Combinatorial Optimization Spring 2018 Convex Sets Instructor: Shaddin Dughmi Outline 1 Convex sets, Affine sets, and Cones 2 Examples of Convex Sets 3 Convexity-Preserving Operations
More informationConvexity: an introduction
Convexity: an introduction Geir Dahl CMA, Dept. of Mathematics and Dept. of Informatics University of Oslo 1 / 74 1. Introduction 1. Introduction what is convexity where does it arise main concepts and
More informationChapter 4 Concepts from Geometry
Chapter 4 Concepts from Geometry An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Line Segments The line segment between two points and in R n is the set of points on the straight line joining
More informationIntroduction to optimization
Introduction to optimization G. Ferrari Trecate Dipartimento di Ingegneria Industriale e dell Informazione Università degli Studi di Pavia Industrial Automation Ferrari Trecate (DIS) Optimization Industrial
More informationExplicit Model Predictive Control by A. Bemporad Controllo di Processo e dei Sistemi di Produzione A.a. 2008/09 1/54
Explicit Model Predictive Control 1/54 KKT Conditions for Optimality When f, g i are convex functions and h j are linear, the condition is also sufficient 2/54 KKT Geometric Interpretation rg 1 (U 1 )
More informationPOLYHEDRAL GEOMETRY. Convex functions and sets. Mathematical Programming Niels Lauritzen Recall that a subset C R n is convex if
POLYHEDRAL GEOMETRY Mathematical Programming Niels Lauritzen 7.9.2007 Convex functions and sets Recall that a subset C R n is convex if {λx + (1 λ)y 0 λ 1} C for every x, y C and 0 λ 1. A function f :
More informationDM545 Linear and Integer Programming. Lecture 2. The Simplex Method. Marco Chiarandini
DM545 Linear and Integer Programming Lecture 2 The Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark Outline 1. 2. 3. 4. Standard Form Basic Feasible Solutions
More informationLecture 2 - Introduction to Polytopes
Lecture 2 - Introduction to Polytopes Optimization and Approximation - ENS M1 Nicolas Bousquet 1 Reminder of Linear Algebra definitions Let x 1,..., x m be points in R n and λ 1,..., λ m be real numbers.
More informationOn the Relationships between Zero Forcing Numbers and Certain Graph Coverings
On the Relationships between Zero Forcing Numbers and Certain Graph Coverings Fatemeh Alinaghipour Taklimi, Shaun Fallat 1,, Karen Meagher 2 Department of Mathematics and Statistics, University of Regina,
More information3 No-Wait Job Shops with Variable Processing Times
3 No-Wait Job Shops with Variable Processing Times In this chapter we assume that, on top of the classical no-wait job shop setting, we are given a set of processing times for each operation. We may select
More informationCS675: Convex and Combinatorial Optimization Spring 2018 The Simplex Algorithm. Instructor: Shaddin Dughmi
CS675: Convex and Combinatorial Optimization Spring 2018 The Simplex Algorithm Instructor: Shaddin Dughmi Algorithms for Convex Optimization We will look at 2 algorithms in detail: Simplex and Ellipsoid.
More informationBilinear Programming
Bilinear Programming Artyom G. Nahapetyan Center for Applied Optimization Industrial and Systems Engineering Department University of Florida Gainesville, Florida 32611-6595 Email address: artyom@ufl.edu
More informationMath 5593 Linear Programming Lecture Notes
Math 5593 Linear Programming Lecture Notes Unit II: Theory & Foundations (Convex Analysis) University of Colorado Denver, Fall 2013 Topics 1 Convex Sets 1 1.1 Basic Properties (Luenberger-Ye Appendix B.1).........................
More informationIntroduction to Optimization Problems and Methods
Introduction to Optimization Problems and Methods wjch@umich.edu December 10, 2009 Outline 1 Linear Optimization Problem Simplex Method 2 3 Cutting Plane Method 4 Discrete Dynamic Programming Problem Simplex
More informationCS599: Convex and Combinatorial Optimization Fall 2013 Lecture 4: Convex Sets. Instructor: Shaddin Dughmi
CS599: Convex and Combinatorial Optimization Fall 2013 Lecture 4: Convex Sets Instructor: Shaddin Dughmi Announcements New room: KAP 158 Today: Convex Sets Mostly from Boyd and Vandenberghe. Read all of
More informationLecture notes on the simplex method September We will present an algorithm to solve linear programs of the form. maximize.
Cornell University, Fall 2017 CS 6820: Algorithms Lecture notes on the simplex method September 2017 1 The Simplex Method We will present an algorithm to solve linear programs of the form maximize subject
More informationPower Set of a set and Relations
Power Set of a set and Relations 1 Power Set (1) Definition: The power set of a set S, denoted P(S), is the set of all subsets of S. Examples Let A={a,b,c}, P(A)={,{a},{b},{c},{a,b},{b,c},{a,c},{a,b,c}}
More informationProgramming, numerics and optimization
Programming, numerics and optimization Lecture C-4: Constrained optimization Łukasz Jankowski ljank@ippt.pan.pl Institute of Fundamental Technological Research Room 4.32, Phone +22.8261281 ext. 428 June
More information60 2 Convex sets. {x a T x b} {x ã T x b}
60 2 Convex sets Exercises Definition of convexity 21 Let C R n be a convex set, with x 1,, x k C, and let θ 1,, θ k R satisfy θ i 0, θ 1 + + θ k = 1 Show that θ 1x 1 + + θ k x k C (The definition of convexity
More informationR n a T i x = b i} is a Hyperplane.
Geometry of LPs Consider the following LP : min {c T x a T i x b i The feasible region is i =1,...,m}. X := {x R n a T i x b i i =1,...,m} = m i=1 {x Rn a T i x b i} }{{} X i The set X i is a Half-space.
More informationMathematical Programming and Research Methods (Part II)
Mathematical Programming and Research Methods (Part II) 4. Convexity and Optimization Massimiliano Pontil (based on previous lecture by Andreas Argyriou) 1 Today s Plan Convex sets and functions Types
More informationConvex Geometry arising in Optimization
Convex Geometry arising in Optimization Jesús A. De Loera University of California, Davis Berlin Mathematical School Summer 2015 WHAT IS THIS COURSE ABOUT? Combinatorial Convexity and Optimization PLAN
More informationOpen and Closed Sets
Open and Closed Sets Definition: A subset S of a metric space (X, d) is open if it contains an open ball about each of its points i.e., if x S : ɛ > 0 : B(x, ɛ) S. (1) Theorem: (O1) and X are open sets.
More informationLinear Optimization. Andongwisye John. November 17, Linkoping University. Andongwisye John (Linkoping University) November 17, / 25
Linear Optimization Andongwisye John Linkoping University November 17, 2016 Andongwisye John (Linkoping University) November 17, 2016 1 / 25 Overview 1 Egdes, One-Dimensional Faces, Adjacency of Extreme
More informationEfficient On-Line Computation of Constrained Optimal Control
Efficient On-Line Computation of Constrained Optimal Control Francesco Borrelli, Mato Baotic, Alberto Bemporad, Manfred Morari Automatic Control Laboratory, ETH Zentrum - ETL, CH-89 Zürich, Switzerland
More informationLinear Machine: A Novel Approach to Point Location Problem
Preprints of the 10th IFAC International Symposium on Dynamics and Control of Process Systems The International Federation of Automatic Control Linear Machine: A Novel Approach to Point Location Problem
More informationSystems of Inequalities
Systems of Inequalities 1 Goals: Given system of inequalities of the form Ax b determine if system has an integer solution enumerate all integer solutions 2 Running example: Upper bounds for x: (2)and
More informationChapter 15 Introduction to Linear Programming
Chapter 15 Introduction to Linear Programming An Introduction to Optimization Spring, 2015 Wei-Ta Chu 1 Brief History of Linear Programming The goal of linear programming is to determine the values of
More informationAdvanced Operations Research Techniques IE316. Quiz 2 Review. Dr. Ted Ralphs
Advanced Operations Research Techniques IE316 Quiz 2 Review Dr. Ted Ralphs IE316 Quiz 2 Review 1 Reading for The Quiz Material covered in detail in lecture Bertsimas 4.1-4.5, 4.8, 5.1-5.5, 6.1-6.3 Material
More information22 Elementary Graph Algorithms. There are two standard ways to represent a
VI Graph Algorithms Elementary Graph Algorithms Minimum Spanning Trees Single-Source Shortest Paths All-Pairs Shortest Paths 22 Elementary Graph Algorithms There are two standard ways to represent a graph
More informationSchool of Management University of Shanghai for Science and Technology Shanghai , China
1 ON THE COMPUTATION OF VIABLE POLYTOPES FOR LINEAR SYSTEMS Yan Gao, John Lygeros School of Management University of Shanghai for Science and Technology Shanghai 293, China email:gaoyan1962@263.net Department
More informationarxiv: v1 [cs.cc] 30 Jun 2017
On the Complexity of Polytopes in LI( Komei Fuuda May Szedlá July, 018 arxiv:170610114v1 [cscc] 30 Jun 017 Abstract In this paper we consider polytopes given by systems of n inequalities in d variables,
More informationConvex Optimization - Chapter 1-2. Xiangru Lian August 28, 2015
Convex Optimization - Chapter 1-2 Xiangru Lian August 28, 2015 1 Mathematical optimization minimize f 0 (x) s.t. f j (x) 0, j=1,,m, (1) x S x. (x 1,,x n ). optimization variable. f 0. R n R. objective
More information2. Convex sets. affine and convex sets. some important examples. operations that preserve convexity. generalized inequalities
2. Convex sets Convex Optimization Boyd & Vandenberghe affine and convex sets some important examples operations that preserve convexity generalized inequalities separating and supporting hyperplanes dual
More informationIn this chapter we introduce some of the basic concepts that will be useful for the study of integer programming problems.
2 Basics In this chapter we introduce some of the basic concepts that will be useful for the study of integer programming problems. 2.1 Notation Let A R m n be a matrix with row index set M = {1,...,m}
More information16.410/413 Principles of Autonomy and Decision Making
16.410/413 Principles of Autonomy and Decision Making Lecture 17: The Simplex Method Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology November 10, 2010 Frazzoli (MIT)
More informationLinear Bilevel Programming With Upper Level Constraints Depending on the Lower Level Solution
Linear Bilevel Programming With Upper Level Constraints Depending on the Lower Level Solution Ayalew Getachew Mersha and Stephan Dempe October 17, 2005 Abstract Focus in the paper is on the definition
More informationResearch Topics (Baotic, Bemporad, Borrelli, Ferrari-Trecate, Geyer, Grieder, Mignone, Torrisi, Morari)
Research Topics (Baotic, Bemporad, Borrelli, Ferrari-Trecate, Geyer, Grieder, Mignone, Torrisi, Morari) Analysis Reachability/Verification Stability Observability Synthesis Control (MPC) Explicit PWA MPC
More information3. The Simplex algorithmn The Simplex algorithmn 3.1 Forms of linear programs
11 3.1 Forms of linear programs... 12 3.2 Basic feasible solutions... 13 3.3 The geometry of linear programs... 14 3.4 Local search among basic feasible solutions... 15 3.5 Organization in tableaus...
More informationGraphs that have the feasible bases of a given linear
Algorithmic Operations Research Vol.1 (2006) 46 51 Simplex Adjacency Graphs in Linear Optimization Gerard Sierksma and Gert A. Tijssen University of Groningen, Faculty of Economics, P.O. Box 800, 9700
More information2. Optimization problems 6
6 2.1 Examples... 7... 8 2.3 Convex sets and functions... 9 2.4 Convex optimization problems... 10 2.1 Examples 7-1 An (NP-) optimization problem P 0 is defined as follows Each instance I P 0 has a feasibility
More informationLecture 2: August 29, 2018
10-725/36-725: Convex Optimization Fall 2018 Lecturer: Ryan Tibshirani Lecture 2: August 29, 2018 Scribes: Adam Harley Note: LaTeX template courtesy of UC Berkeley EECS dept. Disclaimer: These notes have
More information15. Cutting plane and ellipsoid methods
EE 546, Univ of Washington, Spring 2012 15. Cutting plane and ellipsoid methods localization methods cutting-plane oracle examples of cutting plane methods ellipsoid method convergence proof inequality
More information22 Elementary Graph Algorithms. There are two standard ways to represent a
VI Graph Algorithms Elementary Graph Algorithms Minimum Spanning Trees Single-Source Shortest Paths All-Pairs Shortest Paths 22 Elementary Graph Algorithms There are two standard ways to represent a graph
More informationLinear and Integer Programming :Algorithms in the Real World. Related Optimization Problems. How important is optimization?
Linear and Integer Programming 15-853:Algorithms in the Real World Linear and Integer Programming I Introduction Geometric Interpretation Simplex Method Linear or Integer programming maximize z = c T x
More information2. Convex sets. x 1. x 2. affine set: contains the line through any two distinct points in the set
2. Convex sets Convex Optimization Boyd & Vandenberghe affine and convex sets some important examples operations that preserve convexity generalized inequalities separating and supporting hyperplanes dual
More informationA set-based approach to robust control and verification of piecewise affine systems subject to safety specifications
Dipartimento di Elettronica, Informazione e Bioingegneria A set-based approach to robust control and verification of piecewise affine systems subject to safety specifications Maria Prandini maria.prandini@polimi.it
More informationSocial-Network Graphs
Social-Network Graphs Mining Social Networks Facebook, Google+, Twitter Email Networks, Collaboration Networks Identify communities Similar to clustering Communities usually overlap Identify similarities
More informationEfficient implementation of Constrained Min-Max Model Predictive Control with Bounded Uncertainties
Efficient implementation of Constrained Min-Max Model Predictive Control with Bounded Uncertainties D.R. Ramírez 1, T. Álamo and E.F. Camacho2 Departamento de Ingeniería de Sistemas y Automática, Universidad
More informationConvex Optimization. Convex Sets. ENSAE: Optimisation 1/24
Convex Optimization Convex Sets ENSAE: Optimisation 1/24 Today affine and convex sets some important examples operations that preserve convexity generalized inequalities separating and supporting hyperplanes
More informationAlternating Projections
Alternating Projections Stephen Boyd and Jon Dattorro EE392o, Stanford University Autumn, 2003 1 Alternating projection algorithm Alternating projections is a very simple algorithm for computing a point
More informationConstrained Optimal Control of Discrete-Time Linear Hybrid Systems
Constrained Optimal Control of Discrete-Time Linear Hybrid Systems Francesco Borrelli, Mato Baotic, Alberto Bemporad, Manfred Morari Automatic Control Laboratory, ETH Zentrum, ETL K2, CH-892 Zurich, Switzerland
More informationParameterization. Michael S. Floater. November 10, 2011
Parameterization Michael S. Floater November 10, 2011 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to generate from point
More informationConvex Optimization. 2. Convex Sets. Prof. Ying Cui. Department of Electrical Engineering Shanghai Jiao Tong University. SJTU Ying Cui 1 / 33
Convex Optimization 2. Convex Sets Prof. Ying Cui Department of Electrical Engineering Shanghai Jiao Tong University 2018 SJTU Ying Cui 1 / 33 Outline Affine and convex sets Some important examples Operations
More informationLecture 15: The subspace topology, Closed sets
Lecture 15: The subspace topology, Closed sets 1 The Subspace Topology Definition 1.1. Let (X, T) be a topological space with topology T. subset of X, the collection If Y is a T Y = {Y U U T} is a topology
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More information1 Inference for Boolean theories
Scribe notes on the class discussion on consistency methods for boolean theories, row convex constraints and linear inequalities (Section 8.3 to 8.6) Speaker: Eric Moss Scribe: Anagh Lal Corrector: Chen
More informationPRIMAL-DUAL INTERIOR POINT METHOD FOR LINEAR PROGRAMMING. 1. Introduction
PRIMAL-DUAL INTERIOR POINT METHOD FOR LINEAR PROGRAMMING KELLER VANDEBOGERT AND CHARLES LANNING 1. Introduction Interior point methods are, put simply, a technique of optimization where, given a problem
More information9.5 Equivalence Relations
9.5 Equivalence Relations You know from your early study of fractions that each fraction has many equivalent forms. For example, 2, 2 4, 3 6, 2, 3 6, 5 30,... are all different ways to represent the same
More informationHowever, this is not always true! For example, this fails if both A and B are closed and unbounded (find an example).
98 CHAPTER 3. PROPERTIES OF CONVEX SETS: A GLIMPSE 3.2 Separation Theorems It seems intuitively rather obvious that if A and B are two nonempty disjoint convex sets in A 2, then there is a line, H, separating
More informationPolytopes Course Notes
Polytopes Course Notes Carl W. Lee Department of Mathematics University of Kentucky Lexington, KY 40506 lee@ms.uky.edu Fall 2013 i Contents 1 Polytopes 1 1.1 Convex Combinations and V-Polytopes.....................
More informationCSE 20 DISCRETE MATH WINTER
CSE 20 DISCRETE MATH WINTER 2016 http://cseweb.ucsd.edu/classes/wi16/cse20-ab/ Today's learning goals Explain the steps in a proof by (strong) mathematical induction Use (strong) mathematical induction
More information! Linear programming"! Duality "! Smallest enclosing disk"
! Linear programming"! Duality "! Smallest enclosing disk" 14. 24.! Define:" " i types of foods (1!i!d).! " j types of vitamins (1!j!n)." " x i the amount of food of type i." " a ji the amount of vitamin
More informationmaximize c, x subject to Ax b,
Lecture 8 Linear programming is about problems of the form maximize c, x subject to Ax b, where A R m n, x R n, c R n, and b R m, and the inequality sign means inequality in each row. The feasible set
More information15.082J and 6.855J. Lagrangian Relaxation 2 Algorithms Application to LPs
15.082J and 6.855J Lagrangian Relaxation 2 Algorithms Application to LPs 1 The Constrained Shortest Path Problem (1,10) 2 (1,1) 4 (2,3) (1,7) 1 (10,3) (1,2) (10,1) (5,7) 3 (12,3) 5 (2,2) 6 Find the shortest
More informationCS675: Convex and Combinatorial Optimization Spring 2018 Consequences of the Ellipsoid Algorithm. Instructor: Shaddin Dughmi
CS675: Convex and Combinatorial Optimization Spring 2018 Consequences of the Ellipsoid Algorithm Instructor: Shaddin Dughmi Outline 1 Recapping the Ellipsoid Method 2 Complexity of Convex Optimization
More informationConvex sets and convex functions
Convex sets and convex functions Convex optimization problems Convex sets and their examples Separating and supporting hyperplanes Projections on convex sets Convex functions, conjugate functions ECE 602,
More informationLecture 2: August 31
10-725/36-725: Convex Optimization Fall 2016 Lecture 2: August 31 Lecturer: Lecturer: Ryan Tibshirani Scribes: Scribes: Lidan Mu, Simon Du, Binxuan Huang 2.1 Review A convex optimization problem is of
More informationConvex sets and convex functions
Convex sets and convex functions Convex optimization problems Convex sets and their examples Separating and supporting hyperplanes Projections on convex sets Convex functions, conjugate functions ECE 602,
More informationIN the last decade embedded Model Predictive Control
1 Exact Complexity Certification of Active-Set Methods for Quadratic Programming Gionata Cimini, Student Member, IEEE Alberto Bemporad, Fellow, IEEE Abstract Active-set methods are recognized to often
More informationUnconstrained Optimization Principles of Unconstrained Optimization Search Methods
1 Nonlinear Programming Types of Nonlinear Programs (NLP) Convexity and Convex Programs NLP Solutions Unconstrained Optimization Principles of Unconstrained Optimization Search Methods Constrained Optimization
More informationMATH3016: OPTIMIZATION
MATH3016: OPTIMIZATION Lecturer: Dr Huifu Xu School of Mathematics University of Southampton Highfield SO17 1BJ Southampton Email: h.xu@soton.ac.uk 1 Introduction What is optimization? Optimization is
More informationMathematical and Algorithmic Foundations Linear Programming and Matchings
Adavnced Algorithms Lectures Mathematical and Algorithmic Foundations Linear Programming and Matchings Paul G. Spirakis Department of Computer Science University of Patras and Liverpool Paul G. Spirakis
More informationIE 5531: Engineering Optimization I
IE 5531: Engineering Optimization I Lecture 3: Linear Programming, Continued Prof. John Gunnar Carlsson September 15, 2010 Prof. John Gunnar Carlsson IE 5531: Engineering Optimization I September 15, 2010
More informationHW Graph Theory SOLUTIONS (hbovik)
Diestel 1.3: Let G be a graph containing a cycle C, and assume that G contains a path P of length at least k between two vertices of C. Show that G contains a cycle of length at least k. If C has length
More informationApplied Integer Programming
Applied Integer Programming D.S. Chen; R.G. Batson; Y. Dang Fahimeh 8.2 8.7 April 21, 2015 Context 8.2. Convex sets 8.3. Describing a bounded polyhedron 8.4. Describing unbounded polyhedron 8.5. Faces,
More information6 Randomized rounding of semidefinite programs
6 Randomized rounding of semidefinite programs We now turn to a new tool which gives substantially improved performance guarantees for some problems We now show how nonlinear programming relaxations can
More informationLecture 2: August 29, 2018
10-725/36-725: Convex Optimization Fall 2018 Lecturer: Ryan Tibshirani Lecture 2: August 29, 2018 Scribes: Yingjing Lu, Adam Harley, Ruosong Wang Note: LaTeX template courtesy of UC Berkeley EECS dept.
More informationApplied Lagrange Duality for Constrained Optimization
Applied Lagrange Duality for Constrained Optimization Robert M. Freund February 10, 2004 c 2004 Massachusetts Institute of Technology. 1 1 Overview The Practical Importance of Duality Review of Convexity
More informationComputational Methods. Constrained Optimization
Computational Methods Constrained Optimization Manfred Huber 2010 1 Constrained Optimization Unconstrained Optimization finds a minimum of a function under the assumption that the parameters can take on
More informationConic Duality. yyye
Conic Linear Optimization and Appl. MS&E314 Lecture Note #02 1 Conic Duality Yinyu Ye Department of Management Science and Engineering Stanford University Stanford, CA 94305, U.S.A. http://www.stanford.edu/
More informationLECTURE 13: SOLUTION METHODS FOR CONSTRAINED OPTIMIZATION. 1. Primal approach 2. Penalty and barrier methods 3. Dual approach 4. Primal-dual approach
LECTURE 13: SOLUTION METHODS FOR CONSTRAINED OPTIMIZATION 1. Primal approach 2. Penalty and barrier methods 3. Dual approach 4. Primal-dual approach Basic approaches I. Primal Approach - Feasible Direction
More information