Notación axial-angular

El ángulo Plantilla:Math y el vector unitario axial Plantilla:Math definen una rotación, concisamente representada por el vector de rotación Plantilla:Math

En matemáticas, la notación axial-angular de una rotación parametriza una rotación en el espacio euclídeo tridimensional mediante dos elementos: un vector unitario Plantilla:Math que indica la dirección de un eje de rotación, y un ángulo Plantilla:Math que describe la magnitud de la rotación respecto al eje. Solo se necesitan dos números, no tres, para definir la dirección de un vector unitario Plantilla:Math ubicado en el origen, porque Plantilla:Math tiene necesariamente módulo 1, y expresadas dos coordenadas, se puede deducir la tercera (dado que $x^2+y^2+z^2=1$, y por lo tanto, $z=\sqrt{x^2+y^2}$). También es posible utilizar los dos ángulos que definen la elevación y el azimut de Plantilla:Math, suficientes para ubicarlo en cualquier marco de coordenadas cartesianas en particular.

De acuerdo con la fórmula de rotación de Rodrigues, el ángulo y el eje determinan una transformación que determina el giro de vectores tridimensionales. La rotación se produce en el sentido prefijado por la regla de la mano derecha. El eje de rotación a veces se denomina eje de Euler.

Es una de las muchas formalizaciones de la rotación en tres dimensiones. La notación axial-angular se basa en el Teorema de rotación de Euler, que asegura que cualquier rotación o secuencia de rotaciones de un cuerpo rígido en un espacio tridimensional es equivalente a una rotación pura sobre un solo eje fijo.

Vector de rotación

La notación axial-angular es equivalente al más conciso vector de rotación, también llamado vector de Euler. En este caso, tanto el eje de rotación como el ángulo están representados por un vector codireccional con el eje de rotación cuyo módulo (longitud) coincide con el ángulo de rotación Plantilla:Math,

$\mathit{\theta }=\theta 𝐞.$

Se utiliza en las aplicaciones exponenciales y logarítmicas, que involucran este tipo de notación.

Debe tenerse en cuenta que distintos vectores de rotación corresponden a la misma rotación. En particular, un vector de rotación de longitud Plantilla:Math, para cualquier entero Plantilla:Math, codifica exactamente la misma rotación que un vector de rotación de longitud Plantilla:Math. Por lo tanto, existe una infinidad numerable de vectores de rotación correspondientes a una rotación cualquiera dada. Además, todas las rotaciones de valor Plantilla:Math son iguales a una rotación nula, por lo que, para un entero Plantilla:Math dado, todos los vectores de rotación de longitud Plantilla:Math, en todas las direcciones, constituyen una infinidad numerable de vectores de rotación de dos parámetros, que codifican la misma rotación que el vector cero. Estos hechos deben tenerse en cuenta al invertir la aplicación exponencial, es decir, al encontrar un vector de rotación que corresponda a una matriz de rotación dada. En consecuencia, la aplicación exponencial es sobreyectiva, no biyectiva.

Ejemplo

Un observador está situado sobre el suelo, y elige la dirección Plantilla:Math negativa para representar la dirección de la gravedad. Si gira 90° a la izquierda, rotará Plantilla:Math radianes con respecto al eje Plantilla:Math. Viendo la notación axial-angular como un par ordenado, este giro queda representado por

$(\mathrm{e}\mathrm{j}\mathrm{e},\text{ángulo})=(\left[\begin{array}{c}{e}_x\\ e_y\\ e_z\end{array}\right],\theta )=(\left[\begin{array}{c}0\\ 0\\ 1\end{array}\right],\frac{\pi }{2}).$

El ejemplo anterior se puede representar como un vector de rotación con una magnitud de Plantilla:Math apuntando en la dirección Plantilla:Math,

$\left[\begin{array}{c}0\\ 0\\ \frac{\pi }{2}\end{array}\right].$

Uso

La notación axial-angular es conveniente cuando se trata de mecánica del sólido rígido. Es útil tanto para caracterizar movimientos de rotación, como para abordar el estudio del movimiento del sólido rígido mediante transformaciones homogéneas y giros.

Cuando un cuerpo rígido gira alrededor de un eje fijo, sus datos axiales-angulares son una constante y el ángulo de rotación depende continuamente del tiempo.

La conexión entre los tres valores propios 1 y Plantilla:Math y sus tres ejes ortogonales asociados en una representación cartesiana según el teorema de Mercer, es una construcción conveniente de la expresión cartesiana de la matriz de rotación en tres dimensiones.

Rotación de un vector

La fórmula de rotación de Rodrigues, llamada así por Olinde Rodrigues, es un algoritmo eficiente para rotar un vector en un espacio euclídeo dado un eje y un ángulo de rotación. En otras palabras, proporciona un algoritmo para calcular la aplicación exponencial de Plantilla:Math a Plantilla:Math sin calcular la matriz exponencial completa.

Si Plantilla:Math es un vector en Plantilla:Math y Plantilla:Math es un vector unitario definido en el origen que describe un eje de rotación sobre el que Plantilla:Math se gira con un ángulo Plantilla:Math, la fórmula de rotación de Rodrigues para obtener el vector girado es

${𝐯}_{\mathrm{r}\mathrm{o}\mathrm{t}}=(\cos \theta )𝐯+(\mathrm{sen}\theta )(𝐞\times 𝐯)+(1-\cos \theta )(𝐞\cdot 𝐯)𝐞.$

Para la rotación de un solo vector puede ser más eficiente que convertir Plantilla:Math y Plantilla:Math en una matriz de rotación para rotar el vector.

Relación con otras notaciones

Plantilla:VT Hay varias formas de representar una rotación. Es útil comprender cómo se relacionan las diferentes notaciones entre sí, y cómo convertirlas entre ellas. Aquí el vector unitario se denota Plantilla:Math en lugar de Plantilla:Math.

Aplicación exponencial de $𝔰𝔬$ (3) a SO(3)

La aplicación exponencial efectúa una transformación de la notación axial-angular de las rotaciones a matrices de rotación,

$\exp :𝔰𝔬(3)\to \mathrm{S}\mathrm{O}(3).$

Esencialmente, al usar una serie de Taylor, se deduce una relación de forma cerrada entre estas dos notaciones. Dado un vector unidad Plantilla:Math que representa el eje de rotación unitario y un ángulo, Plantilla:Math, se da una matriz de rotación equivalente Plantilla:Math como sigue, donde Plantilla:Math es el producto vectorial de Plantilla:Math, es decir, Plantilla:Math para todos los vectores Plantilla:Math,

$R=\exp (\theta 𝐊)=\sum _{k=0}^{\mathrm{\infty }}\frac{(\theta 𝐊{)}^k}{k!}=I+\theta 𝐊+\frac{1}{2!}(\theta 𝐊{)}^2+\frac{1}{3!}(\theta 𝐊{)}^3+\cdots$

Como Plantilla:Math es asimétrico y la suma de los cuadrados de sus entradas diagonales superiores es 1, el polinomio característico de Plantilla:Math de Plantilla:Math es Plantilla:Math. Dado que, por el teorema de Cayley-Hamilton, Plantilla:Math = 0, esto implica que

${𝐊}^3=-𝐊.$

Como resultado, Plantilla:Math, Plantilla:Math, Plantilla:Math, Plantilla:Math.

Este patrón cíclico continúa indefinidamente, por lo que todas las potencias superiores de Plantilla:Math se pueden expresar en términos de Plantilla:Math y Plantilla:Math. Así, de la ecuación anterior, se deduce que

$R=I+(\theta -\frac{{\theta }^3}{3!}+\frac{{\theta }^5}{5!}-\cdots )𝐊+(\frac{{\theta }^2}{2!}-\frac{{\theta }^4}{4!}+\frac{{\theta }^6}{6!}-\cdots ){𝐊}^2,$

es decir,

$R=I+(\mathrm{sen}\theta )𝐊+(1-\cos \theta ){𝐊}^2.$

Esta es una deducción del álgebra de Lie, en contraste con el concepto geométrico que conduce a la fórmula de rotación de Rodrigues.

Debido a la existencia de la aplicación exponencial mencionada anteriormente, el vector unitario Plantilla:Math que representa el eje de rotación y el ángulo Plantilla:Math, a veces se denominan coordenadas exponenciales de la matriz de rotación Plantilla:Math.

Aplicación logarítmica de SO(3) sobre $𝔰𝔬$(3)

Sea Plantilla:Math, que denota la matriz 3 × 3 definida para efectuar el producto cruzado con el eje de rotación Plantilla:Math: Plantilla:Math para todos los vectores Plantilla:Math.

Para recuperar la notación axial-angular a partir de una matriz de rotación, se debe calcular la traza de la matriz de rotación

$\theta =\arccos \left(\frac{\mathrm{Tr}(R)-1}{2}\right)$

y luego calcular el eje normalizado,

$\mathit{\omega }=\frac{1}{2\mathrm{sen}\theta }\left[\begin{array}{c}R(3,2)-R(2,3)\\ R(1,3)-R(3,1)\\ R(2,1)-R(1,2)\end{array}\right].$

Téngase en cuenta también que el logaritmo de la matriz de rotación Plantilla:Mvar es

$\log R=\{\begin{array}{cc}0& \text{si }\theta =0\\ \frac{\theta }{2\mathrm{sen}\theta }(R-R^{𝖳})& \text{si }\theta \ne 0\text{ and }\theta \in (-\pi ,\pi )\end{array}$

Se produce una excepción cuando Plantilla:Math tiene un autovalor igual a −1. En este caso, el registro no es único. Sin embargo, incluso en el caso donde θ = Π, la norma matricial del logaritmo es

$\Vert \log (R){\Vert }_{\mathrm{F}}=\sqrt{2}|\theta |.$

Dadas las matrices de rotación Plantilla:Mvar y Plantilla:Mvar,

$d_g(A,B):={\Vert \log \left(A^{𝖳}B\right)\Vert }_{\mathrm{F}}$

es la distancia geodésica en la variedad 3D de matrices de rotación.

Para pequeñas rotaciones, el cálculo anterior de Plantilla:Math puede ser numéricamente impreciso, ya que la derivada de la función (arc cos) tiende a infinito cuando Plantilla:Math. En ese caso, los términos fuera del eje proporcionarán una mejor información sobre Plantilla:Math ya que, para ángulos pequeños, Plantilla:Math (esto se debe a que estos son los dos primeros términos de la serie de Taylor para Plantilla:Math).

Esta formulación también tiene problemas numéricos en Plantilla:Math, donde los términos fuera del eje no proporcionan información sobre el eje de rotación (que todavía se define hasta una ambigüedad de signo). En ese caso, se debe reconsiderar la fórmula anterior.

$R=I+𝐊\mathrm{sen}\theta +{𝐊}^2(1-\cos \theta )$

En Plantilla:Math, tenemos

$R=I+2{𝐊}^2=I+2(\mathit{\omega }\otimes \mathit{\omega }-I)=2\mathit{\omega }\otimes \mathit{\omega }-I$

y así se llega a

$B:=\mathit{\omega }\otimes \mathit{\omega }=\frac{1}{2}(R+I),$

y por lo tanto, los términos diagonales de Plantilla:Math son los cuadrados de los elementos de Plantilla:Math y los signos (hasta la ambigüedad de signos) se pueden determinar a partir de los signos de los términos fuera del eje de&nbspPlantilla:Math.

Cuaterniones unitarios

La siguiente expresión transforma las coordenadas de la expresión axial-angular a versores (cuaterniones):

$Q=(\cos \frac{\theta }{2},\mathit{\omega }\mathrm{sen}\frac{\theta }{2})$

Dado un versor Plantilla:Math representado mediante su escalar Plantilla:Math y el vector Plantilla:Math, las coordenadas axiales-angulares se pueden extraer utilizando el procedimiento siguiente:

$\begin{array}{cc}\theta & =2\arccos s\\ [8px]\mathit{\omega }& =\{\begin{array}{cc}\frac{x}{\mathrm{sen}\frac{\theta }{2}},& \text{si }\theta \ne 0\\ 0,& \text{en caso contrario}.\end{array}\end{array}$

Una expresión más estable numéricamente del ángulo de rotación utiliza la función arcotangente de dos parámetros:

$\theta =2\mathrm{atan2}(|x|,s),$

donde Plantilla:Math es la norma vectorial del 3-vector Plantilla:Math.

Véase también

• Coordenadas homogéneas