Infinite-dimensional nonlinear predictive controller design for open-channel hydraulic systems
|
|
- Andrea French
- 6 years ago
- Views:
Transcription
1 Infinite-dimensional nonlinear predictive controller design for open-channel hydraulic systems D. Georges, Control Systems Dept - Gipsa-lab, Grenoble INP Workshop on Irrigation Channels and Related Problems, Salerno, Italy, october 2-4, 28
2 Outline of the talk 1. Control of a single canal pool 2. Modelling of open-channel hydraulic systems: Saint Venant PDE 3. Some backgrounds on predictive control 4. Formulation of the associated optimal control problem 5. Necessary conditions for optimality: the adjoint PDEs 6. Computation of the related two-point boundary value problem 7. Description of the here-proposed predictive control scheme 8. Some simulation results 9. The multi-pool case: a decomposition approach via Lagrangian relaxation 1. Some conclusions and perspectives
3 Control problem Regulation of a single pool of an irrigation canal (pool = canal section delimited by 2 gates) using the upstream regulator gate: B L D h h S h Q P I Upstream Downstream Transversal section Longitudinal section
4 Why to use nonlinear predictive control? 1. In order to take nonlinear transportation and diffusion phenomena into account; 2. Because (finite-dimensional) nonlinear optimal control provides, under some specific assumptions, stabilizing feedbacks (infinite-horizon optimal control, predictive control): possible extension in the infinite-dimensional framework?
5 Open-channel modelling: Saint Venant PDE Rectangular canal section case: S : { B z t + Q x x ( Q2 Bz Q t + = q ) + gbz z x gbz(i J(Q, z)) = kq Q Bz (1) z = relative water level,(m); B=canal width,(m); Q = water flow rate, ( m3 s ); A = wet section, (m2 ); I = canal slope ; g = acc. of gravity, ( m ); J = friction ( m s 2 m ) ; q(x)= withdrawal / length unit ( m2 s ) ; k=, if q >, 1, if q <. { Q(, t) = u(t) (integral control : upstream) B.C. z(l, t) = v(t) (disturbance : downstream) Models of the gates may be introduced (2)
6 Some backgrounds on predictive control Consider a nonlinear system of the form: ẋ = F (x, u), x() = x (3) with F (, ) = (the origin is an equilibrium point). 1. At time t, x(t) = x t is measured or computed from a state observer. 2. Computation of an (open-loop) optimal control problem over a receding control interval [t, t + T ] starting with state x t (so-called associated optimal control problem (AOCP)): min u(.) t+t t L(x(τ), u(τ))dτ (4) s.t. ẋ = F (x, u), x(t) = x t and x(t + T ) E (5) where E denotes a domain of the state space in which x(t + T ) is assigned. 3. Application of the first control u(t) obtained from the sequence {u(t), u(t + t), u(t + 2 t),..., u(t + N t)}. 4. Go back to 1), t + t t.
7 Some backgrounds on predictive control Some sufficient conditions to ensure local asymptotic stability around the origin (Mayne & Michalska, TAC 9): L(.,.) is positive definite in both x and u (in the sense of Hessian matrices w.r.t. x and u) under assumptions (mainly, F is 2 times cont. differentiabble, F is uniformly Lipschitz in x w.r.t. u, local uniform controllability) ensuring existence and uniqueness of the AOCP E = {} V (x, t) = t+t t L(x(τ), u(τ))dτ is a control Lyapunov function. Extension to the infinite-dimensional case?
8 Associated optimal control problem s.t. with T T min m(u(t))dt + u S : B.C. B z t + Q x Q t + x ( Q2 Bz µ = u I.C. = q L l(z(x, t), Q(x, t))dxdt (6) ) + gbz z x gbz(i J) = kq Q Bz Q(x, ) = φ 1 (x), x [, L] z(x, ) = φ 2 (x) µ() = φ 1 () { Q(, t) = µ(t) (control : upstream end) z(l, t) = v(t) (disturbance : downstream end) and constraint on state at final time (terminal constraint): { Q(x, T ) = Q (x), x [, L] T.C. z(x, T ) = z (x) (7) (8) (9) (1)
9 Regulation around an equilibrium state Regularization term: m(u(t)) = r 2 u(t)2, r > (11) State cost function: l(z(x, t), Q(x, t)) = 1 2 [q 1(x)(z(x, t) z (x)) 2 (12) +q 2 (x)(q(x, t) Q (x)) 2 ], q 1 (.), q 2 (.) > where z (x), Q (x), x [, L] = equilibrium state, solution of: S e : { Q x = q x ( Q2 z Bz ) + gbz x gbz(i J) = kq Q Bz with Boundary Conditions: (13) B.C. { Q () = Q o z (L) = z L (14)
10 Some remarks Some other possible control objectives by choosing other m and l functionals: Regulation of the downstream-end level z(x = L, t), with upstream-end boundary control and a downstream-end boundary condition Q(x = L, t) = v(t), where v is the downstream-end water flow acting as a disturbance. Regulation of the downstream-end level z(x = L, t), with downstream-end boundary control and an upstream-end boundary condition Q(x =, t) = v(t), where v is the upstream-end water flow acting as a disturbance. Control problems with withdrawals may be also considered without restriction.
11 Lagrangian formulation L(Q, z, λ 1, λ 2, λ 3, u) = with S : + + Q S 1 = B z t S 2 = Q t + T T L x q x ( Q2 Bz = Q t + x ( Q2 S 3 = µ u {m(u) + λ 3 S 3 }dt (15) {l(z, Q) + λ 1 S 1 + λ 2 S 2 }dxdt ) + gbz z x gbz(i J) kq Q Bz Bz gbz2 ) gbz(i J) kq Q Bz (16)
12 Necessary conditions for optimality: some adjoint PDEs (D. Georges & M.L. Chen, ECC 1) derived from the Lagrangian formulation + variational calculus: In order that both u(t) [, T ] and the trajectory of system: B z = q with and S : are optimal, t + Q x Q t + x ( Q2 Bz µ = u ) + gbz z x gbz(i J) = kq Q Bz Q(x, ) = φ 1 (x), x [, L] I.C. z(x, ) = φ 2 (x) µ() = φ 1 () { Q(, t) = µ(t) (control) B.C. z(l, t) = v(t) (disturbance) T.C. Q(x, T ) = Q (x), x [, L] z(x, T ) = z (x) µ(t ) = Q (x = ) (17) (18) (19) (2)
13 it is necessary that there exists λ = (λ 1 (x, t), λ 2 (x, t), λ 3 (t)), solution of the adjoint system of S: l z B λ 1 t + λ 2 x ( Q2 gbz) λ Bz 2 2 (gb(i J) gbz J z kq Q Bz 2 ) = S adj : l Q λ 1 x λ 2 t 2 λ 2 Q x Bz λ 2(kq 1 J gbz Bz Q ) = Q(, t) λ 1 (, t) + 2λ 2 (, t) Bz(, t) + λ 3 = with and I.T.C. B.C. λ 1 (x, ), λ 1 (x, T ), free, x [, L] λ 2 (x, ), λ 2 (x, T ), free λ 3 (), λ 3 (T ) free (21) { λ2 (, t) =, t [, T ] λ 1 (L, t) + 2λ 2 (L, t) Q(L,t) Bz(L,t) = (22)
14 For all t [, T ], the optimal control u = u, is necessary solution of: m (u(t)) λ 3 (t) = A two-point boundary value problem
15 Computation of the two-point boundary value problem Computation in two stages: 1. Spatial and temporal discretization of the canonical equations: Preissman numerical scheme (1962): applied to both S and S adj ; 2. Solution of the nonlinear algebraic equations derived from the two-dimensional grid via a Newton-Raphson method.
16 Preissman scheme: a finite-difference scheme semi-implicit in time Approximation of fonctions f and their derivatives: f (x, t) = 1 θ 2 [f i+1 + f i ] + θ 2 [f + i+1 + f + i ] f 1 θ (x, t) = x x [f i+1 f i ] + θ x [f + i+1 f + i ] f t (x, t) = 1 2 t [f + i f i + f + i+1 f i+1] (23) where i = spatial index, + = t + t and θ 1 = relaxation coefficient. If θ, 5, we get an unconditionally stable integration scheme.
17 An algebraic set of nonlinear equations For the case when m(u) = 1 2 u2, if N is the number of spatial discretization points and M, the number of temporal discretization points: S + B.C. (2N + 1) M unknown variables (z(x, t), Q(x, t), µ(t)) with (2N + 1) (M 1) equations + 4N + 2 constraints (the states (z, Q, µ) are imposed at t = and t = T ). S adj + B.C. (2N + 1) M unknown variables (λ 1 (x, t), λ 2 (x, t), λ 3 (t)) with (2N + 1) (M 1) equations. 4N M + 2M variables (z(x, t), Q(x, t), µ(t), λ 1 (x, t), λ 2 (x, t), λ 3 (t)) for 4N M + 2M equations: an implicit problem Computation of both state and adjoint state via a Newton-Raphson scheme.
18 The predictive control scheme 1. At time t : computation of the open-loop optimal control problem, with both initial and final states fixed: min u s.t. with S : t+t t m(u(t))dt + B z t + Q x Q t + x ( Q2 Bz µ = u I.C. t+t L t = q l(z(x, t), Q(x, t))dxdt (24) ) + gbz z x gbz(i J) = kq Q Bz Q(x, ) = Q t (x), x [, L] z(x, ) = z t (x) µ() = Q t () (25) (26) Q t (x) et z t (x) = states at time t on [, L], { Q(, τ) = µ(τ), τ [t, t + T ] (control) B.C. z(l, τ) = ˆv(τ) (prediction of downstream end level) (27)
19 The predictive control scheme 2 Apply u(t) defined as the first optimal control input sequence computed on [t, t + T ]: the system reaches state (z(x, t + t), Q(x, t + t)), x [, L]. 3 Go to 1), with t + t t and get z t (x) = z(x, t + t), Q t (x) = Q(x, t + t).
20 Some simulation results Simulation based on Preissman numerical scheme; A 5 km long canal divided into 1 sections of 5 m each; Regulation around a uniform equilibrium state corresponding to a constant relative water level z of 1.5 m, along the pool; Simulation starting from a uniform equilibrium of 1 m; At each time step: complexity of the two-point boundary value problem: N = 11 (spatial discretization) and M = 6 (temporal discretization, t = 1s): 4N M + 2M = 288 equations for 4N M + 2M = 288 unknown variables: computation time < 6 s on a Pentium 1,8 Mhz, 512 Mo Laptop: real-time control is possible.
21 Simulations results 3.9 Water flow rate 3.8 in m3/s Time in mns 1.8 Water level in meters Time in mns
22 Simulations results
23 The multi-pool case A two-pool case: (easy extension to n pools) : Regulator gate x = Pool 1 x = L 1 x + = L 1 Pool 2 x = L 2 Dynamics of the 2 pools: 2 X 2 PDE coupled by a gate model: Q 1 (L 1, t) = Q 2(L + 1, t) = Q(t) and Q such that Q 2 = K 2 α 2 2g(z 1 (L 1, t) z 2(L + 1, t)) F (Q(t), z 1 (L 1, t), z 2(L + 1, t), α(t)) =
24 Two-pool modelling S pool1 : S pool2 : B z 1 t + Q 1 Q 1 t x = q 1 + x ( Q2 1 z Bz 1 ) + gbz 1 1 x gbz Q 1(I J) = kq 1 1 Bz 1 µ = u 1 { (29) Q B.C. 1 (, t) = µ(t) Q 1 (L 1, t) = Q(t) (3) B z 2 t + Q 2 Q 2 t α = u 2 x = q 2 + x ( Q2 2 z Bz 2 ) + gbz 2 2 x gbz Q 2(I J) = kq 2 2 Bz 2 { Q2 (L + B.C. 1, t) = Q(t) z 2 (L 2, t) = v(t) (31) (32) Q(t) 2 = K 2 α 2 2g(z 1 (L 1, t) z 2(L + 1, t)), (33)
25 Optimal control formulation min u 2 [ i=1 T T bi m i (u i )dt + l i (z i (x, t), Q i (x, t))dxdt] (34) a i s.t. S pool1 + B.C. + I.C.+ T.C. and S pool2 + B.C + I.C.+ T.C. and the additional algebraic constraint defined on [, T ]: Q 2 = K 2 α 2 2g(z 1 (L 1, t) z 2(L + 1, t)).
26 Main issues How to reduce the computational complexity? How to take advantage of distributed control architecture (supervisory control and data acquisition: SCADA) used in large-scale water distribution systems? Here-proposed solution: use of a decomposition-coordination algorithm based on Lagrangian relaxation (D. Georges, Workshop NMPC 6)
27 An augmented lagrangian formulation = L c (Q 1, z 1, Q 2, z 2, Q, λ 1 1, λ 1 2, λ 1 3, λ 2 1, λ 2 2, λ 2 3, p, u 1, u 2 ) T T T {[m 1 (u 1 ) + λ 1 3[ Q 1 (, t) u 1 (t)]}dt L 1 T L2 {l 1 (z 1, Q 1 ) + λ 1 1S λ 1 2S 1 2 }dxdt {[m 2 (u 2 ) + λ 2 3[ α u 2 (t)]}dt L + 1 {l 2 (z 2, Q 2 ) + λ 2 1S λ 2 2S 2 2 }dxdt (35)
28 + T [(p + c 2 F (Q(t), z 1(L 1, t), z 2(L + 1, t))) F (Q(t), z 1 (L 1, t), z 2(L + 1, t), α)]dt (36) p is the dual variable associated to the gate constraint F an additional quadratic term c 2 F 2 is introduced An augmented Lagrangian L c
29 Main improvement of augmented Lagrangians: For sufficiently large c >, the augmented Lagrangian admits at least one local saddle-point (no duality gap) even for non convex problems. Convergence of relaxation algorithms is guaranteed
30 Some backgrounds on price-decomposition Augmented Lagrangian: min J(u 1, u 2 ) sous θ(u 1, u 2 ) =. (37) u=(u 1,u 2 ) L c (u, p) = J(u 1, u 2 )+ < p, θ(u 1, u 2 ) > + c 2 θ(u 1, u 2 ) 2 Duality approach: Computation of a saddle-point, solution of: max min L c (u, p) p u A price-decomposition algorithm (UZAWA): 1) Solve min u L c (u, p k ) u k+1 2) p k+1 = p k + ρθ(u1 k+1, u2 k+1 )
31 If both J(u 1, u 2 ) = J 1 (u 1 ) + J 2 (u 2 ) and θ(u 1, u 2 ) = θ 1 (u 1 ) + θ 2 (u 2 ) (separable case) and if c = (ordinary lagrangian formulation), then price decomposition-coordination: 1. Decomposition: Solve two independant subproblems min ui J i (u i )+ < p k, θ i (u i ) > u k+1 i, i = 1, 2 2. Coordination: p k+1 = p k + ρθ(u1 k+1, u2 k+1 )
32 If θ is not separable and c, separability via linearization of c. 2 : ALGORITHM (Cohen, 84) 1. Decomposition : min J i (u i )+ < p k +cθ(u u 1 k, u2 k ), θ i i(u1 k, u2 k ).u i > + 1 2ɛ u i ui k 2, = u k+1 i, where θ i is the gradient of θ w.r.t. u i. 2. Coordination: p k+1 = p k + ρθ(u k+1 1, u k+1 2 ) with < ρ c et < ɛ < 1/2c (for convex problems).
33 Application to the 2-pool AOCP 1. AOCP of pool 1: min u 1 T T m 1 (u 1 )dt + T L 1 l 1 (z 1 (x, t), Q 1 (x, t))dxdt + [< p k + cf k, F k u 1.u 1 > + 1 2ɛ (u 1 u k 1 ) 2 ]dt s.t. S pool1 + B.C. + I.C.+ T.C., = u1 k+1 (.), 2. AOCP of pool 2: min u 2 T m 2 (u 2 )dt + T L2 L + 1 l i (z 2 (x, t), Q 2 (x, t))dxdt + T p k + cf k, F k u 2.u 2 > + 1 2ɛ (u 2 u k 2 ) 2 ]dt, s.t. S pool1 + B.C. + I.C.+ T.C., = u2 k+1 (.), 3. Coordination: T min [< p k + cf k, F Q k.q > + 1 Q 2ɛ (Q Qk ) 2 ]dt, = Q k+1 (.), p k+1 = p k + ρf (Q k+1, z 1 (u k+1 1 ), z 2 (u k+1 2 ), α(u k+1 2 )) [<
34 Some potential advantages Complexity reduction thanks to the decomposition in sub-problems and parallel computation: In our case study, 2 problems of complexity 4(N + 1) M versus one problem of complexity 8(N + 1) M ; Well suited for distributed control application (Networked Control Systems).
35 A distributed predictive control scheme At each instant dt, parallel computation of the two-point boundary value subproblems POOL 1 Computation of a two-point boundary value problem (similar to a single-pool problem) The local variables z1 and u1 at iteration k are sent through the network POOL 2 Computation of a two-point boundary value problem (similar to a single-pool problem) The local variables z2 and u2 at iteration k are sent through the network NETWORK Coordination: The coordination variables are the Lagrangian multiplier p associated to the gate constraint and the flow rate Q Update of p and Q ; p and Q are sent through the network k+1 -> k Broadcast of variables through the communication network until convergence
36 Conclusions and perspectives 1. Practical extension of finite-dimensional predictive control techniques to an infinite-dimensional problem; 2. Computation of the related two-point boundary value problem using a 2D discretization method based on Preissman integration scheme; 3. Theoretical analysis still to be performed (see M. Herty s presentation); 4. Extension to the multi-pool case via a Lagrangian relaxation technique (Decomposition-Coordination) is possible, but has to be validated; 5. For practical implementation: need of a state observer (can be derived by using variational calculus as a dual control problem ).
Robust boundary control of systems of conservation laws: theoretical study with numerical and experimental validations
1 Robust boundary control of systems of conservation laws: theoretical study with numerical and experimental validations LAAS - CNRS Toulouse (France) Workshop on Irrigation Channels and Related Problems
More informationControl Design of a Distributed Parameter. Fixed-Bed Reactor
Workshop on Control of Distributed Parameter Systems p. 1/1 Control Design of a Distributed Parameter Fixed-Bed Reactor Ilyasse Aksikas and J. Fraser Forbes Department of Chemical and Material Engineering
More informationIntroduction to gradient descent
6-1: Introduction to gradient descent Prof. J.C. Kao, UCLA Introduction to gradient descent Derivation and intuitions Hessian 6-2: Introduction to gradient descent Prof. J.C. Kao, UCLA Introduction Our
More informationHybrid Systems Course Lyapunov stability
Hybrid Systems Course Lyapunov stability OUTLINE Focus: stability of an equilibrium point continuous systems decribed by ordinary differential equations (brief review) hybrid automata OUTLINE Focus: stability
More informationAutomatic Control 2. Nonlinear systems. Prof. Alberto Bemporad. University of Trento. Academic year
Automatic Control 2 Nonlinear systems Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011 1 / 18
More informationSuboptimal Open-loop Control Using POD. Stefan Volkwein
Institute for Mathematics and Scientific Computing University of Graz, Austria PhD program in Mathematics for Technology Catania, May 22, 2007 Motivation Optimal control of evolution problems: min J(y,
More informationStochastic dynamical modeling:
Stochastic dynamical modeling: Structured matrix completion of partially available statistics Armin Zare www-bcf.usc.edu/ arminzar Joint work with: Yongxin Chen Mihailo R. Jovanovic Tryphon T. Georgiou
More informationON LYAPUNOV STABILITY OF LINEARISED SAINT-VENANT EQUATIONS FOR A SLOPING CHANNEL. Georges Bastin. Jean-Michel Coron. Brigitte d Andréa-Novel
NETWORKS AND HETEROGENEOUS MEDIA doi:.3934/nhm.9.4.77 c American Institute of Mathematical Sciences Volume 4, Number, June 9 pp. 77 87 ON LYAPUNOV STABILITY OF LINEARISED SAINT-VENANT EQUATIONS FOR A SLOPING
More informationLyapunov stability ORDINARY DIFFERENTIAL EQUATIONS
Lyapunov stability ORDINARY DIFFERENTIAL EQUATIONS An ordinary differential equation is a mathematical model of a continuous state continuous time system: X = < n state space f: < n! < n vector field (assigns
More informationChapter 2 Optimal Control Problem
Chapter 2 Optimal Control Problem Optimal control of any process can be achieved either in open or closed loop. In the following two chapters we concentrate mainly on the first class. The first chapter
More informationUSING NMPC BASED ON A LOW-ORDER MODEL FOR CONTROLLING PRESSURE DURING OIL WELL DRILLING
8th International IFAC Symposium on Dynamics and Control of Process Systems Preprints Vol.1, June 6-8, 27, Cancún, Mexico USING NMPC BASED ON A LOW-ORDER MODEL FOR CONTROLLING PRESSURE DURING OIL WELL
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationminimize x subject to (x 2)(x 4) u,
Math 6366/6367: Optimization and Variational Methods Sample Preliminary Exam Questions 1. Suppose that f : [, L] R is a C 2 -function with f () on (, L) and that you have explicit formulae for
More informationEN Applied Optimal Control Lecture 8: Dynamic Programming October 10, 2018
EN530.603 Applied Optimal Control Lecture 8: Dynamic Programming October 0, 08 Lecturer: Marin Kobilarov Dynamic Programming (DP) is conerned with the computation of an optimal policy, i.e. an optimal
More informationHybrid Systems - Lecture n. 3 Lyapunov stability
OUTLINE Focus: stability of equilibrium point Hybrid Systems - Lecture n. 3 Lyapunov stability Maria Prandini DEI - Politecnico di Milano E-mail: prandini@elet.polimi.it continuous systems decribed by
More informationDouble Smoothing technique for Convex Optimization Problems with Linear Constraints
1 Double Smoothing technique for Convex Optimization Problems with Linear Constraints O. Devolder (F.R.S.-FNRS Research Fellow), F. Glineur and Y. Nesterov Center for Operations Research and Econometrics
More informationStochastic and Adaptive Optimal Control
Stochastic and Adaptive Optimal Control Robert Stengel Optimal Control and Estimation, MAE 546 Princeton University, 2018! Nonlinear systems with random inputs and perfect measurements! Stochastic neighboring-optimal
More informationFlatness based analysis and control of distributed parameter systems Elgersburg Workshop 2018
Flatness based analysis and control of distributed parameter systems Elgersburg Workshop 2018 Frank Woittennek Institute of Automation and Control Engineering Private University for Health Sciences, Medical
More informationRobotics. Control Theory. Marc Toussaint U Stuttgart
Robotics Control Theory Topics in control theory, optimal control, HJB equation, infinite horizon case, Linear-Quadratic optimal control, Riccati equations (differential, algebraic, discrete-time), controllability,
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 20: LMI/SOS Tools for the Study of Hybrid Systems Stability Concepts There are several classes of problems for
More informationTopic # /31 Feedback Control Systems. Analysis of Nonlinear Systems Lyapunov Stability Analysis
Topic # 16.30/31 Feedback Control Systems Analysis of Nonlinear Systems Lyapunov Stability Analysis Fall 010 16.30/31 Lyapunov Stability Analysis Very general method to prove (or disprove) stability of
More informationMathematical Economics. Lecture Notes (in extracts)
Prof. Dr. Frank Werner Faculty of Mathematics Institute of Mathematical Optimization (IMO) http://math.uni-magdeburg.de/ werner/math-ec-new.html Mathematical Economics Lecture Notes (in extracts) Winter
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 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 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 informationOPTIMAL CONTROL. Sadegh Bolouki. Lecture slides for ECE 515. University of Illinois, Urbana-Champaign. Fall S. Bolouki (UIUC) 1 / 28
OPTIMAL CONTROL Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 28 (Example from Optimal Control Theory, Kirk) Objective: To get from
More informationLinear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013
Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013 Abstract As in optimal control theory, linear quadratic (LQ) differential games (DG) can be solved, even in high dimension,
More informationTheory in Model Predictive Control :" Constraint Satisfaction and Stability!
Theory in Model Predictive Control :" Constraint Satisfaction and Stability Colin Jones, Melanie Zeilinger Automatic Control Laboratory, EPFL Example: Cessna Citation Aircraft Linearized continuous-time
More informationPart 5: Penalty and augmented Lagrangian methods for equality constrained optimization. Nick Gould (RAL)
Part 5: Penalty and augmented Lagrangian methods for equality constrained optimization Nick Gould (RAL) x IR n f(x) subject to c(x) = Part C course on continuoue optimization CONSTRAINED MINIMIZATION x
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationOPTIMAL CONTROL AND ESTIMATION
OPTIMAL CONTROL AND ESTIMATION Robert F. Stengel Department of Mechanical and Aerospace Engineering Princeton University, Princeton, New Jersey DOVER PUBLICATIONS, INC. New York CONTENTS 1. INTRODUCTION
More informationDistributed Receding Horizon Control of Cost Coupled Systems
Distributed Receding Horizon Control of Cost Coupled Systems William B. Dunbar Abstract This paper considers the problem of distributed control of dynamically decoupled systems that are subject to decoupled
More informationStability and boundary stabilization of physical networks represented by 1-D hyperbolic balance laws
Stability and boundary stabilization of physical networks represented by 1-D hyperbolic balance laws G. Bastin http://perso.uclouvain.be/georges.bastin/ Georges Bastin Jean Michel Coron Progress in Nonlinear
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 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 informationEE C128 / ME C134 Feedback Control Systems
EE C128 / ME C134 Feedback Control Systems Lecture Additional Material Introduction to Model Predictive Control Maximilian Balandat Department of Electrical Engineering & Computer Science University of
More informationThe Multi-Path Utility Maximization Problem
The Multi-Path Utility Maximization Problem Xiaojun Lin and Ness B. Shroff School of Electrical and Computer Engineering Purdue University, West Lafayette, IN 47906 {linx,shroff}@ecn.purdue.edu Abstract
More informationModule 07 Controllability and Controller Design of Dynamical LTI Systems
Module 07 Controllability and Controller Design of Dynamical LTI Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha October
More informationProper Orthogonal Decomposition in PDE-Constrained Optimization
Proper Orthogonal Decomposition in PDE-Constrained Optimization K. Kunisch Department of Mathematics and Computational Science University of Graz, Austria jointly with S. Volkwein Dynamic Programming Principle
More informationA State-Space Approach to Control of Interconnected Systems
A State-Space Approach to Control of Interconnected Systems Part II: General Interconnections Cédric Langbort Center for the Mathematics of Information CALIFORNIA INSTITUTE OF TECHNOLOGY clangbort@ist.caltech.edu
More informationSelf-Concordant Barrier Functions for Convex Optimization
Appendix F Self-Concordant Barrier Functions for Convex Optimization F.1 Introduction In this Appendix we present a framework for developing polynomial-time algorithms for the solution of convex optimization
More informationDivide-and-combine Strategies in Statistical Modeling for Massive Data
Divide-and-combine Strategies in Statistical Modeling for Massive Data Liqun Yu Washington University in St. Louis March 30, 2017 Liqun Yu (WUSTL) D&C Statistical Modeling for Massive Data March 30, 2017
More informationMACRODISPERSION AND DISPERSIVE TRANSPORT BY UNSTEADY RIVER FLOW UNDER UNCERTAIN CONDITIONS
MACRODISPERSION AND DISPERSIVE TRANSPORT BY UNSTEADY RIVER FLOW UNDER UNCERTAIN CONDITIONS M.L. Kavvas and L.Liang UCD J.Amorocho Hydraulics Laboratory University of California, Davis, CA 95616, USA Uncertainties
More informationLecture 16: Relaxation methods
Lecture 16: Relaxation methods Clever technique which begins with a first guess of the trajectory across the entire interval Break the interval into M small steps: x 1 =0, x 2,..x M =L Form a grid of points,
More informationNonlinear Optimization: What s important?
Nonlinear Optimization: What s important? Julian Hall 10th May 2012 Convexity: convex problems A local minimizer is a global minimizer A solution of f (x) = 0 (stationary point) is a minimizer A global
More informationSubject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)
Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationCache Oblivious Stencil Computations
Cache Oblivious Stencil Computations S. HUNOLD J. L. TRÄFF F. VERSACI Lectures on High Performance Computing 13 April 2015 F. Versaci (TU Wien) Cache Oblivious Stencil Computations 13 April 2015 1 / 19
More informationProximal splitting methods on convex problems with a quadratic term: Relax!
Proximal splitting methods on convex problems with a quadratic term: Relax! The slides I presented with added comments Laurent Condat GIPSA-lab, Univ. Grenoble Alpes, France Workshop BASP Frontiers, Jan.
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 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 information1. Gradient method. gradient method, first-order methods. quadratic bounds on convex functions. analysis of gradient method
L. Vandenberghe EE236C (Spring 2016) 1. Gradient method gradient method, first-order methods quadratic bounds on convex functions analysis of gradient method 1-1 Approximate course outline First-order
More informationInfeasibility Detection and an Inexact Active-Set Method for Large-Scale Nonlinear Optimization
Infeasibility Detection and an Inexact Active-Set Method for Large-Scale Nonlinear Optimization Frank E. Curtis, Lehigh University involving joint work with James V. Burke, University of Washington Daniel
More informationOn variational data assimilation for 1D and 2D fluvial hydraulics
On variational data assimilation for D and D fluvial hydraulics I. Gejadze, M. Honnorat, FX Le Dimet, and J. Monnier LJK - MOISE project-team, Grenoble, France. Contact: Jerome.Monnier@imag.fr Civil Engineering
More informationA Generalization of Barbalat s Lemma with Applications to Robust Model Predictive Control
A Generalization of Barbalat s Lemma with Applications to Robust Model Predictive Control Fernando A. C. C. Fontes 1 and Lalo Magni 2 1 Officina Mathematica, Departamento de Matemática para a Ciência e
More informationLinear Quadratic Zero-Sum Two-Person Differential Games
Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard To cite this version: Pierre Bernhard. Linear Quadratic Zero-Sum Two-Person Differential Games. Encyclopaedia of Systems and Control,
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationConvergence Rate of Nonlinear Switched Systems
Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the
More informationNonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 71 2009 2744 2752 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A nonlinear inequality and applications N.S. Hoang A.G. Ramm
More informationDifferential Equation Types. Moritz Diehl
Differential Equation Types Moritz Diehl Overview Ordinary Differential Equations (ODE) Differential Algebraic Equations (DAE) Partial Differential Equations (PDE) Delay Differential Equations (DDE) Ordinary
More informationFluid Mechanics Prof. S.K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Fluid Mechanics Prof. S.K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture - 42 Flows with a Free Surface Part II Good morning. I welcome you to this session
More informationEE C128 / ME C134 Final Exam Fall 2014
EE C128 / ME C134 Final Exam Fall 2014 December 19, 2014 Your PRINTED FULL NAME Your STUDENT ID NUMBER Number of additional sheets 1. No computers, no tablets, no connected device (phone etc.) 2. Pocket
More informationIntroduction to Alternating Direction Method of Multipliers
Introduction to Alternating Direction Method of Multipliers Yale Chang Machine Learning Group Meeting September 29, 2016 Yale Chang (Machine Learning Group Meeting) Introduction to Alternating Direction
More informationFast Iterative Solution of Saddle Point Problems
Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, GA Acknowledgments NSF (Computational Mathematics) Maxim Olshanskii (Mech-Math, Moscow State U.) Zhen Wang (PhD student,
More informationAlberto Bressan. Department of Mathematics, Penn State University
Non-cooperative Differential Games A Homotopy Approach Alberto Bressan Department of Mathematics, Penn State University 1 Differential Games d dt x(t) = G(x(t), u 1(t), u 2 (t)), x(0) = y, u i (t) U i
More informationSemidefinite Programming Duality and Linear Time-invariant Systems
Semidefinite Programming Duality and Linear Time-invariant Systems Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 2 July 2004 Workshop on Linear Matrix Inequalities in Control LAAS-CNRS,
More informationGideon Weiss University of Haifa. Joint work with students: Anat Kopzon Yoni Nazarathy. Stanford University, MSE, February, 2009
Optimal Finite Horizon Control of Manufacturing Systems: Fluid Solution by SCLP (separated continuous LP) and Fluid Tracking using IVQs (infinite virtual queues) Stanford University, MSE, February, 29
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 informationPDE and Boundary-Value Problems Winter Term 2014/2015
PDE and Boundary-Value Problems Winter Term 2014/2015 Lecture 13 Saarland University 5. Januar 2015 c Daria Apushkinskaya (UdS) PDE and BVP lecture 13 5. Januar 2015 1 / 35 Purpose of Lesson To interpretate
More informationDynamical Systems & Lyapunov Stability
Dynamical Systems & Lyapunov Stability Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Ordinary Differential Equations Existence & uniqueness Continuous dependence
More informationFlood routing. Prof. (Dr.) Rajib Kumar Bhattacharjya Indian Institute of Technology Guwahati
Flood routing Prof. (Dr.) Rajib Kumar Bhattacharjya Indian Institute of Technology Guwahati Guwahati, Assam Email: rkbc@iitg.ernet.in Web: www.iitg.ernet.in/rkbc Visiting Faculty NIT Meghalaya Q (m 3 /sec)
More informationProblem 1 Cost of an Infinite Horizon LQR
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 5243 INTRODUCTION TO CYBER-PHYSICAL SYSTEMS H O M E W O R K # 5 Ahmad F. Taha October 12, 215 Homework Instructions: 1. Type your solutions in the LATEX homework
More information2 The Linear Quadratic Regulator (LQR)
2 The Linear Quadratic Regulator (LQR) Problem: Compute a state feedback controller u(t) = Kx(t) that stabilizes the closed loop system and minimizes J := 0 x(t) T Qx(t)+u(t) T Ru(t)dt where x and u are
More informationALADIN An Algorithm for Distributed Non-Convex Optimization and Control
ALADIN An Algorithm for Distributed Non-Convex Optimization and Control Boris Houska, Yuning Jiang, Janick Frasch, Rien Quirynen, Dimitris Kouzoupis, Moritz Diehl ShanghaiTech University, University of
More informationInterior Point Algorithms for Constrained Convex Optimization
Interior Point Algorithms for Constrained Convex Optimization Chee Wei Tan CS 8292 : Advanced Topics in Convex Optimization and its Applications Fall 2010 Outline Inequality constrained minimization problems
More informationCHAPTER 3 THE MAXIMUM PRINCIPLE: MIXED INEQUALITY CONSTRAINTS. p. 1/73
CHAPTER 3 THE MAXIMUM PRINCIPLE: MIXED INEQUALITY CONSTRAINTS p. 1/73 THE MAXIMUM PRINCIPLE: MIXED INEQUALITY CONSTRAINTS Mixed Inequality Constraints: Inequality constraints involving control and possibly
More informationLecture Note 7: Switching Stabilization via Control-Lyapunov Function
ECE7850: Hybrid Systems:Theory and Applications Lecture Note 7: Switching Stabilization via Control-Lyapunov Function Wei Zhang Assistant Professor Department of Electrical and Computer Engineering Ohio
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 3. Fundamental properties IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Example Consider the system ẋ = f
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 informationLINEAR AND NONLINEAR PROGRAMMING
LINEAR AND NONLINEAR PROGRAMMING Stephen G. Nash and Ariela Sofer George Mason University The McGraw-Hill Companies, Inc. New York St. Louis San Francisco Auckland Bogota Caracas Lisbon London Madrid Mexico
More informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
More informationConstrained optimization: direct methods (cont.)
Constrained optimization: direct methods (cont.) Jussi Hakanen Post-doctoral researcher jussi.hakanen@jyu.fi Direct methods Also known as methods of feasible directions Idea in a point x h, generate a
More informationInexact Newton Methods and Nonlinear Constrained Optimization
Inexact Newton Methods and Nonlinear Constrained Optimization Frank E. Curtis EPSRC Symposium Capstone Conference Warwick Mathematics Institute July 2, 2009 Outline PDE-Constrained Optimization Newton
More informationLegendre Transforms, Calculus of Varations, and Mechanics Principles
page 437 Appendix C Legendre Transforms, Calculus of Varations, and Mechanics Principles C.1 Legendre Transforms Legendre transforms map functions in a vector space to functions in the dual space. From
More informationmoments of inertia at the center of gravity (C.G.) of the first and second pendulums are I 1 and I 2, respectively. Fig. 1. Double inverted pendulum T
Real-Time Swing-up of Double Inverted Pendulum by Nonlinear Model Predictive Control Pathompong Jaiwat 1 and Toshiyuki Ohtsuka 2 Abstract In this study, the swing-up of a double inverted pendulum is controlled
More informationSwitched systems: stability
Switched systems: stability OUTLINE Switched Systems Stability of Switched Systems OUTLINE Switched Systems Stability of Switched Systems a family of systems SWITCHED SYSTEMS SWITCHED SYSTEMS a family
More informationLINEAR-CONVEX CONTROL AND DUALITY
1 LINEAR-CONVEX CONTROL AND DUALITY R.T. Rockafellar Department of Mathematics, University of Washington Seattle, WA 98195-4350, USA Email: rtr@math.washington.edu R. Goebel 3518 NE 42 St., Seattle, WA
More informationConvex Optimization. Newton s method. ENSAE: Optimisation 1/44
Convex Optimization Newton s method ENSAE: Optimisation 1/44 Unconstrained minimization minimize f(x) f convex, twice continuously differentiable (hence dom f open) we assume optimal value p = inf x f(x)
More informationIntelligent Control. Module I- Neural Networks Lecture 7 Adaptive Learning Rate. Laxmidhar Behera
Intelligent Control Module I- Neural Networks Lecture 7 Adaptive Learning Rate Laxmidhar Behera Department of Electrical Engineering Indian Institute of Technology, Kanpur Recurrent Networks p.1/40 Subjects
More informationNumerical investigations of hillslopes with variably saturated subsurface and overland flows
Numerical investigations of hillslopes with variably saturated subsurface and overland flows ARC DYNAS H. Beaugendre, T. Esclaffer, A. Ern and E. Gaume DYNAS Workshop 06-08/12/04 DYNAS Workshop 06-08/12/04
More informationCompositional Safety Analysis using Barrier Certificates
Compositional Safety Analysis using Barrier Certificates Department of Electronic Systems Aalborg University Denmark November 29, 2011 Objective: To verify that a continuous dynamical system is safe. Definition
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 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 informationMath 273a: Optimization Lagrange Duality
Math 273a: Optimization Lagrange Duality Instructor: Wotao Yin Department of Mathematics, UCLA Winter 2015 online discussions on piazza.com Gradient descent / forward Euler assume function f is proper
More informationPontryagin s maximum principle
Pontryagin s maximum principle Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2012 Emo Todorov (UW) AMATH/CSE 579, Winter 2012 Lecture 5 1 / 9 Pontryagin
More informationFrank-Wolfe Method. Ryan Tibshirani Convex Optimization
Frank-Wolfe Method Ryan Tibshirani Convex Optimization 10-725 Last time: ADMM For the problem min x,z f(x) + g(z) subject to Ax + Bz = c we form augmented Lagrangian (scaled form): L ρ (x, z, w) = f(x)
More informationWorst-Case Complexity Guarantees and Nonconvex Smooth Optimization
Worst-Case Complexity Guarantees and Nonconvex Smooth Optimization Frank E. Curtis, Lehigh University Beyond Convexity Workshop, Oaxaca, Mexico 26 October 2017 Worst-Case Complexity Guarantees and Nonconvex
More informationNumerical Optimization Algorithms
Numerical Optimization Algorithms 1. Overview. Calculus of Variations 3. Linearized Supersonic Flow 4. Steepest Descent 5. Smoothed Steepest Descent Overview 1 Two Main Categories of Optimization Algorithms
More informationControl Systems I. Lecture 2: Modeling and Linearization. Suggested Readings: Åström & Murray Ch Jacopo Tani
Control Systems I Lecture 2: Modeling and Linearization Suggested Readings: Åström & Murray Ch. 2-3 Jacopo Tani Institute for Dynamic Systems and Control D-MAVT ETH Zürich September 28, 2018 J. Tani, E.
More informationEXTREMAL ANALYTICAL CONTROL AND GUIDANCE SOLUTIONS FOR POWERED DESCENT AND PRECISION LANDING. Dilmurat Azimov
EXTREMAL ANALYTICAL CONTROL AND GUIDANCE SOLUTIONS FOR POWERED DESCENT AND PRECISION LANDING Dilmurat Azimov University of Hawaii at Manoa 254 Dole Street, Holmes 22A Phone: (88)-956-2863, E-mail: azimov@hawaii.edu
More information