DOCSLIB.ORG

Metrics for 3D Rotations: Comparison and Analysis

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Metrics for 3D Rotations: Comparison and Analysis J Math Imaging Vis (2009) 35: 155–164 DOI 10.1007/s10851-009-0161-2 Metrics for 3D Rotations: Comparison and Analysis Du Q. Huynh Published online: 18 June 2009 © Springer Science+Business Media, LLC 2009 Abstract 3D rotations arise in many computer vision, com- efficiently represent 3D rotations and how to correctly eval- puter graphics, and robotics problems and evaluation of the uate the distance between them. If one of the 3D rotation distance between two 3D rotations is often an essential task. matrices is a true or reference rotation while the other is an This paper presents a detailed analysis of six functions for estimated one, then it is useful also to identify regions where measuring distance between 3D rotations that have been identical error measures occur. A few functions for distance proposed in the literature. Based on the well-developed the- measures between 3D rotations have been proposed in the ory behind 3D rotations, we demonstrate that five of them computer vision literature; however, there has neither been are bi-invariant metrics on SO(3) but that only four of them detailed analysis provided to these functions nor a compari- are boundedly equivalent to each other. We conclude that son of them in the context of the group SO(3) to which 3D it is both spatially and computationally more efficient to use rotations belong. The contributions of this paper are: (1) to quaternions for 3D rotations. Lastly, by treating the two rota- provide this missing information; (2) to analyze and illus- tions as a true and an estimated rotation matrix, we illustrate trate the iso-error contours of a given reference rotation. the geometry associated with iso-error measures. Rotations in 3D space can be represented in various forms. Euler angles are commonly used in robotics appli- Keywords Matrix Lie group · Lie algebra · Quaternions · cations where, because of constraints in the design of the 3D rotations · Distance functions joints of robot arms, rotations often have to be carried out in a certain order (e.g., [1];seealso[19]). For other ap- plications where such a constraint is absent, Euler angles 1 Introduction are less favoured, precisely because the values of these an- gles are dependent on the order of rotations about the three principal axes. For research in computer vision and com- 3D rotations are common entities in many computer vision, puter graphics, 3D rotations are commonly represented as computer graphics, and robotics problems that need to deal unit quaternions (e.g., [5, 15, 17, 18]), rotation axes and an- with the 3D world. Typical applications that involve 3D ro- gles (e.g., [8, 21]), or even as the 3 × 3 rotation matrices tations include the interpolation of trajectory of 3D orienta- themselves (e.g., [3, 4, 6]). The number of computer vision tions, robot kinematics, flight simulation, structure from mo- research papers that involve 3D rotations of any form is far tion, 3D pose recovery of objects, and motion capture. The too many to permit a complete list of citations. The refer- common issues that arise in these applications are how to ences given here are only a very small subset of papers in the literature. In terms of storage, each 3 × 3 rotation matrix requires This research was in part supported by a UWA study leave grant. the space of 9 floating point numbers, whereas in reality the D.Q. Huynh () special orthogonal group is a 3 dimensional object embed- School of Computer Science and Software Engineering, dable in R4. Another concern of the matrix representation The University of Western Australia, Nedlands, WA 6009, Australia is that after several matrix multiplications, round-off errors e-mail: [email protected] within computers can result in “rotation” matrices that are 156 J Math Imaging Vis (2009) 35: 155–164 no longer orthogonal. In that regard, the unit quaternions and noting that√ exp(tA) = I whenever t is an integer mul- are preferred as an alternative way of representing 3D ro- tiple of 2π/ 3. This example, however, suggests at least tations. Although round-off errors may also cause any unit some of the degree of non-uniqueness in the exponential quaternion to have a non-unit magnitude, it is more straight- parametrization. Indeed, the exponential parametrization of forward to renormalize it to unity in comparison with re- SO(3) corresponds exactly to the rotation axis and rota- orthogonalizing a distorted matrix. tion angle formulation. Given a rotation with rotation axis In the following sections, we will first give a brief u = (u1,u2,u3) of unit magnitude and rotation angle θ, overview of the special orthogonal group and the unit the 3 × 3 rotation matrix R can be obtained by the exponen- quaternions. This is then followed by the formal definition tial mapping exp([θu]×), where [ a ]× is chosen to be the of a distance function or metric. Six different functions for matrix defined by [a]×b = a × b. For example, the matrix in 3D rotations will then be studied in turn; their computation (4) corresponds to [(−1, 1, −1)]×. The group SO(3) is cov- complexity will be briefly analyzed. In Sect. 5, a geometrical ered by one-parameter groups (in fact circles) of the form interpretation of one of the functions will be studied further, {exp([θu]×) : θ ∈[−π,π)} and this representation is almost this leads to the discussion on iso-error contours. Finally in unique, since a rotation uniquely specifies its (unit vector) Sect. 6, we conclude the paper. rotation axis up to a multiplication by ±1, and once this is fixed, the angle is specified up to a multiple of 2π. It is straightforward to show that the Rodrigues formula 2 SO(3) and Unit Quaternions: An Overview (see, e.g., [11]) for a rotation matrix R as defined below 3D rotations form the so-called Special Orthogonal Group R = cos θ I + sin θ [u]× + (1 − cos θ)uu , (5) SO(3) of orthogonal matrices with determinant 1. SO(3) is a compact Lie group having the skew-symmetric matrices as is just a simplification of the exponential map given in (3) [ ] × its Lie algebra, so(3). This Lie algebra is a non-associative being applied to θu ×.Givena3 3 rotation matrix R, vector space equipped with a binary operation: the inverse of the exponential map provides a rotation an- gle/axis description of the rotation. Thus, log(R) is the [·, ·]:so(3) × so(3) → so(3), skew-symmetric matrix containing information about the ro- (1) tation axis and angle. Although the inverse process requires [A, B]=AB − BA a choice of rotation axis between the two alternatives, it is which can easily be seen to be a closed operation in so(3). a straightforward procedure to retrieve the rotation axis and The Lie algebra so(3) is the tangent space at the identity ele- angle (see Appendix A). ment of SO(3). The binary operation defined above is known As a unit quaternion, the same 3D rotation matrix [ ] ∈ ≡ ≡ Lie bracket exp( θu ×) SO(3) can be written as q (q0,q1,q2,q3) as the , which satisfies the following properties: (q0, q˜) [ ]=−[ ] =± θ θ A, B B, A , (cos 2 , u sin 2 ) . The unit quaternions are a one-to- (2) one parametrization of the Special Unitary group SU(2) and [A, [B, C]] + [C, [A, B]] + [B, [C, A]] = 0, the q above can be written as a 2 × 2 unitary matrix: for all A, B, C ∈ so(3). It follows immediately from the first q + iq q + iq property that [A, A]=0. 0 1 2 3 . (6) −q + iq q − iq In the theory of Lie groups, the exponential map is a map- 2 3 0 1 ping from the Lie algebra of a Lie group to the group itself. SU(2) provides a double-covering of SO(3); that is, each Such a mapping allows one to recapture the group structure rotation matrix in SO(3) corresponds to two members of from its Lie algebra. As in the general case of matrix Lie SU(2). The group homomorphism from SU(2) to SO(3) : → groups, the exponential map exp so(3) SO(3) is simply: has a two element kernel and this corresponds to the am- biguity in the choice of rotation axis in the logarithmic A2 A3 exp(A) = I + A + + +··· . (3) map. In this regard, we note that SU(2) and SO(3) have 2! 3! the same Lie algebra. The Lie algebra su(2), which con- For the special orthogonal group SO(3), the exponential map sists of the 2×2 skew-hermitian matrices, is isomorphic as a is surjective but not injective. The failure of injectivity is Lie algebra to so(3). Consider representing the rotation axis easily seen by considering a skew-symmetric matrix such as u = (u1,u2,u3) of unit magnitude and the rotation angle ⎡ ⎤ θ as a 2 × 2 skew-hermitian matrix, A, as follows: 011 A = ⎣−101⎦ (4) iθu θu + iθu A = 1 2 3 1 −10 −θu2 + iθu3 −iθu1 J Math Imaging Vis (2009) 35: 155–164 157 iu1 u2 + iu3 ˜ = θ ≡ θA. (7) = −u2 + iu3 −iu1 resp. Φ(R2R1, R3R1) Φ(R2, R3) (12) By applying the exponential mapping to A and noting that for R1, R2, R3 ∈ SO(3). A distance function is bi-invariant A2 =−θ 2I, we obtain a different version of the Rodrigues if it is both left and right invariant. formula: Two distance functions Φ and Ψ are said to be boundedly equivalent if there are positive real numbers a and b such exp : su(2) → SU(2), that sin θ exp(A) = cos θ I + A (8) aΦ(R , R ) ≤ Ψ(R , R ) ≤ bΦ(R , R ) (13) θ 1 2 1 2 1 2 ˜ = cos θ I + sin θ A.
Load more

Recommended publications
  • Quaternions and Cli Ord Geometric Algebras
    Quaternions and Cliord Geometric Algebras Robert Benjamin Easter First Draft Edition (v1) (c) copyright 2015, Robert Benjamin Easter, all rights reserved. Preface As a rst rough draft that has been put together very quickly, this book is likely to contain errata and disorganization. The references list and inline citations are very incompete, so the reader should search around for more references. I do not claim to be the inventor of any of the mathematics found here. However, some parts of this book may be considered new in some sense and were in small parts my own original research. Much of the contents was originally written by me as contributions to a web encyclopedia project just for fun, but for various reasons was inappropriate in an encyclopedic volume. I did not originally intend to write this book. This is not a dissertation, nor did its development receive any funding or proper peer review. I oer this free book to the public, such as it is, in the hope it could be helpful to an interested reader. June 19, 2015 - Robert B. Easter. (v1) [email protected] 3 Table of contents Preface . 3 List of gures . 9 1 Quaternion Algebra . 11 1.1 The Quaternion Formula . 11 1.2 The Scalar and Vector Parts . 15 1.3 The Quaternion Product . 16 1.4 The Dot Product . 16 1.5 The Cross Product . 17 1.6 Conjugates . 18 1.7 Tensor or Magnitude . 20 1.8 Versors . 20 1.9 Biradials . 22 1.10 Quaternion Identities . 23 1.11 The Biradial b/a .
    [Show full text]
  • Understanding Euler Angles 1
    Welcome Guest View Cart (0 items) Login Home Products Support News About Us Understanding Euler Angles 1. Introduction Attitude and Heading Sensors from CH Robotics can provide orientation information using both Euler Angles and Quaternions. Compared to quaternions, Euler Angles are simple and intuitive and they lend themselves well to simple analysis and control. On the other hand, Euler Angles are limited by a phenomenon called "Gimbal Lock," which we will investigate in more detail later. In applications where the sensor will never operate near pitch angles of +/‐ 90 degrees, Euler Angles are a good choice. Sensors from CH Robotics that can provide Euler Angle outputs include the GP9 GPS‐Aided AHRS, and the UM7 Orientation Sensor. Figure 1 ‐ The Inertial Frame Euler angles provide a way to represent the 3D orientation of an object using a combination of three rotations about different axes. For convenience, we use multiple coordinate frames to describe the orientation of the sensor, including the "inertial frame," the "vehicle‐1 frame," the "vehicle‐2 frame," and the "body frame." The inertial frame axes are Earth‐fixed, and the body frame axes are aligned with the sensor. The vehicle‐1 and vehicle‐2 are intermediary frames used for convenience when illustrating the sequence of operations that take us from the inertial frame to the body frame of the sensor. It may seem unnecessarily complicated to use four different coordinate frames to describe the orientation of the sensor, but the motivation for doing so will become clear as we proceed. For clarity, this application note assumes that the sensor is mounted to an aircraft.
    [Show full text]
  • Lecture Notes: Qubit Representations and Rotations
    Phys 711 Topics in Particles & Fields | Spring 2013 | Lecture 1 | v0.3 Lecture notes: Qubit representations and rotations Jeffrey Yepez Department of Physics and Astronomy University of Hawai`i at Manoa Watanabe Hall, 2505 Correa Road Honolulu, Hawai`i 96822 E-mail: [email protected] www.phys.hawaii.edu/∼yepez (Dated: January 9, 2013) Contents mathematical object (an abstraction of a two-state quan- tum object) with a \one" state and a \zero" state: I. What is a qubit? 1 1 0 II. Time-dependent qubits states 2 jqi = αj0i + βj1i = α + β ; (1) 0 1 III. Qubit representations 2 A. Hilbert space representation 2 where α and β are complex numbers. These complex B. SU(2) and O(3) representations 2 numbers are called amplitudes. The basis states are or- IV. Rotation by similarity transformation 3 thonormal V. Rotation transformation in exponential form 5 h0j0i = h1j1i = 1 (2a) VI. Composition of qubit rotations 7 h0j1i = h1j0i = 0: (2b) A. Special case of equal angles 7 In general, the qubit jqi in (1) is said to be in a superpo- VII. Example composite rotation 7 sition state of the two logical basis states j0i and j1i. If References 9 α and β are complex, it would seem that a qubit should have four free real-valued parameters (two magnitudes and two phases): I. WHAT IS A QUBIT? iθ0 α φ0 e jqi = = iθ1 : (3) Let us begin by introducing some notation: β φ1 e 1 state (called \minus" on the Bloch sphere) Yet, for a qubit to contain only one classical bit of infor- 0 mation, the qubit need only be unimodular (normalized j1i = the alternate symbol is |−i 1 to unity) α∗α + β∗β = 1: (4) 0 state (called \plus" on the Bloch sphere) 1 Hence it lives on the complex unit circle, depicted on the j0i = the alternate symbol is j+i: 0 top of Figure 1.
    [Show full text]
  • Arxiv:1001.0240V1 [Math.RA]
    Fundamental representations and algebraic properties of biquaternions or complexified quaternions Stephen J. Sangwine∗ School of Computer Science and Electronic Engineering, University of Essex, Wivenhoe Park, Colchester, CO4 3SQ, United Kingdom. Email: [email protected] Todd A. Ell† 5620 Oak View Court, Savage, MN 55378-4695, USA. Email: [email protected] Nicolas Le Bihan GIPSA-Lab D´epartement Images et Signal 961 Rue de la Houille Blanche, Domaine Universitaire BP 46, 38402 Saint Martin d’H`eres cedex, France. Email: [email protected] October 22, 2018 Abstract The fundamental properties of biquaternions (complexified quaternions) are presented including several different representations, some of them new, and definitions of fundamental operations such as the scalar and vector parts, conjugates, semi-norms, polar forms, and inner and outer products. The notation is consistent throughout, even between representations, providing a clear account of the many ways in which the component parts of a biquaternion may be manipulated algebraically. 1 Introduction It is typical of quaternion formulae that, though they be difficult to find, once found they are immediately verifiable. J. L. Synge (1972) [43, p34] arXiv:1001.0240v1 [math.RA] 1 Jan 2010 The quaternions are relatively well-known but the quaternions with complex components (complexified quaternions, or biquaternions1) are less so. This paper aims to set out the fundamental definitions of biquaternions and some elementary results, which, although elementary, are often not trivial. The emphasis in this paper is on the biquaternions as an applied algebra – that is, a tool for the manipulation ∗This paper was started in 2005 at the Laboratoire des Images et des Signaux (now part of the GIPSA-Lab), Grenoble, France with financial support from the Royal Academy of Engineering of the United Kingdom and the Centre National de la Recherche Scientifique (CNRS).
    [Show full text]
  • The Rolling Ball Problem on the Sphere
    S˜ao Paulo Journal of Mathematical Sciences 6, 2 (2012), 145–154 The rolling ball problem on the sphere Laura M. O. Biscolla Universidade Paulista Rua Dr. Bacelar, 1212, CEP 04026–002 S˜aoPaulo, Brasil Universidade S˜aoJudas Tadeu Rua Taquari, 546, CEP 03166–000, S˜aoPaulo, Brasil E-mail address: [email protected] Jaume Llibre Departament de Matem`atiques, Universitat Aut`onomade Barcelona 08193 Bellaterra, Barcelona, Catalonia, Spain E-mail address: [email protected] Waldyr M. Oliva CAMGSD, LARSYS, Instituto Superior T´ecnico, UTL Av. Rovisco Pais, 1049–0011, Lisbon, Portugal Departamento de Matem´atica Aplicada Instituto de Matem´atica e Estat´ıstica, USP Rua do Mat˜ao,1010–CEP 05508–900, S˜aoPaulo, Brasil E-mail address: [email protected] Dedicated to Lu´ıs Magalh˜aes and Carlos Rocha on the occasion of their 60th birthdays Abstract. By a sequence of rolling motions without slipping or twist- ing along arcs of great circles outside the surface of a sphere of radius R, a spherical ball of unit radius has to be transferred from an initial state to an arbitrary final state taking into account the orientation of the ball. Assuming R > 1 we provide a new and shorter prove of the result of Frenkel and Garcia in [4] that with at most 4 moves we can go from a given initial state to an arbitrary final state. Important cases such as the so called elimination of the spin discrepancy are done with 3 moves only. 1991 Mathematics Subject Classification. Primary 58E25, 93B27. Key words: Control theory, rolling ball problem.
    [Show full text]
  • Vectors, Matrices and Coordinate Transformations
    S. Widnall 16.07 Dynamics Fall 2009 Lecture notes based on J. Peraire Version 2.0 Lecture L3 - Vectors, Matrices and Coordinate Transformations By using vectors and defining appropriate operations between them, physical laws can often be written in a simple form. Since we will making extensive use of vectors in Dynamics, we will summarize some of their important properties. Vectors For our purposes we will think of a vector as a mathematical representation of a physical entity which has both magnitude and direction in a 3D space. Examples of physical vectors are forces, moments, and velocities. Geometrically, a vector can be represented as arrows. The length of the arrow represents its magnitude. Unless indicated otherwise, we shall assume that parallel translation does not change a vector, and we shall call the vectors satisfying this property, free vectors. Thus, two vectors are equal if and only if they are parallel, point in the same direction, and have equal length. Vectors are usually typed in boldface and scalar quantities appear in lightface italic type, e.g. the vector quantity A has magnitude, or modulus, A = |A|. In handwritten text, vectors are often expressed using the −→ arrow, or underbar notation, e.g. A , A. Vector Algebra Here, we introduce a few useful operations which are defined for free vectors. Multiplication by a scalar If we multiply a vector A by a scalar α, the result is a vector B = αA, which has magnitude B = |α|A. The vector B, is parallel to A and points in the same direction if α > 0.
    [Show full text]
  • Chapter 5 ANGULAR MOMENTUM and ROTATIONS
    Chapter 5 ANGULAR MOMENTUM AND ROTATIONS In classical mechanics the total angular momentum L~ of an isolated system about any …xed point is conserved. The existence of a conserved vector L~ associated with such a system is itself a consequence of the fact that the associated Hamiltonian (or Lagrangian) is invariant under rotations, i.e., if the coordinates and momenta of the entire system are rotated “rigidly” about some point, the energy of the system is unchanged and, more importantly, is the same function of the dynamical variables as it was before the rotation. Such a circumstance would not apply, e.g., to a system lying in an externally imposed gravitational …eld pointing in some speci…c direction. Thus, the invariance of an isolated system under rotations ultimately arises from the fact that, in the absence of external …elds of this sort, space is isotropic; it behaves the same way in all directions. Not surprisingly, therefore, in quantum mechanics the individual Cartesian com- ponents Li of the total angular momentum operator L~ of an isolated system are also constants of the motion. The di¤erent components of L~ are not, however, compatible quantum observables. Indeed, as we will see the operators representing the components of angular momentum along di¤erent directions do not generally commute with one an- other. Thus, the vector operator L~ is not, strictly speaking, an observable, since it does not have a complete basis of eigenstates (which would have to be simultaneous eigenstates of all of its non-commuting components). This lack of commutivity often seems, at …rst encounter, as somewhat of a nuisance but, in fact, it intimately re‡ects the underlying structure of the three dimensional space in which we are immersed, and has its source in the fact that rotations in three dimensions about di¤erent axes do not commute with one another.
    [Show full text]
  • Hypercomplex Numbers
    Hypercomplex numbers Johanna R¨am¨o Queen Mary, University of London [email protected] We have gradually expanded the set of numbers we use: first from finger counting to the whole set of positive integers, then to positive rationals, ir- rational reals, negatives and finally to complex numbers. It has not always been easy to accept new numbers. Negative numbers were rejected for cen- turies, and complex numbers, the square roots of negative numbers, were considered impossible. Complex numbers behave like ordinary numbers. You can add, subtract, multiply and divide them, and on top of that, do some nice things which you cannot do with real numbers. Complex numbers are now accepted, and have many important applications in mathematics and physics. Scientists could not live without complex numbers. What if we take the next step? What comes after the complex numbers? Is there a bigger set of numbers that has the same nice properties as the real numbers and the complex numbers? The answer is yes. In fact, there are two (and only two) bigger number systems that resemble real and complex numbers, and their discovery has been almost as dramatic as the discovery of complex numbers was. 1 Complex numbers Complex numbers where discovered in the 15th century when Italian math- ematicians tried to find a general solution to the cubic equation x3 + ax2 + bx + c = 0: At that time, mathematicians did not publish their results but kept them secret. They made their living by challenging each other to public contests of 1 problem solving in which the winner got money and fame.
    [Show full text]
  • Euler Quaternions
    Four different ways to represent rotation my head is spinning... The space of rotations SO( 3) = {R ∈ R3×3 | RRT = I,det(R) = +1} Special orthogonal group(3): Why det( R ) = ± 1 ? − = − Rotations preserve distance: Rp1 Rp2 p1 p2 Rotations preserve orientation: ( ) × ( ) = ( × ) Rp1 Rp2 R p1 p2 The space of rotations SO( 3) = {R ∈ R3×3 | RRT = I,det(R) = +1} Special orthogonal group(3): Why it’s a group: • Closed under multiplication: if ∈ ( ) then ∈ ( ) R1, R2 SO 3 R1R2 SO 3 • Has an identity: ∃ ∈ ( ) = I SO 3 s.t. IR1 R1 • Has a unique inverse… • Is associative… Why orthogonal: • vectors in matrix are orthogonal Why it’s special: det( R ) = + 1 , NOT det(R) = ±1 Right hand coordinate system Possible rotation representations You need at least three numbers to represent an arbitrary rotation in SO(3) (Euler theorem). Some three-number representations: • ZYZ Euler angles • ZYX Euler angles (roll, pitch, yaw) • Axis angle One four-number representation: • quaternions ZYZ Euler Angles φ = θ rzyz ψ φ − φ cos sin 0 To get from A to B: φ = φ φ Rz ( ) sin cos 0 1. Rotate φ about z axis 0 0 1 θ θ 2. Then rotate θ about y axis cos 0 sin θ = ψ Ry ( ) 0 1 0 3. Then rotate about z axis − sinθ 0 cosθ ψ − ψ cos sin 0 ψ = ψ ψ Rz ( ) sin cos 0 0 0 1 ZYZ Euler Angles φ θ ψ Remember that R z ( ) R y ( ) R z ( ) encode the desired rotation in the pre- rotation reference frame: φ = pre−rotation Rz ( ) Rpost−rotation Therefore, the sequence of rotations is concatentated as follows: (φ θ ψ ) = φ θ ψ Rzyz , , Rz ( )Ry ( )Rz ( ) φ − φ θ θ ψ − ψ cos sin 0 cos 0 sin cos sin 0 (φ θ ψ ) = φ φ ψ ψ Rzyz , , sin cos 0 0 1 0 sin cos 0 0 0 1− sinθ 0 cosθ 0 0 1 − − − cφ cθ cψ sφ sψ cφ cθ sψ sφ cψ cφ sθ (φ θ ψ ) = + − + Rzyz , , sφ cθ cψ cφ sψ sφ cθ sψ cφ cψ sφ sθ − sθ cψ sθ sψ cθ ZYX Euler Angles (roll, pitch, yaw) φ − φ cos sin 0 To get from A to B: φ = φ φ Rz ( ) sin cos 0 1.
    [Show full text]
  • Theory of Angular Momentum and Spin
    Chapter 5 Theory of Angular Momentum and Spin Rotational symmetry transformations, the group SO(3) of the associated rotation matrices and the 1 corresponding transformation matrices of spin{ 2 states forming the group SU(2) occupy a very important position in physics. The reason is that these transformations and groups are closely tied to the properties of elementary particles, the building blocks of matter, but also to the properties of composite systems. Examples of the latter with particularly simple transformation properties are closed shell atoms, e.g., helium, neon, argon, the magic number nuclei like carbon, or the proton and the neutron made up of three quarks, all composite systems which appear spherical as far as their charge distribution is concerned. In this section we want to investigate how elementary and composite systems are described. To develop a systematic description of rotational properties of composite quantum systems the consideration of rotational transformations is the best starting point. As an illustration we will consider first rotational transformations acting on vectors ~r in 3-dimensional space, i.e., ~r R3, 2 we will then consider transformations of wavefunctions (~r) of single particles in R3, and finally N transformations of products of wavefunctions like j(~rj) which represent a system of N (spin- Qj=1 zero) particles in R3. We will also review below the well-known fact that spin states under rotations behave essentially identical to angular momentum states, i.e., we will find that the algebraic properties of operators governing spatial and spin rotation are identical and that the results derived for products of angular momentum states can be applied to products of spin states or a combination of angular momentum and spin states.
    [Show full text]
  • Split Quaternions and Spacelike Constant Slope Surfaces in Minkowski 3- Space
    Split Quaternions and Spacelike Constant Slope Surfaces in Minkowski 3- Space Murat Babaarslan and Yusuf Yayli Abstract. A spacelike surface in the Minkowski 3-space is called a constant slope surface if its position vector makes a constant angle with the normal at each point on the surface. These surfaces completely classified in [J. Math. Anal. Appl. 385 (1) (2012) 208-220]. In this study, we give some relations between split quaternions and spacelike constant slope surfaces in Minkowski 3-space. We show that spacelike constant slope surfaces can be reparametrized by using rotation matrices corresponding to unit timelike quaternions with the spacelike vector parts and homothetic motions. Subsequently we give some examples to illustrate our main results. Mathematics Subject Classification (2010). Primary 53A05; Secondary 53A17, 53A35. Key words: Spacelike constant slope surface, split quaternion, homothetic motion. 1. Introduction Quaternions were discovered by Sir William Rowan Hamilton as an extension to the complex number in 1843. The most important property of quaternions is that every unit quaternion represents a rotation and this plays a special role in the study of rotations in three- dimensional spaces. Also quaternions are an efficient way understanding many aspects of physics and kinematics. Many physical laws in classical, relativistic and quantum mechanics can be written nicely using them. Today they are used especially in the area of computer vision, computer graphics, animations, aerospace applications, flight simulators, navigation systems and to solve optimization problems involving the estimation of rigid body transformations. Ozdemir and Ergin [9] showed that a unit timelike quaternion represents a rotation in Minkowski 3-space.
    [Show full text]
  • Quaternion Algebra and Calculus
    Quaternion Algebra and Calculus David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/ or send a letter to Creative Commons, PO Box 1866, Mountain View, CA 94042, USA. Created: March 2, 1999 Last Modified: August 18, 2010 Contents 1 Quaternion Algebra 2 2 Relationship of Quaternions to Rotations3 3 Quaternion Calculus 5 4 Spherical Linear Interpolation6 5 Spherical Cubic Interpolation7 6 Spline Interpolation of Quaternions8 1 This document provides a mathematical summary of quaternion algebra and calculus and how they relate to rotations and interpolation of rotations. The ideas are based on the article [1]. 1 Quaternion Algebra A quaternion is given by q = w + xi + yj + zk where w, x, y, and z are real numbers. Define qn = wn + xni + ynj + znk (n = 0; 1). Addition and subtraction of quaternions is defined by q0 ± q1 = (w0 + x0i + y0j + z0k) ± (w1 + x1i + y1j + z1k) (1) = (w0 ± w1) + (x0 ± x1)i + (y0 ± y1)j + (z0 ± z1)k: Multiplication for the primitive elements i, j, and k is defined by i2 = j2 = k2 = −1, ij = −ji = k, jk = −kj = i, and ki = −ik = j. Multiplication of quaternions is defined by q0q1 = (w0 + x0i + y0j + z0k)(w1 + x1i + y1j + z1k) = (w0w1 − x0x1 − y0y1 − z0z1)+ (w0x1 + x0w1 + y0z1 − z0y1)i+ (2) (w0y1 − x0z1 + y0w1 + z0x1)j+ (w0z1 + x0y1 − y0x1 + z0w1)k: Multiplication is not commutative in that the products q0q1 and q1q0 are not necessarily equal.
    [Show full text]