Theory in Model Predictive Control :" Constraint Satisfaction and Stability!
|
|
- Marjory Booth
- 5 years ago
- Views:
Transcription
1 Theory in Model Predictive Control :" Constraint Satisfaction and Stability Colin Jones, Melanie Zeilinger Automatic Control Laboratory, EPFL
2 Example: Cessna Citation Aircraft Linearized continuous-time model: (at altitude of 5m and a speed of m/sec) ẋ = x y = x 1 Input: elevator angle States: x 1 : angle of attack, x 2 : pitch angle, x 3 : pitch rate, x 4 : altitude Outputs: pitch angle and altitude Constraints: elevator angle ±.262rad (±15º), elevator rate ±.524rad (±6º) pitch angle ±.349 (±3º) horizon Open-loop response is unstable (open-loop poles:,, ±2.29i ).3 17 u V Angle of attack Pitch angle [J. Maciejowski, Predictive Control with constraints, 22] "#$%#&'('$)#*$+( +,$*#$*- 2
3 LQR and Linear MPC with quadratic cost Quadratic performance measure Linear system dynamics x + = Ax + Bu y = Cx + Du Linear constraints on inputs and states LQR: J (x) = min x,u i= x i+1 x x T i Qx i + u T i Ru i = Ax i + Bu i = x MPC: J (x) = min x,u N 1 Q = Q T,R = R T MPC problem can be translated into a quadratic program (QP) i= x i+1 x x T i Cx i + Du i Qx i + u T i Ru i = Ax i + Bu i = x b "#$%#&'('$)#*$+( +,$*#$*- 3
4 Linear MPC with linear costs Linear performance measure (e.g. economically motivated cost, errors) Linear system dynamics x + = Ax + Bu y = Cx + Du Linear constraints on states and inputs Resulting MPC problem: J (x) = min x,u Optimization problem can be translated into a linear program (LP) for p =1/. N 1 i= Qx i p + Ru i p x i+1 = Ax i + Bu i Cx i + Du i x b = x "#$%#&'('$)#*$+( +,$*#$*- 4
5 Example: Cessna Citation Aircraft" LQR with saturation Altitude x 4 (m) Linear quadratic regulator with saturated inputs: At time t = the plane is flying with a deviation of 1m of the desired altitude, i.e. x ()=[; ; ; 1] Elevator angle u (rad) Time (sec).5 Input is saturated Time (sec) Pitch angle x 2 (rad) Problem parameters: Sampling time T s =.25sec, Q=I, R=1 Closed-loop system is unstable Applying LQR control and saturating the controller can lead to instability "#$%#&'('$)#*$+( +,$*#$*- 5
6 Example: Cessna Citation 5 aircraft" MPC with bound constraints on inputs MPC controller with input constraints u i ".262 Problem parameters: Sampling time T s =.25sec, Q=I, R=1, N=1 Altitude x 4 (m) Elevator angle u (rad) Time (sec) Time (sec) Pitch angle x 2 (rad) The MPC controller uses the knowledge that the elevator will saturate, but it does not consider the rate constraints. System does not converge to desired steady-state but to a limit cycle "#$%#&'('$)#*$+( +,$*#$*- 6
7 Example: Cessna Citation Aircraft" MPC with all input constraints MPC controller with input constraints u i ".262 and rate constraints u i.349 approximated by u k u k 1.349T s,u 1 = u prev Altitude x 4 (m) Elevator angle u (rad) Time (sec) Time (sec).2.2 Pitch angle x 2 (rad) Problem parameters: Sampling time T s =.25sec, Q=I, R=1, N=1 The MPC controller considers all constraints on the actuator Closed-loop system is stable Efficient use of the control authority "#$%#&'('$)#*$+( +,$*#$*- 7
8 Example: Cessna Citation Aircraft" Inclusion of state constraints MPC controller with input constraints u i ".262 and rate constraints u i.349 approximated by u k u k 1.349T s,u 1 = u prev Problem parameters: Sampling time T s =.25sec, Q=I, R=1, N=1 Altitude x 4 (m) Elevator angle u (rad) Time (sec).5 Pitch angle #-.9, i.e Time (sec).5.5 Pitch angle x 2 (rad) Increase step: At time t = the plane is flying with a deviation of 1m of the desired altitude, i.e. x ()=[; ; ; 1] Pitch angle too large during transient "#$%#&'('$)#*$+( +,$*#$*- 8
9 Example: Cessna Citation Aircraft" Inclusion of state constraints MPC controller with input constraints u i ".262 and rate constraints u i.349 approximated by u k u k 1.349T s,u 1 = u prev Problem parameters: Sampling time T s =.25sec, Q=I, R=1, N=1 Altitude x 4 (m) Elevator angle u (rad) Time (sec).5 Constraint on pitch angle active Time (sec).4.2 Pitch angle x 2 (rad) Add state constraints for passenger comfort: x 2 ".349 "#$%#&'('$)#*$+( +,$*#$*- 9
10 Example: Cessna Citation Aircraft" Shorter horizon causes loss of stability MPC controller with input constraints u i ".262 and rate constraints u i.349 approximated by u k u k 1.349T s,u 1 = u prev Problem parameters: Sampling time T s =.25sec, Q=I, R=1, N=4 Altitude x 4 (m) Elevator angle u (rad) Time (sec) Time (sec).5 Pitch angle x 2 (rad) Decrease in the prediction horizon causes loss of the stability properties This part of the workshop: How can constraint satisfaction and stability in MPC be guaranteed? "#$%#&'('$)#*$+( +,$*#$*- 1
11 Outline Motivating Example: Cessna Citation Aircraft Constraint satisfaction and stability in MPC Main idea Problem setup for the linear quadratic case How to prove constraint satisfaction and stability in MPC Implementation Theory extensions "#$%#&'('$)#*$+( +,$*#$*- 11
12 Loss of feasibility and stability guarantees What can go wrong with standard MPC approach? No feasibility guarantee, i.e. the MPC problem may not have a solution No stability guarantee, i.e. trajectories may not converge to the origin Definition: Feasible set X N x N X N := {x [u,...,u N 1 ] Cu i + Dx i b, i =1,...,N} min x,u N 1 i= x T i Qx i + u T i Ru i x i+1 = Ax i + Bu i x Cx i + Du i = x b "#$%#&'('$)#*$+( +,$*#$*- 12
13 Example: Loss of feasibility Consider the double integrator: subject to the input constraints.5 u and the state constraints x We choose N =3,Q=,R = 1. Set of initial Time step 1: x ()=[-4 ; 4], u * (x)=-.5 Time step 2: x (1)=[ ; 3], u * (x)=-.5 Time step 3: x (2)=[ 3 ; 2] MPC problem infeasible "#$%#&'('$)#*$+( +,$*#$*- x + = x 2 x u 1 x() feasible states (feasible set) x(1) 5 5 x 1 x(3) 13
14 Example: Loss of stability Consider the unstable system: subject to the input constraints and the state constraints 1 We choose Q = 1 and investigate the stability properties for different horizons N and weights R by solving the finite horizon MPC problem in a receding horizon fashion x + = x u 1 x u "#$%#&'('$)#*$+( +,$*#$*- 14
15 Example: Loss of stability R=1, N=2 R=2, N=3 R=1, N=4 * Initial points leading to trajectories that converge to the origin " Intitial points that diverge 1 x x x 2 x x 1 x 1 Parameters have complex effect on closed-loop trajectory "#$%#&'('$)#*$+( +,$*#$*- 15
16 Feasibility and stability in MPC Main Idea Main idea: Introduce terminal cost and constraints to explicitly ensure stability and feasibility: How to choose P and X f? J (x) = min x,u "#$%#&'('$)#*$+( +,$*#$*- N 1 i= x T i Qx i + u T i Ru i + x T N Px N x i+1 = Ax i + Bu i Cx i + Du i x N x b X f = x Terminal constraint Terminal cost X f t x N x =x 16
17 How to choose the terminal set and cost" Main Idea Problems originate from the use of a short sight strategy Finite horizon causes deviation between the open-loop prediction and the closed-loop system: x 2 5 Closed-loop trajectories x 1 Open-loop predictions x 2 Set of feasible initial states for open-loop prediction Set of initial states leading to feasible closedloop trajectories Ideally we would solve the MPC problem with an infinite horizon, but that is computationally intractable Design finite horizon problem such that it approximates the infinite horizon 5 x 1 "#$%#&'('$)#*$+( +,$*#$*- 17
18 How to choose the terminal cost We can split the infinite horizon problem into two subproblems: Up to time k =N, where the constraints may be active J (x) = min x,u N 1 i= x T i Qx i + u T i Ru i x i+1 = Ax i + Bu i Cx i + Du i x b = x Bound the tail of the infinite horizon cost from N to using the LQR control law u=k LQR x x N T P x N is the corresponding infinite horizon cost P is the solution of the discrete-time algebraic Riccati equation Choice of N such that constraint satisfaction is guaranteed? For k >N, where there are no constraints active + min xi T Qx i + ui T Ru i x,u i=n x i+1 = Ax i + Bu i, + x N T P x N Unconstrained LQR starting from state x N "#$%#&'('$)#*$+( +,$*#$*- 18
19 How to choose the terminal set Terminal constraint provides a sufficient condition for constraint satisfaction : N 1 J (x) = min x Infinite horizon cost i T Qx i + ui T Ru i + xn T Px N x,u starting from x i= N x X f i+1 = Ax i + Bu i Cx t i + Du i b x N X f x =x x = x All input and state constraints are satisfied for the closed-loop system using the LQR control law for x X f Terminal set is often defined by linear or quadratic constraints The bound holds in the terminal set and is used as a terminal cost The terminal set defines the terminal constraint In the following: Show that this problem setup provides feasibility and stability x N "#$%#&'('$)#*$+( +,$*#$*- 19
20 Outline Motivating Example: Cessna Citation Aircraft Constraint satisfaction and stability in MPC Main idea Problem setup for the linear quadratic case How to prove constraint satisfaction and stability in MPC Implementation Theory extensions "#$%#&'('$)#*$+( +,$*#$*- 2
21 Formalize goals: " Definition of feasibility and stability Goal 1: Feasibility at all times Definition: Recursive feasibility The MPC problem is called recursively feasible, if for all feasible initial states feasibility is guaranteed at every state along the closed-loop trajectory. Goal 2: Stability Definition: Lyapunov stability x(k + 1) = Ax(k) +Bκ(x(k)) = f κ (x(k)) X > δ() > x() X % x() δ() x(k) < k N. $ x() "#$%#&'('$)#*$+( +,$*#$*- 21
22 Employed concept for the analysis of feasibility: " Invariant sets Definition: Invariant set O x(k + 1) = f κ (x(k)) x(k) O f κ (x(k)) O, k N. The positively invariant set that contains every closed positively invariant set is called the maximal positively invariant set O. Infeasible after two steps Infeasible after one step O Invariant Recursively feasible "#$%#&'('$)#*$+( +,$*#$*- 22
23 Analysis of Lyapunov stability Lyapunov stability will be analyzed using Lyapunov s direct method: Definition: Lyapunov function X x(k + 1) = f κ (x(k)) V : X R 1 + X x X V (x) > x =,V() =, V (x(k + 1)) V (x(k)) Theorem: (e.g. [Vidyasagar, 1993]) X X 1 For simplicity it is assumed that V(x) is continuous. This assumption can be relaxed by requiring an additional state dependent upper bound on V (x), see e.g. [Rawlings & Mayne, 29] "#$%#&'('$)#*$+( +,$*#$*- 23
24 How to prove feasibility and stability of MPC Main steps: Prove recursive feasibility by showing the existence of a feasible control sequence at all time instants when starting from a feasible initial point NOTE: Recursive feasibility does not imply stability of the closed-loop system Prove stability by showing that the optimal cost function is a Lyapunov function We will discuss two main cases in the following: Terminal constraint at zero: x N = Terminal constraint in some (convex) set: x N & X f For simplicity, we use the more general notation: J (x) = min x,u N 1 i= l(x i,u i )+V f (x N ) stage cost terminal cost ( In the quadratic case: l(x i,u i )=x T Qx i + ui T Ru i,v f (x N )=xn T Px N ) i "#$%#&'('$)#*$+( +,$*#$*- 24
25 Stability of MPC Zero terminal state constraint Terminal constraint x N = : Assume feasibility of x and let [u *, u 1*,...,u N * -1 ] be the optimal control sequence computed at x At x (k + 1) the control sequence [u 1*, u 2*,...,u N * -1, ] is feasible (apply control input to stay at the origin) Recursive feasibility The associated cost function value is given by: We obtain for the optimal solution x 1 x = x x 5 = x 4 x 2 x 3 "#$%#&'('$)#*$+( +,$*#$*- 25
26 Stability of MPC Zero terminal state constraint Terminal constraint x N = : Assume feasibility of x and let [u *, u 1*,...,u N * -1 ] be the optimal control sequence computed at x At x + the control sequence [u 1*, u 2*,...,u N * -1, ] is feasible (apply control input to stay at the origin) Recursive feasibility The associated cost function value is given by: We obtain for the optimal solution x = x x x 5 = x 4 x 1 = x + 2 x 3 "#$%#&'('$)#*$+( +,$*#$*- 26
27 Stability of MPC Zero terminal state constraint Terminal constraint x N = : Assume feasibility of x and let [u *, u 1*,...,u N * -1 ] be the optimal control sequence computed at x At x + the control sequence [u 1*, u 2*,...,u N * -1, ] is feasible (apply control input to stay at the origin) Recursive feasibility The associated cost function value is given by: We obtain for the optimal solution J (x) J(x) J (x + ) J (x) J(x + ) J (x) l(x,u ) J * (x) is a Lyapunov function (Lyapunov) Stability "#$%#&'('$)#*$+( +,$*#$*- J(x + )=J (x) l(x,u )+l(, ) Subtract cost at stage Add cost for staying at x = x x x 5 = x 4 x 1 = x + 2 x 3 27
28 Extension to more general terminal sets Problem: The terminal constraint x N = reduces the feasible set Goal: " Use convex set for X f that contains the origin in its interior Example: 3 Feasible set for x N & X f 2 Feasible set for x N = 1 x "#$%#&'('$)#*$+( +,$*#$*- x 1 How can the proof be generalized to the constraint x N X f? X f Double Integrator: x + = x 5 N =5,Q= x u.5 u.5,r = 1 28
29 Stability of MPC Main result Consider that the following standard assumptions hold: The stage cost is a positive definite function, i.e. it is strictly positive and only zero at the origin The terminal set is invariant under the local control law ' f (x): x + = Ax + Bκ f (x) X f x X f All state and input constraints are satisfied in X f : X f X, κ f (x) U x X f The terminal cost is a continuous Lyapunov function in the terminal set X f : V f (x + ) V f (x) l(x,κ f (x)) x X f The closed-loop system under the MPC control law is stable in the feasible set X N. "#$%#&'('$)#*$+( +,$*#$*- 29
30 Stability of MPC Outline of the proof Assume feasibility of x and let [u *, u 1*,...,u N * -1 ] be the optimal control sequence computed at x At x (k + 1) the control sequence [u 1*, u 2*,...,' f (x N* )] is feasible: x N is in X f ' f (x N* ) is feasible and x N+ =Ax N* +B ' f (x N* ) in X f Terminal constraint provides recursive feasibility The associated cost function value is given by: x = x x 5 X f x 4 x 1 x 2 x 3 V f (x) is a Lyapunov function: " We obtain for the optimal solution : J * (x) is a Lyapunov function (Lyapunov) Stability "#$%#&'('$)#*$+( +,$*#$*- 3
31 Stability of MPC Outline of the proof Assume feasibility of x and let [u *, u 1*,...,u N * -1 ] be the optimal control sequence computed at x At x + the control sequence [u 1*, u 2*,...,' f (x N* )] is feasible: x N is in X f ' f (x N* ) is feasible and x N+ =Ax N* +B ' f (x N* ) in X f Terminal constraint provides recursive feasibility The associated cost function value is given by: x = x x 6 x 5 Xf x 4 x 1 = x + x 2 x 3 V f (x) is a Lyapunov function: " We obtain for the optimal solution : J * (x) is a Lyapunov function (Lyapunov) Stability "#$%#&'('$)#*$+( +,$*#$*- 31
32 Stability of MPC Outline of the proof Assume feasibility of x and let [u *, u 1*,...,u N * -1 ] be the optimal control sequence computed at x At x + the control sequence [u 1*, u 2*,...,' f (x N* )] is feasible: x N is in X f ' f (x N* ) is feasible and x N+ =Ax N* +B ' f (x N* ) in X f Terminal constraint provides recursive feasibility The associated cost function value is given by: x = x x 6 x 5 X f x 4 x 1 = x + x 2 x 3 J(x) =J (x) l(x,u )+V f ( x N+1 ) V f (x N )+l(x N, κ f (x N )) V f (x) is a Lyapunov function: " We obtain for the optimal solution J (x) J(x) : J (x + ) J (x) J(x + ) J (x) l(x,u ) J * (x) is a Lyapunov function (Lyapunov) Stability "#$%#&'('$)#*$+( +,$*#$*- 32
33 Stability of MPC Remarks The terminal set X f and the terminal cost ensure recursive feasibility and stability of the closed-loop system. But: the terminal constraint usually reduces the region of attraction. If the open-loop system is stable, X f can be chosen as the positively invariant set of the system under zero control input, which is feasible. Often no terminal set X f, but N is required to be sufficiently large to ensure recursive feasibility Applied in practice, makes MPC work without terminal constraint Determination of a sufficiently long horizon difficult "#$%#&'('$)#*$+( +,$*#$*- 33
34 Proof of asymptotic stability Definition: Asymptotic stability Extension of Lyapunov s direct method: (see e.g. [Vidyasagar, 1993]) If the continuous Lyapunov function additionally satisfies V (x(k + 1)) V (x(k)) < x then the closed loop system converges to the origin and is hence asymptotically stable. Recall: Decrease of the optimal MPC cost was given by X x(k + 1) = f κ (x(k)) X % X lim k x(k) = x() X J (x(k + 1)) J (x(k)) l(x(k),u ) where the stage cost was assumed to be positive and only at. The closed-loop system under the MPC control law is asymptotically stable. $ x() "#$%#&'('$)#*$+( +,$*#$*- 34
35 Extension to nonlinear MPC Consider the nonlinear system dynamics: Nonlinear MPC problem: J (x) = min x,u N 1 Presented assumptions on the terminal set and cost did not rely on linearity Lyapunov stability is a general framework to analyze stability of nonlinear dynamic systems Results can be directly extended to nonlinear systems. i= l(x i,u i )+V f (x N ) x i+1 = f (x i,u i ) Cx i + Du i b x N X f, x = x, x + = f (x,u) "#$%#&'('$)#*$+( +,$*#$*- 35
36 Outline Motivating Example: Cessna Citation Aircraft Constraint satisfaction and stability in MPC Main idea Problem setup for the linear quadratic case How to prove constraint satisfaction and stability in MPC Implementation Theory extensions "#$%#&'('$)#*$+( +,$*#$*- 36
37 Multi-Parametric Toolbox Modeling Optimal control Analysis Computational geometry MPT Toolbox: M. Kvasnica, P. Grieder and M. Baotic Real-Time Model Predictive Control via Multi-Parametric Programming: Theory and Tools Michal Kvasnica ISBN: "#$%#&'('$)#*$+( +,$*#$*- 37
38 Implementation using Matlab and MPT 1. Compute terminal weight P and LQR control law K by solving the discrete time Riccati equation Matlab [K,P]=dlqr(A,B,Q,R) 2. Compute terminal set X f : Ellipsoidal invariant set of the form Can be written in the form of a Linear Matrix inequality (LMI) [Boyd et al., LMIs in System and Control Theory, 1994] Polytopic invariant set of the form X f := {x x T P E x 1} X f := {x Hx K} Double Integrator: x + = x 5 N =5,Q= x u.5 u.5,r = 1 X f h i "#$%#&'('$)#*$+( +,$*#$*- 38
39 Implementation using Matlab and MPT:" Polytopic Invariant Terminal Set Linear system x + = Ax + Bu LQR controller u = K x Input and state constraints U = {u u min u u max } X = {x x min x x max } M PT X = polytope([i;-i],[xmax; -xmin]) % State constraints U = polytope([k;-k],[umax; -umin]) % Input constraints Compute maximum invariant set that satisfies the constraints O = {x k (A + BK) k x k O, k } M PT P = X & U while 1 Pprev = P; P = P & inv(a+bk)*p; if Pprev == P, break; end end Iteration 2 Iteration 1 Iteration 3 "#$%#&'('$)#*$+( +,$*#$*- 39
40 Implementation using Matlab and MPT:" Polytopic Invariant Terminal Set Linear system x + = Ax + Bu LQR controller u = K x Input and state constraints U = {u u min u u max } X = {x x min x x max } M PT X = polytope([i;-i],[xmax; -xmin]) % State constraints U = polytope([k;-k],[umax; -umin]) % Input constraints Compute maximum invariant set that satisfies the constraints O = {x k (A + BK) k x k O, k } M PT P = X & U while 1 Pprev = P; P = P & inv(a+bk)*p; if Pprev == P, break; end end Iteration 2 Iteration 1 Iteration 3 Simpler procedure M PT X = mpt_infset(a+bk, X&U, 1); "#$%#&'('$)#*$+( +,$*#$*- 4
41 Example: Cessna Citation Aircraft" Revisited MPC controller with input constraints u i ".262 and rate constraints u i.349 approximated by u k u k 1.349T s,u 1 = u prev Altitude x 4 (m) Elevator angle u (rad) Time (sec) Time (sec).5 Pitch angle x 2 (rad) Problem parameters: Sampling time T s =.25sec, Q=I, R=1, N=4 Decrease in the prediction horizon causes loss of the stability properties "#$%#&'('$)#*$+( +,$*#$*- 41
42 Example: Cessna Citation Aircraft" Terminal cost and constraint provide stability guarantee MPC controller with input constraints u i ".262 and rate constraints u i.349 approximated by u k u k 1.349T s,u 1 = u prev Problem parameters: Sampling time T s =.25sec, Q=I, R=1, N=4 Altitude x 4 (m) Elevator angle u (rad) Time (sec) Time (sec).2 Pitch angle x 2 (rad) Inclusion of terminal cost and constraint provides stability "#$%#&'('$)#*$+( +,$*#$*- 42
43 Outline Motivating Example: Cessna Citation Aircraft Constraint satisfaction and stability in MPC Main idea Problem setup for the linear quadratic case How to prove constraint satisfaction and stability in MPC Implementation Theory extensions Reference Tracking Uncertain Systems Soft constrained MPC "#$%#&'('$)#*$+( +,$*#$*- 43
44 Extensions " Reference Tracking Consider the system x + = Ax + Bu Regulation vs. Tracking: y = Cx Regulation: Reject disturbances around one desirable steady-state Tracking: Make output follow a given reference r Steady-state computation: Compute steady-state that yields the desired output [ ][ ] I A B xss C u ss Note: Many solutions possible For a solution to exist the matrix must be full column rank and the steadystate must be feasible, i.e.x ss X,u ss U = [ ] r "#$%#&'('$)#*$+( +,$*#$*- 44
45 Example: Cessna Citation Aircraft" Reference tracking Track given altitude profile: Altitude x 4 (m) Time (sec).5 Pitch angle x 2 (rad) Problem parameters: Sampling time T s =.25sec, Q=I, R=1, N=1 Elevator angle u (rad) Time (sec) "#$%#&'('$)#*$+( +,$*#$*- 45
46 Extensions " Reference Tracking MPC problem for reference tracking: Penalize the deviation from the desired steady-state J (x) = min x,u N 1 (x i x ss ) T Q(x i x ss )+(u i u ss ) T R(u i u ss )+(x N x ss ) T P (x N x ss ) i= x i+1 = Ax i + Bu i Cx i + Du i x b = x Steady-state is computed at each time step Here: Reference assumed constant over horizon. If reference trajectory is known this can be included Preview Same problem structure as the standard MPC problem with additional parameter r Steady-state offset: If the model is not accurate the output will show an offset from the desired reference "#$%#&'('$)#*$+( +,$*#$*- Note: Parameterized terminal set required for stability. 46
47 Extensions " Offset-free reference tracking Goal: Zero steady-state tracking error, i.e. Consider the augmented model: x + = Ax + Bu + B w w y = Cx + C w w with w & R w. [ ] A I Bw Assume size of w =number of states, C w =I and det. C I Then the augmented system is observable. We can correct for the disturbance by choosing x ss, u ss such that [ ][ ] [ ] I A B xss Bw ŵ = C r C w ŵ u ss where ŵ is an estimate of the uncertainty obtained from a state observer. See [Mäder et al., Automatica 29] for more details. y(k) r(k) k "#$%#&'('$)#*$+( +,$*#$*- 47
48 Extensions" Uncertain systems In practice, the nominal system model will not be accurate due to model mismatch or external disturbances that are acting on the system. Consider the uncertain system x + = Ax + Bu+B w w where w &W is a bounded disturbance. Stability: In the presence of uncertainties, asymptotic stability of the origin can often not be achieved. Instead, it can be shown that the trajectories converge to a neighborhood of the origin. Stability can be analyzed in the framework of input-to-state stability. Idea: Effect of the uncertainty on the system is bounded and depends on the size of the uncertainty. Feasibility? "#$%#&'('$)#*$+( +,$*#$*- 48
49 Extensions" Approaches for uncertain systems x terminal set t x + =Ax+Bu+B w w Trajectories differ from what is expected using the nominal system model Uncertainties can lead to infeasibility of the MPC problem Soft constrained MPC: Idea: Tolerate temporary violation of state constraints by constraint relaxation Feasibility of the optimization problem despite disturbances Robust MPC: Idea: Design MPC problem for the worst-case disturbance by constraint tightening Constraint satisfaction and stability of the uncertain system 49 "#$%#&'('$)#*$+( +,$*#$*-
50 Extensions" Soft constrained MPC Idea: Input constraints are hard (e.g. actuator limitations), state constraints may (temporarily) be violated Introduce soft state and terminal constraints by means of slack variables Introduce penalties on the slack variables in the cost J (x) = min x,u N 1 i= x i+1 = Ax i + Bu i Cx i + Du i b + i G N x N x l(x i,u i )+l ( i )+V f (x N )+l ( N ) f N + N = x Slack variables Penalties on the amount of constraint violation Standard soft constrained MPC setup does not provide stability guarantees New soft constrained MPC setup with (robust) stability properties developed in [Zeilinger et al., CDC 21] "#$%#&'('$)#*$+( +,$*#$*- 5
51 Further extensions Reference tracking with recursive feasibility guarantees [Gilbert et al., 1994; Bemporad 1998; Gilbert & Kolmanovsky 1999, 22; Limon et al., 28; Borrelli et al., 29; Ferramosca et al., 29; ] Move blocking: Reduce the computational complexity by fixing the inputs or its derivatives to be constant over several time steps [Li & Xi 27, 29; Cagienard et al. 27; Gondhalekar & Imura, 21; ] Stochastic MPC: Consider uncertainties that are unbounded and/or follow a certain distribution Probabilistic constraints Expected value constraints Expected value cost [Lee et al. 1998; Couchman et al., 25; Cannon et al., 27; Grancharova et al. 27; Hokayem et al., 29; Cinquemani et al., 211; ] "#$%#&'('$)#*$+( +,$*#$*- 51
EE C128 / ME C134 Feedback Control Systems
EE C128 / ME C134 Feedback Control Systems Lecture Additional Material Introduction to Model Predictive Control Maximilian Balandat Department of Electrical Engineering & Computer Science University of
More informationIMPROVED MPC DESIGN BASED ON SATURATING CONTROL LAWS
IMPROVED MPC DESIGN BASED ON SATURATING CONTROL LAWS D. Limon, J.M. Gomes da Silva Jr., T. Alamo and E.F. Camacho Dpto. de Ingenieria de Sistemas y Automática. Universidad de Sevilla Camino de los Descubrimientos
More informationDistributed and Real-time Predictive Control
Distributed and Real-time Predictive Control Melanie Zeilinger Christian Conte (ETH) Alexander Domahidi (ETH) Ye Pu (EPFL) Colin Jones (EPFL) Challenges in modern control systems Power system: - Frequency
More information4F3 - Predictive Control
4F3 Predictive Control - Lecture 3 p 1/21 4F3 - Predictive Control Lecture 3 - Predictive Control with Constraints Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 3 p 2/21 Constraints on
More informationMPC for tracking periodic reference signals
MPC for tracking periodic reference signals D. Limon T. Alamo D.Muñoz de la Peña M.N. Zeilinger C.N. Jones M. Pereira Departamento de Ingeniería de Sistemas y Automática, Escuela Superior de Ingenieros,
More informationPostface to Model Predictive Control: Theory and Design
Postface to Model Predictive Control: Theory and Design J. B. Rawlings and D. Q. Mayne August 19, 2012 The goal of this postface is to point out and comment upon recent MPC papers and issues pertaining
More informationStochastic Tube MPC with State Estimation
Proceedings of the 19th International Symposium on Mathematical Theory of Networks and Systems MTNS 2010 5 9 July, 2010 Budapest, Hungary Stochastic Tube MPC with State Estimation Mark Cannon, Qifeng Cheng,
More informationDecentralized and distributed control
Decentralized and distributed control Centralized control for constrained discrete-time systems M. Farina 1 G. Ferrari Trecate 2 1 Dipartimento di Elettronica, Informazione e Bioingegneria (DEIB) Politecnico
More informationImproved MPC Design based on Saturating Control Laws
Improved MPC Design based on Saturating Control Laws D.Limon 1, J.M.Gomes da Silva Jr. 2, T.Alamo 1 and E.F.Camacho 1 1. Dpto. de Ingenieria de Sistemas y Automática. Universidad de Sevilla, Camino de
More informationEfficient robust optimization for robust control with constraints Paul Goulart, Eric Kerrigan and Danny Ralph
Efficient robust optimization for robust control with constraints p. 1 Efficient robust optimization for robust control with constraints Paul Goulart, Eric Kerrigan and Danny Ralph Efficient robust optimization
More informationOutline. 1 Linear Quadratic Problem. 2 Constraints. 3 Dynamic Programming Solution. 4 The Infinite Horizon LQ Problem.
Model Predictive Control Short Course Regulation James B. Rawlings Michael J. Risbeck Nishith R. Patel Department of Chemical and Biological Engineering Copyright c 217 by James B. Rawlings Outline 1 Linear
More informationMATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem
MATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem Pulemotov, September 12, 2012 Unit Outline Goal 1: Outline linear
More informationFINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES. Danlei Chu, Tongwen Chen, Horacio J. Marquez
FINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES Danlei Chu Tongwen Chen Horacio J Marquez Department of Electrical and Computer Engineering University of Alberta Edmonton
More informationFurther results on Robust MPC using Linear Matrix Inequalities
Further results on Robust MPC using Linear Matrix Inequalities M. Lazar, W.P.M.H. Heemels, D. Muñoz de la Peña, T. Alamo Eindhoven Univ. of Technology, P.O. Box 513, 5600 MB, Eindhoven, The Netherlands,
More informationCourse on Model Predictive Control Part III Stability and robustness
Course on Model Predictive Control Part III Stability and robustness Gabriele Pannocchia Department of Chemical Engineering, University of Pisa, Italy Email: g.pannocchia@diccism.unipi.it Facoltà di Ingegneria,
More informationLearning Model Predictive Control for Iterative Tasks: A Computationally Efficient Approach for Linear System
Learning Model Predictive Control for Iterative Tasks: A Computationally Efficient Approach for Linear System Ugo Rosolia Francesco Borrelli University of California at Berkeley, Berkeley, CA 94701, USA
More informationNonlinear Model Predictive Control Tools (NMPC Tools)
Nonlinear Model Predictive Control Tools (NMPC Tools) Rishi Amrit, James B. Rawlings April 5, 2008 1 Formulation We consider a control system composed of three parts([2]). Estimator Target calculator Regulator
More informationIEOR 265 Lecture 14 (Robust) Linear Tube MPC
IEOR 265 Lecture 14 (Robust) Linear Tube MPC 1 LTI System with Uncertainty Suppose we have an LTI system in discrete time with disturbance: x n+1 = Ax n + Bu n + d n, where d n W for a bounded polytope
More informationRobustly stable feedback min-max model predictive control 1
Robustly stable feedback min-max model predictive control 1 Eric C. Kerrigan 2 and Jan M. Maciejowski Department of Engineering, University of Cambridge Trumpington Street, Cambridge CB2 1PZ, United Kingdom
More information4F3 - Predictive Control
4F3 Predictive Control - Lecture 2 p 1/23 4F3 - Predictive Control Lecture 2 - Unconstrained Predictive Control Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 2 p 2/23 References Predictive
More informationPrinciples of Optimal Control Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 6.33 Principles of Optimal Control Spring 8 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 6.33 Lecture 6 Model
More informationModel Predictive Control Short Course Regulation
Model Predictive Control Short Course Regulation James B. Rawlings Michael J. Risbeck Nishith R. Patel Department of Chemical and Biological Engineering Copyright c 2017 by James B. Rawlings Milwaukee,
More informationIntroduction to Model Predictive Control. Dipartimento di Elettronica e Informazione
Introduction to Model Predictive Control Riccardo Scattolini Riccardo Scattolini Dipartimento di Elettronica e Informazione Finite horizon optimal control 2 Consider the system At time k we want to compute
More informationROBUST CONSTRAINED PREDICTIVE CONTROL OF A 3DOF HELICOPTER MODEL WITH EXTERNAL DISTURBANCES
ABCM Symposium Series in Mechatronics - Vol 3 - pp19-26 Copyright c 28 by ABCM ROBUST CONSTRAINED PREDICTIVE CONTROL OF A 3DOF HELICOPTER MODEL WITH EXTERNAL DISTURBANCES Marcelo Handro Maia, handro@itabr
More informationOn the Inherent Robustness of Suboptimal Model Predictive Control
On the Inherent Robustness of Suboptimal Model Predictive Control James B. Rawlings, Gabriele Pannocchia, Stephen J. Wright, and Cuyler N. Bates Department of Chemical & Biological Engineering Computer
More informationFast Model Predictive Control with Soft Constraints
European Control Conference (ECC) July 7-9,, Zürich, Switzerland. Fast Model Predictive Control with Soft Constraints Arthur Richards Department of Aerospace Engineering, University of Bristol Queens Building,
More informationCourse on Model Predictive Control Part II Linear MPC design
Course on Model Predictive Control Part II Linear MPC design Gabriele Pannocchia Department of Chemical Engineering, University of Pisa, Italy Email: g.pannocchia@diccism.unipi.it Facoltà di Ingegneria,
More informationPart II: Model Predictive Control
Part II: Model Predictive Control Manfred Morari Alberto Bemporad Francesco Borrelli Unconstrained Optimal Control Unconstrained infinite time optimal control V (x 0 )= inf U k=0 [ x k Qx k + u k Ru k]
More informationMODEL PREDICTIVE SLIDING MODE CONTROL FOR CONSTRAINT SATISFACTION AND ROBUSTNESS
MODEL PREDICTIVE SLIDING MODE CONTROL FOR CONSTRAINT SATISFACTION AND ROBUSTNESS Yizhou Wang, Wenjie Chen, Masayoshi Tomizuka Department of Mechanical Engineering University of California Berkeley, California
More informationNonlinear Reference Tracking with Model Predictive Control: An Intuitive Approach
onlinear Reference Tracking with Model Predictive Control: An Intuitive Approach Johannes Köhler, Matthias Müller, Frank Allgöwer Abstract In this paper, we study the system theoretic properties of a reference
More informationStability and feasibility of MPC for switched linear systems with dwell-time constraints
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Stability and feasibility of MPC for switched linear systems with dwell-time constraints Bridgeman, L.; Danielson, C.; Di Cairano, S. TR016-045
More informationOn the Inherent Robustness of Suboptimal Model Predictive Control
On the Inherent Robustness of Suboptimal Model Predictive Control James B. Rawlings, Gabriele Pannocchia, Stephen J. Wright, and Cuyler N. Bates Department of Chemical and Biological Engineering and Computer
More informationOn robustness of suboptimal min-max model predictive control *
Manuscript received June 5, 007; revised Sep., 007 On robustness of suboptimal min-max model predictive control * DE-FENG HE, HAI-BO JI, TAO ZHENG Department of Automation University of Science and Technology
More informationOutline. Linear regulation and state estimation (LQR and LQE) Linear differential equations. Discrete time linear difference equations
Outline Linear regulation and state estimation (LQR and LQE) James B. Rawlings Department of Chemical and Biological Engineering 1 Linear Quadratic Regulator Constraints The Infinite Horizon LQ Problem
More informationAdaptive Nonlinear Model Predictive Control with Suboptimality and Stability Guarantees
Adaptive Nonlinear Model Predictive Control with Suboptimality and Stability Guarantees Pontus Giselsson Department of Automatic Control LTH Lund University Box 118, SE-221 00 Lund, Sweden pontusg@control.lth.se
More informationA new low-and-high gain feedback design using MPC for global stabilization of linear systems subject to input saturation
A new low-and-high gain feedbac design using MPC for global stabilization of linear systems subject to input saturation Xu Wang 1 Håvard Fjær Grip 1; Ali Saberi 1 Tor Arne Johansen Abstract In this paper,
More informationLinear-Quadratic Optimal Control: Full-State Feedback
Chapter 4 Linear-Quadratic Optimal Control: Full-State Feedback 1 Linear quadratic optimization is a basic method for designing controllers for linear (and often nonlinear) dynamical systems and is actually
More informationA Stable Block Model Predictive Control with Variable Implementation Horizon
American Control Conference June 8-,. Portland, OR, USA WeB9. A Stable Block Model Predictive Control with Variable Implementation Horizon Jing Sun, Shuhao Chen, Ilya Kolmanovsky Abstract In this paper,
More informationEnlarged terminal sets guaranteeing stability of receding horizon control
Enlarged terminal sets guaranteeing stability of receding horizon control J.A. De Doná a, M.M. Seron a D.Q. Mayne b G.C. Goodwin a a School of Electrical Engineering and Computer Science, The University
More informationContinuous-time linear MPC algorithms based on relaxed logarithmic barrier functions
Preprints of the 19th World Congress The International Federation of Automatic Control Continuous-time linear MPC algorithms based on relaxed logarithmic barrier functions Christian Feller Christian Ebenbauer
More informationA Globally Stabilizing Receding Horizon Controller for Neutrally Stable Linear Systems with Input Constraints 1
A Globally Stabilizing Receding Horizon Controller for Neutrally Stable Linear Systems with Input Constraints 1 Ali Jadbabaie, Claudio De Persis, and Tae-Woong Yoon 2 Department of Electrical Engineering
More informationA SIMPLE TUBE CONTROLLER FOR EFFICIENT ROBUST MODEL PREDICTIVE CONTROL OF CONSTRAINED LINEAR DISCRETE TIME SYSTEMS SUBJECT TO BOUNDED DISTURBANCES
A SIMPLE TUBE CONTROLLER FOR EFFICIENT ROBUST MODEL PREDICTIVE CONTROL OF CONSTRAINED LINEAR DISCRETE TIME SYSTEMS SUBJECT TO BOUNDED DISTURBANCES S. V. Raković,1 D. Q. Mayne Imperial College London, London
More informationLecture 9: Discrete-Time Linear Quadratic Regulator Finite-Horizon Case
Lecture 9: Discrete-Time Linear Quadratic Regulator Finite-Horizon Case Dr. Burak Demirel Faculty of Electrical Engineering and Information Technology, University of Paderborn December 15, 2015 2 Previous
More informationOptimized Robust Control Invariant Tubes. & Robust Model Predictive Control
Optimized Robust Control Invariant Tubes & Robust Model Predictive Control S. V. Raković (sasa.rakovic@imperial.ac.uk) Imperial College, London Imperial College, 17. November 2004 Imperial College ICM
More informationLINEAR TIME VARYING TERMINAL LAWS IN MPQP
LINEAR TIME VARYING TERMINAL LAWS IN MPQP JA Rossiter Dept of Aut Control & Systems Eng University of Sheffield, Mappin Street Sheffield, S1 3JD, UK email: JARossiter@sheffieldacuk B Kouvaritakis M Cannon
More informationNonlinear Model Predictive Control for Periodic Systems using LMIs
Marcus Reble Christoph Böhm Fran Allgöwer Nonlinear Model Predictive Control for Periodic Systems using LMIs Stuttgart, June 29 Institute for Systems Theory and Automatic Control (IST), University of Stuttgart,
More informationARTICLE IN PRESS Journal of Process Control xxx (2011) xxx xxx
Journal of Process Control xxx (2011) xxx xxx Contents lists available at ScienceDirect Journal of Process Control ournal homepage: www.elsevier.com/locate/procont Stable MPC for tracking with maximal
More informationPrashant Mhaskar, Nael H. El-Farra & Panagiotis D. Christofides. Department of Chemical Engineering University of California, Los Angeles
HYBRID PREDICTIVE OUTPUT FEEDBACK STABILIZATION OF CONSTRAINED LINEAR SYSTEMS Prashant Mhaskar, Nael H. El-Farra & Panagiotis D. Christofides Department of Chemical Engineering University of California,
More informationGiulio Betti, Marcello Farina and Riccardo Scattolini
1 Dipartimento di Elettronica e Informazione, Politecnico di Milano Rapporto Tecnico 2012.29 An MPC algorithm for offset-free tracking of constant reference signals Giulio Betti, Marcello Farina and Riccardo
More informationRobust Learning Model Predictive Control for Uncertain Iterative Tasks: Learning From Experience
Robust Learning Model Predictive Control for Uncertain Iterative Tasks: Learning From Experience Ugo Rosolia, Xiaojing Zhang, and Francesco Borrelli Abstract We present a Robust Learning Model Predictive
More informationReal Time Economic Dispatch for Power Networks: A Distributed Economic Model Predictive Control Approach
Real Time Economic Dispatch for Power Networks: A Distributed Economic Model Predictive Control Approach Johannes Köhler, Matthias A. Müller, Na Li, Frank Allgöwer Abstract Fast power fluctuations pose
More informationRobust Adaptive MPC for Systems with Exogeneous Disturbances
Robust Adaptive MPC for Systems with Exogeneous Disturbances V. Adetola M. Guay Department of Chemical Engineering, Queen s University, Kingston, Ontario, Canada (e-mail: martin.guay@chee.queensu.ca) Abstract:
More informationOptimizing Economic Performance using Model Predictive Control
Optimizing Economic Performance using Model Predictive Control James B. Rawlings Department of Chemical and Biological Engineering Second Workshop on Computational Issues in Nonlinear Control Monterey,
More informationTime-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2015
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 15 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
More informationarxiv: v2 [math.oc] 15 Jan 2014
Stability and Performance Guarantees for MPC Algorithms without Terminal Constraints 1 Jürgen Pannek 2 and Karl Worthmann University of the Federal Armed Forces, 85577 Munich, Germany, juergen.pannek@googlemail.com
More informationMPC Infeasibility Handling
MPC Handling Thomas Wiese, TU Munich, KU Leuven supervised by H.J. Ferreau, Prof. M. Diehl (both KUL) and Dr. H. Gräb (TUM) October 9, 2008 1 / 42 MPC General MPC Strategies 2 / 42 Linear Discrete-Time
More informationEconomic model predictive control with self-tuning terminal weight
Economic model predictive control with self-tuning terminal weight Matthias A. Müller, David Angeli, and Frank Allgöwer Abstract In this paper, we propose an economic model predictive control (MPC) framework
More informationHomework Solution # 3
ECSE 644 Optimal Control Feb, 4 Due: Feb 17, 4 (Tuesday) Homework Solution # 3 1 (5%) Consider the discrete nonlinear control system in Homework # For the optimal control and trajectory that you have found
More informationEconomic and Distributed Model Predictive Control: Recent Developments in Optimization-Based Control
SICE Journal of Control, Measurement, and System Integration, Vol. 10, No. 2, pp. 039 052, March 2017 Economic and Distributed Model Predictive Control: Recent Developments in Optimization-Based Control
More informationOptimal control and estimation
Automatic Control 2 Optimal control and estimation Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011
More informationECE7850 Lecture 9. Model Predictive Control: Computational Aspects
ECE785 ECE785 Lecture 9 Model Predictive Control: Computational Aspects Model Predictive Control for Constrained Linear Systems Online Solution to Linear MPC Multiparametric Programming Explicit Linear
More informationDecentralized and distributed control
Decentralized and distributed control Constrained distributed control for discrete-time systems M. Farina 1 G. Ferrari Trecate 2 1 Dipartimento di Elettronica e Informazione (DEI) Politecnico di Milano,
More informationarxiv: v1 [cs.sy] 28 May 2013
From Parametric Model-based Optimization to robust PID Gain Scheduling Minh H.. Nguyen a,, K.K. an a a National University of Singapore, Department of Electrical and Computer Engineering, 3 Engineering
More informationRegional Solution of Constrained LQ Optimal Control
Regional Solution of Constrained LQ Optimal Control José DeDoná September 2004 Outline 1 Recap on the Solution for N = 2 2 Regional Explicit Solution Comparison with the Maximal Output Admissible Set 3
More informationTime-Invariant Linear Quadratic Regulators!
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 17 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
More informationNumerical Methods for Model Predictive Control. Jing Yang
Numerical Methods for Model Predictive Control Jing Yang Kongens Lyngby February 26, 2008 Technical University of Denmark Informatics and Mathematical Modelling Building 321, DK-2800 Kongens Lyngby, Denmark
More informationCONSTRAINED MODEL PREDICTIVE CONTROL ON CONVEX POLYHEDRON STOCHASTIC LINEAR PARAMETER VARYING SYSTEMS. Received October 2012; revised February 2013
International Journal of Innovative Computing, Information and Control ICIC International c 2013 ISSN 1349-4198 Volume 9, Number 10, October 2013 pp 4193 4204 CONSTRAINED MODEL PREDICTIVE CONTROL ON CONVEX
More informationQuadratic Stability of Dynamical Systems. Raktim Bhattacharya Aerospace Engineering, Texas A&M University
.. Quadratic Stability of Dynamical Systems Raktim Bhattacharya Aerospace Engineering, Texas A&M University Quadratic Lyapunov Functions Quadratic Stability Dynamical system is quadratically stable if
More informationExam. 135 minutes, 15 minutes reading time
Exam August 6, 208 Control Systems II (5-0590-00) Dr. Jacopo Tani Exam Exam Duration: 35 minutes, 5 minutes reading time Number of Problems: 35 Number of Points: 47 Permitted aids: 0 pages (5 sheets) A4.
More informationUCLA Chemical Engineering. Process & Control Systems Engineering Laboratory
Constrained Innite-time Optimal Control Donald J. Chmielewski Chemical Engineering Department University of California Los Angeles February 23, 2000 Stochastic Formulation - Min Max Formulation - UCLA
More informationarzelier
COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.1 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS STABILITY ANALYSIS Didier HENRION www.laas.fr/ henrion henrion@laas.fr Denis ARZELIER www.laas.fr/
More informationRobust output feedback model predictive control of constrained linear systems
Automatica 42 (2006) 1217 1222 Brief paper Robust output feedback model predictive control of constrained linear systems D.Q. Mayne a, S.V. Raković a,, R. Findeisen b, F. Allgöwer b a Department of Electrical
More informationMODEL PREDICTIVE CONTROL FUNDAMENTALS
Nigerian Journal of Technology (NIJOTECH) Vol 31, No 2, July, 2012, pp 139 148 Copyright 2012 Faculty of Engineering, University of Nigeria ISSN 1115-8443 MODEL PREDICTIVE CONTROL FUNDAMENTALS PE Orukpe
More informationOn the design of Robust tube-based MPC for tracking
Proceedings of the 17th World Congress The International Federation of Automatic Control On the design of Robust tube-based MPC for tracking D. Limon I. Alvarado T. Alamo E. F. Camacho Dpto. de Ingeniería
More informationTube Model Predictive Control Using Homothety & Invariance
Tube Model Predictive Control Using Homothety & Invariance Saša V. Raković rakovic@control.ee.ethz.ch http://control.ee.ethz.ch/~srakovic Collaboration in parts with Mr. Mirko Fiacchini Automatic Control
More informationOptimal Control. Lecture 18. Hamilton-Jacobi-Bellman Equation, Cont. John T. Wen. March 29, Ref: Bryson & Ho Chapter 4.
Optimal Control Lecture 18 Hamilton-Jacobi-Bellman Equation, Cont. John T. Wen Ref: Bryson & Ho Chapter 4. March 29, 2004 Outline Hamilton-Jacobi-Bellman (HJB) Equation Iterative solution of HJB Equation
More informationThe ϵ-capacity of a gain matrix and tolerable disturbances: Discrete-time perturbed linear systems
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 11, Issue 3 Ver. IV (May - Jun. 2015), PP 52-62 www.iosrjournals.org The ϵ-capacity of a gain matrix and tolerable disturbances:
More informationA FAST, EASILY TUNED, SISO, MODEL PREDICTIVE CONTROLLER. Gabriele Pannocchia,1 Nabil Laachi James B. Rawlings
A FAST, EASILY TUNED, SISO, MODEL PREDICTIVE CONTROLLER Gabriele Pannocchia, Nabil Laachi James B. Rawlings Department of Chemical Engineering Univ. of Pisa Via Diotisalvi 2, 5626 Pisa (Italy) Department
More informationOptimal Control. Lecture 3. Optimal Control of Discrete Time Dynamical Systems. John T. Wen. January 22, 2004
Optimal Control Lecture 3 Optimal Control of Discrete Time Dynamical Systems John T. Wen January, 004 Outline optimization of a general multi-stage discrete time dynamical systems special case: discrete
More informationPiecewise-affine functions: applications in circuit theory and control
Piecewise-affine functions: applications in circuit theory and control Tomaso Poggi Basque Center of Applied Mathematics Bilbao 12/04/2013 1/46 Outline 1 Embedded systems 2 PWA functions Definition Classes
More informationESTIMATES ON THE PREDICTION HORIZON LENGTH IN MODEL PREDICTIVE CONTROL
ESTIMATES ON THE PREDICTION HORIZON LENGTH IN MODEL PREDICTIVE CONTROL K. WORTHMANN Abstract. We are concerned with model predictive control without stabilizing terminal constraints or costs. Here, our
More informationMulti-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures
Preprints of the 19th World Congress The International Federation of Automatic Control Multi-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures Eric Peterson Harry G.
More informationA CHANCE CONSTRAINED ROBUST MPC AND APPLICATION TO HOT ROLLING MILLS
A CHANCE CONSTRAINED ROBUST MPC AND APPLICATION TO HOT ROLLING MILLS Il Seop Choi Anthony Rossiter Peter Fleming University of Sheffield, Department of Automatic Control & Systems Engineering, Sheffield
More informationLinear Systems. Manfred Morari Melanie Zeilinger. Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich
Linear Systems Manfred Morari Melanie Zeilinger Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich Spring Semester 2016 Linear Systems M. Morari, M. Zeilinger - Spring
More informationFinite horizon robust model predictive control with terminal cost constraints
Finite horizon robust model predictive control with terminal cost constraints Danlei Chu, Tongwen Chen and Horacio J Marquez Department of Electrical & Computer Engineering, University of Alberta, Canada,
More informationMPC: Tracking, Soft Constraints, Move-Blocking
MPC: Tracking, Soft Constraints, Move-Blocking M. Morari, F. Borrelli, C. Jones, M. Zeilinger Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich UC Berkeley EPFL Spring
More informationRobust Explicit MPC Based on Approximate Multi-parametric Convex Programming
43rd IEEE Conference on Decision and Control December 4-7, 24 Atlantis, Paradise Island, Bahamas WeC6.3 Robust Explicit MPC Based on Approximate Multi-parametric Convex Programming D. Muñoz de la Peña
More informationLinear Offset-Free Model Predictive Control
Linear Offset-Free Model Predictive Control Urban Maeder a,, Francesco Borrelli b, Manfred Morari a a Automatic Control Lab, ETH Zurich, CH-892 Zurich, Switzerland b Department of Mechanical Engineering,
More informationTechnical work in WP2 and WP5
Technical work in WP2 and WP5 UNIZG-FER Mato Baotić, Branimir Novoselnik, Jadranko Matuško, Mario Vašak, Andrej Jokić Aachen, October 9, 2015 This project has received funding from the European Union s
More informationPrediktivno upravljanje primjenom matematičkog programiranja
Prediktivno upravljanje primjenom matematičkog programiranja Doc. dr. sc. Mato Baotić Fakultet elektrotehnike i računarstva Sveučilište u Zagrebu www.fer.hr/mato.baotic Outline Application Examples PredictiveControl
More informationMOST control systems are designed under the assumption
2076 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 53, NO. 9, OCTOBER 2008 Lyapunov-Based Model Predictive Control of Nonlinear Systems Subject to Data Losses David Muñoz de la Peña and Panagiotis D. Christofides
More informationECE7850 Lecture 8. Nonlinear Model Predictive Control: Theoretical Aspects
ECE7850 Lecture 8 Nonlinear Model Predictive Control: Theoretical Aspects Model Predictive control (MPC) is a powerful control design method for constrained dynamical systems. The basic principles and
More informationPenalty and Barrier Methods General classical constrained minimization problem minimize f(x) subject to g(x) 0 h(x) =0 Penalty methods are motivated by the desire to use unconstrained optimization techniques
More informationESC794: Special Topics: Model Predictive Control
ESC794: Special Topics: Model Predictive Control Discrete-Time Systems Hanz Richter, Professor Mechanical Engineering Department Cleveland State University Discrete-Time vs. Sampled-Data Systems A continuous-time
More informationFast model predictive control based on linear input/output models and bounded-variable least squares
7 IEEE 56th Annual Conference on Decision and Control (CDC) December -5, 7, Melbourne, Australia Fast model predictive control based on linear input/output models and bounded-variable least squares Nilay
More informationIndirect Adaptive Model Predictive Control for Linear Systems with Polytopic Uncertainty
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Indirect Adaptive Model Predictive Control for Linear Systems with Polytopic Uncertainty Di Cairano, S. TR2016-043 July 2016 Abstract We develop
More informationEconomic MPC using a Cyclic Horizon with Application to Networked Control Systems
Economic MPC using a Cyclic Horizon with Application to Networked Control Systems Stefan Wildhagen 1, Matthias A. Müller 1, and Frank Allgöwer 1 arxiv:1902.08132v1 [cs.sy] 21 Feb 2019 1 Institute for Systems
More informationLecture 13: Constrained optimization
2010-12-03 Basic ideas A nonlinearly constrained problem must somehow be converted relaxed into a problem which we can solve (a linear/quadratic or unconstrained problem) We solve a sequence of such problems
More informationCourse Outline. Higher Order Poles: Example. Higher Order Poles. Amme 3500 : System Dynamics & Control. State Space Design. 1 G(s) = s(s + 2)(s +10)
Amme 35 : System Dynamics Control State Space Design Course Outline Week Date Content Assignment Notes 1 1 Mar Introduction 2 8 Mar Frequency Domain Modelling 3 15 Mar Transient Performance and the s-plane
More informationAn Introduction to Model-based Predictive Control (MPC) by
ECE 680 Fall 2017 An Introduction to Model-based Predictive Control (MPC) by Stanislaw H Żak 1 Introduction The model-based predictive control (MPC) methodology is also referred to as the moving horizon
More information