Main Content

dcm2rod

Convert direction cosine matrix to Euler-Rodrigues vector

collapse all in page

Syntax

R = dcm2rod(dcm)
R = dcm2rod(dcm,action)
R = dcm2rod(dcm,action,tolerance)

Description

R = dcm2rod(dcm) function calculates the Euler-Rodrigues vector (R) from the direction cosine matrix. This function applies only to direction cosine matrices that are orthogonal with determinant +1. The direction cosine matrix input and resulting Euler-Rodrigues vector represent a right-hand passive transformation from frame A to frame B.

R = dcm2rod(dcm,action) performs action if the direction cosine matrix is invalid (not orthogonal).

R = dcm2rod(dcm,action,tolerance) uses a tolerance level to evaluate if the direction cosine matrix, n, is valid (orthogonal).

example

Examples

collapse all

Open Live Script

This example shows how to determine the Rodrigues vector from the direction cosine matrix.

DCM = [0.433 0.75 0.5;-0.25 -0.433 0.866;0.866 -0.5 0.0];
r = dcm2rod(DCM)
r = 1×3

    1.3660    0.3660    1.0000
Open Live Script

Determine the Rodrigues vector from the direction cosine matrix validated within tolerance.

DCM = [0.433 0.75 0.5;-0.25 -0.433 0.866;0.866 -0.5 0.0];
r = dcm2rod(DCM,'Warning',0.1)
r = 1×3

    1.3660    0.3660    1.0000

Input Arguments

collapse all

3-by-3-by-M containing M direction cosine matrices.

Data Types: double

Function behavior when direction cosine matrix is invalid (not orthogonal).

  • Warning — Displays warning and indicates that the direction cosine matrix is invalid.

  • Error — Displays error and indicates that the direction cosine matrix is invalid.

  • None — Does not display warning or error (default).

Data Types: char | string

Tolerance of direction cosine matrix validity, specified as a scalar. The function considers the direction cosine matrix valid if these conditions are true:

  • The transpose of the direction cosine matrix times itself equals 1 within the specified tolerance (transpose(n)*n == 1±tolerance)

  • The determinant of the direction cosine matrix equals 1 within the specified tolerance (det(n) == 1±tolerance).

Data Types: double

Output Arguments

collapse all

M-by-3 matrix containing M Rodrigues vectors.

Data Types: double

Algorithms

An Euler-Rodrigues vector b⇀ represents a rotation by integrating a direction cosine of a rotation axis with the tangent of half the rotation angle as follows:

b→=[bxbybz]

where:

bx=tan(12θ)sx,by=tan(12θ)sy,bz=tan(12θ)sz

are the Rodrigues parameters. Vector s⇀ represents a unit vector around which the rotation is performed. Due to the tangent, the rotation vector is indeterminate when the rotation angle equals ±pi radians or ±180 deg. Values can be negative or positive.

References

[1] Dai, J.S. "Euler-Rodrigues formula variations, quaternion conjugation and intrinsic connections." Mechanism and Machine Theory, 92, 144-152. Elsevier, 2015.

Version History

Introduced in R2017a

See Also

| | | |