|
|
| Line 1: |
Line 1: |
| {{About|the Euler angles used in [[mathematics]]|the use of the term in [[physics]] and [[aerospace engineering]]|Rigid body dynamics}}
| | Read this article and learn these tips to prevent water damage at home. During those situations remember that the first thing to do is to approach a professional water or fire restoration services. Morgan Vawter is the managing director of which is the leader in market of Water Damage Control, catering to the needs of Clients all across the globe. These professional will find roots and structural damage in you sewer line and will give you with an accurate report on your sewer form. A key lock can be a simple lock affixed to a door (deadlock, handle locks, etc. <br><br>When you have almost any questions regarding in which along with the best way to work with [http://tupalo.com/en/plano-texas/dallas-roofing-consultants Plano roofer], you possibly can email us in our own internet site. It's always important to practice everything in moderation and most importantly, to take care of your body and your health. Clean and well organized kitchen is always found in a successful man's life. " Extract standing water in the building and dry wet surfaces and items as quickly as possible to stop mold growth. The process of cleaning up homes damaged by water and restoring them to their pre-loss condition is called water damage restoration and when facing such circumstances it is most often in your best interest to hire a professional water damage restoration company. Another technique is cover the item in toothpaste that has baking soda and let it stand for 15 minutes, soften your toothbrush with water and brush the jewelry item gently, continue to dip your brush into water until you get a good lather. <br><br>A bronze substrate, for example, may need acrylic-based coatings that do not create a galvanic reaction with the bronze substrate. Fabric panels do not meet the impact certification test and will not. If they can't, it's important to remember that towns come and go, but extinct species are forever. We have to be prepared when misfortunes like this occurs. Yes, it is time consuming but the adjuster will be able to see the reasoning behind the line item. <br><br>Some say wash it everyday for healthy strands - others say wash every other day. Going online for putting your cherished photographs to canvas is a great idea. Side note: Be sure to call a plumber to stop the source of the water from flooding your home. They've their modern facilities and devices also for a far more prosperous process. ' Air out drench locations or areas as soon as possible. <br><br>Wipe out drawers and cabinets and blot them as well. If you must keep one or more servers up during a long-term power outage, consider a generator. Make sure that your work authorization form includes a section in it that would allow the insurance company to be able to pay you directly. These water restoration teams are very willing to help you out concerning this matter. Due to its best selling price and excellent worldwide customer support i - Pod recovery software has become popular in international market. |
| | |
| The '''Euler angles''' are three angles introduced by [[Leonhard Euler]] to describe the [[Orientation (geometry)|orientation]] of a [[rigid body]].<ref>Novi Commentarii academiae scientiarum Petropolitanae 20, 1776, pp. 189–207 (E478) [http://math.dartmouth.edu/~euler/docs/originals/E478.pdf pdf]</ref> To describe such an orientation in [[dimension|3-dimensional]] [[Euclidean space]] three parameters are required. They can be given in several ways, Euler angles being one of them; see [[charts on SO(3)]] for others. Euler angles are also used to represent the orientation of a [[frame of reference]] (typically, a [[coordinate system]] or [[Basis (linear algebra)|basis]]) relative to another. They are typically denoted as α, β, γ, or <math>\varphi, \theta, \psi</math>.
| |
| | |
| Euler angles also represent a sequence of three ''elemental'' [[rotation]]s, i.e. rotations about the axes of a [[coordinate system]]. For instance, a first rotation about ''z'' by an angle α, a second rotation about ''x'' by an angle β, and a last rotation again about ''z'', by an angle γ. These rotations start from a known standard orientation. In physics, this standard initial orientation is typically represented by a motionless (''fixed'', ''global'', or ''world'') [[coordinate system]]. In linear algebra, by a [[standard basis]].
| |
| | |
| Any orientation can be achieved by composing three elemental rotations. The elemental rotations can either occur about the axes of the fixed coordinate system ([[#Extrinsic rotations|extrinsic rotations]]) or about the axes of a rotating coordinate system, which is initially aligned with the fixed one, and modifies its orientation after each elemental rotation ([[#Intrinsic rotations|intrinsic rotations]]). The rotating coordinate system may be imagined to be rigidly attached to a rigid body. In this case, it is sometimes called a ''local'' coordinate system. Without considering the possibility of using two different conventions for the definition of the rotation axes (intrinsic or extrinsic), there exist twelve possible sequences of rotation axes, divided in two groups:
| |
| * '''Euler angles''' (''z-x-z, x-y-x, y-z-y, z-y-z, x-z-x, y-x-y'')
| |
| * '''Tait–Bryan angles''' (''x-y-z, y-z-x, z-x-y, x-z-y, z-y-x, y-x-z'').
| |
| Tait–Bryan angles are also called '''Cardan angles''', '''nautical angles''', '''heading, elevation, and bank''', or '''yaw, pitch, and roll'''. Sometimes, both kinds of sequences are called "Euler angles". In that case, the sequences of the first group are called ''proper'' or ''classic'' Euler angles.
| |
| | |
| ==Proper Euler angles==
| |
| [[Image:Eulerangles.svg|right|thumb|300px|Proper Euler angles representing rotations about ''z''-, ''N''-, and ''Z''-axis. The ''xyz'' (original) system is shown in blue, the ''XYZ'' (rotated) system is shown in red. The [[line of nodes]] (''N'') is shown in green.]]
| |
| | |
| ===Classic definition===
| |
| Euler angles are a means of representing the spatial orientation of any [[Frame of reference|reference frame]] ([[coordinate system]] or [[Basis (linear algebra)|basis]]) as a composition of three elemental rotations starting from a known standard orientation, represented by another frame (sometimes referred to as the ''original'' or ''fixed'' reference frame, or [[standard basis]]). The reference orientation can be imagined to be an initial orientation from which the frame virtually rotates to reach its actual orientation. In the following, the axes of the original frame are denoted as ''x'',''y'',''z'' and the axes of the rotated frame are denoted as ''X'',''Y'',''Z''. In [[geometry]] and [[physics]], the rotated coordinate system is often imagined to be rigidly attached to a [[rigid body]]. In this case, it is called a "local" coordinate system, and it is meant to represent both the position and the orientation of the body.
| |
| | |
| The geometrical definition (referred sometimes as static) of the Euler angles is based on the axes of the above-mentioned (original and rotated) reference frames and an additional axis called the [[line of nodes]]. The line of nodes (''N'') is defined as the intersection of the ''xy'' and the ''XY'' coordinate planes. In other words, it is a line passing through the origin of both frames, and perpendicular to the ''zZ'' plane, on which both ''z'' and ''Z'' lie. The three Euler angles are defined as follows:
| |
| *α (or <math>\varphi</math>) is the angle between the ''x''-axis and the ''N''-axis.
| |
| *β (or <math>\theta</math>) is the angle between the ''z''-axis and the ''Z''-axis.
| |
| *γ (or <math>\psi</math>) is the angle between the ''N''-axis and the ''X''-axis.
| |
| | |
| This definition implies that:
| |
| *α represents a rotation around the ''z''-axis,
| |
| *β represents a rotation around the ''N''-axis,
| |
| *γ represents a rotation around the ''Z''-axis.
| |
| | |
| If β is zero, there is no rotation about ''N''. As a consequence, ''Z'' coincides with ''z'', α and γ represent rotations about the same axis (''z''), and the final orientation can be obtained with a single rotation about ''z'', by an angle equal to α+γ.
| |
| | |
| ===Alternative definition===
| |
| The rotated frame ''XYZ'' may be imagined to be initially aligned with ''xyz'', before undergoing the three elemental rotations represented by Euler angles. Its successive orientations may be denoted as follows:
| |
| * ''x-y-z'', or ''x''<sub>0</sub>-''y''<sub>0</sub>-''z''<sub>0</sub> (initial)
| |
| * ''x’-y’-z’'', or ''x''<sub>1</sub>-''y''<sub>1</sub>-''z''<sub>1</sub> (after first rotation)
| |
| * ''x″-y″-z″'', or ''x''<sub>2</sub>-''y''<sub>2</sub>-''z''<sub>2</sub> (after second rotation)
| |
| * ''X-Y-Z'', or ''x''<sub>3</sub>-''y''<sub>3</sub>-''z''<sub>3</sub> (final)
| |
| | |
| For the above-listed sequence of rotations, the [[line of nodes]] ''N'' can be simply defined as the orientation of ''X'' after the first elemental rotation. Hence, ''N'' can be simply denoted ''x’''. Moreover, since the third elemental rotation occurs about ''Z'', it does not change the orientation of ''Z''. Hence ''Z'' coincides with ''z″''. This allows us to simplify the definition of the Euler angles as follows:
| |
| *α (or <math>\varphi</math>) represents a rotation around the ''z''-axis,
| |
| *β (or <math>\theta</math>) represents a rotation around the ''x’''-axis,
| |
| *γ (or <math>\psi</math>) represents a rotation around the ''z″''-axis.
| |
| | |
| ===Conventions===
| |
| Different authors may use different sets of rotation axes to define Euler angles, or different names for the same angles. Therefore any discussion employing Euler angles should always be preceded by their definition.<ref>[http://mathworld.wolfram.com/EulerAngles.html Mathworld] does a good job describing this issue</ref> Unless otherwise stated, this article will use the convention described above.
| |
| | |
| The three elemental rotations may occur either about the axes ''xyz'' of the original coordinate system, which is assumed to remain motionless ([[#Extrinsic rotations|extrinsic rotations]]), or about the axes of the rotating coordinate system ''XYZ'', which changes its orientation after each elemental rotation ([[#Intrinsic rotations|intrinsic rotations]]). The definition above uses intrinsic rotations.
| |
| | |
| There are six possibilities of choosing the rotation axes for proper Euler angles. In all of them, the first and third rotation axis are the same. The six possible sequences are:
| |
| # ''z-x’-z″'' (intrinsic rotations) or ''z-x-z'' (extrinsic rotations)
| |
| # ''x-y’-x″'' (intrinsic rotations) or ''x-y-x'' (extrinsic rotations)
| |
| # ''y-z’-y″'' (intrinsic rotations) or ''y-z-y'' (extrinsic rotations)
| |
| # ''z-y’-z″'' (intrinsic rotations) or ''z-y-z'' (extrinsic rotations)
| |
| # ''x-z’-x″'' (intrinsic rotations) or ''x-z-x'' (extrinsic rotations)
| |
| # ''y-x’-y″'' (intrinsic rotations) or ''y-x-y'' (extrinsic rotations)
| |
| | |
| Euler angles between two reference frames are defined only if both frames have the same [[Orientation (mathematics)|handedness]].
| |
| | |
| ===Signs and ranges===
| |
| Angles are commonly defined according to the [[right hand rule]]. Namely, they have positive values when they represent a rotation that appears counter-clockwise when observed from a point lying on the positive part of the rotation axis, and negative values when the rotation appears clockwise. The opposite convention (left hand rule) is less frequently adopted.
| |
| | |
| About the ranges:
| |
| *for α and γ, the range is defined [[Modular arithmetic|modulo]] 2π [[radian]]s. A valid range could be [−π, π].
| |
| *for β, the range covers π radians (but can't be said to be modulo π). For example could be [0, π] or [−π/2, π/2].
| |
| | |
| The angles α, β and γ are uniquely determined except for the singular case that the ''xy'' and the ''XY'' planes are identical, the ''z'' axis and the ''Z'' axis having the same or opposite directions. Indeed, if the ''z''-axis and the Z-axis are the same, β = 0 and only (α + γ) is uniquely defined (not the individual values), and, similarly, if the ''z''-axis and the ''Z''-axis are opposite, β = π and only (α − γ) is uniquely defined (not the individual values). These ambiguities are known as [[gimbal lock]] in applications.
| |
| | |
| ===Geometric derivation===
| |
| [[Image:EulerProjections.svg|right|thumb|200px|Projections of Z vector.]]
| |
| [[Image:EulerProjections2.svg|right|thumb|200px|Projections of Y vector.]]
| |
| The fastest way to get the Euler Angles of a given frame is to write the three given vectors as columns of a matrix and compare it with the expression of the theoretical matrix (see later table of matrices). Hence the three Euler Angles can be calculated. Nevertheless, the same result can be reached avoiding matrix algebra, which is more geometrical. Assuming a frame with unitary vectors (''X'', ''Y'', ''Z'') as in the main diagram, it can be seen that:
| |
| | |
| :<math>\cos (\beta) = Z_3.</math>
| |
| | |
| And, since
| |
| :<math>\sin^2 x = 1 - \cos^2 x,</math>
| |
| we have
| |
| :<math>\sin (\beta) = \sqrt {1 - Z_3^2}.</math>
| |
| | |
| As <math>Z_2</math> is the double projection of a unitary vector,
| |
| :<math>\cos(\alpha) \cdot \sin(\beta) = -Z_2,</math>
| |
| :<math>\cos(\alpha) = -Z_2 / \sqrt{1 - Z_3^2}.</math>
| |
| | |
| There is a similar construction for <math>Y_3</math>, projecting it first over the plane defined by the axis ''z'' and the line of nodes. As the angle between the planes is <math>\pi/2 - \beta</math> and <math>\cos(\pi/2 - \beta) = \sin(\beta)</math>, this leads to:
| |
| | |
| :<math>\sin(\beta) \cdot \cos(\gamma) = Y_3,</math>
| |
| :<math>\cos(\gamma) = Y_3 / \sqrt{1 - Z_3^2},</math>
| |
| | |
| and finally, using the [[inverse trigonometric functions|inverse cosine]] function,
| |
| | |
| :<math>\alpha = \arccos(-Z_2 / \sqrt{1 - Z_3^2}),</math>
| |
| :<math>\beta = \arccos(Z_3),</math>
| |
| :<math>\gamma = \arccos(Y_3 / \sqrt{1 - Z_3^2}).</math>
| |
| | |
| It is interesting to note that the inverse cosine function yields two possible values for the argument. In this geometrical description only one of the solutions is valid. When Euler Angles are defined as a sequence of rotations, all the solutions can be valid, but there will be only one inside the angle ranges. This is because the sequence of rotations to reach the target frame is not unique if the ranges are not previously defined.<ref>[http://www.soi.city.ac.uk/~sbbh653/publications/euler.pdf Gregory G. Slabaugh, Computing Euler angles from a rotation matrix]</ref>
| |
| | |
| For computational purposes, it may be useful to represent the angles using [[atan2]](y,x):
| |
| | |
| :<math>\alpha = \operatorname{atan2}(Z_1 , -Z_2),</math>
| |
| :<math>\gamma = \operatorname{atan2}(X_3 , Y_3).</math>
| |
| | |
| ==Tait–Bryan angles==
| |
| [[Image:taitbrianzyx.svg|thumb|left|200px|Tait–Bryan angles. ''z-y’-x″'' sequence (intrinsic rotations; ''N'' coincides with ''y’'')]]
| |
| | |
| The second type of formalism is called '''Tait–Bryan angles''', after [[Peter Guthrie Tait]] and [[George H. Bryan]].
| |
| | |
| The definitions and notations used for Tait-Bryan angles are similar to those described above for proper Euler angles ([[#Classic definition|Classic definition]], [[#Alternative definition|Alternative definition]]). The only difference is that Tait–Bryan angles represent rotations about three distinct axes (e.g. ''x-y-z'', or ''x-y’-z″''), while proper Euler angles use the same axis for both the first and third elemental rotations (e.g., ''z-x-z'', or ''z-x’-z″'').
| |
| | |
| This implies a different definition for the [[line of nodes]]. In the first case it was defined as the intersection between two homologous Cartesian planes (parallel when Euler angles are zero; e.g. ''xy'' and ''XY''). In the second one, it is defined as the intersection of two non-homologous planes (perpendicular when Euler angles are zero; e.g. ''xy'' and ''YZ'').
| |
| | |
| [[Image:taitbrianangles.svg|thumb|right|325px|Tait–Bryan angles. ''z-x’-y″'' sequence (intrinsic rotations; ''N'' coincides with ''x’'')]]
| |
| | |
| ===Conventions===
| |
| The three elemental rotations may occur either about the axes of the original coordinate system, which remains motionless ([[#Extrinsic rotations|extrinsic rotations]]), or about the axes of the rotating coordinate system, which changes its orientation after each elemental rotation ([[#Intrinsic rotations|intrinsic rotations]]).
| |
| | |
| There are six possibilities of choosing the rotation axes for Tait–Bryan angles. The six possible sequences are:
| |
| # ''x-y’-z″'' (intrinsic rotations) or ''x-y-z'' (extrinsic rotations)
| |
| # ''y-z’-x″'' (intrinsic rotations) or ''y-z-x'' (extrinsic rotations)
| |
| # ''z-x’-y″'' (intrinsic rotations) or ''z-x-y'' (extrinsic rotations)
| |
| # ''x-z’-y″'' (intrinsic rotations) or ''x-z-y'' (extrinsic rotations)
| |
| # ''z-y’-x″'' (intrinsic rotations) or ''z-y-x'' (extrinsic rotations)
| |
| # ''y-x’-z″'' (intrinsic rotations) or ''y-x-z'' (extrinsic rotations)
| |
| | |
| ===Alternative names===
| |
| [[Image:Yaw Axis Corrected.svg|thumb|right|250px|The [[Aircraft principal axes|principal axes]] of an aircraft]]
| |
| Tait-Bryan angles are also known as '''nautical angles''', because they can be used to describe the orientation of a ship or aircraft, or '''Cardan angles''', after the Italian mathematician and physicist [[Gerolamo Cardano]] ({{lang-fr|Jérôme Cardan}}; 24 September 1501 – 21 September 1576) who first described in detail the [[Cardan suspension]] and the [[Cardan joint]].
| |
| | |
| They are also called '''heading, elevation and bank''', or '''yaw, pitch and roll'''. Notice that the second set of terms is also used for the three [[aircraft principal axes]].
| |
| | |
| == Relationship with physical motions ==
| |
| :See also [[Givens rotation]]s
| |
| | |
| ===Intrinsic rotations===
| |
| | |
| {{multiple image
| |
| | width1 = 170
| |
| | footer = Any target orientation can be reached, starting from a known reference orientation, using a specific sequence of intrinsic rotations, whose magnitudes are the Euler Angles of the target orientation. This example uses the z-x'-z" sequence.
| |
| | image1 = euler2a.gif
| |
| | alt1 =
| |
| | width2 = 150
| |
| | image2 = intermediateframes.svg
| |
| | alt2 =
| |
| | caption2 =
| |
| }}
| |
| | |
| Intrinsic rotations are elemental rotations that occur about the axes of the rotating coordinate system ''XYZ'', which changes its orientation after each elemental rotation. The ''XYZ'' system rotates, while ''xyz'' is fixed. Starting with ''XYZ'' overlapping ''xyz'', a composition of three intrinsic rotations can be used to reach any target orientation for ''XYZ''. The Euler or Tait Bryan angles (α, β, γ) are the amplitudes of these elemental rotations. For instance, the target orientation can be reached as follows:
| |
| | |
| *The ''XYZ''-system rotates by α about the ''Z''-axis (which coincides with the ''z''-axis). The ''X''-axis now lies on the line of nodes.
| |
| *The ''XYZ''-system rotates about the now rotated ''X''-axis by β. The ''Z''-axis is now in its final orientation, and the ''X''-axis remains on the line of nodes.
| |
| *The ''XYZ''-system rotates a third time about the new ''Z''-axis by γ.
| |
| | |
| The above-mentioned [[#Alternative definition|notation]] allows us to summarize this as follows: the three elemental rotations of the XYZ-system occur about ''z'', ''x''’ and ''z''″. Indeed, this sequence is often denoted ''z-x’-z″''. Sets of rotation axes associated with both proper Euler angles and Tait Bryan angles are commonly named using this notation (see above for details). Sometimes, the same sequence is simply called ''z-x-z'', ''Z-X-Z'', or ''3-1-3'', but this notation may be ambiguous as it may be identical to that used for extrinsic rotations. In this case, it becomes necessary to separately specify whether the rotations are intrinsic or extrinsic.
| |
| | |
| [[Rotation matrix|Rotation matrices]] can be used to represent a sequence of intrinsic rotations. For instance,
| |
| :<math>R = X(\alpha) Y(\beta) Z(\gamma)</math>
| |
| represents a composition of intrinsic rotations about axes ''x-y’-z″'', if used to pre-multiply [[column vector]]s, while
| |
| :<math>R = Z(\gamma) Y(\beta) X(\alpha) </math>
| |
| represents exactly the same composition when used to post-multiply [[row vector]]s. See [[Rotation matrix#Ambiguities|Ambiguities in the definition of rotation matrices]] for more details.
| |
| | |
| ===Extrinsic rotations===
| |
| | |
| Extrinsic rotations are elemental rotations that occur about the axes of the fixed coordinate system ''xyz''. The ''XYZ'' system rotates, while ''xyz'' is fixed. Starting with ''XYZ'' overlapping ''xyz'', a composition of three extrinsic rotations can be used to reach any target orientation for ''XYZ''. The Euler or Tait Bryan angles (α, β, γ) are the amplitudes of these elemental rotations. For instance, the target orientation can be reached as follows:
| |
| | |
| *The ''XYZ''-system rotates about the ''z''-axis by α. The ''X''-axis is now at angle α with respect to the ''x''-axis.
| |
| *The ''XYZ''-system rotates again about the ''x''-axis by β. The ''Z''-axis is now at angle β with respect to the ''z''-axis.
| |
| *The ''XYZ''-system rotates a third time about the ''z''-axis by γ.
| |
| | |
| In sum, the three elemental rotations occur about ''z'', ''x'' and ''z''. Indeed, this sequence is often denoted ''z-x-z'' (or 3-1-3). Sets of rotation axes associated with both proper Euler angles and Tait Bryan angles are commonly named using this notation (see above for details).
| |
| | |
| [[Rotation matrix|Rotation matrices]] can be used to represent a sequence of extrinsic rotations. For instance,
| |
| :<math>R = Z(\gamma) Y(\beta) X(\alpha) </math>
| |
| represents a composition of extrinsic rotations about axes ''x-y-z'', if used to pre-multiply [[column vector]]s, while
| |
| :<math>R = X(\alpha) Y(\beta) Z(\gamma)</math>
| |
| represents exactly the same composition when used to post-multiply [[row vector]]s. See [[Rotation matrix#Ambiguities|Ambiguities in the definition of rotation matrices]] for more details.
| |
| | |
| ===Conversion between intrinsic and extrinsic rotations===
| |
| [[Image:EulerG.png|thumb|300px|A rotation represented by Euler angles (α, β, γ) = (−60°, 30°, 45°), using ''z-x’-z″'' intrinsic rotations]]
| |
| [[Image:EulerX.png|thumb|300px|The same rotation represented by (γ, β, α) = (45°, 30°, −60°), using ''z-x-z'' extrinsic rotations]]
| |
| | |
| Any extrinsic rotation is equivalent to an intrinsic rotation by the same angles but with inverted order of elemental rotations, and vice-versa. For instance, the intrinsic rotations ''x-y’-z″'' by angles α, β, γ are equivalent to the extrinsic rotations ''z-y-x'' by angles γ, β, α. Both are represented by a matrix
| |
| | |
| :<math>R = X(\alpha) Y(\beta) Z(\gamma)</math>
| |
| if ''R'' is used to pre-multiply [[column vector]]s, and by a matrix
| |
| :<math>R = Z(\gamma) Y(\beta) X(\alpha) </math>
| |
| if ''R'' is used to post-multiply [[row vector]]s. See [[Rotation matrix#Ambiguities|Ambiguities in the definition of rotation matrices]] for more details.
| |
| | |
| ===Euler rotations===
| |
| [[Image:Praezession.svg|thumb|170px|left|Euler rotations of the Earth. Intrinsic (green), Precession (blue) and Nutation (red)]]
| |
| Euler rotations are defined as the movement obtained by changing one of the Euler angles while leaving the other two constant. Euler rotations are never expressed in terms of the external frame, or in terms of the co-moving rotated body frame, but in a mixture. They constitute a '''mixed axes of rotation''' system, where the first angle moves the line of nodes around the external axis z, the second rotates around the line of nodes and the third one is an intrinsic rotation around an axis fixed in the body that moves.
| |
| | |
| These rotations are called [[precession]], [[nutation]], and [[rotation|intrinsic rotation]] (spin). As an example, consider a [[top]]. The top spins around its own axis of symmetry; this corresponds to its intrinsic rotation. It also rotates around its pivotal axis, with its center of mass orbiting the pivotal axis; this rotation is a precession. Finally, the top can wobble up and down; the inclination angle is the nutation angle. While all three are rotations when applied over individual frames, only precession is valid as a rotation operator, and only precession can be expressed in general as a matrix in the basis of the space.
| |
| | |
| ====Gimbal analogy====
| |
| {{multiple image
| |
| | image1 = gimbaleuler.svg
| |
| | width1 = 180
| |
| | image2 = gimbaleuler2.svg
| |
| | width2 = 160
| |
| | footer = '''Left:''' A three axes z-x-z-[[gimbal]] where the external frame and external axis 'x' are not shown and axes 'Y' are perpendicular to each gimbal ring. '''Right:''' A simple diagram showing the Euler angles and where the axes 'Y' of intermediate frames are located.
| |
| }}
| |
| | |
| If we suppose a set of frames, able to move each with respect to the former according to just one angle, like a gimbal, there will exist an external fixed frame, one final frame and two frames in the middle, which are called "intermediate frames". The two in the middle work as two gimbal rings that allow the last frame to reach any orientation in space. | |
| | |
| ====Intermediate frames====
| |
| The gimbal rings indicate some intermediate frames. They can be defined statically too. Taking some vectors '''i''', '''j''' and '''k''' over the axes x, y and z, and vectors '''I''', '''J''', '''K''' over X, Y and Z, and a vector '''N''' over the line of nodes, some intermediate frames can be defined using the vector [[cross product]], as following:
| |
| | |
| *origin: ['''i''','''j''','''k'''] (where '''k''' = '''i''' × '''j''')
| |
| *first: ['''N''','''k''' × '''N''','''k''']
| |
| *second: ['''N''','''K''' × '''N''','''K''']
| |
| *final: ['''I''','''J''','''K''']
| |
| | |
| These intermediate frames are equivalent to those of the gimbal. They are such that they differ from the previous one in just a single elemental rotation. This proves that:
| |
| | |
| *Any target frame can be reached from the reference frame just composing three rotations.
| |
| *The values of these three rotations are exactly the Euler angles of the target frame.
| |
| | |
| ==Relationship to other representations==
| |
| {{Main|Rotation formalisms in three dimensions#Conversion formulae between formalisms}}
| |
| | |
| Euler angles are one way to represent orientations. There are others, and it is possible to change to and from other conventions.
| |
| | |
| ===Rotation matrix===
| |
| | |
| Any orientation can be achieved by composing three elemental rotations, starting from a known standard orientation. Equivalently, any [[rotation matrix]] ''R'' can be [[Matrix decomposition|decomposed]] as a product of three elemental rotation matrices. For instance:
| |
| :<math>R = X(\alpha) Y(\beta) Z(\gamma)</math>
| |
| is a rotation matrix that may be used to represent a composition of [[Euler angles#Intrinsic rotations|intrinsic rotations]] about axes ''x-y’-z″''. However, both the definition of the elemental rotation matrices ''X'', ''Y'', ''Z'', and their multiplication order depend on the choices taken by the user about the definition of both rotation matrices and Euler angles (see, for instance, [[Rotation matrix#Ambiguities|Ambiguities in the definition of rotation matrices]]). Unfortunately, different sets of conventions are adopted by users in different contexts. The following table was built according to this set of conventions:
| |
| # Each matrix is meant to operate by pre-multiplying [[column vector]]s (see [[Rotation matrix#Ambiguities|Ambiguities in the definition of rotation matrices]])
| |
| # Each matrix is meant to represent an [[Active and passive transformation|active rotation]] (the composing and composed matrices are supposed to act on the coordinates of vectors defined in the initial fixed reference frame and give as a result the coordinates of a rotated vector defined in the same reference frame).
| |
| # Each matrix is meant to represent the composition of [[intrinsic rotations]] (around the axes of the rotating reference frame).
| |
| # [[Right hand rule|Right handed]] reference frames are adopted, and the [[right hand rule]] is used to determine the sign of the angles α, β, γ.
| |
| For the sake of simplicity, the following table uses the following nomenclature:
| |
| # 1, 2, 3 represent the angles α, β, γ.
| |
| # ''X, Y, Z'' are the matrices representing the elemental rotations about the axes ''x, y, z'' of the fixed frame (e.g., ''X''<sub>1</sub> represents a rotation about ''x'' by an angle α).
| |
| # ''s'' and ''c'' represent sine and cosine (e.g., ''s''<sub>1</sub> represents the sine of α).
| |
| # Each matrix is denoted by the formula used to calculate it. If <math>R = Z_1 X_2 Z_3</math>, we name it <math>Z_1 X_2 Z_3</math>.
| |
| | |
| :{| class="wikitable" style="background-color:white;font-weight:bold"
| |
| |-
| |
| ! Proper Euler angles !! Tait-Bryan angles
| |
| |-
| |
| |<math>X_1 Z_2 X_3 = \begin{bmatrix}
| |
| c_2 & - c_3 s_2 & s_2 s_3 \\
| |
| c_1 s_2 & c_1 c_2 c_3 - s_1 s_3 & - c_3 s_1 - c_1 c_2 s_3 \\
| |
| s_1 s_2 & c_1 s_3 + c_2 c_3 s_1 & c_1 c_3 - c_2 s_1 s_3
| |
| \end{bmatrix}</math>
| |
| |<math>X_1 Z_2 Y_3 = \begin{bmatrix}
| |
| c_2 c_3 & - s_2 & c_2 s_3 \\
| |
| s_1 s_3 + c_1 c_3 s_2 & c_1 c_2 & c_1 s_2 s_3 - c_3 s_1 \\
| |
| c_3 s_1 s_2 - c_1 s_3 & c_2 s_1 & c_1 c_3 + s_1 s_2 s_3
| |
| \end{bmatrix}</math>
| |
| |-
| |
| |<math>X_1 Y_2 X_3 = \begin{bmatrix}
| |
| c_2 & s_2 s_3 & c_3 s_2 \\
| |
| s_1 s_2 & c_1 c_3 - c_2 s_1 s_3 & - c_1 s_3 - c_2 c_3 s_1 \\
| |
| - c_1 s_2 & c_3 s_1 + c_1 c_2 s_3 & c_1 c_2 c_3 - s_1 s_3
| |
| \end{bmatrix}</math>
| |
| |<math>X_1 Y_2 Z_3 = \begin{bmatrix}
| |
| c_2 c_3 & - c_2 s_3 & s_2 \\
| |
| c_1 s_3 + c_3 s_1 s_2 & c_1 c_3 - s_1 s_2 s_3 & - c_2 s_1 \\
| |
| s_1 s_3 - c_1 c_3 s_2 & c_3 s_1 + c_1 s_2 s_3 & c_1 c_2
| |
| \end{bmatrix}</math>
| |
| |-
| |
| |<math>Y_1 X_2 Y_3 = \begin{bmatrix}
| |
| c_1 c_3 - c_2 s_1 s_3 & s_1 s_2 & c_1 s_3 + c_2 c_3 s_1 \\
| |
| s_2 s_3 & c_2 & - c_3 s_2 \\
| |
| - c_3 s_1 - c_1 c_2 s_3 & c_1 s_2 & c_1 c_2 c_3 - s_1 s_3
| |
| \end{bmatrix}</math>
| |
| |<math>Y_1 X_2 Z_3 = \begin{bmatrix}
| |
| c_1 c_3 + s_1 s_2 s_3 & c_3 s_1 s_2 - c_1 s_3 & c_2 s_1 \\
| |
| c_2 s_3 & c_2 c_3 & - s_2 \\
| |
| c_1 s_2 s_3 - c_3 s_1 & s_1 s_3 + c_1 c_3 s_2 & c_1 c_2
| |
| \end{bmatrix}</math>
| |
| |-
| |
| |<math>Y_1 Z_2 Y_3 = \begin{bmatrix}
| |
| c_1 c_2 c_3 - s_1 s_3 & - c_1 s_2 & c_3 s_1 + c_1 c_2 s_3 \\
| |
| c_3 s_2 & c_2 & s_2 s_3 \\
| |
| - c_1 s_3 - c_2 c_3 s_1 & s_1 s_2 & c_1 c_3 - c_2 s_1 s_3
| |
| \end{bmatrix}</math>
| |
| |<math>Y_1 Z_2 X_3 = \begin{bmatrix}
| |
| c_1 c_2 & s_1 s_3 - c_1 c_3 s_2 & c_3 s_1 + c_1 s_2 s_3 \\
| |
| s_2 & c_2 c_3 & - c_2 s_3 \\
| |
| - c_2 s_1 & c_1 s_3 + c_3 s_1 s_2 & c_1 c_3 - s_1 s_2 s_3
| |
| \end{bmatrix}</math>
| |
| |-
| |
| |<math>Z_1 Y_2 Z_3 = \begin{bmatrix}
| |
| c_1 c_2 c_3 - s_1 s_3 & - c_3 s_1 - c_1 c_2 s_3 & c_1 s_2 \\
| |
| c_1 s_3 + c_2 c_3 s_1 & c_1 c_3 - c_2 s_1 s_3 & s_1 s_2 \\
| |
| - c_3 s_2 & s_2 s_3 & c_2
| |
| \end{bmatrix}</math>
| |
| |<math>Z_1 Y_2 X_3 = \begin{bmatrix}
| |
| c_1 c_2 & c_1 s_2 s_3 - c_3 s_1 & s_1 s_3 + c_1 c_3 s_2 \\
| |
| c_2 s_1 & c_1 c_3 + s_1 s_2 s_3 & c_3 s_1 s_2 - c_1 s_3 \\
| |
| - s_2 & c_2 s_3 & c_2 c_3
| |
| \end{bmatrix}</math>
| |
| |-
| |
| |<math>Z_1 X_2 Z_3 = \begin{bmatrix}
| |
| c_1 c_3 - c_2 s_1 s_3 & - c_1 s_3 - c_2 c_3 s_1 & s_1 s_2 \\
| |
| c_3 s_1 + c_1 c_2 s_3 & c_1 c_2 c_3 - s_1 s_3 & - c_1 s_2 \\
| |
| s_2 s_3 & c_3 s_2 & c_2
| |
| \end{bmatrix}</math>
| |
| |<math>Z_1 X_2 Y_3 = \begin{bmatrix}
| |
| c_1 c_3 - s_1 s_2 s_3 & - c_2 s_1 & c_1 s_3 + c_3 s_1 s_2 \\
| |
| c_3 s_1 + c_1 s_2 s_3 & c_1 c_2 & s_1 s_3 - c_1 c_3 s_2 \\
| |
| - c_2 s_3 & s_2 & c_2 c_3
| |
| \end{bmatrix}</math>
| |
| |}
| |
| | |
| To change the formulas for the opposite direction of rotation, change the signs of the sine functions. To change the formulas for [[Active and passive transformation|passive rotations]], transpose the matrices (then each matrix transforms the initial coordinates of a vector remaining fixed to the coordinates of the same vector measured in the rotated reference system; same rotation axis, same angles, but now the coordinate system rotates, rather than the vector).
| |
| | |
| ===Quaternions===
| |
| | |
| [[quaternions|Unit quaternions]], also known as [[Euler–Rodrigues parameters]], provide [[Quaternions and spatial rotation|another mechanism]] for representing 3D rotations. This is equivalent to the special unitary group description.
| |
| | |
| Expressing rotations in 3D as unit quaternions instead of matrices has some advantages:
| |
| * Concatenating rotations is computationally faster and numerically more stable.
| |
| * Extracting the angle and axis of rotation is simpler.
| |
| * Interpolation is more straightforward. See for example [[slerp]].
| |
| | |
| ===Geometric algebra===
| |
| Other representation comes from the [[Geometric algebra]](GA). GA is a higher level abstraction, in which the quaternions are an even subalgebra. The principal tool in GA is the rotor <math> \bold\R=[\cos(\theta/2)-I u \sin(\theta/2)] </math> where <math>\bold\theta =</math>[[angle of rotation]], <math>\bold(u) =</math>rotation axis (unitary vector) and <math>\bold(I)=</math>pseudoscalar (trivector in <math>\mathbb{R}^3</math>)
| |
| | |
| == Properties ==
| |
| {{see also|Charts on SO(3)|Quaternions and spatial rotation}}
| |
| The Euler angles form a [[chart (topology)|chart]] on all of [[SO(3)]], the [[special orthogonal group]] of rotations in 3D space. The chart is smooth except for a polar coordinate style singularity along β=0. See [[charts on SO(3)]] for a more complete treatment.
| |
| | |
| The space of rotations is called in general "The [[Quaternions and spatial rotation#Quaternion rotation operations|Hypersphere of rotations]]", though this is a misnomer: the group [[Spin group|Spin(3)]] is [[isometric embedding|isometric]] to the hypersphere ''S''<sup>3</sup>, but the rotation space SO(3) is instead isometric to the [[real projective space]] '''RP'''<sup>3</sup> which is a 2-fold [[quotient space]] of the hypersphere. This 2-to-1 ambiguity is the mathematical origin of [[Spin (physics)|spin in physics]].
| |
| | |
| A similar three angle decomposition applies to [[SU(2)]], the [[special unitary group]] of rotations in complex 2D space, with the difference that β ranges from 0 to 2π. These are also called Euler angles.
| |
| | |
| The [[Haar measure]] for Euler angles has the simple form sin(β).dα.dβ.dγ, usually normalized by a factor of 1/8π².
| |
| | |
| For example, to generate uniformly randomized orientations, let α and γ be uniform from 0 to 2π, let ''z'' be uniform from −1 to 1, and let β = arccos(''z'').
| |
| | |
| ==Higher dimensions==
| |
| | |
| It is possible to define parameters analogous to the Euler angles in dimensions higher than three.<ref>{{it icon}} [http://ansi.altervista.org A generalization of Euler Angles to ''n''-dimensional real spaces]</ref>
| |
| | |
| The number of degrees of freedom of a rotation matrix is always less than the dimension of the matrix squared. That is, the elements of a rotation matrix are not all completely independent. For example, the rotation matrix in dimension 2 has only one degree of freedom, since all four of its elements depend on a single angle of rotation. A rotation matrix in dimension 3 (which has nine elements) has three degrees of freedom, corresponding to each independent rotation, for example by its three Euler angles or a magnitude one (unit) quaternion.
| |
| | |
| In [[SO(4)]] the rotation matrix [[Rotation (mathematics)#In four dimensions|is defined by two quaternions]], and is therefore 6-parametric (three degrees of freedom for every quaternion). The 4x4 rotation matrices have therefore 6 out of 16 independent components.
| |
| | |
| Any set of 6 parameters that define the rotation matrix could be considered an extension of Euler angles to dimension 4.
| |
| | |
| In general, the number of euler angles in dimension D is quadratic in D; since any one rotation consists of choosing two dimensions to rotate between, the total number of rotations available in dimension <math>D</math> is <math>N_{rot}=\binom{D}{2}=D(D-1)/2</math>, which for <math>D=2,3,4</math> yields <math>N_{rot}=1,3,6</math>.
| |
| | |
| == Applications ==
| |
| [[image:Gyroscope operation.gif|thumb|right|A [[gyroscope]] keeps its rotation axis constant. Therefore, angles measured in this frame are equivalent to angles measured in the lab frame]]
| |
| | |
| ===Vehicles and moving frames===
| |
| {{main|rigid body}}
| |
| {{see also|axes conventions}}
| |
| | |
| Their main advantage over other orientation descriptions is that they are directly measurable from a gimbal mounted in a vehicle. As gyroscopes keep their rotation axis constant, angles measured in a gyro frame are equivalent to angles measured in the lab frame. Therefore gyros are used to know the actual orientation of moving spacecraft, and Euler angles are directly measurable. Intrinsic rotation angle cannot be read from a single gimbal, so there has to be more than one gimbal in a spacecraft. Normally there are at least three for redundancy. There is also a relation to the well-known [[gimbal lock]] problem of [[Mechanical Engineering]] <ref>The relation between the Euler angles and the Cardan suspension is explained in chap. 11.7 of the following textbook: U. Krey, A. Owen, ''Basic Theoretical Physics – A Concise Overview, New York, London, Berlin, Heidelberg, Springer (2007) .</ref> .
| |
| | |
| [[Image:plane.svg|thumb|left|200px|Heading, elevation and bank for an aircraft with axes [[Deutsches Institut für Normung|DIN]] 9300]]
| |
| The most popular application is to describe aircraft attitudes, normally using a Tait–Bryan convention so that zero degrees elevation represents the horizontal attitude. Tait–Bryan angles represent the orientation of the aircraft respect a reference axis system (''world frame'') with three angles which in the context of an aircraft are normally called Heading, Elevation and Bank. When dealing with vehicles, different [[axes conventions]] are possible.
| |
| | |
| When studying rigid bodies in general, one calls the ''xyz'' system ''space coordinates'', and the ''XYZ'' system ''body coordinates''. The space coordinates are treated as unmoving, while the body coordinates are considered embedded in the moving body. Calculations involving [[acceleration]], [[angular acceleration]], [[angular velocity]], [[angular momentum]], and [[kinetic energy]] are often easiest in body coordinates, because then the moment of inertia tensor does not change in time. If one also diagonalizes the rigid body's moment of inertia tensor (with nine components, six of which are independent), then one has a set of coordinates (called the principal axes) in which the moment of inertia tensor has only three components.
| |
| | |
| The angular velocity of a rigid body takes a [[Angular velocity#Components_from_Euler_angles|simple form]] using Euler angles in the moving frame. Also the [[Euler's equations (rigid body dynamics)|Euler's rigid body equations]] are simpler because the inertia tensor is constant in that frame.
| |
| | |
| [[File:Automation of foundry with robot.jpg|thumb|right|Industrial robot operating in a foundry.]]
| |
| | |
| ===Others===
| |
| Euler angles, normally in the Tait–Bryan convention, are also used in [[Industrial robot|robotics]] for speaking about the degrees of freedom of a [[Robotic arm|wrist]]. They are also used in [[Electronic stability control]] in a similar way.
| |
| | |
| Gun fire control systems require corrections to gun-order angles (bearing and elevation) to compensate for deck tilt (pitch and roll). In traditional systems, a stabilizing gyroscope with a vertical spin axis corrects for deck tilt, and stabilizes the optical sights and radar antenna. However, gun barrels point in a direction different from the line of sight to the target, to anticipate target movement and fall of the projectile due to gravity, among other factors. Gun mounts roll and pitch with the deck plane, but also require stabilization. Gun orders include angles computed from the vertical gyro data, and those computations involve Euler angles.
| |
| | |
| Euler angles are also used extensively in the quantum mechanics of angular momentum. In quantum mechanics, explicit descriptions of the representations of ''SO(3)'' are very important for calculations, and almost all the work has been done using Euler angles. In the early history of quantum mechanics, when physicists and chemists had a sharply negative reaction towards abstract group theoretic methods (called the ''Gruppenpest''), reliance on Euler angles was also essential for basic theoretical work.
| |
| | |
| In materials science, crystallographic [[texture (crystalline)|texture]] (or preferred orientation) can be described using Euler angles. In texture analysis, the Euler angles provide the necessary mathematical depiction of the orientation of individual crystallites within a polycrystalline material, allowing for the quantitative description of the macroscopic material.<ref>{{Citation
| |
| |last1 = Kocks
| |
| |first1 = U.F.
| |
| |last2 = Tomé
| |
| |first2 = C.N.
| |
| |last3 = Wenk
| |
| |first3 = H.-R.
| |
| |title = Texture and Anisotropy: Preferred Orientations in Polycrystals and their effect on Materials Properties
| |
| |publisher = [[Cambridge]]
| |
| |year = 2000
| |
| |isbn = 978-0-521-79420-6}}</ref>
| |
| The most common definition of the angles is due to Bunge and corresponds to the ZXZ convention. It is important to note, however, that the application generally involves axis transformations of tensor quantities, i.e. passive rotations. Thus the matrix that corresponds to the Bunge Euler angles is the transpose of that shown in the table above.<ref>{{Citation
| |
| |last = Bunge
| |
| |first = H.
| |
| |title = Texture Analysis in Materials Science: Mathematical Methods
| |
| |publisher = [[CUVILLIER VERLAG]]
| |
| |year = 1993
| |
| |asin = B0014XV9HU}}</ref>
| |
| | |
| Many mobile computing devices contain [[accelerometer]]s which can determine these devices' Euler angles with respect to the earth's gravitational attraction. These are used in applications such as games, [[bubble level]] simulations, and [[kaleidoscope]]s.{{citation needed|reason=accelerometers can measure orientation, yes; is there any evidence they use Euler angles for this?|date=May 2011}}
| |
| | |
| == See also ==
| |
| * [[3D projection]]
| |
| * [[Axis-angle representation]]
| |
| * [[Conversion between quaternions and Euler angles]]
| |
| * [[Euler's rotation theorem]]
| |
| * [[Quaternion]]
| |
| * [[Quaternions and spatial rotation]]
| |
| * [[Rotation formalisms in three dimensions]]
| |
| * [[Spherical coordinate system]]
| |
| | |
| ==References==
| |
| {{Reflist}}<!--added under references heading by script-assisted edit-->
| |
| | |
| ==Bibliography==
| |
| * {{citation
| |
| |last1= Biedenharn
| |
| |first1= L. C.
| |
| |last2= Louck
| |
| |first2= J. D.
| |
| |title= Angular Momentum in Quantum Physics
| |
| |publisher= [[Addison–Wesley]]
| |
| |place= Reading, MA
| |
| |year= 1981
| |
| |isbn= 978-0-201-13507-7
| |
| }}
| |
| * {{citation
| |
| |last1= Goldstein
| |
| |first1= Herbert
| |
| |author1-link= Herbert Goldstein
| |
| |title= Classical Mechanics
| |
| |edition= 2nd
| |
| |publisher= Addison–Wesley
| |
| |place= Reading, MA
| |
| |year= 1980
| |
| |isbn= 978-0-201-02918-5
| |
| }}
| |
| * {{citation
| |
| |last1= Gray
| |
| |first1= Andrew
| |
| |title= A Treatise on Gyrostatics and Rotational Motion
| |
| |publisher= [[Macmillan Publishers|Macmillan]]
| |
| |place= London
| |
| |year= 1918
| |
| |publication-date= 2007
| |
| |isbn= 978-1-4212-5592-7
| |
| }}
| |
| * {{citation
| |
| |last1= Rose
| |
| |first1= M. E.
| |
| |title= Elementary Theory of Angular Momentum
| |
| |publisher= [[John Wiley & Sons]]
| |
| |place= New York, NY
| |
| |year= 1957
| |
| |publication-date= 1995
| |
| |isbn= 978-0-486-68480-2
| |
| }}
| |
| * {{citation
| |
| |last1= Symon
| |
| |first1= Keith
| |
| |title= Mechanics
| |
| |publisher= Addison-Wesley
| |
| |place= Reading, MA
| |
| |year= 1971
| |
| |isbn= 978-0-201-07392-8
| |
| }}
| |
| * {{citation
| |
| |last1= Landau
| |
| |first1= L.D.
| |
| |author1-link=Lev Landau
| |
| |last2= Lifshitz
| |
| |first2= E. M.
| |
| |author2-link= Evgeny Lifshitz
| |
| |title= Mechanics
| |
| |edition= 3rd
| |
| |publisher= Butterworth-Heinemann
| |
| |year= 1996
| |
| |isbn= 978-0-7506-2896-9
| |
| |location= Oxford
| |
| }}
| |
| | |
| ==External links==
| |
| {{Commons|Euler angles}}
| |
| | |
| *{{MathWorld|urlname=EulerAngles|title=Euler Angles}}
| |
| * Java applet for the simulation of Euler angles available at http://www.parallemic.org/Java/EulerAngles.html.
| |
| * [https://itunes.apple.com/us/app/eulerangles/id584911325 EulerAngles] - An iOS app for visualizing in 3D the three rotations associated with Euler angles.
| |
| * http://sourceforge.net/projects/orilib – A collection of routines for rotation / orientation manipulation, including special tools for crystal orientations.
| |
| * Online tool to compose rotation matrices available at http://www.vectoralgebra.info/eulermatrix.html
| |
| | |
| {{DEFAULTSORT:Euler Angles}}
| |
| [[Category:Rotation in three dimensions]]
| |
| [[Category:Euclidean symmetries]]
| |
| [[Category:Angle]]
| |
| [[Category:Analytic geometry]]
| |