PlanetMath (more info)
 Math for the people, by the people. Sponsor PlanetMath
Encyclopedia | Requests | Forums | Docs | Wiki | Random | RSS  
Login
create new user
name:
pass:
forget your password?
Main Menu
Owner confidence rating: High Entry average rating: No information on entry rating
[parent] proof of Rodrigues' rotation formula (Proof)

Let $[ \bx, \by, \bz]$ be a frame of right-handed orthonormal vectors in $\real^3$ , and let $\bv = a\bx + b\by + c\bz$ (with $a, b, c \in \real$ ) be any vector to be rotated on the $\bz$ axis, by an angle $\theta$ counter-clockwise.

\includegraphics{rodrigues.eps}

The image vector $\bv'$ is the vector $\bv$ with its component in the $\bx,\by$ plane rotated, so we can write $$ \bv' = a\bx' + b\by' + c\bz\,,\\ $$ where $\bx'$ and $\by'$ are the rotations by angle $\theta$ of the $\bx$ and $\by$ vectors in the $\bx,\by$ plane. By the rotation formula in two dimensions, we have

$\displaystyle \mathbf{x}'$ $\displaystyle = \cos \theta \mathbf{x}+ \sin \theta \mathbf{y} ,$    
$\displaystyle \mathbf{y}'$ $\displaystyle = -\sin \theta \mathbf{x}+ \cos \theta \mathbf{y} .$    

So $$ \bv' = \cos \theta (a\bx + b\by) + \sin \theta (a \by - b\bx) + c\bz\,. $$ The vector $a \bx + b\by$ is the projection of $\bv$ onto the $\bx,\by$ plane, and $a \by - b\bx$ is its rotation by $90^\circ$ . So these two vectors form an orthogonal frame in the $\bx,\by$ plane, although they are not necessarily unit vectors. Alternate expressions for these vectors are easily derived -- especially with the help of the picture:
$\displaystyle \mathbf{v}- (\mathbf{v}\cdot \mathbf{z}) \mathbf{z}$ $\displaystyle = \mathbf{v}- c\mathbf{z}= a\mathbf{x}+ b\mathbf{y} ,$    
$\displaystyle \mathbf{z}\times \mathbf{v}$ $\displaystyle = a (\mathbf{z}\times \mathbf{x}) + b (\mathbf{z}\times \mathbf{y}) + c (\mathbf{z}\times \mathbf{z}) = a\mathbf{y}- b\mathbf{x} .$    

Substituting these into the expression for $\bv'$ : $$ \bv' = \cos \theta (\bv - (\bv \cdot \bz) \bz) + \sin\theta (\bz \times \bv) + c\bz\,, $$ which could also have been derived directly if we had first considered the frame $[\bv - (\bv \cdot \bz), \bz \times \bv]$ instead of $[\bx, \by]$ .

We attempt to simplify further: $$ \bv' = \bv + \sin\theta (\bz \times \bv) + (\cos \theta - 1)(\bv - (\bv \cdot \bz)\bz )\,. $$ Since $\bz \times \bv$ is linear in $\bv$ , this transformation is represented by a linear operator $A$ . Under a right-handed orthonormal basis, the matrix representation of $A$ is directly computed to be $$ A\bv = \bz \times \bv = \begin{bmatrix} 0 & -z_3 & z_2 \\ z_3 & 0 & -z_1 \\ -z_2 & z_1 & 0 \end{bmatrix} \, \begin{bmatrix} v_1 \\ v_2 \\ v_3 \end{bmatrix} \,. $$ We also have

$\displaystyle -(\mathbf{v}- (\mathbf{v}\cdot \mathbf{z}) \mathbf{z})$ $\displaystyle = -a\mathbf{x}- b\mathbf{y}$ (rotate $ a\mathbf{x}+ b\mathbf{y}$ by $ 180^\circ$)    
  $\displaystyle = \mathbf{z}\times (a\mathbf{y}- b\mathbf{x})$ (rotate $ a\mathbf{y}- b\mathbf{x}$ by $ 90^\circ$)    
  $\displaystyle = \mathbf{z}\times (\mathbf{z}\times (a\mathbf{x}+ b\mathbf{y}+c\mathbf{z}))$    
  $\displaystyle = A^2 \mathbf{v} .$    

So $$ \bv' = I\bv + \sin \theta \, A\bv + (1- \cos\theta) A^2 \, \bv\,, $$ proving Rodrigues' rotation formula.

Relation with the matrix exponential

Here is a curious fact. Notice that the matrix 1 $A$ is skew-symmetric. This is not a coincidence -- for any skew-symmetric matrix $B$ , we have $\transpose{(e^B)} = e^{\transpose{B}} = e^{-B} = (e^B)^{-1}$ , and $\det e^B = e^{\trace B} = e^0 = 1$ , so $e^B$ is always a rotation. It is in fact the case that:
$\displaystyle I + \sin \theta A + (1- \cos\theta) A^2 = e^{\theta A}$    

for the matrix $A$ we had above! To prove this, observe that powers of $A$ cycle like so:
$\displaystyle I, A, A^2, -A, -A^2, A, A^2, -A, -A^2, \dotsc$    

Then
$\displaystyle \sin \theta A$ $\displaystyle = \sum_{k=0}^\infty \frac{(-1)^k \theta^{2k+1} A}{(2k+1)!} = \sum... ...heta^{2k+1} A^{2k+1}}{(2k+1)!} = \sum_{k \textrm{ odd}} \frac{(\theta A)^k}{k!}$    
$\displaystyle (1 - \cos \theta) A^2$ $\displaystyle = \sum_{k=1}^\infty \frac{(-1)^{k-1} \theta^{2k} A^2}{(2k)!} = \s... ...^{2k} A^{2k}}{(2k)!} = \sum_{k\geq 2 \textrm{ even}} \frac{(\theta A)^k}{k!} .$    

Adding $\sin \theta \, A$ , $(1-\cos \theta) A^2$ and $I$ together, we obtain the power series for $e^{\theta A}$ .

Second proof: If we regard $\theta$ as time, and differentiate the equation $\bv' = a\bx' + b\by' + c\bz$ with respect to $\theta$ , we obtain $d\bv'/d\theta = a\by' - b\bx' = \bz \times \bv' = A \bv'$ , whence the solution (to this linear ODE) is $\bv' = e^{\theta A} \bv$ .

Remark: The operator $e^{\theta A}$ , as $\theta$ ranges over $\real$ , is a one-parameter subgroup of $\mathrm{SO}(3)$ . In higher dimensions $n$ , every rotation in $\mathrm{SO}(n)$ is of the form $e^A$ for a skew-symmetric $A$ , and the second proof above can be modified to prove this more general fact.



Footnotes

...http://planetmath.org/encyclopedia/Order7.html 1
If we want to use coordinate-free language, then in this section, ``matrix'' should be replaced by ``linear operator'' and transposes should be replaced by the adjoint operation.



"proof of Rodrigues' rotation formula" is owned by stevecheng.
(view preamble | get metadata)

View style:

See Also: dimension of the special orthogonal group


This object's parent.
Log in to rate this entry.
(view current ratings)

Cross-references: one-parameter subgroup, ranges, operator, ODE, solution, equation, differentiate, proof, power series, cycle, powers, skew-symmetric matrix, skew-symmetric, operation, adjoint, transposes, section, matrix, Rodrigues rotation formula, matrix representation, orthonormal basis, linear operator, transformation, expressions, unit vectors, orthogonal, onto, projection, dimensions, formula, rotations, plane, component, image, angle, axis, vectors, orthonormal, right-handed, frame

This is version 11 of proof of Rodrigues' rotation formula, born on 2005-07-12, modified 2006-06-16.
Object id is 7221, canonical name is ProofOfRodriguesRotationFormula.
Accessed 7324 times total.

Classification:
AMS MSC15-00 (Linear and multilinear algebra; matrix theory :: General reference works )
 51-00 (Geometry :: General reference works )

Pending Errata and Addenda
None.
[ View all 1 ]
Discussion
Style: Expand: Order:
forum policy

No messages.

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