Constrained Optimal Control I
|
|
- Marilynn Gladys Morris
- 5 years ago
- Views:
Transcription
1 Optimal Control, Guidance and Estimation Lecture 34 Constrained Optimal Control I Pro. Radhakant Padhi Dept. o Aerospace Engineering Indian Institute o Science - Bangalore opics Motivation Brie Summary o Unconstrained Optimal Control Pontryagin Minimum Principle ime Optimal Control o LI Systems ime Optimal Control o Double-Integral System Fuel Optimal Control Energy Optimal Control State Constrained optimal control 2
2 Motivation Pro. Radhakant Padhi Dept. o Aerospace Engineering Indian Institute o Science - Bangalore Motivation Physical systems are always restricted by constraints on control and state variables. Examples: hrust delection o the rocket engine cannot not exceed a certain designed value Control surace delections are constrained by hard bounds Aircrats cannot climb beyond a certain altitude (else, they will loose lit because o low dynamic pressure) Robotic arms are constrained by physical limits on angular movements Speed o electric motors should not increase beyond a limit (to prevent wear and tear) Current in a circuit must not increase beyond a limit. Else, some component may burn out. 4
3 Motivation Question: Can these constraints be explicitly handled in the control design? Answer: YES! (optimal control ramework allows that) Ways to handle Sot constrainormulations Hard constrainormulation Problem classiication Control constrained problems State constrained problems Mixed state and control constrained problems 5 Pioneers o Optimal Control 1700s Bernoulli, Newton Euler (Student o Bernoulli) Lagrange years later s Pontryagin Bellman Kalman Bernoulli Pontryagin Euler Bellman Lagrange Kalman Newton 6
4 Lev Semyonovich Pontryagin Lev Semyonovich Pontryagin (September 3,1908: May 3, 1988) Moscow Russia. Lost his eyesight when he was about 14 years old due to an explosion. Entered Moscow State University (1925). In 1930s & 1940s signiicant contribution to topology which was translated into several Languages. As a head o Steklov Mathematical Institute he ocused on general theory o singularly perturbed systems o ordinary dierential equations and maximum principle in optimal control theory. 7 L. S. Pontryagin In 1955, he ormulated a general time-optimal control problem or a ith-order dynamical system describing optimal maneuvers o an aircrat with bounded control unctions. o invent a new calculus o variation he spent three consecutive sleepless nights and came up with the idea o the Hamiltonian ormulation or the problem and the adjoint dierential equations. His other contributions include "singular perturbation theory" and "dierential game theory". He and his co-workers were awarded "Lenin Prize" in
5 Brie Summary o Unconstrained Optimal Control Pro. Radhakant Padhi Dept. o Aerospace Engineering Indian Institute o Science - Bangalore Objective o ind an "admissible" time history o control variable U t, t,, which: ( ) 1) Causes the system governed by Xɺ = t, X, U to ollow an admissible trajectory ( ) 2) Optimizes (minimizes/maximizes) a "meaningul" perormance index = ϕ (, ) + (,, ) J t X L t X U dt 3) Forces the system to satisy "proper boundary conditions". 10
6 Optimal Control Problem Perormance Index (to minimize / maximize): Path Constraint: Xɺ = t X U Boundary Conditions: X = ϕ (, ) + (,, ) J t X L t X U dt (,, ) t ( ) t 0 = X : Speciied 0 ( ) : Fixed, X t : Free 11 Necessary Conditions o Optimality Augmented PI Hamiltonian First Variation ( ɺ ) J = ϕ + L + λ X dt ( + λ ) H L t ( ɺ ) δ J = δϕ + δ H λ X dt ( ɺ ) = δϕ + δ H λ X dt 12
7 Necessary Conditions o Optimality First Variation Individual terms δϕ ( t, X ) ( δ X ) ( ɺ ɺ ) δ J = δϕ + δ H δλ X λ δ X dt ϕ = X δ H H H H ( t, X, U, λ ) ( δ X ) ( U ) ( ) X δ U δλ = + + λ 13 Necessary Conditions o Optimality d ( δ X ) λ δ Xɺ dt λ dt ( ) = dt, δ X dλ = λ δ X, δ X 0 t dt ( ) ɺ 0 X ( X ) ɺ dt δ X dt = λ δ X λ δ X δ X λ dt = λ δ δ λ 14
8 Necessary Conditions o Optimality First Variation ϕ δ J = ( δ X ) ( δ X ) λ X H H H + ( δ X ) + ( δu ) + ( δλ ) dt X U λ ( δ ) ɺ λ ( δλ ) + X dt Xɺ dt 15 Necessary Conditions o Optimality First Variation ϕ δ J = ( δ X ) λ X = 0 ( δ ) ɺ λ ( δ ) H H + X + dt + U dt X U H + ( δλ ) X ɺ dt λ 16
9 Necessary Conditions o Optimality: Summary State Equation Costate Equation H Xɺ = = t X U λ ɺ H λ = X (,, ) Optimal Control Equation H U = 0 Boundary Condition λ ϕ = ( ) X X t X Fixed = : Necessary Conditions o Optimality: Some Comments State and Costate equations are dynamic equations. I one is stable, the other turns out to be unstable! Optimal control equation is a stationary equation Boundary conditions are split: it leads to wo-point- Boundary-Value Problem (PBVP) State equation develops orward whereas Costate equation develops backwards. It is known as Curse o Complexity in optimal control raditionally, PBVPs demand computationally-intensive iterative numerical procedures, which lead to open-loop control structure. 18
10 Control Constrained Problems: Pontryagin Minimum Principle Reerence: D. S. Naidu: Optimal Control Systems, CRC Press, Pontryagin Pro. Radhakant Padhi Dept. o Aerospace Engineering Indian Institute o Science - Bangalore Objective o ind an "admissible" time history o control variable U t, t, t, + ( or, component wise, j j j ) where U ( t) U U u ( t) U,which: ( ) 1) Causes the system governed by Xɺ = t, X, U to ollow an admissible trajectory 2) Optimizes (minimizes/maximizes) a "meanigul" perormance index = ϕ (, ) + (,, ) J t X L t X U dt 3) Forces the system to satisy "proper boundary conditions". ( ) 20
11 Optimum o Control Functional u, u ( ) ( ) ( ) Variation: δ u t = u t u t u ( t) δu( t) u ( t) : Optimum control a b t 21 Optimum o Control Functional ( ) ( ) u( t) u( t) u ( t) A unctional u t is said to have a relative optimum at u t, i ε > 0 such thaor all unctions Ω which satisy < ε, the increment o J has the "same sign". ( ) ( ) ( ) 1) I J = J u J u 0, then J u is a relative (local) "Minimum". ( ) ( ) ( ) 2) I J = J u J u 0, then J u is a relative (local) "Maximum". Note: I the above relationships are satisied or arbitrarily large ε > 0, then ( ) J u is a "global optimum". 22
12 Pontryagin Minimum Principle U U δu With variations in control = +, ( δ ) ( ) ( ) = δ J ( U δu ) J U, U = J U J U 0 or Minimum, + HO J δu (Neglecting HO) U However, when U ( t) U, δu is no longer arbitrary or all t, 23 Control Constrained Problems Reerence: D. S. Naidu: Optimal Control Systems, CRC Press, Note: he condition δ J = 0 is valid only i u (t) lies within the boundary (i.e. it has no constraint) or the entire time interval t, 0 t 24
13 Pontryagin Minimum Principle U ( δ ) Necessary Condition: δ J U ( t), U ( t) 0 First Variation: ( t) : Optimal solution δu t U t ( ) : Allowable variation about ( ) H H H δ J = + ɺ λ δ X + δ U + Xɺ δλ dt X U λ ϕ + λ X δ X 25 Pontryagin Minimum Principle 1) In control constrained problems, variations in costates δλ( t) can be arbitrary. his gives H Xɺ = 0 λ H Xɺ = = ( t, X, U ) (state equation) λ 2) I the costate λ( t) is selected such that the coeicient o δ X ( t) is 0 (i.e. variations in states can be arbitrary), then ɺ H λ + = 0 X ɺ H λ = (costate equation) X 26
14 Pontryagin Minimum Principle 3) Boundary conditions are not eected by the control constraints. Hence, the ollowing ransversality condition still holds good. ϕ λ = X With the above observations, the necessary condition becomes H δ J ( U, δu ) = δu dt U (, δ, λ ) (,, λ ) = H X U U H X U + dt 0 admissible δu arbitrarily small 27 Pontryagin Minimum Principle Since δu ( t) is arbitrarily small, the integrand 0. his gives us H X, U + δu, λ H X, U, U i. e. (, U, λ ) H ( X, U, λ ) H X U ( t) U ( λ ) "Necessary condition" or constrained optimal control U is given by min H ( X, U, λ ) = H ( X, U, λ ) i. e. the optimal control should minimize the Hamiltonian his is known as the "Pontryagin Minimum Principle". 28
15 Solution Procedure o a given Problem Hamiltonian : Necessary Conditions : (,, λ) = (, ) + λ (, ) H X U L X U X U H (i) State Equation: Xɺ = = ( t, X, U ) λ H (ii) Costate Equation: ɺ λ = X (ii i) Optimal Control Equation: Minimize H with repect to U ( t) U ( 0) ( λ) H ( X U λ ) i.e. H X, U,,, (iv) Boundary conditions: X ϕ = Speciied, λ = X 29 Some Important Observations 1) he optimality condition (,, λ ) H ( X, U, λ ) H X U is valid or both constrained and unconstrained control system, ( H U ) whereas the control relation / = 0 is valid or unconstrained systems only. 2) he results given above provide the necessary conditions only. 3) he suicient condition or unconstrained control problem is that 2 H 2 should be positive deinite matrix t, U ( X, U, λ ) 30
16 A Simple Scalar Algebraic Example Problem : 2 Minimize the unction H = u 6u + 7 subject to the constraint relation u 2, i. e. 2 u + 2 Solution : Using the relation or unconstrained control, H = 2u 6 = 0 u u = 3 (Not admissible!) 31 Plot o Hamiltonian Reerence: D. S. Naidu: Optimal Control Systems, CRC Press,
17 A Simple Scalar Example In this case, the admissible optimal value is u = + 2 can also be obtained rom static optimization results using Kharush-Kuhn-ucker conditions. Note : I the constraint had been u 3, i. e. 3 u + 3, then either o the relation could be used and obtain the optimal value as u = 3. However, unortunately many practical constraints do not admit such solutions! 33 Additional Necessary Conditions (Due to Pontryagin & Co-workers) 1) I the inal time t is " ixed" and the Hamiltonian H does not depend on time t explicitly, then the Hamiltonian must be constant along the optimal trajectory, i.e. H t t t = Constant 0, 2) I the inal time t is " ree" and the Hamiltonian does not depend on time t explicitly, then the the Hamiltonian must be identically zero along the optimal trajectory, i.e. H = 0 t, t 34
18 Proo or Unconstrained Problem heorem: I the Hamiltonian H is not an expliciunction o time, then H is constant along the optimal path. Proo: dh H H H H H X U But X and X X = + ɺ + ɺ + ɺ λ = ɺ ɺ λ ɺ = ɺ ɺ λ dt t X U λ λ H H H = + Xɺ + ɺ λ + Uɺ t X U 0 0 dh H = ( on optimal path) dt t = 0 i H is not an expliciunction o t. Hence, the result! ( ) 35 Conclusions Physical systems are always restricted by constraints on control and state variables. In this class we studied about Brie Summary o Unconstrained Optimal Control Pontryagin Minimum Principle or Control Constrained Optimal Control (in a generic sense) In the next two classes, we will study about ime Optimal Control o LI Systems Fuel Optimal Control Energy Optimal Control State Constrained optimal control 36
19 Reerences D. S. Naidu: Optimal Control Systems, CRC Press, L. M. Hocking: Optimal Control: An Introduction to heory and Applications, Oxord University Press, New York, NY, L. S. Pontryagin, V. G. Boltyanskii, R. V. Gamkrelidze and E. F. Mishchenko: he Mathematical heory o Optimal Processes, Wiley-Intersciences, New York, NY, 1962 (ranslated rom Russian). 37 hanks or the Attention.!! 38
Constrained Optimal Control. Constrained Optimal Control II
Optimal Control, Guidance and Estimation Lecture 35 Constrained Optimal Control II Prof. Radhakant Padhi Dept. of Aerospace Engineering Indian Institute of Science - Bangalore opics: Constrained Optimal
More informationLinear Quadratic Regulator (LQR) I
Optimal Control, Guidance and Estimation Lecture Linear Quadratic Regulator (LQR) I Pro. Radhakant Padhi Dept. o Aerospace Engineering Indian Institute o Science - Bangalore Generic Optimal Control Problem
More informationLinear Quadratic Regulator (LQR) II
Optimal Control, Guidance and Estimation Lecture 11 Linear Quadratic Regulator (LQR) II Pro. Radhakant Padhi Dept. o Aerospace Engineering Indian Institute o Science - Bangalore Outline Summary o LQR design
More informationLinear Quadratic Regulator (LQR) Design I
Lecture 7 Linear Quadratic Regulator LQR) Design I Dr. Radhakant Padhi Asst. Proessor Dept. o Aerospace Engineering Indian Institute o Science - Bangalore LQR Design: Problem Objective o drive the state
More informationOptimal Control. with. Aerospace Applications. James M. Longuski. Jose J. Guzman. John E. Prussing
Optimal Control with Aerospace Applications by James M. Longuski Jose J. Guzman John E. Prussing Published jointly by Microcosm Press and Springer 2014 Copyright Springer Science+Business Media New York
More informationSWEEP METHOD IN ANALYSIS OPTIMAL CONTROL FOR RENDEZ-VOUS PROBLEMS
J. Appl. Math. & Computing Vol. 23(2007), No. 1-2, pp. 243-256 Website: http://jamc.net SWEEP METHOD IN ANALYSIS OPTIMAL CONTROL FOR RENDEZ-VOUS PROBLEMS MIHAI POPESCU Abstract. This paper deals with determining
More informationOptimal Control Design
Optimal Control Design Prof. Radhakant Padhi Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Acknowledgement: Indian Institute of Science Founded in 1909 more than 100 years old
More informationClassical Numerical Methods to Solve Optimal Control Problems
Lecture 26 Classical Numerical Methods to Solve Optimal Control Problems Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Necessary Conditions
More information8. INTRODUCTION TO STATISTICAL THERMODYNAMICS
n * D n d Fluid z z z FIGURE 8-1. A SYSTEM IS IN EQUILIBRIUM EVEN IF THERE ARE VARIATIONS IN THE NUMBER OF MOLECULES IN A SMALL VOLUME, SO LONG AS THE PROPERTIES ARE UNIFORM ON A MACROSCOPIC SCALE 8. INTRODUCTION
More informationA Simple Explanation of the Sobolev Gradient Method
A Simple Explanation o the Sobolev Gradient Method R. J. Renka July 3, 2006 Abstract We have observed that the term Sobolev gradient is used more oten than it is understood. Also, the term is oten used
More informationThe Ascent Trajectory Optimization of Two-Stage-To-Orbit Aerospace Plane Based on Pseudospectral Method
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 00 (014) 000 000 www.elsevier.com/locate/procedia APISAT014, 014 Asia-Paciic International Symposium on Aerospace Technology,
More informationDynamic Inversion Design II
Lecture 32 Dynamic Inversion Design II Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Topics Summary of Dynamic Inversion Design Advantages Issues
More informationUsing Lyapunov Theory I
Lecture 33 Stability Analysis of Nonlinear Systems Using Lyapunov heory I Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Outline Motivation Definitions
More informationEE291E/ME 290Q Lecture Notes 8. Optimal Control and Dynamic Games
EE291E/ME 290Q Lecture Notes 8. Optimal Control and Dynamic Games S. S. Sastry REVISED March 29th There exist two main approaches to optimal control and dynamic games: 1. via the Calculus of Variations
More informationThree-Dimensional Trajectory Optimization in Constrained Airspace
Three-Dimensional Traectory Optimization in Constrained Airspace Ran Dai * and John E. Cochran, Jr. Auburn University, Auburn, Alabama 36849 The operational airspace o aerospace vehicles, including airplanes
More informationStatic and Dynamic Optimization (42111)
Static and Dynamic Optimization (421) Niels Kjølstad Poulsen Build. 0b, room 01 Section for Dynamical Systems Dept. of Applied Mathematics and Computer Science The Technical University of Denmark Email:
More informationMinimum-Time Trajectory Optimization of Multiple Revolution Low-Thrust Earth-Orbit Transfers
Minimum-Time Trajectory Optimization o Multiple Revolution Low-Thrust Earth-Orbit Transers Kathryn F. Graham Anil V. Rao University o Florida Gainesville, FL 32611-625 Abstract The problem o determining
More informationFeedback Optimal Control of Low-thrust Orbit Transfer in Central Gravity Field
Vol. 4, No. 4, 23 Feedback Optimal Control of Low-thrust Orbit Transfer in Central Gravity Field Ashraf H. Owis Department of Astronomy, Space and Meteorology, Faculty of Science, Cairo University Department
More informationPRECISION ZEM/ZEV FEEDBACK GUIDANCE ALGORITHM UTILIZING VINTI S ANALYTIC SOLUTION OF PERTURBED KEPLER PROBLEM
AAS 16-345 PRECISION ZEM/ZEV FEEDBACK GUIDANCE ALGORITHM UTILIZING VINTI S ANALYTIC SOLUTION OF PERTURBED KEPLER PROBLEM Jaemyung Ahn, * Yanning Guo, and Bong Wie A new implementation o a zero-eort-miss/zero-eort-velocity
More informationFeedback Optimal Control for Inverted Pendulum Problem by Using the Generating Function Technique
(IJACSA) International Journal o Advanced Computer Science Applications Vol. 5 No. 11 14 Feedback Optimal Control or Inverted Pendulum Problem b Using the Generating Function echnique Han R. Dwidar Astronom
More informationOptimal Control, Guidance and Estimation. Lecture 16. Overview of Flight Dynamics II. Prof. Radhakant Padhi. Prof. Radhakant Padhi
Optimal Control, Guidance and Estimation Lecture 16 Overview of Flight Dynamics II Prof. Radhakant Padhi Dept. of erospace Engineering Indian Institute of Science - Bangalore Point Mass Dynamics Prof.
More informationLinear Quadratic Regulator (LQR) Design II
Lecture 8 Linear Quadratic Regulator LQR) Design II Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Outline Stability and Robustness properties
More information18-660: Numerical Methods for Engineering Design and Optimization
8-66: Numerical Methods or Engineering Design and Optimization Xin Li Department o ECE Carnegie Mellon University Pittsburgh, PA 53 Slide Overview Linear Regression Ordinary least-squares regression Minima
More informationModule 2 Selection of Materials and Shapes. IIT, Bombay
Module Selection o Materials and Shapes Lecture Selection o Materials - I Instructional objectives By the end o this lecture, the student will learn (a) what is a material index and how does it help in
More informationSpatial Vector Algebra
A Short Course on The Easy Way to do Rigid Body Dynamics Roy Featherstone Dept. Inormation Engineering, RSISE The Australian National University Spatial vector algebra is a concise vector notation or describing
More informationLecture 25: Heat and The 1st Law of Thermodynamics Prof. WAN, Xin
General Physics I Lecture 5: Heat and he 1st Law o hermodynamics Pro. WAN, Xin xinwan@zju.edu.cn http://zimp.zju.edu.cn/~xinwan/ Latent Heat in Phase Changes Latent Heat he latent heat o vaporization or
More informationYURI LEVIN AND ADI BEN-ISRAEL
Pp. 1447-1457 in Progress in Analysis, Vol. Heinrich G W Begehr. Robert P Gilbert and Man Wah Wong, Editors, World Scientiic, Singapore, 003, ISBN 981-38-967-9 AN INVERSE-FREE DIRECTIONAL NEWTON METHOD
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 informationIntroduction to Multicopter Design and Control
Introduction to Multicopter Design and Control Lesson 10 Stability and Controllability Quan Quan, Associate Proessor qq_buaa@buaa.edu.cn BUAA Reliable Flight Control Group, http://rly.buaa.edu.cn/ Beihang
More informationFeedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 2016.
Feedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 216 Ravi N Banavar banavar@iitbacin September 24, 216 These notes are based on my readings o the two books Nonlinear Control
More information1 Relative degree and local normal forms
THE ZERO DYNAMICS OF A NONLINEAR SYSTEM 1 Relative degree and local normal orms The purpose o this Section is to show how single-input single-output nonlinear systems can be locally given, by means o a
More informationGeneral solution of the Inhomogeneous Div-Curl system and Consequences
General solution o the Inhomogeneous Div-Curl system and Consequences Briceyda Berenice Delgado López Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional March 07 2017 CINVESTAV
More information2.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics. References - geostatistics. References geostatistics (cntd.
.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics Spline interpolation was originally developed or image processing. In GIS, it is mainly used in visualization o spatial
More informationAn Ensemble Kalman Smoother for Nonlinear Dynamics
1852 MONTHLY WEATHER REVIEW VOLUME 128 An Ensemble Kalman Smoother or Nonlinear Dynamics GEIR EVENSEN Nansen Environmental and Remote Sensing Center, Bergen, Norway PETER JAN VAN LEEUWEN Institute or Marine
More informationLecture : Feedback Linearization
ecture : Feedbac inearization Niola Misovic, dipl ing and Pro Zoran Vuic June 29 Summary: This document ollows the lectures on eedbac linearization tought at the University o Zagreb, Faculty o Electrical
More informationChapter 2 Lecture 7 Longitudinal stick fixed static stability and control 4 Topics
hapter 2 Lecture 7 Longitudinal stick ixed static stability and control 4 Topics 2.4.6 Revised expression or mcgt 2.4.7 mαt in stick-ixed case 2.5 ontributions o uselage to mcg and mα 2.5.1 ontribution
More informationAC&ST AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS. Claudio Melchiorri
C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI)
More informationg(x,y) = c. For instance (see Figure 1 on the right), consider the optimization problem maximize subject to
1 of 11 11/29/2010 10:39 AM From Wikipedia, the free encyclopedia In mathematical optimization, the method of Lagrange multipliers (named after Joseph Louis Lagrange) provides a strategy for finding the
More informationChapter 6 Reliability-based design and code developments
Chapter 6 Reliability-based design and code developments 6. General Reliability technology has become a powerul tool or the design engineer and is widely employed in practice. Structural reliability analysis
More informationSupplementary material for Continuous-action planning for discounted infinite-horizon nonlinear optimal control with Lipschitz values
Supplementary material or Continuous-action planning or discounted ininite-horizon nonlinear optimal control with Lipschitz values List o main notations x, X, u, U state, state space, action, action space,
More informationAnalysis of the regularity, pointwise completeness and pointwise generacy of descriptor linear electrical circuits
Computer Applications in Electrical Engineering Vol. 4 Analysis o the regularity pointwise completeness pointwise generacy o descriptor linear electrical circuits Tadeusz Kaczorek Białystok University
More informationDIFFERENTIAL POLYNOMIALS GENERATED BY SECOND ORDER LINEAR DIFFERENTIAL EQUATIONS
Journal o Applied Analysis Vol. 14, No. 2 2008, pp. 259 271 DIFFERENTIAL POLYNOMIALS GENERATED BY SECOND ORDER LINEAR DIFFERENTIAL EQUATIONS B. BELAÏDI and A. EL FARISSI Received December 5, 2007 and,
More informationGlobal Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance
Global Weak Solution o Planetary Geostrophic Equations with Inviscid Geostrophic Balance Jian-Guo Liu 1, Roger Samelson 2, Cheng Wang 3 Communicated by R. Temam) Abstract. A reormulation o the planetary
More informationAUGMENTED POLYNOMIAL GUIDANCE FOR TERMINAL VELOCITY CONSTRAINTS
AUGMENTED POLYNOMIAL GUIDANCE FOR TERMINAL VELOCITY CONSTRAINTS Gun-Hee Moon*, Sang-Woo Shim*, and Min-Jea Tah* *Korea Advanced Institute Science and Technology Keywords: Polynomial guidance, terminal
More informationA new Control Strategy for Trajectory Tracking of Fire Rescue Turntable Ladders
Proceedings o the 17th World Congress The International Federation o Automatic Control Seoul, Korea, July 6-11, 28 A new Control Strategy or Trajectory Tracking o Fire Rescue Turntable Ladders N. Zimmert,
More informationMath 248B. Base change morphisms
Math 248B. Base change morphisms 1. Motivation A basic operation with shea cohomology is pullback. For a continuous map o topological spaces : X X and an abelian shea F on X with (topological) pullback
More informationLecture 6 Classical Control Overview IV. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore
Lecture 6 Classical Control Overview IV Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Lead Lag Compensator Design Dr. Radhakant Padhi Asst.
More informationClassification of effective GKM graphs with combinatorial type K 4
Classiication o eective GKM graphs with combinatorial type K 4 Shintarô Kuroki Department o Applied Mathematics, Faculty o Science, Okayama Uniervsity o Science, 1-1 Ridai-cho Kita-ku, Okayama 700-0005,
More informationIonosphere Ray-Tracing of RF Signals and Solution Sensitivities to Model Parameters
Ionosphere Ray-Tracing o RF Signals and Solution Sensitivities to Model Parameters Mark L. Psiaki Virginia Tech, Blacksburg, VA 2461, USA BIOGRAPHY Mark L. Psiaki is Proessor and Kevin Croton Faculty Chair
More informationENSC327 Communications Systems 2: Fourier Representations. School of Engineering Science Simon Fraser University
ENSC37 Communications Systems : Fourier Representations School o Engineering Science Simon Fraser University Outline Chap..5: Signal Classiications Fourier Transorm Dirac Delta Function Unit Impulse Fourier
More informationMicro-canonical ensemble model of particles obeying Bose-Einstein and Fermi-Dirac statistics
Indian Journal o Pure & Applied Physics Vol. 4, October 004, pp. 749-757 Micro-canonical ensemble model o particles obeying Bose-Einstein and Fermi-Dirac statistics Y K Ayodo, K M Khanna & T W Sakwa Department
More informationPulling by Pushing, Slip with Infinite Friction, and Perfectly Rough Surfaces
1993 IEEE International Conerence on Robotics and Automation Pulling by Pushing, Slip with Ininite Friction, and Perectly Rough Suraces Kevin M. Lynch Matthew T. Mason The Robotics Institute and School
More informationRESOLUTION MSC.362(92) (Adopted on 14 June 2013) REVISED RECOMMENDATION ON A STANDARD METHOD FOR EVALUATING CROSS-FLOODING ARRANGEMENTS
(Adopted on 4 June 203) (Adopted on 4 June 203) ANNEX 8 (Adopted on 4 June 203) MSC 92/26/Add. Annex 8, page THE MARITIME SAFETY COMMITTEE, RECALLING Article 28(b) o the Convention on the International
More informationMaximum Flow. Reading: CLRS Chapter 26. CSE 6331 Algorithms Steve Lai
Maximum Flow Reading: CLRS Chapter 26. CSE 6331 Algorithms Steve Lai Flow Network A low network G ( V, E) is a directed graph with a source node sv, a sink node tv, a capacity unction c. Each edge ( u,
More informationNumerical Methods - Lecture 2. Numerical Methods. Lecture 2. Analysis of errors in numerical methods
Numerical Methods - Lecture 1 Numerical Methods Lecture. Analysis o errors in numerical methods Numerical Methods - Lecture Why represent numbers in loating point ormat? Eample 1. How a number 56.78 can
More informationPhysics 5153 Classical Mechanics. Solution by Quadrature-1
October 14, 003 11:47:49 1 Introduction Physics 5153 Classical Mechanics Solution by Quadrature In the previous lectures, we have reduced the number o eective degrees o reedom that are needed to solve
More informationUnderwater vehicles: a surprising non time-optimal path
Underwater vehicles: a surprising non time-optimal path M. Chyba Abstract his paper deals with the time-optimal problem for a class of underwater vehicles. We prove that if two configurations at rest can
More informationEstimation of Sample Reactivity Worth with Differential Operator Sampling Method
Progress in NUCLEAR SCIENCE and TECHNOLOGY, Vol. 2, pp.842-850 (2011) ARTICLE Estimation o Sample Reactivity Worth with Dierential Operator Sampling Method Yasunobu NAGAYA and Takamasa MORI Japan Atomic
More informationNumerical Analysis II. Problem Sheet 8
P. Grohs S. Hosseini Ž. Kereta Spring Term 15 Numerical Analysis II ETH Zürich D-MATH Problem Sheet 8 Problem 8.1 The dierential equation Autonomisation and Strang-Splitting ẏ = A(t)y, y(t ) = y, A(t)
More informationMultibody simulation
Multibody simulation Dynamics of a multibody system (Newton-Euler formulation) Dimitar Dimitrov Örebro University June 8, 2012 Main points covered Newton-Euler formulation forward dynamics inverse dynamics
More informationSOME CHARACTERIZATIONS OF HARMONIC CONVEX FUNCTIONS
International Journal o Analysis and Applications ISSN 2291-8639 Volume 15, Number 2 2017, 179-187 DOI: 10.28924/2291-8639-15-2017-179 SOME CHARACTERIZATIONS OF HARMONIC CONVEX FUNCTIONS MUHAMMAD ASLAM
More informationFluctuationlessness Theorem and its Application to Boundary Value Problems of ODEs
Fluctuationlessness Theorem and its Application to Boundary Value Problems o ODEs NEJLA ALTAY İstanbul Technical University Inormatics Institute Maslak, 34469, İstanbul TÜRKİYE TURKEY) nejla@be.itu.edu.tr
More informationOn a Closed Formula for the Derivatives of e f(x) and Related Financial Applications
International Mathematical Forum, 4, 9, no. 9, 41-47 On a Closed Formula or the Derivatives o e x) and Related Financial Applications Konstantinos Draais 1 UCD CASL, University College Dublin, Ireland
More informationScattered Data Approximation of Noisy Data via Iterated Moving Least Squares
Scattered Data Approximation o Noisy Data via Iterated Moving Least Squares Gregory E. Fasshauer and Jack G. Zhang Abstract. In this paper we ocus on two methods or multivariate approximation problems
More informationX-ray Diffraction. Interaction of Waves Reciprocal Lattice and Diffraction X-ray Scattering by Atoms The Integrated Intensity
X-ray Diraction Interaction o Waves Reciprocal Lattice and Diraction X-ray Scattering by Atoms The Integrated Intensity Basic Principles o Interaction o Waves Periodic waves characteristic: Frequency :
More informationComments on Magnetohydrodynamic Unsteady Flow of A Non- Newtonian Fluid Through A Porous Medium
Comments on Magnetohydrodynamic Unsteady Flow o A Non- Newtonian Fluid Through A Porous Medium Mostaa A.A.Mahmoud Department o Mathematics, Faculty o Science, Benha University (358), Egypt Abstract The
More informationIntegrator Backstepping using Barrier Functions for Systems with Multiple State Constraints
Integrator Backstepping using Barrier Functions for Systems with Multiple State Constraints Khoi Ngo Dep. Engineering, Australian National University, Australia Robert Mahony Dep. Engineering, Australian
More informationGROWTH PROPERTIES OF COMPOSITE ANALYTIC FUNCTIONS OF SEVERAL COMPLEX VARIABLES IN UNIT POLYDISC FROM THE VIEW POINT OF THEIR NEVANLINNA L -ORDERS
Palestine Journal o Mathematics Vol 8209 346 352 Palestine Polytechnic University-PPU 209 GROWTH PROPERTIES OF COMPOSITE ANALYTIC FUNCTIONS OF SEVERAL COMPLEX VARIABLES IN UNIT POLYDISC FROM THE VIEW POINT
More informationChapter 5. Pontryagin s Minimum Principle (Constrained OCP)
Chapter 5 Pontryagin s Minimum Principle (Constrained OCP) 1 Pontryagin s Minimum Principle Plant: (5-1) u () t U PI: (5-2) Boundary condition: The goal is to find Optimal Control. 2 Pontryagin s Minimum
More information2. ETA EVALUATIONS USING WEBER FUNCTIONS. Introduction
. ETA EVALUATIONS USING WEBER FUNCTIONS Introduction So ar we have seen some o the methods or providing eta evaluations that appear in the literature and we have seen some o the interesting properties
More informationProbabilistic Observations and Valuations (Extended Abstract) 1
Replace this ile with prentcsmacro.sty or your meeting, or with entcsmacro.sty or your meeting. Both can be ound at the ENTCS Macro Home Page. Probabilistic Observations and Valuations (Extended Abstract)
More informationStrong Lyapunov Functions for Systems Satisfying the Conditions of La Salle
06 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 49, NO. 6, JUNE 004 Strong Lyapunov Functions or Systems Satisying the Conditions o La Salle Frédéric Mazenc and Dragan Ne sić Abstract We present a construction
More informationLegendre Transforms, Calculus of Varations, and Mechanics Principles
page 437 Appendix C Legendre Transforms, Calculus of Varations, and Mechanics Principles C.1 Legendre Transforms Legendre transforms map functions in a vector space to functions in the dual space. From
More informationCircuit Complexity / Counting Problems
Lecture 5 Circuit Complexity / Counting Problems April 2, 24 Lecturer: Paul eame Notes: William Pentney 5. Circuit Complexity and Uniorm Complexity We will conclude our look at the basic relationship between
More informationAnalysis Scheme in the Ensemble Kalman Filter
JUNE 1998 BURGERS ET AL. 1719 Analysis Scheme in the Ensemble Kalman Filter GERRIT BURGERS Royal Netherlands Meteorological Institute, De Bilt, the Netherlands PETER JAN VAN LEEUWEN Institute or Marine
More informationNumerical Solution of Ordinary Differential Equations in Fluctuationlessness Theorem Perspective
Numerical Solution o Ordinary Dierential Equations in Fluctuationlessness Theorem Perspective NEJLA ALTAY Bahçeşehir University Faculty o Arts and Sciences Beşiktaş, İstanbul TÜRKİYE TURKEY METİN DEMİRALP
More informationTour Planning for an Unmanned Air Vehicle under Wind Conditions
Tour Planning or an Unmanned Air Vehicle under Wind Conditions Rachelle L. McNeely and Ram V. Iyer Texas Tech University, Lubbock, TX 79409-1042 and Phillip R. Chandler U. S. Air Force Research Laboratory,
More informationEXISTENCE OF SOLUTIONS TO SYSTEMS OF EQUATIONS MODELLING COMPRESSIBLE FLUID FLOW
Electronic Journal o Dierential Equations, Vol. 15 (15, No. 16, pp. 1 8. ISSN: 17-6691. URL: http://ejde.math.txstate.e or http://ejde.math.unt.e tp ejde.math.txstate.e EXISTENCE OF SOLUTIONS TO SYSTEMS
More informationProbabilistic Optimisation applied to Spacecraft Rendezvous on Keplerian Orbits
Probabilistic Optimisation applied to pacecrat Rendezvous on Keplerian Orbits Grégory aive a, Massimiliano Vasile b a Université de Liège, Faculté des ciences Appliquées, Belgium b Dipartimento di Ingegneria
More informationMinimum Energy Control of Fractional Positive Electrical Circuits with Bounded Inputs
Circuits Syst Signal Process (26) 35:85 829 DOI.7/s34-5-8-7 Minimum Energy Control o Fractional Positive Electrical Circuits with Bounded Inputs Tadeusz Kaczorek Received: 24 April 25 / Revised: 5 October
More informationRobust and Optimal Control of Uncertain Dynamical Systems with State-Dependent Switchings Using Interval Arithmetic
Robust and Optimal Control o Uncertain Dynamical Systems with State-Dependent Switchings Using Interval Arithmetic Andreas Rauh Chair o Mechatronics, University o Rostock, D-1859 Rostock, Germany Andreas.Rauh@uni-rostock.de
More informationPercentile Policies for Inventory Problems with Partially Observed Markovian Demands
Proceedings o the International Conerence on Industrial Engineering and Operations Management Percentile Policies or Inventory Problems with Partially Observed Markovian Demands Farzaneh Mansouriard Department
More informationCategories and Natural Transformations
Categories and Natural Transormations Ethan Jerzak 17 August 2007 1 Introduction The motivation or studying Category Theory is to ormalise the underlying similarities between a broad range o mathematical
More informationOBSERVER/KALMAN AND SUBSPACE IDENTIFICATION OF THE UBC BENCHMARK STRUCTURAL MODEL
OBSERVER/KALMAN AND SUBSPACE IDENTIFICATION OF THE UBC BENCHMARK STRUCTURAL MODEL Dionisio Bernal, Burcu Gunes Associate Proessor, Graduate Student Department o Civil and Environmental Engineering, 7 Snell
More informationApplication of singular perturbation theory in modeling and control of flexible robot arm
Research Article International Journal of Advanced Technology and Engineering Exploration, Vol 3(24) ISSN (Print): 2394-5443 ISSN (Online): 2394-7454 http://dx.doi.org/10.19101/ijatee.2016.324002 Application
More informationFLUID MECHANICS. Lecture 7 Exact solutions
FLID MECHANICS Lecture 7 Eact solutions 1 Scope o Lecture To present solutions or a ew representative laminar boundary layers where the boundary conditions enable eact analytical solutions to be obtained.
More informationSolution of the Synthesis Problem in Hilbert Spaces
Solution o the Synthesis Problem in Hilbert Spaces Valery I. Korobov, Grigory M. Sklyar, Vasily A. Skoryk Kharkov National University 4, sqr. Svoboda 677 Kharkov, Ukraine Szczecin University 5, str. Wielkopolska
More informationNumerical approximation for optimal control problems via MPC and HJB. Giulia Fabrini
Numerical approximation for optimal control problems via MPC and HJB Giulia Fabrini Konstanz Women In Mathematics 15 May, 2018 G. Fabrini (University of Konstanz) Numerical approximation for OCP 1 / 33
More informationMathematical Theory of Control Systems Design
Mathematical Theory of Control Systems Design by V. N. Afarias'ev, V. B. Kolmanovskii and V. R. Nosov Moscow University of Electronics and Mathematics, Moscow, Russia KLUWER ACADEMIC PUBLISHERS DORDRECHT
More informationGain Scheduling and Dynamic Inversion
Lecture 31 Gain Scheduling and Dynamic Inversion Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Topics A Brief Overview of Gain Scheduling Philosophy
More informationExtreme Values of Functions
Extreme Values o Functions When we are using mathematics to model the physical world in which we live, we oten express observed physical quantities in terms o variables. Then, unctions are used to describe
More informationLecture 5 Classical Control Overview III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore
Lecture 5 Classical Control Overview III Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore A Fundamental Problem in Control Systems Poles of open
More informationNONPARAMETRIC PREDICTIVE INFERENCE FOR REPRODUCIBILITY OF TWO BASIC TESTS BASED ON ORDER STATISTICS
REVSTAT Statistical Journal Volume 16, Number 2, April 2018, 167 185 NONPARAMETRIC PREDICTIVE INFERENCE FOR REPRODUCIBILITY OF TWO BASIC TESTS BASED ON ORDER STATISTICS Authors: Frank P.A. Coolen Department
More informationMinimum-Fuel Trajectory Optimization of Many Revolution Low-Thrust Earth-Orbit Transfers
Minimum-Fuel Trajectory Optimization o Many Revolution Low-Thrust Earth-Orbit Transers Kathryn F. Graham 1, Anil V. Rao 2 University o Florida, Gainesville, FL 32611-625 Abstract The problem o determining
More informationRobin-Robin preconditioned Krylov methods for fluid-structure interaction problems
MOX Report No. 19/2008 Robin-Robin preconditioned Krylov methods or luid-structure interaction problems Santiago Badia, Fabio Nobile, Christian Vergara MOX, Dipartimento di Matematica F. Brioschi Politecnico
More informationExtremal Trajectories for Bounded Velocity Differential Drive Robots
Extremal Trajectories for Bounded Velocity Differential Drive Robots Devin J. Balkcom Matthew T. Mason Robotics Institute and Computer Science Department Carnegie Mellon University Pittsburgh PA 523 Abstract
More informationUltra Fast Calculation of Temperature Profiles of VLSI ICs in Thermal Packages Considering Parameter Variations
Ultra Fast Calculation o Temperature Proiles o VLSI ICs in Thermal Packages Considering Parameter Variations Je-Hyoung Park, Virginia Martín Hériz, Ali Shakouri, and Sung-Mo Kang Dept. o Electrical Engineering,
More informationOptimal Control of Nonlinear Systems using RBF Neural Network and Adaptive Extended Kalman Filter
9 American Control Conerence Hyatt Regency Riverront, St. Louis, MO, USA June -, 9 WeA. Optimal Control o Nonlinear Systems using RBF Neural Network and Adaptive Extended Kalman Filter Peda V. Medagam,
More informationSymplectic Möbius integrators for LQ optimal control problems
Symplectic Möbius integrators or LQ optimal control problems Jason Frank Abstract The paper presents symplectic Möbius integrators or Riccati equations arising in linear quadratic optimal control problems.
More informationFeedback Linearization
Feedback Linearization Peter Al Hokayem and Eduardo Gallestey May 14, 2015 1 Introduction Consider a class o single-input-single-output (SISO) nonlinear systems o the orm ẋ = (x) + g(x)u (1) y = h(x) (2)
More information