PlanetPhysics/Direction Cosine Matrix

From testwiki
Jump to navigation Jump to search

A direction cosine matrix (DCM) is a transformation [[../Matrix/|matrix]] that transforms one coordinate [[../CosmologicalConstant2/|reference frame]] to another. If we extend the [[../PreciseIdea/|concept]] of how the three dimensional [[../DirectionCosines/|direction cosines]] locate a [[../Vectors/|vector]], then the DCM locates three [[../PureState/|unit vectors]] that describe a coordinate reference frame. Using the notation in equation 1, we need to find the matrix elements that correspond to the correct transformation matrix.

DCM=[A11A12A13A21A22A23A31A32A33]

The first unit vector of the second coordinate frame can be located in the first frame by normal vector notation. See figure 1 for relationship.

y^1=A11x^1+A12x^2+A13x^3

\medskip \begin{figure} \includegraphics[scale=0.78]{DCM.eps} \end{figure} \medskip

Similarily, the other two unit vectors can be described by

y^2=A21x^1+A22x^2+A23x^3 y^3=A31x^1+A32x^2+A33x^3

It is easy to see how equation 1 works as a transformation matrix through simple [[../Matrix/|matrix multiplication]].

[y^1y^2y^3]=[A11A12A13A21A22A23A31A32A33][x^1x^2x^3]

Once this transformation matrix is found, it can be used to transform vectors from the second frame to the first frame and vice versa. Equation 2 transforms the x frame to the y frame and can be denoted as R12. In order to get R21, which transforms the y frame to the x frame, we use a property of transformation matrices of orthonormal reference frames (a frame that is described by unit vectors and are perpindicular to each other). See the entry on a transformation matrix for more info on its properties. We use the properties that

R121=R12T=R21 R12R12T=[100010001]

so using these properties and rearranging equation 2 y^=R12x^ yields

R121y^=R121R12x^

giving the transformation of the y frame to the x frame

x^=R21y^

So to extend this concept to transform vectors from one frame to another a closer examination of a vector being represented in both frames is needed. If we denote the second frame as the prime () frame, then a vector expressed in each of these is given by

v=v1x^1+v2x^2+v3x^3
v=v1y^1+v2y^2+v3y^3

Since both equations describe the same vector, let us set them equal to each other so

v1x^1+v2x^2+v3x^3=v1y^1+v2y^2+v3y^3

This notation is clumsy so we want to represent it in matrix notation. This is simple enough if you have an understanding of multiplying a column vector by a row vector. This allows us to describe equations 3 and 4 by

v=[v1v2v3][x^1x^2x^3] v=[v1v2v3][y^1y^2y^3]

Setting them equal and substituting equation 2 in for the second coordinate frame yields

v=[v1v2v3][x^1x^2x^3]=[v1v2v3][A11A12A13A21A22A23A31A32A33][x^1x^2x^3]

Then by inspection (or go through the matrix manipulation to cancel the x frame)

[v1v2v3]=[v1v2v3][A11A12A13A21A22A23A31A32A33]

Representing the transformation matrix as R12 as the transformation from the first frame to the second frame and transposing the previous equation gives

[v1v2v3]=([v1v2v3]R12)T

Performing the transposition and using a transposition property for two matrices A and B such that

(AB)T=BTAT

leads to the relationship

[v1v2v3]=R12T[v1v2v3]

Finally giving us the ability to transform a vector from the second (prime) frame to the first frame.

v=R21v

Much much more can be found in the general entry about the Transformation matrix.

Template:CourseCat