Sensitivity computation based on CasADi for dynamic optimization problems
|
|
- Claire Jennings
- 5 years ago
- Views:
Transcription
1 Sensitivity computation based on CasADi for dynamic optimization problems 16 th European Workshop on Automatic Differentiation 8 9 December 2014, Jena, Germany Evgeny Lazutkin, Abebe Geletu, Siegbert Hopfgarten, Pu Li Technische Universität Ilmenau Simulation and Optimal Processes Group evgeny.lazutkin@tu-ilmenau.de Page 1
2 Motivation Content Combined Multiple Shooting and Collocation (CMSC) method Application of the symbolic computations Sensitivity analysis First Order Sensitivity (FOS) Path constraints Summary and current work Page 2
3 1. Motivation Why? Safe operations Fast and efficient What? startup / shutdown High dimensional, dynamical, and nonlinear systems with constraints Where? Power plants Robotics Petrochemical industry Power plant¹ Petrochemical industry³ Robotics² ¹ ² ³ Page 3
4 1. Motivation Easy application Object-oriented modelling for optimization Solution methods Advanced numerical algorithms Nonlinear programming Combined multiple shooting and collocation Computational implementation Automatic and symbolic differentiation Complex control system FOS are required Fast analysis and solution Functions, sparsity, etc. Sensitivities Approximation Analytical Control system Offline solution Online solution (Nonlinear Model Predictive Control (NMPC)) Example Modelica models⁴ Example: Automatic differentiation⁵ ⁴ ⁵ Page 4
5 2. CMSC method [1,2,5] Original NOCP Multiple Shooting* min x(t),u(t) s.t. J(x t, u(t)) x t 0 = x 0, x t = f x t, u t, t, x mmm x t x mmm, u mmm u t u mmm. Collocation** Ipopt Sensitivity calculation Solution *Bock, H. G., Plitt, K. J. A Multiple Shooting Algorithm for Direct Solution of Optimal Control Problems, Prepr. 9th IFAC World Congress, Budapest, 1984, pp **Cuthrell, J. E., Biegler, L. T.: Simultaneous optimization and solution methods for batch reactor control profiles. Comput. Chem. Eng. 13(1989),pp Page 5
6 2. CMSC method [1,2,5] Original NOCP Multiple Shooting* min x(t),u(t) s.t. J(x t, u(t)) x t 0 = x 0, x t = f x t, u t, t, x mmm x t x mmm, u mmm u t u mmm. Collocation** Ipopt Sensitivity calculation Solution *Bock, H. G., Plitt, K. J. A Multiple Shooting Algorithm for Direct Solution of Optimal Control Problems, Prepr. 9th IFAC World Congress, Budapest, 1984, pp **Cuthrell, J. E., Biegler, L. T.: Simultaneous optimization and solution methods for batch reactor control profiles. Comput. Chem. Eng. 13(1989),pp Page 6
7 2. CMSC method [1,2,5] Original NOCP Multiple Shooting* Collocation** Ipopt Sensitivity calculation Solution *Bock, H. G., Plitt, K. J. A Multiple Shooting Algorithm for Direct Solution of Optimal Control Problems, Prepr. 9th IFAC World Congress, Budapest, 1984, pp **Cuthrell, J. E., Biegler, L. T.: Simultaneous optimization and solution methods for batch reactor control profiles. Comput. Chem. Eng. 13(1989),pp Page 7
8 2. CMSC method [1,2,5] Original NOCP Multiple Shooting* Collocation** Ipopt Sensitivity calculation Solution *Bock, H. G., Plitt, K. J. A Multiple Shooting Algorithm for Direct Solution of Optimal Control Problems, Prepr. 9th IFAC World Congress, Budapest, 1984, pp **Cuthrell, J. E., Biegler, L. T.: Simultaneous optimization and solution methods for batch reactor control profiles. Comput. Chem. Eng. 13(1989),pp Page 8
9 2. CMSC method [1,2,5] Original NOCP Multiple Shooting* Collocation** Ipopt Sensitivity calculation Solution *Bock, H. G., Plitt, K. J. A Multiple Shooting Algorithm for Direct Solution of Optimal Control Problems, Prepr. 9th IFAC World Congress, Budapest, 1984, pp **Cuthrell, J. E., Biegler, L. T.: Simultaneous optimization and solution methods for batch reactor control profiles. Comput. Chem. Eng. 13(1989),pp Page 9
10 2. CMSC method [1,2,5] Original NOCP Multiple Shooting* Collocation** Ipopt Sensitivity calculation Solution *Bock, H. G., Plitt, K. J. A Multiple Shooting Algorithm for Direct Solution of Optimal Control Problems, Prepr. 9th IFAC World Congress, Budapest, 1984, pp **Cuthrell, J. E., Biegler, L. T.: Simultaneous optimization and solution methods for batch reactor control profiles. Comput. Chem. Eng. 13(1989),pp Page 10
11 2. CMSC method ADVANTAGES Combines advantages of multiple shooting method and collocation method Reduction of the NLP variables in optimizer (Ipopt) Precision approximation of the state trajectories Can be applied for NMPC Allows parallel computation due to independence of each interval The symbolic and automatic pre-calculation of the derivatives Use Modelica models Implementation on JModelica.org and CasADi platform for wider public availability Page 11
12 3. Application of symbolic computations Original model formulation Modelica model formulation model batch Real x1(start=1, fixed=true, min=0,max=1); Real x2(start=0, fixed=true, min=0,max=1); input Real u1; equation der(x1) = -(u *u1*u1)*x1; der(x2) = u1*x1; end batch; optimization batch_ocp(objective= -x2(finaltime), starttime=0, finaltime=1) extends batch(u1(free=true,min=0,max=5)); end batch_ocp; Modelica model can be transferred into Python through JModelica.org software for further manipulations Page 12
13 3. Symbolic computations with CasADi Transferred ODE system in symbolic form SX_ODE = [(der(x1)+((u1+((0.5*u1)*u1))*x1)), (der(x2)-(u1*x1))] Inputs for ODE system INPUTS = [Matrix<SX>(der(x1)), Matrix<SX>(der(x2)), Matrix<SX>(x1), Matrix<SX>(x2), Matrix<SX>(u1)] SXFunction instance SX_ODE_FUNCTION = SXFunction([vertcat(INPUTS)], [SX_ODE]) SX_ODE_FUNCTION.init() Features Evaluation and AD Symbolically Numerically Required variables States at collocation points Parameterized states Parameterized controls state = ssym("s",n_int*(n_d N_CP)) => [s_0,s_1,s_2,s_3,s_4,s_5] p_state = ssym("ps",(n_int+1) N_D) p_ctrl = ssym("c",n_int N_C) Advantages Check results Problem independent Simple to implement Page 13
14 3. Application of symbolic computations Applied method for SXFunction: eval RES = SX_ODE_FUNCTION.eval([vertcat([der_expr, diff, ctrl])])[0] Time discretization by CMSC method: 1 intervals, 3 collocation points [((((( *s_0)+(1.0328*s_2))+( *s_4))+( *ps_0))+((c+((0.5*c)*c))*s_0)), ((((( *s_1)+(1.0328*s_3))+( *s_5))+( *ps_1))-(c*s_0))] [((((( *s_0)+(1.7746*s_2))+( *s_4))+(2.6619*ps_0))+((c+((0.5*c)*c))*s_2)), ((((( *s_1)+(1.7746*s_3))+( *s_5))+(2.6619*ps_1))-(c*s_2))] [((((( *s_0)+( *s_2))+( *s_4))+( *ps_0))+((c+((0.5*c)*c))*s_4)), ((((( *s_1)+( *s_3))+( *s_5))+( *ps_1))-(c*s_4))] Pack this system to new symbolic function G and create new SXFunction with only symbolical inputs The same system will be solved in a Newton solver to obtain state trajectories NewtonImplicitSolver KinsolSolver How to calculate required FOS? Page 14
15 4. Sensitivity analysis: FOS X ii X ii X ii = G X ii Sensitivity w.r.t. parameterized state Matrix equations X ii X ii U = G U Sensitivity w.r.t. parameterized control,, can be constructed using AD! X ii X ii U Function = SXFunction([vertcat([state, p_state, p_ctrl])],[vertcat([g])]) Function.init() full_jacobian=function.jac() dgdx_sym = SXMatrix(SYSTEM_INDEX,SYSTEM_INDEX) dgdx_sym = full_jacobian[range(system_index), range(system_index)] dgdx_sx = SXFunction([c_s,vertcat([p_s, p_c])],[dgdx_sym]) dgdx_sx.init() Step 0 Step 1.1 A = msym('a', dgdx_sym.sparsity()) Step 2.1 X.1 means matrix number SYSTEM_INDEX is a number of collocation variables per interval Page 15
16 4. Sensitivity analysis: FOS X ii X ii X ii = G X ii Sensitivity w.r.t. parameterized state Matrix equations,, can be constructed using AD! X ii X ii U X ii X ii U = G U Sensitivity w.r.t. parameterized control dgdxp_sym = SXMatrix(SYSTEM_INDEX,N_D) self.dgdxp_sym = full_jacobian[system_index:system_index+n_d]] dgdxp_sx = SXFunction([d_s,vertcat([p_s, p_c])],[dgdxp_sym]) dgdxp_sx.init() Step 1.2 B = msym( B', dgdxp_sym.sparsity()) SOLVER_X = CSparse(A.sparsity(), shape(b)[1]) SOLVER_X.init() SOLVER_X_SOL = SOLVER_X.solve(A,B) Step 2.2 Step 3 Create symbolical MXFunction to obtain sensitivities Page 16
17 4. Sensitivity analysis: Path constraints X ii = G X ii X ii X ii Sensitivity w.r.t. parameterized state Matrix equations x 1 2 x 0,1, x 0,2, u + x 2 2 x 0,1, x 0,2, u 1 0 X ii X ii U = G U Sensitivity w.r.t. parameterized control How to calculate sensitivity for optimizer automatically? pc_function = SXFunction([vertcat([ORIGINAL_STATES])],[vertcat([PATH])]) pc_function.init() GRADINET=pc_function.grad() x 1 [2 x 1 x 0,1, x 0,2, u, 2 x 2 (x 0,1, x 0,2, u)] x 0,1 x 2 AND x 1 u Automatic approach via CasADi (even for complicated path constraints) x 0,1 Page 17
18 5. Summary and current work Summary Pre-calculated derivatives and their symbolic representation Sensitivities are automatically available Using Modelica model formulation (limited to CasADi 1.8.2) CasADi package (not only for AD) Can handle complicated expressions Symbolical and numerical evaluation / Solvers Documentation via API, examples, Google group, manual Limited to explicit form of sensitivities Implementation time Current work Investigation of the CPU/GPU parallelization Analytic Hessian matrix Page 18
19 References [1] Tamimi, J.: Development of the Efficient Algorithms for Model Predictive Control of Fast Systems. PhD Thesis, Technische Universität Ilmenau, VDI Verlag, [2] Tamimi, J., Li, P.: A Combined approach to nonlinear model predictive control of fast systems. J. Process Control, 20(2010)9, pp [3] Andersson, J., Åkesson, J., Diehl, M.: CasADi: A Symbolic Package for Automatic Differentiation and Optimal Control. Lecture Notes in Computational Science and Engineering, Vol. 87, Springer, 2012, pp [4] Wächter, A., Biegler, L.T.: On the Implementation of a Primal-Dual Interior Point Filter Line Search Algorithm for Large-Scale Nonlinear Programming, Mathematical Programming, Ser. A, 106(2006)1, pp [5] E. Lazutkin, A. Geletu, S. Hopfgarten, P. Li: Modified Multiple Shooting Combined with Collocation Method in JModelica.org with Symbolic Calculations, Proceedings of the 10th International Modelica Conference, DOI /ECP , March, 2014, Lund, Sweden, pp Page 19
20 Acknowledgments This work has been supported by Model Driven Physical Systems Operation project (MODRIO) by ITEA2, No , and by the German BMBF (BMBF Förderkennzeichen: 01IS12022H). Page 20
21 DISCUSSION Thank you very much for your attention! Page 21
Model-Based Dynamic Optimization with OpenModelica and CasADi
Model-Based Dynamic Optimization with OpenModelica and CasADi Alachew Shitahun PELAB Programming Environment Lab, Dept. Computer Science Linköping University, SE-581 83 Linköping, Sweden Vitalij Ruge Mathematics
More informationNonlinear Model Predictive Control Using Multiple Shooting Combined with Collocation on Finite Elements
Nonlinear Model Predictive Control Using Multiple Shooting Combined with Collocation on Finite Elements Jasem Tamimi and Pu Li Simulation and Optimal Processes Group, Institute of Automation and Systems
More informationCasADi tutorial Introduction
Lund, 6 December 2011 CasADi tutorial Introduction Joel Andersson Department of Electrical Engineering (ESAT-SCD) & Optimization in Engineering Center (OPTEC) Katholieke Universiteit Leuven OPTEC (ESAT
More informationR. Franke, PSP-A24, OpenModelica Annual Workshop 2015 Mathematical Optimization of Dynamic Systems with OpenModelica. ABB February 2, 2015
R. Frane, PSP-A24, 02.02.2015 OpenModelica Annual Worshop 2015 Mathematical Optimization of Dynamic Systems with OpenModelica Slide 1 Overview Introduction to ABB OPTIMAX DrumBoiler example 1. Setpoint
More informationNMPC Application using JModelica.org: Features and Performance
NMPC Application using JModelica.org: Features and Performance Christian Hartlep 1 Toivo Henningsson 2 1 Siemens AG, Germany,christian.hartlep.ext@siemens.com 2 Modelon AB, Sweden,toivo.henningsson@modelon.com
More informationOptimal Control Techniques for Dynamic Walking
Optimal Control Techniques for Dynamic Walking Optimization in Robotics & Biomechanics IWR, University of Heidelberg Presentation partly based on slides by Sebastian Sager, Moritz Diehl and Peter Riede
More informationMultiple-Shooting Optimization using the JModelica.org Platform
Multiple-Shooting Optimization using the JModelica.org Platform Jens Rantil* Johan Åkesson Claus Führer* Magnus Gäfvert ** Department of Numerical Analysis, Lund University Department of Automatic Control,
More informationDynamic Optimization in JModelica.org
Dynamic Optimization in JModelica.org Magnusson, Fredrik; Åkesson, Johan Published in: Processes DOI: 10.3390/pr3020471 2015 Document Version: Publisher's PDF, also known as Version of record Link to publication
More informationMagnusson, Fredrik; Berntorp, Karl; Olofsson, Björn; Åkesson, Johan
Symbolic Transformations of Dynamic Optimization Problems Magnusson, Fredrik; Berntorp, Karl; Olofsson, Björn; Åkesson, Johan Published in: Proceedings of the 1th International Modelica Conference 214
More informationACADO TOOLKIT - AUTOMATIC CONTROL AND DYNAMIC OPTIMIZATION
ACADO Toolkit - Automatic Control and Dynamic Optimization p. 1/24 ACADO TOOLKIT - AUTOMATIC CONTROL AND DYNAMIC OPTIMIZATION Moritz Diehl, Boris Houska, Hans Joachim Ferreau Optimization in Engineering
More informationOutline 19/3/ Introduction. 2. Current Status of qpoases. 3. Future Developments
Hans Joachim Ferreau, ABB Corporate Research, Switzerland Past, Present, Future Outline 1. Introduction 2. Current Status of 3. Future Developments Slide 2 1 Outline 1. Introduction 2. Current Status of
More informationReal-Time Optimization for Fast Nonlinear MPC: Algorithms, Theory, and Applications
Real-Time Optimization for Fast Nonlinear MPC: Algorithms, Theory, and Applications Moritz Diehl Optimization in Engineering Center OPTEC & ESAT, K.U. Leuven Joint work with H. J. Ferreau*, B. Houska*,
More informationModeling and Optimization with Modelica and Optimica Using the JModelica.org Open Source Platform
Modeling and Optimization with Modelica and Optimica Using the JModelica.org Open Source Platform Johan Åkesson a,b ToveBergdahl a MagnusGäfvert a HubertusTummescheit a a ModelonAB, Sweden b DepartmentofAutomaticControl,
More informationNumerical and Symbolic Methods for Dynamic Optimization
Numerical and Symbolic Methods for Dynamic Optimization Magnusson, Fredrik Published: 2016-11-18 Document Version Publisher's PDF, also known as Version of record Link to publication Citation for published
More informationModeling, Calibration and Control of a Paper Machine Dryer Section
Modeling, Calibration and Control of a Paper Machine Dryer Section Modeling, Calibration and Control of a Paper Machine Dryer Section JohanÅkesson 1 OlaSlätteke 2 1 DepartmentofAutomaticControl,FacultyofEngineering,LundUniversity,Sweden
More informationSolving Optimal Control Problems with ACADO Toolkit
Solving Optimal Control Problems with ACADO Toolkit Boris Houska, Hans Joachim Ferreau, Moritz Diehl Electrical Engineering Department K.U. Leuven OPTEC Seminar, 2/9/2009 Electrical Engineering Department
More informationWORHP Lab The Graphical User Interface for Optimisation and Optimal Control
WORHP Lab The Graphical User Interface for Optimisation and Optimal Control M. Knauer, C. Büskens Zentrum für Universität Bremen 3rd European Optimisation in Space Engineering 17 th - 18 th September 2015
More informationIntegration of CasADi and JModelica.org
Integration of CasADi and JModelica.org Joel Andersson c Johan Åkesson a,b Francesco Casella d Moritz Diehl c a Department of Automatic Control, Lund University, Sweden b Modelon AB, Sweden c Department
More informationINTERIOR POINT METHOD BASED CONTACT ALGORITHM FOR STRUCTURAL ANALYSIS OF ELECTRONIC DEVICE MODELS
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationUnderstanding Concepts of Optimization and Optimal Control with WORHP Lab
Understanding Concepts of Optimization and Optimal Control with WORHP Lab M. Knauer, C. Büskens Zentrum für Universität Bremen 6th International Conference on Astrodynamics Tools and Techniques 14 th -
More informationUser Documentation for CasADi v Joel Andersson Joris Gillis Moritz Diehl
User Documentation for CasADi v2.0.0 Joel Andersson Joris Gillis Moritz Diehl July 25, 2014 2 Contents 1 Introduction 5 1.1 What CasADi is and what it is not...................... 5 1.2 Help and support................................
More informationpyomo.dae: A Modeling and Automatic Discretization Framework for Optimization with Differential and Algebraic Equations
Noname manuscript No. (will be inserted by the editor) pyomo.dae: A Modeling and Automatic Discretization Framework for Optimization with Differential and Algebraic Equations Bethany Nicholson John D.
More informationGrey-Box Building Models for Model Order Reduction and Control
Grey-Box Building Models for Model Order Reduction and Control Roel De Coninck a,b, Fredrik Magnusson c, Johan Åkesson c,d, Lieve Helsen b a 3E nv., 1000 Brussels, Belgium, b KU Leuven, Department of Mechanical
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 1 3.1 Linearization and Optimization of Functions of Vectors 1 Problem Notation 2 Outline 3.1.1 Linearization 3.1.2 Optimization of Objective Functions 3.1.3 Constrained
More informationPerformance Evaluation of an Interior Point Filter Line Search Method for Constrained Optimization
6th WSEAS International Conference on SYSTEM SCIENCE and SIMULATION in ENGINEERING, Venice, Italy, November 21-23, 2007 18 Performance Evaluation of an Interior Point Filter Line Search Method for Constrained
More informationA New Tool Providing an Integrated Framework for Process Optimization
EngOpt 2008 - International Conference on Engineering Optimization Rio de Janeiro, Brazil, 01-05 June 2008. A New Tool Providing an Integrated Framework for Process Optimization E.C. do Valle; R.P. Soares;
More informationNew features in dynamic optimization with EcosimPro
Journal of Engineering and Technology for Industrial Applications, 2018. Edition. 13.Vol: 04 https://www.itegam-jetia.org ISSN ONLINE: 2447-0228 DOI: https://dx.doi.org/10.5935/2447-0228.20180001 New features
More informationOlofsson, Björn; Nilsson, Henrik; Robertsson, Anders; Åkesson, Johan
Optimal Tracking and Identification of Paths for Industrial Robots Olofsson, Björn; Nilsson, Henrik; Robertsson, Anders; Åkesson, Johan Published in: Proceedings of the 8th IFAC World Congress, 2 2 Document
More informationFast Direct Multiple Shooting Algorithms for Optimal Robot Control
Fast Direct Multiple Shooting Algorithms for Optimal Robot Control Moritz Diehl, Hans Georg Bock, Holger Diedam, Pierre-Brice Wieber To cite this version: Moritz Diehl, Hans Georg Bock, Holger Diedam,
More informationParallel Multiple-Shooting and Collocation Optimization with OpenModelica
arallel Multiple-Shooting and Collocation Optimization with OpenModelica Bernhard Bachmann 3, Lennart Ochel 3, Vitalij Ruge 3, Mahder Gebremedhin 1, eter Fritzson 1, Vaheed Nezhadali, Lars Eriksson, Martin
More informationThis is the published version of a paper presented at IEEE PES General Meeting 2013.
http://www.diva-portal.org This is the published version of a paper presented at IEEE PES General Meeting 2013. Citation for the original published paper: Vanfretti, L., Li, W., Bogodorova, T., Panciatici,
More informationEfficient implementation of Constrained Min-Max Model Predictive Control with Bounded Uncertainties
Efficient implementation of Constrained Min-Max Model Predictive Control with Bounded Uncertainties D.R. Ramírez 1, T. Álamo and E.F. Camacho2 Departamento de Ingeniería de Sistemas y Automática, Universidad
More informationComparison of Interior Point Filter Line Search Strategies for Constrained Optimization by Performance Profiles
INTERNATIONAL JOURNAL OF MATHEMATICS MODELS AND METHODS IN APPLIED SCIENCES Comparison of Interior Point Filter Line Search Strategies for Constrained Optimization by Performance Profiles M. Fernanda P.
More informationGPC AND NEURAL GENERALIZED PREDICTIVE CONTROL
S. Chidrawar Nikhil Bidwai, L. Waghmare and B. M. Patre MGM s College of Engineering, Nanded (MS) 43 60, India SGGS Institute of Engineering and Technology, Nanded (MS) 43 606, India sadhana_kc@rediffmail.com
More informationCURRICULUM VITAE: DI MARTIN BERNAUER. DI Martin Bernauer. Chemnitz University of Technology Faculty of Mathematics D Chemnitz Germany
CURRICULUM CURRICULUM VITAE: DI MARTIN BERNAUER Chemnitz, March 15, 2010 DI Martin Bernauer Chemnitz University of Technology Faculty of Mathematics D-09107 Chemnitz Germany +49-371-531-37497 martin.bernauer@mathematik.tu-chemnitz.de
More informationIncremental Single Shooting - A Robust Method for the Estimation of Parameters in Dynamical Systems
Aachen Institute for Advanced Study in Computational Engineering Science Preprint: AICES-2009-8 28/January/2009 Incremental Single Shooting - A Robust Method for the Estimation of Parameters in Dynamical
More informationParametric Sensitivity Analysis of NLP Problems and its Applications to Real-Time Controller Design
SADCO - Sensitivity Analysis for Deterministic Controller Design Parametric Sensitivity Analysis of NLP Problems and its Applications to Real-Time Controller Design Part 1: Nonlinear Real-Time Optimization
More informationOnline algorithms for clustering problems
University of Szeged Department of Computer Algorithms and Artificial Intelligence Online algorithms for clustering problems Summary of the Ph.D. thesis by Gabriella Divéki Supervisor Dr. Csanád Imreh
More informationDeveloping Optimization Algorithms for Real-World Applications
Developing Optimization Algorithms for Real-World Applications Gautam Ponnappa PC Training Engineer Viju Ravichandran, PhD Education Technical Evangelist 2015 The MathWorks, Inc. 1 2 For a given system,
More informationFlight Systems are Cyber-Physical Systems
Flight Systems are Cyber-Physical Systems Dr. Christopher Landauer Software Systems Analysis Department The Aerospace Corporation Computer Science Division / Software Engineering Subdivision 08 November
More informationACADO Code Generation tool
ACADO Code Generation tool Rien Quirynen July 31, 2015 Introduction Automatic Code Generation Real-Time Iterations Application examples ACADO demo Outline Introduction Automatic Code Generation Real-Time
More informationReal-Time Simulation of Modelica-based Models
Real-Time Simulation of Modelica-based Models Torsten Blochwitz Thomas Beutlich ITI GmbH Webergasse 1, 01067 Dresden, Germany {blochwitz,beutlich}@iti.de Abstract This paper shows the various steps a simulation
More informationModern Methods of Data Analysis - WS 07/08
Modern Methods of Data Analysis Lecture XV (04.02.08) Contents: Function Minimization (see E. Lohrmann & V. Blobel) Optimization Problem Set of n independent variables Sometimes in addition some constraints
More informationCGT 581 G Geometric Modeling Curves
CGT 581 G Geometric Modeling Curves Bedrich Benes, Ph.D. Purdue University Department of Computer Graphics Technology Curves What is a curve? Mathematical definition 1) The continuous image of an interval
More informationExperience with Industrial In-House Application of FMI
Kilian Link 1 Leo Gall 2 Monika Mühlbauer 3 Stephanie Gallardo-Yances 4 1,3,4 Siemens AG, Germany,{kilian.link, monika.muehlbauer, stephanie.gallardo}@siemens.com 2 LTX Simulation GmbH, Germany,leo.gall@ltx.de
More informationLECTURE NOTES Non-Linear Programming
CEE 6110 David Rosenberg p. 1 Learning Objectives LECTURE NOTES Non-Linear Programming 1. Write out the non-linear model formulation 2. Describe the difficulties of solving a non-linear programming model
More informationInverse KKT Motion Optimization: A Newton Method to Efficiently Extract Task Spaces and Cost Parameters from Demonstrations
Inverse KKT Motion Optimization: A Newton Method to Efficiently Extract Task Spaces and Cost Parameters from Demonstrations Peter Englert Machine Learning and Robotics Lab Universität Stuttgart Germany
More informationOPTIMIZATION METHODS
D. Nagesh Kumar Associate Professor Department of Civil Engineering, Indian Institute of Science, Bangalore - 50 0 Email : nagesh@civil.iisc.ernet.in URL: http://www.civil.iisc.ernet.in/~nagesh Brief Contents
More informationNonlinear State Estimation for Robotics and Computer Vision Applications: An Overview
Nonlinear State Estimation for Robotics and Computer Vision Applications: An Overview Arun Das 05/09/2017 Arun Das Waterloo Autonomous Vehicles Lab Introduction What s in a name? Arun Das Waterloo Autonomous
More informationDynamic Analysis of Manipulator Arm for 6-legged Robot
American Journal of Mechanical Engineering, 2013, Vol. 1, No. 7, 365-369 Available online at http://pubs.sciepub.com/ajme/1/7/42 Science and Education Publishing DOI:10.12691/ajme-1-7-42 Dynamic Analysis
More informationIlmenau University of Technology
Ilmenau University of Technology Computer Science and Automation Faculty Institute for Automation and Systems Engineering Simulation and Optimal Processes MASTER THESIS Comparison of hierarchical and distributed
More informationEllipsoid Algorithm :Algorithms in the Real World. Ellipsoid Algorithm. Reduction from general case
Ellipsoid Algorithm 15-853:Algorithms in the Real World Linear and Integer Programming II Ellipsoid algorithm Interior point methods First polynomial-time algorithm for linear programming (Khachian 79)
More informationOptimized Control of Hot-Gas Cycle for Solar Thermal Power Plants
Optimized Control of Hot-Gas Cycle for Solar Thermal Power Plants Jan Gall Dirk Abel IRT Institut für Regelungstechnik, RWTH Aachen, Steinbachstr. 54, D-52074 Aachen, Germany j.gall@irt.rwth-aachen.de
More informationSymbolic Evaluation of Sums for Parallelising Compilers
Symbolic Evaluation of Sums for Parallelising Compilers Rizos Sakellariou Department of Computer Science University of Manchester Oxford Road Manchester M13 9PL United Kingdom e-mail: rizos@csmanacuk Keywords:
More informationOutline. Robust MPC and multiparametric convex programming. A. Bemporad C. Filippi. Motivation: Robust MPC. Multiparametric convex programming
Robust MPC and multiparametric convex programming A. Bemporad C. Filippi D. Muñoz de la Peña CC Meeting Siena /4 September 003 Outline Motivation: Robust MPC Multiparametric convex programming Kothares
More informationFramework for dynamic optimization of district heating systems using Optimica Compiler Toolkit
Framework for dynamic optimization of district heating systems using Optimica Compiler Toolkit Gerald Schweiger 1 Håkan Runvik 2 Fredrik Magnusson 3,2 Per-Ola Larsson 2 Stéphane Velut 2 1 AEE INTEC, 8200
More informationLecture 23: Starting to put it all together #2... More 2-Point Boundary value problems
Lecture 23: Starting to put it all together #2... More 2-Point Boundary value problems Outline 1) Our basic example again: -u'' + u = f(x); u(0)=α, u(l)=β 2) Solution of 2-point Boundary value problems
More informationA PACKAGE FOR DEVELOPMENT OF ALGORITHMS FOR GLOBAL OPTIMIZATION 1
Mathematical Modelling and Analysis 2005. Pages 185 190 Proceedings of the 10 th International Conference MMA2005&CMAM2, Trakai c 2005 Technika ISBN 9986-05-924-0 A PACKAGE FOR DEVELOPMENT OF ALGORITHMS
More informationOlivier Brüls. Department of Aerospace and Mechanical Engineering University of Liège
Fully coupled simulation of mechatronic and flexible multibody systems: An extended finite element approach Olivier Brüls Department of Aerospace and Mechanical Engineering University of Liège o.bruls@ulg.ac.be
More informationPETRI NET ANALYSIS OF BATCH RECIPES
Presented at FOCAPO 98, Snowbird, USA. PETRI NET ANALYSIS OF BATCH RECIPES STRUCTURED WITH GRAFCHART Charlotta Johnsson and Karl-Erik Årzén Department of Automatic Control, Lund Institute of Technology,
More informationVirtual Manikin Controller Calculating the movement of a human model Master s thesis in Complex Adaptive Systems DANIEL GLEESON
Virtual Manikin Controller Calculating the movement of a human model Master s thesis in Complex Adaptive Systems DANIEL GLEESON Department of Applied Mechanics Division of Vehicle Engineering and Autonomous
More informationPreface. and Its Applications 81, ISBN , doi: / , Springer Science+Business Media New York, 2013.
Preface This book is for all those interested in using the GAMS technology for modeling and solving complex, large-scale, continuous nonlinear optimization problems or applications. Mainly, it is a continuation
More informationHardware Tailored Linear Algebra for Implicit Integrators in Embedded NMPC.
Downloaded from orbit.dtu.dk on: Jan 12, 2019 Hardware Tailored Linear Algebra for Implicit Integrators in Embedded NMPC. Frison, Gianluca; Quirynen, Rien; Zanelli, Andrea; Diehl, Moritz; Jørgensen, John
More informationIntroduction to Optimization Problems and Methods
Introduction to Optimization Problems and Methods wjch@umich.edu December 10, 2009 Outline 1 Linear Optimization Problem Simplex Method 2 3 Cutting Plane Method 4 Discrete Dynamic Programming Problem Simplex
More informationInternal-growth Strategies for Dynamic Process Optimization with Differential-algebraic Equations
Preprints of the 19th World Congress The International Federation of Automatic Control Internal-growth Strategies for Dynamic Process Optimization with Differential-algebraic Equations Zhiqiang Wang, Zhijiang
More informationTrajectory Optimization
C H A P T E R 12 Trajectory Optimization So far, we have discussed a number of ways to solve optimal control problems via state space search (e.g., Dijkstra s and Dynamic Programming/Value Iteration).
More informationCompressive Sensing Based Image Reconstruction using Wavelet Transform
Compressive Sensing Based Image Reconstruction using Wavelet Transform Sherin C Abraham #1, Ketki Pathak *2, Jigna J Patel #3 # Electronics & Communication department, Gujarat Technological University
More informationDLR FF-DR-ER, Technical Report TR R101-93, March
DLR FF-DR-ER, Technical Report TR R101-93, March 1993. 1 The ANDECS Simulation Environment DSSIM Martin Otter Institute for Robotics and System Dynamics German Aerospace Research Establishment (DLR), Oberpfaenhofen
More informationParameter Estimation in Differential Equations: A Numerical Study of Shooting Methods
Parameter Estimation in Differential Equations: A Numerical Study of Shooting Methods Franz Hamilton Faculty Advisor: Dr Timothy Sauer January 5, 2011 Abstract Differential equation modeling is central
More informationUsing Local Trajectory Optimizers To Speed Up Global. Christopher G. Atkeson. Department of Brain and Cognitive Sciences and
Using Local Trajectory Optimizers To Speed Up Global Optimization In Dynamic Programming Christopher G. Atkeson Department of Brain and Cognitive Sciences and the Articial Intelligence Laboratory Massachusetts
More informationautomatic differentiation techniques used in jump
automatic differentiation techniques used in jump Miles Lubin, Iain Dunning, and Joey Huchette June 22, 2016 MIT Operations Research Center 1 / 36 - Solver-independent, fast, extensible, open-source algebraic
More informationFrom Theory to Application (Optimization and Optimal Control in Space Applications)
From Theory to Application (Optimization and Optimal Control in Space Applications) Christof Büskens Optimierung & Optimale Steuerung 02.02.2012 The paradox of mathematics If mathematics refer to reality,
More informationParallel Architectures and Algorithms for Large-Scale Nonlinear Programming
Parallel Architectures and Algorithms for Large-Scale Nonlinear Programming Carl D. Laird Associate Professor, School of Chemical Engineering, Purdue University Faculty Fellow, Mary Kay O Connor Process
More informationEfficient Synthesis of Production Schedules by Optimization of Timed Automata
Efficient Synthesis of Production Schedules by Optimization of Timed Automata Inga Krause Institute of Automatic Control Engineering Technische Universität München inga.krause@mytum.de Joint Advanced Student
More informationUsing SPARK as a Solver for Modelica. Michael Wetter Philip Haves Michael A. Moshier Edward F. Sowell. July 30, 2008
Using SPARK as a Solver for Modelica Michael Wetter Philip Haves Michael A. Moshier Edward F. Sowell July 30, 2008 1 Overview Overview of SPARK, Modelica, OpenModelica, Dymola Problem reduction SPARK integration
More informationA software environment for simultaneous dynamic optimization
Computers and Chemical Engineering 31 (2007) 931 942 A software environment for simultaneous dynamic optimization Y.-D. Lang, L.T. Biegler Department of Chemical Engineering, Carnegie Mellon University,
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. XVI - Control Reconfiguration - Jan Lunze
CONTROL RECONFIGURATION Jan Lunze The Institute of Automation and Computer Control, Ruhr University Bochum, Germany Keywords: Model-matching, Reconfigurability, Reconfiguration, Recoverability. Contents
More informationRecent developments in simulation, optimization and control of flexible multibody systems
Recent developments in simulation, optimization and control of flexible multibody systems Olivier Brüls Department of Aerospace and Mechanical Engineering University of Liège o.bruls@ulg.ac.be Katholieke
More informationProblem characteristics. Dynamic Optimization. Examples. Policies vs. Trajectories. Planning using dynamic optimization. Dynamic Optimization Issues
Problem characteristics Planning using dynamic optimization Chris Atkeson 2004 Want optimal plan, not just feasible plan We will minimize a cost function C(execution). Some examples: C() = c T (x T ) +
More informationAn Efficient Method for Solving the Direct Kinematics of Parallel Manipulators Following a Trajectory
An Efficient Method for Solving the Direct Kinematics of Parallel Manipulators Following a Trajectory Roshdy Foaad Abo-Shanab Kafr Elsheikh University/Department of Mechanical Engineering, Kafr Elsheikh,
More informationMOSAIC: an Online Platform for Combined Process Model and Measurement Data Management
MOSAIC: an Online Platform for Combined Process Model and Measurement Data Management Erik Esche*, David Müller**, Robert Kraus***, Sandra Fillinger ****, Victor Alejandro Merchan Restrepo*****, Günter
More informationEXPRESSING REQUIREMENTS IN MODELICA
EXPRESSING REQUIREMENTS IN MODELICA Lena Buffoni and Peter Fritzson Linköping University SE-581 83 Linköping Sweden ABSTRACT As cyber-physical systems grow increasingly complex, the need for methodologies
More informationSystem Identification Algorithms and Techniques for Systems Biology
System Identification Algorithms and Techniques for Systems Biology by c Choujun Zhan A Thesis submitted to the School of Graduate Studies in partial fulfillment of the requirements for the degree of Doctor
More informationMBS MODELLING WITH SIMMECHANICS: CASE STUDIES IN RESEARCH AND EDUCATION
MBS MODELLING WITH SIMMECHANICS: CASE STUDIES IN RESEARCH AND EDUCATION Grepl, R., Lee, B., Singule, V., Švejda, P., Vlachý, D., Zezula, P. Laboratory of mechatronics, FME, Brno University of Technology
More informationSimple shooting-projection method for numerical solution of two-point Boundary Value Problems
Simple shooting-projection method for numerical solution of two-point Boundary Value Problems Stefan M. Filipo Ivan D. Gospodino a Department of Programming and Computer System Application, University
More informationIntroduction to Optimization
Introduction to Optimization Second Order Optimization Methods Marc Toussaint U Stuttgart Planned Outline Gradient-based optimization (1st order methods) plain grad., steepest descent, conjugate grad.,
More informationFast-Lipschitz Optimization
Fast-Lipschitz Optimization DREAM Seminar Series University of California at Berkeley September 11, 2012 Carlo Fischione ACCESS Linnaeus Center, Electrical Engineering KTH Royal Institute of Technology
More informationqpoases - Online Active Set Strategy for Fast Linear MPC
qpoases - Online Active Set Strategy for Fast Linear MPC Moritz Diehl, Hans Joachim Ferreau, Lieboud Vanden Broeck, Jan Swevers Dept. ESAT and Center of Excellence for Optimization in Engineering OPTEC
More informationWavelet-Galerkin Solutions of One and Two Dimensional Partial Differential Equations
VOL 3, NO0 Oct, 202 ISSN 2079-8407 2009-202 CIS Journal All rights reserved http://wwwcisjournalorg Wavelet-Galerkin Solutions of One and Two Dimensional Partial Differential Equations Sabina, 2 Vinod
More informationLOAD SHEDDING AN EFFICIENT USE OF LTC TRANSFORMERS
LOAD SHEDDING AN EFFICIENT USE OF LTC TRANSFORMERS Luciano V. Barboza André A. P. Lerm Roberto S. Salgado Catholic University of Pelotas Federal Center for Technological Education of Pelotas Federal University
More informationTechnische Universität München Zentrum Mathematik
Technische Universität München Zentrum Mathematik Prof. Dr. Dr. Jürgen Richter-Gebert, Bernhard Werner Projective Geometry SS 28 https://www-m.ma.tum.de/bin/view/lehre/ss8/pgss8/webhome Solutions for Worksheet
More informationCan we quantify the hardness of learning manipulation? Kris Hauser Department of Electrical and Computer Engineering Duke University
Can we quantify the hardness of learning manipulation? Kris Hauser Department of Electrical and Computer Engineering Duke University Robot Learning! Robot Learning! Google used 14 identical robots 800,000
More informationLaGO. Ivo Nowak and Stefan Vigerske. Humboldt-University Berlin, Department of Mathematics
LaGO a Branch and Cut framework for nonconvex MINLPs Ivo Nowak and Humboldt-University Berlin, Department of Mathematics EURO XXI, July 5, 2006 21st European Conference on Operational Research, Reykjavik
More informationResearch on time optimal trajectory planning of 7-DOF manipulator based on genetic algorithm
Acta Technica 61, No. 4A/2016, 189 200 c 2017 Institute of Thermomechanics CAS, v.v.i. Research on time optimal trajectory planning of 7-DOF manipulator based on genetic algorithm Jianrong Bu 1, Junyan
More informationDIRECT SEQUENTIAL DYNAMIC OPTIMIZATION WITH AUTOMATIC SWITCHING STRUCTURE DETECTION. Martin Schlegel, Wolfgang Marquardt 1
DIRECT SEQUENTIAL DYNAMIC OPTIMIZATION WITH AUTOMATIC SWITCHING STRUCTURE DETECTION Martin Schlegel, Wolfgang Marquardt 1 Lehrstuhl für Prozesstechnik, RWTH Aachen University D 52056 Aachen, Germany Abstract:
More informationLECTURE 6: INTERIOR POINT METHOD. 1. Motivation 2. Basic concepts 3. Primal affine scaling algorithm 4. Dual affine scaling algorithm
LECTURE 6: INTERIOR POINT METHOD 1. Motivation 2. Basic concepts 3. Primal affine scaling algorithm 4. Dual affine scaling algorithm Motivation Simplex method works well in general, but suffers from exponential-time
More informationA Framework for the Simulation and Validation of Distributed Control Architectures for Technical Systems of Systems
epartment of Biochemical and Chemical Engineering Process ynamics Group () A Framework for the Simulation and Validation of istributed Control Architectures for Technical Systems of Systems S. azari *,
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationAone-phaseinteriorpointmethodfornonconvexoptimization
Aone-phaseinteriorpointmethodfornonconvexoptimization Oliver Hinder, Yinyu Ye Department of Management Science and Engineering Stanford University April 28, 2018 Sections 1 Motivation 2 Why a one-phase
More informationPrediction-based diagnosis and loss prevention using qualitative multi-scale models
European Symposium on Computer Arded Aided Process Engineering 15 L. Puigjaner and A. Espuña (Editors) 2005 Elsevier Science B.V. All rights reserved. Prediction-based diagnosis and loss prevention using
More information