Physics Library
 An open source physics library
Encyclopedia | Forums | Docs | Random |  
Login
create new user
Username:
Password:
forget your password?
Main Menu
Sections

Meta

Talkback

Downloads

Information
direction cosine matrix to Euler 232 angles (Definition)

Starting with a direction cosine matrix (DCM), we need to determine the three Euler Angles. The connection is made by comparing the DCM elements with the combined Euler 232 sequence. It is important to note that the 12 combinations of Euler angles for a given sequence can be found from a given DCM. The DCM matrix is

         ⌊                ⌋
           A11  A12   A13
DCM   =  ⌈ A21  A22   A23 ⌉
           A31  A32   A33
(1)

The Euler 232 sequence is

                    ⌊                                       ⌋
                       cψc 𝜃cϕ − sψs ϕ  cψs𝜃  − cψc 𝜃sϕ − s ψcϕ
R  (ψ )R (𝜃)R  (ϕ) = ⌈     − s c        c         s s        ⌉
  2     3    2               𝜃 ϕ        𝜃         𝜃 ϕ
                       sψc𝜃cϕ + cψs ϕ sψs𝜃  − sψc𝜃sϕ + cψc ϕ
(2)

If we examine the element in coloum 2 row 2, then by inspection

A22 = cos(𝜃)

Solving for 𝜃 yields

        −1
𝜃 = cos   (A22)
(3)

Care must now be taken when evaluating the inverse cosine. It is a multivalued function. Analytically, the convention is to choose the principle value such that

0 𝜃 π

If a numerical program is used, a function acos() usually does this for us. The next step is to analyze the ratio A23-
A
 21. Using these values from the Euler sequence we get

A23
----
A21 =  sin(𝜃)sin(ϕ )
--------------
− sin(𝜃)cos(ϕ)

Rearranging and using the tangent yields

tan(ϕ) =  A
--23--
− A21

Solving the quadrant ambiquity caused by the inverse tangent is done by examining the signs of the numerator and denominator. Denoting y as the numerator and x as the denominator, then the quadrant is chosen by:

quadrant 1 [        ]
   y  x
  +   +

so if x and y are both positive, then 0 ϕ π∕2. Similarily for the other quadrant possibilites

quadrant 2 [  y  x  ]

  +   −

quadrant 3 [        ]
   y  x
  −   −

quadrant 4 [        ]
   y  x
  −   +

Ofcourse, it is much simplier to use a calculator or numerical program that uses the atan2() function which will choose the principle value. In a similar fashion we can find the final Euler angle by looking at the ratio A32-
A12 which gives the relation A32-
A12 = sin(ψ)sin(𝜃)-
cos(ψ)sin(𝜃)

Canceling terms and rearrangeing gives us

ψ = tan1(A32-
A12)

which uses the same method to resolove quadrant ambiguity as above. To summarize, we will give the formulas for the conversion from a direction cosine matrix to the Euler 232 angles in Matlab syntax. Be careful how you implement this in other numerical programs. In Matlab it goes atan2(y,x) and in Mathematica it is ArcTan[x, y].

ϕ = atan2 (A23,− A21)

𝜃 = acos(A22)

ψ = atan2 (A32,A12 )

"direction cosine matrix to Euler 232 angles" is owned by bloftin.
(view preamble)
View style:

Cross-references: formulas, relation, program, function, matrix, Euler 232 sequence, Euler Angles, direction cosine matrix

This is version 4 of direction cosine matrix to Euler 232 angles, born on 2005-08-17, modified 2005-08-25.
Object id is 70, canonical name is DirectionCosineMatrixToEuler232Sequence.
Accessed 2740 times total.

Classification:
Physics Classification45.40.-f (Dynamics and kinematics of rigid bodies)
Pending Errata and Addenda
None.
Discussion
Style: Expand: Order:

No messages.

Interact
rate | post | correct | update request | add derivation | add example | add (any)