Quaternions. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Quaternions Semester 1, / 40
|
|
- Clarissa Barker
- 5 years ago
- Views:
Transcription
1 Quaternions Basilio Bona DAUIN Politecnico di Torino Semester 1, B. Bona (DAUIN) Quaternions Semester 1, / 40
2 Introduction Complex numbers with unit norm can be used as rotation operators in the 2D plane. Quaternions were discovered in 1843 by Sir William Hamilton, who was trying to find new hyper-complex numbers as rotation operators in 3D space. Figure: Sir William Rowland Hamilton ( ) and the plaque on Broom Bridge, where he discovered the quaternions. B. Bona (DAUIN) Quaternions Semester 1, / 40
3 Definitions The generic quaternion will be indicated as q. Quaternions are elements of a 4D linear space H(R), defined on the real numbers field F = R, with base {1 i j k}. i, j and k are hypercomplex numbers that satisfy the following anticommutative multiplication rules: i 2 = j 2 = k 2 = ijk = 1 ij = ji = k jk = kj = i ki = ik = j Notice the similarity of quaternions i, j and k with the basis vectors i, j and k B. Bona (DAUIN) Quaternions Semester 1, / 40
4 Definitions A quaternion q H is defined as a linear combination expressed in the base {1 i j k}: q = q q 1 i + q 2 j + q 3 k where the coefficients {q i } 3 i=0 are real. Another way to represent a quaternion is to define it as a quadruple of reals, (q 0,q 1,q 2,q 3 ) q = (q 0,q 1,q 2,q 3 ) in analogy with complex numbers c = a + jb, where c is represented by a couple of reals, c = (a,b), B. Bona (DAUIN) Quaternions Semester 1, / 40
5 Definitions Quaternions are also defined as hypercomplex numbers, i.e., those complex numbers having complex coefficients: where c 1 = q 0 + kq 3 e c 2 = q 2 + kq 1. q = c 1 + jc 2, Therefore, considering multiplication rules, it results: q = c 1 + jc 2 = q 0 + kq 3 + jq 2 + jkq 1 = q q 1 i + q 2 j + q 3 k B. Bona (DAUIN) Quaternions Semester 1, / 40
6 Definitions In analogy with complex numbers that are the sum of a real part and an imaginary part, quaternions are the sum of a real part and a vectorial part. The real part q r is defined as q r = q 0, and the vectorial part q v is defined as q v = q 1 i + q 2 j + q 3 k. We write q = (q r, q v ) or q = q r + q v ; note that the vectorial part is not transposed since the conventional definition for the vectorial part of a quaternion assumes a row representation. Using the convention that defines vectors as column vectors, we can write q = (q r, q T v ). B. Bona (DAUIN) Quaternions Semester 1, / 40
7 Definitions Quaternions are mathematical objects that include real numbers r = (r, 0, 0, 0), r R complex numbers a + ib = (a, b, 0, 0), a,b R real vectors in R 3 (with some caution, since not all vectorial parts represents vectors) v = (0, v 1, v 2, v 3 ), v i R. In this last case, elements { i j k } are to be understood as unit vectors { i j k } forming an orthonormal base in a cartesian right-handed reference frame. B. Bona (DAUIN) Quaternions Semester 1, / 40
8 Definitions Multiplication rules between elements i, j, k have the same properties of the cross product between unit vectors i,j,k: ij = k i j = k ji = k j i = k etc. In the following we will use all the possible alternative notations to indicate quaternions q = q q 1 i + q 2 j + q 3 k = (q r,q v ) = q r + q v = (q 0,q 1,q 2,q 3 ) i.e., a) a hypercomplex number; b) the sum of a real part and a vectorial part; c) a quadruple of reals. B. Bona (DAUIN) Quaternions Semester 1, / 40
9 Definitions Another way to write a quaternion is using 2 2 complex matrices q = q q 1 i + q 2 j + q 3 k where 1 = [ ] [ ] 1 0 i 0 ; i = ; j = i [ ] 0 1 ; k = 1 0 [ ] 0 i ; i 0 and i 2 = 1. Hence every matrix is of the form [ ] c d d c ; These matrices are called Cayley matrices [see slide on Rotations - Pauli spin matrices]. B. Bona (DAUIN) Quaternions Semester 1, / 40
10 Quaternions algebra Given a quaternion q = q q 1 i + q 2 j + q 3 k = q r + q v = (q 0,q 1,q 2,q 3 ), the following properties hold: null or zero 0 quaternion 0 = i + 0j + 0k = (0,0) = = (0,0,0,0) conjugate quaternion q, having the same real part and the opposite vectorial part: q = q 0 (q 1 i + q 2 j + q 3 k) = (q r, q v ) = q r q v = (q 0, q 1, q 2, q 3 ) Conjugate quaternion satisfies (q ) = q B. Bona (DAUIN) Quaternions Semester 1, / 40
11 Quaternions algebra non-negative function, called quaternion norm q, defined as q 2 = qq = q q = 3 l=0 q 2 l = q q T v q v A quaternion with unit norm q = 1 is called unit quaternion. Quaternion q and its conjugate q have the same norm The quaternion q = q q v = 01 + q 1 i + q 2 j + q 3 k = (0,q v ) = 0 + q v = (0,q 1,q 2,q 3 ), that has a zero real part is called pure quaternion or vector. The conjugate of a pure quaternion q v is the opposite of the original pure quaternion q v = q v B. Bona (DAUIN) Quaternions Semester 1, / 40
12 Quaternions algebra Given two quaternions h = h h 1 i + h 2 j + h 3 k = (h r,h v ) = h r + h v = (h 0,h 1,h 2,h 3 ) and g = g g 1 i + g 2 j + g 3 k = (g r,g v ) = g r + g v = (g 0,g 1,g 2,g 3 ) the following operations are defined B. Bona (DAUIN) Quaternions Semester 1, / 40
13 Sum Sum or addition h + g h + g = (h 0 + g 0 )1 + (h 1 + g 1 )i + (h 2 + g 2 )j + (h 3 + g 3 )k = ((h r + g r ), (h v + g v )) = (h r + g r ) + (h v + g v ) = (h 0 + g 0,h 1 + g 1,h 2 + g 2,h 3 + g 3 ) B. Bona (DAUIN) Quaternions Semester 1, / 40
14 Difference Difference or subtraction h g = (h 0 g 0 )1 + (h 1 g 1 )i + (h 2 g 2 )j + (h 3 g 3 )k = ((h r g r ), (h v g v )) = (h r g r ) + (h v g v ) = (h 0 g 0,h 1 g 1,h 2 g 2,h 3 g 3 ) B. Bona (DAUIN) Quaternions Semester 1, / 40
15 Product Product hg = (h 0 g 0 h 1 g 1 h 2 g 2 h 3 g 3 )1 + (h 1 g 0 + h 0 g 1 h 3 g 2 + h 2 g 3 )i + (h 2 g 0 + h 3 g 1 + h 0 g 2 h 1 g 3 )j + (h 3 g 0 h 2 g 1 + h 1 g 2 + h 0 g 3 )k = (h r g r h v g v, h r g v + g r h v + h v g v ) where h v g v is the scalar product h v g v = h vi g vi = h T v g v = g T v h v i defined in R n, and h v g v is the cross product (defined only in R 3 ) h v g v = h 2g 3 h 3 g 2 h 3 g 1 h 1 g 3 = S(h v )g v h 1 g 2 h 2 g 1 B. Bona (DAUIN) Quaternions Semester 1, / 40
16 Product The quaternion product is anti-commutative, since, being it follows that g v h v = h v g v gh = (h r g r h v g v, h r g v + g r h v h v g v ) hg Notice that the real part remains the same, while the vectorial part changes. Product commutes only if h v g v = 0, i.e., when the vectorial parts are parallel. The conjugate of a quaternion product satisfies this relation (gh) = h g The product norm satisfies this relation gh = hg = g h = h g B. Bona (DAUIN) Quaternions Semester 1, / 40
17 Product properties The quaternion product has the following properties associative multiplication by the unit scalar (gh)p = g(hp) 1q = q1 = (1,0)(q r,q v ) = (1q r,1q v ) = (q r,q v ) multiplication by the real λ bilinearity, with real λ 1,λ 2 λq = (λ,0)(q r,q v ) = (λq r,λq v ) g(λ 1 h 1 + λ 2 h 2 ) = λ 1 gh 1 + λ 2 gh 2 (λ 1 g 1 + λ 2 g 2 )h = λ 1 g 1 h + λ 2 g 2 h B. Bona (DAUIN) Quaternions Semester 1, / 40
18 Product Alternative forms Quaternion product may be written as matrix product forms: h 0 h 1 h 2 h 3 g 0 [ ] g 0 h hg = 1 h 0 h 3 h 2 g 1 h0 h h 2 h 3 h 0 h 1 g = T v g 1 2 h v h 0 I + S(h v ) g 2 h 3 h 2 h 1 h 0 g 3 g 3 = F L (h)g or, alternatively, as g 0 g 1 g 2 g 3 h 0 [ ] h 0 g hg = 1 g 0 g 3 g 2 h 1 g0 g g 2 g 3 g 0 g 1 h = T v h 1 2 g v g 0 I S(g v ) h 2 g 3 g 2 g 1 g 0 h 3 h 3 = F R (g)h B. Bona (DAUIN) Quaternions Semester 1, / 40
19 Quotient Since the quaternion product is anti-commutative we must distinguish between the left and the right quotient or division. Given two quaternions h e p, we define the left quotient of p by h the quaternion q l that satisfies hq l = p while we define the right quotient of p by h the quaternion q r that satisfies q r h = p Hence q l = h h 2 p; q r = p h h 2 B. Bona (DAUIN) Quaternions Semester 1, / 40
20 Inverse Given a quaternion q, in principle one must define the right q 1 r left inverse q 1 l as and the qq 1 l = 1 = (1,0,0,0); q 1 r q = 1 = (1,0,0,0) Since qq = q q = q 2 = q q, one can write q q q = q q q q = 1 = (1,0,0,0) q It follows that right inverse and left inverse are equal q 1 l = q 1 r = q 1 = q q 2 similar to c 1 = c c 2 B. Bona (DAUIN) Quaternions Semester 1, / 40
21 Inverse For a unit quaternion u, u = 1, inverse and conjugate coincide u 1 = u, u = 1 and for a pure unit quaternion q = (0,q v ), q = 1, i.e., a unit vector Inverse satisfies q 1 v = q v = q v. (q 1 ) 1 = q; (pq) 1 = q 1 p 1 B. Bona (DAUIN) Quaternions Semester 1, / 40
22 Unit quaternions Before starting to illustrate the relations between quaternions and rotations, we look closer to the properties of the unit quaternions, that we indicate with the symbol u. A unit quaternion u with u = 1, has an unit inverse and the product of two unit quaternions is still a unit quaternion. We assume that a unit quaternion is represented by a sum of two trigonometric functions u = cosθ + usinθ = (cosθ,usinθ) where u is a unit norm vector and θ a generic angle. B. Bona (DAUIN) Quaternions Semester 1, / 40
23 Unit quaternions Notice the analogy of the formula above with the unit complex number expression c = cosθ + jsinθ The analogy applies also to the exponential expression c = e jθ ; substituting x with uθ in the series expansion of e x and recalling that uu = 1, we have e uθ = cosθ + usinθ = u The relation above shows a formal identity between a unit quaternion and the exponential of a unit vector multiplied by a scalar θ Notice the similarity between u = e uθ and R(u,θ) = e S(u)θ B. Bona (DAUIN) Quaternions Semester 1, / 40
24 Unit quaternions From the previous relation one obtains the power p of a unit quaternion as u p = (cosθ + usinθ) p = e uθp = cos(θp) + usin(θp) and the logarithm of a unit quaternion logu = log(cosθ + usinθ) = log(e uθ ) = uθ Notice that the anti-commutativity of the quaternion product inhibits to use the standard identities between exponential and logarithms. For instance, e u 1 e u 2 it is not necessarily equal to e u 1+u 2, and log(u 1 u 2 ) is not necessarily equal to log(u 1 ) + log(u 2 ). B. Bona (DAUIN) Quaternions Semester 1, / 40
25 Quaternions and Rotations Now we relate rotations and unit quaternions: the unit quaternion u = (u 0,u 1,u 2,u 3 ) = (u 0,u) = cosθ + usinθ represents the rotation Rot(u,2θ) around the axis specified by the unit vector u = [ u 1 u 2 u 3 ] T The converse is also true, i.e., given a rotation Rot(u,θ) of an angle θ around the axis specified by the unit vector u = [ u 1 u 2 u 3 ] T, the unit quaternion ( u = cos θ 2,u 1 sin θ 2,u 2 sin θ 2,u 3 sin θ ) ( = cos θ 2 2, usin θ ) = 2 represents the same rotations. cos θ 2 + usin θ 2 B. Bona (DAUIN) Quaternions Semester 1, / 40
26 Quaternions and Rotations We know that a rigid rotation in R 3 is represented by a rotation (orthonormal) matrix R SO(3) R 3 3. We can associate to every rotation matrix R a unit quaternion u and viceversa, indicating this relation as R(u) u. B. Bona (DAUIN) Quaternions Semester 1, / 40
27 Quaternions and Rotations To compute the rotation matrix R(u) given a unit quaternion u = (u 0,u), we use the following relation R(u) = (u0 2 u T u)i + 2uu T 2u 0 S(u) = u0 2 + u2 1 u2 2 u2 3 2(u 1 u 2 u 3 u 0 ) 2(u 1 u 3 + u 2 u 0 ) 2(u 1 u 2 + u 3 u 0 ) u0 2 u2 1 + u2 2 u2 3 2(u 2 u 3 u 1 u 0 ) 2(u 1 u 3 u 2 u 0 ) 2(u 2 u 3 + u 1 u 0 ) u 2 0 u2 1 u2 2 + u2 3 where S(u) is an skew-symmetric matrix. B. Bona (DAUIN) Quaternions Semester 1, / 40
28 Quaternions and Rotations To compute the unit quaternion u = (u 0,u) given a rotation matrix R we use the following relations u 0 = ± 1 2 (1 + r11 + r 22 + r 33 ) u 1 = 1 4u 0 (r 32 r 23 ) u 2 = 1 4u 0 (r 13 r 31 ) u 3 = 1 4u 0 (r 21 r 12 ) B. Bona (DAUIN) Quaternions Semester 1, / 40
29 Quaternions and Rotations Another alternative set of relations is u 0 = 1 2 (1 + r11 + r 22 + r 33 ) u 1 = 1 2 sgn(r 32 r 23 ) (1 + r 11 r 22 r 33 ) u 2 = 1 2 sgn(r 13 r 31 ) (1 r 11 + r 22 r 33 ) u 3 = 1 2 sgn(r 21 r 12 ) (1 r 11 r 22 + r 33 ) where sgn(x) is the sign function, with sgn(0) = 0. This relation is never singular compared with the previous one that is singular for u 0 = 0 B. Bona (DAUIN) Quaternions Semester 1, / 40
30 Quaternions and Rotations Elementary rotations around the three principal axes R(i, α), R(j, β) and R(k, γ), correspond to the following elementary quaternions R(i,α) u x = (cos α 2, sin α ) 2 (cos, 0, 0 β2, 0, sin β2 ), 0 R(j,β) u y = R(k,γ) u z = ( cos γ 2, 0, 0, sin γ ) 2 It is easy to see that the vectorial base of quaternions correspond to elementary rotations of π around the principal axes i = (0,1,0,0) R(i,π) j = (0,0,1,0) R(j,π) k = (0,0,0,1) R(k,π) B. Bona (DAUIN) Quaternions Semester 1, / 40
31 Quaternions and Rotations Notice an important fact: while the product among the unit base quaternions gives ii = jj = kk = ijk = ( 1, 0, 0, 0), the product among the associated rotation matrices gives: R(i,π)R(i,π) = R(j,π)R(j,π) = R(k,π)R(k,π) = R(i,π)R(j,π)R(k,π) = I Since I represents a rotation that leaves the vectors unchanged, it seems natural to associate it to positive unit scalar, i.e., the quaternion (1,0,0,0). The apparent discrepancy between the two results can be explained only introducing a most basic mathematical quantity, called spinor, not discussed in the present context B. Bona (DAUIN) Quaternions Semester 1, / 40
32 Quaternions and Rotations We shall see now some correspondence between quaternion operations and matrix operations Rotation product Given n rotations R 1, R 2,, R n and the corresponding unit quaternions u 1, u 2,, u n, the product corresponds to the product in the shown order. R(u) = R(u 1 )R(u 2 ) R(u n ) u = u 1 u 2 u n B. Bona (DAUIN) Quaternions Semester 1, / 40
33 Quaternions and Rotations Transpose matrix Given the rotation R(u) and its corresponding unit quaternion u, the transpose matrix (i.e., the inverse rotation) R T corresponds to the conjugate unit quaternion u (i.e., the inverse quaternion) R u R T u B. Bona (DAUIN) Quaternions Semester 1, / 40
34 Quaternions and Rotations Vector rotation Given a generic vector x, and the corresponding pure quaternion q x = (0,x) = (0,x 1,x 2,x 3 ) and given a rotation matrix R(u) with its corresponding unit quaternion u, the rotated vector y = R(u)x is given by the vectorial part of the quaternion obtained as q y = (0,y) = uq x u where q y = q x B. Bona (DAUIN) Quaternions Semester 1, / 40
35 The map q x q y = uq x u y = R(u)x that transforms a pure vector into its rotated counterpart, in quaternion form, is called conjugation by u. Notice that the transpose map is equivalent to exchange the order of the conjugation q x q y = u q x u y = R T (u)x Quaternions u and u = u 1 are called antipodal, because they represent opposite points on the 3-sphere of unit quaternions. B. Bona (DAUIN) Quaternions Semester 1, / 40
36 Quaternions and Rotations If we use the homogeneous coordinates to express the vector x x = [ wx 1 wx 2 wx 3 w ] T and the quaternion x is defined as x = (w,wx) the product uxu provides the quaternion y, defined as y = (w,wr(u)x) The resulting vector y, in homogeneous coordinates is therefore ỹ = [ wy 1 wy 2 wy 3 w ] T y = R(u)x B. Bona (DAUIN) Quaternions Semester 1, / 40
37 Quaternions and Rotations Product matrices Since the bi-linearity property holds for the product between two quaternions, this can be represented by linear operators (i.e., matrices). We recall that the product qp can be expressed in matrix form as qp = F L (q)p This can be interpreted as the left product of q by p. Similarly for the right product pq, expressed as pq = F R (q)p. B. Bona (DAUIN) Quaternions Semester 1, / 40
38 Quaternions and Rotations We can obtain pq as F R (q )p and where the matrix Q R 4 4 is qpq = F L (q)f R (q )p = Qp Q = F L (q)f R (q ) = q0 2 + q2 1 q2 2 q2 3 2(q 1 q 2 q 3 q 0 ) 2(q 1 q 3 + q 2 q 0 ) 0 2(q 1 q 2 + q 3 q 0 ) q0 2 q2 1 + q2 2 q2 3 2(q 2 q 3 q 1 q 0 ) 0 2(q 1 q 3 q 2 q 0 ) 2(q 2 q 3 + q 1 q 0 ) q0 2 q2 1 q2 2 + q q 2 We observe that the upper left 3 3 matrix of Q equals the matrix R(q) previously defined B. Bona (DAUIN) Quaternions Semester 1, / 40
39 Infinitesimal rotations and quaternions Quaternions allow to compute angles directly, if they are small. Assume a small angles rotation Rot(dα) and the associated quaternion h = ( h 0 h 1 h 2 h 3 ). We have R(dα) = I + S(dα) = 1 dα 3 dα 2 dα 3 1 dα 1 = r 11 r 12 r 13 r 21 r 22 r 23 dα 2 dα 1 1 r 31 r 32 r 33 hence, considering h 0 = cos dθ 2 1, dα 1 = 1 2 (r 32 r 23 ) = 2h 1 h 0 2h 1 dα 2 = 1 2 (r 13 r 31 ) = 2h 2 h 0 2h 2 dα 3 = 1 2 (r 21 r 12 ) = 2h 3 h 0 2h 3 B. Bona (DAUIN) Quaternions Semester 1, / 40
40 Quaternions in Aerospace literature In space applications quaternions are used when dealing with satellite orientation control; unfortunately they may be organized in a different way wrt to our conventions: often the real part is the last element of quaternions q = (q 1,q 2,q 3,q 0 ) B. Bona (DAUIN) Quaternions Semester 1, / 40
Geometric Fundamentals in Robotics Quaternions
Geometric Fundamentals in Robotics Quaternions Basilio Bona DAUIN-Politecnico di Torino July 2009 Basilio Bona (DAUIN-Politecnico di Torino) Quaternions July 2009 1 / 47 Introduction Quaternions were discovered
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 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 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 / 44 Definitions Definition A matrix is a set of N real or complex
More informationECS 178 Course Notes QUATERNIONS
ECS 178 Course Notes QUATERNIONS Kenneth I. Joy Institute for Data Analysis and Visualization Department of Computer Science University of California, Davis Overview The quaternion number system was discovered
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 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 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 informationMSMS Vectors and Matrices
MSMS Vectors and Matrices Basilio Bona DAUIN Politecnico di Torino Semester 1, 2015-2016 B. Bona (DAUIN) MSMS-Vectors and matrices Semester 1, 2015-2016 1 / 39 Introduction Most of the topics introduced
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 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 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 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 informationQuaternions. R. J. Renka 11/09/2015. Department of Computer Science & Engineering University of North Texas. R. J.
Quaternions R. J. Renka Department of Computer Science & Engineering University of North Texas 11/09/2015 Definition A quaternion is an element of R 4 with three operations: addition, scalar multiplication,
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 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 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 informationLeandra Vicci. Microelectronic Systems Laboratory. Department of Computer Science. University of North Carolina at Chapel Hill. 27 April 2001.
Quaternions and Rotations in 3-Space: The Algebra and its Geometric Interpretation Leandra Vicci Microelectronic Systems Laboratory Department of Computer Science University of North Carolina at Chapel
More informationQuaternions and their applications
Quaternions and their applications Ramachandran Subramanian February 20, 2014 Jack B. Kuipers. Quaternions and Rotation Sequences. university press, New Jersey, 1998 Princeton Outline 1 Background 2 Introduction
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 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 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 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 informationClifford Algebras and Spin Groups
Clifford Algebras and Spin Groups Math G4344, Spring 2012 We ll now turn from the general theory to examine a specific class class of groups: the orthogonal groups. Recall that O(n, R) is the group of
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 informationFoundations of Matrix Analysis
1 Foundations of Matrix Analysis In this chapter we recall the basic elements of linear algebra which will be employed in the remainder of the text For most of the proofs as well as for the details, the
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 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 informationKinematics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Kinematics Semester 1, / 15
Kinematics Basilio Bona DAUIN Politecnico di Torino Semester 1, 2016-17 B. Bona (DAUIN) Kinematics Semester 1, 2016-17 1 / 15 Introduction The kinematic quantities used to represent a body frame are: position
More informationQuaternions and Octonions
Quaternions and Octonions Alberto Elduque Universidad de Zaragoza 1 Real and complex numbers 2 Quaternions 3 Rotations in three-dimensional space 4 Octonions 2 / 32 1 Real and complex numbers 2 Quaternions
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 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 informationChapter 2 Math Fundamentals
Chapter 2 Math Fundamentals Part 5 2.8 Quaternions 1 Outline 2.8.1 Representations and Notation 2.7.2 Quaternion Multiplication 2.7.3 Other Quaternion Operations 2.7.4 Representing 3D Rotations 2.7.5 Attitude
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 informationMatrices and Vectors
Matrices and Vectors James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 11, 2013 Outline 1 Matrices and Vectors 2 Vector Details 3 Matrix
More informationA FIRST COURSE IN LINEAR ALGEBRA. An Open Text by Ken Kuttler. Matrix Arithmetic
A FIRST COURSE IN LINEAR ALGEBRA An Open Text by Ken Kuttler Matrix Arithmetic Lecture Notes by Karen Seyffarth Adapted by LYRYX SERVICE COURSE SOLUTION Attribution-NonCommercial-ShareAlike (CC BY-NC-SA)
More informationLecture 10: A (Brief) Introduction to Group Theory (See Chapter 3.13 in Boas, 3rd Edition)
Lecture 0: A (Brief) Introduction to Group heory (See Chapter 3.3 in Boas, 3rd Edition) Having gained some new experience with matrices, which provide us with representations of groups, and because symmetries
More informationMatrices A brief introduction
Matrices A brief introduction Basilio Bona DAUIN Politecnico di Torino September 2013 Basilio Bona (DAUIN) Matrices September 2013 1 / 74 Definitions Definition A matrix is a set of N real or complex numbers
More informationWhat is the Matrix? Linear control of finite-dimensional spaces. November 28, 2010
What is the Matrix? Linear control of finite-dimensional spaces. November 28, 2010 Scott Strong sstrong@mines.edu Colorado School of Mines What is the Matrix? p. 1/20 Overview/Keywords/References Advanced
More informationLinear Equations and Matrix
1/60 Chia-Ping Chen Professor Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Gaussian Elimination 2/60 Alpha Go Linear algebra begins with a system of linear
More informationChapter 2. Linear Algebra. rather simple and learning them will eventually allow us to explain the strange results of
Chapter 2 Linear Algebra In this chapter, we study the formal structure that provides the background for quantum mechanics. The basic ideas of the mathematical machinery, linear algebra, are rather simple
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 informationCS 246 Review of Linear Algebra 01/17/19
1 Linear algebra In this section we will discuss vectors and matrices. We denote the (i, j)th entry of a matrix A as A ij, and the ith entry of a vector as v i. 1.1 Vectors and vector operations A vector
More informationRotation outline of talk
Rotation outline of talk Properties Representations Hamilton s Quaternions Rotation as Unit Quaternion The Space of Rotations Photogrammetry Closed Form Solution of Absolute Orientation Division Algebras,
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 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 information1 Matrices and vector spaces
Matrices and vector spaces. Which of the following statements about linear vector spaces are true? Where a statement is false, give a counter-example to demonstrate this. (a) Non-singular N N matrices
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 informationGENERALIZED QUATERNIONS AND ROTATION IN 3-SPACE E 3 αβ
TWMS J. Pure Appl. Math. V.6, N., 015, pp.4-3 GENERALIZED QUATERNIONS AND ROTATION IN 3-SPACE E 3 αβ MEHDI JAFARI 1, YUSUF YAYLI Abstract. The paper explains how a unit generalized quaternion is used to
More informationMath Problem set # 7
Math 128 - Problem set # 7 Clifford algebras. April 8, 2004 due April 22 Because of disruptions due the the Jewish holidays and surgery on my knee there will be no class next week (April 13 or 15). (Doctor
More informationLinear Algebra March 16, 2019
Linear Algebra March 16, 2019 2 Contents 0.1 Notation................................ 4 1 Systems of linear equations, and matrices 5 1.1 Systems of linear equations..................... 5 1.2 Augmented
More informationA Review of Matrix Analysis
Matrix Notation Part Matrix Operations Matrices are simply rectangular arrays of quantities Each quantity in the array is called an element of the matrix and an element can be either a numerical value
More informationKinematics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Kinematics Semester 1, / 15
Kinematics Basilio Bona DAUIN Politecnico di Torino Semester 1, 2014-15 B. Bona (DAUIN) Kinematics Semester 1, 2014-15 1 / 15 Introduction The kinematic quantities used are: position r, linear velocity
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 informationElementary maths for GMT
Elementary maths for GMT Linear Algebra Part 2: Matrices, Elimination and Determinant m n matrices The system of m linear equations in n variables x 1, x 2,, x n a 11 x 1 + a 12 x 2 + + a 1n x n = b 1
More informationOn complexified quantum mechanics and space-time
On complexified quantum mechanics and space-time Dorje C. Brody Mathematical Sciences Brunel University, Uxbridge UB8 3PH dorje.brody@brunel.ac.uk Quantum Physics with Non-Hermitian Operators Dresden:
More informationChapter 3: Complex Numbers
Chapter 3: Complex Numbers Daniel Chan UNSW Semester 1 2018 Daniel Chan (UNSW) Chapter 3: Complex Numbers Semester 1 2018 1 / 48 Philosophical discussion about numbers Q In what sense is 1 a number? DISCUSS
More informationDUAL SPLIT QUATERNIONS AND SCREW MOTION IN MINKOWSKI 3-SPACE * L. KULA AND Y. YAYLI **
Iranian Journal of Science & echnology, ransaction A, Vol, No A Printed in he Islamic Republic of Iran, 6 Shiraz University DUAL SPLI UAERNIONS AND SCREW MOION IN MINKOWSKI -SPACE L KULA AND Y YAYLI Ankara
More informationCourse Name : Physics I Course # PHY 107
Course Name : Physics I Course # PHY 107 Lecture-2 : Representation of Vectors and the Product Rules Abu Mohammad Khan Department of Mathematics and Physics North South University http://abukhan.weebly.com
More informationQuaternion 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 informationQM and Angular Momentum
Chapter 5 QM and Angular Momentum 5. Angular Momentum Operators In your Introductory Quantum Mechanics (QM) course you learned about the basic properties of low spin systems. Here we want to review that
More informationThe scope of the midterm exam is up to and includes Section 2.1 in the textbook (homework sets 1-4). Below we highlight some of the important items.
AMS 10: Review for the Midterm Exam The scope of the midterm exam is up to and includes Section 2.1 in the textbook (homework sets 1-4). Below we highlight some of the important items. Complex numbers
More informationDynamics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Dynamics Semester 1, / 18
Dynamics Basilio Bona DAUIN Politecnico di Torino Semester 1, 2016-17 B. Bona (DAUIN) Dynamics Semester 1, 2016-17 1 / 18 Dynamics Dynamics studies the relations between the 3D space generalized forces
More informationBasic Concepts in Linear Algebra
Basic Concepts in Linear Algebra Grady B Wright Department of Mathematics Boise State University February 2, 2015 Grady B Wright Linear Algebra Basics February 2, 2015 1 / 39 Numerical Linear Algebra Linear
More information1 Matrices and matrix algebra
1 Matrices and matrix algebra 1.1 Examples of matrices A matrix is a rectangular array of numbers and/or variables. For instance 4 2 0 3 1 A = 5 1.2 0.7 x 3 π 3 4 6 27 is a matrix with 3 rows and 5 columns
More informationMatrix Algebra & Elementary Matrices
Matrix lgebra & Elementary Matrices To add two matrices, they must have identical dimensions. To multiply them the number of columns of the first must equal the number of rows of the second. The laws below
More informationReview of Basic Concepts in Linear Algebra
Review of Basic Concepts in Linear Algebra Grady B Wright Department of Mathematics Boise State University September 7, 2017 Math 565 Linear Algebra Review September 7, 2017 1 / 40 Numerical Linear Algebra
More informationEXERCISES. a b = a + b l aq b = ab - (a + b) + 2. a b = a + b + 1 n0i) = oii + ii + fi. A. Examples of Rings. C. Ring of 2 x 2 Matrices
/ rings definitions and elementary properties 171 EXERCISES A. Examples of Rings In each of the following, a set A with operations of addition and multiplication is given. Prove that A satisfies all the
More informationPrepared by: M. S. KumarSwamy, TGT(Maths) Page
Prepared by: M. S. KumarSwamy, TGT(Maths) Page - 50 - CHAPTER 3: MATRICES QUICK REVISION (Important Concepts & Formulae) MARKS WEIGHTAGE 03 marks Matrix A matrix is an ordered rectangular array of numbers
More informationThe groups SO(3) and SU(2) and their representations
CHAPTER VI The groups SO(3) and SU() and their representations Two continuous groups of transformations that play an important role in physics are the special orthogonal group of order 3, SO(3), and the
More informationKinematics. Basilio Bona. October DAUIN - Politecnico di Torino. Basilio Bona (DAUIN - Politecnico di Torino) Kinematics October / 15
Kinematics Basilio Bona DAUIN - Politecnico di Torino October 2013 Basilio Bona (DAUIN - Politecnico di Torino) Kinematics October 2013 1 / 15 Introduction The kinematic quantities used are: position r,
More informationA = 3 B = A 1 1 matrix is the same as a number or scalar, 3 = [3].
Appendix : A Very Brief Linear ALgebra Review Introduction Linear Algebra, also known as matrix theory, is an important element of all branches of mathematics Very often in this course we study the shapes
More informationSystems of Linear Equations and Matrices
Chapter 1 Systems of Linear Equations and Matrices System of linear algebraic equations and their solution constitute one of the major topics studied in the course known as linear algebra. In the first
More informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
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 informationSAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear Algebra
1.1. Introduction SAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear algebra is a specific branch of mathematics dealing with the study of vectors, vector spaces with functions that
More informationChapter XI Novanion rings
Chapter XI Novanion rings 11.1 Introduction. In this chapter we continue to provide general structures for theories in physics. J. F. Adams proved in 1960 that the only possible division algebras are at
More informationThe Quaternions. The Quaternions. John Huerta. Department of Mathematics UC Riverside. Cal State Stanislaus
John Huerta Department of Mathematics UC Riverside Cal State Stanislaus The Complex Numbers The complex numbers C form a plane. Their operations are very related to two dimensional geometry. In particular,
More informationA VERY BRIEF LINEAR ALGEBRA REVIEW for MAP 5485 Introduction to Mathematical Biophysics Fall 2010
A VERY BRIEF LINEAR ALGEBRA REVIEW for MAP 5485 Introduction to Mathematical Biophysics Fall 00 Introduction Linear Algebra, also known as matrix theory, is an important element of all branches of mathematics
More informationRobotics & Automation. Lecture 03. Representation of SO(3) John T. Wen. September 3, 2008
Robotics & Automation Lecture 03 Representation of SO(3) John T. Wen September 3, 2008 Last Time Transformation of vectors: v a = R ab v b Transformation of linear transforms: L a = R ab L b R ba R SO(3)
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 informationTABLE OF CONTENTS. Our aim is to give people math skillss in a very simple way Raymond J. Page 2 of 29
TABLE OF CONTENTS Topic.Page# 1. Numbers..04 2. Ratio, Profit & Loss 06 3. Angles......06 4. Interest..06 5. Algebra..07 6. Quadratic Equations.08 7. Logarithms...09 8. Series..10 9. Sphere.11 10. Coordinate
More informationSystems of Linear Equations and Matrices
Chapter 1 Systems of Linear Equations and Matrices System of linear algebraic equations and their solution constitute one of the major topics studied in the course known as linear algebra. In the first
More informationMATH Topics in Applied Mathematics Lecture 12: Evaluation of determinants. Cross product.
MATH 311-504 Topics in Applied Mathematics Lecture 12: Evaluation of determinants. Cross product. Determinant is a scalar assigned to each square matrix. Notation. The determinant of a matrix A = (a ij
More informationROBOTICS 01PEEQW. Basilio Bona DAUIN Politecnico di Torino
ROBOTICS 01PEEQW Basilio Bona DAUIN Politecnico di Torino Kinematic Functions Kinematic functions Kinematics deals with the study of four functions(called kinematic functions or KFs) that mathematically
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 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 information1.3 LECTURE 3. Vector Product
12 CHAPTER 1. VECTOR ALGEBRA Example. Let L be a line x x 1 a = y y 1 b = z z 1 c passing through a point P 1 and parallel to the vector N = a i +b j +c k. The equation of the plane passing through the
More informationLinear Algebra and Matrix Inversion
Jim Lambers MAT 46/56 Spring Semester 29- Lecture 2 Notes These notes correspond to Section 63 in the text Linear Algebra and Matrix Inversion Vector Spaces and Linear Transformations Matrices are much
More informationMathematical Fundamentals
Mathematical Fundamentals Ming-Hwa Wang, Ph.D. COEN 148/290 Computer Graphics COEN 166/366 Artificial Intelligence COEN 396 Interactive Multimedia and Game Programming Department of Computer Engineering
More informationCLASS 12 ALGEBRA OF MATRICES
CLASS 12 ALGEBRA OF MATRICES Deepak Sir 9811291604 SHRI SAI MASTERS TUITION CENTER CLASS 12 A matrix is an ordered rectangular array of numbers or functions. The numbers or functions are called the elements
More informationMatrix representation of a linear map
Matrix representation of a linear map As before, let e i = (0,..., 0, 1, 0,..., 0) T, with 1 in the i th place and 0 elsewhere, be standard basis vectors. Given linear map f : R n R m we get n column vectors
More informationA Tutorial on Euler Angles and Quaternions
A Tutorial on Euler Angles and Quaternions Moti Ben-Ari Department of Science Teaching Weimann Institute of Science http://www.weimann.ac.il/sci-tea/benari/ Version.0.1 c 01 17 b Moti Ben-Ari. This work
More informationRigid body dynamics. Basilio Bona. DAUIN - Politecnico di Torino. October 2013
Rigid body dynamics Basilio Bona DAUIN - Politecnico di Torino October 2013 Basilio Bona (DAUIN - Politecnico di Torino) Rigid body dynamics October 2013 1 / 16 Multiple point-mass bodies Each mass is
More informationPhysically Based Rendering ( ) Geometry and Transformations
Phsicall Based Rendering (6.657) Geometr and Transformations 3D Point Specifies a location Origin 3D Point Specifies a location Represented b three coordinates Infinitel small class Point3D { public: Coordinate
More informationMatrices over Hyperbolic Split Quaternions
Filomat 30:4 (2016), 913 920 DOI 102298/FIL1604913E Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://wwwpmfniacrs/filomat Matrices over Hyperbolic Split
More informationDr. Tess J. Moon, (office), (fax),
THE UNIVERSITY OF TEXAS AT AUSTIN DEPARTMENT OF MECHANICAL ENGINEERING ME 380Q ENGINEERING ANALYSIS: ANALYTICAL METHODS Dr Tess J Moon 471-0094 (office) 471-877 (fax) tmoon@mailutexasedu Mathematical Preliminaries
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationLecture Note 1: Background
ECE5463: Introduction to Robotics Lecture Note 1: Background Prof. Wei Zhang Department of Electrical and Computer Engineering Ohio State University Columbus, Ohio, USA Spring 2018 Lecture 1 (ECE5463 Sp18)
More informationNumerical Methods I Solving Square Linear Systems: GEM and LU factorization
Numerical Methods I Solving Square Linear Systems: GEM and LU factorization Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 18th,
More information