Index Theory and Periodic Solution of Delay Differential Systems
|
|
- Rolf Shepherd
- 5 years ago
- Views:
Transcription
1 Index Theory and Periodic Solution of Delay Differential Systems Liu Chungen School of Mathematics, Nankai University workshop on nonlinear PDE and calculus of variation, Chern Mathematical Institute
2 Outline 1 M-boundary problem of a Hamiltonian system 2 Background of delay differential equations First order delay Hamiltonian systems 6 Second order delay Hamiltonian systems
3 1 M-boundary problem of a Hamiltonian system For a skew-symmetric non-degenerate 2N 2N matrix J = (a ij ), it can define a symplectic structure on R 2N by ω(v, w) = v T J 1 w or ω = 1 2 a ij dx i dx j with J 1 = (a ij ). A 2N 2N matrix M is called i,j J -symplectic if there holds M T J 1 M = J 1. We denote the set of all J -symplectic matrices by Sp ( J (2N). The ) usual symplectic group Sp(2N) is the 0 IN special case of J = J N =, i.e., Sp(2N) = Sp I N 0 JN (2N). Here I N is the N N identity matrix. We will write J for J N if the dimension 2N is clear from the text
4 For a J -symplectic matrix M with M k = I 2N, and a function H C 2 (R R 2n ) with H(t + τ, Mz) = H(t, z), we consider kτ-periodic solution of the following M-boundary value problem { ż(t) = J H(t, z(t)) (1) z(τ) = Mz(0). The corresponding functional is defined in E = W 1/2,2 (S 1, R 2n ) with S 1 = R/(kτZ) by ϕ(z) = 1 2 kτ 0 (J 1 ż(t), z(t))dt kτ 0 H(t, z(t))dt. (2) The critical point of ϕ in E is a kτ-periodic solution of the nonlinear Hamiltonian system in (1). In order to solve the problem (1), we define a group action σ on E by σz(t) = Mz(t τ)
5 It is clear that σ k = id and ϕ is σ-invariant, i.e., there holds ϕ(σz) = ϕ(z). (3) Setting E σ = {z E σz = z} = fix(σ), by the well known Palais symmetric principal(see [27]), a critical point of ϕ in E σ is a solution of the boundary problem (1). For z E σ, there holds ϕ(z) = k 2 τ 0 (J 1 ż(t), z(t))dt k τ 0 H(t, z(t))dt
6 The linearized system along a solution z(t) of the nonlinear Hamiltonian system in (1) is the following linear Hamiltonian system ẏ(t) = J H (t, z(t))y(t). Its fundamental solution γ z (t) with γ z (0) = I 2N should satisfy γ z (t) = J H (t, z(t))γ z (t). The following result is well known. Lemma 1.1 γ z is an J -symplectic path, i.e., γ z (t) T J 1 γ z (t) = J 1 for all t R. For simplicity we take τ = 1. For a symmetric continuous matrix function B(t) satisfying M T B(t + 1)M = B(t), suppose γ B (t) is the fundamental solution of the linear Hamiltonian system ẏ(t) = J B(t)y(t)
7 Definition 1.1 The (J, M)-nullity of the symmetric matrix function B is defined by ν J M (B) = dim C ker C (γ B (1) M). For the standard case of J = J N and a matrix P Sp JN (2N), the (J, P )-nullity and Maslov-type index of a symmetric matrix function B was defined first in 2006 by algebra method. We will define the (J, M)-index i J M (B) below via analytic method. So the index pair (i J M (B), νj M (B)) Z {0, 1,, 2N} makes sense for all symmetric continuous matrix function B(t) satisfying M T B(t + 1)M = B(t)
8 Definition 1.2 We define the relative index by I(A M, A M B M ) = m d (P n(a M B M )P n ) m d (P na M P n ), n n, where n > 0 is a constant large enough such that the difference becomes a constant independent of n n. In this case, we define i J M ( B) = I(A M, A M B M ), ν J M ( B) = m 0 (A M B M ). Here we denote by m d (A) the total multiplicity of eigenvalues of the adjointoperator A in the interval (, d], and m 0 (A) = dim ker(a). We recall that P n : E σ En σ is the projection maps and {En} σ is a Galerkin approximation sequence. A M and B M are operators in E σ defined respectively by (A M z, z) = 1 0 (J 1 ż(t), z(t))dt, (B M z, z) = 1 0 (B(t)z(t), z(t))dt. Up to a sign, the relative index I(A M, A M B M ) is exactly the spectral flow of the operator family A M sb M. Precisely there holds I(A M, A M B M ) = sf(a M sb M ). Here sf(a) is the spectral flow of the bounded self-adjoint operator A
9 2 Background of delay differential equations In 1974, Kaplan and Yorke in [Kaplan-Yorke 1974] studied the autonomous delay differential equation and introduced a new technique for establishing the existence of periodic solutions. More precisely, Kaplan and Yorke solved the periodic solutions of the following kinds of delay differential equations and with odd function f. ẋ(t) = f(x(t 1)) ẋ(t) = f(x(t 1)) + f(x(t 2))
10 They turned their problems into the problems of periodic solution of autonomous Hamiltonian system and under some twisted condition on the origin and infinity for the function f, it was proved that there exists at least one periodic solution. Since then many papers (see [G.Fei 1995,2006], [Z.Guo-J.Yu 2005, 2011], [J.Li-X.He 1998,1999] and the references therein) used Kaplan and Y- orke s original idea to search for periodic solutions of more general differential delay equations of the following form ẋ(t) = f(x(t 1)) + f(x(t 2)) + + f(x(t m + 1)). The existence of periodic solutions of above delay differential equation has been investigated by Nussbaum in [Nussbaum 1978] using different techniques
11 3 For simplicity, we first consider the 4τ-periodic solutions of the following delay differential systems x (t) = V (t, x(t τ)), (4) where the function V C 2 (R R n, R) is τ-periodic in variable t and is even in variables x. To find 4τ-periodic solution x(t), we only need to find solution with x(t + 2τ) = x(t). If x(t) is so a solution, let x 1 (t) = x(t), x 2 (t) = x(t τ) and z(t) = (x 1 (t), x 2 (t)) T, then there holds { x 1 (t) = V (t, x 2 (t)) x (5) 2(t) = V (t, x 1 (t))
12 Set H(t, x 1, x 2 ) = V (t, x 1 ) + V (t, x 2 ), then we can rewrite (5) as ( ) 0 In ż(t) = J 1 H(t, z(t)), J 1 =. (6) I n 0 Moreover, if z(t) = (x 1 (t), x 2 (t)) T is a 4τ-periodic solution of (6) with z(t) = σz(t) for the 4-periodic action σz(t) = J 1 z(t + τ), (7) then x(t) = x 1 (t) is a solution of (4) with x(t + 2τ) = x(t). The condition (7) is equivalent to z(τ) = J1 1 z(0). (8) So the problem can be transformed to the problem (1) with J = J 1 and M = J
13 In general, for a function V C 2 (R R n, R) with period τ in variable t and even in variables x, we consider the 2mτ-periodic solutions of the following delay differential system x (t) = V (t, x(t τ))+ V (t, x(t 2τ))+ + V (t, x(t (m 1)τ)). (9)
14 If we get a solution x(t) with x(t mτ) = x(t), then by setting x 1 (t) = x(t), x 2 (t) = x(t τ),, x m (t) = x(t (m 1)τ) and H(t, x 1,, x m ) = V (t, x 1 ) + + V (t, x m ), z = (x 1,, x m ) T, we rewrite the system (9) as ż(t) = A m H(t, z(t)), (10) where the mn mn skew symmetric matrix A m is defined by 0 I n I n A m = I n 0 I n.... I n I n 0 We see that det A m 0 if m 2N and det A m = 0 if m 2N
15 Setting 0 I n 0 0 ( T m = 0 0 I n = I n 0 I n In(m 1) ). If z(t) is a solution of (10) with z(τ) = Tm 1 z(0) then x(t) = x 1 (t) is a 2mτperiodic solution of equation (9) for m 2N. So we also turn our problem into the problem (1) with J = A m and M = Tm 1. In order to understand the case of m 2N + 1, we first recall the notation about
16 4 Any k k skew symmetric matrix A determine a on R k. For any functions F, H C (R k, R), the defined by {F, H} = ( F ) T A H C (R k, R). (11) We recall the definition of a {, } on a manifold M. For any two functions F, H C (M), the Poisson bracket {F, H} C (M) is defined by the following properties (1) {c 1 F 1 + c 2 F 2, H} = c 1 {F 1, H} + c 2 {F 2, H}, c 1, c 2 R, (2) {F, H} = {H, F }, (3) {{F, H}, P } + {{P, F }, H} + {{H, P }, F } = 0, (4) {F, H P } = {F, H}P + H{F, P }
17 A differential manifold M with a is called Poisson manifold. Furthermore if there holds (5) {F, H} = 0 for any function F implies H c, then the {, } determines a symplectic structure, and M is a symplectic manifold. For example, when k = 2m and A = J m = then it defines a on R 2m by m ( F H {F, H} = F ) H. q i p i p i q i i=1 ( 0 Im I m 0 This on R 2m is the standard and it determines the standard symplectic structure on R 2m. ),
18 Let M be a Poisson manifold and H : M R be a smooth function. The Hamiltonian vector field associated with H is the unique smooth vector field V H on M satisfying V H (F ) = {F, H} for every smooth function F : M R. If the on R k is defined by (11), the Hamiltonian vector field is V H (x) = A H(x). The Hamiltonian equation becomes ẋ(t) = A H(x(t)). For a time depending function H t (x) = H(t, x), the Hamiltonian vector field makes sense as V Ht (x) = A x H(t, x) and the Hamiltonian equation becomes ẋ(t) = A x H(t, x(t)). The following Darboux Theorem will be very useful
19 Theorem 4.1 (Darboux Theorem) Suppose the rank of the k k matrix A in (11) is 2m with k = 2m + l, then there is a coordinates transformation y = Bx such that m ( f h {F, H}(x) = f ) h (y) = f(y) T J h(y), q i p i p i q i i=1 where f(y) = F (B 1 y), h(y) = H(B 1 y) and ( y ) = (p 1,, p m, q 1,, q m, z 1,, z l ) T. That is to say J = BAB T Jm 0 =
20 Particularly, for the matrix A m defined in (10) with m 2N + 1, and the matrix B m defined by I n I n 0 0 B m = I n 0, I n I n I n I n there holds B m A m B T m = ( ) Am So by choosing matrix ( C m 1 satisfying ) C m 1 A m 1 Cm 1 T = J nl with ( l = m 1 2) and taking B m Cm 1 0 = there holds B 0 I m mb m A m (B mb m ) T Jnl 0 =. ( ) 0 0 ỹ By taking the coordinates transformation y = = B m z and h(t, y) = H(t, z), the system (10) becomes y m ỹ(t) = A m 1 ỹh(t, y), y m c. (12)
21 From the coordinates transformation, we see that y i = x i for 1 i m 1 and y m = x 1 x 2 + x m 1 + x m. So system (12) becomes as z(t) = A m 1 H(t, z(t)), (13) where H(t, z) = H(t, z, c x 1 +x 2 +x m 1 ) and z = (x 1,, x m 1 ). When choosing c = 0, any solution z(t) of (13) with z(0) = B m 1 z(1) determines a solution of system (10) with x m (t) = x 1 (t) + x 2 (t) + x m 1 (t), so x(t) = x 1 (t) is a solution of the delay system (9), where 0 I n I n 0 B m 1 = I n. I n I n I n I n So in this case, we also transform the problem into the problem (1) with J = 1 A m 1 and M = B m
22 Now by some arguments of index theory and variational method, we have the following statement: Theorem 4.2 Suppose V satisfies the following conditions (V1) There exists a constant C > 0 such that V (t, x) C, (t, x) [0, 1] R n. (V ) There exists a continuous symmetric matrix function C (t) such that V (t, x) = C (t)x + o( x ) uniformly in t R as x. Let C 0 (t) = V (t, 0). For m 2N C 0 (t) 0 0 B 0 (t) = 0 C 0 (t) 0..., 0 0 C 0 (t) and C (t) 0 0 B (t) = 0 C (t) C (t) Set (i 0, ν 0 ) = (i J M (B 0), ν J M (B 0)), (i, ν ) = (i J M (B ), ν J M (B )) with J = A m, M = T 1 m
23 For m 2N + 1, 2C 0 (t) C 0 (t) C 0 (t) C 0 (t) C 0 (t) B 0 (t) = C 0 (t) 2C 0 (t) C 0 (t)... C 0 (t). C 0 (t). C 0 (t) C 0 (t) C 0 (t) C 0 (t) 2C 0 (t) and B (t) = 2C (t) C (t) C (t) C (t) C (t) C (t) 2C (t) C (t) C (t) C (t)..... C (t) C (t) C (t) C (t) 2C (t). Set (i 0, ν 0 ) = (i J M (B 0), ν J M (B 0)), (i, ν ) = (i J M (B ), ν J M (B )) with J = 1 A m 1, M = B m
24 If i / [i 0, i 0 + ν 0 ], the delay differential system (9) possesses at least one nontrivial 2mτ-periodic solution x with x(t mτ) = x(t). Furthermore, if the trivial solution z 0 = 0 of problem (10), with H(t, z) = V (t, x 1 ) + + V (t, x m ) for even m and H(t, z) = V (t, x 1 ) + + V (t, x m 1 ) + V (t, x 1 + x 2 + x m 1 ) for odd m, is not pseudo-degenerated, and i / [i 0 2N, i 0 + ν 0 + 2N], the delay differential system (9) possesses at least two nontrivial 2mτ-periodic solutions as above
25 5 First order delay Hamiltonian systems For a function G C 2 (R R 2n, R) with G(t + τ, x) = G(t, x), we consider the 2τ periodic solutions of following first order delay Hamiltonian system ẋ(t) = J n G(t, x(t τ)). (14) For a 2τ period solution x(t) of (14), by setting x 1 (t) = x(t), x 2 (t) = x(t τ) and z = (x 1, x 2 ) T, the delay Hamiltonian system (14) is read as ż(t) = J 2n H(t, z(t)) (15) ( ) with H(t, z) = H(t, x 1, x 2 ) = G(t, x 1 ) + G(t, x 2 ) and J 0 Jn 2n =. ( J ) n 0 0 In If z(t) is a solution of (15) with z(0) = P 2n z(τ), P 2n =, then x(t) = I n 0 x 1 (t) is a 2τ-periodic solution of (14)
26 In general, we consider the following first order delay Hamiltonian systems ẋ(t) = J n ( G(t, x(t τ))+ G(t, x(t 2τ))+ + G(t, x(t (m 1)τ))). (16) For an mτ-periodic solution x(t) of (16), by setting x 1 (t) = x(t), x 2 (t) = x(t τ),, x m (t) = x(t (m 1)τ) and z = (x 1, x 2,, x m ) T, the delay Hamiltonian system (16) is read as ż(t) = J n,m H(t, z(t)) (17) with H(t, z) = G(t, x 1 ) + + G(t, x m ) and 0 J n J n J n J n,m = J n 0 J n J n..... J n J n J n
27 Conversely, if z(t) is a solution of (17) with z(0) = P n,m z(τ), where 0 I n 0 0 ( ) P n,m = 0 0 I n = In(m 1), I n 0 I n then x(t) = x 1 (t) is a mτ-periodic solution of (16)
28 Theorem 5.1 Suppose G satisfies the conditions (V1) and (V ) in Theorem4.2. Then the system (16) possesses an mτ-periodic solution x 0. Suppose z 0 (t) is the solution of (17) corresponding to x 0. Let B 0 (t) = H (t, z 0 (t)) and B (t) = C (t) C (t) C (t) Set (i 0, ν 0 ) = (i J M (B 0), ν J M (B 0)), (i, ν ) = (i J M (B ), ν J M (B )) with J = J n,m, M = Pn,m. 1 If i / [i 0, i 0 + ν 0 ] holds, the system (16) possesses at least two mτ-periodic solutions. Furthermore, if z 0 is not pseudo-degenerated, and i / [i 0 2N, i 0 + ν 0 + 2N] holds, then the system (16) possesses at least three an mτ-periodic solutions
29 6 Second order delay Hamiltonian systems For a function V C 2 (R R n, R) with V (t + τ, x) = V (t, x), we consider the periodic solutions of following second order delay Hamiltonian system ẍ(t) + V (t, x(t 2τ)) = 0. We can turn it into a first order delay Hamiltonian system as (14) with H(t, x, y) = y 2 2 V (t, x) and y(t) = ẋ(t + τ). In general, we consider the mτ-periodic solutions of the following second order delay system ẍ(t) = [ V (t, x(t τ)) + V (t, x(t 2τ)) + + V (t, x(t (m 1)τ))]. (18)
30 Let x 1 (t) = x(t), x 2 (t) = x(t τ), x m (t) = x(t (m 1)τ) and z(t) = (x 1 (t),, x m (t)) T, then by x(t + mτ) = x(t), there holds z(t) = A m H(t, z(t)), where H(t, z) = V (t, x 1 ) + + V (t, x m ) and 0 I n I n A m = I n 0 I n.... I n I n 0 It is easy to see that det A m 0. By taking y(t) = A 1 m ż(t) and w = (z, y) T, there holds ẇ(t) = J mn H(t, w(t)), (19) where H(t, w) = 1 2 (A my, y) H(t, z)
31 A solution ( of this Hamiltonian ) system with the boundary value condition Pn,m 0 w(0) = 0 A 1 w(τ) determine an mτ-periodic solution x(t) = m P n,m A m x 1 (t) of the second order delay system. We also turn our problem into the problem (1) with J = J mn and M = Pn,m. 1 By the same arguments, we also obtained an existence and multiplicity result as above theorems for this case
32 Thank you! liucg@nankai.edu.cn
Symplectic and Poisson Manifolds
Symplectic and Poisson Manifolds Harry Smith In this survey we look at the basic definitions relating to symplectic manifolds and Poisson manifolds and consider different examples of these. We go on to
More informationLecture 11. Andrei Antonenko. February 26, Last time we studied bases of vector spaces. Today we re going to give some examples of bases.
Lecture 11 Andrei Antonenko February 6, 003 1 Examples of bases Last time we studied bases of vector spaces. Today we re going to give some examples of bases. Example 1.1. Consider the vector space P the
More informationRandom and Deterministic perturbations of dynamical systems. Leonid Koralov
Random and Deterministic perturbations of dynamical systems Leonid Koralov - M. Freidlin, L. Koralov Metastability for Nonlinear Random Perturbations of Dynamical Systems, Stochastic Processes and Applications
More informationChap. 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 informationYORK UNIVERSITY. Faculty of Science Department of Mathematics and Statistics MATH M Test #2 Solutions
YORK UNIVERSITY Faculty of Science Department of Mathematics and Statistics MATH 3. M Test # Solutions. (8 pts) For each statement indicate whether it is always TRUE or sometimes FALSE. Note: For this
More informationx R d, λ R, f smooth enough. Goal: compute ( follow ) equilibrium solutions as λ varies, i.e. compute solutions (x, λ) to 0 = f(x, λ).
Continuation of equilibria Problem Parameter-dependent ODE ẋ = f(x, λ), x R d, λ R, f smooth enough. Goal: compute ( follow ) equilibrium solutions as λ varies, i.e. compute solutions (x, λ) to 0 = f(x,
More informationStabilization 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 informationEcon Lecture 14. Outline
Econ 204 2010 Lecture 14 Outline 1. Differential Equations and Solutions 2. Existence and Uniqueness of Solutions 3. Autonomous Differential Equations 4. Complex Exponentials 5. Linear Differential Equations
More informationBACKGROUND 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 informationMore chapter 3...linear dependence and independence... vectors
More chapter 3...linear dependence and independence... vectors It is important to determine if a set of vectors is linearly dependent or independent Consider a set of vectors A, B, and C. If we can find
More informationPeriod Domains. Carlson. June 24, 2010
Period Domains Carlson June 4, 00 Carlson - Period Domains Period domains are parameter spaces for marked Hodge structures. We call Γ\D the period space, which is a parameter space of isomorphism classes
More informationNORMS ON SPACE OF MATRICES
NORMS ON SPACE OF MATRICES. Operator Norms on Space of linear maps Let A be an n n real matrix and x 0 be a vector in R n. We would like to use the Picard iteration method to solve for the following system
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationSymplectic Geometry versus Riemannian Geometry
Symplectic Geometry versus Riemannian Geometry Inês Cruz, CMUP 2010/10/20 - seminar of PIUDM 1 necessarily very incomplete The aim of this talk is to give an overview of Symplectic Geometry (there is no
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationLecture 4: Harmonic forms
Lecture 4: Harmonic forms Jonathan Evans 29th September 2010 Jonathan Evans () Lecture 4: Harmonic forms 29th September 2010 1 / 15 Jonathan Evans () Lecture 4: Harmonic forms 29th September 2010 2 / 15
More informationu n 2 4 u n 36 u n 1, n 1.
Exercise 1 Let (u n ) be the sequence defined by Set v n = u n 1 x+ u n and f (x) = 4 x. 1. Solve the equations f (x) = 1 and f (x) =. u 0 = 0, n Z +, u n+1 = u n + 4 u n.. Prove that if u n < 1, then
More informationSpectral properties of Toeplitz+Hankel Operators
Spectral properties of Toeplitz+Hankel Operators Torsten Ehrhardt University of California, Santa Cruz ICM Satellite Conference on Operator Algebras and Applications Cheongpung, Aug. 8-12, 2014 Overview
More informationINSTABILITY OF RAPIDLY-OSCILLATING PERIODIC SOLUTIONS FOR DISCONTINUOUS DIFFERENTIAL DELAY EQUATIONS
Differential and Integral Equations Volume 15, Number 1, January 2002, Pages 53 90 INSTABILITY OF RAPIDLY-OSCILLATING PERIODIC SOLUTIONS FOR DISCONTINUOUS DIFFERENTIAL DELAY EQUATIONS Marianne Akian INRIA,
More informationExercises in Geometry II University of Bonn, Summer semester 2015 Professor: Prof. Christian Blohmann Assistant: Saskia Voss Sheet 1
Assistant: Saskia Voss Sheet 1 1. Conformal change of Riemannian metrics [3 points] Let (M, g) be a Riemannian manifold. A conformal change is a nonnegative function λ : M (0, ). Such a function defines
More information1 Contact structures and Reeb vector fields
1 Contact structures and Reeb vector fields Definition 1.1. Let M be a smooth (2n + 1)-dimensional manifold. A contact structure on M is a maximally non-integrable hyperplane field ξ T M. Suppose that
More informationLECTURE 2: SYMPLECTIC VECTOR BUNDLES
LECTURE 2: SYMPLECTIC VECTOR BUNDLES WEIMIN CHEN, UMASS, SPRING 07 1. Symplectic Vector Spaces Definition 1.1. A symplectic vector space is a pair (V, ω) where V is a finite dimensional vector space (over
More informationSome recent results on controllability of coupled parabolic systems: Towards a Kalman condition
Some recent results on controllability of coupled parabolic systems: Towards a Kalman condition F. Ammar Khodja Clermont-Ferrand, June 2011 GOAL: 1 Show the important differences between scalar and non
More informationx 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7
Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)
More informationA proof of the Jordan normal form theorem
A proof of the Jordan normal form theorem Jordan normal form theorem states that any matrix is similar to a blockdiagonal matrix with Jordan blocks on the diagonal. To prove it, we first reformulate it
More informationFlat Bi-Hamiltonian Structures and Invariant Densities
Lett Math Phys (2016) 106:1415 1427 DOI 10.1007/s11005-016-0875-1 Flat Bi-Hamiltonian Structures and Invariant Densities ANTON IZOSIMOV Department of Mathematics, University of Toronto, 40 St. George Street,
More informationLinear Ordinary Differential Equations
MTH.B402; Sect. 1 20180703) 2 Linear Ordinary Differential Equations Preliminaries: Matrix Norms. Denote by M n R) the set of n n matrix with real components, which can be identified the vector space R
More informationDirac structures. Henrique Bursztyn, IMPA. Geometry, mechanics and dynamics: the legacy of J. Marsden Fields Institute, July 2012
Dirac structures Henrique Bursztyn, IMPA Geometry, mechanics and dynamics: the legacy of J. Marsden Fields Institute, July 2012 Outline: 1. Mechanics and constraints (Dirac s theory) 2. Degenerate symplectic
More informationMcGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA
McGill University Department of Mathematics and Statistics Ph.D. preliminary examination, PART A PURE AND APPLIED MATHEMATICS Paper BETA 17 August, 2018 1:00 p.m. - 5:00 p.m. INSTRUCTIONS: (i) This paper
More informationA Z N -graded generalization of the Witt algebra
A Z N -graded generalization of the Witt algebra Kenji IOHARA (ICJ) March 5, 2014 Contents 1 Generalized Witt Algebras 1 1.1 Background............................ 1 1.2 A generalization of the Witt algebra..............
More informationAbstract. 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 informationNon-associative Deformations of Geometry in Double Field Theory
Non-associative Deformations of Geometry in Double Field Theory Michael Fuchs Workshop Frontiers in String Phenomenology based on JHEP 04(2014)141 or arxiv:1312.0719 by R. Blumenhagen, MF, F. Haßler, D.
More informationAvailable online at ScienceDirect. Procedia IUTAM 19 (2016 ) IUTAM Symposium Analytical Methods in Nonlinear Dynamics
Available online at www.sciencedirect.com ScienceDirect Procedia IUTAM 19 (2016 ) 11 18 IUTAM Symposium Analytical Methods in Nonlinear Dynamics A model of evolutionary dynamics with quasiperiodic forcing
More informationEcon 204 Differential Equations. 1 Existence and Uniqueness of Solutions
Econ 4 Differential Equations In this supplement, we use the methods we have developed so far to study differential equations. 1 Existence and Uniqueness of Solutions Definition 1 A differential equation
More informationTHE WAVE EQUATION. d = 1: D Alembert s formula We begin with the initial value problem in 1 space dimension { u = utt u xx = 0, in (0, ) R, (2)
THE WAVE EQUATION () The free wave equation takes the form u := ( t x )u = 0, u : R t R d x R In the literature, the operator := t x is called the D Alembertian on R +d. Later we shall also consider the
More informationEXERCISES 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 informationIntegrating exact Courant Algebroids
Integrating exact Courant Algebroids Rajan Mehta (Joint with Xiang Tang) Smith College September 28, 2013 The Courant bracket The Courant bracket is a bracket on Γ(TM T M), given by Twisted version: [X
More informationComputations of Critical Groups at a Degenerate Critical Point for Strongly Indefinite Functionals
Journal of Mathematical Analysis and Applications 256, 462 477 (2001) doi:10.1006/jmaa.2000.7292, available online at http://www.idealibrary.com on Computations of Critical Groups at a Degenerate Critical
More informationSturm-Liouville Matrix Differential Systems with Singular Leading Coefficients
Sturm-Liouville Matrix Differential Systems with Singular Leading Coefficients Iva Drimalová, Werner Kratz und Roman Simon Hilscher Preprint Series: 2016-05 Fakultät für Mathematik und Wirtschaftswissenschaften
More informationLECTURE 16: LIE GROUPS AND THEIR LIE ALGEBRAS. 1. Lie groups
LECTURE 16: LIE GROUPS AND THEIR LIE ALGEBRAS 1. Lie groups A Lie group is a special smooth manifold on which there is a group structure, and moreover, the two structures are compatible. Lie groups are
More informationLECTURE 25-26: CARTAN S THEOREM OF MAXIMAL TORI. 1. Maximal Tori
LECTURE 25-26: CARTAN S THEOREM OF MAXIMAL TORI 1. Maximal Tori By a torus we mean a compact connected abelian Lie group, so a torus is a Lie group that is isomorphic to T n = R n /Z n. Definition 1.1.
More informationLECTURE 26: THE CHERN-WEIL THEORY
LECTURE 26: THE CHERN-WEIL THEORY 1. Invariant Polynomials We start with some necessary backgrounds on invariant polynomials. Let V be a vector space. Recall that a k-tensor T k V is called symmetric if
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationOn periodic solutions of superquadratic Hamiltonian systems
Electronic Journal of Differential Equations, Vol. 22(22), No. 8, pp. 1 12. ISSN: 172-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) On periodic solutions
More informationEntrance 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 information11 a 12 a 21 a 11 a 22 a 12 a 21. (C.11) A = The determinant of a product of two matrices is given by AB = A B 1 1 = (C.13) and similarly.
C PROPERTIES OF MATRICES 697 to whether the permutation i 1 i 2 i N is even or odd, respectively Note that I =1 Thus, for a 2 2 matrix, the determinant takes the form A = a 11 a 12 = a a 21 a 11 a 22 a
More informationTHE GEOMETRY OF B-FIELDS. Nigel Hitchin (Oxford) Odense November 26th 2009
THE GEOMETRY OF B-FIELDS Nigel Hitchin (Oxford) Odense November 26th 2009 THE B-FIELD IN PHYSICS B = i,j B ij dx i dx j flux: db = H a closed three-form Born-Infeld action: det(g ij + B ij ) complexified
More informationOSCILLATION AND SPECTRAL THEORY FOR LINEAR HAMILTONIAN SYSTEMS WITH NONLINEAR DEPENDENCE ON THE SPECTRAL PARAMETER
OSCILLATION AND SPECTRAL THEORY FOR LINEAR HAMILTONIAN SYSTEMS WITH NONLINEAR DEPENDENCE ON THE SPECTRAL PARAMETER Martin Bohner, Werner Kratz und Roman Simon Hilscher Preprint Series: 2011-09 Fakultät
More informationProblem set 2. Math 212b February 8, 2001 due Feb. 27
Problem set 2 Math 212b February 8, 2001 due Feb. 27 Contents 1 The L 2 Euler operator 1 2 Symplectic vector spaces. 2 2.1 Special kinds of subspaces....................... 3 2.2 Normal forms..............................
More informationOn implicit Lagrangian differential systems
ANNALES POLONICI MATHEMATICI LXXIV (2000) On implicit Lagrangian differential systems by S. Janeczko (Warszawa) Bogdan Ziemian in memoriam Abstract. Let (P, ω) be a symplectic manifold. We find an integrability
More informationBOUNDED WEAK SOLUTION FOR THE HAMILTONIAN SYSTEM. Q-Heung Choi and Tacksun Jung
Korean J. Math. 2 (23), No., pp. 8 9 http://dx.doi.org/.568/kjm.23.2..8 BOUNDED WEAK SOLUTION FOR THE HAMILTONIAN SYSTEM Q-Heung Choi and Tacksun Jung Abstract. We investigate the bounded weak solutions
More informationREVIEW 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 informationMath 115A: Homework 5
Math 115A: Homework 5 1 Suppose U, V, and W are finite-dimensional vector spaces over a field F, and that are linear a) Prove ker ST ) ker T ) b) Prove nullst ) nullt ) c) Prove imst ) im S T : U V, S
More informationGraphs With Many Valencies and Few Eigenvalues
Graphs With Many Valencies and Few Eigenvalues By Edwin van Dam, Jack H. Koolen, Zheng-Jiang Xia ELA Vol 28. 2015 Vicente Valle Martinez Math 595 Spectral Graph Theory Apr 14, 2017 V. Valle Martinez (Math
More informationEACAT7 at IISER, Mohali Homology of 4D universe for every 3-manifold
EACAT7 at IISER, Mohali Homology of 4D universe for every 3-manifold Akio Kawauchi Osaka City University Advanced Mathematical Institute December 2, 2017 1. Types I, II, full and punctured 4D universes
More informationStationary trajectories, singular Hamiltonian systems and ill-posed Interconnection
Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection S.C. Jugade, Debasattam Pal, Rachel K. Kalaimani and Madhu N. Belur Department of Electrical Engineering Indian Institute
More informationHYPERKÄHLER MANIFOLDS
HYPERKÄHLER MANIFOLDS PAVEL SAFRONOV, TALK AT 2011 TALBOT WORKSHOP 1.1. Basic definitions. 1. Hyperkähler manifolds Definition. A hyperkähler manifold is a C Riemannian manifold together with three covariantly
More informationEquivalence of dynamical systems by bisimulation
Equivalence of dynamical systems by bisimulation Arjan van der Schaft Department of Applied Mathematics, University of Twente P.O. Box 217, 75 AE Enschede, The Netherlands Phone +31-53-4893449, Fax +31-53-48938
More informationOn the Virtual Fundamental Class
On the Virtual Fundamental Class Kai Behrend The University of British Columbia Seoul, August 14, 2014 http://www.math.ubc.ca/~behrend/talks/seoul14.pdf Overview Donaldson-Thomas theory: counting invariants
More informationVector fields in the presence of a contact structure
Vector fields in the presence of a contact structure Valentin Ovsienko To cite this version: Valentin Ovsienko. Vector fields in the presence of a contact structure. Preprint ICJ. 10 pages. 2005.
More informationAn 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 informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 1. Linear systems Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Linear systems 1.1 Differential Equations 1.2 Linear flows 1.3 Linear maps
More informationThe dynamical rigid body with memory
The dynamical rigid body with memory Ion Doru Albu, Mihaela Neamţu and Dumitru Opriş Abstract. In the present paper we describe the dynamics of the revised rigid body, the dynamics of the rigid body with
More informationREMARKS ON THE TIME-OPTIMAL CONTROL OF A CLASS OF HAMILTONIAN SYSTEMS. Eduardo D. Sontag. SYCON - Rutgers Center for Systems and Control
REMARKS ON THE TIME-OPTIMAL CONTROL OF A CLASS OF HAMILTONIAN SYSTEMS Eduardo D. Sontag SYCON - Rutgers Center for Systems and Control Department of Mathematics, Rutgers University, New Brunswick, NJ 08903
More informationRepresentations and Linear Actions
Representations and Linear Actions Definition 0.1. Let G be an S-group. A representation of G is a morphism of S-groups φ G GL(n, S) for some n. We say φ is faithful if it is a monomorphism (in the category
More informationMañé s Conjecture from the control viewpoint
Mañé s Conjecture from the control viewpoint Université de Nice - Sophia Antipolis Setting Let M be a smooth compact manifold of dimension n 2 be fixed. Let H : T M R be a Hamiltonian of class C k, with
More informationA Model of Evolutionary Dynamics with Quasiperiodic Forcing
paper presented at Society for Experimental Mechanics (SEM) IMAC XXXIII Conference on Structural Dynamics February 2-5, 205, Orlando FL A Model of Evolutionary Dynamics with Quasiperiodic Forcing Elizabeth
More informationLie groupoids and Lie algebroids
Lie groupoids and Lie algebroids Songhao Li University of Toronto, Canada University of Waterloo July 30, 2013 Li (Toronto) Groupoids and algebroids Waterloo 2013 1 / 22 Outline Lie groupoid Lie groupoid
More informationSYMPLECTIC MANIFOLDS, GEOMETRIC QUANTIZATION, AND UNITARY REPRESENTATIONS OF LIE GROUPS. 1. Introduction
SYMPLECTIC MANIFOLDS, GEOMETRIC QUANTIZATION, AND UNITARY REPRESENTATIONS OF LIE GROUPS CRAIG JACKSON 1. Introduction Generally speaking, geometric quantization is a scheme for associating Hilbert spaces
More informationLecture 5. Numerical continuation of connecting orbits of iterated maps and ODEs. Yu.A. Kuznetsov (Utrecht University, NL)
Lecture 5 Numerical continuation of connecting orbits of iterated maps and ODEs Yu.A. Kuznetsov (Utrecht University, NL) May 26, 2009 1 Contents 1. Point-to-point connections. 2. Continuation of homoclinic
More informationSymmetric and self-adjoint matrices
Symmetric and self-adjoint matrices A matrix A in M n (F) is called symmetric if A T = A, ie A ij = A ji for each i, j; and self-adjoint if A = A, ie A ij = A ji or each i, j Note for A in M n (R) that
More informationSEMILINEAR ELLIPTIC EQUATIONS WITH DEPENDENCE ON THE GRADIENT
Electronic Journal of Differential Equations, Vol. 2012 (2012), No. 139, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SEMILINEAR ELLIPTIC
More informationUniformly accurate averaging numerical schemes for oscillatory evolution equations
Uniformly accurate averaging numerical schemes for oscillatory evolution equations Philippe Chartier University of Rennes, INRIA Joint work with M. Lemou (University of Rennes-CNRS), F. Méhats (University
More informationLeast Squares Based Self-Tuning Control Systems: Supplementary Notes
Least Squares Based Self-Tuning Control Systems: Supplementary Notes S. Garatti Dip. di Elettronica ed Informazione Politecnico di Milano, piazza L. da Vinci 32, 2133, Milan, Italy. Email: simone.garatti@polimi.it
More informationThe spectral zeta function
The spectral zeta function Bernd Ammann June 4, 215 Abstract In this talk we introduce spectral zeta functions. The spectral zeta function of the Laplace-Beltrami operator was already introduced by Minakshisundaram
More informationHAMILTONIAN ACTIONS IN GENERALIZED COMPLEX GEOMETRY
HAMILTONIAN ACTIONS IN GENERALIZED COMPLEX GEOMETRY TIMOTHY E. GOLDBERG These are notes for a talk given in the Lie Groups Seminar at Cornell University on Friday, September 25, 2009. In retrospect, perhaps
More informationNonlinear Systems Theory
Nonlinear Systems Theory Matthew M. Peet Arizona State University Lecture 2: Nonlinear Systems Theory Overview Our next goal is to extend LMI s and optimization to nonlinear systems analysis. Today we
More informationMATH 263: PROBLEM SET 1: BUNDLES, SHEAVES AND HODGE THEORY
MATH 263: PROBLEM SET 1: BUNDLES, SHEAVES AND HODGE THEORY 0.1. Vector Bundles and Connection 1-forms. Let E X be a complex vector bundle of rank r over a smooth manifold. Recall the following abstract
More informationDefinition 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 informationRepresentations of Sp(6,R) and SU(3) carried by homogeneous polynomials
Representations of Sp(6,R) and SU(3) carried by homogeneous polynomials Govindan Rangarajan a) Department of Mathematics and Centre for Theoretical Studies, Indian Institute of Science, Bangalore 560 012,
More informationOn a generalized Sturm theorem
On a generalized Sturm theorem Alessandro Portaluri May 21, 29 Abstract Sturm oscillation theorem for second order differential equations was generalized to systems and higher order equations with positive
More informationTHE EULER CHARACTERISTIC OF A LIE GROUP
THE EULER CHARACTERISTIC OF A LIE GROUP JAY TAYLOR 1 Examples of Lie Groups The following is adapted from [2] We begin with the basic definition and some core examples Definition A Lie group is a smooth
More informationInvariant 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 informationInstability of Periodic Orbits of Some Rhombus and Parallelogram Four Body Problems
Instability of Periodic Orbits of Some Rhombus and Parallelogram Four Body Problems by Abdalla Mansur A thesis submitted to the Department of Mathematics & Statistics in conformity with the requirements
More informationRIESZ BASES AND UNCONDITIONAL BASES
In this paper we give a brief introduction to adjoint operators on Hilbert spaces and a characterization of the dual space of a Hilbert space. We then introduce the notion of a Riesz basis and give some
More informationOn the dynamics of strongly tridiagonal competitive-cooperative system
On the dynamics of strongly tridiagonal competitive-cooperative system Chun Fang University of Helsinki International Conference on Advances on Fractals and Related Topics, Hong Kong December 14, 2012
More informationA short introduction with a view toward examples. Short presentation for the students of: Dynamical Systems and Complexity Summer School Volos, 2017
A short introduction with a view toward examples Center of Research and Applications of Nonlinear (CRANS) Department of Mathematics University of Patras Greece sanastassiou@gmail.com Short presentation
More informationEstimates for the resolvent and spectral gaps for non self-adjoint operators. University Bordeaux
for the resolvent and spectral gaps for non self-adjoint operators 1 / 29 Estimates for the resolvent and spectral gaps for non self-adjoint operators Vesselin Petkov University Bordeaux Mathematics Days
More informationIntroduction to the Baum-Connes conjecture
Introduction to the Baum-Connes conjecture Nigel Higson, John Roe PSU NCGOA07 Nigel Higson, John Roe (PSU) Introduction to the Baum-Connes conjecture NCGOA07 1 / 15 History of the BC conjecture Lecture
More informationCanonical Forms for BiHamiltonian Systems
Canonical Forms for BiHamiltonian Systems Peter J. Olver Dedicated to the Memory of Jean-Louis Verdier BiHamiltonian systems were first defined in the fundamental paper of Magri, [5], which deduced the
More informationAn Invitation to Geometric Quantization
An Invitation to Geometric Quantization Alex Fok Department of Mathematics, Cornell University April 2012 What is quantization? Quantization is a process of associating a classical mechanical system to
More informationA Problem of Hsiang-Palais-Terng on Isoparametric Submanifolds
arxiv:math/0312251v1 [math.dg] 12 Dec 2003 A Problem of Hsiang-Palais-Terng on Isoparametric Submanifolds Haibao Duan Institute of Mathematics, Chinese Academy of Sciences, Beijing 100080, dhb@math.ac.cn
More informationMATH 23a, FALL 2002 THEORETICAL LINEAR ALGEBRA AND MULTIVARIABLE CALCULUS Solutions to Final Exam (in-class portion) January 22, 2003
MATH 23a, FALL 2002 THEORETICAL LINEAR ALGEBRA AND MULTIVARIABLE CALCULUS Solutions to Final Exam (in-class portion) January 22, 2003 1. True or False (28 points, 2 each) T or F If V is a vector space
More informationMath 594, HW2 - Solutions
Math 594, HW2 - Solutions Gilad Pagi, Feng Zhu February 8, 2015 1 a). It suffices to check that NA is closed under the group operation, and contains identities and inverses: NA is closed under the group
More informationPeriods of meromorphic quadratic differentials and Goldman bracket
Periods of meromorphic quadratic differentials and Goldman bracket Dmitry Korotkin Concordia University, Montreal Geometric and Algebraic aspects of integrability, August 05, 2016 References D.Korotkin,
More informationMASLOV INDEX. Contents. Drawing from [1]. 1. Outline
MASLOV INDEX J-HOLOMORPHIC CURVE SEMINAR MARCH 2, 2015 MORGAN WEILER 1. Outline 1 Motivation 1 Definition of the Maslov Index 2 Questions Demanding Answers 2 2. The Path of an Orbit 3 Obtaining A 3 Independence
More informationHolomorphic symplectic fermions
Holomorphic symplectic fermions Ingo Runkel Hamburg University joint with Alexei Davydov Outline Investigate holomorphic extensions of symplectic fermions via embedding into a holomorphic VOA (existence)
More informationHamiltonian Dynamics
Hamiltonian Dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS Feb. 10, 2009 Joris Vankerschaver (CDS) Hamiltonian Dynamics Feb. 10, 2009 1 / 31 Outline 1. Introductory concepts; 2. Poisson brackets;
More informationTwisted Poisson manifolds and their almost symplectically complete isotropic realizations
Twisted Poisson manifolds and their almost symplectically complete isotropic realizations Chi-Kwong Fok National Center for Theoretical Sciences Math Division National Tsing Hua University (Joint work
More informationAubry Mather Theory from a Topological Viewpoint
Aubry Mather Theory from a Topological Viewpoint III. Applications to Hamiltonian instability Marian Gidea,2 Northeastern Illinois University, Chicago 2 Institute for Advanced Study, Princeton WORKSHOP
More information