Two Worked out Examples of Rotations Using Quaternions

Two Worked out Examples of Rotations Using Quaternions

TWO WORKED OUT EXAMPLES OF ROTATIONS USING QUATERNIONS This note is an attachment to the article \Rotations and Quaternions" which in turn is a companion to the video of the talk by the same title. Example 1. Determine the image of the point (1; −1; 2) under the rotation by an angle of 60◦ about an axis in the yz-plane that is inclined at an angle of 60◦ to the positive y-axis. p ◦ ◦ 1 3 Solution: The unit vector u in the direction of the axis of rotation is cos 60 j + sin 60 k = 2 j + 2 k. The quaternion (or vector) corresponding to the point p = (1; −1; 2) is of course p = i − j + 2k. To find −1 θ θ the image of p under the rotation, we calculate qpq where q is the quaternion cos 2 + sin 2 u and θ the angle of rotation (60◦ in this case). The resulting quaternion|if we did the calculation right|would have no constant term and therefore we can interpret it as a vector. That vector gives us the answer. p p p p p We have q = 3 + 1 u = 3 + 1 j + 3 k = 1 (2 3 + j + 3k). Since q is by construction a unit quaternion, 2 2 2 4 p4 4 p −1 1 its inverse is its conjugate: q = 4 (2 3 − j − 3k). Now, computing qp in the routine way, we get 1 p p p p qp = ((1 − 2 3) + (2 + 3 3)i − 3j + (4 3 − 1)k) 4 and then another long but routine computation gives 1 p p p qpq−1 = ((10 + 4 3)i + (1 + 2 3)j + (14 − 3 3)k) 8 The point corresponding to the vector on the right hand side in the above equation is the image of (1; −1; 2) under the given rotation. Explicitly, that point is p p p 10 + 4 3 1 + 2 3 14 − 3 3 ( ; ; ) 8 8 8 Example 2. Consider the rotation by 60◦ about an axis in the xz-plane inclined at an angle of 60◦ to the positive x-axis. Determine the composition of this rotation with that in the earlier example: the one in the earlier example acts first. Solution: We need only compute q0q where q is as in the previous example and q0 = cos 30◦ + sin 30◦v where v is the unit vector in the direction of the axis of the rotation specified in this example. We have p p p p p ◦ ◦ 1 0 3 1 1 v = cos 60 i + sin 60 k = 2 (i + 3k). So q = 2 + 4 (i + 3k) = 4 (2 3 + i + 3k). Now a long but routine calculation gives: 1 p p 1 p p q0q = (2 3 + i + 3k) · (2 3 + j + 3k) 4 4 1 p p = (9 + 3i + 3j + 13k) 16 We now need to determine an angle θ and a unit vector v so that θ θ 1 p p cos + sin v = (9 + 3i + 3j + 13k) 2 2 16 Let us assume 0 ≤ θ ≤ 180◦, so that the value of cos θ (which is 9 in our case) determines θ . Then sin θ 2 2 p 16 2 2 θ 5 7 is non-negative and in the present case we get sin 2 = 16 . Now we determine v: 1 1 p p 1 p p v = · ( 3i + 3j + 13k) = p ( 3i + 3j + 13k) θ 16 sin 2 5 7 1 p p p = ( 21i + 21j + 13 7k) 35 The composite rotation is by the angle θ around the axis determined by the unit vector v. 1 The association p 7! qpq−1, where p is a point in 3-space interpreted as a vector and q a unit quaternion, is what enables us to attach rotations to unit quaternions. Since the association does not change if we replace q by −q, it follows that both q and −q are attached to the same rotation. Moreover, q and −q are the only ones that have the same association. Indeed, if qpq−1 = q0pq0−1 for all p, then qq0−1 commutes with all quaternions, and so has no imaginary component; being a unit quaternion as well, it can only be ±1; thus q0 = ±q. 2.

View Full Text

Details

  • File Type
    pdf
  • Upload Time
    -
  • Content Languages
    English
  • Upload User
    Anonymous/Not logged-in
  • File Pages
    2 Page
  • File Size
    -

Download

Channel Download Status
Express Download Enable

Copyright

We respect the copyrights and intellectual property rights of all users. All uploaded documents are either original works of the uploader or authorized works of the rightful owners.

  • Not to be reproduced or distributed without explicit permission.
  • Not used for commercial purposes outside of approved use cases.
  • Not used to infringe on the rights of the original creators.
  • If you believe any content infringes your copyright, please contact us immediately.

Support

For help with questions, suggestions, or problems, please contact us