Optimal control of nonstructured nonlinear descriptor systems
|
|
- Hilary Stanley
- 6 years ago
- Views:
Transcription
1 Optimal control of nonstructured nonlinear descriptor systems TU Berlin DFG Research Center Institut für Mathematik MATHEON Workshop Elgersburg joint work with Peter Kunkel
2 Overview Applications Why DAEs and not ODEs A crash course in DAE theory Optimal control for ODEs Examples where the theory does not work for DAEs Optimality conditions for linear DAEs Differential-algebraic Riccati equations General nonlinear problems Numerical example
3 Optimal control of DAEs Optimal control problem: t J (x, u) = M(x(t)) + K(t, x, u) dt = min! t subject to a DAE constraint x state, u input, y output. F(t, x, u, ẋ) = 0, x(t) = x.
4 Linear quadratic optimal control Cost functional: J (x, u) = 1 2 x(t)t Mx(t) t t (x T Wx + 2x T Su + u T Ru) dt, W = W T C 0 (I, R n,n ), S C 0 (I, R n,l ), R = R T C 0 (I, R l,l ), M = M T R n,n. Constraint: E(t)ẋ = A(t)x + B(t)u + f, x(t) = x, E C 1 (I, R n,n ), A C 0 (I, R n,n ), B C 0 (I, R n,l ), f C 0 (I, R n ), x R n. Here: Determine optimal controls u U = C 0 (I, R l )., More general spaces and also nonsquare E, A are possible.
5 Drop size distributions in stirred liquid/liquid systems with M. Kraume from Chemical Engineering (S. Schlauch)
6 Drop size distributions in stirred liquid/liquid systems with M. Kraume from Chemical Engineering (S. Schlauch) simulation & control
7 Technological Application, Tasks Chemical industry: pearl polymerization and extraction processes Modelling of coalescence and breakage in turbulent flow. Numerical methods for simulation of coupled system of population balance equations/fluid flow equations. Development of optimal control methods for large scale coupled systems Model reduction and observer design. Control of real configurations via stirrer speed. Ultimate goal: Achieve specified average drop diameter and small standard deviation for distribution by real time-control of stirrer-speed.
8 Mathematical system components Navier Stokes equation (flow field) ( Sonjas Film). Population balance equation (drop size distribution). One or two way coupling. Initial and boundary conditions. Space discretization leads to an extremely large control system of nonlinear DAEs.
9 Active flow control with F. Tröltzsch (M. Navier Stokes equations (3D + boundary conditions speaker microphones input output Schmidt)
10 Technological Application, Tasks Control of detached turbulent flow on airline wing Test case (backward step to compare experiment/numerics.) Modelling of turbulent flow. Development of control methods for large scale coupled systems. Model reduction and observer design. Optimal control of real configurations via blowing and sucking of air in wing. Ultimate goal: Force detached flow back to wing.
11 Simulated flow
12 Controlled flow Movement of recirculation bubble following reference curve. 9 8 x ref (t) x r (t) time t 4 2 u(t) time t
13 Model based control of high-tech automatic gearboxes with P. Hamann Daimler/Chrysler ( DC Film)
14 Technological Application, Tasks Decrease full consumption, improve smoothness of switching Modelling of multi-model system, including elasticity, hydraulics Development of control methods for coupled system. Model reduction and observer design. Real time optimal control of real configuration. Ultimate goal: Increase efficiency, reduce costs Space discretization leads to a large hybrid, switched control system of nonlinear DAEs.
15 DAE control systems After space discretization all these problems are DAE control systems or in the linear case F(t, x, ẋ, u) = 0, G(t, x, u) = y, E(t)ẋ(t) = A(t)x(t) + B(t)u(t) + f (t), y(t) = C(t)x(t) + D(t)u(t) + g(t). Using a behavior approach, i.e., forming z(t) = (x, u, y) we obtain general non-square DAEs F(t, z, ż) = 0, E(t)ż = A(t)z. The behavior approach allows a uniform mathematical treatment of simulation and control problems!
16 Why DAEs and not ODEs? DAEs provide a unified framework for the analysis, simulation and control of coupled dynamical systems (continuous and discrete time). Automatic modelling leads to DAEs. (Constraints at interfaces). Conservation laws lead to DAEs. (Conservation of mass, energy, momentum). Coupling of solvers leads to DAEs (discrete time). Control problems are DAEs (behavior).
17 How does one solve such complex problems today? Simplified models. Space discretization with very coarse meshes. Identification and realization of black box models. Model reduction (mostly based on heuristic methods). Coupling of simulation packages. Use of standard optimal control techniques for simplified mathematical model. But do they work for these models?
18 Is there anything to do? Why not just apply the Pontryagin maximum principle? For linear ODEs the initial value problem has a unique solution x C 1 (I, R n ) for every u U, every f C 0 (I, R n ), and every initial value x R n. DAEs, where E(t) may be singular, may not be (uniquely) solvable for any u U and also the initial conditons may be restricted. Furthermore, we need solutions x X, where X usually is a larger space than C 1 (I, R n ).
19 Previous work Linear constant coefficient index 1 case, Bender/Laub 87, Campbell 87 M. 91, Geerts 93. Regularization to index 1, Bunse-Gerstner/M./Nichols 92, 94, Byers/Geerts/M. 97, Byers/Kunkel/M. 97. Linear variable coefficients index 1 case, Kunkel./M. 97. Semi-explicit nonlinear index 1 case, maximum principle, De Pinho/Vinter 97, Devdariani/Ledyaev 99. Semi-explicit index 2, 3 case Roubicek/Valasek 02. Linear index 1, 2 case with properly stated leading term, Balla/März, 02,04, Balla/Linh 05, Kurina/März 04, Backes 06. Multibody systems (structured index 3), Büskens/Gerdts 00, Gerdts 03,06.
20 A crash course in DAE Theory For numerical methods and for the design of controllers, we use derivative arrays (Campbell 1989). We assume that derivatives of original functions are available or can be obtained via computer algebra or automatic differentiation. Linear case: We put E(t)ẋ = A(t)x + f (t) and its derivatives up to order µ into a large DAE for z k = (x, ẋ,..., x (k) ). M 2 = M k (t)ż k = N k (t)z k + g k (t), k N 0 E 0 0 A Ė E 0 Ȧ 2Ë A Ė E, N 2 = A 0 0 Ȧ 0 0 Ä 0 0, z 2 = ẋ x ẍ.
21 Theorem, Kunkel/M Under some constant rank assumptions, for a linear DAE there exist integers µ, a, d and v such that: 1. corank M µ+1 (t) corank M µ (t) = v. 2. rank M µ (t) = (µ + 1)m a v on I, and there exists a smooth matrix function Z 2,3 (left nullspace of M µ ) with Z T 2,3 M µ(t) = The projection Z 2,3 can be partitioned into two parts: Z 2 (left nullspace of (M µ, N µ )) so that the first block column Â2 of Z 2 N µ(t) has full rank a and Z 3 N µ(t) = 0. Let T 2 be a smooth matrix function such that Â2T 2 = 0, (right nullspace of Â2). 4. rank E(t)T 2 = d = l a v and there exists a smooth matrix function Z 1 of size (n, d) with rank Ê1 = d, where Ê1 = Z T 1 E.
22 Reduced problem The quantity µ is called the strangeness-index. It describes the smoothness requirements for forcing or input functions. It generalizes the usual differentiation index to general DAEs (and counts slightly differently). We obtain a numerically computable condensed form Ê 1 (t)ẋ = Â1(t)x + ˆf 1 (t), ddifferential equations 0 = Â2(t)x + ˆf 2 (t), a algebraic equations 0 = ˆf3 (t), v consistency equations where Â1 = Z T 1 A, ˆf 1 = Z T 1 f, and ˆf 2 = Z T 2 g µ, ˆf 3 = Z T 3 g µ. The reduced system has the same solution set as the orignal problem but now it has strangeness-index 0. Remodeling!
23 Calculus of variations for ODEs (E=I) Introduce Lagrange multiplier function λ(t) and couple constraint into cost function, i.e. minimize J (x, u, λ) = 1 2 x(t)t Mx(t) t + λ T (ẋ Ax + Bu + f ) dt. t (x T Wx + 2x T Su + u T Ru) Consider x + δx, u + δu and λ + δλ. For a minimum the cost function has to go up in the neighborhood, so we get optimality conditions (Euler-Lagrange equations):
24 Optimality system Theorem If (x, u) is a solution to the optimal control problem, then there exists a Lagrange multiplier function λ C 1 (I, R n ), such that (x, λ, u) satisfy the optimality boundary value problem (a) ẋ = Ax + Bu + f, x(t) = x, (b) λ = Wx + Su A T λ, λ(t) = Mx(t), (c) 0 = S T x + Ru B T λ.
25 Naive Idea for DAEs Replace the identity in front of x by E(t) and then do the analysis in the same way. For DAEs the formal optimality system then could be (a) Eẋ = Ax + Bu + f, x(t) = x d (b) dt (E T λ) = Wx + Su A T λ, (E T λ)(t) = Mx(t), (b) 0 = S T x + Ru B T λ. This works if the system has strangeness-index µ = 0 as a free system with u = 0 but not in general.
26 What are the difficulties In the proof one has to guarantee that the resulting adjoint equation for λ has a unique solution. One needs a density argument in the solution space. But, the formal adjoint equation may not have a (unique) solution. The formal boundary conditions may not be consistent. The solution of the optimality systems may not exist or may not be unique.
27 Example Consider J (x, u) = (x u 2 )dt = min! subject to the differential-algebraic system [ ] [ ] [ ] [ ] [ 0 1 ẋ1 1 0 x1 1 = ẋ2 0 1 x 2 0 A simple calculation yields the optimal solution ] [ ] f1 u +. f 2 x 1 = u = λ 1 = 1 2 (f 1 + ḟ2), x 2 = f 2, λ 2 = 0.
28 In the formal optimality system we get x 1 = u = λ 1 = 1 2 (f 1 + ḟ2), x 2 = f 2, λ 2 = 1 2 (ḟ1 + f 2 ) without using the initial condition λ 1 (1) = 0. Depending on the data, this initial condition may be consistent or not. This initial condition should not be present. Moreover, λ 2 requires more smoothness of the inhomogeneity than in the optimal solution. Further examples, see Dissertation Backes 06.
29 Solution space To derive optimality conditions for DAEs, we need the right solution space for x: X = C 1 E + E(I, R n ) = { x C 0 (I, R n ) E + Ex C 1 (I, R n ) }, where E + denotes the Moore-Penrose inverse of the matrix valued function, which is the unique matrix function that satisfies the Penrose axioms. EE + E = E, E + EE + = E +, (EE + ) T = EE +, (E + E) T = E + E U is usually a set of piecewise continuous functions or even distributions.
30 Necessary optimality condition Theorem Consider the linear quadratic optimal control problem with a consistent initial condition. Suppose that the system has µ = 0 as a behavior system and that Mx(t) cokernel E(t). If (x, u) X U is a solution to this optimal control problem, then there exists a Lagrange multiplier function λ C 1 E + E (I, Rn ), such that (x, λ, u) satisfy the optimality boundary value problem E d dt (E + Ex) = (A + E d dt (E + E))x + Bu + f, (E + Ex)(t) = x, E T d dt (EE + λ) = Wx + Su (A + EE +Ė)T λ, (EE + λ)(t) = E + (t) T Mx(t), 0 = S T x + Ru B T λ.
31 Sufficient condition For this we need some kind of convexity. Theorem Consider the optimal control problem with a consistent initial condition and suppose that in the cost functional we have that [ ] W S S T, M R are (pointwise) positive semidefinite. If (x, u, λ) satisfies the formal optimality system then for any (x, u) satisfying the constraint we have J (x, u) J (x, u ).
32 Remarks If a minimum exists, then it satisfies the optimality system. If a unique solution to the formal optimality system exists, then x, u are the same as from the optimality system, λ may be different. To get a sufficient condition, we need some extra convexity. The optimality DAE may have µ > 0. Then it is numerically difficult to solve and further consistency conditions or smoothness requirements arise. The condition that the original system has µ = 0 as a behavior system is not necessary if the cost function is chosen appropriately, so that the resulting optimality system has µ = 0.
33 Differential-algebraic Riccati equations If R in the cost functional is invertible, and if the system has µ = 0 as a free system with u = 0, then one can (at least in theory) apply the usual Riccati approach to E d dt (E + Ex) = (A + E d dt (E + E))x + Bu + f, (E + Ex)(t) = x, E T d dt (EE + λ) = Wx + Su (A + EE +Ė)T λ, (EE + λ)(t) = E + (t) T Mx(t), 0 = S T x + Ru B T λ. If µ > 0 or R is singular, then the Riccati approach may not work, even if the boundary value problem has a unique solution.
34 Nonlinear problems For nonlinear systems F(t, x, ẋ) = 0 one considers nonlinear derivative arrays: We set 0 = F k (t, x, ẋ,..., x (k+1) ) = F (t, x, ẋ) F(t, x, ẋ)... d dt d k dt k F (t, x, ẋ). M k (t, x, ẋ,..., x (k+1) ) = F k;ẋ,...,x (k+1)(t, x, ẋ,..., x (k+1) ), N k (t, x, ẋ,..., x (k+1) ) = (F k;x (t, x, ẋ,..., x (k+1) ), 0,..., 0), z k = (t, x, ẋ,..., x (k+1) ).
35 Hypothesis: There exist integers µ, r, a, d, and v such that L = Fµ 1 ({0}). We have rank F µ;t,x,ẋ,...,x (µ+1) = rank F µ;x,ẋ,...,x (µ+1) = r, in a neighborhood of L such that there exists an equivalent system F(z µ ) = 0 with a Jacobian of full row rank r. On L we have 1. corank F µ;x,ẋ,...,x (µ+1) corank F µ 1;x,ẋ,...,x (µ+1) = v. 2. corank F x,ẋ,...,x (µ+1) = a and there exist smooth matrix functions Z 2 (left nullspace of M µ ) and T 2 (right nullspace of Â2 = F x ) with Z2 T F x,ẋ,...,x (µ+1) = 0 and Z2 T Â2T 2 = rank FẋT 2 = d, d = m a v, and there exists a smooth matrix function Z 1 with rank Z1 T F ẋ = d.
36 Theorem Kunkel/M The solution set L forms a (smooth) manifold of dimension (µ + 2)n + 1 r. The DAE can locally be transformed (by application of the implicit function theorem) to a reduced DAE of the form ẋ 1 = G 1 (t, x 1, x 3 ), (d differential equations), x 2 = G 2 (t, x 1, x 3 ), (a algebraic equations), 0 = 0 (v redundant equations). The variables x 3 represent undetermined components (controls).
37 Optimality conditions Assume that µ = 0 for the system in behavior form with z = (x, u), then in terms of the reduced DAE, the local optimality system is (a) ẋ 1 = L(t, x 1, u), x 1 (t) = x 1, (b) x 2 = R(t, x 1, u), (c) λ1 = K x1 (t, x 1, x 2, u) T L x1 (t, x 1, x 2, u) T λ 1 R x1 (t, x 1, u) T λ 1, λ 1 (t) = M x1 (x 1 (t), x 2 (t)) T (d) 0 = K x2 (t, x 1, x 2, u) T + λ 2, (e) 0 = K u (t, x 1, x 2, u) T L u (t, x 1, u) T λ 1 R u (t, x 1, u) T λ 2, (f) γ = λ 1 (t) Here λ 1, λ 2 are Lagrange multipliers associated with x 1, x 2 and γ is associated with the initial value constraint.
38 Remarks These are local results. All the results can be generalized to general nonsquare nonlinear systems. End point conditions for x can be included. Input and state constraints can be included to give a maximum principle.
39 Numerical Methods The linear optimality system that we would have to solve numerically has the form (a) Ê 1 ẋ = Â1x + ˆB 1 u + ˆf 1, (Ê + 1 Ê1x)(t) = x (b) 0 = Â2x + ˆB 2 u + ˆf 2, d (c) dt (Ê 1 T λ 1) = Wx + Su ÂT 1 λ 1 ÂT 2 λ 2, λ 1 (t) = [ Ê + 1 (t)t 0 ]Mx(t), (d) 0 = S T x + Ru ˆB 1 T λ 1 ˆB 2 T λ 2. where E i, A i, B i are obtained by projection with smooth orthogonal projections Z i from the derivative array. An analogous structure arises locally in the nonlinear case.
40 Numerical Problems In the implementation of the numerical integration codes we use nonsmooth projectors Z T 1, Z T 2. It would be too expensive to carry smooth projectors along. For numerical forward (in time) simulation, it is enough that we know the existence of smooth projectors. Integration methods like Runge-Kutta or BDF don t see the nonsmooth behavior. But the adjoint variables (Lagrange multipliers) depend on these projections and their derivatives. However, even if Z T 1, Z T 2 are nonsmooth, Z 1Z T 1 and Z 2Z T 2 are smooth.
41 Smooth optimality system Choose Ê T 1 λ 1 = E T Z 1 λ 1 = E T Z 1 Z T 1 Z 1 λ 1 = E T Z 1 Z T 1 ˆλ 1. With ˆλ 1 = Z 1 λ 1 we obtain smooth coefficients for ˆλ 1. However, we have to add the condition that ˆλ 1 range Z 1 to the system. If Z i completes Z i to a full orthogonal matrix (we compute these anyway when doing a QR or SVD computation) then these conditions can be expressed as Z i T ˆλi = 0, i = 1, 2
42 New linear optimality system For the numerical solution we use the optimality system. (a) Ê 1 ẋ = Â1x + ˆB 1 u + ˆf 1, (Ê + 1 Ê1x)(t) = x, (b) 0 = Â2x + ˆB 2 u + ˆf 2, (c) d dt (E T Z 1 Z1 T ˆλ 1 ) = Wx + Su A T ˆλ1 [ I n ]Nµ T ˆλ 2, (Z1 T ˆλ 1 )(t) = [ Ê + 1 (t)t 0 ]Mx(t), (d) 0 = S T x + Ru B T ˆλ1 [ 0 I l ]Nµ T ˆλ 2 (e) 0 = Z 1 T ˆλ1, (f) 0 = Z 2 T ˆλ 2. All quantities are available for all time steps. An analoguous system can be derived for each Gauss-Newton step in the nonlinear case.
43 Numerical Example A motor controlled pendulum with a motor in the origin shall be driven into its equilibrium with minimal costs, ex. from Büskens/Gerdts y m 3 l x mg
44 Model problem The model problem for this reads J(x, u) = 3 0 u(t) 2 dt = min! s.t. ẋ 1 = x 3, x 1 (0) = 2 1 2, g = 9.81 ẋ 2 = x 4, x 2 (0) = 2 1 2, ẋ 3 = 2x 1 x 5 + x 2 u, x 3 (0) = 0, ẋ 4 = g 2x 2 x 5 x 1 u, x 4 (0) = 0, 0 = x1 2 + x 2 2 1, x 5(0) = 1 2 gx 2(0).
45 Properties The DAE satisfies the Hypothesis with µ = 2, a = 3, d = 2, and v = 0. Hence, only two scalar initial values are sufficient to describe the initial state. We chose them to be x 2 (0) = and x3 (0) = 0. Similarly, x 1 (3) = 0 and x 3 (0) = 0 are sufficient to describe the equilibrium at the end point. As initial trajectory we took x 1 (t) = t, x3 (t) = 0, x 2 (t) = 1 x 1 (t) 2, x 4 (t) = 0, x 5 (t) = 1 2 gx 2(t). with all other unknowns set to zero on a equidistant grid of 60 intervals.
46 Gauss-Newton method Tolerance for the Gauß-Newton method was Let k count the iterations and w k 2 denote the Euclidian norm of the corresponding Gauß-Newton correction. k w k D D D D D-11 Initial bad convergence is due to a bad initial guess. Final value of cost function is J opt = 3.82 which is correct up to discretization and roundoff errors.
47 What has been achieved? Theoretical analysis (solvability) for general over and under-determined linear and nonlinear DAEs of arbitrary index (blue book) Optimal control and maximum principle. Model verification, model reduction and removal of redundancies is possible in a numerically stable way. General nonlinear (regular) DAE boundary value problem solver, Kunkel/M./Stöver 2004.
48 Thank you very much for your attention.
Modelling, Simulation and Control in Key Technologies
Modelling, Simulation and Control in Key Technologies Volker Mehrmann TU Berlin, Institut für Mathematik DFG Research Center MATHEON Mathematics for key technologies 09.02.11, Delft Overview 1 Introduction
More informationNumerical Treatment of Unstructured. Differential-Algebraic Equations. with Arbitrary Index
Numerical Treatment of Unstructured Differential-Algebraic Equations with Arbitrary Index Peter Kunkel (Leipzig) SDS2003, Bari-Monopoli, 22. 25.06.2003 Outline Numerical Treatment of Unstructured Differential-Algebraic
More informationDifferential-algebraic equations. Control and Numerics I
Differential-algebraic equations. Control and Numerics I Volker Mehrmann Institut für Mathematik Technische Universität Berlin P. Kunkel DFG Research Center MATHEON Mathematics for key technologies 8.3.10
More informationStability analysis of differential algebraic systems
Stability analysis of differential algebraic systems Volker Mehrmann TU Berlin, Institut für Mathematik with Vu Hoang Linh and Erik Van Vleck DFG Research Center MATHEON Mathematics for key technologies
More informationNUMERICAL SOLUTION OF HYBRID SYSTEMS OF DIFFERENTIAL-ALGEBRAIC EQUATIONS
NUMERICAL SOLUTION OF HYBRID SYSTEMS OF DIFFERENTIAL-ALGEBRAIC EQUATIONS PETER HAMANN AND VOLKER MEHRMANN Abstract. We present a mathematical framework for general over- and underdetermined hybrid (switched)
More informationDefect-based a-posteriori error estimation for implicit ODEs and DAEs
1 / 24 Defect-based a-posteriori error estimation for implicit ODEs and DAEs W. Auzinger Institute for Analysis and Scientific Computing Vienna University of Technology Workshop on Innovative Integrators
More informationBohl exponent for time-varying linear differential-algebraic equations
Bohl exponent for time-varying linear differential-algebraic equations Thomas Berger Institute of Mathematics, Ilmenau University of Technology, Weimarer Straße 25, 98693 Ilmenau, Germany, thomas.berger@tu-ilmenau.de.
More informationOn linear quadratic optimal control of linear time-varying singular systems
On linear quadratic optimal control of linear time-varying singular systems Chi-Jo Wang Department of Electrical Engineering Southern Taiwan University of Technology 1 Nan-Tai Street, Yungkung, Tainan
More informationModeling & Simulation 2018 Lecture 12. Simulations
Modeling & Simulation 2018 Lecture 12. Simulations Claudio Altafini Automatic Control, ISY Linköping University, Sweden Summary of lecture 7-11 1 / 32 Models of complex systems physical interconnections,
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationTechnische Universität Berlin Institut für Mathematik
Technische Universität Berlin Institut für Mathematik A flow-on-manifold formulation of differential-algebraic equations Ann-Kristin Baum Preprint 24-37 Preprint-Reihe des Instituts für Mathematik Technische
More informationObservability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u)
Observability Dynamic Systems Lecture 2 Observability Continuous time model: Discrete time model: ẋ(t) = f (x(t), u(t)), y(t) = h(x(t), u(t)) x(t + 1) = f (x(t), u(t)), y(t) = h(x(t)) Reglerteknik, ISY,
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 3 3.4 Differential Algebraic Systems 3.5 Integration of Differential Equations 1 Outline 3.4 Differential Algebraic Systems 3.4.1 Constrained Dynamics 3.4.2 First and Second
More informationCHAPTER 10: Numerical Methods for DAEs
CHAPTER 10: Numerical Methods for DAEs Numerical approaches for the solution of DAEs divide roughly into two classes: 1. direct discretization 2. reformulation (index reduction) plus discretization Direct
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 informationNumerical Optimal Control Overview. Moritz Diehl
Numerical Optimal Control Overview Moritz Diehl Simplified Optimal Control Problem in ODE path constraints h(x, u) 0 initial value x0 states x(t) terminal constraint r(x(t )) 0 controls u(t) 0 t T minimize
More informationTheory and Applications of Constrained Optimal Control Proble
Theory and Applications of Constrained Optimal Control Problems with Delays PART 1 : Mixed Control State Constraints Helmut Maurer 1, Laurenz Göllmann 2 1 Institut für Numerische und Angewandte Mathematik,
More informationOptimal Linear Feedback Control for Incompressible Fluid Flow
Optimal Linear Feedback Control for Incompressible Fluid Flow Miroslav K. Stoyanov Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of the
More informationLecture 6. Regularized least-squares and minimum-norm methods 6 1
Regularized least-squares and minimum-norm methods 6 1 Lecture 6 Regularized least-squares and minimum-norm methods EE263 Autumn 2004 multi-objective least-squares regularized least-squares nonlinear least-squares
More informationME751 Advanced Computational Multibody Dynamics
ME751 Advanced Computational Multibody Dynamics Review: Elements of Linear Algebra & Calculus September 9, 2016 Dan Negrut University of Wisconsin-Madison Quote of the day If you can't convince them, confuse
More informationOptimal control problems with PDE constraints
Optimal control problems with PDE constraints Maya Neytcheva CIM, October 2017 General framework Unconstrained optimization problems min f (q) q x R n (real vector) and f : R n R is a smooth function.
More informationThe Additional Dynamics of the Least Squares Completions of Linear Differential Algebraic Equations. by Irfan Okay
ABSTRACT OKAY, IRFAN The Additional Dynamics of the Least Squares Completions of Linear Differential Algebraic Equations ( Under the direction of Dr Stephen L Campbell ) Differential equations of the form
More information10.34: Numerical Methods Applied to Chemical Engineering. Lecture 19: Differential Algebraic Equations
10.34: Numerical Methods Applied to Chemical Engineering Lecture 19: Differential Algebraic Equations 1 Recap Differential algebraic equations Semi-explicit Fully implicit Simulation via backward difference
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 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 informationLinear Algebra Review (Course Notes for Math 308H - Spring 2016)
Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Dr. Michael S. Pilant February 12, 2016 1 Background: We begin with one of the most fundamental notions in R 2, distance. Letting (x 1,
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More informationA Gauss Lobatto quadrature method for solving optimal control problems
ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes
More informationA relaxation of the strangeness index
echnical report from Automatic Control at Linköpings universitet A relaxation of the strangeness index Henrik idefelt, orkel Glad Division of Automatic Control E-mail: tidefelt@isy.liu.se, torkel@isy.liu.se
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationOrdinary Differential Equation Theory
Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit
More informationDelay Differential-Algebraic Equations
ECHNISCHE UNIVERSIÄ BERLIN Analysis and Numerical Solution of Linear Delay Differential-Algebraic Equations Ha Phi and Volker Mehrmann Preprint 214/42 Preprint-Reihe des Instituts für Mathematik echnische
More informationThe Direct Transcription Method For Optimal Control. Part 2: Optimal Control
The Direct Transcription Method For Optimal Control Part 2: Optimal Control John T Betts Partner, Applied Mathematical Analysis, LLC 1 Fundamental Principle of Transcription Methods Transcription Method
More informationOn Linear-Quadratic Control Theory of Implicit Difference Equations
On Linear-Quadratic Control Theory of Implicit Difference Equations Daniel Bankmann Technische Universität Berlin 10. Elgersburg Workshop 2016 February 9, 2016 D. Bankmann (TU Berlin) Control of IDEs February
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 informationThe Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications
MAX PLANCK INSTITUTE Elgersburg Workshop Elgersburg February 11-14, 2013 The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications Timo Reis 1 Matthias Voigt 2 1 Department
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS
More informationSteady State Kalman Filter
Steady State Kalman Filter Infinite Horizon LQ Control: ẋ = Ax + Bu R positive definite, Q = Q T 2Q 1 2. (A, B) stabilizable, (A, Q 1 2) detectable. Solve for the positive (semi-) definite P in the ARE:
More informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
More informationIterative Methods for Solving A x = b
Iterative Methods for Solving A x = b A good (free) online source for iterative methods for solving A x = b is given in the description of a set of iterative solvers called templates found at netlib: http
More informationControl, Stabilization and Numerics for Partial Differential Equations
Paris-Sud, Orsay, December 06 Control, Stabilization and Numerics for Partial Differential Equations Enrique Zuazua Universidad Autónoma 28049 Madrid, Spain enrique.zuazua@uam.es http://www.uam.es/enrique.zuazua
More informationCS 257: Numerical Methods
CS 57: Numerical Methods Final Exam Study Guide Version 1.00 Created by Charles Feng http://www.fenguin.net CS 57: Numerical Methods Final Exam Study Guide 1 Contents 1 Introductory Matter 3 1.1 Calculus
More informationNumerical Methods I Solving Square Linear Systems: GEM and LU factorization
Numerical Methods I Solving Square Linear Systems: GEM and LU factorization Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 18th,
More informationVerification of analog and mixed-signal circuits using hybrid systems techniques
FMCAD, November 2004, Austin Verification of analog and mixed-signal circuits using hybrid systems techniques Thao Dang, Alexandre Donze, Oded Maler VERIMAG Grenoble, France Plan 1. Introduction 2. Verification
More informationLinear Algebra and Robot Modeling
Linear Algebra and Robot Modeling Nathan Ratliff Abstract Linear algebra is fundamental to robot modeling, control, and optimization. This document reviews some of the basic kinematic equations and uses
More informationModel order reduction of electrical circuits with nonlinear elements
Model order reduction of electrical circuits with nonlinear elements Andreas Steinbrecher and Tatjana Stykel 1 Introduction The efficient and robust numerical simulation of electrical circuits plays a
More informationLinear Multistep Methods I: Adams and BDF Methods
Linear Multistep Methods I: Adams and BDF Methods Varun Shankar January 1, 016 1 Introduction In our review of 5610 material, we have discussed polynomial interpolation and its application to generating
More informationTrajectory Planning from Multibody System Dynamics
Trajectory Planning from Multibody System Dynamics Pierangelo Masarati Politecnico di Milano Dipartimento di Ingegneria Aerospaziale Manipulators 2 Manipulator: chain of
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 informationLecture 6. Numerical methods. Approximation of functions
Lecture 6 Numerical methods Approximation of functions Lecture 6 OUTLINE 1. Approximation and interpolation 2. Least-square method basis functions design matrix residual weighted least squares normal equation
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
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 informationNonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1
Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems p. 1/1 p. 2/1 Converse Lyapunov Theorem Exponential Stability Let x = 0 be an exponentially stable equilibrium
More informationRegularity and approximations of generalized equations; applications in optimal control
SWM ORCOS Operations Research and Control Systems Regularity and approximations of generalized equations; applications in optimal control Vladimir M. Veliov (Based on joint works with A. Dontchev, M. Krastanov,
More informationDirect Methods. Moritz Diehl. Optimization in Engineering Center (OPTEC) and Electrical Engineering Department (ESAT) K.U.
Direct Methods Moritz Diehl Optimization in Engineering Center (OPTEC) and Electrical Engineering Department (ESAT) K.U. Leuven Belgium Overview Direct Single Shooting Direct Collocation Direct Multiple
More informationTransient Sensitivity Analysis CASA Day 13th Nov 2007 Zoran Ilievski. Zoran Ilievski Transient Sensitivity Analysis
CASA Day 13th Nov 2007 Talk Structure Talk Structure Introduction Talk Structure Introduction Recap Sensitivity Talk Structure Introduction Recap Sensitivity Examples and Results Talk Structure Introduction
More informationENGI9496 Lecture Notes State-Space Equation Generation
ENGI9496 Lecture Notes State-Space Equation Generation. State Equations and Variables - Definitions The end goal of model formulation is to simulate a system s behaviour on a computer. A set of coherent
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 informationNumerical Optimal Control Part 3: Function space methods
Numerical Optimal Control Part 3: Function space methods SADCO Summer School and Workshop on Optimal and Model Predictive Control OMPC 2013, Bayreuth Institute of Mathematics and Applied Computing Department
More informationSuppose that we have a specific single stage dynamic system governed by the following equation:
Dynamic Optimisation Discrete Dynamic Systems A single stage example Suppose that we have a specific single stage dynamic system governed by the following equation: x 1 = ax 0 + bu 0, x 0 = x i (1) where
More informationDifferential-algebraic equations. Control and Numerics V
Differential-algebraic equations. Control and Numerics V Volker Mehrmann Institut für Mathematik Technische Universität Berlin with V.H. Linh, P. Kunkel, E. Van Vleck DFG Research Center MATHEON Mathematics
More informationConvergence of a Gauss Pseudospectral Method for Optimal Control
Convergence of a Gauss Pseudospectral Method for Optimal Control Hongyan Hou William W. Hager Anil V. Rao A convergence theory is presented for approximations of continuous-time optimal control problems
More informationLecture 4 Continuous time linear quadratic regulator
EE363 Winter 2008-09 Lecture 4 Continuous time linear quadratic regulator continuous-time LQR problem dynamic programming solution Hamiltonian system and two point boundary value problem infinite horizon
More informationFinal Examination. CS 205A: Mathematical Methods for Robotics, Vision, and Graphics (Fall 2013), Stanford University
Final Examination CS 205A: Mathematical Methods for Robotics, Vision, and Graphics (Fall 2013), Stanford University The exam runs for 3 hours. The exam contains eight problems. You must complete the first
More information1.2 Derivation. d p f = d p f(x(p)) = x fd p x (= f x x p ). (1) Second, g x x p + g p = 0. d p f = f x g 1. The expression f x gx
PDE-constrained optimization and the adjoint method Andrew M. Bradley November 16, 21 PDE-constrained optimization and the adjoint method for solving these and related problems appear in a wide range of
More informationLösning: Tenta Numerical Analysis för D, L. FMN011,
Lösning: Tenta Numerical Analysis för D, L. FMN011, 090527 This exam starts at 8:00 and ends at 12:00. To get a passing grade for the course you need 35 points in this exam and an accumulated total (this
More informationOptimal Control. McGill COMP 765 Oct 3 rd, 2017
Optimal Control McGill COMP 765 Oct 3 rd, 2017 Classical Control Quiz Question 1: Can a PID controller be used to balance an inverted pendulum: A) That starts upright? B) That must be swung-up (perhaps
More informationNumerical Methods for Embedded Optimization and Optimal Control. Exercises
Summer Course Numerical Methods for Embedded Optimization and Optimal Control Exercises Moritz Diehl, Daniel Axehill and Lars Eriksson June 2011 Introduction This collection of exercises is intended to
More informationNumerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationStabilität differential-algebraischer Systeme
Stabilität differential-algebraischer Systeme Caren Tischendorf, Universität zu Köln Elgersburger Arbeitstagung, 11.-14. Februar 2008 Tischendorf (Univ. zu Köln) Stabilität von DAEs Elgersburg, 11.-14.02.2008
More informationApproximation of spectral intervals and associated leading directions for linear DAEs via smooth SVDs*
Approximation of spectral intervals and associated leading directions for linear DAEs via smooth SVDs* (supported by Alexander von Humboldt Foundation) Faculty of Mathematics, Mechanics and Informatics
More informationOPTIMAL FUSION OF SENSOR DATA FOR DISCRETE KALMAN FILTERING Z. G. FENG, K. L. TEO, N. U. AHMED, Y. ZHAO, AND W. Y. YAN
Dynamic Systems and Applications 16 (2007) 393-406 OPTIMAL FUSION OF SENSOR DATA FOR DISCRETE KALMAN FILTERING Z. G. FENG, K. L. TEO, N. U. AHMED, Y. ZHAO, AND W. Y. YAN College of Mathematics and Computer
More informationEnsemble prediction and strategies for initialization: Tangent Linear and Adjoint Models, Singular Vectors, Lyapunov vectors
Ensemble prediction and strategies for initialization: Tangent Linear and Adjoint Models, Singular Vectors, Lyapunov vectors Eugenia Kalnay Lecture 2 Alghero, May 2008 Elements of Ensemble Forecasting
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationThe Plan. Initial value problems (ivps) for ordinary differential equations (odes) Review of basic methods You are here: Hamiltonian systems
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit V Solution of Differential Equations Part 2 Dianne P. O Leary c 2008 The Plan
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationUpon successful completion of MATH 220, the student will be able to:
MATH 220 Matrices Upon successful completion of MATH 220, the student will be able to: 1. Identify a system of linear equations (or linear system) and describe its solution set 2. Write down the coefficient
More informationRobotics 1 Inverse kinematics
Robotics 1 Inverse kinematics Prof. Alessandro De Luca Robotics 1 1 Inverse kinematics what are we looking for? direct kinematics is always unique; how about inverse kinematics for this 6R robot? Robotics
More informationA Tutorial on Recursive methods in Linear Least Squares Problems
A Tutorial on Recursive methods in Linear Least Squares Problems by Arvind Yedla 1 Introduction This tutorial motivates the use of Recursive Methods in Linear Least Squares problems, specifically Recursive
More informationFormula Sheet for Optimal Control
Formula Sheet for Optimal Control Division of Optimization and Systems Theory Royal Institute of Technology 144 Stockholm, Sweden 23 December 1, 29 1 Dynamic Programming 11 Discrete Dynamic Programming
More informationNumerical Integration of Equations of Motion
GraSMech course 2009-2010 Computer-aided analysis of rigid and flexible multibody systems Numerical Integration of Equations of Motion Prof. Olivier Verlinden (FPMs) Olivier.Verlinden@fpms.ac.be Prof.
More informationECEN 615 Methods of Electric Power Systems Analysis Lecture 18: Least Squares, State Estimation
ECEN 615 Methods of Electric Power Systems Analysis Lecture 18: Least Squares, State Estimation Prof. om Overbye Dept. of Electrical and Computer Engineering exas A&M University overbye@tamu.edu Announcements
More informationKalman Filter and Parameter Identification. Florian Herzog
Kalman Filter and Parameter Identification Florian Herzog 2013 Continuous-time Kalman Filter In this chapter, we shall use stochastic processes with independent increments w 1 (.) and w 2 (.) at the input
More informationSolving Linear Systems of Equations
Solving Linear Systems of Equations Gerald Recktenwald Portland State University Mechanical Engineering Department gerry@me.pdx.edu These slides are a supplement to the book Numerical Methods with Matlab:
More informationNumerical Methods. Elena loli Piccolomini. Civil Engeneering. piccolom. Metodi Numerici M p. 1/??
Metodi Numerici M p. 1/?? Numerical Methods Elena loli Piccolomini Civil Engeneering http://www.dm.unibo.it/ piccolom elena.loli@unibo.it Metodi Numerici M p. 2/?? Least Squares Data Fitting Measurement
More informationModelling Physical Phenomena
Modelling Physical Phenomena Limitations and Challenges of the Differential Algebraic Equations Approach Olaf Trygve Berglihn Department of Chemical Engineering 30. June 2010 2 Outline Background Classification
More informationThe integrating factor method (Sect. 1.1)
The integrating factor method (Sect. 1.1) Overview of differential equations. Linear Ordinary Differential Equations. The integrating factor method. Constant coefficients. The Initial Value Problem. Overview
More informationA THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS
A THEORETICAL INTRODUCTION TO NUMERICAL ANALYSIS Victor S. Ryaben'kii Semyon V. Tsynkov Chapman &. Hall/CRC Taylor & Francis Group Boca Raton London New York Chapman & Hall/CRC is an imprint of the Taylor
More informationTwo-Point Boundary Value Problem and Optimal Feedback Control based on Differential Algebra
Two-Point Boundary Value Problem and Optimal Feedback Control based on Differential Algebra Politecnico di Milano Department of Aerospace Engineering Milan, Italy Taylor Methods and Computer Assisted Proofs
More informationECE7850 Lecture 7. Discrete Time Optimal Control and Dynamic Programming
ECE7850 Lecture 7 Discrete Time Optimal Control and Dynamic Programming Discrete Time Optimal control Problems Short Introduction to Dynamic Programming Connection to Stabilization Problems 1 DT nonlinear
More informationMIT Final Exam Solutions, Spring 2017
MIT 8.6 Final Exam Solutions, Spring 7 Problem : For some real matrix A, the following vectors form a basis for its column space and null space: C(A) = span,, N(A) = span,,. (a) What is the size m n of
More informationCHAPTER 5: Linear Multistep Methods
CHAPTER 5: Linear Multistep Methods Multistep: use information from many steps Higher order possible with fewer function evaluations than with RK. Convenient error estimates. Changing stepsize or order
More informationMinimization with Equality Constraints
Path Constraints and Numerical Optimization Robert Stengel Optimal Control and Estimation, MAE 546, Princeton University, 2015 State and Control Equality Constraints Pseudoinverse State and Control Inequality
More informationLinear Algebra Exercises
9. 8.03 Linear Algebra Exercises 9A. Matrix Multiplication, Rank, Echelon Form 9A-. Which of the following matrices is in row-echelon form? 2 0 0 5 0 (i) (ii) (iii) (iv) 0 0 0 (v) [ 0 ] 0 0 0 0 0 0 0 9A-2.
More informationAn Implicit Runge Kutta Solver adapted to Flexible Multibody System Simulation
An Implicit Runge Kutta Solver adapted to Flexible Multibody System Simulation Johannes Gerstmayr 7. WORKSHOP ÜBER DESKRIPTORSYSTEME 15. - 18. March 2005, Liborianum, Paderborn, Germany Austrian Academy
More informationSuboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids
Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität
More informationDirect and indirect methods for optimal control problems and applications in engineering
Direct and indirect methods for optimal control problems and applications in engineering Matthias Gerdts Computational Optimisation Group School of Mathematics The University of Birmingham gerdtsm@maths.bham.ac.uk
More informationTrajectory Tracking Control of Bimodal Piecewise Affine Systems
25 American Control Conference June 8-1, 25. Portland, OR, USA ThB17.4 Trajectory Tracking Control of Bimodal Piecewise Affine Systems Kazunori Sakurama, Toshiharu Sugie and Kazushi Nakano Abstract This
More information