Keywords: Kinematics; Principal planes; Special orthogonal matrices; Trajectory optimization
|
|
- Alan McCormick
- 5 years ago
- Views:
Transcription
1 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: Kinematics; Principal planes; Special orthogonal matrices; Trajectory optimization Corresponding Author: Dr. Andrew Sinclair, Corresponding Author's Institution: Auburn University First Author: Andrew Sinclair Order of Authors: Andrew Sinclair; John E Hurtado; John L Junkins Abstract: The group of special (or proper) orthogonal matrices, SO(N), is used throughout engineering mechanics in the analysis and representation of mechanical systems. In this paper, a solution is presented for the optimal transformation between two elements of SO(N). The transformation is assumed to occur during a specified finite time, and a cost function that penalizes the transformation rates is utilized. The optimal transformation is found as a constant-rate rotation in each of the principal planes relating the two elements. Although the kinematics of SO(N) are nonlinear and governed by Poisson's equation, the solution is found to be a linear function of the generalized principal angles. This is made possible by the extension of principal-rotation kinematics from three-dimensional rotations to the general SO(N) group. This extension relates the N-dimensional angular velocity to the derivatives of the principal angles. The cost of the optimal transformation, the square root of the sum of the principal angles squared, also provides a useful measure for the angular distance between two elements of SO(N). Suggested Reviewers: Shane Ross Assistant Professor, Engineering Science and Mechanics, Virginia Polytechnic Institute and State University sdross@vt.edu
2 Craig Woolsey Associate Professor, Aerospace & Ocean Engineering, Virginia Polytechnic Institute and State University Maruthi Akella Associate Professor, Aerospace Engineering & Engineering Mechanics, University of Texas at Austin
3 Manuscript Click here to download Manuscript: SHJ_NonlinDyn.ps Minimum Distance and Optimal Transformations on SO(N) ANDREW J. SINCLAIR 1,, JOHN E. HURTADO 2, and JOHN L. JUNKINS 2 1 Department of Aerospace Engineering, Auburn University, 211 Aerospace Engineering Building, Auburn, AL , USA; 2 Department of Aerospace Engineering, Texas A&M University, TAMU-3141, College Station, TX , USA; Author for correspondence ( sinclair@auburn.edu; fax: ; phone: ) Abstract The group of special (or proper) orthogonal matrices, SO(N), is used throughout engineering mechanics in the analysis and representation of mechanical systems. In this paper, a solution is presented for the optimal transformation between two elements of SO(N). The transformation is assumed to occur during a specified finite time, and a cost function that penalizes the transformation rates is utilized. The optimal transformation is found as a constant-rate rotation in each of the principal planes relating the two elements. Although the kinematics of SO(N) are nonlinear and governed by Poisson s equation, the solution is found to be a linear function of the generalized principal angles. This is made possible by the extension of principal-rotation kinematics from three-dimensional rotations to the general SO(N) group. This extension relates the N-dimensional angular velocity to the derivatives of the principal angles. The cost of the optimal transformation, the square root of the sum of the principal angles squared, 1
4 also provides a useful measure for the angular distance between two elements of SO(N). Keywords: Kinematics, Principal planes, Special orthogonal matrices, Trajectory optimization 1 Introduction Special (or proper) orthogonal matrices, C SO(N), arise frequently in the analysis and representation of dynamic systems. Most commonly the group SO(3) is used to represent the orientation of rigid bodies. By analogy, however, the group SO(N) can be used to describe the configuration of M-degree of freedom systems, where M = N(N 1)/2 [11, 13]. Additionally, the group SO(N) is used frequently in similarity transformations and other decompositions of fundamental quantities related to engineering mechanics. The variety of uses for these matrices motivates the investigation of the following two problems. First, what is the angular distance between two elements of SO(N)? Second, what is the best transition between two orientations belonging to SO(N)? In this paper, these questions will be investigated as the solution to a nonlinear trajectoryoptimization problem subject to the kinematic behavior of the SO(N) group. Of course, many properties of the SO(N) group are generalizations of properties of threedimensional rotations, and it will be shown that using descriptions of principal rotations leads to a solution that is linear in the principal angles. For three-dimensional rotations, Euler s theorem describes any general orientation 2
5 in terms of a single principal rotation. The principal rotation concept also extends to N-dimensional rotations [1, 2, 7, 9, 10]; however, for higher-dimensional spaces a general orientation requires N/2 principal rotations for even N and (N 1)/2 principal rotations for odd N. In general the number of required rotations can be expressed as L = N/2. These rotations take place on completely orthogonal planes, called the principal planes. For even dimensions these planes completely occupy the space. For odd dimensions, however, one axis is left out of the rotational motion and is referred to as a principal axis. The following section reviews some of the methods that have been developed to describe SO(N) in terms of principal rotations. This view of principal rotations is then extended to develop new results for the kinematic behavior of the rotational group. Using these results, a solution is found for the optimal transition between two elements of SO(N). The resulting optimal cost function is a convenient definition of the angular distance between two elements of SO(N). The two main contributions of the paper are the principal-rotation kinematics and the measure of distance between elements of SO(N). As will be described, the canonical form decomposes the various representation matrices into a block-diagonal form with a 2 2 block associated with each principal plane and a 1 1 block associated with the principal axis if it exists. In illustrating these blocks, a subscript k will be used to indicate the kth 2 2 block on the diagonal of the indicated matrix. 3
6 2 Review of N-dimensional rotations Much of the description of N-dimensional orientations in this section was given by Mortari [7] and Bauer [2]. The current work attempts to follow their notation and conventions as closely as possible with one exception. In discussing a rotated frame both of the above authors define representations of the rotation as the mapping from the rotated frame back to a reference frame. Here the convention will be to give the mapping from the reference frame to the rotated frame. Therefore, many of the definitions given below correspond to the transpose of the rotation matrices given by Mortari and Bauer. 2.1 Rotation matrix The transformation of an N-dimensional vector by a proper orthogonal matrix, C, describes a rotation in N-dimensional space [3]. The following equation describes the transformation from a column matrix parameterizing a vector in a reference coordinate system, the n frame, to a column matrix parameterizing the vector in a rotated coordinate system, the b frame. [r] b = [C] [r] n (1) This matrix C is called a rotation matrix and is the most fundamental representation of N-dimensional rotations. An important similarity transformation of the rotation matrix is the block-diagonal skew-symmetric canonical form for normal matrices [5], referred to here as simply the 4
7 canonical representation. C = P T C P; C = PCP T (2) Here, P is a proper orthogonal matrix, and C is a block-diagonal proper orthogonal matrix. The rows of P are the coordinatization of the principal coordinate vectors in the b frame. [ [P] = [p 1 ] b [p 2 ] b... [p N ] b ] T (3) The matrix C is related to the principal angles. The kth block on the diagonal of C has the following form. C k = cos (φ k + 2πn k ) sin (φ k + 2πn k ) sin (φ k + 2πn k ) cos (φ k + 2πn k ) (4) Here, each angle π φ k π is the value of rotation in the kth principal plane, and the values n k can be any integer. For odd N, the (N, N) element of C forms a 1 1 block and is equal to positive one. The matrix P is itself an N-dimensional rotation matrix that describes the transformation from the rotated frame to a third frame, the principal frame. This coordinate system has coordinate vectors {p 1, p 2,...,p N } which are aligned with the principal planes of the rotation described by C: (p 1, p 2 ), (p 3, p 4 ), etc. The similarity transformation in the second of Eqs. (2) can be described as coordinatizing or viewing C in this principal frame. Note that consistent with his convention mentioned earlier, Bauer defines P as the transpose of the definition given here; thus it is the mapping from the principal to rotated frame. 5
8 2.2 Euler matrix Another representation of N-dimensional rotations that will be useful for the current purposes is the Euler matrix. Whereas this matrix has been used tangentially in previous works [1,7], it will be developed more fully here. The Euler matrix is a skewsymmetric matrix, E, that can be related to the rotation matrix using properties of the matrix exponential and determinant [4]. exp (E) (exp (E)) T = exp (E) exp ( E T) = exp ( E + E T) = exp (0) = I (5) det (exp (E)) = exp (Tr (E)) = exp (0) = +1 (6) Because exp (E) is proper orthogonal, the following relationship to the rotation matrix can be considered the definition of the Euler matrix. C = exp (E) (7) Based on this definition it is possible to relate the Euler matrix to the principal rotations and solve for E. Although it is tempting to simply write E = ln(c), it will be shown that this mapping is not unique because infinitely many solutions for E correspond to any particular orientation. Additionally, the matrix logarithm can suffer from limited range of convergence [4]. Because E is also normal, it can be expressed in the canonical representation. E = P T E P; E = PEP T (8) The first of Eqs. (8) can be substituted into Eq. (7). C = P T exp (E ) P (9) 6
9 Next, consider the following form for the blocks on the diagonal of E. E k = 0 φ k + 2πn k φ k 2πn k 0 = (φ k + 2πn k )J (10) The exponential of this form is the canonical form of the rotation matrix. exp (E k) = sin (φ k + 2πn k )J + cos (φ k + 2πn k ) I = C k (11) Therefore, the canonical form of E is made of blocks as shown in Eq. (10), the matrix P is the canonical transformation for both C and E, and C = exp(e ). The matrix E is shown below for odd N. [E ] = 0 φ 1 + 2πn φ 1 2πn φ L + 2πn L φ L 2πn L (12) For even N, the form is identical except the omission of the row and column of zeros. Again, the sign convention above is chosen to be consistent with C in Eq. (4). Because of the ambiguity in the integers n, infinitely many values of E and E exist which correspond to any particular C. 3 Kinematics of principal rotations In the previous section several methods to describe N-dimensional rotations in terms of the principal planes and angles were reviewed. In this section new results are presented 7
10 that relate the N-dimensional angular velocity and the derivatives of the principal planes and angles. This results in kinematic differential equations for P and φ k [12], but here the focus will be on φ k. Traditionally, the kinematic evolution of N-dimensional rotations has not been related to the principal rotations. Instead, Poisson s equation for the derivative of C is used directly. Ċ = ΩC ; Ω = CC T (13) Here, Ω is the N-dimensional skew-symmetric angular-velocity matrix. Whereas Poisson s equation holds for any value of N, for three-dimensional rotations there also exists relationships between the angular velocity and the derivatives of the principal angle, φ, and principal axis, â [6]. Following an analogous approach, the derivatives of the principal angles and principal planes of higher dimensions can also be related to the angular velocity [12]. For the current purposes, it will be sufficient to examine the principal angles. The angular velocity and the principal-angle derivatives can be related using the canonical form of C. The derivative of the canonical form is taken as follows. Ċ = P T C P + P T d dt (C ) P + P T C P (14) As mentioned, P itself is a rotation matrix describing the transformation from the rotated to principal frame. Its derivative is related to the angular velocity of the principal frame relative to the rotated frame. This skew-symmetric matrix is defined as Ψ = PP T, and P satisfies the Poisson equation P = ΨP. This is used to 8
11 rewrite the derivative of C. ( ) d Ċ = P T dt (C ) + ΨC C Ψ P (15) The angular velocity Ω can now be written in terms of the canonical form C and its derivative d dt (C ). ( ) Ω = CC d T = P T dt (C ) + ΨC C Ψ PP T C T P = P T ( d dt (C ) C T + Ψ C ΨC T ) P (16) The blocks on the diagonal of C have the form shown in Eq. (4), and the derivative of C will of course also be block diagonal with blocks of the following form. [ ] d dt (C ) = φ sin (φ k + 2πn k ) cos (φ k + 2πn k ) k (17) k cos (φ k + 2πn k ) sin (φ k + 2πn k ) For odd N, the (N, N) element of d dt (C ) forms a 1 1 block and is equal to zero. Because both C and d dt (C ) are block diagonal, the product d dt (C ) C T will also be block diagonal with blocks of the following form. [ ] d dt (C ) C T k = [ ] d dt (C ) k [ ] C T = k 0 φk φ k 0 (18) Again, for odd N the (N, N) element of the product is zero. The remaining terms in parentheses on the right-hand side of Eq. (16) are skew symmetric but in general not block diagonal. To develop the relationship between the angular velocity and the principal angle derivatives, however, the blocks on the diagonal of these terms will be investigated. Because C is block diagonal, the third term has the following form where 9
12 the 2πn k terms in C have been dropped for convenience. [ C ΨC T] = k C kψ k C k T cos (φ k ) sin (φ k ) 0 Ψ 2k 1,2k cos (φ k ) sin (φ k ) = sin (φ k ) cos (φ k ) Ψ 2k 1,2k 0 sin (φ k ) cos (φ k ) = 0 Ψ 2k 1,2k Ψ 2k 1,2k 0 (19) The blocks on the diagonal of C ΨC T are identical to the blocks on the diagonal of Ψ. Therefore, the only contribution to the blocks on the diagonal of Eq. (16) comes from Eq. (18). [ PΩP T ] k = 0 φk φ k 0 (20) Thus, the blocks on the diagonal of the angular velocity viewed in the principal frame are simply related to the derivative of the principal angles. This relation is a new result and is a generalization of the three-dimensional concept φ = a T ω [6]. The minus signs in these components are artifacts of the convention chosen to define the angular-velocity matrix as Ω = CC T. These minus signs are entirely equivalent to the convention in defining the (1, 2) component of the three-dimensional angular-velocity matrix as ω 3. Of course, this choice is made to make angular-velocity matrix multiplication equivalent to the angular-velocity vector cross product. These conventions are maintained for N-dimensions even though the angular-velocity vector and cross product lose their physical significance. In the following section, this result for the components of angular velocity lying in the principal planes will be used to study the optimal transition between two elements of SO(N). 10
13 4 Optimal kinematic maneuvers Using the results of the previous two sections, the optimal transformation between two elements of SO(N) is now investigated. It is well known that the minimum angular distance between two orientations in three dimensions is the principal angle associated with the rotation matrix relating them. A rigid-body rotational maneuver about the corresponding Euler axis through the principal angle is called an eigenaxis rotation. Clearly, a related solution is anticipated in the N-dimensional case. Of course, infinitely many transformations can relate two boundary conditions in SO(N). To obtain a single transformation, a cost function can be defined to determine an optimal solution. A cost function that penalizes the transformational rates used in the maneuver is considered. Ċdt = C (t + dt) C (t) dc (21) J 1 = T 0 dc = T 0 Ċ dt (22) Here, indicates the Frobenius matrix norm: A tr (AA T ). This function can be used to evaluate the cost or distance of a particular transitional path between two orientations. In Eq. (22), the initial time is zero, and the final time is T. Without loss of generality, the boundary conditions are chosen such that the rotated frame is initially aligned with the reference frame, C (0) = I, and an arbitrary final orientation is chosen, C (T) = F. The matrix Ċ obeys Poisson s equation. Ċ = ΩC (23) 11
14 This kinematic equation can be used in Eq. (22) to obtain the following. J 1 = T 0 ΩC dt = T 0 tr(ωωt )dt (24) The minimum distance between the two orientations is now given as the minimization of Eq. (24) subject to the boundary conditions and kinematic equations, Eq. (23). It is convenient to disregard the square root function in Eq. (24) and investigate a slightly modified problem. J 2 = 1 2 T 0 tr ( ΩΩ T) dt (25) The minimization of Eq. (25) is now sought subject to the kinematic equations and boundary conditions. Because Eqs. (24) and (25) are related via a monotone transformation of the integrand, they share the same minimizer. The Hamiltonian for the problem is easily found by introducing a matrix of costates, λ. H = 1 2 tr( ΩΩ T) tr ( λ T ΩC ) (26) From the Hamiltonian, the first-order necessary conditions are computed using properties of the trace. Ċ = H λ = ΩC (27) λ = H C = ΩT λ (28) 0 = H Ω = 1 2 Ω Ω λct (29) The third condition gives the following. Ω = λc T (30) 12
15 Next, consider the derivative of Eq. (30) together with the first two conditions, Eqs. (27) and (28). Ω = λc T + λċt = Ω T λc T λc T Ω T (31) Equation (30) itself can now be used, as well as the skew-symmetry of Ω. Ω = Ω T Ω ΩΩ T = Ω T Ω Ω T Ω = 0 (32) Therefore, the optimal angular velocity is a constant. Equations (26 32) demonstrate the optimality of constant angular velocity, relying only on kinematic principles. Alternative demonstrations focus on kinetic principles, like Hamilton s principle of motion. Equation (25) can be recognized as the action integral for an N-dimensional rigid body with identity mass matrix. According to Hamilton s principle, the equations of motion, given by Ratiu [8], give a stationary value for the action integral. For identity mass matrix, these equations simplify to Ω = 0. For constant angular velocity, the Poisson kinematic equations have the following solution. C (t) = exp ( Ωt) C (0) (33) The value of the optimal angular velocity can be related to the boundary conditions of the problem. F = exp ( ΩT) (34) The implication of this is that along the optimal solution, the matrix ΩT is an Euler matrix of the final orientation. It was observed earlier that infinitely many Euler 13
16 matrices exist for any particular orientation. The optimal angular velocity can be written using the canonical form of the Euler matrix, which was presented in Eqs. (8) and (12). For odd N, this is shown below. [Ω] = 1 T [P]T 0 φ 1 + 2πn φ 1 2πn [P] (35) φ L + 2πn L φ L 2πn L For even N, the form is identical omitting the row and columns of zeros. Therefore, infinitely many angular velocities satisfy the first-order optimality conditions, and it remains to be shown that selecting n k = 0 gives the global minimum of the cost function. Because the optimal angular velocity has been found to be a constant, the cost function evaluated for the optimal solution can be rewritten as follows. J 2 = T 2 tr( Ω T Ω ) (36) In the previous section it was observed that applying the canonical transformation P of the rotation matrix to the angular velocity does not in general produce a block-diagonal form. Because the optimal angular velocity is proportional to the Euler matrix of F, however, Eq. (35) shows that in this case PΩP T will be block diagonal. This will be defined as Ω. J 2 = T 2 tr( P T Ω T Ω P ) = T 2 tr ( Ω T Ω ) (37) 14
17 The product Ω T Ω, however, is a diagonal matrix, and its trace can be used to write the cost function as shown below. J 2 = 1 2T [ (φ1 + 2πn 1 ) (φ L + 2πn L ) 2] (38) Therefore, finding the optimal-angular velocity is reduced to minimizing each of the parenthetical terms above. This is clearly done by selecting each n k equal to zero. This gives the final result for the optimal angular velocity. [Ω] = 1 T [P]T 0 φ φ [P] (39) φ L φ L Considering the results from the previous section, the optimal maneuver is a rotation in each of the principal planes relating the initial and final orientations with a rotational rate of φ k /T. It is the coordinatization in the principal planes that allows the solution to the nonlinear trajectory-optimization problem to be expressed in this form, which is linear in the principal angles. As mentioned, this solution for the minimization of J 2 must also minimize J 1. The optimal cost for J 1 can be evaluated using the found angular velocity. J 1 = T tr (ΩΩ T ) = φ φ 2 L (40) This demonstrates that the minimum angular distance between two N-dimensional orientations is indeed the square root of the sum of the squares of the principal angles. 15
18 5 Conclusion In this paper, the kinematic interpretation of principal rotations was extended to higher dimensional elements of SO(N). This result led to a solution for the optimal transition between two elements of this group. The optimal cost of the resulting solution provides a convenient, intuitive measure for the angular distance between such matrices. It is anticipated that these results will prove useful in a wide variety of mechanics fields where elements of SO(N) are used in the analysis and comparison of engineering systems. Critically, the posed trajectory-optimization problem was subject to the nonlinear kinematics of the SO(N) group described by Poisson s equation. Judicious selection for the representation and coordinatization of the control variables, however, resulted in a physically meaningful solution that is linear in the orientation variables. This was made possible through the use of the Euler matrix, canonical form, and principal plane concepts. These transformations allowed the presentation of a solution of linear form to the nonlinear problem. References [1] Itzhack Y. Bar-Itzhack. Extension of Euler s theorem to n-dimensional spaces. IEEE Transactions on Aerospace and Electronic Systems, 25(6): , [2] Robert Bauer. Euler s theorem on rigid body displacements generalized to n dimensions. The Journal of the Astronautical Sciences, 50(3): ,
19 [3] O. Bottema and B. Roth. Theoretical Kinematics. North-Holland Publishing Company, Amsterdam, Chapters 1 and 2. [4] Morton L. Curtis. Matrix Groups. Springer-Verlag, New York, [5] Roger A. Horn and Charles R. Johnson. Matrix Analysis. Cambridge University Press, pp [6] Peter C. Hughes. Spacecraft Attitude Dynamics. John Wiley & Sons, New York, pp [7] Daniele Mortari. On the rigid rotation concept in n-dimensional spaces. Journal of the Astronautical Sciences, 49(3): , [8] T. Ratiu. The motion of the free n-dimensional rigid body. Indiana University Mathematics Journal, 29(4): , [9] Hanspeter Schaub, Panagiotis Tsiotras, and John L. Junkins. Principal rotation representations of proper n n orthogonal matrices. International Journal of Engineering Science, 33(15): , [10] P. H. Schoute. Le déplacement le plus général dans l espace à n dimensions. Annales de l Ecole Polytechnique de Delft, 7: , [11] A. J. Sinclair and J. E. Hurtado. Cayley kinematics and the Cayley form of dynamic equations. Proceedings of the Royal Society of London Series A, 461(2055): ,
20 [12] A. J. Sinclair, J. E. Hurtado, and J. L. Junkins. Kinematics of N-dimensional principal rotations. In AAS/AIAA Spaceflight Mechanics Meeting, Copper Mountain, Colorado, January Paper AAS [13] A. J. Sinclair, J. E. Hurtado, and J. L. Junkins. Application of the Cayley form to general spacecraft motion. Journal of Guidance, Control, and Dynamics, 29(2): ,
Minimum-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 informationLinear 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 informationApplication of the Cayley Form to General Spacecraft Motion
JOURNAL OF GUIDANCE, CONTROL, AND DYNAMICS Vol. 29, No. 2, March April 2006 Application of the Cayley Form to General Spacecraft Motion Andrew J. Sinclair, John E. Hurtado, and John L. Junkins Texas A&M
More information9 th AAS/AIAA Astrodynamics Specialist Conference Breckenridge, CO February 7 10, 1999
AAS 99-139 American Institute of Aeronautics and Astronautics MATLAB TOOLBOX FOR RIGID BODY KINEMATICS Hanspeter Schaub and John L. Junkins 9 th AAS/AIAA Astrodynamics Specialist Conference Breckenridge,
More informationApplication of the Cayley Form to General. Spacecraft Motion
Application of the Cayley Form to General Spacecraft Motion Andrew J. Sinclair, John E. Hurtado, and John L. Junkins Texas A&M University, College Station, Texas 77843 Abstract The study of N-dimensional
More informationUSING CARLEMAN EMBEDDING TO DISCOVER A SYSTEM S MOTION CONSTANTS
(Preprint) AAS 12-629 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
More informationNew Results for Time-Optimal Three-Axis. Reorientation of a Rigid Spacecraft
New Results for Time-Optimal Three-Axis Reorientation of a Rigid Spacecraft Xiaoli Bai, John L. Junkins Texas A&M University, College Station, Texas 77843 Abstract Some new results for a classical minimum,
More informationTracking Rigid Body Motion Using Thrusters and Momentum. Wheels
JAS 199 Tracking Rigid Body Motion Using Thrusters and Momentum Wheels Christopher D. Hall, Panagiotis Tsiotras, and Haijun Shen Abstract Tracking control laws are developed for a rigid spacecraft using
More informationPrincipal Rotation Representations of Proper NxN Orthogonal Matrices
Principal Rotation Representations of Proper NxN Orthogonal Matrices Hanspeter Schaub, Panagiotis siotras and John L. Junkins Simulated Reprint from International Journal of Engineering Science Vol., No.,
More informationAnalytical Mechanics. of Space Systems. tfa AA. Hanspeter Schaub. College Station, Texas. University of Colorado Boulder, Colorado.
Analytical Mechanics of Space Systems Third Edition Hanspeter Schaub University of Colorado Boulder, Colorado John L. Junkins Texas A&M University College Station, Texas AIM EDUCATION SERIES Joseph A.
More information*School of Aeronautics and Astronautics, Purdue University, West Lafayette, IN
NEW CONTROL LAWS FOR THE ATTITUDE STABILIZATION OF RIGID BODIES PANAGIOTIS TSIOTRAS *School of Aeronautics and Astronautics, Purdue University, West Lafayette, IN 7907. Abstract. This paper introduces
More information1 - z w = 1 + z. 1 - w z = 1 + w S 1
Journal of Guidance, Control, and Dynamics Vol. 0, No., May-June 997 Higher Order Cayley Transforms with Applications to Attitude Representations Panagiotis Tsiotras University of Virginia, Charlottesville,
More informationMixed Control Moment Gyro and Momentum Wheel Attitude Control Strategies
AAS03-558 Mixed Control Moment Gyro and Momentum Wheel Attitude Control Strategies C. Eugene Skelton II and Christopher D. Hall Department of Aerospace & Ocean Engineering Virginia Polytechnic Institute
More informationLecture Notes - Modeling of Mechanical Systems
Thomas Bak Lecture Notes - Modeling of Mechanical Systems February 19, 2002 Aalborg University Department of Control Engineering Fredrik Bajers Vej 7C DK-9220 Aalborg Denmark 2 Table of Contents Table
More informationOrtho-Skew and Ortho-Sym Matrix Trigonometry 1
Abstract Ortho-Skew and Ortho-Sym Matrix Trigonometry Daniele Mortari Department of Aerospace Engineering, Texas A&M University, College Station, TX 77843-34 This paper introduces some properties of two
More informationAngular Velocity Determination Directly from Star Tracker Measurements
Angular Velocity Determination Directly from Star Tracker Measurements John L. Crassidis Introduction Star trackers are increasingly used on modern day spacecraft. With the rapid advancement of imaging
More informationDifferential of the Exponential Map
Differential of the Exponential Map Ethan Eade May 20, 207 Introduction This document computes ɛ0 log x + ɛ x ɛ where and log are the onential mapping and its inverse in a Lie group, and x and ɛ are elements
More informationGEOMETRICAL TRANSFORMATIONS IN HIGHER DIMENSIONAL EUCLIDEAN SPACES
GEOMETRICAL TRANSFORMATIONS IN HIGHER DIMENSIONAL EUCLIDEAN SPACES A Thesis by AMIT KUMAR SANYAL Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements
More informationOssama Abdelkhalik and Daniele Mortari Department of Aerospace Engineering, Texas A&M University,College Station, TX 77843, USA
Two-Way Orbits Ossama Abdelkhalik and Daniele Mortari Department of Aerospace Engineering, Texas A&M University,College Station, TX 77843, USA November 17, 24 Abstract. This paper introduces a new set
More informationIAC-11-C1.5.9 INERTIA-FREE ATTITUDE CONTROL OF SPACECRAFT WITH UNKNOWN TIME-VARYING MASS DISTRIBUTION
6nd International Astronautical Congress, Cape Town, SA. Copyright by the International Astronautical Federation. All rights reserved IAC--C.5.9 INERTIA-FREE ATTITUDE CONTROL OF SPACECRAFT WITH UNKNOWN
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 informationMethods in Computer Vision: Introduction to Matrix Lie Groups
Methods in Computer Vision: Introduction to Matrix Lie Groups Oren Freifeld Computer Science, Ben-Gurion University June 14, 2017 June 14, 2017 1 / 46 Definition and Basic Properties Definition (Matrix
More informationGeneralization of the Euler Angles
The Journal of the Astronautical Sciences, Vol. 51, No. 2, AprilJune 2003, pp. 123123 Generalization of the Euler Angles Malcolm D. Shuster 1 and F. Landis Markley 2 Abstract It is shown that the Euler
More informationKINEMATIC EQUATIONS OF NONNOMINAL EULER AXIS/ANGLE ROTATION
IAA-AAS-DyCoSS -14-10-01 KINEMATIC EQUATIONS OF NONNOMINAL EULER AXIS/ANGLE ROTATION Emanuele L. de Angelis and Fabrizio Giulietti INTRODUCTION Euler axis/angle is a useful representation in many attitude
More informationOn Sun-Synchronous Orbits and Associated Constellations
On Sun-Synchronous Orbits and Associated Constellations Daniele Mortari, Matthew P. Wilkins, and Christian Bruccoleri Department of Aerospace Engineering, Texas A&M University, College Station, TX 77843,
More informationGENERAL ARTICLE Realm of Matrices
Realm of Matrices Exponential and Logarithm Functions Debapriya Biswas Debapriya Biswas is an Assistant Professor at the Department of Mathematics, IIT- Kharagpur, West Bengal, India. Her areas of interest
More informationMoving Mass Control for Underwater Vehicles
Moving Mass Control for Underwater Vehicles C. A. Woolsey 1 & N. E. Leonard 2 1 Aerospace & Ocean Engineering Virginia Tech Blacksburg, VA 2461 cwoolsey@vt.edu 2 Mechanical & Aerospace Engineering Princeton
More informationUnit quaternion observer based attitude stabilization of a rigid spacecraft without velocity measurement
Proceedings of the 45th IEEE Conference on Decision & Control Manchester Grand Hyatt Hotel San Diego, CA, USA, December 3-5, 6 Unit quaternion observer based attitude stabilization of a rigid spacecraft
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 informationRotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep.
Rotational Motion Chapter 4 P. J. Grandinetti Chem. 4300 Sep. 1, 2017 P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. 1, 2017 1 / 76 Angular Momentum The angular momentum of a particle with respect
More informationSynodic and Relative Flower Constellations with Applications to Planetary Explorations
AAS 5-5 Synodic and Relative Flower Constellations with Applications to Planetary Explorations Daniele Mortari, Ossama Abdelkhalik, and Christian Bruccoleri Texas A&M University, College Station, Texas
More informationA Numerical Integration Scheme For The Dynamic Motion Of Rigid Bodies Using The Euler Parameters
International Journal of Research in Engineering and Science (IJRES) ISSN (Online): 30-9364, ISSN (Print): 30-9356 Volume Issue 8 ǁ August 014 ǁ PP30-37 A Numerical Integration Scheme For The Dynamic Motion
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 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 informationStabilization of Angular Velocity of Asymmetrical Rigid Body. Using Two Constant Torques
Stabilization of Angular Velocity of Asymmetrical Rigid Body Using Two Constant Torques Hirohisa Kojima Associate Professor Department of Aerospace Engineering Tokyo Metropolitan University 6-6, Asahigaoka,
More informationSpace Surveillance with Star Trackers. Part II: Orbit Estimation
AAS -3 Space Surveillance with Star Trackers. Part II: Orbit Estimation Ossama Abdelkhalik, Daniele Mortari, and John L. Junkins Texas A&M University, College Station, Texas 7783-3 Abstract The problem
More informationGlobal Formulations of Lagrangian and Hamiltonian Dynamics on Embedded Manifolds
1 Global Formulations of Lagrangian and Hamiltonian Dynamics on Embedded Manifolds By Taeyoung Lee, Melvin Leok, and N. Harris McClamroch Mechanical and Aerospace Engineering, George Washington University,
More informationAttitude Control of a Bias Momentum Satellite Using Moment of Inertia
I. INTRODUCTION Attitude Control of a Bias Momentum Satellite Using Moment of Inertia HYOCHOONG BANG Korea Advanced Institute of Science and Technology HYUNG DON CHOI Korea Aerospace Research Institute
More informationResearch Article Shorter Path Design and Control for an Underactuated Satellite
Hindawi International Journal of Aerospace Engineering Volume 27 Article ID 8536732 9 pages https://doi.org/.55/27/8536732 Research Article Shorter Path Design and Control for an Underactuated Satellite
More informationRADIALLY ADAPTIVE EVALUATION OF THE SPHERICAL HARMONIC GRAVITY SERIES FOR NUMERICAL ORBITAL PROPAGATION
AAS 15-440 RADIALLY ADAPTIVE EVALUATION OF THE SPHERICAL HARMONIC GRAVITY SERIES FOR NUMERICAL ORBITAL PROPAGATION Austin B. Probe, * Brent Macomber,* Julie I. Read,* Robyn M. Woollands,* and John L. Junkins
More informationLecture 38: Equations of Rigid-Body Motion
Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can
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 informationAutonomous Underwater Vehicles: Equations of Motion
Autonomous Underwater Vehicles: Equations of Motion Monique Chyba - November 18, 2015 Departments of Mathematics, University of Hawai i at Mānoa Elective in Robotics 2015/2016 - Control of Unmanned Vehicles
More informationRigid Body Equations of Motion for Modeling and Control of Spacecraft Formations - Part 1: Absolute Equations of Motion.
Rigid ody Equations of Motion for Modeling and Control of Spacecraft Formations - Part 1: Absolute Equations of Motion. Scott R. Ploen, Fred Y. Hadaegh, and Daniel P. Scharf Jet Propulsion Laboratory California
More informationSpacecraft Attitude Dynamics for Undergraduates
Session 1123 Spacecraft Attitude Dynamics for Undergraduates Dr. Rachel Shinn Embry Riddle Aeronautical University, Prescott, AZ Abstract Teaching spacecraft attitude dynamics to undergraduate students
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 informationSpacecraft Attitude Control with RWs via LPV Control Theory: Comparison of Two Different Methods in One Framework
Trans. JSASS Aerospace Tech. Japan Vol. 4, No. ists3, pp. Pd_5-Pd_, 6 Spacecraft Attitude Control with RWs via LPV Control Theory: Comparison of Two Different Methods in One Framework y Takahiro SASAKI,),
More informationLecture 38: Equations of Rigid-Body Motion
Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can
More information(3.1) a 2nd-order vector differential equation, as the two 1st-order vector differential equations (3.3)
Chapter 3 Kinematics As noted in the Introduction, the study of dynamics can be decomposed into the study of kinematics and kinetics. For the translational motion of a particle of mass m, this decomposition
More informationPhysics 5153 Classical Mechanics. Canonical Transformations-1
1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant
More information(1) for all (2) for all and all
8. Linear mappings and matrices A mapping f from IR n to IR m is called linear if it fulfills the following two properties: (1) for all (2) for all and all Mappings of this sort appear frequently in the
More informationAppendix Composite Point Rotation Sequences
Appendix Composite Point Rotation Sequences A. Euler Rotations In Chap. 6 we considered composite Euler rotations comprising individual rotations about the x, y and z axes such as R γ,x R β,y R α,z and
More informationSpacecraft Angular Rate Estimation Algorithms For Star Tracker-Based Attitude Determination
AAS 3-191 Spacecraft Angular Rate Estimation Algorithms For Star Tracker-Based Attitude Determination Puneet Singla John L. Crassidis and John L. Junkins Texas A&M University, College Station, TX 77843
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 information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationV Control Distance with Dynamics Parameter Uncertainty
V Control Distance with Dynamics Parameter Uncertainty Marcus J. Holzinger Increasing quantities of active spacecraft and debris in the space environment pose substantial risk to the continued safety of
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 informationProjective Bivector Parametrization of Isometries Part I: Rodrigues' Vector
Part I: Rodrigues' Vector Department of Mathematics, UACEG International Summer School on Hyprcomplex Numbers, Lie Algebras and Their Applications, Varna, June 09-12, 2017 Eulerian Approach The Decomposition
More information2 S. I. Kruglov, V. Barzda
Annales de la Fondation Louis de Broglie, Volume 4, numéro, 017 1 Modified Gibbs s representation of rotation matrix S. I. Kruglov a, and V. Barzda b Department of Physics, University of Toronto, 60 St.
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationOptimal Fault-Tolerant Configurations of Thrusters
Optimal Fault-Tolerant Configurations of Thrusters By Yasuhiro YOSHIMURA ) and Hirohisa KOJIMA, ) ) Aerospace Engineering, Tokyo Metropolitan University, Hino, Japan (Received June st, 7) Fault tolerance
More informationGlobal Trajectory Design for Position and Attitude Control of an Underactuated Satellite *
Trans. Japan Soc. Aero. Space Sci. Vol. 59, No. 3, pp. 7 4, 26 Global Trajectory Design for Position and Attitude Control of an Underactuated Satellite * Yasuhiro YOSHIMURA, ) Takashi MATSUNO, 2) and Shinji
More informationPhysics 312, Winter 2007, Practice Final
Physics 312, Winter 2007, Practice Final Time: Two hours Answer one of Question 1 or Question 2 plus one of Question 3 or Question 4 plus one of Question 5 or Question 6. Each question carries equal weight.
More informationCOUPLED ORBITAL AND ATTITUDE CONTROL SIMULATION
COUPLED ORBITAL AND ATTITUDE CONTROL SIMULATION Scott E. Lennox AOE 5984: Advanced Attitude Spacecraft Dynamics and Control December 12, 23 INTRODUCTION In the last few years the space industry has started
More informationRotational motion of a rigid body spinning around a rotational axis ˆn;
Physics 106a, Caltech 15 November, 2018 Lecture 14: Rotations The motion of solid bodies So far, we have been studying the motion of point particles, which are essentially just translational. Bodies with
More informationABSOLUTELY FLAT IDEMPOTENTS
ABSOLUTELY FLAT IDEMPOTENTS JONATHAN M GROVES, YONATAN HAREL, CHRISTOPHER J HILLAR, CHARLES R JOHNSON, AND PATRICK X RAULT Abstract A real n-by-n idempotent matrix A with all entries having the same absolute
More informationInstitute of Geometry, Graz, University of Technology Mobile Robots. Lecture notes of the kinematic part of the lecture
Institute of Geometry, Graz, University of Technology www.geometrie.tugraz.at Institute of Geometry Mobile Robots Lecture notes of the kinematic part of the lecture Anton Gfrerrer nd Edition 4 . Contents
More informationSpace Surveillance using Star Trackers. Part I: Simulations
AAS 06-231 Space Surveillance using Star Trackers. Part I: Simulations Iohan Ettouati, Daniele Mortari, and Thomas Pollock Texas A&M University, College Station, Texas 77843-3141 Abstract This paper presents
More informationLeonhard Euler ( September 1783)
LEONHARD EULER (1707 - SEPTEMBER 1783) Leonhard Euler (1707 - September 1783) BEYOND EQUATIONS Leonhard Euler was born in Basle, Switzerland; he was in fact a born mathematician, who went on to become
More informationSATELLITE ATTITUDE CONTROL SYSTEM DESIGN WITH NONLINEAR DYNAMICS AND KINEMTICS OF QUATERNION USING REACTION WHEELS
SATELLITE ATTITUDE CONTROL SYSTEM DESIGN WITH NONLINEAR DYNAMICS AND KINEMTICS OF QUATERNION USING REACTION WHEELS Breno Braga Galvao Maria Cristina Mendes Faustino Luiz Carlos Gadelha de Souza breno.braga.galvao@gmail.com
More informationMathematics. EC / EE / IN / ME / CE. for
Mathematics for EC / EE / IN / ME / CE By www.thegateacademy.com Syllabus Syllabus for Mathematics Linear Algebra: Matrix Algebra, Systems of Linear Equations, Eigenvalues and Eigenvectors. Probability
More informationAttitude Estimation Based on Solution of System of Polynomials via Homotopy Continuation
Attitude Estimation Based on Solution of System of Polynomials via Homotopy Continuation Yang Cheng Mississippi State University, Mississippi State, MS 3976 John L. Crassidis University at Buffalo, State
More informationPhys 201. Matrices and Determinants
Phys 201 Matrices and Determinants 1 1.1 Matrices 1.2 Operations of matrices 1.3 Types of matrices 1.4 Properties of matrices 1.5 Determinants 1.6 Inverse of a 3 3 matrix 2 1.1 Matrices A 2 3 7 =! " 1
More information3 Rigid Spacecraft Attitude Control
1 3 Rigid Spacecraft Attitude Control Consider a rigid spacecraft with body-fixed frame F b with origin O at the mass centre. Let ω denote the angular velocity of F b with respect to an inertial frame
More informationConvexity of the Joint Numerical Range
Convexity of the Joint Numerical Range Chi-Kwong Li and Yiu-Tung Poon October 26, 2004 Dedicated to Professor Yik-Hoi Au-Yeung on the occasion of his retirement. Abstract Let A = (A 1,..., A m ) be an
More informationStabilization of a Specified Equilibrium in the Inverted Equilibrium Manifold of the 3D Pendulum
Proceedings of the 27 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 27 ThA11.6 Stabilization of a Specified Equilibrium in the Inverted Equilibrium
More informationReview of Linear Algebra
Review of Linear Algebra Definitions An m n (read "m by n") matrix, is a rectangular array of entries, where m is the number of rows and n the number of columns. 2 Definitions (Con t) A is square if m=
More informationLecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More information7 Pendulum. Part II: More complicated situations
MATH 35, by T. Lakoba, University of Vermont 60 7 Pendulum. Part II: More complicated situations In this Lecture, we will pursue two main goals. First, we will take a glimpse at a method of Classical Mechanics
More informationTIME-OPTIMAL CONTROL OF AXI-SYMMETRIC RIGID SPACECRAFT
AIAA-98-438 TIME-OPTIMAL CONTROL OF AXI-SYMMETRIC RIGID SPACECRAFT Haijun Shen and Panagiotis Tsiotras University of Virginia, Charlottesville, VA 93-44 Abstract In this paper, we consider the minimum-time
More informationMP463 QUANTUM MECHANICS
MP463 QUANTUM MECHANICS Introduction Quantum theory of angular momentum Quantum theory of a particle in a central potential - Hydrogen atom - Three-dimensional isotropic harmonic oscillator (a model of
More informationPARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL
More informationGoldstein Problem 2.17 (3 rd ed. # 2.18)
Goldstein Problem.7 (3 rd ed. #.8) The geometry of the problem: A particle of mass m is constrained to move on a circular hoop of radius a that is vertically oriented and forced to rotate about the vertical
More informationThe kinetic equation (Lecture 11)
The kinetic equation (Lecture 11) January 29, 2016 190/441 Lecture outline In the preceding lectures we focused our attention on a single particle motion. In this lecture, we will introduce formalism for
More informationEigenstructure Assignment for Helicopter Hover Control
Proceedings of the 17th World Congress The International Federation of Automatic Control Eigenstructure Assignment for Helicopter Hover Control Andrew Pomfret Stuart Griffin Tim Clarke Department of Electronics,
More information2 Quantization of the Electromagnetic Field
2 Quantization of the Electromagnetic Field 2.1 Basics Starting point of the quantization of the electromagnetic field are Maxwell s equations in the vacuum (source free): where B = µ 0 H, D = ε 0 E, µ
More informationElectrons in a periodic potential
Chapter 3 Electrons in a periodic potential 3.1 Bloch s theorem. We consider in this chapter electrons under the influence of a static, periodic potential V (x), i.e. such that it fulfills V (x) = V (x
More informationMATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2
MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a
More informationRELATIVE NAVIGATION FOR SATELLITES IN CLOSE PROXIMITY USING ANGLES-ONLY OBSERVATIONS
(Preprint) AAS 12-202 RELATIVE NAVIGATION FOR SATELLITES IN CLOSE PROXIMITY USING ANGLES-ONLY OBSERVATIONS Hemanshu Patel 1, T. Alan Lovell 2, Ryan Russell 3, Andrew Sinclair 4 "Relative navigation using
More informationMATHEMATICAL METHODS INTERPOLATION
MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU
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 informationSingular Value Decomposition
Chapter 6 Singular Value Decomposition In Chapter 5, we derived a number of algorithms for computing the eigenvalues and eigenvectors of matrices A R n n. Having developed this machinery, we complete our
More informationStudy Notes on Matrices & Determinants for GATE 2017
Study Notes on Matrices & Determinants for GATE 2017 Matrices and Determinates are undoubtedly one of the most scoring and high yielding topics in GATE. At least 3-4 questions are always anticipated from
More informationA note on stability of spacecrafts and underwater vehicles
A note on stability of spacecrafts and underwater vehicles arxiv:1603.00343v1 [math-ph] 1 Mar 2016 Dan Comănescu Department of Mathematics, West University of Timişoara Bd. V. Pârvan, No 4, 300223 Timişoara,
More informationFormation Dynamics of Coulomb Satellites
Formation Dynamics of Coulomb Satellites Hyunsik Joe a, Hanspeter Schaub a and Gordon G. Parker b a Virginia Polytechnic Institute, Blacksburg, USA b Michigan Technological University, Houghton, USA Abstract
More informationECE 275A Homework #3 Solutions
ECE 75A Homework #3 Solutions. Proof of (a). Obviously Ax = 0 y, Ax = 0 for all y. To show sufficiency, note that if y, Ax = 0 for all y, then it must certainly be true for the particular value of y =
More informationSeminar 6: COUPLED HARMONIC OSCILLATORS
Seminar 6: COUPLED HARMONIC OSCILLATORS 1. Lagrangian Equations of Motion Let consider a system consisting of two harmonic oscillators that are coupled together. As a model, we will use two particles attached
More informationI. Perturbation Theory and the Problem of Degeneracy[?,?,?]
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Chemistry 5.76 Spring 19 THE VAN VLECK TRANSFORMATION IN PERTURBATION THEORY 1 Although frequently it is desirable to carry a perturbation treatment to second or third
More informationAnalytical Mechanics for Relativity and Quantum Mechanics
Analytical Mechanics for Relativity and Quantum Mechanics Oliver Davis Johns San Francisco State University OXPORD UNIVERSITY PRESS CONTENTS Dedication Preface Acknowledgments v vii ix PART I INTRODUCTION:
More informationLeast-squares Solutions of Linear Differential Equations
1 Least-squares Solutions of Linear Differential Equations Daniele Mortari dedicated to John Lee Junkins arxiv:1700837v1 [mathca] 5 Feb 017 Abstract This stu shows how to obtain least-squares solutions
More information