Articulo de referencia

Fórmula de Euler-Rodrigues

En matemáticas y mecánica , la fórmula de Euler - Rodrigues describe la rotación de un vector en tres dimensiones. Se basa en la fórmula de rotación de Rodrigues , pero utiliza ...

En matemáticas y mecánica , la fórmula de Euler - Rodrigues describe la rotación de un vector en tres dimensiones. Se basa en la fórmula de rotación de Rodrigues , pero utiliza una parametrización diferente.

La rotación se describe mediante cuatro parámetros de Euler, debido a Leonhard Euler . La fórmula de rotación de Rodrigues (que lleva el nombre de Olinde Rodrigues ), un método para calcular la posición de un punto rotado, se utiliza en algunas aplicaciones de software, como simuladores de vuelo y videojuegos .

Definición

Una rotación alrededor del origen está representada por cuatro números reales, a , b , c , d , tales que

a2+b2+do2+d2=1.{\displaystyle a^{2}+b^{2}+c^{2}+d^{2}=1.}

Cuando se aplica la rotación, un punto en la posiciónincógnita{\displaystyle {\vec {x}}}gira a su nueva posición, [ 1 ]

incógnita=(a2+b2do2d22(bdoad)2(bd+ado)2(bdo+ad)a2+do2b2d22(dodab)2(bdado)2(dod+ab)a2+d2b2do2)incógnita.{\displaystyle {\vec {x}}'={\begin{pmatrix}a^{2}+b^{2}-c^{2}-d^{2}&2(bc-ad)&2(bd+ac)\\2(bc+ad)&a^{2}+c^{2}-b^{2}-d^{2}&2(cd-ab)\\2(bd-ac)&2(cd+ab)&a^{2}+d^{2}-b^{2}-c^{2}\end{pmatrix}}{\vec {x}}.}

Formulación vectorial

El parámetro a puede denominarse parámetro escalar yω=(b,do,d){\displaystyle {\vec {\omega }}=(b,c,d)}el parámetro vectorial . En notación vectorial estándar, la fórmula de rotación de Rodrigues toma la forma compacta [ 2 ] [ 3 ]

incógnita=incógnita+2a(ω×incógnita)+2(ω×(ω×incógnita)){\displaystyle {\vec {x}}'={\vec {x}}+2a({\vec {\omega }}\times {\vec {x}})+2\left({\vec {\omega }}\times ({\vec {\omega }}\times {\vec {x}})\right)}

Simetría

Los parámetros ( a , b , c , d ) y ( −a , −b , −c , −d ) describen la misma rotación. Aparte de esta simetría, cada conjunto de cuatro parámetros describe una rotación única en el espacio tridimensional.

Composición de rotaciones

La composición de dos rotaciones es en sí misma una rotación. Sean ( a₁ , b₁, c₁, d₁) y (a₂, b₂, c₂, d₂ ) los parámetros de Euler de dos rotaciones . Los parámetros para la rotación compuesta ( rotación 2 después de la rotación 1) son los siguientes:

a=a1a2b1b2do1do2d1d2;b=a1b2+b1a2do1d2+d1do2;do=a1do2+do1a2d1b2+b1d2;d=a1d2+d1a2b1do2+do1b2.{\displaystyle {\begin{aligned}a&=a_{1}a_{2}-b_{1}b_{2}-c_{1}c_{2}-d_{1}d_{2};\\b&=a_{1}b_{2}+b_{1}a_{2}-c_{1}d_{2}+d_{1}c_{2};\\c&=a_{1}c_{2}+c_{1}a_{2}-d_{1}b_{2}+b_{1}d_{2};\\d&=a_{1}d_{2}+d_{1}a_{2}-b_{1}c_{2}+c_{1}b_{2}.\end{aligned}}}

Es sencillo, aunque tedioso, comprobar que + + + = 1. (Esto es esencialmente la identidad de Euler de cuatro cuadrados ) .

Ángulo de rotación y eje de rotación

Cualquier rotación central en tres dimensiones está determinada unívocamente por su eje de rotación (representado por un vector unitario k = ( ​​k x , k y , k z ) ) y el ángulo de rotación φ . Los parámetros de Euler para esta rotación se calculan de la siguiente manera:

a=porqueφ2;b=kincógnitapecadoφ2;do=kypecadoφ2;d=kzpecadoφ2.{\displaystyle {\begin{aligned}a&=\cos {\frac {\varphi }{2}};\\b&=k_{x}\sin {\frac {\varphi }{2}};\\c&=k_{y}\sin {\frac {\varphi }{2}};\\d&=k_{z}\sin {\frac {\varphi }{2}}.\end{aligned}}}

Nótese que si φ aumenta en una rotación completa de 360 ​​grados, los argumentos del seno y el coseno solo aumentan en 180 grados. Los parámetros resultantes son opuestos a los valores originales, (−a , −b , −c , −d ) ; representan la misma rotación.

En particular, la transformación identidad (rotación nula, φ = 0 ) corresponde a valores de parámetros ( a , b , c , d ) = (±1, 0, 0, 0) . Las rotaciones de 180 grados alrededor de cualquier eje dan como resultado a = 0 .

Conexión con los cuaterniones

Los parámetros de Euler pueden verse como los coeficientes de un cuaternión ; el parámetro escalar a es la parte real, los parámetros vectoriales b , c , d son las partes imaginarias. Por lo tanto, tenemos el cuaternión

q=a+bi+doj+dk,{\displaystyle q=a+bi+cj+dk,}

que es un cuaternión de longitud unitaria (o versor ) ya que

q2=a2+b2+do2+d2=1.{\displaystyle \left\|q\right\|^{2}=a^{2}+b^{2}+c^{2}+d^{2}=1.}

Lo más importante es que las ecuaciones anteriores para la composición de rotaciones son precisamente las ecuaciones para la multiplicación de cuaterniones.q=q2q1{\displaystyle q=q_{2}\,q_{1}}En otras palabras, el grupo de cuaterniones unitarios con multiplicación, módulo el signo negativo, es isomorfo al grupo de rotaciones con composición.

Una rotación en 3D puede representarse, por lo tanto, mediante un cuaternión unitario q :

q=porqueφ2+kpecadoφ2=a+ω,{\displaystyle q=\cos {\frac {\varphi }{2}}+\mathbf {k} \sin {\frac {\varphi }{2}}=a+{\vec {\omega }},}

dónde:

  • porqueφ2{\displaystyle \cos {\tfrac {\varphi }{2}}}es la parte escalar (real),
  • kpecadoφ2{\displaystyle \mathbf {k} \sin {\tfrac {\varphi }{2}}}es la parte vectorial (imaginaria),
  • k{\displaystyle \mathbf {k} }es un vector unitario que representa el eje de rotación.

Para un vector 3Dv{\displaystyle \mathbf {v} }identificado como un cuaternión puro, el vector rotadov{\displaystyle \mathbf {v} '}está dado por:

v=qvq1,{\displaystyle \mathbf {v} '=q\mathbf {v} q^{-1},}

dóndeq1=porqueφ2kpecadoφ2=aω.{\displaystyle q^{-1}=\cos {\tfrac {\varphi }{2}}-\mathbf {k} \sin {\tfrac {\varphi }{2}}=a-{\vec {\omega }}.}

Se puede demostrar que la ecuación anterior es equivalente a la ecuación de Euler-Rodrigues aprovechando la relación de conmutación.

q2q1=q1q22k1×k2{\displaystyle q_{2}q_{1}=q_{1}q_{2}-2\mathbf {k_{1}} \times \mathbf {k_{2}} }

y el hecho de queqq1=1{\displaystyle qq^{-1}=1}como sigue:

qvq1=q(q1v2(ω)×v){\displaystyle q\mathbf {v} q^{-1}=q(q^{-1}\mathbf {v} -2(-{\vec {\omega }})\times \mathbf {v} )}
qvq1=qq1v+(a+ω)q(2ω×v){\displaystyle {\phantom {q\mathbf {v} q^{-1}}}={\cancel {qq^{-}1}}\mathbf {v} +\underbrace {(a+{\vec {\omega }})} _{q}(2{\vec {\omega }}\times \mathbf {v} )}
qvq1=v+2a(ω×v)+2(ω(ω×v)+ω×(ω×v)){\displaystyle {\phantom {q\mathbf {v} q^{-1}}}=\mathbf {v} +2a({\vec {\omega }}\times \mathbf {v} )+2({\cancel {-{\vec {\omega }}\cdot ({\vec {\omega }}\times \mathbf {v} )}}+{\vec {\omega }}\times ({\vec {\omega }}\times \mathbf {v} ))}
qvq1=v+2a(ω×v)+2(ω×(ω×v)){\displaystyle {\phantom {q\mathbf {v} q^{-1}}}=\mathbf {v} +2a({\vec {\omega }}\times \mathbf {v} )+2({\vec {\omega }}\times ({\vec {\omega }}\times \mathbf {v} ))}

donde, en la penúltima igualdad, el producto escalar de dos vectores ortogonales es cero.

Conexión con matrices de espín SU(2)

El grupo de Lie SU(2) se puede utilizar para representar rotaciones tridimensionales en matrices complejas de 2 × 2. La matriz SU(2) correspondiente a una rotación, en términos de sus parámetros de Euler, es

U=( adidobidobia+di).{\displaystyle U={\begin{pmatrix}\ a-di&-c-bi\\c-bi&a+di\end{pmatrix}}.}

que se puede escribir como la suma

U=a (1001)ib (0110)ido (0ii0)id (1001)=aIibσincógnitaidoσyidσz,{\displaystyle {\begin{aligned}U&=a\ {\begin{pmatrix}1&0\\0&1\end{pmatrix}}-ib\ {\begin{pmatrix}0&1\\1&0\end{pmatrix}}-ic\ {\begin{pmatrix}0&-i\\i&0\end{pmatrix}}-id\ {\begin{pmatrix}1&0\\0&-1\end{pmatrix}}\\&=a\,I-ib\,\sigma _{x}-ic\,\sigma _{y}-id\,\sigma _{z},\end{aligned}}}

donde las σ i son las matrices de espín de Pauli .

La rotación viene dada porincógnita(incógnita1σincógnita+incógnita2σy+incógnita3σz)=UincógnitaU=(aIibσincógnitaidoσyidσz)(incógnita1σincógnita+incógnita2σy+incógnita3σz)(aI+ibσincógnita+idoσy+idσz){\displaystyle X^{\prime }\equiv (x_{1}^{\prime }\sigma _{x}+x_{2}^{\prime }\sigma _{y}+x_{3}^{\prime }\sigma _{z})=U\;X\;U^{\dagger }=(a\,I-ib\,\sigma _{x}-ic\,\sigma _{y}-id\,\sigma _{z})(x_{1}\sigma _{x}+x_{2}\sigma _{y}+x_{3}\sigma _{z})(a\,I+ib\,\sigma _{x}+ic\,\sigma _{y}+id\,\sigma _{z})}, lo cual se puede confirmar multiplicando se obtiene la fórmula de Euler-Rodrigues como se indicó anteriormente.

Así, los parámetros de Euler son las coordenadas real e imaginaria en una matriz SU(2) correspondiente a un elemento del grupo de espín Spin(3), que se transforma mediante una doble cobertura en una rotación del grupo ortogonal SO(3). Esto realizaR3{\displaystyle \mathbb {R} ^{3}}como la única representación irreducible tridimensional del grupo de Lie SU(2) ≈ Spin(3).

Parámetros de Cayley-Klein

Los elementos de la matrizU{\displaystyle U}se conocen como parámetros de Cayley - Klein , en honor a los matemáticos Arthur Cayley y Felix Klein , [ a ]

α=adiβ=dobiγ=dobiδ= a+di{\displaystyle {\begin{aligned}\alpha &=a-di&\beta &=-c-bi\\\gamma &=c-bi&\delta &=\ a+di\end{aligned}}}

En términos de estos parámetros, la fórmula de Euler-Rodrigues también se puede escribir [ 4 ] [ 8 ] [ a ].

incógnita=(12(α2γ2+δ2β2)12i(γ2α2+δ2β2)γδαβ12i(α2+γ2β2δ2)12(α2+γ2+β2+δ2)i(αβ+γδ)βδαγi(αγ+βδ)αδ+βγ)incógnita.{\displaystyle {\vec {x}}'={\begin{pmatrix}{\frac {1}{2}}(\alpha ^{2}-\gamma ^{2}+\delta ^{2}-\beta ^{2})&{\frac {1}{2}}i(\gamma ^{2}-\alpha ^{2}+\delta ^{2}-\beta ^{2})&\gamma \delta -\alpha \beta \\{\frac {1}{2}}i(\alpha ^{2}+\gamma ^{2}-\beta ^{2}-\delta ^{2})&{\frac {1}{2}}(\alpha ^{2}+\gamma ^{2}+\beta ^{2}+\delta ^{2})&-i(\alpha \beta +\gamma \delta )\\\beta \delta -\alpha \gamma &i(\alpha \gamma +\beta \delta )&\alpha \delta +\beta \gamma \end{pmatrix}}{\vec {x}}.}

Klein y Arnold Sommerfeld utilizaron ampliamente los parámetros en relación con las transformaciones de Möbius y las relaciones cruzadas en su análisis de la dinámica de los giroscopios . [ 5 ] [ 9 ]

Véase también

Notas

  1. 1 2 Goldstein (1980) [ 4 ] considera una transformación pasiva ( contravariante o "alias"), en lugar de la transformación activa (covariante o "coartada") aquí. Su matrizA{\displaystyle A}por lo tanto corresponde a la transpuesta de la matriz de Euler-Rodrigues dada al inicio de este artículo, o, equivalentemente, a la matriz de Euler-Rodrigues para una rotación activa deφ{\displaystyle -\varphi }en vez deφ{\displaystyle \varphi }. Teniendo esto en cuenta, es evidente que sumi1{\displaystyle e_{1}},mi2{\displaystyle e_{2}}, ymi3{\displaystyle e_{3}}en la ecuación 4-67 (pág. 153) son iguales ab{\displaystyle b},do{\displaystyle c}, yd{\displaystyle d}aquí. Sin embargo suα{\displaystyle \alpha },β{\displaystyle \beta },γ{\displaystyle \gamma }, yδ{\displaystyle \delta }, los elementos de su matrizQ{\displaystyle Q}, corresponden a los elementos de la matrizU{\displaystyle U^{\dagger }}aquí, en lugar de la matrizU{\displaystyle U}. Esto da entonces su parametrización
    α=a+diβ=do+biγ=do+biδ=adi{\displaystyle {\begin{aligned}\alpha &=\;\;a+di&\beta &=c+bi\\\gamma &=-c+bi&\delta &=a-di\end{aligned}}}
    En consecuencia, mientras que su fórmula (4-64) es idéntica símbolo por símbolo a la matriz de transformación dada aquí, utilizando sus definiciones paraα{\displaystyle \alpha },β{\displaystyle \beta },γ{\displaystyle \gamma }, yδ{\displaystyle \delta }le da su matrizA{\displaystyle A}, mientras que las definiciones basadas en la matrizU{\displaystyle U}arriba conducen a la matriz de Euler-Rodrigues (activa) que se presenta aquí. Pennestrì et al (2016) [ 5 ] definen de manera similar suα{\displaystyle \alpha },β{\displaystyle \beta },γ{\displaystyle \gamma }, yδ{\displaystyle \delta }en términos de la matriz pasivaQ{\displaystyle Q}en lugar de la matriz activaU{\displaystyle U}. La parametrización aquí concuerda con la utilizada en, por ejemplo, Sakurai y Napolitano (2020), [ 6 ] p. 165, y Altmann (1986), [ 7 ] ec. 5 p. 113 / ec. 9 p. 117.

Lecturas adicionales

  • Cartan, Élie (1981). La teoría de los espinores . Dover. ISBN 0-486-64070-1.
  • Hamilton, WR (1899). Elementos de cuaterniones . Cambridge University Press.
  • Haug, EJ (1984). Análisis y optimización asistidos por computadora de la dinámica de sistemas mecánicos . Springer-Verlag.
  • Garza, Eduardo; Pacheco Quintanilla, ME (junio de 2011). "Benjamin Olinde Rodrigues, matemático y filántropo, y su influencia en la Física Mexicana" . Revista Mexicana de Física E (en español): 109– 113.
  • Shuster, Malcolm D. (1993). "Un estudio de las representaciones de actitudes" (PDF) . Journal of the Astronautical Sciences . 41 (4): 439– 517.
  • Dai, Jian S. (octubre de 2015). "Variaciones de la fórmula de Euler-Rodrigues, conjugación de cuaterniones y conexiones intrínsecas" . Mechanism and Machine Theory . 92 : 144–152 . doi : 10.1016/j.mechmachtheory.2015.03.004 .

Referencias

  1. p. ej. Felix Klein (1897), La teoría matemática de la peonza , Nueva York: Scribner. pág. 4
  2. Valdenebro, Angel G. (2016). "Visualizing rotations and compositions of rotations with the Rodrigues vector". Eur. J. Phys . 37 (6) 065001. arXiv : 1607.05999 . doi : 10.1088/0143-0807/37/6/065001 ., ecuación (2)
  3. Siminovitch, David J. (1997). "Rotaciones en RMN: Parte I. Parámetros de Euler-Rodríguez y cuaterniones". Concep. Magn. Res . 9 (3): 117– 193. doi : 10.1002/(SICI)1099-0534(1997)9:3 < 149::AID-CMR3 > 3.0.CO ; 2-Z .ecuación [37]
  4. 1 2 Goldstein, H. (1980), "Los parámetros de Cayley-Klein y cantidades relacionadas". §4-5 en Mecánica clásica , 2.ª ed. Reading, MA: Addison-Wesley. pág. 153
  5. 1 2 E. Pennestrì, PP Valentini, G. Figliolini, J. Angeles (2016), "Parámetros duales de Cayley-Klein y transformada de Möbius: Teoría y aplicaciones",Mechanism and Machine Theory 106 (enero): 50-67. doi : 10.1016/j.mechmachtheory.2016.08.008.pdf disponible a través de ResearchGate
  6. Sakurai, JJ ; Napolitano, Jim (2020). Mecánica cuántica moderna (3.ª ed.). Cambridge. ISBN  978-1-108-47322-4OCLC 1202949320 .​ {{cite book}}: CS1 mantenimiento: falta el editor de ubicación ( enlace )
  7. Altmann, S. (1986), Rotaciones, cuaterniones y grupos dobles . Oxford: Clarendon Press. ISBN 0-19-855372-2
  8. Weisstein, Eric W. , Parámetros de Cayley-Klein , MathWorld . Consultado el 10 de mayo de 2024.
  9. Felix Klein y Arnold Sommerfeld , Über die Theorie des Kreisels , vol 1. (Teubner, 1897). Traducido (2008) como: La teoría de la cima , vol 1. Boston: Birkhauser. ISBN 0817647201