Quaternions - Revision history

Yash at 10:45, 23 February 2018

2018-02-23T10:45:44Z

Yash at 22:03, 19 February 2018

2018-02-19T22:03:19Z

Yash at 19:05, 19 February 2018

2018-02-19T19:05:16Z

Yash at 16:24, 4 February 2018

2018-02-04T16:24:18Z

Yash at 16:22, 4 February 2018

2018-02-04T16:22:36Z

Yash at 16:22, 4 February 2018

2018-02-04T16:22:16Z

Yash at 20:47, 1 February 2018

2018-02-01T20:47:55Z

Yash: Created page with "''“There are two types of geniuses. Ordinary geniuses do great things, but they leave you room to believe that you could do the same if only you worked hard enough. Then the..."

2018-01-25T09:44:12Z

Created page with "''“There are two types of geniuses. Ordinary geniuses do great things, but they leave you room to believe that you could do the same if only you worked hard enough. Then the..."

New page

''“There are two types of geniuses. Ordinary geniuses do great things, but they leave you room to believe that you could do the same if only you worked hard enough. Then there are magicians, and you can have no idea how they do it.”''
----
It is surprising how simplified mathematical modelling of rotations has become with the advent of quaternions. One can’t stop admiring how each of the properties of quaternions has a subtle role to play in accomplishing this task of rotation. And these seem to have been defined with a magical foresight that they would all somehow fit into place later. 
Quaternions were constructed in attempt to extend the idea of rotations in a complex plane to 3D. Each quaternion is a set of four parameters. Why four? Well, any rigid body rotation can be done about a unique axis. So, 3 parameters to specify the axis and 1 for the angle rotated about it. But we have an endless choice of vectors along the axis to represent it. Here comes an applicational constraint from the fact that we want a rotation that doesn’t scale. So we restrict to unit quaternions (quaternions of unit magnitude). In retrospection, this has also simplified the mathematical design for rotation using quaternions. 
So we’ll start with vector notation for rotation and make an attempt to build a mathematical construct for unit quaternions.
[[File:Quaternion2.png|frame|center]]
In the picture above, x is being rotated to x’ about n in anticlockwise sense. This is equivalent to keeping its component parallel to n preserved and rotating the perpendicular component. Mathematically this is given as,
[[File:Equation13.png|frame|center]]
Now, a general quaternion is written as a + b i + c j + d k (or equivalently (a,'''v''')). ‘a’ is called the scalar part and the rest is the vector part (think of i,j,k to be similar to the unit orthogonal vectors of 3D space). 
To proceed, we’ll define few quaternion properties, 

* '''Addition'''
[[File:Equation14.png|frame|center]]

* '''Magnitude'''
[[File:Equation15.png|frame|center]]
(where q* is the conjugate quaternion of q. It’s scalar part is the same but the vector part is multiplied with (-1) )

* '''Multiplication'''
[[File:Equation16.JPG|frame|center]]
This is done defining ij=k ; jk=i ; ki=j ; ii=jj=kk= -1; So, ijk=-1. Notice, this is defined in accordance with the right- handed cross product of vectors. This was deliberate (as will be evident when we use quaternions to represent the vector rotation equation above). 
In short,
[[File:Equation17.png|frame|center]]
An useful observation here is that if above r1 =r2 =0 and v1 and v2 are perpendicular then 

(0, v1)(0, v2)= (0, v1 X v2). {we are doing vector cross product in v1 X v2} 
We’ll use this to update the vector rotation equation to quaternion notation. To extend vectors to quaternions, we’ll simply put the scalar part to zero and retain the vector notation in the vector part.
[[File:Equation18.JPG|frame|center]]

[[File:Equation19.JPG|frame|center]]
So, back to the rotation equation:
[[File:Equation20.JPG|frame|center]]
For unit quaternion multiplication with vector (in quaternion form),
[[File:Equation21.JPG|frame|center]]
The above can be directly verified from the rules of multiplication. 
So finally,
[[File:Equation22.JPG|frame|center]]
The above equation is used to perform rotations using unit quaternions. All along, the requirement of preserving vector magnitude never explicitly forced us to choose unit quaternions. However, this has been ensured because we started with a vector equation that took this into account and unit quaternion was an output of subsequent manipulations of this equation. 
Formulation of quaternions has made two major accomplishments. Firstly, it was able to give Euler’s rotation theorem (which says any rotation of a ‘rigid body’ is equivalent to a single rotation about a unique axis) a more solid mathematical application. Secondly, it was able to “mathematically” resolve the problem of gimbal lock ([https://www.youtube.com/watch?v=zc8b2Jo7mno here's a nice video on gimbal lock]) that euler angles face. It should be pointed out now that though quaternions provided a ‘mathematical’ alternative to euler angles, they didn’t fully dislodge them. For instance, there are still mechanical systems that perform successive rotations in a predefined manner to realise a final rotation. It is natural and convenient to describe these using euler rotation angles. Quaternions in this case can’t help circumvent the gimbal lock issue because the euler angles are a physical requirement. Quaternions are, hence, only a ‘mathematical’ antidote to gimbal lock. Nevertheless, they are a huge simplification of their existing mathematical counterparts. 
If you are done reading this page, you can go back to [[Attitude Determination and Control Subsystem]]

@@ Line 4: / Line 4: @@
 Quaternions were constructed in attempt to extend the idea of rotations in a complex plane to 3D. Each quaternion is a set of four parameters. Why four? Well, any rigid body rotation can be done about a unique axis. So, 3 parameters to specify the axis and 1 for the angle rotated about it. But we have an endless choice of vectors along the axis to represent it. Here comes an applicational constraint from the fact that we want a rotation that doesn’t scale. So we restrict to unit quaternions (quaternions of unit magnitude). In retrospection, this has also simplified the mathematical design for rotation using quaternions. <br \>
 So we’ll start with vector notation for rotation and make an attempt to build a mathematical construct for unit quaternions.
-[[File:Quaternions.png|thumb|center|500px|Image inspired from [http://www.mlahanas.de/Math/orientation.htm here]]]
+[[File:Quaternions.png|thumb|center|500px|Decomposing a given point x into components parallel and perpendicular to the rotation axis. Image inspired from [http://www.mlahanas.de/Math/orientation.htm here]]]
 In the picture above, x is being rotated to x’ about n in anticlockwise sense. This is equivalent to keeping its component parallel to n preserved and rotating the perpendicular component. Mathematically this is given as,
 [[File:Equation13.png|frame|center]]

@@ Line 4: / Line 4: @@
 Quaternions were constructed in attempt to extend the idea of rotations in a complex plane to 3D. Each quaternion is a set of four parameters. Why four? Well, any rigid body rotation can be done about a unique axis. So, 3 parameters to specify the axis and 1 for the angle rotated about it. But we have an endless choice of vectors along the axis to represent it. Here comes an applicational constraint from the fact that we want a rotation that doesn’t scale. So we restrict to unit quaternions (quaternions of unit magnitude). In retrospection, this has also simplified the mathematical design for rotation using quaternions. <br \>
 So we’ll start with vector notation for rotation and make an attempt to build a mathematical construct for unit quaternions.
-[[File:Quaternions.png|frame|center| Image inspired from [http://www.mlahanas.de/Math/orientation.htm here]]]
+[[File:Quaternions.png|thumb|center|500px|Image inspired from [http://www.mlahanas.de/Math/orientation.htm here]]]
 In the picture above, x is being rotated to x’ about n in anticlockwise sense. This is equivalent to keeping its component parallel to n preserved and rotating the perpendicular component. Mathematically this is given as,
 [[File:Equation13.png|frame|center]]

@@ Line 4: / Line 4: @@
 Quaternions were constructed in attempt to extend the idea of rotations in a complex plane to 3D. Each quaternion is a set of four parameters. Why four? Well, any rigid body rotation can be done about a unique axis. So, 3 parameters to specify the axis and 1 for the angle rotated about it. But we have an endless choice of vectors along the axis to represent it. Here comes an applicational constraint from the fact that we want a rotation that doesn’t scale. So we restrict to unit quaternions (quaternions of unit magnitude). In retrospection, this has also simplified the mathematical design for rotation using quaternions. <br \>
 So we’ll start with vector notation for rotation and make an attempt to build a mathematical construct for unit quaternions.
-[[File:Quaternion2.png|frame|center| Image reproduced from [http://www.mlahanas.de/Math/orientation.htm here]]]
+[[File:Quaternions.png|frame|center| Image inspired from [http://www.mlahanas.de/Math/orientation.htm here]]]
 In the picture above, x is being rotated to x’ about n in anticlockwise sense. This is equivalent to keeping its component parallel to n preserved and rotating the perpendicular component. Mathematically this is given as,
 [[File:Equation13.png|frame|center]]

← Older revision		Revision as of 16:22, 4 February 2018
Line 40:		Line 40:
	----		----
	If you are done reading this page, you can go back to [[Attitude Determination and Control Subsystem]]		If you are done reading this page, you can go back to [[Attitude Determination and Control Subsystem]]
		+
		+	== References ==

@@ Line 8: / Line 8: @@
 [[File:Equation13.png|frame|center]]
 Now, a general quaternion is written as a + b i + c j + d k (or equivalently (a,'''v''')). ‘a’ is called the scalar part and the rest is the vector part (think of i,j,k to be similar to the unit orthogonal vectors of 3D space). <br \>
-To proceed, we’ll define few quaternion properties, <br \>
+To proceed, we’ll define few quaternion properties, <ref>https://en.wikipedia.org/wiki/Quaternion </ref><br \>
 * '''Addition'''

@@ Line 38: / Line 38: @@
 The above equation is used to perform rotations using unit quaternions. All along, the requirement of preserving vector magnitude never explicitly forced us to choose unit quaternions. However, this has been ensured because we started with a vector equation that took this into account and unit quaternion was an output of subsequent manipulations of this equation. <br \>
 Formulation of quaternions has made two major accomplishments. Firstly, it was able to give Euler’s rotation theorem (which says any rotation of a ‘rigid body’ is equivalent to a single rotation about a unique axis) a more solid mathematical application. Secondly, it was able to “mathematically” resolve the problem of gimbal lock ([https://www.youtube.com/watch?v=zc8b2Jo7mno here's a nice video on gimbal lock]) that euler angles face. It should be pointed out now that though quaternions provided a ‘mathematical’ alternative to euler angles, they didn’t fully dislodge them. For instance, there are still mechanical systems that perform successive rotations in a predefined manner to realise a final rotation. It is natural and convenient to describe these using euler rotation angles. Quaternions in this case can’t help circumvent the gimbal lock issue because the euler angles are a physical requirement. Quaternions are, hence, only a ‘mathematical’ antidote to gimbal lock. Nevertheless, they are a huge simplification of their existing mathematical counterparts. <br \>
 If you are done reading this page, you can go back to [[Attitude Determination and Control Subsystem]]