Efficient On-Line Computation of Constrained Optimal Control
|
|
- Meryl Thornton
- 5 years ago
- Views:
Transcription
1 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 tel: , fax: +--6, borrelli,baotic,bemporad,morari@aut.ee.ethz.ch Dip. Ingegneria dell Informazione, Università di Siena, Via Roma 56, I-5 Siena, Italy tel: , fax: , bemporad@dii.unisi.it Abstract For discrete-time linear time-invariant systems with constraints on inputs and outputs, the constrained finite time optimal controller can be obtained explicitly as a piecewise affine function of the initial state via multiparametric programming [, ]. Exploiting the properties of the value function, we present two algorithms that efficiently perform the online evaluation of the explicit optimal control law both in terms of storage demands and computational complexity.the algorithms are particularly effective when used for Model Predictive Control (MPC) where an open loop constrained finite time optimal control problem has to be solved at each sampling time. Introduction Recently in [] and in [] the authors have shown how to compute the solution of the constrained finite time optimal control (CFTOC) problem as a explicit piecewise affine function of the initial state.such a function is computed off-line by using a multiparametric programming solver [, ], which divides the state space into polyhedral regions, and for each region determines the linear gain which produces the optimal control action. This method reveals its effectiveness when applied to Model Predictive Control (MPC) [8, 6].MPC requires to solve at each sampling time an open-loop CFTOC problem.the optimal command signal is applied to the process only during the following sampling interval.at the next time step a new optimal control problem based on new measurements of the state is solved over a shifted horizon.clearly the explicit solution reduces the on-line computation of the MPC control law to the evaluation of a piecewise affine function instead of the on-line solution of a quadratic or linear program. The only drawback of such explicit optimal control laws is that the number of polyhedral regions could grow exponentially with the size of the optimal control problem. In this paper we focus on efficient on-line methods for the evaluation of this piecewise affine control law.the simplest algorithm would require: (i) the storage of the list of polyhedral regions and of the corresponding linear control laws, (ii) a sequential search through the list of polyhedra for the i-th polyhedron that contains the current state in order to implement the i-th control law. By exploiting the properties of the value function, for CFTOC based on LP and QP, we propose two new algorithms that avoid storing the polyhedral regions.the new algorithms significantly reduce the on-line storage demands and computational complexity of the evaluation of the explicit CFTOC control law and consequently of the explicit MPC control law. CFTOC, MPC and Their Explicit Solution. CFTOC problem formulation Consider the linear time-invariant system { x t+ = Ax t + Bu t () y t = Cx t subject to the constraints y min y t y max, u min u t u max () at all time instants t. In () (), x t R n, u t R m, and y t R p are the state, input, and output vector respectively, y min y max (u min u max ) are bounds on outputs (inputs), and the pair (A, B) is stabilizable. Let x(t) be the state at the generic time t and consider the constrained finite time optimal control problem N J (x(t)) = min U xt+n P p + (x t+k) Q p + u t+k R p k= subj. to y min y t+k y max, k =,...,N u min u t+k u max, k =,...,N x t+k+ = Ax t+k + Bu t+k, k () In()wedenotewithU [u t,...,u t+n ] the optimization vector and with x M p the p-norm of the vector x weighted with the matrix M, in particular x M p = x Mx for p =, x M p = Mx p for p =,. In the following, we will assume that Q = Q, R = R, P, for p =,andthatq, R, P are full column rank matrices for p =,. Depending on the norm p used in the cost function, the optimization problem () can be translated into an LP Although the form () is very common in practical implementations of CFTOC, the results of this paper also hold for the more general mixed constraints Ex t + Lu t M arising, for example, from constraints on the input rate u t u t u t.
2 (p =orp =+ ) orintoaqp(p = ).The optimizer U (t) ={u t,...,u t+n } of problem () is a function of the initial state x(t).it can be computed by solving aqporanlponcex(t) is fixed or it can be computed explicitly for all x(t) within a given range of values as explained in the following.. MPC formulation The solution U (t) ={u t,...,u t+n } of problem () is an open loop optimal control trajectory over a finite horizon.mpc uses it together with a receding horizon strategy to obtain an infinite horizon feedback control law. Consider the problem of regulating to the origin the discrete-time linear time-invariant system () while fulfilling the constraints ().MPC solves such a constrained regulation problem in the following way.assume that a full measurement of the state x(t) is available at the current time t.then, the CFTOC problem () is solved at each time t for t. Let U (t) ={u t,...,u t+n } be the optimal solution of ().Then at time t u(t) =u t () is applied as input to system (). The two main issues regarding this policy are the feasibility of the optimization problem () for all t and the stability of the resulting closed-loop system. We will assume that the matrices Q, R, P, the horizon length N and the terminal constraints in () have been chosen to guarantee the stability and the feasibility of MPC control law ()-().For a detailed discussion see, e.g. [,, 6].. Explicit Solution Based on LP Consider the case of p =orp =+ in problem (). Is well known (e.g. []) that by introducing the vector z {ε x,...,εx N,εu,...,εu N,u t,...,u t+n } R s, s (m +)N, ε x k Qx t+k, ε x N Px t+n, ε u k Ru t+k, and substituting x t+k = A k x(t) + k j= Aj Bu t+k j in (), the optimization problem () can be rewritten in the form J (x) =min f T z z (5) subj.to Gz S + Fx where x = x(t), and the matrices f R s, G R q s, F R q n, S R q are easily obtained from Q, R, P and (), as explained in []. Because the problem depends on x = x(t) the implementation of MPC can be performed in two different ways: solving the LP (5) on-line at each time step or by solving problem (5) off-line for all x within a given range of values, i.e., by considering (5) as a multi-parametric Linear Program (mp-lp). The mp-lp problem consists of computing the optimizer z (x) and the value function J (x) for all possible vectors x in a given compact set X.The solution of mp-lp problems can be approached as proposed in [] or [5]. Once the multi-parametric problem (5) has been solved off-line for a polyhedral set X R n of states, the CFTOC explicit solution z (x) of (5) is available as a piecewise affine function of x, and the model predictive controller ()-() is also available explicitly, as the optimal input u(t) consists simply of m components of z (x) u(t) =[... I m... ]z (x). (6) In [5] the following results about the properties of the solution are proved: Theorem Consider the multi-parametric linear program (5). Then the set of feasible parameters X f is convex,the optimizer (if unique) z (x) :X f R s is continuous and piecewise affine,and the optimal solution J (x) : X f R is continuous,convex and piecewise linear. Corollary The MPC control law (6) u(t) = f(x), f : R n R m,defined by the optimization problem () and () is continuous and piecewise affine.. Explicit Solution Based on QP Consider the case p = in problem ().By substituting x t+k t = A k x(t) + k j= Aj Bu t+k j in (), the optimization problem () can be rewritten as the QP J (x) = x Yx+ min U U HU + x FU subj.to GU W + Ex (7) where x = x(t), the column vector U [u t,...,u t+n ] R s, s mn, is the optimization vector, H = H, and H, F, Y, G, W, E are easily obtained from Q, R, P and ().As only the optimizer U is needed, the term involving Y is usually removed from (7). As in the -norm or -norm case, because the problem depends on x(t), the implementation of MPC can be performed in two different ways: solving the QP (7) on-line at each time step or, as recently proposed in [], by solving problem (7) off-line for all x(t) within a given range of values, i.e., by considering (7) as a multi-parametric Quadratic Program (mp-qp). The mp-qp problem consists of computing the optimizer U (x) and the value function J (x) for all possible vectors x in a given set X.Once the multi-parametric problem (7) has been solved off line, i.e., the CFTOC solution U (x) of (7) has been found, the explicit MPC law is simply u(t) =[I m... ]U (x). (8) In [] the authors prove the following results about the properties of the solution Theorem Consider the multi-parametric quadratic program (7) and let H. Then the set of feasible parameters X f is convex,the optimizer U (x) :X f R s is continuous and piecewise affine,and the optimal solution J (x) :X f R is continuous,convex and piecewise quadratic.
3 Corollary The MPC control law (8) u(t) = f(x), f : R n R m,defined by the optimization problem () and () is continuous and piecewise affine. Corollary states that by using an mp-qp solver the computational effort for MPC is reduced to a piecewise affine function evaluation.in the next section we propose a method to efficiently evaluate such a piecewise affine function without storing the polyhedral regions of the feasible domain X f. Efficient On-Line Algorithms Let the explicit optimal control law be: u (x) =F i x + G i, x P i, i =,...,N r (9) { where F i R m n, G i R m, and} P i = x R n H i x K i,h i R N i c n,k i R N i c, i =,...,N r is a polyhedral partition of X f.in the following H j i denotes the j-row of the matrix H i.the online implementation of the control law (9) can be simply executed according to the following steps: Algorithm. Measure the current state x. Search for the j-th polyhedron that contains x, i.e. H j x K j. Implement the j-th control law, i.e. u (x) =F j (x)+g j In Algorithm, step (.) is critical and it is the only step whose efficiency can be improved.a trivial implementation of step (.) would consist of searching for the polyhedral region that contains the state x as in the following algorithm Algorithm. i =, notfound=;. while i N r and notfound.. j =, feasible=.. while j Nc i and feasible... if H j i x>kj i then feasible=... else j = j +.. end.. if feasible= then notfound=. end Algorithm requires the storage of all polyhedra P i,i.e., (n +)N C real numbers and in the worst case (the state is contained in the last region of the list) it will give a solution after nn C multiplications, (n )N C sums and N C comparisons, where N C N r i= N c i. By using the properties of the value function we show how Algorithm can be replaced by a more efficient algorithm that avoids storing the polyhedral regions P i, i =,...,N r and is computationally faster. In the following we need to distinguish between optimal control based on LP and optimal control based on QP. fx () 5 f f f P P P P x f Figure : Example for Algorithm in one-dimensional case. For a given point x P (x =5)wehave f (x) =max(f (x),f (x),f (x),f (x)).. Efficient Implementation - LP case From Theorem, the value function J (x) corresponding to the explicit solution of CFTOC () based on LP is convex and piecewise affine: J (x) =T i x + V i, x P i, i =,...,N r. () The value function can be used to avoid storing the polyhedral regions P i.from the equivalence of the representations of piecewise linear convex functions [9], the function J (x) in equation () can be represented alternatively as J (x) =max{t i x + V i,i=,...,n} for x X f () Exploiting the equivalence of () and (), the polyhedral region P j containing the state x can simply be identified by searching for the maximum among the list {T i x + V i, } Nr i= : x P j T jx + V j =max{t i x + V i, i =,...,N r }. () Therefore, instead of searching for the polyhedron j that contains the point x, we store the expression of the piecewise-linear value function and identify the region j by searching for the maximum in the list of numbers composed of the single affine function T i x+v i evaluated at x (see Figure ). Algorithm is then substituted by the following: Algorithm. Compute the vector n R Nr, where n i = T i x + V i. Find j such that n j =max(n) Algorithm does not require the storage of any polyhedra P i, but it only requires the storage of the value function,i.e.,(n +)N r real numbers, and it will give a solution after nn r multiplications, (n )N r sums, and N r comparisons.in Table we compare the complexity of Algorithm against Algorithm in terms of storage demand and number of flops. Table : Complexity of Algorithm and Algorithm Algorithm Algorithm Storage (real numbers) (n +)N C (n +)N r Flops (worst case) nn C nn r
4 . Efficient Implementation - QP case Consider the explicit solution of CFTOC () based on QP.From Theorem, the value function J (x) isconvex and piecewise quadratic.the simple Algorithm described in the previous subsection cannot be applied. Therefore in a new algorithm is proposed. Definition Two polyhedra P i, P j are called neighboring polyhedra if they share a facet. Let {P i } Nr i= be the polyhedral partition obtained by solving the mp-qp (7) and denote by C i = {j P j is a neighbor of P i, j =,...,N r,j i} the list of all neighboring polyhedra of P i. Theorem Let f(x) be a continuous real-valued PWA function f(x) ={f i(x) =A ix + B i x P i, i =,...,N r}, () with A i A j, j C i, i =,...,N r. () Define the list O i (x) = { o i j (x) j C i},where { o i j(x) = + f i (x) f j (x) f i (x) <f j (x) Then O i (x) has the following properties: (i) O i (x) =S i = const, x P i, i =,...,N r. (5) (ii) O i (x) S i, x / P i, i =,...,N r. Proof: Let F = P i P j be the common facet of P i and P j.define the linear function gj i (x) =f i(x) f j (x). (6) From the continuity of f(x) it follows that gj i(x) =, x F.From convexity of P i and P j it follows that gj i(x) is the separating hyperplane between P i and P j. (i) Because gj i (x) = is a separating hyperplane, the function gj i(x) does not change its sign for all x P i, i.e., o i j (x) =si j, x P i with s i j =+orsi j =.The same reasoning can be applied to all neighbors of P i to get S i = {s i j, j C i, x P i }.(ii) x / P i, j C i such that H j i x>kj i.since gi j (x) = is a separating hyperplane, o i j (x) = si j. Equivalently, Theorem states the following property x P i O i (x) =S i (7) where S i uniquely characterizes P i.therefore to check on-line if the polyhedral region i contains the state x it is sufficient to compute the binary vector O i (x) and compare it with S i.list S i = {s i j j C i} is calculated off-line, by comparing the values of f i (x) andf j (x) for j C i in a point x belonging to P i, for instance, the Chebychev center of P i. Figure illustrates the procedure with regions.the list of neighboring regions C i and the list S i can be constructed by simply looking at the figure: C = {}, C = {, }, C = {, }, C = {}, S = { }, S = {, }, S = {, }, S = { }.The point x = is in region and we have O (x) = {, } = S, while O (x) ={, } S, O (x) = fx () 5 f f x f f P P P P Figure : Example for Algorithm in one-dimension. {} S, O (x) ={} S.Obviously, if we are in the wrong region, the checking procedure gives us information about the right search direction(s).the solution can be found by searching in the direction where constraint(s) are violated, i.e., we should check the neighboring region P j for which o i j (x) si j. The overall procedure is composed of two parts:. (off-line) Construction of the PWA function f(x) in () satisfying () and computation of the list of neighbors C i and the list S i,. (on-line) Execution of the following algorithm Algorithm. i =, notfound=;. while notfound.. compute O i (x).. if O i (x) =S i then notfound=.. else i = q, where o i q (x) si q, q C i. end Algorithm does not require the storage of the polyhedra P i, but it only requires the storage of one linear function f i (x) per polyhedron, i.e., (n +)N r real numbers and the list of neighbors C i which requires N C integers.in the worst case Algorithm terminates after nn r multiplications, (n )N r sums and N C comparisons.in Table we compare the complexity of Algorithm against Algorithm in terms of storage demand and number of flops. Remark Note that the computation of O i (x) inalgorithm requires the evaluation of Nc i linear functions, but the overall computation never exceeds N r linear function evaluations.consequently, Algorithm will outperform Algorithm, since typically N C N r. Table : Complexity of Algorithm and Algorithm Algorithm Algorithm Storage (real n.) (n +)N C (n +)N r Flops (worst case) nn C (n )N r + N C The PWA function f(x) in () satisfying () will be referred to as PWA descriptor function.in the rest of the section we propose a way to construct a PWA descriptor function from the PWQ value function.in the following we will assume that the mp-qp (7) is not degenerate (refer to [7] for a detailed discussion on the degeneracy of mp-qp).
5 j qx j( ) qx i( ) (a) i j qx j( ) (b) qx i( ) Figure : Two convex piecewise quadratic functions: (a) not differentiable one, (b) differentiable one. Let J (x) be the convex and piecewise quadratic (CPWQ) value function corresponding to the explicit solution of CFTOC () based on QP: } J (x) = {q i (x) x Q i x + T i x + V i, for x P i, i =,...,N r, (8) We will prove that the J (x) isac () function and we will obtain a PWA descriptor function by differentiating it.before going further we recall the following result [, 7]: Theorem Consider the set of active constraints at the optimum of QP (7) and assume there are no degeneracies: A (x) ={i G i U (x) =W i + S i x}, (9) then. A (x) is constant x P i, i =,...,N r,i.e. A (x) A i x P i. If P i and P j are neighboring polyhedra then A i A j or A j A i. Theorem 5 Consider the value function J (x) in (8). Let P i, P j be two neighboring polyhedra then Q i Q j or Q i Q j and Q i Q j () and Q i Q j iff A i A j () Proof: The proof is given in []. We recall a property of convex piecewise quadratic functions [9]: Proposition Consider the value function J (x) in (8) satisfying () and its quadratic expression q i (x) and q j (x) on two neighboring polyhedra P i, P j with common boundary a x = b. Then there exist constants γ R/{} and b such that q i (x) =q j (x)+(a x b)(γa x b) () Equation () states that the neighboring expressions q i (x) andq j (x) of a CPWQ function satisfying () either intersect on two parallel hyperplanes a x = b and γa x = b if b γb (see Figure (a)), or they are tangent in one hyperplane a x = b if b = γb (see Figure (b)). Theorem 6 Under the hypothesis of Theorem,the value function J (x) in (8) is C (). i Proof: We will prove that b and b in () satisfy b = γb by contradiction.suppose there exists i {,...,N r } and j C i such that b γb.without loss of generality we assume that (i) Q i Q j and (ii) P i is on the side a x b of the common boundary.let F ij be the common facet between P i and P j and F ij its interior. From (i) and from (), γ<andγ =iffq i Q j =. Take x F ij, for sufficiently small ε, the point x = x aε belongs to P i. Let J (ε) =J (x aε) then q i (ε) =q j (ε)+(a aε)(γa aε ( b γb)) () From convexity of J (ε), J (ε) J + (ε)wherej (ε) and J + (ε) are the left and right derivatives of J (ε). This implies q j (ε) q i (ε) whereq j (ε) andq i (ε) arethe derivatives of q j (ε) andq i (ε).from () the latter is true if and only if (a a)( b γb)) γ(a a) ε,that implies b <γbsince γ<andε>. From () q i (ε) <q j (ε) for all ε (, b γb γa a ). Thus there exists x P i with q j (x) <q i (x).this is a contradiction since from Theorem, A i A j. Note that in case of degeneracy the value function J (x) in (8) is not C (), counterexamples are given in [5]. Combining Theorem 6 and (), we obtain q i (x) =q j (x)+γ(a x b) () Theorem 7 Consider the value function J (x) in (8). Let m(x) J (x),be the gradient vector of J (x), i.e., m(x) ={m i (x) q i (x) =Q i x+t i x P i,i=,...,n r },where m : x R n m(x) R n satisfies the following properties (i) m(x) is continuous; (ii) w R n such that for each i and j with P i, P j neighboring polyhedra gj i(x) w (m i (x) m j (x)) for any x which is not an element of their common hyperplane. Proof: (i) Follows easily from Theorem 6.(ii) Consider two neighboring polyhedra P i, P j and their common hyperplane a x = b, a.from equation () we get gj(x) i =γw a(a x b). (5) Since γ,a,anda x b it follows that gj i(x) = if and only if w a = (i.e., vectors w and a are orthogonal).in order to have gj i (x) forall neighboring polyhedra P i and P j, the vector w must not be parallel to any of the common hyperplanes.as the number of common hyperplanes defining the polyhedral partition is finite, such a vector w exists and can be easily computed. It follows from Theorem 7 that an appropriate PWA descriptor function f(x) can be calculated from the gradient of the value function.it is sufficient to find off-line
6 a vector w that is not parallel to any hyperplane of the polyhedral partition.then f i (x) = A i x + B i, where A i =w Q i and B i = w T i. Remark Note that Algorithm applies also to the LP case by using the value function J (x) in () as descriptor function.algorithm is less simple than Algorithm but it requires the evaluation of N r linear functions only in the worst case. Example We compare the performance of Algorithm, Algorithm and Algorithm on an CFTOC problem for the linear system y(t) = u(t), (6) s or its equivalent discrete-time (T s = ) state-space representation x(t +) = x(t)+.5 u(t) [ ] y(t) =.8... x(t) (7) subject to the input constraint u(k) (8) and the output constraint y(k) (9). Explicit MPC based on LP To regulate (7), we design an MPC controller based on the optimization problem () where p =, N =, Q = diag{[5 ]}, R =.8, P =.The explicit solution of the mp-lp problem consist of a polyhedral partition of the state-space consist in 6 regions.in Table we report the comparison between the complexity of Algorithm and Algorithm for this example. The average on-line MPC computation for a set of random points in the state space is 59 flops (Algorithm ), and 88 flops (Algorithm ).The implicit solution with Matlab Optimization Toolbox LP solver takes 559 flops on average. Table : Complexity comparison of Algorithm andalgorithm for example in Section. Algorithm Algorithm Storage (real numbers) Flops (worst case) Explicit MPC based on QP To regulate (7), we design an MPC controller based on the optimization problem () where p =,N =7, Q = diag{[]}, R =., P =.The explicit solution of the mp-qp problem consists of a polyhedral partition of the state-space consist in regions.for this example the choice of w =[] is satisfactory.in Table we report the comparison between Algorithm and Algorithm for this example. The average on-line MPC computation for a set of random points in the state space is flops (Algorithm ), and 75 flops (Algorithm ).The implicit solution with Matlab s QP solver takes 5 flops on average. Table : Complexity comparison of Algorithm andalgorithm for example in Section. Algorithm Algorithm Storage (real numbers) Flops (worst case) Conclusion By exploiting the properties of the value function of MPC, we presented two algorithms that significantly improve the efficiency of the on-line explicit MPC (both based on LP and QP) in terms of storage demands and computational complexity.the following improvements are achieved.there is no need to store the polyhedral partition of the state space for computing the optimal control law..in the worst case, the optimal control law is computed after the evaluation of one linear function per polyhedron. References [] A. Bemporad, F. Borrelli, and M. Morari. Explicit solution of constrained / -norm model predictive control. In Proc. 9th IEEE Conf. on Decision and Control, December. [] A. Bemporad, M. Morari, V. Dua, and E. N. Pistikopoulos. The explicit linear quadratic regulator for constrained systems. Automatica, toappear. [] F. Borrelli, M. Baotic, A. Bemporad, and M. Morari. Efficient on-line computation of explicit model predictive control. Technical Report AUT-5, Automatic Control Laboratory, ETH Zurich, Switzerland, August. [] F. Borrelli, A. Bemporad, and M. Morari. A geometric algorithm for multi-parametric linear programming. Technical Report AUT-6, Automatic Control Laboratory, ETH Zurich, Switzerland, February. [5] T. Gal. Postoptimal Analyses, Parametric Programming, and Related Topics. de Gruyter, Berlin, nd edition, 995. [6] D.Q. Mayne, J.B. Rawlings, C.V. Rao, and P.O.M. Scokaert. Constrained model predictive control: Stability and optimality. Automatica, 6(6):789 8, June. [7] P. Tøndel, T.A. Johansen, and A. Bemporad. An algorithm for multi-parametric quadratic programming and explicit MPC solutions. In Proc. th IEEE Conf. on Decision and Control, December. [8] S.J. Qin and T.A. Badgwell. An overview of industrial model predictive control technology. In Chemical Process Control -V, volume 9, no. 6, pages 56. AIChE Symposium Series - American Institute of Chemical Engineers, 997. [9] M. Schechter. Polyhedral functions and multiparametric linear programming. Journal of Optimization Theory and Applications, 5():69 8, May 987.
Efficient 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 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 informationAn 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 informationComputation of the Constrained Infinite Time Linear Quadratic Optimal Control Problem
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
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 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 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 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 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 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 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 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 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 informationEvaluation of Piecewise Affine Control via Binary Search Tree
Evaluation of Piecewise Affine Control via Binary Search Tree P. Tøndel a T. A. Johansen a A. Bemporad b a Department of Engineering Cybernetics, Norwegian University of Science and Technology, N-7491
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 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 informationTechnical Notes and Correspondence
1600 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 8, NO. 9, SEPTEMBER 2003 Technical Notes and Correspondence Min Max Control of Constrained Uncertain Discrete-Time Linear Systems Alberto Bemporad, Francesco
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 informationComputation of Voronoi Diagrams and Delaunay Triangulation via Parametric Linear Programming
Computation of Voronoi Diagrams and Delaunay Triangulation via Parametric Linear Programming Saša V. Raković, Pascal Grieder and Colin Jones September 14, 2004 Abstract This note illustrates how Voronoi
More informationReverse Search for Parametric Linear Programming
Reverse Search for Parametric Linear Programg Colin N. Jones Automatic Control Laboratory Swiss Federal Institute of Technology Physikstrasse 3, CH-8092 Zurich Switzerland colin.jones@cantab.net Jan M.
More informationImproving Reliability of Partition Computation in Explicit MPC with MPT Toolbox
Improving Reliability of Partition Computation in Explicit MPC with MPT Toolbox Samo Gerkšič Dept. of Systems and Control, Jozef Stefan Institute, Jamova 39, Ljubljana, Slovenia (e-mail: samo.gerksic@ijs.si).
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 informationAn algorithm for multi-parametric quadratic programming and explicit MPC solutions
Available online at www.sciencedirect.com Automatica 39 (2003) 489 497 www.elsevier.com/locate/automatica Brief Paper An algorithm for multi-parametric quadratic programming and explicit MPC solutions
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 informationEfficient Mode Enumeration of Compositional Hybrid Models
Efficient Mode Enumeration of Compositional Hybrid Models F. D. Torrisi, T. Geyer, M. Morari torrisi geyer morari@aut.ee.ethz.ch, http://control.ethz.ch/~hybrid Automatic Control Laboratory Swiss Federal
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 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 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 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 informationA Survey on Explicit Model Predictive Control
A Survey on Explicit Model Predictive Control Alessandro Alessio and Alberto Bemporad Abstract. Explicit model predictive control (MPC) addresses the problem of removing one of the main drawbacks of MPC,
More informationLecture 3. Corner Polyhedron, Intersection Cuts, Maximal Lattice-Free Convex Sets. Tepper School of Business Carnegie Mellon University, Pittsburgh
Lecture 3 Corner Polyhedron, Intersection Cuts, Maximal Lattice-Free Convex Sets Gérard Cornuéjols Tepper School of Business Carnegie Mellon University, Pittsburgh January 2016 Mixed Integer Linear Programming
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 informationEmbedded Optimization for Mixed Logic Dynamical Systems
Embedded Optimization for Mixed Logic Dynamical Systems Alexander Domahidi Joint work with Damian Frick and Manfred Morari EMBOPT Workshop IMT Lucca, Italy September 8, 2014 Application: Optimal Traffic
More informationShiqian Ma, MAT-258A: Numerical Optimization 1. Chapter 2. Convex Optimization
Shiqian Ma, MAT-258A: Numerical Optimization 1 Chapter 2 Convex Optimization Shiqian Ma, MAT-258A: Numerical Optimization 2 2.1. Convex Optimization General optimization problem: min f 0 (x) s.t., f i
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 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 informationReal-time MPC Stability through Robust MPC design
Real-time MPC Stability through Robust MPC design Melanie N. Zeilinger, Colin N. Jones, Davide M. Raimondo and Manfred Morari Automatic Control Laboratory, ETH Zurich, Physikstrasse 3, ETL I 28, CH 8092
More informationLinear programming and duality theory
Linear programming and duality theory Complements of Operations Research Giovanni Righini Linear Programming (LP) A linear program is defined by linear constraints, a linear objective function. Its variables
More informationDiscrete Optimization. Lecture Notes 2
Discrete Optimization. Lecture Notes 2 Disjunctive Constraints Defining variables and formulating linear constraints can be straightforward or more sophisticated, depending on the problem structure. The
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 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 informationApproximate nonlinear explicit MPC based on reachability analysis
Approximate nonlinear explicit MPC based on reachability analysis D.M. Raimondo 1, M. Schulze Darup 2, M. Mönnigmann 2 Università degli studi di Pavia Ruhr-Universität Bochum The system Introduction: Nonlinear
More informationComputation and Approximation of Piecewise Affine Control Laws
Computation and Approximation of Piecewise Affine Control Laws via Binary Search Trees P Tøndel TA Johansen and A Bemporad Abstract We present an algorithm for generating a binary search tree that allows
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 informationPerformance Driven Reachability Analysis for Optimal Scheduling and Control of Hybrid Systems
Performance Driven Reachability Analysis for Optimal Scheduling and Control of Hybrid Systems A. Bemporad,L.Giovanardi,F.D.Torrisi Automatic Control Laboratory, Swiss Federal Institute of Technology -
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 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 informationIntroduction to Modern Control Systems
Introduction to Modern Control Systems Convex Optimization, Duality and Linear Matrix Inequalities Kostas Margellos University of Oxford AIMS CDT 2016-17 Introduction to Modern Control Systems November
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 informationMulti-Parametric Toolbox 3.0
Multi-Parametric Toolbox 3.0 Martin Herceg Manfred Morari Michal Kvasnica Colin N. Jones ETH Zürich STU Bratislava EPFL Lausanne What is MPT? Matlab toolbox for application of explicit MPC high-speed implementation
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 informationPredictive control for hybrid systems. Implications of polyhedral pre-computations
Nonlinear Analysis: Hybrid Systems ( ) www.elsevier.com/locate/nahs Predictive control for hybrid systems. Implications of polyhedral pre-computations Sorin Olaru a,, Didier Dumur a, Jean Thomas b, Marius
More informationLecture 6: Faces, Facets
IE 511: Integer Programming, Spring 2019 31 Jan, 2019 Lecturer: Karthik Chandrasekaran Lecture 6: Faces, Facets Scribe: Setareh Taki Disclaimer: These notes have not been subjected to the usual scrutiny
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 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 informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,900 116,000 120M Open access books available International authors and editors Downloads Our
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 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 informationChapter II. Linear Programming
1 Chapter II Linear Programming 1. Introduction 2. Simplex Method 3. Duality Theory 4. Optimality Conditions 5. Applications (QP & SLP) 6. Sensitivity Analysis 7. Interior Point Methods 1 INTRODUCTION
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 informationSimplification of Explicit MPC Feedback Laws via Separation Functions
Simplification of Explicit MPC Feedback Laws via Separation Functions Michal Kvasnica,1, Ivana Rauová, and Miroslav Fikar Institute of Automation, Information Engineering and Mathematics, Slovak University
More informationLP Geometry: outline. A general LP. minimize x c T x s.t. a T i. x b i, i 2 M 1 a T i x = b i, i 2 M 3 x j 0, j 2 N 1. where
LP Geometry: outline I Polyhedra I Extreme points, vertices, basic feasible solutions I Degeneracy I Existence of extreme points I Optimality of extreme points IOE 610: LP II, Fall 2013 Geometry of Linear
More informationChapter 9. Linear algebra applications in geometry
Chapter 9. Linear algebra applications in geometry C.O.S. Sorzano Biomedical Engineering August 25, 2013 9. Linear algebra applications in geometry August 25, 2013 1 / 73 Outline 9 Linear algebra applications
More informationA robust optimization based approach to the general solution of mp-milp problems
21 st European Symposium on Computer Aided Process Engineering ESCAPE 21 E.N. Pistikopoulos, M.C. Georgiadis and A. Kokossis (Editors) 2011 Elsevier B.V. All rights reserved. A robust optimization based
More informationExplicit model predictive control of gas-liquid. separation plant via orthogonal search tree. partitioning Trondheim, Norway
Explicit model predictive control of gas-liquid separation plant via orthogonal search tree partitioning Alexandra Grancharova,,, Tor A. Johansen, 3, Juš Kocijan Department of Engineering Cybernetics,
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 informationModel Predictive Control Design: New Trends and Tools
Proceedings of the 45th IEEE Conference on Decision & Control Manchester Grand Hyatt Hotel San Diego, CA, USA, December 13-15, 2006 Model Predictive Control Design: New Trends and Tools Alberto Bemporad
More informationConvex Sets (cont.) Convex Functions
Convex Sets (cont.) Convex Functions Optimization - 10725 Carlos Guestrin Carnegie Mellon University February 27 th, 2008 1 Definitions of convex sets Convex v. Non-convex sets Line segment definition:
More informationEvaluation ofpiecewise ane control via binary search tree
Available online at www.sciencedirect.com Automatica 39 (2003) 945 950 www.elsevier.com/locate/automatica Technical Communique Evaluation ofpiecewise ane control via binary search tree P. THndel a;, T.A.
More informationExplicit MPC in Mechatronics Industry:
European Control Conference, July 8 th, 23 Zurich, CH MITSUBISHI ELECTRIC RESEARCH LABORATORIES Cambridge, Massachusetts Explicit MPC in Mechatronics Industry: Technology Transfer Potential and Limitations
More informationOptimization Methods. Final Examination. 1. There are 5 problems each w i t h 20 p o i n ts for a maximum of 100 points.
5.93 Optimization Methods Final Examination Instructions:. There are 5 problems each w i t h 2 p o i n ts for a maximum of points. 2. You are allowed to use class notes, your homeworks, solutions to homework
More informationPiecewise Quadratic Optimal Control
EECE 571M/491M, Spring 2007 Lecture 15 Piecewise Quadratic Optimal Control Meeko Oishi, Ph.D. Electrical and Computer Engineering University of British Columbia, BC http://www.ece.ubc.ca/~elec571m.html
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 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 informationMax Min Control Problems for Constrained Discrete Time Systems
Proceedings of the 47th IEEE Conference on Decision and Control Cancun, Mexico, Dec. 9-11, 2008 Max Min Control Problems for Constrained Discrete Time Systems Saša V. Raković, Miroslav Barić and Manfred
More informationModeling and Control of Hybrid Systems
Modeling and Control of Hybrid Systems Alberto Bemporad Dept. of Information Engineering University of Siena, Italy bemporad@unisi.it MTNS 2004 Leuven, July 8 COHES Group Control and Optimization of Hybrid
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 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 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 informationConvexization in Markov Chain Monte Carlo
in Markov Chain Monte Carlo 1 IBM T. J. Watson Yorktown Heights, NY 2 Department of Aerospace Engineering Technion, Israel August 23, 2011 Problem Statement MCMC processes in general are governed by non
More informationGlobal Minimization via Piecewise-Linear Underestimation
Journal of Global Optimization,, 1 9 (2004) c 2004 Kluwer Academic Publishers, Boston. Manufactured in The Netherlands. Global Minimization via Piecewise-Linear Underestimation O. L. MANGASARIAN olvi@cs.wisc.edu
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 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 informationLecture 4: Convexity
10-725: Convex Optimization Fall 2013 Lecture 4: Convexity Lecturer: Barnabás Póczos Scribes: Jessica Chemali, David Fouhey, Yuxiong Wang Note: LaTeX template courtesy of UC Berkeley EECS dept. Disclaimer:
More informationDiscrete Optimization 2010 Lecture 5 Min-Cost Flows & Total Unimodularity
Discrete Optimization 2010 Lecture 5 Min-Cost Flows & Total Unimodularity Marc Uetz University of Twente m.uetz@utwente.nl Lecture 5: sheet 1 / 26 Marc Uetz Discrete Optimization Outline 1 Min-Cost Flows
More informationSHIP heading control, also known as course keeping, has
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 20, NO. 1, JANUARY 2012 257 Disturbance Compensating Model Predictive Control With Application to Ship Heading Control Zhen Li, Member, IEEE, and Jing
More informationCS 372: Computational Geometry Lecture 10 Linear Programming in Fixed Dimension
CS 372: Computational Geometry Lecture 10 Linear Programming in Fixed Dimension Antoine Vigneron King Abdullah University of Science and Technology November 7, 2012 Antoine Vigneron (KAUST) CS 372 Lecture
More informationCS 473: Algorithms. Ruta Mehta. Spring University of Illinois, Urbana-Champaign. Ruta (UIUC) CS473 1 Spring / 36
CS 473: Algorithms Ruta Mehta University of Illinois, Urbana-Champaign Spring 2018 Ruta (UIUC) CS473 1 Spring 2018 1 / 36 CS 473: Algorithms, Spring 2018 LP Duality Lecture 20 April 3, 2018 Some of the
More informationOptimality certificates for convex minimization and Helly numbers
Optimality certificates for convex minimization and Helly numbers Amitabh Basu Michele Conforti Gérard Cornuéjols Robert Weismantel Stefan Weltge May 10, 2017 Abstract We consider the problem of minimizing
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 informationSlovak University of Technology in Bratislava Institute of Information Engineering, Automation, and Mathematics PROCEEDINGS
Slovak University of Technology in Bratislava Institute of Information Engineering, Automation, and Mathematics PROCEEDINGS of the 18 th International Conference on Process Control Hotel Titris, Tatranská
More informationSequential Coordinate-wise Algorithm for Non-negative Least Squares Problem
CENTER FOR MACHINE PERCEPTION CZECH TECHNICAL UNIVERSITY Sequential Coordinate-wise Algorithm for Non-negative Least Squares Problem Woring document of the EU project COSPAL IST-004176 Vojtěch Franc, Miro
More informationqpoases - Online Active Set Strategy for Fast Linear MPC
qpoases - Online Active Set Strategy for Fast Linear MPC Moritz Diehl, Hans Joachim Ferreau, Lieboud Vanden Broeck, Jan Swevers Dept. ESAT and Center of Excellence for Optimization in Engineering OPTEC
More informationOptimality certificates for convex minimization and Helly numbers
Optimality certificates for convex minimization and Helly numbers Amitabh Basu Michele Conforti Gérard Cornuéjols Robert Weismantel Stefan Weltge October 20, 2016 Abstract We consider the problem of minimizing
More informationCMU-Q Lecture 9: Optimization II: Constrained,Unconstrained Optimization Convex optimization. Teacher: Gianni A. Di Caro
CMU-Q 15-381 Lecture 9: Optimization II: Constrained,Unconstrained Optimization Convex optimization Teacher: Gianni A. Di Caro GLOBAL FUNCTION OPTIMIZATION Find the global maximum of the function f x (and
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 informationCS 473: Algorithms. Ruta Mehta. Spring University of Illinois, Urbana-Champaign. Ruta (UIUC) CS473 1 Spring / 29
CS 473: Algorithms Ruta Mehta University of Illinois, Urbana-Champaign Spring 2018 Ruta (UIUC) CS473 1 Spring 2018 1 / 29 CS 473: Algorithms, Spring 2018 Simplex and LP Duality Lecture 19 March 29, 2018
More informationA Learning Algorithm for Piecewise Linear Regression
A Learning Algorithm for Piecewise Linear Regression Giancarlo Ferrari-Trecate 1, arco uselli 2, Diego Liberati 3, anfred orari 1 1 nstitute für Automatik, ETHZ - ETL CH 8092 Zürich, Switzerland 2 stituto
More information11 Linear Programming
11 Linear Programming 11.1 Definition and Importance The final topic in this course is Linear Programming. We say that a problem is an instance of linear programming when it can be effectively expressed
More information