From Local to Global Control

Size: px
Start display at page:

Download "From Local to Global Control"

Transcription

1 Proceeings of the 47th IEEE Conference on Decision an Control Cancun, Mexico, Dec. 9-, 8 ThB. From Local to Global Control Stephen P. Banks, M. Tomás-Roríguez. Automatic Control Engineering Department, The University of Sheffiel Sheffiel S3 JD Unite Kingom. s.banks@sheffiel.ac.uk Mathematics an Engineering School, City University. Lonon ECV HB. Unite Kingom. Maria.Tomas-Rorigue.@city.ac.uk Abstract In this paper the authors stuy the problem of the existence of multiple local operating points in control systems. In particular, they consier a metho of going from local to global control, i.e. given a number of local, linearize systems, from which global system o they come, an, can global controllers be etermine in this case? I. INTRODUCTION In many practical applications of the control of nonlinear systems, the ynamics are usually linearise about operating points such as trim conitions in aircraft. These local moels are then use to esign local controllers an some gain scheuling proceure is use to switch between the controllers at ifferent operating points. In fact, in many inustrial plants or aerospace systems, only a small number of operating conitions are known an the global nonlinear ynamical system is not known. Gain scheuling methos are wiely use an are object of many research papers, see [], [], [3] an references within as example. In this paper, the problem of going from some known local moels to a global one will be consiere. Furthermore, the obtaine global moel can then be use instea of gain scheuling for control purposes. The contains of this paper are as follows: Section II recalls the basics of nonlinear systems linearisation proceure, operating conitions an trackability properties. Section III introuces an algebraic metho that allows the reconstruction of an unknown nonlinear system taking as starting points its linear representations at some given operating conitions. Section IV presents a global control metho that can be use instea of the well-known of gain scheuling once the nonlinear system is known. Section V tackles the problem of going from the local moels to the topology of the manifol on which the system is efine. Section VI contains a summary of the ieas here presente. II. NONLINEAR SYSTEMS, OPERATING CONDITIONS AND TRACKABILITY Consier the following nonlinear system: ẋ fx,u efine on R n xr m. This is in fact, a kin of local moel for systems on manifols - this point will be aresse later. Intuitively, by an operating point it is unerstoo a pair x t,u t R consisting of an open loop control an a function x t which satisfies equation when the control u t is applie, i.e.: ẋ t fx t,u t On the other han, the esire function x t is trackable if there exists a control u t such that hols. If x t is trackable for, the local variables aroun the operating point yt an vt, are efine as: yt xt x t vt ut u t Then by applying Taylor s theorem: an ẏt ẋt ẋ t fx,u fx t,u t AtytBtvt for small yt an vt. At fx t,u t x Bt fx t,u t If At,Bt is a controllable pair, then the control /8/$5. 8 IEEE 446

2 47th IEEE CDC, Cancun, Mexico, Dec. 9-, 8 of system aroun the operating conition is achievable by using the control: ut u tvt5 In general, there coul be a number i of these operating points: x i,ui C [, ];R nm, i K 6 associate to a corresponent number i of local moels of the nonlinear system : ẏ it A ity itb itv it, 7 The problem to be consiere now is the inverse one: How to go from a set of i local moels 7 to the global one?. For simplicity, it will be assume that all the operating points are constant, so that the local moels are linear, time-invariant systems of the form: ẏ i A i y itb i v it, i K 8 Also, in this case, becomes an algebraic conition of the form: fx i,ui 9 if the values of x i are constant, i K. fx i,ui B i, 4 for i k, an moreover, f must satisfy: fx i,ui Consier the equations: with an N i N n i n, i k. 5 f p x l,ul 6 M j M m j m a p i...i n, j... j m x l i u l j l k, p n. 7 A l pq N i N q i q N n i n M j M m j m ap i...i n, j... j m i q x l i q u l j 8 B l pq N i N u i n M j M q j q M m j m ap i...i n, j... j m j q x l i u l j q i i,,i n, ThB. 9 Having obtaine a global moel from the i local ones, a global controller which will rive the system from one operating point to another has to be etermine. This will involve an application of an iteration scheme [4] which replaces a nonlinear not necessarily quaratic optimal control problem which is linear time-varying an quaratic, which can be solve by classical methos. Finally, the authors will consier the tracking global problem of etermining to some egree the topology of the manifol on which a nonlinear system is efine from a knowlege of the local representatives assuming that the local systems are complete in the sense that their efining neighborhoos cover the manifol. III. FROM LOCAL TO GLOBAL Suppose the set of k unknown constant operating points: x i,ui, i k at which there exist k linearisations of the form: ẏ it A i y itb i v it, i k of some unknown nonlinear system on R n. If the unknown system is of the form then ẋ fx,u fx i,ui A i, 3 x 447 an j j,, j m q,,, }{{},,, q in the n i N i m j M j variables a p i, j, f is assume to be able to be approximate by a polynomial function. The number of equations is an they can be written in the form: nn nmk LV W L is a linear operator, V is a vector of the unknown parameters a p i, j an W contains the known local representations A i,b i. L maps R α into R β, : an The system is: α ΠN i ΠM j β nn nmk 3 i Overetermine if α < β, ii Determine if α β an L is invertible, iii Uneretermine if α > β. Therefore, accoring to the above classification, there will be a unique solution in case ii, a solution in case i if W RangeL an in case iii if RankL β.

3 47th IEEE CDC, Cancun, Mexico, Dec. 9-, 8 ThB. IV. GLOBAL CONTROL In the previous section, a way to fin a global moel from the local moels aroun some operating conitions has been presente. Now, a global controller which rives the system from one operating conition to another is foun. Suppose the existence of two istinct operating points x,u an,u, so that in the obtaine global moel of the form: ẋ fx,u4 this will be, fx i,u i, i,. 5 Now the control objective is to rive the global nonlinear system 4 from x to, this is, to seek a esire trajectory x t such that: x x, x t f Therefore, the question is: Does there exist a control u t so that: ẋ t fx t,u t6 In orer to answer this question, the notion of projection fiel shoul be introuce: This is efine as a section of the bunle of projection operators on R n of rank m. Thus, a projection fiel on R n associates a projection operator P x : R n R m to each point x R n such that the function x P x is smooth. The main result can be summarize as: Theorem: Given a esire trajectory x : R R n, there is a control u : R R m satisfying 6, if there exist a projection fiel x P x such that the function gx,u P x fx,u satisfies: gx t,u t for each t R an ẋ t RP x t for all t R, RP is the range of the projection P. Proof: Since ẋ RP x, then: Hence, P x tẋ P x t fx t,u t gx t,u t ẋ gx t,u t an since gx t,u t the results follow from the implicit function theorem. Suppose there exists a control u t satisfying the above theorem, then from 4 an 6, it can be written that: ẋ ẋ fx,u fx,u fx x x,u u u fx,u gx x,u u,t g,,t, t. Hence, if y x x an v u u, then the system 7 can be written as: ẏ gy,v,t7 therefore, the regulator problem for v can be solve, an then write u vu. To solve this problem, a well establishe technique of linear, time-varying approximations to the problem see [4] is applie: It is assume that the system 7 can be written on the form: ẏt Ay,tytBy,tvt, 8 In the case the system is not affine in the control, nonlinear control terms can be inclue in Ay,t an By,t. Now, equation 8 is replace by the following sequence of LTV systems: ẏ [i] t Ay [i ] ty [i] tby [i ] tv [i] t9 an to each of the equations 9, the following quaratic cost functional is applie: J x[i]t t f Fx [i] t f t f [ ] x [i]t tqx [i] tu [i]t tru [i] t t 3 The problems 9 an 3 can then be solve by stanar methos see [5] or [7]. Example Let ξ,ν, ξ,ν R be two points which satisfy: ν ξ ξ ξ 3 ν ξ ξ 3, an suppose the two local systems: ẋ ξ i x These systems are local versions of: ẋ ξ an if: Then the control: x x t ξ tξ t t ẋ x x 3 u x ẋ 3 ẋ x will erive the system between the two points. u, i,. 3 u 3 448

4 47th IEEE CDC, Cancun, Mexico, Dec. 9-, 8 ThB. The system 7 becomes: ẏ y x y x x ẏ y v y 33 an this generates a close-loop control so that the following expression of ut, u vu, 34 will regulate the system between the two given operating values. V. GLOBAL SYSTEMS ON MANIFOLDS In this section, a more general problem is consiere; that of going from local moels to the topology of the ambient manifol on which the system is efine. Hence, consier a compact n-imensional ifferentiable manifol M an let Xu be a parameterize vector fiel on M, u R m. Let φ i,u i, i K be a finite covering of M by local parameters, each φ i : U i R n is a homeomorphism. In each coorinate neighborhoo, the local form of the ynamical system corresponing to the vector fiel Xu ui can be written as: ẋ i f ix i,u i. 35 It will be assume that, in each neighborhoo U i there is an operating point x i φ iu i. This, as before, means that there is a constant open loop control u i such that f ix i,ui 36 Then the new local coorinates are efine, an the new control so the local moel is of the form: y i x i x i 37 v i u i u i y 38 y i g iy i,v i 39 g iy i,v i f iy i x i,vi u i Note that g,. The question is now: Knowing the local moels 39 an the operating points x i,ui, what can be sai about the topology of M?. It is assume that complete information is known, in the sense that the coorinate patches on which the local systems are efine cover the unknown ambient manifol. So the problem is, how o they fit together?. Taking zero controls from 39: ẏ i g iy i,, i K each of these systems is efine on some region V i, say of R n. Now, it is sai that two local systems ẏ i g iy i,, ẏ j g jy j,4 are compatible if there exist a iffeomorphism φ i j : ˆV i V j from some nonempty subset ˆV i of V i onto some subset ˆV j of V j such that the system is a topologically conjugate on ˆV j, i.e. φ i j maps trajectories of systems i on those of system j. Note that, if these two systems are compatible, then y it φ i j y it an so i.e., ẏ j φ i j ẏi y i g jy i, φ i j y i giy i,. Hence the matrices φ i j y i gi from the transmition matrices for this tangent bunle of a manifol. The transition matrices are enote by γ i j φ i j y iyi Note that they obey the stanar cocycle conitions: γ i j γ jk γ ik on U i U j U k, γ ii I, γ i j γ ji I on U i U j. For example, for a sphere S, regare as S C { }: g : U U GLI,C g z z n for each integer n an U, U are neighborhoos of an respectively. This efines a complex line bunle on S, usually enote by H n. A connection on the vector bunle efine by transition functions γ i j is a collection of ifferential operators w i efine on U i such that w i γ i j γi j γ i j w j γi j on U i U j is the exterior erivative. These glue together to give a global map: for a bunle E. A : Ω E Ω E4 The curvature of A is A an is represente locally by a matrix K i of two-forms for which K i γ i j K j γ i j 4 The K th characteristic class of the bunle is efine by τ k E [τ k A] H k M;R43 449

5 47th IEEE CDC, Cancun, Mexico, Dec. 9-, 8 ThB. τa trace [ π i K j k]. The first Chern class c E τ E. Then, < c τm,[m] > π M KA XM by the Gauss-Bonnet theorem, so the Chern class states if the system is trivial or not. Example The systems: ẋ x u 44 ẏ ẏ y y u 45 Generate a system in S, since the transition functions look like g z z in local coorinates. It can be seen that < c,τm,[m] > XM VI. CONCLUSIONS In this paper the authors have consiere the problem of piecing together a set of given local systems to form a global one on some manifol. If the information is incomplete as in most practical cases interpolation methos can be use an if it is complete, then topological manifol theory will be use to escribe the ambient manifol. Using an approximation metho for the obtaine nonlinear system, global controllers can be now esigne in orer to rive the system from one local operating point to another. REFERENCES [] Hunt, K.J., Johasen,T.A. Design an analysis of gain-scheule control using local controller networks. International Journal of Control, 66, [] Hye, R.A., Glover, K.A. Comparison of Different Scheuling Techniques for H Controllers. Proceeings of the Institute of Measurement an Control Symposium on Robust Control, Cambrige. 99. [3] Hye, R.A., Glover, K.A. The Application of Scheule H Controllers to a VSTOL Aircraft. IEEE Transactions on Automatic Control, 38, [4] M. Tomas-Roriguez an S.P. Banks, Linear Approximations to Nonlinear Dynamical Systems with Applications to Stability an Spectral Theory. IMA J. Math. Cont Inf.,, 3, pp [5] M. Tomas-Roriguez, M., Navarro Hernanez, C. an Banks, S.P., Parametric Approach to Optimal Nonlinear Control Problem using Orthogonal Expansions, In proceeings of 6th IFAC Worl Congress, Prague, July 5, Cz. [6] Navarro Hernanez, C, Banks, S P an Aleen. Observer Design for Nonlinear Systems using Linear Approximations. IMA J. Math. Cont Inf, [7] Tayfun Cimen an Stephen P. Banks. Nonlinear optimal tracking control with application to super-tankers for autopilot esign. Automatica, Volume 4, Issue, November 4, Pages

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the

More information

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements

Adaptive Gain-Scheduled H Control of Linear Parameter-Varying Systems with Time-Delayed Elements Aaptive Gain-Scheule H Control of Linear Parameter-Varying Systems with ime-delaye Elements Yoshihiko Miyasato he Institute of Statistical Mathematics 4-6-7 Minami-Azabu, Minato-ku, okyo 6-8569, Japan

More information

DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10

DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10 DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10 5. Levi-Civita connection From now on we are intereste in connections on the tangent bunle T X of a Riemanninam manifol (X, g). Out main result will be a construction

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

ELEC3114 Control Systems 1

ELEC3114 Control Systems 1 ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.

More information

Approximate reduction of dynamic systems

Approximate reduction of dynamic systems Systems & Control Letters 57 2008 538 545 www.elsevier.com/locate/sysconle Approximate reuction of ynamic systems Paulo Tabuaa a,, Aaron D. Ames b, Agung Julius c, George J. Pappas c a Department of Electrical

More information

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain Nonlinear Aaptive Ship Course Tracking Control Base on Backstepping an Nussbaum Gain Jialu Du, Chen Guo Abstract A nonlinear aaptive controller combining aaptive Backstepping algorithm with Nussbaum gain

More information

Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity

Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Preprints of the 9th Worl Congress The International Feeration of Automatic Control Cape Town, South Africa. August -9, Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Zhengqiang

More information

Slide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13)

Slide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13) Slie10 Haykin Chapter 14: Neuroynamics (3r E. Chapter 13) CPSC 636-600 Instructor: Yoonsuck Choe Spring 2012 Neural Networks with Temporal Behavior Inclusion of feeback gives temporal characteristics to

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching Time Optimization in Discretized Hybrid Dynamical Systems Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set

More information

Multi-agent Systems Reaching Optimal Consensus with Time-varying Communication Graphs

Multi-agent Systems Reaching Optimal Consensus with Time-varying Communication Graphs Preprints of the 8th IFAC Worl Congress Multi-agent Systems Reaching Optimal Consensus with Time-varying Communication Graphs Guoong Shi ACCESS Linnaeus Centre, School of Electrical Engineering, Royal

More information

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy, NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

Interpolated Rigid-Body Motions and Robotics

Interpolated Rigid-Body Motions and Robotics Interpolate Rigi-Boy Motions an Robotics J.M. Selig Faculty of Business, Computing an Info. Management. Lonon South Bank University, Lonon SE AA, U.K. seligjm@lsbu.ac.uk Yaunquing Wu Dept. Mechanical Engineering.

More information

Tensors, Fields Pt. 1 and the Lie Bracket Pt. 1

Tensors, Fields Pt. 1 and the Lie Bracket Pt. 1 Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 PHYS 500 - Southern Illinois University September 8, 2016 PHYS 500 - Southern Illinois University Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 September 8,

More information

On state representations of time-varying nonlinear systems

On state representations of time-varying nonlinear systems On state representations of time-varying nonlinear systems Paulo Sérgio Pereira a Silva a, Simone Batista a, a University of São Paulo, Escola Politécnicca PTC Av. Luciano Gualberto trav. 03, 158, 05508-900

More information

Darboux s theorem and symplectic geometry

Darboux s theorem and symplectic geometry Darboux s theorem an symplectic geometry Liang, Feng May 9, 2014 Abstract Symplectic geometry is a very important branch of ifferential geometry, it is a special case of poisson geometry, an coul also

More information

A global Implicit Function Theorem without initial point and its applications to control of non-affine systems of high dimensions

A global Implicit Function Theorem without initial point and its applications to control of non-affine systems of high dimensions J. Math. Anal. Appl. 313 (2006) 251 261 www.elsevier.com/locate/jmaa A global Implicit Function Theorem without initial point an its applications to control of non-affine systems of high imensions Weinian

More information

arxiv: v1 [math.dg] 1 Nov 2015

arxiv: v1 [math.dg] 1 Nov 2015 DARBOUX-WEINSTEIN THEOREM FOR LOCALLY CONFORMALLY SYMPLECTIC MANIFOLDS arxiv:1511.00227v1 [math.dg] 1 Nov 2015 ALEXANDRA OTIMAN AND MIRON STANCIU Abstract. A locally conformally symplectic (LCS) form is

More information

Rank, Trace, Determinant, Transpose an Inverse of a Matrix Let A be an n n square matrix: A = a11 a1 a1n a1 a an a n1 a n a nn nn where is the jth col

Rank, Trace, Determinant, Transpose an Inverse of a Matrix Let A be an n n square matrix: A = a11 a1 a1n a1 a an a n1 a n a nn nn where is the jth col Review of Linear Algebra { E18 Hanout Vectors an Their Inner Proucts Let X an Y be two vectors: an Their inner prouct is ene as X =[x1; ;x n ] T Y =[y1; ;y n ] T (X; Y ) = X T Y = x k y k k=1 where T an

More information

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS Francesco Bullo Richar M. Murray Control an Dynamical Systems California Institute of Technology Pasaena, CA 91125 Fax : + 1-818-796-8914 email

More information

Linear Algebra- Review And Beyond. Lecture 3

Linear Algebra- Review And Beyond. Lecture 3 Linear Algebra- Review An Beyon Lecture 3 This lecture gives a wie range of materials relate to matrix. Matrix is the core of linear algebra, an it s useful in many other fiels. 1 Matrix Matrix is the

More information

Linear and quadratic approximation

Linear and quadratic approximation Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

arxiv: v2 [math.dg] 16 Dec 2014

arxiv: v2 [math.dg] 16 Dec 2014 A ONOTONICITY FORULA AND TYPE-II SINGULARITIES FOR THE EAN CURVATURE FLOW arxiv:1312.4775v2 [math.dg] 16 Dec 2014 YONGBING ZHANG Abstract. In this paper, we introuce a monotonicity formula for the mean

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

The Non-abelian Hodge Correspondence for Non-Compact Curves

The Non-abelian Hodge Correspondence for Non-Compact Curves 1 Section 1 Setup The Non-abelian Hoge Corresponence for Non-Compact Curves Chris Elliott May 8, 2011 1 Setup In this talk I will escribe the non-abelian Hoge theory of a non-compact curve. This was worke

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

Approximate Reduction of Dynamical Systems

Approximate Reduction of Dynamical Systems Proceeings of the 4th IEEE Conference on Decision & Control Manchester Gran Hyatt Hotel San Diego, CA, USA, December 3-, 6 FrIP.7 Approximate Reuction of Dynamical Systems Paulo Tabuaa, Aaron D. Ames,

More information

u t v t v t c a u t b a v t u t v t b a

u t v t v t c a u t b a v t u t v t b a Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying

More information

Math 11 Fall 2016 Section 1 Monday, September 19, Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v

Math 11 Fall 2016 Section 1 Monday, September 19, Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v Math Fall 06 Section Monay, September 9, 06 First, some important points from the last class: Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v passing through

More information

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Permanent vs. Determinant

Permanent vs. Determinant Permanent vs. Determinant Frank Ban Introuction A major problem in theoretical computer science is the Permanent vs. Determinant problem. It asks: given an n by n matrix of ineterminates A = (a i,j ) an

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

More information

Laplacian Cooperative Attitude Control of Multiple Rigid Bodies

Laplacian Cooperative Attitude Control of Multiple Rigid Bodies Laplacian Cooperative Attitue Control of Multiple Rigi Boies Dimos V. Dimarogonas, Panagiotis Tsiotras an Kostas J. Kyriakopoulos Abstract Motivate by the fact that linear controllers can stabilize the

More information

State observers and recursive filters in classical feedback control theory

State observers and recursive filters in classical feedback control theory State observers an recursive filters in classical feeback control theory State-feeback control example: secon-orer system Consier the riven secon-orer system q q q u x q x q x x x x Here u coul represent

More information

On classical orthogonal polynomials and differential operators

On classical orthogonal polynomials and differential operators INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS A: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 38 2005) 6379 6383 oi:10.1088/0305-4470/38/28/010 On classical orthogonal polynomials an ifferential

More information

Left-invariant extended Kalman filter and attitude estimation

Left-invariant extended Kalman filter and attitude estimation Left-invariant extene Kalman filter an attitue estimation Silvere Bonnabel Abstract We consier a left-invariant ynamics on a Lie group. One way to efine riving an observation noises is to make them preserve

More information

1 Heisenberg Representation

1 Heisenberg Representation 1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.

More information

A new approach to explicit MPC using self-optimizing control

A new approach to explicit MPC using self-optimizing control 28 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June -3, 28 WeA3.2 A new approach to explicit MPC using self-optimizing control Henrik Manum, Sriharakumar Narasimhan an Sigur

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

ON THE RIEMANN EXTENSION OF THE SCHWARZSCHILD METRICS

ON THE RIEMANN EXTENSION OF THE SCHWARZSCHILD METRICS ON THE RIEANN EXTENSION OF THE SCHWARZSCHILD ETRICS Valerii Dryuma arxiv:gr-qc/040415v1 30 Apr 004 Institute of athematics an Informatics, AS R, 5 Acaemiei Street, 08 Chisinau, olova, e-mail: valery@ryuma.com;

More information

The canonical controllers and regular interconnection

The canonical controllers and regular interconnection Systems & Control Letters ( www.elsevier.com/locate/sysconle The canonical controllers an regular interconnection A.A. Julius a,, J.C. Willems b, M.N. Belur c, H.L. Trentelman a Department of Applie Mathematics,

More information

Minimum-time constrained velocity planning

Minimum-time constrained velocity planning 7th Meiterranean Conference on Control & Automation Makeonia Palace, Thessaloniki, Greece June 4-6, 9 Minimum-time constraine velocity planning Gabriele Lini, Luca Consolini, Aurelio Piazzi Università

More information

Implicit Lyapunov control of closed quantum systems

Implicit Lyapunov control of closed quantum systems Joint 48th IEEE Conference on Decision an Control an 28th Chinese Control Conference Shanghai, P.R. China, December 16-18, 29 ThAIn1.4 Implicit Lyapunov control of close quantum systems Shouwei Zhao, Hai

More information

APPPHYS 217 Thursday 8 April 2010

APPPHYS 217 Thursday 8 April 2010 APPPHYS 7 Thursay 8 April A&M example 6: The ouble integrator Consier the motion of a point particle in D with the applie force as a control input This is simply Newton s equation F ma with F u : t q q

More information

Curvature, Conformal Mapping, and 2D Stationary Fluid Flows. Michael Taylor

Curvature, Conformal Mapping, and 2D Stationary Fluid Flows. Michael Taylor Curvature, Conformal Mapping, an 2D Stationary Flui Flows Michael Taylor 1. Introuction Let Ω be a 2D Riemannian manifol possibly with bounary). Assume Ω is oriente, with J enoting counterclockwise rotation

More information

model considered before, but the prey obey logistic growth in the absence of predators. In

model considered before, but the prey obey logistic growth in the absence of predators. In 5.2. First Orer Systems of Differential Equations. Phase Portraits an Linearity. Section Objective(s): Moifie Preator-Prey Moel. Graphical Representations of Solutions. Phase Portraits. Vector Fiels an

More information

BLOW-UP FORMULAS FOR ( 2)-SPHERES

BLOW-UP FORMULAS FOR ( 2)-SPHERES BLOW-UP FORMULAS FOR 2)-SPHERES ROGIER BRUSSEE In this note we give a universal formula for the evaluation of the Donalson polynomials on 2)-spheres, i.e. smooth spheres of selfintersection 2. Note that

More information

Some Remarks on the Boundedness and Convergence Properties of Smooth Sliding Mode Controllers

Some Remarks on the Boundedness and Convergence Properties of Smooth Sliding Mode Controllers International Journal of Automation an Computing 6(2, May 2009, 154-158 DOI: 10.1007/s11633-009-0154-z Some Remarks on the Bouneness an Convergence Properties of Smooth Sliing Moe Controllers Wallace Moreira

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

Characteristic classes of vector bundles

Characteristic classes of vector bundles Characteristic classes of vector bunles Yoshinori Hashimoto 1 Introuction Let be a smooth, connecte, compact manifol of imension n without bounary, an p : E be a real or complex vector bunle of rank k

More information

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables*

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables* 51st IEEE Conference on Decision an Control December 1-13 212. Maui Hawaii USA Total Energy Shaping of a Class of Uneractuate Port-Hamiltonian Systems using a New Set of Close-Loop Potential Shape Variables*

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

COUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS

COUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS VI International Conference on Computational Methos for Couple Problems in Science an Engineering COUPLED PROBLEMS 15 B. Schrefler, E. Oñate an M. Paparakakis(Es) COUPLING REQUIREMENTS FOR WELL POSED AND

More information

PDE Notes, Lecture #11

PDE Notes, Lecture #11 PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =

More information

Multirate Feedforward Control with State Trajectory Generation based on Time Axis Reversal for Plant with Continuous Time Unstable Zeros

Multirate Feedforward Control with State Trajectory Generation based on Time Axis Reversal for Plant with Continuous Time Unstable Zeros Multirate Feeforwar Control with State Trajectory Generation base on Time Axis Reversal for with Continuous Time Unstable Zeros Wataru Ohnishi, Hiroshi Fujimoto Abstract with unstable zeros is known as

More information

LECTURE 1: BASIC CONCEPTS, PROBLEMS, AND EXAMPLES

LECTURE 1: BASIC CONCEPTS, PROBLEMS, AND EXAMPLES LECTURE 1: BASIC CONCEPTS, PROBLEMS, AND EXAMPLES WEIMIN CHEN, UMASS, SPRING 07 In this lecture we give a general introuction to the basic concepts an some of the funamental problems in symplectic geometry/topology,

More information

Polynomial Inclusion Functions

Polynomial Inclusion Functions Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl

More information

Energy-preserving affine connections

Energy-preserving affine connections 2 A. D. Lewis Enery-preservin affine connections Anrew D. Lewis 28/07/1997 Abstract A Riemannian affine connection on a Riemannian manifol has the property that is preserves the kinetic enery associate

More information

Math 300 Winter 2011 Advanced Boundary Value Problems I. Bessel s Equation and Bessel Functions

Math 300 Winter 2011 Advanced Boundary Value Problems I. Bessel s Equation and Bessel Functions Math 3 Winter 2 Avance Bounary Value Problems I Bessel s Equation an Bessel Functions Department of Mathematical an Statistical Sciences University of Alberta Bessel s Equation an Bessel Functions We use

More information

Invariant Extended Kalman Filter: Theory and application to a velocity-aided estimation problem

Invariant Extended Kalman Filter: Theory and application to a velocity-aided estimation problem Invariant Extene Kalman Filter: Theory an application to a velocity-aie estimation problem S. Bonnabel (Mines ParisTech) Joint work with P. Martin (Mines ParisTech) E. Salaun (Georgia Institute of Technology)

More information

An algebraic expression of stable inversion for nonminimum phase systems and its applications

An algebraic expression of stable inversion for nonminimum phase systems and its applications Proceeings of the 17th Worl Congress The International Feeration of Automatic Control An algebraic expression of stable inversion for nonminimum phase systems an its applications Takuya Sogo Chubu University,

More information

Practical Computation of Flat Outputs for Nonlinear Control Systems

Practical Computation of Flat Outputs for Nonlinear Control Systems IOP onference Series: Materials Science an Engineering PPER OPEN ESS Practical omputation of Flat Outputs for Nonlinear ontrol Systems To cite this article: Joe Imae et al 25 IOP onf Ser: Mater Sci Eng

More information

WUCHEN LI AND STANLEY OSHER

WUCHEN LI AND STANLEY OSHER CONSTRAINED DYNAMICAL OPTIMAL TRANSPORT AND ITS LAGRANGIAN FORMULATION WUCHEN LI AND STANLEY OSHER Abstract. We propose ynamical optimal transport (OT) problems constraine in a parameterize probability

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

Summary: Differentiation

Summary: Differentiation Techniques of Differentiation. Inverse Trigonometric functions The basic formulas (available in MF5 are: Summary: Differentiation ( sin ( cos The basic formula can be generalize as follows: Note: ( sin

More information

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes Differentiation of vectors 12.5 Introuction The area known as vector calculus is use to moel mathematically a vast range of engineering phenomena incluing electrostatics, electromagnetic fiels, air flow

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

THE DIVERGENCE. 1. Interpretation of the divergence 1 2. The divergence theorem Gauss law 6

THE DIVERGENCE. 1. Interpretation of the divergence 1 2. The divergence theorem Gauss law 6 THE DIVERGENCE CRISTIAN E. GUTIÉRREZ SEPTEMBER 16, 2008 CONTENTS 1. Interpretation of the ivergence 1 2. The ivergence theorem 4 2.1. Gauss law 6 1. INTERPRETATION OF THE DIVERGENCE Let F(, t) : R n R

More information

SYMPLECTIC GEOMETRY: LECTURE 3

SYMPLECTIC GEOMETRY: LECTURE 3 SYMPLECTIC GEOMETRY: LECTURE 3 LIAT KESSLER 1. Local forms Vector fiels an the Lie erivative. A vector fiel on a manifol M is a smooth assignment of a vector tangent to M at each point. We think of M as

More information

Convergence rates of moment-sum-of-squares hierarchies for optimal control problems

Convergence rates of moment-sum-of-squares hierarchies for optimal control problems Convergence rates of moment-sum-of-squares hierarchies for optimal control problems Milan Kora 1, Diier Henrion 2,3,4, Colin N. Jones 1 Draft of September 8, 2016 Abstract We stuy the convergence rate

More information

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an

More information

Optimal Control of Spatially Distributed Systems

Optimal Control of Spatially Distributed Systems Optimal Control of Spatially Distribute Systems Naer Motee an Ali Jababaie Abstract In this paper, we stuy the structural properties of optimal control of spatially istribute systems. Such systems consist

More information

II. First variation of functionals

II. First variation of functionals II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent

More information

Applications of the Wronskian to ordinary linear differential equations

Applications of the Wronskian to ordinary linear differential equations Physics 116C Fall 2011 Applications of the Wronskian to orinary linear ifferential equations Consier a of n continuous functions y i (x) [i = 1,2,3,...,n], each of which is ifferentiable at least n times.

More information

Exam 2 Review Solutions

Exam 2 Review Solutions Exam Review Solutions 1. True or False, an explain: (a) There exists a function f with continuous secon partial erivatives such that f x (x, y) = x + y f y = x y False. If the function has continuous secon

More information

Optimization of a point-mass walking model using direct collocation and sequential quadratic programming

Optimization of a point-mass walking model using direct collocation and sequential quadratic programming Optimization of a point-mass walking moel using irect collocation an sequential quaratic programming Chris Dembia June 5, 5 Telescoping actuator y Stance leg Point-mass boy m (x,y) Swing leg x Leg uring

More information

Robustness of the Moore-Greitzer Compressor Model s Surge Subsystem with New Dynamic Output Feedback Controllers

Robustness of the Moore-Greitzer Compressor Model s Surge Subsystem with New Dynamic Output Feedback Controllers Preprints of the 19th Worl Congress The International Feeration of Automatic Control Cape Town, South Africa. August 4-9, 14 Robustness of the Moore-Greiter Compressor Moel s Surge Subsystem with New Dynamic

More information

Characterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations

Characterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations Characterizing Real-Value Multivariate Complex Polynomials an Their Symmetric Tensor Representations Bo JIANG Zhening LI Shuzhong ZHANG December 31, 2014 Abstract In this paper we stuy multivariate polynomial

More information

u!i = a T u = 0. Then S satisfies

u!i = a T u = 0. Then S satisfies Deterministic Conitions for Subspace Ientifiability from Incomplete Sampling Daniel L Pimentel-Alarcón, Nigel Boston, Robert D Nowak University of Wisconsin-Maison Abstract Consier an r-imensional subspace

More information

Relation between the propagator matrix of geodesic deviation and the second-order derivatives of the characteristic function

Relation between the propagator matrix of geodesic deviation and the second-order derivatives of the characteristic function Journal of Electromagnetic Waves an Applications 203 Vol. 27 No. 3 589 60 http://x.oi.org/0.080/0920507.203.808595 Relation between the propagator matrix of geoesic eviation an the secon-orer erivatives

More information

Energy behaviour of the Boris method for charged-particle dynamics

Energy behaviour of the Boris method for charged-particle dynamics Version of 25 April 218 Energy behaviour of the Boris metho for charge-particle ynamics Ernst Hairer 1, Christian Lubich 2 Abstract The Boris algorithm is a wiely use numerical integrator for the motion

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the

More information

Function Spaces. 1 Hilbert Spaces

Function Spaces. 1 Hilbert Spaces Function Spaces A function space is a set of functions F that has some structure. Often a nonparametric regression function or classifier is chosen to lie in some function space, where the assume structure

More information

A check digit system over a group of arbitrary order

A check digit system over a group of arbitrary order 2013 8th International Conference on Communications an Networking in China (CHINACOM) A check igit system over a group of arbitrary orer Yanling Chen Chair of Communication Systems Ruhr University Bochum

More information

LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form

LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1 We consier ifferential equations of the form y + a()y = b(), (1) y( 0 ) = y 0, where a() an b() are functions. Observe that this class of equations inclues equations

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Topic 7: Convergence of Random Variables

Topic 7: Convergence of Random Variables Topic 7: Convergence of Ranom Variables Course 003, 2016 Page 0 The Inference Problem So far, our starting point has been a given probability space (S, F, P). We now look at how to generate information

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information