USING CARLEMAN EMBEDDING TO DISCOVER A SYSTEM S MOTION CONSTANTS
|
|
- Brent Holland
- 5 years ago
- Views:
Transcription
1 (Preprint) AAS USING CARLEMAN EMBEDDING TO DISCOVER A SYSTEM S MOTION CONSTANTS John E. Hurtado and Andrew J. Sinclair INTRODUCTION Although the solutions with respect to time are commonly sought when analyzing the behavior of dynamic systems, there are other expressions that can be as revealing and important. Specifically, motion constants, which are time independent algebraic or transcendental equations involving the system states, can provide a broad view of a system s motion. Locally, at least, the motion constants of autonomous systems exist. Nevertheless, for finite dimensional, autonomous, nonlinear systems, the motion constants are generally difficult to identify. The motion constants for every finite dimensional, autonomous, linear system, regardless of its dimension, however, have recently been found. This paper attempts to discover the motion constants of nonlinear systems through the motion constants of a linear representation. The linear representation is achieved through a Carleman embedding, which provides a procedure to transform a finite dimensional system of analytic ordinary differential equations into an infinite system of linear differential equations The solutions of dynamic systems are commonly presented in terms of collections of integrals. These integrals can be either explicitly time dependent or time independent. Time-independent integrals are able to describe the shape of the state-space trajectories and provide a broad picture of the motion irrespective of time. A complete set of time-independent integrals reduces the solution process to finding the speed along these trajectories [Contopoulas 2002]. Here, time-independent integrals are referred to as motion constants. These are algebraic or transcendental equations involving the system states. Sinclair and Hurtado [2012] have recently found exhaustive solutions for the motion constants of every finite dimensional, autonomous, linear system. Significantly, a central feature of their process is the single motion constant for a 2 2 system. For an n-dimensional linear system, the n 1 motion constants are constructed through a coordinate transformation to a Jordan canonical form. The diagonal block structure of the Jordan form allows one to identify motion constants that are internal to each block by using patterns learned from simple 2 2 groups. The remaining motion constants of the entire system are found by relating the time behavior of separate blocks to one another. With the exception of the case of incommensurate oscillators, the procedure provides a complete set of motion constants. Carleman showed in 1932 that a finite dimensional system of analytic ordinary differential equations can be embedded within an infinite system of linear differential equations [Steeb 1989, Kowalski and Steeb 1991]. The Carleman embedding process is well understood, and researchers have Associate Professor, Aerospace Engineering Department, Texas A&M University, 3141 TAMU, College Station, TX Associate Professor, Aerospace Engineering Department, Auburn University, 211 Davis Hall, Auburn, AL
2 answered important questions like, what is the relationship between the solutions of a finite dimensional nonlinear system and the solutions of its infinite dimensional linear representation? Here, the Carleman embedding process is investigated further yet, by seeking to understand the relationship between the motion constants of a finite dimensional nonlinear system and the motion constants of its infinite dimensional linear representation. INTEGRALS AND MOTION CONSTANTS OF DYNAMIC SYSTEMS Nonlinear Autonomous Systems An nth-order autonomous dynamic system describes the evolution of a state vector x(t) R n. ẋ = f(x) (1) This system of differential equations may be solved by finding n functions of the states and possibly time that are constant along the trajectory associated with an initial condition x 0. Such functions are called integrals [Szebehely 1967]. φ k (x(t), (t)) = φ k (x 0, 0) for t R, k = 1,..., n (2) Although the value of each φ k is constant along the motion, evaluating these integrals may require knowledge of the current states and time. Separately, an nth-order autonomous dynamic system can have up to n 1 time-independent integrals. ψ k (x(t)) = ψ k (x 0 ) for t R, k = 1,..., n 1 (3) Because each ψ k is constant, its Lie derivative along the direction of motion is zero. ψ k = ( ψ k / x) T ẋ = ( ψ k / x) T f(x) = 0 (4) This equation represents a partial differential equation for an unknown time independent function ψ k ; the ψ k are called motion constants. Although Frobenius theorem conveniently guarantees the existence of n 1 local solutions, determining the motion constants for nonlinear autonomous systems can be challenging. Linear Autonomous Systems Linear autonomous systems are a special class of eq. (1). ẋ = Ax (5) These systems are completely integrable in the sense of possessing n independent analytic integrals. φ(x, t) = exp At x (6) But these systems also possess motion constants ψ k (x), and their explicit forms have recently been discovered [Sinclair and Hurtado 2012]. A discussion of these forms can not be adequately addressed in this note; nevertheless, the important aspects are that the motion constants of n- dimensional linear systems can be assembled in a systematic way using Jordan canonical forms and a thorough comprehension of the single motion constant of 2 2 systems. 2
3 CARLEMAN EMBEDDING In simple terms, Carleman embedding is a process wherein a finite dimensional nonlinear system is mapped to a higher (perhaps infinitely higher) dimensional linear representation via coordinate transformations. A popular illustrative example is the rate equation [Kowalski and Steeb 1991]. Suppose new coordinates v are introduced. ẋ = x + x 2 (7) v 1 = x, v 2 = x 2,..., v n = x n (8) The differential equations that govern the new coordinates can be computed. v n = n x n 1 ẋ = n x n 1 ( x + x 2 ) = n x n + n x n+1 = n v n + n v n+1, n = 1, (9) The single nonlinear autonomous system has been mapped to an infinite collection of linear autonomous systems. The behavior and characteristics of the new infinite dimensional linear system must reflect the behavior and characteristics of the original finite dimensional nonlinear system in some way. Here, we pursue an understanding of the nonlinear system through the motion constants of the new linear representation; this pursuit is hopeful because the motion constants of every autonomous linear systems are known. At its heart, the Carleman embedding technique represents an overparameterization of the original system. We suspect that motion constants of the original system are captured in the overparameterization; more specifically, we suspect that the motion constants are captured at some finite truncation. This thinking leads one to ponder over the meaning of any extraneous motion constants. This brings to mind that one quality of overparameterizations is that natural constraint relationships are created between the new states [Hurtado and Sinclair 2011], thus we suspect that the extraneous motion constants reflect these natural constraints. A ROTATIONAL KINEMATICS INVESTIGATION The idea of discovering the motion constants of a nonlinear autonomous system through a coordinate embedding is examined through an example. Consider the rotational kinematics of a rigid body wherein the body-fixed components of the angular velocity vector are known constants. This system takes the form of a nonlinear autonomous system if the classic Rodrigues parameters are selected as the attitude coordinates. ẋ = f(x) ẋ = x2 1 x 1 x 2 x 3 x 1 x 3 + x 2 x 2 x 1 + x x 2 2 x 2 x 3 x 1 ω (10) 2 x 3 x 1 x 2 x 3 x 2 + x x 2 3 Consider an overparameterization that maps the three current coordinates x to four new coordinates β. [β 0, β 1, β 2, β 3 ] = [1, x 1, x 2, x 3 ]/ 1 + x T x (11) The kinematics in terms of the new coordinates are linear. β = Aβ where A = ω 1 ω 2 ω 3 ω 1 0 ω 3 ω 2 ω 2 ω 3 0 ω 1 ω 3 ω 2 ω 1 0 (12) 3
4 This state matrix has four eigenvalues, λ 1,2 = ±ı ω T ω/2, λ 3,4 = ±ı ω T ω/2, with four independent eigenvectors. The substitution ω 0 ω T ω is helpful in the ensuing equations. Following the process laid out by Sinclair and Hurtado, a Jordan canonical form is formed via A = MJM 1 where M is composed from the eigenvectors and J contains 2 2 blocks of zeros and the imaginary parts of the eigenvalues on the sub and super diagonals. M = ω 0 ω 2 /(ω1 2 + ω2 2 ) ω 1ω 3 /(ω1 2 + ω2 2 ) ω 0ω 1 /(ω1 2 + ω2 2 ) ω 2ω 3 /(ω1 2 + ω2 2 ) ω 0 ω 1 /(ω1 2 + ω2 2 ) ω 2ω 3 /(ω1 2 + ω2 2 ) ω 0ω 2 /(ω1 2 + ω2 2 ) ω 1ω 3 /(ω1 2 + ω2 2 ) J = ω ω ω ω 0 0 (13) (14) The matrix M is used to define Jordan coordinates y = M 1 β which are governed by ẏ = Jy. The Jordan system is partitioned so that motion constants are easily constructed. For this example, the three independent motion constants in terms of the four Jordan coordinates are given as follows. ψ 1 = y y 2 2 ; ψ 2 = y y 2 4 ; ψ 3 = y 2y 3 y 1 y 4 y 1 y 3 + y 2 y 4 (15) Another motion constant that is not independent of these three is ψ 0 = y T M T My. This motion constant is evident from eq. (11) which produces 1 = β T β. Consequently, any three of the four ψ k motion constants form an independent set. Judiciously selecting ψ 0, ψ 1, and ψ 2 and writing these in terms of the overparameterized β coordinates gives the following. ψ 0 = β β β β 2 3 (16) ψ 1 = ω 2 0β (ω 2 β 0 ω 3 β 1 + ω 1 β 3 ) 2 (17) ψ 2 = ω 2 0β (ω 1 β 0 + ω 3 β 2 ω 2 β 3 ) 2 (18) These can also be written in terms of the original (minimal) x coordinates. ψ 0 = 1 (19) ψ 1 = ω2 0 x2 2 + ω2 1 x2 3 + ω2 2 + ω2 3 x2 1 2ω 1ω 3 x 1 x 3 2ω 2 ω 3 x 1 + 2ω 1 ω 2 x 3 ω 2 0 (1 + x2 1 + x2 2 + x2 3 ) (20) ψ 2 = ω2 0 x2 1 + ω2 1 + ω2 2 x2 3 + ω2 3 x2 2 2ω 2ω 3 x 2 x 3 2ω 1 ω 2 x 3 + 2ω 1 ω 3 x 2 ω 2 0 (1 + x2 1 + x2 2 + x2 3 ) (21) One can confirm that functions ψ 1 and ψ 2 are constant along the trajectories of the dynamic system given by eq. (10), ( ψ k / x) T f(x) = 0. These equations illustrate that three independent motion constants for the overparameterized fourcoordinate system reflect two independent motion constants for the minimal three-coordinate system plus an extraneous motion constant that reflects the natural constraint involving the excess parameters. Two things to note, however, are that the coordinate embedding in this case is exact (i.e., a linear representation is exactly achieved with a finite overparameterization), and the natural constraint (1 = β T β) was not precisely found to be a motion constant of the linear representation, although it is functionally dependent on the extracted motion constants. 4
5 A ROTATIONAL DYNAMICS The examination of discovering the motion constants of a nonlinear autonomous system through a coordinate embedding is continued through a second example. A classic nonlinear autonomous system is the torque-free rotational dynamics of a rigid body. ω = f(ω) ω 1 = a 1 ω 2 ω 3 ; ω 2 = a 2 ω 1 ω 3 ; ω 3 = a 3 ω 1 ω 2 (22) Here, ω are the body-fixed components of the angular velocity vector and a k represent the rigid body inertia parameters. a 1 = (I 2 I 3 )/I 1 ; a 2 = (I 3 I 1 )/I 2 ; a 3 = (I 1 I 2 )/I 3 (23) It is well-known that this system has two time-independent integrals, which are the rotational kinetic energy and the magnitude of the angular momentum. T = 1 2 I 1ω I 2ω I 3ω 2 3 ; h 2 = I 2 1ω I 2 2ω I 2 3ω 2 3 (24) The Carleman embedding process to convert the three-dimensional nonlinear system to an infinitedimensional, overparameterized, linear system begins with the definition of new coordinates. x (ijk) = ω i 1 ω j 2 ωk 3, i, j, k = 1,..., (25) According to this definition, there are three x coordinates that are linear in the ω states; and six x coordinates that are quadratic in the ω states; and ten x coordinates that are cubic in the ω states; or (p + 1)(p + 2)/2 coordinates that are order p in the ω states. The new coordinates are governed by a system of linear autonomous equations. ẋ (ijk) = i a 1 x (i 1 j+1 k+1) + j a 2 x (i+1 j 1 k+1) + k a 3 x (i+1 j+1 k 1) (26) Our approach to discover two motion constants for the original nonlinear system is outlined in the following steps: 1. Truncate the length of the new coordinates at some point, say m. The collection of truncated states is denoted as x. 2. Retain the upper m m block of A. This is the truncated state matrix A. 3. Create a Jordan canonical form via A = MJM 1 where M is composed from the eigenvectors and J contains eigenvalues on the diagonal and sometimes ones on the subdiagonal. 4. Use M to define Jordan coordinates y = M 1 x which are governed by ẏ = Jy. 5. Use the methods discussed in Sinclair and Hurtado (2012) to determine the motion constants for the Jordan coordinates, ψ = ψ(y). 6. Use y = M 1 x to write the motion constants in terms of the truncated linear states x. 7. Use x = x(ω) x to write the motion constants in terms of the original coordinates ω. At this point, specific inertia values are used to facilitate the computations; we select a tri-inertial case I 1 = 1, I 2 = 2, and I 3 = 3. Regardless of the truncation limit, an eigenanalysis of the matrix A reveals that all of its eigenvalues equal zero. Some of the zero eigenvalues generate an independent eigenvector whereas others generate an eigenvector that is linearly dependent on the 5
6 others. Sinclair and Hurtado (2012) have shown that a Jordan coordinate associated to a 1 1 Jordan block with a zero eigenvalue is itself a motion constant. Subsequently using y = M 1 x for each of these special Jordan coordinates gives an expression for the motion constant in terms of the coordinates x. In this example we ll look for true motion constants of the original nonlinear system within these special Jordan coordinates. For example, suppose we truncate the linear representation to the coordinates x that are no more than quadratic in the ω states. This creates a nine element vector x : x 100, x 010, x 001, x 110, x 101, x 011, x 200, x 020, x 002 The 9 9 matrix A has three 1 1 Jordan blocks with a zero eigenvalue. The Jordan form is arranged so that the special Jordan coordinates that are constant for all time are y 1, y 2, and y 3. These three Jordan coordinates are motion constants, and using M they can be mapped to motion constants in terms of x, which are assumed to honor the Carleman coordinate definitions in eq. (25); this then gives motion constants in terms of ω. ψ 1 = y 1 = x x 200 = ω 1 ω 3 + ω 2 1 (27) ψ 2 = y 2 = x x 020 = ω ω 2 2 (28) ψ 3 = y 3 = x x 001 = ω ω 2 3 (29) Computing the Lie derivative shows that these expressions in ω are not motion constants for the true nonlinear system, so this truncation level does not capture the true motion constants of the original system. Moreover, unlike the previous example, none of the eight motion constants at this truncation level reflect the coordinate constraints that are part of the overparameterized definitions (see eq. (25) which can be interpreted to read as x (ijk) = x i 100 xj 010 xk 001 ). Truncating the linear representation to the coordinates x that are no more than quintic in the ω states gives a total of fifty-five elements for x. The matrix A has fourteen 1 1 Jordan blocks with a zero eigenvalue. The Jordan form is arranged so that the special Jordan coordinates that are constant for all time are y 1 through y 14. These fourteen Jordan coordinates are motion constants, and they can be mapped to motion constants in terms of x or ω using M. Two of these fourteen motion constants read as follows: ψ 1 = x 020 3x x x 022 x x 220 (30) = ω 2 2 3ω ω ω 2 2ω 2 3 ω 2 1ω ω2 1ω 2 2 (31) ψ 2 = x x x x x 400 (32) = ω ω ω ω 2 1ω ω 4 1 (33) The Lie derivative of these expressions in ω confirm that these are motion constants for the true nonlinear system. Moreover, the rotational kinetic energy T and magnitude of angular momentum h 2, see eq. (24), are not independent functions from ψ 1 and ψ 2. Therefore, this truncation level captures the true motion constants of the original system. Although motion constants for the true nonlinear system have been identified, the extraneous motion constants do not seem to reflect the coordinate constraints that are part of the overparameterized definitions. This is because each level of truncation neglects the influence of higher order states in the governing differential equations. 6
7 ω 1 ω 2 ω 3 x 100 x 010 x 001 ψ 1 ψ Figure 1. Graphs of a representative simulation. The solid lines show ω from a nonlinear simulation; the dashed lines show x 100, x 010, x 001 for a truncated linear representation; and the dotted lines show the two motion constants when evaluated using the x ( ijk) states. Interestingly, finding the exact motion constants for the nonlinear system through the Carleman embedding does not mean the time evolution is exactly matched. An illustration is shown in Figure 1. The solid lines show a two second evolution of ω from simulating the nonlinear equations of motion. The initial conditions are ω = [1.1, 1.2, 1.3] T. The dashed lines show the evolution of x 100, x 010, x 001 for a 55 state linear representation. Neglecting the affect of the truncation, these three states represent ω in the linear representation. Notice that the dashed lines diverge from the solid lines before much time has elasped. Also shown as dotted lines are the two motion constants when evaluated using the evolution of the linear states x ( ijk). That is, the dotted lines are plots of eqs. (30) and (32). CONCLUSIONS An idea for discovering the motion constants of a nonlinear autonomous system is presented in this paper. The idea uses an overparameterization, or Carleman embedding method, to map a finite dimensional nonlinear system to a higher dimensional linear representation via coordinate transformations. The idea further uses the result that the motion constants for every finite dimensional, autonomous, linear system, regardless of its dimension, have recently been found. We have shown through examples that the idea has merit: in each case we were able to discover exact motion constants for the original nonlinear system using overparameterization. When the overparameterization or embedding is exact, then the extraneous motion constants reflect the coordinate constraints that arise through the overparameterized coordinate definitions. This is as it should be because the governing equations are exact differentials of the overparameterized definitions. When the overparameterization or embedding is inexact (because of truncation), then the extraneous motion constants do not necessarily reflect the coordinate constraints because the truncation produces an approximation of the true governing equations. Thus the governing equations are no longer exact differentials of 7
8 the overparameterized definitions. REFERENCES [1] G. Contopoulos, Order and chaos in dynamical astronomy. Springer, [2] J.E. Hurtado, A.J. Sinclair, Lagrangian mechanics of overparameterized systems, Nonlinear Dynamics, Vol. 66, [3] K. Kowalski and W.-H. Steeb, Nonlinear dynamical systems and Carleman linearization. World Scientific Publishing Co Inc, [4] A.J. Sinclair and J.E. Hurtado, The time-independent integrals of linear autonomous dynamical systems, in preparation. [5] W.-H. Steeb, A note on Carleman linearization, Physics Letters A, Vol. 4, No. 6, [6] V. Szebehely, Theory of Orbits, Sects. 1.3, 1.8 & 2.1. Academic Press,
Linear Feedback Control Using Quasi Velocities
Linear Feedback Control Using Quasi Velocities Andrew J Sinclair Auburn University, Auburn, Alabama 36849 John E Hurtado and John L Junkins Texas A&M University, College Station, Texas 77843 A novel approach
More informationKeywords: Kinematics; Principal planes; Special orthogonal matrices; Trajectory optimization
Editorial Manager(tm) for Nonlinear Dynamics Manuscript Draft Manuscript Number: Title: Minimum Distance and Optimal Transformations on SO(N) Article Type: Original research Section/Category: Keywords:
More informationTWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY P. Yu 1,2
More informationMinimum-Parameter Representations of N-Dimensional Principal Rotations
Minimum-Parameter Representations of N-Dimensional Principal Rotations Andrew J. Sinclair and John E. Hurtado Department of Aerospace Engineering, Texas A&M University, College Station, Texas, USA Abstract
More informationConstruction of Lyapunov functions by validated computation
Construction of Lyapunov functions by validated computation Nobito Yamamoto 1, Kaname Matsue 2, and Tomohiro Hiwaki 1 1 The University of Electro-Communications, Tokyo, Japan yamamoto@im.uec.ac.jp 2 The
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More informationPhysics 221A Fall 2018 Notes 22 Bound-State Perturbation Theory
Copyright c 2018 by Robert G. Littlejohn Physics 221A Fall 2018 Notes 22 Bound-State Perturbation Theory 1. Introduction Bound state perturbation theory applies to the bound states of perturbed systems,
More informationJUST THE MATHS UNIT NUMBER 9.9. MATRICES 9 (Modal & spectral matrices) A.J.Hobson
JUST THE MATHS UNIT NUMBER 9.9 MATRICES 9 (Modal & spectral matrices) by A.J.Hobson 9.9. Assumptions and definitions 9.9.2 Diagonalisation of a matrix 9.9.3 Exercises 9.9.4 Answers to exercises UNIT 9.9
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationRigid bodies - general theory
Rigid bodies - general theory Kinetic Energy: based on FW-26 Consider a system on N particles with all their relative separations fixed: it has 3 translational and 3 rotational degrees of freedom. Motion
More informationPhysics 342 Lecture 27. Spin. Lecture 27. Physics 342 Quantum Mechanics I
Physics 342 Lecture 27 Spin Lecture 27 Physics 342 Quantum Mechanics I Monday, April 5th, 2010 There is an intrinsic characteristic of point particles that has an analogue in but no direct derivation from
More informationJohn E Hurtado Associate Professor of Aerospace Engineering, Texas A&M University, College Station, TX
John E Hurtado Associate Professor of Aerospace Engineering, Texas A&M University, College Station, TX http://aero.tamu.edu/faculty/hurtado/ Education Expertise Mentorship Teaching Awards Impact Scholarship
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear
More informationMath 1270 Honors ODE I Fall, 2008 Class notes # 14. x 0 = F (x; y) y 0 = G (x; y) u 0 = au + bv = cu + dv
Math 1270 Honors ODE I Fall, 2008 Class notes # 1 We have learned how to study nonlinear systems x 0 = F (x; y) y 0 = G (x; y) (1) by linearizing around equilibrium points. If (x 0 ; y 0 ) is an equilibrium
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
More informationAn Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems
An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems Scott Zimmerman MATH181HM: Dynamical Systems Spring 2008 1 Introduction The Hartman-Grobman and Poincaré-Bendixon Theorems
More information21 Linear State-Space Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
More informationDiagonalization of Matrices
LECTURE 4 Diagonalization of Matrices Recall that a diagonal matrix is a square n n matrix with non-zero entries only along the diagonal from the upper left to the lower right (the main diagonal) Diagonal
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationFinding eigenvalues for matrices acting on subspaces
Finding eigenvalues for matrices acting on subspaces Jakeniah Christiansen Department of Mathematics and Statistics Calvin College Grand Rapids, MI 49546 Faculty advisor: Prof Todd Kapitula Department
More information/. 3D. a very abstract sort of problem, this research is a step on the way to applicable knowledge.
/. 3D Corwin Atwood-Stone Physics Department, The College of Wooster, Wooster, Ohio 691, USA (Dated: May 7, 009) A simulation of the mutual orbital motion of a point and a stick of equal masses in three
More informationThe goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T
1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.
More information1. Matrix multiplication and Pauli Matrices: Pauli matrices are the 2 2 matrices. 1 0 i 0. 0 i
Problems in basic linear algebra Science Academies Lecture Workshop at PSGRK College Coimbatore, June 22-24, 2016 Govind S. Krishnaswami, Chennai Mathematical Institute http://www.cmi.ac.in/~govind/teaching,
More informationTranscendental cases in stability problem. Hamiltonian systems
of Hamiltonian systems Boris S. Bardin Moscow Aviation Institute (Technical University) Faculty of Applied Mathematics and Physics Department of Theoretical Mechanics Hamiltonian Dynamics and Celestial
More informationOn 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 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 informationSolutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.
Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the
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 informationRobotics. Control Theory. Marc Toussaint U Stuttgart
Robotics Control Theory Topics in control theory, optimal control, HJB equation, infinite horizon case, Linear-Quadratic optimal control, Riccati equations (differential, algebraic, discrete-time), controllability,
More informationPhase portraits in two dimensions
Phase portraits in two dimensions 8.3, Spring, 999 It [ is convenient to represent the solutions to an autonomous system x = f( x) (where x x = ) by means of a phase portrait. The x, y plane is called
More informationPHYSICS 110A : CLASSICAL MECHANICS
PHYSICS 110A : CLASSICAL MECHANICS 1. Introduction to Dynamics motion of a mechanical system equations of motion : Newton s second law ordinary differential equations (ODEs) dynamical systems simple 2.
More informationPhysics 106b/196b Problem Set 9 Due Jan 19, 2007
Physics 06b/96b Problem Set 9 Due Jan 9, 2007 Version 3: January 8, 2007 This problem set focuses on dynamics in rotating coordinate systems (Section 5.2), with some additional early material on dynamics
More informationElectric and Magnetic Forces in Lagrangian and Hamiltonian Formalism
Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Benjamin Hornberger 1/26/1 Phy 55, Classical Electrodynamics, Prof. Goldhaber Lecture notes from Oct. 26, 21 Lecture held by Prof. Weisberger
More informationJordan normal form notes (version date: 11/21/07)
Jordan normal form notes (version date: /2/7) If A has an eigenbasis {u,, u n }, ie a basis made up of eigenvectors, so that Au j = λ j u j, then A is diagonal with respect to that basis To see this, let
More informationIntroduction to Group Theory
Chapter 10 Introduction to Group Theory Since symmetries described by groups play such an important role in modern physics, we will take a little time to introduce the basic structure (as seen by a physicist)
More informationInvariant Manifolds of Dynamical Systems and an application to Space Exploration
Invariant Manifolds of Dynamical Systems and an application to Space Exploration Mateo Wirth January 13, 2014 1 Abstract In this paper we go over the basics of stable and unstable manifolds associated
More informationLight-Cone Quantization of Electrodynamics
Light-Cone Quantization of Electrodynamics David G. Robertson Department of Physics, The Ohio State University Columbus, OH 43210 Abstract Light-cone quantization of (3+1)-dimensional electrodynamics is
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
More information1/30. Rigid Body Rotations. Dave Frank
. 1/3 Rigid Body Rotations Dave Frank A Point Particle and Fundamental Quantities z 2/3 m v ω r y x Angular Velocity v = dr dt = ω r Kinetic Energy K = 1 2 mv2 Momentum p = mv Rigid Bodies We treat a rigid
More information22.2. Applications of Eigenvalues and Eigenvectors. Introduction. Prerequisites. Learning Outcomes
Applications of Eigenvalues and Eigenvectors 22.2 Introduction Many applications of matrices in both engineering and science utilize eigenvalues and, sometimes, eigenvectors. Control theory, vibration
More informationSystems of Algebraic Equations and Systems of Differential Equations
Systems of Algebraic Equations and Systems of Differential Equations Topics: 2 by 2 systems of linear equations Matrix expression; Ax = b Solving 2 by 2 homogeneous systems Functions defined on matrices
More informationA NONLINEARITY MEASURE FOR ESTIMATION SYSTEMS
AAS 6-135 A NONLINEARITY MEASURE FOR ESTIMATION SYSTEMS Andrew J. Sinclair,JohnE.Hurtado, and John L. Junkins The concept of nonlinearity measures for dynamical systems is extended to estimation systems,
More informationM3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011
M3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011 Text for the course: Professor Darryl D Holm 25 October 2011 Imperial College London d.holm@ic.ac.uk http://www.ma.ic.ac.uk/~dholm/ Geometric Mechanics
More informationJohn E Hurtado Associate Professor of Aerospace Engineering, Texas A&M University, College Station, TX
John E Hurtado Associate Professor of Aerospace Engineering, Texas A&M University, College Station, TX http://aero.tamu.edu/faculty/hurtado/ Education Expertise Mentorship Teaching Awards Impact Scholarship
More informationarxiv: v4 [quant-ph] 9 Jun 2016
Applying Classical Geometry Intuition to Quantum arxiv:101.030v4 [quant-ph] 9 Jun 016 Spin Dallin S. Durfee and James L. Archibald Department of Physics and Astronomy, Brigham Young University, Provo,
More informationOPTIMAL SPACECRAF1 ROTATIONAL MANEUVERS
STUDIES IN ASTRONAUTICS 3 OPTIMAL SPACECRAF1 ROTATIONAL MANEUVERS JOHNL.JUNKINS Texas A&M University, College Station, Texas, U.S.A. and JAMES D.TURNER Cambridge Research, Division of PRA, Inc., Cambridge,
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling
1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples
More informationMechanics. In the Science Program, Mechanics contributes to the following program goals described in the Exit Profile:
Mechanics Objectives: 00UR Discipline: Physics Ponderation: 3-2-3 Course Code: 203-NYA-05 Prerequisite: Sec. V Physics 534, Mathematics 536 (or equivalent) Course Credit: 2 2/3 Corequisite: 00UP (Calculus
More informationThe Liapunov Method for Determining Stability (DRAFT)
44 The Liapunov Method for Determining Stability (DRAFT) 44.1 The Liapunov Method, Naively Developed In the last chapter, we discussed describing trajectories of a 2 2 autonomous system x = F(x) as level
More informationPhysics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top
Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem
More information9.1 Eigenanalysis I Eigenanalysis II Advanced Topics in Linear Algebra Kepler s laws
Chapter 9 Eigenanalysis Contents 9. Eigenanalysis I.................. 49 9.2 Eigenanalysis II................. 5 9.3 Advanced Topics in Linear Algebra..... 522 9.4 Kepler s laws................... 537
More informationROTATIONAL STABILITY OF A CHARGED DIELEC- TRIC RIGID BODY IN A UNIFORM MAGNETIC FIELD
Progress In Electromagnetics Research Letters, Vol. 11, 103 11, 009 ROTATIONAL STABILITY OF A CHARGED DIELEC- TRIC RIGID BODY IN A UNIFORM MAGNETIC FIELD G.-Q. Zhou Department of Physics Wuhan University
More informationKinematics. Chapter Multi-Body Systems
Chapter 2 Kinematics This chapter first introduces multi-body systems in conceptual terms. It then describes the concept of a Euclidean frame in the material world, following the concept of a Euclidean
More informationGetting Started with Communications Engineering
1 Linear algebra is the algebra of linear equations: the term linear being used in the same sense as in linear functions, such as: which is the equation of a straight line. y ax c (0.1) Of course, if we
More informationDynamical Systems. August 13, 2013
Dynamical Systems Joshua Wilde, revised by Isabel Tecu, Takeshi Suzuki and María José Boccardi August 13, 2013 Dynamical Systems are systems, described by one or more equations, that evolve over time.
More informationLinear Systems of ODE: Nullclines, Eigenvector lines and trajectories
Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 6, 203 Outline
More informationMath 396. Quotient spaces
Math 396. Quotient spaces. Definition Let F be a field, V a vector space over F and W V a subspace of V. For v, v V, we say that v v mod W if and only if v v W. One can readily verify that with this definition
More informationPhysics 221A Fall 2017 Notes 27 The Variational Method
Copyright c 2018 by Robert G. Littlejohn Physics 221A Fall 2017 Notes 27 The Variational Method 1. Introduction Very few realistic problems in quantum mechanics are exactly solvable, so approximation methods
More informationAn Exactly Solvable 3 Body Problem
An Exactly Solvable 3 Body Problem The most famous n-body problem is one where particles interact by an inverse square-law force. However, there is a class of exactly solvable n-body problems in which
More informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More informationMATHEMATICAL PHYSICS
MATHEMATICAL PHYSICS Third Year SEMESTER 1 015 016 Classical Mechanics MP350 Prof. S. J. Hands, Prof. D. M. Heffernan, Dr. J.-I. Skullerud and Dr. M. Fremling Time allowed: 1 1 hours Answer two questions
More informationLinear Systems of ODE: Nullclines, Eigenvector lines and trajectories
Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 6, 2013 Outline
More informationThe quantum state as a vector
The quantum state as a vector February 6, 27 Wave mechanics In our review of the development of wave mechanics, we have established several basic properties of the quantum description of nature:. A particle
More informationMATH 415, WEEKS 14 & 15: 1 Recurrence Relations / Difference Equations
MATH 415, WEEKS 14 & 15: Recurrence Relations / Difference Equations 1 Recurrence Relations / Difference Equations In many applications, the systems are updated in discrete jumps rather than continuous
More informationProblem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension
105A Practice Final Solutions March 13, 01 William Kelly Problem 1: Lagrangians and Conserved Quantities Consider the following action for a particle of mass m moving in one dimension S = dtl = mc dt 1
More informationAN ALTERNATING MINIMIZATION ALGORITHM FOR NON-NEGATIVE MATRIX APPROXIMATION
AN ALTERNATING MINIMIZATION ALGORITHM FOR NON-NEGATIVE MATRIX APPROXIMATION JOEL A. TROPP Abstract. Matrix approximation problems with non-negativity constraints arise during the analysis of high-dimensional
More informationarxiv: v1 [hep-ph] 5 Sep 2017
A First Step Towards Effectively Nonperturbative Scattering Amplitudes in the Perturbative Regime Neil Christensen, Joshua Henderson, Santiago Pinto, and Cory Russ Department of Physics, Illinois State
More informationMaps and differential equations
Maps and differential equations Marc R. Roussel November 8, 2005 Maps are algebraic rules for computing the next state of dynamical systems in discrete time. Differential equations and maps have a number
More informationWEIGHTING MATRICES DETERMINATION USING POLE PLACEMENT FOR TRACKING MANEUVERS
U.P.B. Sci. Bull., Series D, Vol. 75, Iss. 2, 2013 ISSN 1454-2358 WEIGHTING MATRICES DETERMINATION USING POLE PLACEMENT FOR TRACKING MANEUVERS Raluca M. STEFANESCU 1, Claudiu L. PRIOROC 2, Adrian M. STOICA
More informationAdvanced Engineering Mathematics Prof. Pratima Panigrahi Department of Mathematics Indian Institute of Technology, Kharagpur
Advanced Engineering Mathematics Prof. Pratima Panigrahi Department of Mathematics Indian Institute of Technology, Kharagpur Lecture No. #07 Jordan Canonical Form Cayley Hamilton Theorem (Refer Slide Time:
More information13 Path Planning Cubic Path P 2 P 1. θ 2
13 Path Planning Path planning includes three tasks: 1 Defining a geometric curve for the end-effector between two points. 2 Defining a rotational motion between two orientations. 3 Defining a time function
More informationSTEP Support Programme. STEP 2 Matrices Topic Notes
STEP Support Programme STEP 2 Matrices Topic Notes Definitions............................................. 2 Manipulating Matrices...................................... 3 Transformations.........................................
More informationState observers for invariant dynamics on a Lie group
State observers for invariant dynamics on a Lie group C. Lageman, R. Mahony, J. Trumpf 1 Introduction This paper concerns the design of full state observers for state space systems where the state is evolving
More informationAcceleration of Levenberg-Marquardt method training of chaotic systems fuzzy modeling
ISSN 746-7233, England, UK World Journal of Modelling and Simulation Vol. 3 (2007) No. 4, pp. 289-298 Acceleration of Levenberg-Marquardt method training of chaotic systems fuzzy modeling Yuhui Wang, Qingxian
More informationarxiv: v1 [math.ds] 18 Nov 2008
arxiv:0811.2889v1 [math.ds] 18 Nov 2008 Abstract Quaternions And Dynamics Basile Graf basile.graf@epfl.ch February, 2007 We give a simple and self contained introduction to quaternions and their practical
More informationSimplest Chaotic Flows with Involutional Symmetries
International Journal of Bifurcation and Chaos, Vol. 24, No. 1 (2014) 1450009 (9 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127414500096 Simplest Chaotic Flows with Involutional Symmetries
More informationJordan Canonical Form Homework Solutions
Jordan Canonical Form Homework Solutions For each of the following, put the matrix in Jordan canonical form and find the matrix S such that S AS = J. [ ]. A = A λi = λ λ = ( λ) = λ λ = λ =, Since we have
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationSection 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d
Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d July 6, 2009 Today s Session Today s Session A Summary of This Session: Today s Session A Summary of This
More informationA plane autonomous system is a pair of simultaneous first-order differential equations,
Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium
More informationarxiv: v2 [nlin.cd] 8 Sep 2012
An analytical limitation for time-delayed feedback control in autonomous systems Edward W. Hooton 1 1, 2, and Andreas Amann 1 School of Mathematical Sciences, University College Cork, Ireland 2 Tyndall
More informationOne Dimensional Dynamical Systems
16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with one-dimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar
More informationJim Lambers MAT 610 Summer Session Lecture 2 Notes
Jim Lambers MAT 610 Summer Session 2009-10 Lecture 2 Notes These notes correspond to Sections 2.2-2.4 in the text. Vector Norms Given vectors x and y of length one, which are simply scalars x and y, the
More informationPhysics 209 Fall 2002 Notes 5 Thomas Precession
Physics 209 Fall 2002 Notes 5 Thomas Precession Jackson s discussion of Thomas precession is based on Thomas s original treatment, and on the later paper by Bargmann, Michel, and Telegdi. The alternative
More informationVideo 8.1 Vijay Kumar. Property of University of Pennsylvania, Vijay Kumar
Video 8.1 Vijay Kumar 1 Definitions State State equations Equilibrium 2 Stability Stable Unstable Neutrally (Critically) Stable 3 Stability Translate the origin to x e x(t) =0 is stable (Lyapunov stable)
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationBasic Theory of Dynamical Systems
1 Basic Theory of Dynamical Systems Page 1 1.1 Introduction and Basic Examples Dynamical systems is concerned with both quantitative and qualitative properties of evolution equations, which are often ordinary
More informationWhat s Eigenanalysis? Matrix eigenanalysis is a computational theory for the matrix equation
Eigenanalysis What s Eigenanalysis? Fourier s Eigenanalysis Model is a Replacement Process Powers and Fourier s Model Differential Equations and Fourier s Model Fourier s Model Illustrated What is Eigenanalysis?
More informationTheory of Vibrations in Stewart Platforms
Theory of Vibrations in Stewart Platforms J.M. Selig and X. Ding School of Computing, Info. Sys. & Maths. South Bank University London SE1 0AA, U.K. (seligjm@sbu.ac.uk) Abstract This article develops a
More informationCOSMOS: Making Robots and Making Robots Intelligent Lecture 3: Introduction to discrete-time dynamics
COSMOS: Making Robots and Making Robots Intelligent Lecture 3: Introduction to discrete-time dynamics Jorge Cortés and William B. Dunbar June 3, 25 Abstract In this and the coming lecture, we will introduce
More information7.6 Additional Problems
7in x 10in Felder c07_online.tex V3 - January 21, 2015 10:11 A.M. Page 26 26 Chapter 7 Linear Algebra II (nline) 7.6 Additional Problems 7.137 The picture below shows a rubber sheet with three small stones
More informationQuaternion-Based Tracking Control Law Design For Tracking Mode
A. M. Elbeltagy Egyptian Armed forces Conference on small satellites. 2016 Logan, Utah, USA Paper objectives Introduction Presentation Agenda Spacecraft combined nonlinear model Proposed RW nonlinear attitude
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationColumbia University Department of Physics QUALIFYING EXAMINATION
Columbia University Department of Physics QUALIFYING EXAMINATION Wednesday, January 13, 2016 3:10PM to 5:10PM Modern Physics Section 4. Relativity and Applied Quantum Mechanics Two hours are permitted
More informationApplications of Diagonalization
Applications of Diagonalization Hsiu-Hau Lin hsiuhau@phys.nthu.edu.tw Apr 2, 200 The notes cover applications of matrix diagonalization Boas 3.2. Quadratic curves Consider the quadratic curve, 5x 2 4xy
More informationModal Decomposition and the Time-Domain Response of Linear Systems 1
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING.151 Advanced System Dynamics and Control Modal Decomposition and the Time-Domain Response of Linear Systems 1 In a previous handout
More informationNOTES ON LINEAR ALGEBRA CLASS HANDOUT
NOTES ON LINEAR ALGEBRA CLASS HANDOUT ANTHONY S. MAIDA CONTENTS 1. Introduction 2 2. Basis Vectors 2 3. Linear Transformations 2 3.1. Example: Rotation Transformation 3 4. Matrix Multiplication and Function
More informationCOMPLETE ALL ROUGH WORKINGS IN THE ANSWER BOOK AND CROSS THROUGH ANY WORK WHICH IS NOT TO BE ASSESSED.
BSc/MSci EXAMINATION PHY-304 Time Allowed: Physical Dynamics 2 hours 30 minutes Date: 28 th May 2009 Time: 10:00 Instructions: Answer ALL questions in section A. Answer ONLY TWO questions from section
More informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
More information