CDS 110b: Lecture 1-2 Introduction to Optimal Control

Similar documents
The Linear Quadratic Regulator

CDS 110b: Lecture 2-1 Linear Quadratic Regulators

Formal Methods for Deriving Element Equations

Partial Differential Equations with Applications

Optimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications

Robust Tracking and Regulation Control of Uncertain Piecewise Linear Hybrid Systems

Stability of Model Predictive Control using Markov Chain Monte Carlo Optimisation

STABILIZATIO ON OF LONGITUDINAL AIRCRAFT MOTION USING MODEL PREDICTIVE CONTROL AND EXACT LINEARIZATION

Linear and Nonlinear Model Predictive Control of Quadruple Tank Process

Move Blocking Strategies in Receding Horizon Control

Adaptive Fault-tolerant Control with Control Allocation for Flight Systems with Severe Actuator Failures and Input Saturation

Symbol R R + C Table : Notation Meaning Set of all real nmbers Set of positive real nmbers Set of all complex nmbers A(s) T, conjgate transpose A(s) Λ

Dynamic Optimization of First-Order Systems via Static Parametric Programming: Application to Electrical Discharge Machining

M 21 M 22 1 S. . x. x u B D A C

Lecture 13: Duality Uses and Correspondences

Robot Dynamics Fixed Wing UAS: Control and Solar UAS

Mixed Type Second-order Duality for a Nondifferentiable Continuous Programming Problem

Designing MIPS Processor

Nonlinear parametric optimization using cylindrical algebraic decomposition

AN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS

AN ALTERNATIVE DECOUPLED SINGLE-INPUT FUZZY SLIDING MODE CONTROL WITH APPLICATIONS

A Model-Free Adaptive Control of Pulsed GTAW

Fast Algorithms for Restoration of Color Wireless Capsule Endoscopy Images

Discontinuous Fluctuation Distribution for Time-Dependent Problems

FRTN10 Exercise 12. Synthesis by Convex Optimization

Online Solution of State Dependent Riccati Equation for Nonlinear System Stabilization

Safe Manual Control of the Furuta Pendulum

Concepts Introduced. Digital Electronics. Logic Blocks. Truth Tables

1. LQR formulation 2. Selection of weighting matrices 3. Matlab implementation. Regulator Problem mm3,4. u=-kx

Chapter 4 Supervised learning:

Approximate Solution of Convection- Diffusion Equation by the Homotopy Perturbation Method

Linear System Theory (Fall 2011): Homework 1. Solutions

Review of course Nonlinear control Lecture 5. Lyapunov based design Torkel Glad Lyapunov design Control Lyapunov functions Control Lyapunov function

Hybrid modelling and model reduction for control & optimisation

EE2 Mathematics : Functions of Multiple Variables

Lecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2

Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers

This is a repository copy of Enhanced Bounded Integral Control of Input-to-State Stable Nonlinear Systems.

Gradient Projection Anti-windup Scheme on Constrained Planar LTI Systems. Justin Teo and Jonathan P. How

Pulses on a Struck String

Restricted Three-Body Problem in Different Coordinate Systems

Affine Invariant Total Variation Models

Nonparametric Identification and Robust H Controller Synthesis for a Rotational/Translational Actuator

Model Predictive Control Based Energy Management Algorithm for a Hybrid Excavator

Chapter 2 Difficulties associated with corners

A Computational Study with Finite Element Method and Finite Difference Method for 2D Elliptic Partial Differential Equations

OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1

Outline. Model Predictive Control: Current Status and Future Challenges. Separation of the control problem. Separation of the control problem

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.

Chapter 2 Introduction to the Stiffness (Displacement) Method. The Stiffness (Displacement) Method

Andrew W. Moore Professor School of Computer Science Carnegie Mellon University

1. State-Space Linear Systems 2. Block Diagrams 3. Exercises

Simplified Identification Scheme for Structures on a Flexible Base

Lecture 8: September 26

08.06 Shooting Method for Ordinary Differential Equations

Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS

A sixth-order dual preserving algorithm for the Camassa-Holm equation

Bumpless transfer for discrete-time switched systems

Section 7.4: Integration of Rational Functions by Partial Fractions

An extremum seeking approach to parameterised loop-shaping control design

FRÉCHET KERNELS AND THE ADJOINT METHOD

BREGMAN ITERATIVE ALGORITHMS FOR l 1 -MINIMIZATION WITH APPLICATIONS TO COMPRESSED SENSING

New Fourth Order Explicit Group Method in the Solution of the Helmholtz Equation Norhashidah Hj. Mohd Ali, Teng Wai Ping

2E1252 Control Theory and Practice

A Macroscopic Traffic Data Assimilation Framework Based on Fourier-Galerkin Method and Minimax Estimation

3. Several Random Variables

Adaptive Predictive Control of Laboratory Heat Exchanger

Essentials of optimal control theory in ECON 4140

DEPTH CONTROL OF THE INFANTE AUV USING GAIN-SCHEDULED REDUCED-ORDER OUTPUT FEEDBACK 1. C. Silvestre A. Pascoal

Step-Size Bounds Analysis of the Generalized Multidelay Adaptive Filter

Math 273b: Calculus of Variations

Assignment Fall 2014

FUZZY BOUNDARY ELEMENT METHODS: A NEW MULTI-SCALE PERTURBATION APPROACH FOR SYSTEMS WITH FUZZY PARAMETERS

FREQUENCY DOMAIN FLUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS

Deakin Research Online

REINFORCEMENT LEARNING AND OPTIMAL ADAPTIVE CONTROL

Execution time certification for gradient-based optimization in model predictive control

Control Performance Monitoring of State-Dependent Nonlinear Processes

A Note on Irreducible Polynomials and Identity Testing

Discussion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli

Convergence analysis of ant colony learning

A Characterization of the Domain of Beta-Divergence and Its Connection to Bregman Variational Model

MULTIOBJECTIVE OPTIMAL STRUCTURAL CONTROL OF THE NOTRE DAME BUILDING MODEL BENCHMARK *

NEURAL CONTROL OF NONLINEAR SYSTEMS: A REFERENCE GOVERNOR APPROACH

Chap. 17 Optimization with constrained conditions...

An Investigation into Estimating Type B Degrees of Freedom

Linearly Solvable Markov Games

3.4-Miscellaneous Equations

AMS 212B Perturbation Methods Lecture 05 Copyright by Hongyun Wang, UCSC

CDS 101/110a: Lecture 8-1 Frequency Domain Design

Conceptual Questions. Problems. 852 CHAPTER 29 Magnetic Fields

RESGen: Renewable Energy Scenario Generation Platform

Integration of Basic Functions. Session 7 : 9/23 1

Distribution Network Planning Based on Entropy Fuzzy Comprehensive

Bridging the Gap Between Multigrid, Hierarchical, and Receding-Horizon Control

A Single Species in One Spatial Dimension

On Multiobjective Duality For Variational Problems

A NEW EVALUATION PLATFORM FOR NAVIGATION SYSTEMS

An optimizing reduced order FDS for the tropical Pacific Ocean reduced gravity model

Geometry of Span (continued) The Plane Spanned by u and v

Transcription:

CDS 110b: Lectre 1-2 Introdction to Optimal Control Richard M. Mrray 4 Janary 2006 Goals: Introdce the problem of optimal control as method of trajectory generation State the maimm principle and give eamples of its sage Reading: Lewis and Syrmos Section 3.1-3.2 available on web page RMM corse notes available on web page Homework #1 Apply the maimm principle to some simple problems De Wed 11 Jan by 5 pm in bo otside 109 Steele 4 Jan 06 R. M. Mrray Caltech 1

Trajectory Generation Trajectory Generation Controller Process Estimator System Dynamics: Trajectory Generation: r d d Tracking Controller: = d + K - d Approach: optimal control Find trajectory that minimies desired cost fnction and satisfies the nonlinear dynamics Use linear state feedback K that minimies qadratic criteria Modern variation: receding horion control Resolve trajectory generation in real-time based on crrent state 4 Jan 06 R. M. Mrray Caltech 2

Review: Optimiation * Optimiation of fnctions: given F:R n R find * R n sch that F * F for all R n Necessary conditions: Remarks: Basic idea: if gradient is non-ero move in that direction to increase cost Conditions are not sfficient conditions: all marked points satisfy conditions only one is global maimm To find the minimm of a fnction F find the maimm of the fnction -F conditions are nchanged so we will switch back and forth MATLAB: X = FMINSEARCHFUNX0 starts at X0 and attempts to find a local minimier X of the fnction FUN 4 Jan 06 R. M. Mrray Caltech 3

Constrained Optimiation * G i = 0 F Optimiation with constraints: given cost fnction F:R n R and constraints G i :R n R i = 1 k find * R n sch that G i * = 0 satisfies constraints and F * F for all sch that G i = 0. Necessary conditions Lagrange mltipliers: i Geometric interpretation: free variables that cancel the gradient of F in directions normal to the constraint Algebraic interpretation: minimie the fnction with respect to leaving free. Variables can take any vale need to choose sch that G = 0 otherwise we can choose to generate a large cost MATLAB: se Optimiation Toolbo 4 Jan 06 R. M. Mrray Caltech 4

Optimal Control Problem Statement: Given nonlinear system Find trajectory * * that satisfies the dynamics and minimies cost Integrated cost Terminal cost Remarks Usally assme we are trying to take system to the origin = 0; easy to generalie by rewriting dynamics in terms of write e = - r Integrated cost L is sed to tradeoff state error from inpt cost Typical cost fnction: qadratic cost 4 Jan 06 R. M. Mrray Caltech 5

Pontryagin s Maimm Principle System: Cost: Hamiltonian: Terminal constraint optional Theorem Pontryagin: If * * is optimal then there eists * t and * sch that and Bondary conditions 2n total Remark: this is a very general and sefl set of conditions 4 Jan 06 R. M. Mrray Caltech 6

Using the Maimm Principle 1. Formlate problem in standard form: 2. Constrct Hamiltonian: 3. Compte necessary conditions: 4. Find the optimal inpt 5. Solve for the optimal trajectory: Sbstitte optimal inpt into necessary conditions and solved bondary vale problem. In general this is hard to do in closed form Good news: can convert this to a comptational problem 4 Jan 06 R. M. Mrray Caltech 7

Eample: Bang-Bang Control 1. Problem formlation: move to origin in minimm amont of time Linear dynamics bonded inpt Minimm time make T small 0 Terminal constraint 2. Constrct Hamiltonian: 3 4. Compte necessary conditions and optimal inpt: Retrns a copy of the dynamics Evoltion of Lagrange mltipliers costates Optimal inpt Remarks: Form of the soltion is easy: apply ma or min inpt at all times Finding actal trajectory is hard: need to search over switching times 4 Jan 06 R. M. Mrray Caltech 8

Application: Dcted Fan flight control Nonlinear design global nonlinearities inpt satration state space constraints Local design ref Trajectory Generation d d noise Plant P Local Control otpt Op Maps NTG MPC+CLF Use real-time trajectory generation to constrct sboptimal feasible trajectories Use model predictive control for reconfigrable tracking & robst performance Etension to mlti-vehicle systems performing cooperative tasks 4 Jan 06 R. M. Mrray Caltech 9

Approach: Receding Horion Optimiation T T t+ T tt+ T = arg min L d + V t+ T [ ] = t = t+ T 0 t f & = f g 0 d Online control cstomiation System: f Constraints/environment: g Mission: L Update in real-time to achieve reconfigrable operation 4 Jan 06 R. M. Mrray Caltech 10

4 Jan 06 R. M. Mrray Caltech 11 Dcted Fan Terrain Avoidance Real-Time Trajectory Generation / Optimiation Collocation Flatness Qasi-collocation f = & i i i i i t t f t t = = & = = = = t i i q q p K & K & K & 0 p q y y y y y y h y K & K & = = =

Eperimental Reslts: Caltech Dcted Fan NTG + MPC + CLF Pitch Control dspace RTOS 4 Jan 06 R. M. Mrray Caltech 12

Eperiments: Caltech Dcted Fan Milam Fran Haser 2002 ACC Dnbar Milam Mrray 2002 IFAC Real-Time MPC on Caltech Dcted Fan Ag 01 NTG with qasi-flat otpts + Lyapnov CLF Average comptation time of ~100 msec Inner pitch loop closed sing local control law; MPC for position variables Inner/oter tradeoff: how mch can be pshed into optimiation 4 Jan 06 R. M. Mrray Caltech 13

Smmary Trajectory Generation Controller Process Estimator Optimal Control for Trajectory Generation Find feasible optimal trajectories for a nonlinear control system Necessary conditions: Pontryagin Maimm Principle Homework Work throgh several eamples of maimm principle Net week: application of optimal control to design state feedback Note: NO CLASS MONDAY - net reglar class on Wed 11 Jan 1:30 pm Corse Project information meeting: Friday 2-3 pm 125 Steele 4 Jan 06 R. M. Mrray Caltech 14