Quaternions. R. J. Renka 11/09/2015. Department of Computer Science & Engineering University of North Texas. R. J.
|
|
- Susanna Parsons
- 6 years ago
- Views:
Transcription
1 Quaternions R. J. Renka Department of Computer Science & Engineering University of North Texas 11/09/2015
2 Definition A quaternion is an element of R 4 with three operations: addition, scalar multiplication, and quaternion multiplication. The first two operations are identical to those on R 4. Multiplication is defined in terms of basis elements, 1, i, j, k, where 1 is the multiplicative identity element: a quaternion can be written as a linear combination q = w + xi + yj + zk, where we have suppressed 1 = (1, 0, 0, 0), and i 2 = j 2 = k 2 = ijk = 1, and multiplication is associative, so that ij = k, jk = i, ki = j, ji = k, kj = i, ik = j. This is a four-dimensional analog of complex numbers. Notice that i, j, and k are unit quaternions that behave like vectors in R 3 with vector cross product, except that cross product is not associative.
3 Quaternion Operations It is convenient to write a quaternion as a scalar/vector pair Then it is easily verified that q = (w, v), for v = (x, y, z). q 1 q 2 = (w 1 w 2 v 1, v 2, w 1 v 2 + w 2 v 1 + v 1 v 2 ). Note that quaternion multiplication is commutative if and only if the vector parts are linearly dependent. We could have used the above expression to define multiplication, and then derived i 2 = (0, e 1 )(0, e 1 ) = ( 1, 0) = 1, ij = (0, e 1 )(0, e 2 ) = (0, e 3 ) = k, etc. We also have q 1 + q 2 = (w 1, v 1 ) + (w 2, v 2 ) = (w 1 + w 2, v 1 + v 2 ) and α(w, v) = (αw, αv) α R.
4 Conjugate, Norm, and Reciprocal The conjugate of q = (w,v) is q = (w, v), and the norm of q is q = w 2 + v, v, so that q q = (w 2 + v, v, 0) = w 2 + v, v = q 2. The reciprocal of q is q 1 = q/ q 2. Since qq = q q, q 1 is both a left inverse and right inverse. However, pq 1 q 1 p, and the expression p/q is therefore ambiguous. A unit quaternion is one whose norm is 1. Its reciprocal is its conjugate. It is straightforward to show that q 1 q 2 = q 1 q 2 and (q 1 q 2 ) 1 = q 1 2 q 1 1. The subset of purely imaginary quaternions with zero scalar parts includes the zero quaternion and is closed under linear combinations (addition and scalar multiplication). It is therefore a three-dimensional linear space over the reals, isomorphic to R 3.
5 Unit Quaternions and Rotation The complex number z = x + iy = ρe iφ = ρ(cos φ + i sin φ) can be rotated through angle θ by multiplying it by the unit complex number e iθ : e iθ z = ρe i(φ+θ) = cos θ x sin θ y + i(sin θ x + cos θ y). Similarly, a unit quaternion represents a rotation about an arbitrary pole (unit vector defining an axis direction) in R 3. Theorem A vector p R 3 can be rotated through angle θ about a pole n, n = 1, by forming the unit quaternion q = (cos θ 2, sin θ 2 n), and computing qpq 1, where p = (0,p) is the purely imaginary quaternion associated with vector p.
6 Proof of Theorem qp = = qpq 1 = = (cos θ2, sin θ2 ) n (0, p) ( sin θ2 nt p, sin θ2 (n p) + cos θ2 ) p ( sin θ2 nt p, sin θ2 (n p) + cos θ2 p ) (cos θ2, sin θ2 n ) ( sin θ 2 cos θ 2 nt p + sin θ 2 cos θ 2 nt p, sin 2 θ 2 (n p n) cos θ 2 sin θ (p n) 2 + sin 2 θ ( ) n T p n + cos θ 2 2 sin θ 2 (n p) + θ ) cos2 2 p
7 Proof of Theorem continued qp q = ( 0, sin 2 θ 2 p + 2 θ sin2 2 + cos 2 θ ) 2 p ( ) n T p n + 2 cos θ 2 sin θ 2 = (0, cos θp + (1 cos θ)(n T p)n + sin θ(n p)) = (0, Rp) for rotation matrix R, where we have used the trig identites ( θ cos θ = cos 2 + θ ) = cos 2 θ 2 2 θ sin2 2, 1 cos θ = 2 sin 2 θ 2, and ( θ sin θ = sin 2 + θ ) = 2 cos θ 2 2 sin θ 2. (n p)
8 Rotation Matrix R for Pole n p (n T p)n n φ 0 p p (n T p)n θ n p + (n T p)n Let p = Rp. Then p (n T p)n = cos θ(p (n T p)n) + sin θ (n p), where n p = p sin φ = p (n T p)n. Hence p = (cos θ)p + (1 cos θ) (n T p)n + sin θ(n p).
9 Alternative Expression for R Since the vector cross product is bilinear, n p can be computed by applying a linear transformation L to p. Denote the (skew symmetric) matrix representing L by ˆn. Then it is easily shown that ˆn = 0 n 3 n 2 n 3 0 n 1 n 2 n 1 0 We now have the following expression for R:. R = cos θi + (1 cos θ)nn T + sin θˆn. Also, ˆn 2 p = n (n p) = n p n = n 2 p+ (n T p)n = p + nn T p. Hence ˆn 2 = I + nn T, and nn T = I + ˆn 2, so that R = cos θi + (1 cos θ)(i + ˆn 2 ) + sin θˆn = I + sin θˆn + (1 cos θ)ˆn 2.
10 Properties of Cross Product Operator For u R 3, define Then û = ûv = u v v R 3. û T = û. 0 u 3 u 2 u 3 0 u 1 u 2 u 1 0. αu = αû and û + v = û + ˆv for scalars α and vectors u, v. û (t) = û (t), and u (t) u(t) = û (t)û(t) û(t)û (t). û 3 = û, u(t). û 4 = û 2, and û(t)û (t)û(t) = 0 for unit vectors
11 R in Terms of the Components of q = (w,x,y,z) In terms of the components of q = (w, v) = (cos(θ/2), sin(θ/2)n), we have where ˆv = so that R = I + 2 cos θ 2 sin θ 2 ˆn + 2 θ sin2 2 ˆn2 = I + 2w sin θ ( 2 ˆn + 2 sin θ ) 2 2 ˆn 0 z y z 0 x y x 0 R = = I + 2wˆv + 2ˆv 2,, ˆv 2 = y 2 z 2 xy xz xy x 2 z 2 yz xz yz x 2 y 2 1 2y 2 2z 2 2xy 2wz 2xz + 2wy 2xy + 2wz 1 2x 2 2z 2 2yz 2wx 2xz 2wy 2yz + 2wx 1 2x 2 2y 2.,
12 Algebraic Derivation of R Let a, b, and n be orthonormal unit vectors defining a right hand coordinate system, and let c = cos θ, s = sin θ. Then the matrix R that rotates through angle θ (CCW as viewed toward the origin) about n satisfies Ra = ca + sb, Rb = -sa + cb, and Rn = n, or R(a b n) = (a b n) c s 0 s c ; i.e., RP = PR z for rotations P and R z. Hence R = PR z P T = ( a b a T n )R z b T n T = ( a b n ) = a(ca T sb T ) + b(sa T + cb T ) + nn T = c(aa T + bb T ) + s(ba T ab T ) + nn T. ca T sb T sa T + cb T n T
13 Algebraic Derivation of R continued We will rewrite R in terms of n, ˆn, and ˆn 2 instead of a, b, and n. Any vector p in R 3 may be written p = αa + βb + γn for α = a T p, β = b T p, and γ = n T p. Then ˆnp = n p = αb βa = ba T p ab T p, and hence ˆn = ba T ab T. Similarly, ˆn 2 p = αa βb = aa T p bb T p, and I p =αa + βb + γn = aa T p + bb T p + nn T p giving I = aa T + bb T + nn T, and Thus, we have ˆn 2 = aa T bb T = nn T I. R = cˆn 2 + sˆn + ˆn 2 + I = I + sˆn + (1 c)ˆn 2.
14 Properties of Unit Quaternions The product of unit quaternions is a unit quaternion. (cos(θ/2), sin(θ/2) n)(cos(φ/2), sin(φ/2)n) = (cos(θ/2 + φ/2), sin(θ/2 +φ/2)n) for n = 1; i.e., rotation through angle φ followed by rotation through angle θ is equivalent to rotation through angle θ + φ. More generally, unit quaternions commute iff the poles coincide. Unit quaternions i, j, and k represent 180-degree axis rotations, and i 2 = -1 = (-1,0) represents a 360-degree rotation (a special case of the following property). In fact every purely imaginary unit quaternion q = (0, n) represents a 180-degree rotation R = 2nn T I and satisfies q 2 = 1. q and -q represent the same rotation: q(0, p)q 1 = q(0, p)( q) 1 for all p. Rotation by q 1 followed by rotation by q 2 is equivalent to rotation by q 2 q 1.
15 Euler s Theorem Theorem Any orientation defined by a rotation matrix (order-3 orthogonal matrix R with det(r) = 1) can be represented by a rotation about a single Euler pole n; i.e., Rn = n for some nonzero vector n (and R has 1 as an eigenvalue). In particular, any product of rotations in R 3 can be represented by a rotation about a single Euler pole. proof It suffices to prove that 1 is an eigenvalue of R. det(i R) = det((i R) T ) = det(i R T ) = det( R T (I R)) = det(r T ) det(i R) = det(i R), and hence det(i R) = 0.
16 Euler Angles Any rigid-body orientation can be represented by three Euler angles or, equivalently, by a product of three rotations. An interactive user must be able to incrementally rotate an object about the object axes to change heading (yaw), pitch, or bank (roll), for example, or to smoothly rotate the camera or eye position around a scene. In OpenGL this could be accomplished by left multiplying the modelview matrix M by an axis rotation, but there are two problems with that. First, it would be inefficient because the OpenGL rotate function right multiplies M by the rotation. Second, we would be accumulating roundoff error, requiring periodic re-orthogonalization. The solution is to represent the viewing transformation by a single pole and angle, convert the axis rotation to a pole and angle, and compute a product of unit quaternions to update the viewing transformation.
17 Viewing Direction Problem: Given an arbitrary viewing direction u = (E O)/ E O, for eye position E and object point O, construct a rotation (pole n and angle θ) that maps u to e 3 so that the viewing direction is along the z axis. Solution: n = u e 3 = (u 2, u 1, 0), θ = cos 1 u, e 3 = cos 1 u 3.
18 Rigid Body Animation Given key frames associated with different orientations, we need to construct in-between frames for animation by interpolation. If the orientations are defined by pairs of axis rotations or triples of Euler angles, it is not clear how to choose in-between frames. In terms of a single pole and angle, or quaternion, however, we have a one dimensional parameter space, and we simply use a sequence of angles along with the fixed pole. The more general problem is to produce constant-speed motion along a unit-radius arc of a great circle in two or more dimensions. This is referred to as spherical linear interpolation (SLERP). Given unit quaternions (points on the unit hypersphere in R 4 ) representing an initial orientation and a final orientation, we seek to construct a sequence of unit quaternions representing intermediate orientations on the arc of great circle joining the given pair.
19 Spherical Linear Interpolation (SLERP) Given a pair of unit vectors q 0 and q 1 in R n, n 2, separated by angle θ = cos 1 (q 0 q 1 ) with 0 < θ < π, we wish to construct a constant-speed parameterization {q(t) : t [0, 1]} of the geodesic that joins q 0 and q 1. The solution is characterized as follows. 1 q(t) is a linear combination of q 0 and q 1, 2 q(t) = 1, and 3 the angular distance from q 0 to q(t) is tθ for all t [0, 1]. Let q(t) = α 0 q 0 + α 1 q 1. Then q(t) 2 = α α 0 α 1 cos θ + α 2 1 = 1 and q 0 q(t) = α 0 + α 1 cos θ = cos tθ α 0 = cos tθ α 1 cos θ, so that 1 = cos 2 tθ 2α 1 cos θ cos tθ + α1 2 cos 2 θ + 2α 1 cos θ (cos tθ α 1 cos θ) + α1 2 = cos 2 tθ α1 2 cos 2 θ + α1. 2
20 Spherical Linear Interpolation (SLERP) continued We thus have 1 cos 2 tθ = α 2 1(1 cos 2 θ) sin 2 tθ = α 2 1 sin 2 θ, and hence α 1 = sin tθ/ sin θ = sin tθ/ sin θ since 0 < θ < π. By choosing α 1 0, α 0 and α 1 will be in the interval [0, 1], and q(t) will lie on the shorter of the two great circle arcs defined by q 0 and q 1. We now have sin tθ sin θ cos tθ cos θ sin tθ α 0 = cos tθ cos θ = sin θ sin θ i.e., = sin(θ tθ) ; sin θ q(t) = slerp(t; q 0, q 1 ) = sin[(1 t)θ] sin tθ q 0 + sin θ sin θ q 1. Note the symmetry: slerp(t;q 0, q 1 ) = slerp(1-t;q 1, q 0 ). Also, in the limit as θ 0, slerp(t;q 0, q 1 ) (1 t)q 0 + tq 1.
21 Spherical Linear Interpolation (SLERP) continued An alternative derivation of the expression for q(t) is based on geometry (and trigonometry). Let q0 denote the 90-degree rotation of q 0 in the direction toward q 1. Then q(t) = (cos tθ)q 0 + (sin tθ)q0, and q 1 = (cos θ)q 0 + (sin θ)q0 implying that q0 = (q 1 (cos θ)q 0 )/ sin θ. Hence ( ) q1 (cos θ)q 0 q(t) = (cos tθ)q 0 + sin tθ sinθ = (sin θ cos tθ cos θ sin tθ)q 0 + (sin tθ)q 1 sin θ sin[θ(1 t)] sin tθ = q 0 + sin θ sin θ q 1. It is a straightforward exercise to show that the derivative satisfies q (t) = θ t, so that q (t)/θ is a unit tangent to the unit hypersphere at q(t) for all t.
22 Spherical Linear Interpolation (SLERP) continued There remains an additional issue in the application of SLERP to unit quaternions the double coverage of the set of rotations. Each rotation is represented by two unit quaternions, q and q. Thus, slerp(t;q 0, q 1 ) and slerp(t;q 0, q 1 ) both correspond to great circle arcs associated with the same pair of rotations, but only the shorter of the two arcs is the geodesic. We therefore begin by reversing the sign of either q 0 or q 1 if q 0 q 1 < 0. Then we have θ [0, π/2].
23 Time-varying Quaternion Given a unit quaternion q(t) representing the orientation of a rigid body at time t, we need to evaluate the derivative q (t). We require the following. Lemma: Given a time-dependent angle θ(t) and unit vector u(t) defining a rotation R(t) = I + sin(θ(t))û(t) + (1 cos(θ(t)))û(t) 2 with derivative the angular velocity is R (t) = ˆω(t)R(t), ω(t) = θu + sin(θ) u + (cos(θ) 1)( u u).
24 Proof of Lemma Differentiating R(t) = I + sin(θ(t))û(t) + (1 cos(θ(t)))û(t) 2, we get Ṙ = θ cos(θ)û + sin(θ) û + (1 cos(θ))(û û + ûû) + θ sin(θ)û 2, and R T = I sin(θ)û + (1 cos(θ))û 2. Using the properties of the cross product operator for unit vectors, and a bit of algebra, we obtain ṘR T = θû + sin(θ) û + (1 cos(θ))(û û ûû) = ˆω.
25 Derivative of Time-varying Unit Quaternion Theorem: q (t) = (1/2)(0, ω(t))q(t). proof: Let q(t) = γ + σu for γ = cos(θ(t)/2), σ = sin(θ(t)/2), and u = (0, u(t)) for u(t) = 1 for all t. Then Also, Thus q (t) = ( θ/2)σ + σ u + ( θ/2)γu = ( θ/2)(γu σ) + σ u = ( θ/2)uq + σ u = (1/2)( θu + 2σ uq 1 )q. uq 1 = u(γ σu) = γ u σ uu = γ u σ( u T u + (0, u u)) = γ u σ(0, u u). q (t) = (1/2)[ θu + 2σγ u 2σ 2 (0, u u)]q = (1/2)[ θu + sin(θ) u + (cos(θ) 1)(0, u u)]q = (1/2)(0, ω)q.
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
More informationCourse 2BA1: Hilary Term 2007 Section 8: Quaternions and Rotations
Course BA1: Hilary Term 007 Section 8: Quaternions and Rotations David R. Wilkins Copyright c David R. Wilkins 005 Contents 8 Quaternions and Rotations 1 8.1 Quaternions............................ 1 8.
More informationWhy Transforms? Want to animate objects and camera Translations Rotations Shears And more.. Want to be able to use projection transforms
Why Transforms? Want to animate objects and camera Translations Rotations Shears And more.. Want to be able to use projection transforms ITCS 3050:Game Engine Programming 1 Geometric Transformations Implementing
More informationBe fruitful and multiply Genesis 1:28
Lecture 15: Quaternions: Multiplication in the Space of Mass-Points Be fruitful and multiply Genesis 1:8 15.1 Vector Spaces and Division Algebras The space of mass-points is a 4-dimensional vector space:
More informationBe fruitful and multiply Genesis 1:28
Chapter 15: Quaternions: Multiplication in the Space of Mass-Points Be fruitful and multiply Genesis 1:8 1. Vector Spaces and Division Algebras The space of mass-points is a 4-dimensional vector space:
More informationA Linear Algebraic Approach to Quaternions
A Linear Algebraic Approach to Quaternions David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International
More informationQUATERNIONS AND ROTATIONS
QUATERNIONS AND ROTATIONS SVANTE JANSON 1. Introduction The purpose of this note is to show some well-known relations between quaternions and the Lie groups SO(3) and SO(4) (rotations in R 3 and R 4 )
More informationCourse MA2C02, Hilary Term 2010 Section 4: Vectors and Quaternions
Course MA2C02, Hilary Term 2010 Section 4: Vectors and Quaternions David R. Wilkins Copyright c David R. Wilkins 2000 2010 Contents 4 Vectors and Quaternions 47 4.1 Vectors...............................
More informationIntroduction to quaternions
. Introduction Introduction to uaternions Invented and developed by William Hamilton in 843, uaternions are essentially a generalization of complex numbers to four dimensions (one real dimension, three
More informationComplex Numbers and Quaternions for Calc III
Complex Numbers and Quaternions for Calc III Taylor Dupuy September, 009 Contents 1 Introduction 1 Two Ways of Looking at Complex Numbers 1 3 Geometry of Complex Numbers 4 Quaternions 5 4.1 Connection
More information1 Vectors and Tensors
PART 1: MATHEMATICAL PRELIMINARIES 1 Vectors and Tensors This chapter and the next are concerned with establishing some basic properties of vectors and tensors in real spaces. The first of these is specifically
More informationAlgebra I Fall 2007
MIT OpenCourseWare http://ocw.mit.edu 18.701 Algebra I Fall 007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 18.701 007 Geometry of the Special Unitary
More informationQuaternions 2 AUI Course Denbigh Starkey
Quaternions 2 AUI Course Denbigh Starkey 1. Background 2 2. Some Basic Quaternion Math 4 3. The Justification of the Quaternion Rotation Formula 5 4. Interpolation between two Unit Quaternions SLERP vs.
More informationPHYS 705: Classical Mechanics. Rigid Body Motion Introduction + Math Review
1 PHYS 705: Classical Mechanics Rigid Body Motion Introduction + Math Review 2 How to describe a rigid body? Rigid Body - a system of point particles fixed in space i r ij j subject to a holonomic constraint:
More informationLinear Algebra Review. Vectors
Linear Algebra Review 9/4/7 Linear Algebra Review By Tim K. Marks UCSD Borrows heavily from: Jana Kosecka http://cs.gmu.edu/~kosecka/cs682.html Virginia de Sa (UCSD) Cogsci 8F Linear Algebra review Vectors
More informationQuaternions. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Quaternions Semester 1, / 40
Quaternions Basilio Bona DAUIN Politecnico di Torino Semester 1, 2016-2017 B. Bona (DAUIN) Quaternions Semester 1, 2016-2017 1 / 40 Introduction Complex numbers with unit norm can be used as rotation operators
More informationVectors. September 2, 2015
Vectors September 2, 2015 Our basic notion of a vector is as a displacement, directed from one point of Euclidean space to another, and therefore having direction and magnitude. We will write vectors in
More information:25 1. Rotations. A rotation is in general a transformation of the 3D space with the following properties:
2011-02-17 12:25 1 1 Rotations in general Rotations A rotation is in general a transformation of the 3D space with the following properties: 1. does not change the distances between positions 2. does not
More informationSemester University of Sheffield. School of Mathematics and Statistics
University of Sheffield School of Mathematics and Statistics MAS140: Mathematics (Chemical) MAS152: Civil Engineering Mathematics MAS152: Essential Mathematical Skills & Techniques MAS156: Mathematics
More informationLecture 7. Quaternions
Matthew T. Mason Mechanics of Manipulation Spring 2012 Today s outline Motivation Motivation have nice geometrical interpretation. have advantages in representing rotation. are cool. Even if you don t
More informationsin(α + θ) = sin α cos θ + cos α sin θ cos(α + θ) = cos α cos θ sin α sin θ
Rotations in the 2D Plane Trigonometric addition formulas: sin(α + θ) = sin α cos θ + cos α sin θ cos(α + θ) = cos α cos θ sin α sin θ Rotate coordinates by angle θ: 1. Start with x = r cos α y = r sin
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 informationUnderstanding Quaternions: Rotations, Reflections, and Perspective Projections. Ron Goldman Department of Computer Science Rice University
Understanding Quaternions: Rotations, Reflections, and Perspective Projections Ron Goldman Department of Computer Science Rice University The invention of the calculus of quaternions is a step towards
More informationMatrix Representation
Matrix Representation Matrix Rep. Same basics as introduced already. Convenient method of working with vectors. Superposition Complete set of vectors can be used to express any other vector. Complete set
More informationGroup, Rings, and Fields Rahul Pandharipande. I. Sets Let S be a set. The Cartesian product S S is the set of ordered pairs of elements of S,
Group, Rings, and Fields Rahul Pandharipande I. Sets Let S be a set. The Cartesian product S S is the set of ordered pairs of elements of S, A binary operation φ is a function, S S = {(x, y) x, y S}. φ
More informationEXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. 1. Determinants
EXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. Determinants Ex... Let A = 0 4 4 2 0 and B = 0 3 0. (a) Compute 0 0 0 0 A. (b) Compute det(2a 2 B), det(4a + B), det(2(a 3 B 2 )). 0 t Ex..2. For
More informationLinear Algebra Notes. Lecture Notes, University of Toronto, Fall 2016
Linear Algebra Notes Lecture Notes, University of Toronto, Fall 2016 (Ctd ) 11 Isomorphisms 1 Linear maps Definition 11 An invertible linear map T : V W is called a linear isomorphism from V to W Etymology:
More information5. Vector Algebra and Spherical Trigonometry (continued)
MA232A Euclidean and Non-Euclidean Geometry School of Mathematics, Trinity College Michaelmas Term 2017 Section 5: Vector Algebra and Spherical Trigonometry David R. Wilkins 5.1. Vectors in Three-Dimensional
More informationMathematical Foundations: Intro
Mathematical Foundations: Intro Graphics relies on 3 basic objects: 1. Scalars 2. Vectors 3. Points Mathematically defined in terms of spaces: 1. Vector space 2. Affine space 3. Euclidean space Math required:
More informationBSc (Hons) in Computer Games Development. vi Calculate the components a, b and c of a non-zero vector that is orthogonal to
1 APPLIED MATHEMATICS INSTRUCTIONS Full marks will be awarded for the correct solutions to ANY FIVE QUESTIONS. This paper will be marked out of a TOTAL MAXIMUM MARK OF 100. Credit will be given for clearly
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 informationVector Calculus Æ Z 4 { 9
Vector Calculus 2008ôÆ 9 4ß ii 1 Much of the material presented here and some figures are copied from Vector Calculus of J. Marsden and A. Tromba for teaching purpose only. Please do not distribute This
More informationRigid Geometric Transformations
Rigid Geometric Transformations Carlo Tomasi This note is a quick refresher of the geometry of rigid transformations in three-dimensional space, expressed in Cartesian coordinates. 1 Cartesian Coordinates
More informationRotations. sin θ cos θ 0. a a 2
Rotations Ari Grant Sugih Jamin EECS 487 Winter 010 1. rotation of a vector around an arbitrary axis Given a vector b, an axis a, and an angle α the goal is to find a formula (or function or operator)
More informationLecture 1 Complex Numbers. 1 The field of complex numbers. 1.1 Arithmetic operations. 1.2 Field structure of C. MATH-GA Complex Variables
Lecture Complex Numbers MATH-GA 245.00 Complex Variables The field of complex numbers. Arithmetic operations The field C of complex numbers is obtained by adjoining the imaginary unit i to the field R
More informationContents. D. R. Wilkins. Copyright c David R. Wilkins
MA232A Euclidean and Non-Euclidean Geometry School of Mathematics, Trinity College Michaelmas Term 2017 Section 5: Vector Algebra and Spherical Trigonometry D. R. Wilkins Copyright c David R. Wilkins 2015
More informationFurther Mathematical Methods (Linear Algebra) 2002
Further Mathematical Methods (Linear Algebra) 22 Solutions For Problem Sheet 3 In this Problem Sheet, we looked at some problems on real inner product spaces. In particular, we saw that many different
More informationSINGLE MATHEMATICS B : Vectors Summary Notes
Preprint typeset in JHEP style - HYPER VERSION SINGLE MATHEMATICS B : Vectors Summary Notes Ruth Gregory Abstract: These notes sum up all you need to know about the mathematics of vectors at this stage.
More informationDual unitary matrices and unit dual quaternions
Dual unitary matrices and unit dual quaternions Erhan Ata and Yusuf Yayli Abstract. In this study, the dual complex numbers defined as the dual quaternions have been considered as a generalization of complex
More informationThe Matrix Representation of a Three-Dimensional Rotation Revisited
Physics 116A Winter 2010 The Matrix Representation of a Three-Dimensional Rotation Revisited In a handout entitled The Matrix Representation of a Three-Dimensional Rotation, I provided a derivation of
More informationChapter 8. Rigid transformations
Chapter 8. Rigid transformations We are about to start drawing figures in 3D. There are no built-in routines for this purpose in PostScript, and we shall have to start more or less from scratch in extending
More informationMATH 135: COMPLEX NUMBERS
MATH 135: COMPLEX NUMBERS (WINTER, 010) The complex numbers C are important in just about every branch of mathematics. These notes 1 present some basic facts about them. 1. The Complex Plane A complex
More informationCalculus Vector Principia Mathematica. Lynne Ryan Associate Professor Mathematics Blue Ridge Community College
Calculus Vector Principia Mathematica Lynne Ryan Associate Professor Mathematics Blue Ridge Community College Defining a vector Vectors in the plane A scalar is a quantity that can be represented by a
More informationSPLIT QUATERNIONS and CANAL SURFACES. in MINKOWSKI 3-SPACE
INTERNATIONAL JOURNAL OF GEOMETRY Vol. 5 (016, No., 51-61 SPLIT QUATERNIONS and CANAL SURFACES in MINKOWSKI 3-SPACE SELAHATTIN ASLAN and YUSUF YAYLI Abstract. A canal surface is the envelope of a one-parameter
More informationMatrices A brief introduction
Matrices A brief introduction Basilio Bona DAUIN Politecnico di Torino Semester 1, 2014-15 B. Bona (DAUIN) Matrices Semester 1, 2014-15 1 / 41 Definitions Definition A matrix is a set of N real or complex
More information5.6. PSEUDOINVERSES 101. A H w.
5.6. PSEUDOINVERSES 0 Corollary 5.6.4. If A is a matrix such that A H A is invertible, then the least-squares solution to Av = w is v = A H A ) A H w. The matrix A H A ) A H is the left inverse of A and
More informationMA232A Euclidean and Non-Euclidean Geometry School of Mathematics, Trinity College Michaelmas Term 2017 Vector Algebra and Spherical Trigonometry
MA232A Euclidean and Non-Euclidean Geometry School of Mathematics, Trinity College Michaelmas Term 2017 Vector Algebra and Spherical Trigonometry David R. Wilkins 3.1. Vectors in Three-Dimensional Euclidean
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 informationRigid Geometric Transformations
Rigid Geometric Transformations Carlo Tomasi This note is a quick refresher of the geometry of rigid transformations in three-dimensional space, expressed in Cartesian coordinates. 1 Cartesian Coordinates
More informationMath 489AB Exercises for Chapter 1 Fall Section 1.0
Math 489AB Exercises for Chapter 1 Fall 2008 Section 1.0 1.0.2 We want to maximize x T Ax subject to the condition x T x = 1. We use the method of Lagrange multipliers. Let f(x) = x T Ax and g(x) = x T
More informationIntroduction to Lie Groups
Introduction to Lie Groups MAT 4144/5158 Winter 2015 Alistair Savage Department of Mathematics and Statistics University of Ottawa This work is licensed under a Creative Commons Attribution-ShareAlike
More informationLie Algebra of Unit Tangent Bundle in Minkowski 3-Space
INTERNATIONAL ELECTRONIC JOURNAL OF GEOMETRY VOLUME 12 NO. 1 PAGE 1 (2019) Lie Algebra of Unit Tangent Bundle in Minkowski 3-Space Murat Bekar (Communicated by Levent Kula ) ABSTRACT In this paper, a one-to-one
More informationRobotics, Geometry and Control - Rigid body motion and geometry
Robotics, Geometry and Control - Rigid body motion and geometry Ravi Banavar 1 1 Systems and Control Engineering IIT Bombay HYCON-EECI Graduate School - Spring 2008 The material for these slides is largely
More informationExam 1 Review SOLUTIONS
1. True or False (and give a short reason): Exam 1 Review SOLUTIONS (a) If the parametric curve x = f(t), y = g(t) satisfies g (1) = 0, then it has a horizontal tangent line when t = 1. FALSE: To make
More informationKnowledge Discovery and Data Mining 1 (VO) ( )
Knowledge Discovery and Data Mining 1 (VO) (707.003) Review of Linear Algebra Denis Helic KTI, TU Graz Oct 9, 2014 Denis Helic (KTI, TU Graz) KDDM1 Oct 9, 2014 1 / 74 Big picture: KDDM Probability Theory
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 informationPosition and orientation of rigid bodies
Robotics 1 Position and orientation of rigid bodies Prof. Alessandro De Luca Robotics 1 1 Position and orientation right-handed orthogonal Reference Frames RF A A z A p AB B RF B z B x B y A rigid body
More informationNOTES on LINEAR ALGEBRA 1
School of Economics, Management and Statistics University of Bologna Academic Year 207/8 NOTES on LINEAR ALGEBRA for the students of Stats and Maths This is a modified version of the notes by Prof Laura
More information88 CHAPTER 3. SYMMETRIES
88 CHAPTER 3 SYMMETRIES 31 Linear Algebra Start with a field F (this will be the field of scalars) Definition: A vector space over F is a set V with a vector addition and scalar multiplication ( scalars
More informationRidig Body Motion Homogeneous Transformations
Ridig Body Motion Homogeneous Transformations Claudio Melchiorri Dipartimento di Elettronica, Informatica e Sistemistica (DEIS) Università di Bologna email: claudio.melchiorri@unibo.it C. Melchiorri (DEIS)
More information3. Interpret the graph of x = 1 in the contexts of (a) a number line (b) 2-space (c) 3-space
MA2: Prepared by Dr. Archara Pacheenburawana Exercise Chapter 3 Exercise 3.. A cube of side 4 has its geometric center at the origin and its faces parallel to the coordinate planes. Sketch the cube and
More informationSpecial Lecture - The Octionions
Special Lecture - The Octionions March 15, 2013 1 R 1.1 Definition Not much needs to be said here. From the God given natural numbers, we algebraically build Z and Q. Then create a topology from the distance
More informationREVIEW - Vectors. Vectors. Vector Algebra. Multiplication by a scalar
J. Peraire Dynamics 16.07 Fall 2004 Version 1.1 REVIEW - Vectors By using vectors and defining appropriate operations between them, physical laws can often be written in a simple form. Since we will making
More informationQuaternion kinematics for the error-state Kalman filter
Quaternion kinematics for the error-state Kalman filter Joan Solà October 12, 2017 Contents 1 Quaternion definition and properties 4 1.1 Definition of quaternion............................ 4 1.1.1 Alternative
More informationChapter 8 Vectors and Scalars
Chapter 8 193 Vectors and Scalars Chapter 8 Vectors and Scalars 8.1 Introduction: In this chapter we shall use the ideas of the plane to develop a new mathematical concept, vector. If you have studied
More informationPart IA. Vectors and Matrices. Year
Part IA Vectors and Matrices Year 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2018 Paper 1, Section I 1C Vectors and Matrices For z, w C define the principal value of z w. State de Moivre s
More informationQuantum Computing Lecture 2. Review of Linear Algebra
Quantum Computing Lecture 2 Review of Linear Algebra Maris Ozols Linear algebra States of a quantum system form a vector space and their transformations are described by linear operators Vector spaces
More informationp 2 p 3 p y p z It will not be considered in the present context; the interested reader can find more details in [05].
1. Geometrical vectors A geometrical vector p represents a point P in space. The point P is an abstraction that often, but not always, requires a representation. Vector representations are given wrt a
More informationREPRESENTATION THEORY WEEK 7
REPRESENTATION THEORY WEEK 7 1. Characters of L k and S n A character of an irreducible representation of L k is a polynomial function constant on every conjugacy class. Since the set of diagonalizable
More informationGeometry of the Special Unitary Group
Geometry of the Special Unitary Group The elements of SU 2 are the unitary 2 2 matrices with determinant 1. It is not hard to see that they have the form a 1) ) b, b ā with āa+bb = 1. This is the transpose
More information6 Inner Product Spaces
Lectures 16,17,18 6 Inner Product Spaces 6.1 Basic Definition Parallelogram law, the ability to measure angle between two vectors and in particular, the concept of perpendicularity make the euclidean space
More informationJim Lambers MAT 610 Summer Session Lecture 1 Notes
Jim Lambers MAT 60 Summer Session 2009-0 Lecture Notes Introduction This course is about numerical linear algebra, which is the study of the approximate solution of fundamental problems from linear algebra
More informationIntroduction to Geometry
Introduction to Geometry it is a draft of lecture notes of H.M. Khudaverdian. Manchester, 18 May 211 Contents 1 Euclidean space 3 1.1 Vector space............................ 3 1.2 Basic example of n-dimensional
More informationMAT 1339-S14 Class 8
MAT 1339-S14 Class 8 July 28, 2014 Contents 7.2 Review Dot Product........................... 2 7.3 Applications of the Dot Product..................... 4 7.4 Vectors in Three-Space.........................
More informationIntroduction to Algebraic and Geometric Topology Week 14
Introduction to Algebraic and Geometric Topology Week 14 Domingo Toledo University of Utah Fall 2016 Computations in coordinates I Recall smooth surface S = {f (x, y, z) =0} R 3, I rf 6= 0 on S, I Chart
More informationClifford Analysis, Homework 1
Clifford Analysis, Homework 1 November 1, 017 1 Let w v 1... v k, for vectors v j. Show that ŵ is the result of negating the vectors: ŵ ( v 1 )... ( v k ). Show that w is the result of changing the order
More informationLinear Algebra: Lecture notes from Kolman and Hill 9th edition.
Linear Algebra: Lecture notes from Kolman and Hill 9th edition Taylan Şengül March 20, 2019 Please let me know of any mistakes in these notes Contents Week 1 1 11 Systems of Linear Equations 1 12 Matrices
More informationChapter y. 8. n cd (x y) 14. (2a b) 15. (a) 3(x 2y) = 3x 3(2y) = 3x 6y. 16. (a)
Chapter 6 Chapter 6 opener A. B. C. D. 6 E. 5 F. 8 G. H. I. J.. 7. 8 5. 6 6. 7. y 8. n 9. w z. 5cd.. xy z 5r s t. (x y). (a b) 5. (a) (x y) = x (y) = x 6y x 6y = x (y) = (x y) 6. (a) a (5 a+ b) = a (5
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 information2) If a=<2,-1> and b=<3,2>, what is a b and what is the angle between the vectors?
CMCS427 Dot product review Computing the dot product The dot product can be computed via a) Cosine rule a b = a b cos q b) Coordinate-wise a b = ax * bx + ay * by 1) If a b, a and b all equal 1, what s
More informationLecture 13: Matrix Representations for Affine and Projective Transformations
Lecture 13: Matrix Representations for Affine and Projective Transformations thou shalt set apart unto the LORD all that openeth the matrix Exodus 13:12 1 Matrix Representations for Affine Transformations
More informationClosed-Form Solution Of Absolute Orientation Using Unit Quaternions
Closed-Form Solution Of Absolute Orientation Using Unit Berthold K. P. Horn Department of Computer and Information Sciences November 11, 2004 Outline 1 Introduction 2 3 The Problem Given: two sets of corresponding
More information1.1 Single Variable Calculus versus Multivariable Calculus Rectangular Coordinate Systems... 4
MATH2202 Notebook 1 Fall 2015/2016 prepared by Professor Jenny Baglivo Contents 1 MATH2202 Notebook 1 3 1.1 Single Variable Calculus versus Multivariable Calculus................... 3 1.2 Rectangular Coordinate
More informationLecture 5. Complex Numbers and Euler s Formula
Lecture 5. Complex Numbers and Euler s Formula University of British Columbia, Vancouver Yue-Xian Li March 017 1 Main purpose: To introduce some basic knowledge of complex numbers to students so that they
More informationMCE/EEC 647/747: Robot Dynamics and Control. Lecture 2: Rigid Motions and Homogeneous Transformations
MCE/EEC 647/747: Robot Dynamics and Control Lecture 2: Rigid Motions and Homogeneous Transformations Reading: SHV Chapter 2 Mechanical Engineering Hanz Richter, PhD MCE503 p.1/22 Representing Points, Vectors
More informationCOMPLEX NUMBERS AND QUADRATIC EQUATIONS
Chapter 5 COMPLEX NUMBERS AND QUADRATIC EQUATIONS 5. Overview We know that the square of a real number is always non-negative e.g. (4) 6 and ( 4) 6. Therefore, square root of 6 is ± 4. What about the square
More informationInterpolation of Rigid Motions in 3D
Interpolation of Rigid Motions in 3D David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License.
More information3 + 4i 2 + 3i. 3 4i Fig 1b
The introduction of complex numbers in the 16th century was a natural step in a sequence of extensions of the positive integers, starting with the introduction of negative numbers (to solve equations of
More informationELA
LEFT EIGENVALUES OF 2 2 SYMPLECTIC MATRICES E. MACíAS-VIRGÓS AND M.J. PEREIRA-SÁEZ Abstract. A complete characterization is obtained of the 2 2 symplectic matrices that have an infinite number of left
More informationMatrix Theory and Differential Equations Homework 6 Solutions, 10/5/6
Matrix Theory and Differential Equations Homework 6 Solutions, 0/5/6 Question Find the general solution of the matrix system: x 3y + 5z 8t 5 x + 4y z + t Express your answer in the form of a particulaolution
More information8 Velocity Kinematics
8 Velocity Kinematics Velocity analysis of a robot is divided into forward and inverse velocity kinematics. Having the time rate of joint variables and determination of the Cartesian velocity of end-effector
More informationON A NEW SPECIES OF IMAGINARY QUANTITIES CONNECTED WITH A THEORY OF QUATERNIONS. William Rowan Hamilton
ON A NEW SPECIES OF IMAGINARY QUANTITIES CONNECTED WITH A THEORY OF QUATERNIONS By William Rowan Hamilton (Proceedings of the Royal Irish Academy, 2 (1844), 424 434.) Edited by David R. Wilkins 1999 On
More informationMathematics for Graphics and Vision
Mathematics for Graphics and Vision Steven Mills March 3, 06 Contents Introduction 5 Scalars 6. Visualising Scalars........................ 6. Operations on Scalars...................... 6.3 A Note on
More informationComplex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,
Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )
More informationA Review of Linear Algebra
A Review of Linear Algebra Gerald Recktenwald Portland State University Mechanical Engineering Department gerry@me.pdx.edu These slides are a supplement to the book Numerical Methods with Matlab: Implementations
More informationQuaternion and Rotation
Quaternion and Rotation Com S 477/577 Notes) Yan-Bin Jia Sep 5, 017 1 Introduction Up until now we have learned that a rotation in R 3 about an axis through the origin can be represented by a 3 3 orthogonal
More informationUnit Generalized Quaternions in Spatial Kinematics
Vol (07) - Pages 7-4 Unit Generalized Quaternions in Spatial Kinematics Mehdi Jafari, echnical Vocational University, Urmia, Iran, mj_msc@yahoo.com Abstract- After a review of some fundamental properties
More informationEquality: Two matrices A and B are equal, i.e., A = B if A and B have the same order and the entries of A and B are the same.
Introduction Matrix Operations Matrix: An m n matrix A is an m-by-n array of scalars from a field (for example real numbers) of the form a a a n a a a n A a m a m a mn The order (or size) of A is m n (read
More informationMath 3108: Linear Algebra
Math 3108: Linear Algebra Instructor: Jason Murphy Department of Mathematics and Statistics Missouri University of Science and Technology 1 / 323 Contents. Chapter 1. Slides 3 70 Chapter 2. Slides 71 118
More informationA (Mostly) Linear Algebraic Introduction to Quaternions
A (Mostly) Linear Algebraic Introduction to Quaternions Joe McMahon Program in Applied Mathematics University of Arizona Fall 23 1 Some History 1.1 Hamilton s Discovery and Subsequent Vandalism Having
More information