The local structure of affine systems

Size: px
Start display at page:

Download "The local structure of affine systems"

Transcription

1 The local structure of affine systems César Aguilar and Andrew D. Lewis Department of Mathematics and Statistics, Queen s University 05/06/2008 C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

2 Algebra in geometry D TM a distribution, x 0 M. π Dx0 : T x0 M T x0 M/D x0 the canonical projection. Define B D (x 0 ) 2 (D x 0 ) T x0 M/D x0 by B D (x 0 )(u, v) = π Dx0 ([U, V](x 0 )), where U and V are vector fields extending u, v D x0. Frobenius: If D is analytic or smooth with constant rank, it is integrable if and only if B D (x) = 0 for all x M. C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

3 Feedback invariance Many controllability theorems are not feedback invariant. Example On R m R n m consider a system with governing equations ẋ 1 = u, ẋ 2 = Q(x 1 ), where Q is a R n m -valued homogeneous polynomial of degree 2. Write Q(x 1 ) = (B 1 (x 1, x 1 ),..., B n m (x 1, x 1 )) for quadratic forms B 1,..., B n m. C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

4 Feedback invariance Example (cont d) By the generalised Hermes condition, (Sussmann 1 ) the system is STLC from (0, 0) if the diagonal entries in the matrices B 1,..., B n m are zero. Is this condition necessary? No. Is this condition invariant under feedback transformations of the form u Pu for P GL(m; R)? No. However... the system is STLC from (0, 0) if and only if there exists P GL(m; R) such that the diagonal entries of the matrices P T B 1 P,..., P T B n m P are zero. 1 SIAM J. Control Optim., 25(1), C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

5 Feedback invariance Control-affine systems ẋ(t) = f 0 (x(t)) + ẋ(t) = g 0 (x(t)) + m u a (t)f a (x(t)), a=1 n v b (t)g b (x(t)) on a manifold M are feedback equivalent when they have the same trajectories (essentially). Have a relation like g 0 (x) = m λ a (x)f a (x), g b (x) = a=0 b=1 m Λ a b (x)f a(x), a=1 b {1,..., n}. C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

6 Feedback invariance Since trajectories are of most interest, one s approach to control theory should be feedback invariant, i.e., identify systems that are feedback equivalent. Working with a specific f 0, f 1,..., f m is like working with coordinates in differential geometry: Sometimes it is necessary, but it is best avoided if comprehension is what is sought. Two possible approaches: 1 Choose f 0, f 1,..., f m, then show results are independent of this choice. 2 Develop an approach that does not involve this choice. C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

7 Feedback invariance Example (cont d) Consider again: ẋ 1 = u, ẋ 2 = Q(x 1 ). TFAE: 1 STLC from (0, 0); 2 there exists P GL(m; R) such that the diagonal entries of the matrices P T B 1 P,..., P T B n m P are zero; 3 0 int(conv(image(q))). C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

8 Feedback invariance Condition 2 represents the choose generators and show independence on these approach. Condition 3 represents generator independent approach. Generalisation of the first approach seems hopeless. Generalisation of the second approach... how does one do it? Let us choose ignorance over hopelessness. C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

9 Affine systems Replace the data {f 0, f 1,..., f m } with a subset A TM such that, in a neighbourhood of any point x 0 M, there exist vector fields X 0, X 1,..., X k such that A x A T x M = { X 0 (x) + k u a X a (x) u R k}. The object A is a locally finitely generated affine distribution on M. The (not uniquely defined) vector fields X 0, X 1,..., X k are local generators. a=1 C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

10 Affine systems Typically controls for a control-affine system ẋ(t) = f 0 (x(t)) + m u a (t)f a (x(t)) are restricted to lie in a subset U R m. We generalise this by... An affine system in an affine distribution A assigns to each point x M a subset A(x) A x. Require the nondegeneracy condition that aff(a(x)) = Ax and require some fussy smoothness conditions that I will not state here. a=1 A trajectory of A is a locally absolutely continuous curve γ : I M such that γ (t) A(γ(t)) for a.e. t I. An A-vector field is a map ξ : M TM such that ξ(x) A(x). C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

11 Problems for affine systems We aim to study the local properties of affine systems. To what end? 1 Understand controllability. 2 Understand stabilisability. 3 Understand relationships between controllability and stabilisability. What geometric properties of affine distributions come into play to address these problems? We sort of think Lie brackets are involved, but how? Today: Understand this through controllability. C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

12 Controllability of affine systems Have the expected notions of controllability for an affine system A. Denote R A (x 0, T) = {γ(t) γ is a trajectory on [0, T] such that γ(0) = x 0 } and R A (x 0, T) = t [0,T] R A (x 0, t). Two flavours of controllability from x 0 M: 1 Accessibility: int(r A (x 0, T)) ; 2 Small-time local controllability (STLC): x 0 int(r A (x 0, T)). C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

13 Controllability of affine systems Bang bang type theorems = it suffices to consider piecewise constant controls. Thus need to consider trajectories of the form Φ ξ 1 t 1 Φ ξp t p (x 0 ) for A-vector fields ξ 1,..., ξ p and t 1,..., t p R 0. Baker Campbell Hausdorff: For k Z >0 there exists BCH k (η 1,..., η p ) such that Φ ξ 1 t 1 Have Φ ξp t p (x 0 ) = Φ BCH k(t 1 ξ 1,...,t pξ p) 1 (x 0 ) + O((t t p ) k+1 ). BCH k (η 1,..., η p ) = η η p j,k {1,...,p} j<k [η j, η k ] + C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

14 Controllability of affine systems Standard technique in controllability: Find curves in the reachable set so the (possibly higher-order) tangent vector of this curve at x 0 is a reachable direction. Construction of such tangent vectors: 1 Abbreviate Φ ξ1 t 1 Φ ξp t p (x 0 ) = Φ ξ x 0 (t 1,..., t p ). 2 Consider the map R p 0 t Φξ x 0 (t) M. 3 Consider a curve R 0 s τ (s) R p 0 with τ (0) = 0. 4 The composition Φ ξ x 0 τ is a curve in the reachable set. 5 Suppose the first k 1 derivatives of Φ ξ x 0 τ vanish. 6 Then dk ds k s=0 Φ ξ x 0 τ (s) (assume nonzero) is a kth-order tangent vector to the reachable set. 7 Denote ord x0 (ξ, τ ) = k. C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

15 Controllability of affine systems Example Take ξ = ( Y, X, Y, X) and τ (s) = (s, s, s, s). Then ord x0 (ξ, τ ) = 2 and d 2 ds 2 Φ ξ x 0 τ (s) = [X, Y](x 0 ). s=0 This is nothing new. C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

16 Algebraic constructions using jets For analytic systems data is contained in the Taylor expansions at x 0, i.e., the infinite jets at x 0. Jet bundle notation and properties. 1 For manifolds M and N, J k (x,y)(m; N) denotes the k-jets of maps for which x y. J k (x,0)(m; R) is a R-algebra. Elements of J k (x,y)(m; N) are homomorphisms of the R-algebras J k (y,0) (N; R) and Jk (x,0)(m; R). 2 For a vector bundle π : E M, J k xπ denotes the k-jets of sections of E at x. J k xπ is a R-vector space. C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

17 Algebraic constructions using jets If ord x0 (ξ, τ ) = k we are interested in the k-jet of s Φ ξ x 0 τ (s). We have j k (Φ ξ x 0 τ )(0) = j k τ (0) j k Φ ξ x 0 (0) (composition of algebra homomorphisms). j k τ (0) J k (0,0) (R; Rp ) is a rather canonical object. What about j k Φ ξ x 0 (0)? C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

18 Algebraic constructions using jets Some notation: 1 Denote T k x M = J k (x,0) (M; R) and (Rp ) k = J k (0,0) (Rp ; R). 2 π p TM : TMp M denotes the p-fold Whitney sum. 3 For a R-vector space V, S j (V) denotes the symmetric tensors of degree j. 4 S k (V) = k j=1 Sj (V). 5 Define k : V S k (V) v v (v v) (v v). 6 L(U; V) denotes the linear maps between R-vector spaces U and V. 7 Hom(A; B) denotes the homomorphisms of R-algebras A and B. C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

19 Algebraic constructions using jets Theorem For each k, p Z >0 there exists a unique map such that T k p (x 0 ) L(S k (J k 1 x 0 π p TM ); L(T k x 0 M; (R p ) k )) T k p (x 0 )( k (j k 1 ξ(x 0 ))) = j k Φ ξ x 0 (0) for every family ξ = (ξ 1,..., ξ p ) of C -vector fields. Moreover, the diagram 1 (J 0 x 0 π p TM ) 2 (J 1 x 0 π p TM ) 3 (J 2 x 0 π p TM ) Tp 1 (x 0 ) Hom(T 1 x M; (R p ) 1 ) 0 T 2 p (x 0 ) Hom(T 2 x 0 M; (R p ) 2 ) T 3 p (x 0 ) Hom(T 3 x M; (R p ) 3 ) 0 commutes, where the horizontal arrows are the canonical projections. C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

20 Algebraic constructions using jets Important points: 1 The diagram implies proj lim k exists. 2 The map Tp k (x 0 ) can be computed (and the proof of its existence is made) using Baker Campbell Hausdorff. 3 The map Tp k (x 0 ) is system independent: It depends only on k, p, and dim(m). How is T k p (x 0 ) used for controllability? C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

21 Controllability of affine systems (again) Fix k 2. Neutralisation: Find A-vector fields ξ 1,..., ξ p and τ : R 0 R p 0 such that j k 1 τ (0) T k 1 p ( k 1 (j k 2 ξ(x 0 ))) = 0. Bad news: This is an intractable nonlinear equation for j k 1 τ (0) and j k 2 ξ(x 0 ). Good news: The equation is algebraic and it can be written down! The point: If j k 1 τ (0) and j k 2 ξ(x 0 ) satisfy the neutralisability condition, then j k τ (0) T k p ( k (j k 1 ξ(x 0 ))) is a kth-order tangent vector to the reachable set. Note: The system only comes in through its jets. C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

22 Conclusions and future work Conclusion: The local structure of an analytic affine system about x 0 is captured by the canonical algebraic objects Tp k (x 0 ) restricted to system data. Future work: 1 Must do examples = need symbolic software to do computations. 2 See what kind of controllability theorems come from this framework. 3 Think about stabilisability = connection with Control Lyapunov Functions? 4. C. Aguilar, A. Lewis (Queen s University) The local structure of affine systems 05/06/ / 22

Jacobian linearisation in a geometric setting

Jacobian linearisation in a geometric setting Jacobian linearisation in a geometric setting David R. Tyner 12/12/2003 Slide 0 Department of Mathematics and Statistics, Queen s University Email: andrew.lewis@queensu.ca URL: http://penelope.queensu.ca/~andrew/

More information

On small-time local controllability

On small-time local controllability On small-time local controllability Saber Jafarpour 1 Queen s University 7th Meeting on System and Control Theory, Kingston, 12 May 2016 1 Department of Mathematics and Statistics, Queen s University 1

More information

Hamiltonian Systems of Negative Curvature are Hyperbolic

Hamiltonian Systems of Negative Curvature are Hyperbolic Hamiltonian Systems of Negative Curvature are Hyperbolic A. A. Agrachev N. N. Chtcherbakova Abstract The curvature and the reduced curvature are basic differential invariants of the pair: Hamiltonian system,

More information

CALCULUS ON MANIFOLDS

CALCULUS ON MANIFOLDS CALCULUS ON MANIFOLDS 1. Introduction. Crash course on the Multivariate Calculus 1.1. Linear algebra. Field of real numbers R and its completeness. Real line R 1 and real spaces R n. Linear functionals,

More information

Chap. 3. Controlled Systems, Controllability

Chap. 3. Controlled Systems, Controllability Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :

More information

NOTES ON FORMAL NEIGHBORHOODS AND JET BUNDLES

NOTES ON FORMAL NEIGHBORHOODS AND JET BUNDLES NOTES ON FORMAL NEIGHBORHOODS AND JET BUNDLES SHILIN U ABSTRACT. The purpose of this note is to review the construction of smooth and holomorphic jet bundles and its relation to formal neighborhood of

More information

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,

More information

Problems in Linear Algebra and Representation Theory

Problems in Linear Algebra and Representation Theory Problems in Linear Algebra and Representation Theory (Most of these were provided by Victor Ginzburg) The problems appearing below have varying level of difficulty. They are not listed in any specific

More information

A Sufficient Condition for Local Controllability

A Sufficient Condition for Local Controllability A Sufficient Condition for Local Controllability of Nonlinear Systems Along Closed Orbits Kwanghee Nam and Aristotle Arapostathis Abstract We present a computable sufficient condition to determine local

More information

THEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)

THEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.) 4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M

More information

Definition 5.1. A vector field v on a manifold M is map M T M such that for all x M, v(x) T x M.

Definition 5.1. A vector field v on a manifold M is map M T M such that for all x M, v(x) T x M. 5 Vector fields Last updated: March 12, 2012. 5.1 Definition and general properties We first need to define what a vector field is. Definition 5.1. A vector field v on a manifold M is map M T M such that

More information

Vectors. January 13, 2013

Vectors. January 13, 2013 Vectors January 13, 2013 The simplest tensors are scalars, which are the measurable quantities of a theory, left invariant by symmetry transformations. By far the most common non-scalars are the vectors,

More information

Atiyah classes and homotopy algebras

Atiyah classes and homotopy algebras Atiyah classes and homotopy algebras Mathieu Stiénon Workshop on Lie groupoids and Lie algebroids Kolkata, December 2012 Atiyah (1957): obstruction to existence of holomorphic connections Rozansky-Witten

More information

1 Smooth manifolds and Lie groups

1 Smooth manifolds and Lie groups An undergraduate approach to Lie theory Slide 1 Andrew Baker, Glasgow Glasgow, 12/11/1999 1 Smooth manifolds and Lie groups A continuous g : V 1 V 2 with V k R m k open is called smooth if it is infinitely

More information

where m is the maximal ideal of O X,p. Note that m/m 2 is a vector space. Suppose that we are given a morphism

where m is the maximal ideal of O X,p. Note that m/m 2 is a vector space. Suppose that we are given a morphism 8. Smoothness and the Zariski tangent space We want to give an algebraic notion of the tangent space. In differential geometry, tangent vectors are equivalence classes of maps of intervals in R into the

More information

Convergence Rate of Nonlinear Switched Systems

Convergence Rate of Nonlinear Switched Systems Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the

More information

STABILITY OF PLANAR NONLINEAR SWITCHED SYSTEMS

STABILITY OF PLANAR NONLINEAR SWITCHED SYSTEMS LABORATOIRE INORMATIQUE, SINAUX ET SYSTÈMES DE SOPHIA ANTIPOLIS UMR 6070 STABILITY O PLANAR NONLINEAR SWITCHED SYSTEMS Ugo Boscain, régoire Charlot Projet TOpModel Rapport de recherche ISRN I3S/RR 2004-07

More information

MILNOR SEMINAR: DIFFERENTIAL FORMS AND CHERN CLASSES

MILNOR SEMINAR: DIFFERENTIAL FORMS AND CHERN CLASSES MILNOR SEMINAR: DIFFERENTIAL FORMS AND CHERN CLASSES NILAY KUMAR In these lectures I want to introduce the Chern-Weil approach to characteristic classes on manifolds, and in particular, the Chern classes.

More information

Mathematical Systems Theory: Advanced Course Exercise Session 5. 1 Accessibility of a nonlinear system

Mathematical Systems Theory: Advanced Course Exercise Session 5. 1 Accessibility of a nonlinear system Mathematical Systems Theory: dvanced Course Exercise Session 5 1 ccessibility of a nonlinear system Consider an affine nonlinear control system: [ ẋ = f(x)+g(x)u, x() = x, G(x) = g 1 (x) g m (x) ], where

More information

OR MSc Maths Revision Course

OR MSc Maths Revision Course OR MSc Maths Revision Course Tom Byrne School of Mathematics University of Edinburgh t.m.byrne@sms.ed.ac.uk 15 September 2017 General Information Today JCMB Lecture Theatre A, 09:30-12:30 Mathematics revision

More information

Sheaves of Lie Algebras of Vector Fields

Sheaves of Lie Algebras of Vector Fields Sheaves of Lie Algebras of Vector Fields Bas Janssens and Ori Yudilevich March 27, 2014 1 Cartan s first fundamental theorem. Second lecture on Singer and Sternberg s 1965 paper [3], by Bas Janssens. 1.1

More information

Hybrid Control and Switched Systems. Lecture #11 Stability of switched system: Arbitrary switching

Hybrid Control and Switched Systems. Lecture #11 Stability of switched system: Arbitrary switching Hybrid Control and Switched Systems Lecture #11 Stability of switched system: Arbitrary switching João P. Hespanha University of California at Santa Barbara Stability under arbitrary switching Instability

More information

2. Intersection Multiplicities

2. Intersection Multiplicities 2. Intersection Multiplicities 11 2. Intersection Multiplicities Let us start our study of curves by introducing the concept of intersection multiplicity, which will be central throughout these notes.

More information

Mathematics for Control Theory

Mathematics for Control Theory Mathematics for Control Theory Geometric Concepts in Control Involutivity and Frobenius Theorem Exact Linearization Hanz Richter Mechanical Engineering Department Cleveland State University Reading materials

More information

The existence of light-like homogeneous geodesics in homogeneous Lorentzian manifolds. Sao Paulo, 2013

The existence of light-like homogeneous geodesics in homogeneous Lorentzian manifolds. Sao Paulo, 2013 The existence of light-like homogeneous geodesics in homogeneous Lorentzian manifolds Zdeněk Dušek Sao Paulo, 2013 Motivation In a previous project, it was proved that any homogeneous affine manifold (and

More information

Generalized complex geometry and topological sigma-models

Generalized complex geometry and topological sigma-models Generalized complex geometry and topological sigma-models Anton Kapustin California Institute of Technology Generalized complex geometry and topological sigma-models p. 1/3 Outline Review of N = 2 sigma-models

More information

TIME VARYING STABILISING FEEDBACK DESIGN FOR BILINEAR SYSTEMS

TIME VARYING STABILISING FEEDBACK DESIGN FOR BILINEAR SYSTEMS IME VARYING SABILISING FEEDBACK DESIGN FOR BILINEAR SYSEMS HANNAH MICHALSKA and MIGUEL A. ORRES-ORRII Department of Electrical Engineering, McGill University 348 University Street, Montréal, H3A 2A7 Canada

More information

3.1. Derivations. Let A be a commutative k-algebra. Let M be a left A-module. A derivation of A in M is a linear map D : A M such that

3.1. Derivations. Let A be a commutative k-algebra. Let M be a left A-module. A derivation of A in M is a linear map D : A M such that ALGEBRAIC GROUPS 33 3. Lie algebras Now we introduce the Lie algebra of an algebraic group. First, we need to do some more algebraic geometry to understand the tangent space to an algebraic variety at

More information

Choice of Riemannian Metrics for Rigid Body Kinematics

Choice of Riemannian Metrics for Rigid Body Kinematics Choice of Riemannian Metrics for Rigid Body Kinematics Miloš Žefran1, Vijay Kumar 1 and Christopher Croke 2 1 General Robotics and Active Sensory Perception (GRASP) Laboratory 2 Department of Mathematics

More information

REVIEW OF DIFFERENTIAL CALCULUS

REVIEW OF DIFFERENTIAL CALCULUS REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be

More information

Affine Connections: Part 2

Affine Connections: Part 2 Affine Connections: Part 2 Manuscript for Machine Learning Reading Group Talk R. Simon Fong Abstract Note for online manuscript: This is the manuscript of a one hour introductory talk on (affine) connections.

More information

Stabilization of Nonlinear Systems in the Plane

Stabilization of Nonlinear Systems in the Plane Stabilization of Nonlinear Systems in the Plane Matthias Kawski 1 Department of Mathematics Arizona State University Tempe, Arizona 85287 Abstract. It is shown that every small-time locally controllable

More information

Chap. 1. Some Differential Geometric Tools

Chap. 1. Some Differential Geometric Tools Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U

More information

EXERCISES IN POISSON GEOMETRY

EXERCISES IN POISSON GEOMETRY EXERCISES IN POISSON GEOMETRY The suggested problems for the exercise sessions #1 and #2 are marked with an asterisk. The material from the last section will be discussed in lecture IV, but it s possible

More information

Lecture 2: Controllability of nonlinear systems

Lecture 2: Controllability of nonlinear systems DISC Systems and Control Theory of Nonlinear Systems 1 Lecture 2: Controllability of nonlinear systems Nonlinear Dynamical Control Systems, Chapter 3 See www.math.rug.nl/ arjan (under teaching) for info

More information

LECTURE 14: LIE GROUP ACTIONS

LECTURE 14: LIE GROUP ACTIONS LECTURE 14: LIE GROUP ACTIONS 1. Smooth actions Let M be a smooth manifold, Diff(M) the group of diffeomorphisms on M. Definition 1.1. An action of a Lie group G on M is a homomorphism of groups τ : G

More information

Some results on energy shaping feedback

Some results on energy shaping feedback Some results on energy shaping feedback Andrew D. Lewis Bahman Gharesifard 13/08/2007 Department of Mathematics and Statistics, Queen s University Email: andrew@mast.queensu.ca URL: http://www.mast.queensu.ca/~andrew/

More information

Quadratic forms. Here. Thus symmetric matrices are diagonalizable, and the diagonalization can be performed by means of an orthogonal matrix.

Quadratic forms. Here. Thus symmetric matrices are diagonalizable, and the diagonalization can be performed by means of an orthogonal matrix. Quadratic forms 1. Symmetric matrices An n n matrix (a ij ) n ij=1 with entries on R is called symmetric if A T, that is, if a ij = a ji for all 1 i, j n. We denote by S n (R) the set of all n n symmetric

More information

Linear Algebra and Robot Modeling

Linear Algebra and Robot Modeling Linear Algebra and Robot Modeling Nathan Ratliff Abstract Linear algebra is fundamental to robot modeling, control, and optimization. This document reviews some of the basic kinematic equations and uses

More information

Abstract. Jacobi curves are far going generalizations of the spaces of \Jacobi

Abstract. Jacobi curves are far going generalizations of the spaces of \Jacobi Principal Invariants of Jacobi Curves Andrei Agrachev 1 and Igor Zelenko 2 1 S.I.S.S.A., Via Beirut 2-4, 34013 Trieste, Italy and Steklov Mathematical Institute, ul. Gubkina 8, 117966 Moscow, Russia; email:

More information

On common linear/quadratic Lyapunov functions for switched linear systems

On common linear/quadratic Lyapunov functions for switched linear systems On common linear/quadratic Lyapunov functions for switched linear systems M. Seetharama Gowda Department of Mathematics and Statistics University of Maryland, Baltimore County Baltimore, Maryland, USA

More information

2.10 Saddles, Nodes, Foci and Centers

2.10 Saddles, Nodes, Foci and Centers 2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one

More information

A TUTORIAL ON GEOMETRIC CONTROL THEORY

A TUTORIAL ON GEOMETRIC CONTROL THEORY A TUTORIAL ON GEOMETRIC CONTROL THEORY SLOBODAN N. SIMIĆ This tutorial consists of three parts: 1. Basic concepts, 2. Basic results, 3. Steering with piecewise constant inputs. The goal is to present only

More information

Control, Stabilization and Numerics for Partial Differential Equations

Control, Stabilization and Numerics for Partial Differential Equations Paris-Sud, Orsay, December 06 Control, Stabilization and Numerics for Partial Differential Equations Enrique Zuazua Universidad Autónoma 28049 Madrid, Spain enrique.zuazua@uam.es http://www.uam.es/enrique.zuazua

More information

Invariant Lagrangian Systems on Lie Groups

Invariant Lagrangian Systems on Lie Groups Invariant Lagrangian Systems on Lie Groups Dennis Barrett Geometry and Geometric Control (GGC) Research Group Department of Mathematics (Pure and Applied) Rhodes University, Grahamstown 6140 Eastern Cape

More information

MATRIX LIE GROUPS AND LIE GROUPS

MATRIX LIE GROUPS AND LIE GROUPS MATRIX LIE GROUPS AND LIE GROUPS Steven Sy December 7, 2005 I MATRIX LIE GROUPS Definition: A matrix Lie group is a closed subgroup of Thus if is any sequence of matrices in, and for some, then either

More information

Hamiltonian flows, cotangent lifts, and momentum maps

Hamiltonian flows, cotangent lifts, and momentum maps Hamiltonian flows, cotangent lifts, and momentum maps Jordan Bell jordan.bell@gmail.com Department of Mathematics, University of Toronto April 3, 2014 1 Symplectic manifolds Let (M, ω) and (N, η) be symplectic

More information

Lecture 7. Lie brackets and integrability

Lecture 7. Lie brackets and integrability Lecture 7. Lie brackets and integrability In this lecture we will introduce the Lie bracket of two vector fields, and interpret it in several ways. 7.1 The Lie bracket. Definition 7.1.1 Let X and Y by

More information

fy (X(g)) Y (f)x(g) gy (X(f)) Y (g)x(f)) = fx(y (g)) + gx(y (f)) fy (X(g)) gy (X(f))

fy (X(g)) Y (f)x(g) gy (X(f)) Y (g)x(f)) = fx(y (g)) + gx(y (f)) fy (X(g)) gy (X(f)) 1. Basic algebra of vector fields Let V be a finite dimensional vector space over R. Recall that V = {L : V R} is defined to be the set of all linear maps to R. V is isomorphic to V, but there is no canonical

More information

Deterministic Dynamic Programming

Deterministic Dynamic Programming Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0

More information

Review: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function.

Review: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function. Taylor Series (Sect. 10.8) Review: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function. Review: Power series define functions Remarks:

More information

Lecture 8. Connections

Lecture 8. Connections Lecture 8. Connections This lecture introduces connections, which are the machinery required to allow differentiation of vector fields. 8.1 Differentiating vector fields. The idea of differentiating vector

More information

10. Smooth Varieties. 82 Andreas Gathmann

10. Smooth Varieties. 82 Andreas Gathmann 82 Andreas Gathmann 10. Smooth Varieties Let a be a point on a variety X. In the last chapter we have introduced the tangent cone C a X as a way to study X locally around a (see Construction 9.20). It

More information

Riemannian Geometry. Eckhard Meinrenken. Lecture Notes, University of Toronto, Spring 2002

Riemannian Geometry. Eckhard Meinrenken. Lecture Notes, University of Toronto, Spring 2002 Riemannian Geometry Eckhard Meinrenken Lecture Notes, University of Toronto, Spring 2002 Contents 1. Manifolds 4 2. Examples of manifolds 6 3. Submanifolds 7 4. Tangent spaces 9 5. Tangent map 11 6. Tangent

More information

The Atiyah bundle and connections on a principal bundle

The Atiyah bundle and connections on a principal bundle Proc. Indian Acad. Sci. (Math. Sci.) Vol. 120, No. 3, June 2010, pp. 299 316. Indian Academy of Sciences The Atiyah bundle and connections on a principal bundle INDRANIL BISWAS School of Mathematics, Tata

More information

QUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1)

QUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1) QUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1) Each of the six questions is worth 10 points. 1) Let H be a (real or complex) Hilbert space. We say

More information

Manifolds and tangent bundles. Vector fields and flows. 1 Differential manifolds, smooth maps, submanifolds

Manifolds and tangent bundles. Vector fields and flows. 1 Differential manifolds, smooth maps, submanifolds MA 755 Fall 05. Notes #1. I. Kogan. Manifolds and tangent bundles. Vector fields and flows. 1 Differential manifolds, smooth maps, submanifolds Definition 1 An n-dimensional C k -differentiable manifold

More information

BACKGROUND IN SYMPLECTIC GEOMETRY

BACKGROUND IN SYMPLECTIC GEOMETRY BACKGROUND IN SYMPLECTIC GEOMETRY NILAY KUMAR Today I want to introduce some of the symplectic structure underlying classical mechanics. The key idea is actually quite old and in its various formulations

More information

Modeling & Control of Hybrid Systems Chapter 4 Stability

Modeling & Control of Hybrid Systems Chapter 4 Stability Modeling & Control of Hybrid Systems Chapter 4 Stability Overview 1. Switched systems 2. Lyapunov theory for smooth and linear systems 3. Stability for any switching signal 4. Stability for given switching

More information

Reminder on basic differential geometry

Reminder on basic differential geometry Reminder on basic differential geometry for the mastermath course of 2013 Charts Manifolds will be denoted by M, N etc. One should think of a manifold as made out of points (while the elements of a vector

More information

The Toda Lattice. Chris Elliott. April 9 th, 2014

The Toda Lattice. Chris Elliott. April 9 th, 2014 The Toda Lattice Chris Elliott April 9 th, 2014 In this talk I ll introduce classical integrable systems, and explain how they can arise from the data of solutions to the classical Yang-Baxter equation.

More information

Cost-extended control systems on SO(3)

Cost-extended control systems on SO(3) Cost-extended control systems on SO(3) Ross M. Adams Mathematics Seminar April 16, 2014 R.M. Adams (RU) Cost-extended control systems on SO(3) RU Maths Seminar 1 / 34 Outline 1 Introduction 2 Control systems

More information

Stabilization and Passivity-Based Control

Stabilization and Passivity-Based Control DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive

More information

Relativistic Mechanics

Relativistic Mechanics Physics 411 Lecture 9 Relativistic Mechanics Lecture 9 Physics 411 Classical Mechanics II September 17th, 2007 We have developed some tensor language to describe familiar physics we reviewed orbital motion

More information

Varieties of Signature Tensors Lecture #1

Varieties of Signature Tensors Lecture #1 Lecture #1 Carlos Enrique Améndola Cerón (TU Munich) Summer School in Algebraic Statistics Tromsø, September 218 Setup: Signatures Let X : [, 1] R d be a piecewise differentiable path. Coordinate functions:

More information

Chapter 3. Riemannian Manifolds - I. The subject of this thesis is to extend the combinatorial curve reconstruction approach to curves

Chapter 3. Riemannian Manifolds - I. The subject of this thesis is to extend the combinatorial curve reconstruction approach to curves Chapter 3 Riemannian Manifolds - I The subject of this thesis is to extend the combinatorial curve reconstruction approach to curves embedded in Riemannian manifolds. A Riemannian manifold is an abstraction

More information

WHY POLYNOMIALS? PART 1

WHY POLYNOMIALS? PART 1 WHY POLYNOMIALS? PART 1 The proofs of finite field Kakeya and joints and short and clean, but they also seem strange to me. They seem a little like magic tricks. In particular, what is the role of polynomials

More information

LECTURE 3: SHEAVES ON R 2, CONSTRUCTIBLE WITH RESPECT TO THE STRATIFICATION BY THE AXES

LECTURE 3: SHEAVES ON R 2, CONSTRUCTIBLE WITH RESPECT TO THE STRATIFICATION BY THE AXES LECTURE 3: SHEAVES ON R 2, CONSTRUCTIBLE WITH RESPECT TO THE STRATIICATION BY THE AXES VIVEK SHENDE Today we will do a calculation in some detail of the microsupports of sheaves on R 2 which are constructible

More information

2 Two-Point Boundary Value Problems

2 Two-Point Boundary Value Problems 2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x

More information

PARALLEL TRANSPORT AND MULTI-TEMPORAL CONTROLLABILITY

PARALLEL TRANSPORT AND MULTI-TEMPORAL CONTROLLABILITY U.P.B. Sci. Bull., Series A, Vol. 77, Iss. 2, 25 ISSN 22-727 PARALLEL TRANSPORT AND MULTI-TEMPORAL CONTROLLABILITY Cristian Ghiu, Constantin Udrişte 2 Our idea was to make a bridge between the parallel

More information

CHAPTER 1 PRELIMINARIES

CHAPTER 1 PRELIMINARIES CHAPTER 1 PRELIMINARIES 1.1 Introduction The aim of this chapter is to give basic concepts, preliminary notions and some results which we shall use in the subsequent chapters of the thesis. 1.2 Differentiable

More information

Jacobian linearisation in a geometric setting

Jacobian linearisation in a geometric setting Proceedings of the 42nd IEEE Conference on Decision and Control 6084 6089, December 2003 doi: 10.1109/CDC.2003.1272230 2003 IEEE. Personal use of this material is permitted. However, permission to reprint/republish

More information

These notes are incomplete they will be updated regularly.

These notes are incomplete they will be updated regularly. These notes are incomplete they will be updated regularly. LIE GROUPS, LIE ALGEBRAS, AND REPRESENTATIONS SPRING SEMESTER 2008 RICHARD A. WENTWORTH Contents 1. Lie groups and Lie algebras 2 1.1. Definition

More information

TEST CODE: MMA (Objective type) 2015 SYLLABUS

TEST CODE: MMA (Objective type) 2015 SYLLABUS TEST CODE: MMA (Objective type) 2015 SYLLABUS Analytical Reasoning Algebra Arithmetic, geometric and harmonic progression. Continued fractions. Elementary combinatorics: Permutations and combinations,

More information

10 Transfer Matrix Models

10 Transfer Matrix Models MIT EECS 6.241 (FALL 26) LECTURE NOTES BY A. MEGRETSKI 1 Transfer Matrix Models So far, transfer matrices were introduced for finite order state space LTI models, in which case they serve as an important

More information

LOW-ORDER CONTROLLABILITY AND KINEMATIC REDUCTIONS FOR AFFINE CONNECTION CONTROL SYSTEMS

LOW-ORDER CONTROLLABILITY AND KINEMATIC REDUCTIONS FOR AFFINE CONNECTION CONTROL SYSTEMS LOW-ORDER CONTROLLABILITY AND KINEMATIC REDUCTIONS FOR AFFINE CONNECTION CONTROL SYSTEMS FRANCESCO BULLO AND ANDREW D. LEWIS Abstract. Controllability and kinematic modelling notions are investigated for

More information

Problems and Challenges in Control Theory (for Mechanical Systems)

Problems and Challenges in Control Theory (for Mechanical Systems) Problems and Challenges in Control Theory (for Mechanical Systems) 16/03/2005 Professor, Department of Mathematics and Statistics, Queen s University, Kingston, ON K7L 3N6, Canada Email: andrew.lewis@queensu.ca,

More information

class # MATH 7711, AUTUMN 2017 M-W-F 3:00 p.m., BE 128 A DAY-BY-DAY LIST OF TOPICS

class # MATH 7711, AUTUMN 2017 M-W-F 3:00 p.m., BE 128 A DAY-BY-DAY LIST OF TOPICS class # 34477 MATH 7711, AUTUMN 2017 M-W-F 3:00 p.m., BE 128 A DAY-BY-DAY LIST OF TOPICS [DG] stands for Differential Geometry at https://people.math.osu.edu/derdzinski.1/courses/851-852-notes.pdf [DFT]

More information

Dirac Structures on Banach Lie Algebroids

Dirac Structures on Banach Lie Algebroids DOI: 10.2478/auom-2014-0060 An. Şt. Univ. Ovidius Constanţa Vol. 22(3),2014, 219 228 Dirac Structures on Banach Lie Algebroids Vlad-Augustin VULCU Abstract In the original definition due to A. Weinstein

More information

Math 273, Final Exam Solutions

Math 273, Final Exam Solutions Math 273, Final Exam Solutions 1. Find the solution of the differential equation y = y +e x that satisfies the condition y(x) 0 as x +. SOLUTION: y = y H + y P where y H = ce x is a solution of the homogeneous

More information

On the role of regularity in mathematical control theory

On the role of regularity in mathematical control theory On the role of regularity in mathematical control theory by Saber Jafarpour A thesis submitted to the Department of Mathematics and Statistics in conformity with the requirements for the degree of Doctor

More information

(df (ξ ))( v ) = v F : O ξ R. with F : O ξ O

(df (ξ ))( v ) = v F : O ξ R. with F : O ξ O Math 396. Derivative maps, parametric curves, and velocity vectors Let (X, O ) and (X, O) be two C p premanifolds with corners, 1 p, and let F : X X be a C p mapping. Let ξ X be a point and let ξ = F (ξ

More information

Converse Lyapunov theorem and Input-to-State Stability

Converse Lyapunov theorem and Input-to-State Stability Converse Lyapunov theorem and Input-to-State Stability April 6, 2014 1 Converse Lyapunov theorem In the previous lecture, we have discussed few examples of nonlinear control systems and stability concepts

More information

Entrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.

Entrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4. Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the

More information

An Introduction to Differential Geometry through Computation

An Introduction to Differential Geometry through Computation Utah State University DigitalCommons@USU Textbooks Open Texts Winter 1-26-2016 An Introduction to Differential Geometry through Computation Mark E. Fels Utah State University Follow this and additional

More information

The second-order 1D wave equation

The second-order 1D wave equation C The second-order D wave equation C. Homogeneous wave equation with constant speed The simplest form of the second-order wave equation is given by: x 2 = Like the first-order wave equation, it responds

More information

CONTROLLABILITY OF NILPOTENT SYSTEMS

CONTROLLABILITY OF NILPOTENT SYSTEMS GEOMETRY IN NONLINEAR CONTROL AND DIFFERENTIAL INCLUSIONS BANACH CENTER PUBLICATIONS, VOLUME 32 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 1995 CONTROLLABILITY OF NILPOTENT SYSTEMS VICTOR

More information

12. Linear systems Theorem Let X be a scheme over a ring A. (1) If φ: X P n A is an A-morphism then L = φ O P n

12. Linear systems Theorem Let X be a scheme over a ring A. (1) If φ: X P n A is an A-morphism then L = φ O P n 12. Linear systems Theorem 12.1. Let X be a scheme over a ring A. (1) If φ: X P n A is an A-morphism then L = φ O P n A (1) is an invertible sheaf on X, which is generated by the global sections s 0, s

More information

B 1 = {B(x, r) x = (x 1, x 2 ) H, 0 < r < x 2 }. (a) Show that B = B 1 B 2 is a basis for a topology on X.

B 1 = {B(x, r) x = (x 1, x 2 ) H, 0 < r < x 2 }. (a) Show that B = B 1 B 2 is a basis for a topology on X. Math 6342/7350: Topology and Geometry Sample Preliminary Exam Questions 1. For each of the following topological spaces X i, determine whether X i and X i X i are homeomorphic. (a) X 1 = [0, 1] (b) X 2

More information

DS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.

DS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra. DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1

More information

Feedback stabilisation of locally controllable systems 1

Feedback stabilisation of locally controllable systems 1 Feedback stabilisation of locally controllable systems 1 Pantelis Isaiah 2 2012/09/24 1 PhD thesis, Department of Mathematics and Statistics, Queen s University 2 Postdoctoral Fellow, Aerospace Engineering,

More information

Linear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems

Linear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems Linear ODEs p. 1 Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Autonomous linear systems Linear ODEs Definition (Linear ODE) A linear ODE is a differential equation taking the form

More information

Three Descriptions of the Cohomology of Bun G (X) (Lecture 4)

Three Descriptions of the Cohomology of Bun G (X) (Lecture 4) Three Descriptions of the Cohomology of Bun G (X) (Lecture 4) February 5, 2014 Let k be an algebraically closed field, let X be a algebraic curve over k (always assumed to be smooth and complete), and

More information

Dimensional Reduction of Invariant Fields and Differential Operators. II. Reduction of Invariant Differential Operators

Dimensional Reduction of Invariant Fields and Differential Operators. II. Reduction of Invariant Differential Operators Dimensional Reduction of Invariant Fields and Differential Operators. II. Reduction of Invariant Differential Operators Petko A. Nikolov and Nikola P. Petrov Abstract. In the present paper, which is a

More information

On the convergence of normal forms for analytic control systems

On the convergence of normal forms for analytic control systems Chapter 1 On the convergence of normal forms for analytic control systems Wei Kang Department of Mathematics, Naval Postgraduate School Monterey, CA 93943 wkang@nps.navy.mil Arthur J. Krener Department

More information

Manipulating exponential products Instead of working with complicated concatenations of flows

Manipulating exponential products Instead of working with complicated concatenations of flows Manipulating exponential products Instead of working with complicated concatenations of flows z(t) = e t9 t2 t8 (f +f 1 +f 2 ) dt... e t1 (f +f 1 f 2 ) dt t1 e (f +f 1 +f 2 ) dt (p) it is desirable to

More information

On the Stability of the Best Reply Map for Noncooperative Differential Games

On the Stability of the Best Reply Map for Noncooperative Differential Games On the Stability of the Best Reply Map for Noncooperative Differential Games Alberto Bressan and Zipeng Wang Department of Mathematics, Penn State University, University Park, PA, 68, USA DPMMS, University

More information

Observability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u)

Observability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u) Observability Dynamic Systems Lecture 2 Observability Continuous time model: Discrete time model: ẋ(t) = f (x(t), u(t)), y(t) = h(x(t), u(t)) x(t + 1) = f (x(t), u(t)), y(t) = h(x(t)) Reglerteknik, ISY,

More information

On local normal forms of completely integrable systems

On local normal forms of completely integrable systems On local normal forms of completely integrable systems ENCUENTRO DE Răzvan M. Tudoran West University of Timişoara, România Supported by CNCS-UEFISCDI project PN-II-RU-TE-2011-3-0103 GEOMETRÍA DIFERENCIAL,

More information

An Optimal Control Problem for Rigid Body Motions in Minkowski Space

An Optimal Control Problem for Rigid Body Motions in Minkowski Space Applied Mathematical Sciences, Vol. 5, 011, no. 5, 559-569 An Optimal Control Problem for Rigid Body Motions in Minkowski Space Nemat Abazari Department of Mathematics, Ardabil Branch Islamic Azad University,

More information