Cuaterniones y rotación en el espacio

De Wikipedia, la enciclopedia libre

[editar] I. Introducción

Se recuerda la versión geométrica del producto de dos cuaterniones, q = (a, u) y q´= (a', u'), donde a y a´ son las partes reales, u y u´ son las partes imaginarias, también vistas como vectores del espacio tridimensional R ³  : q.q' = (aa' - u·u', au' + a'u + u×u'). u·u´ designa el producto escalar, y u×u´ el producto vectorial. Notaremos q^~ el cuaternión conjugado de q: q^~ = (a, -u).

Para permanecer en el espacio tridimensional, hace falta hacer desaparecer las partes reales. Tomemos a = a´= 0.
Entonces q . q^~ = (0, u)(0, u') = ( - u·u' , u×u´ ).
Bien es sabido que el producto vectorial está relacionado con la rotación en el espacio. Por lo tanto, a base de productos, debe ser posible expresar cualquier rotation tridimensional. El objetivo es obtener una fórmula parecida a la expresión compleja de la rotación en el plano:

z´= wz, con w = e ^iθ cuando se gira alrededor del origen, y z´ = w(z-c) + c si se rota alrededor del punto c.

[editar] II. Descubriendo la fórmula

Tomemos el ejemplo más sencillo: ¿Cómo expresar analíticamente la rotación alrededor de eje de los x,(O,i) con un ángulo de 90 grados?

El vector i debe tener un papel. Miremos a la multiplicación por i por la izquierda: ij = k (como i×j =k) e ik = -j (como i×k =-j).

Parece por lo tanto que la función f: q → iq es la respuesta a la pregunta. En el plano (j, k) esa función rota de 90 grados. ¿Pero qué pasa en el resto del espacio? f(i) = ii = -1 y f(1) = i. Por linealidad, nos damos cuenta que hace girar el plano (1, i) de 90 grados también, y esto ¡no lo queremos! El punto i tiene que permanecer inmóvil, y la función f no tiene que enviar ni un punto del espacio usual iR + jR + kR en la cuarta dimensión (aquí, en 1R).

Como sabemos que la multiplicación no es conmutativa en H, cuerpo de los cuaterniones, miremos al producto por la derecha, por i:

g : q → qi. q(j)= ji = - k; q(k) = ki = j lo que corresponde a la rotación inversa en el plano (j, k). Pero q(1) = i e q(i) = - 1 da la misma rotación parásita que f sobre (i,j). Si tomamos la función opuesta (e inversa) : h = -g : q → -qi = q(-i) nos damos cuenta que h tiene la misma acción sobre el plano (j, k) que f pero la acción opuesta sobre (1, i). Entonces f y h se compensan en (1, i), pero se acumulan en (j,k), y la función compuesta: m = foh (=hof): q → iq(-i) deja el plano (1, i) quieto, pero gira el plano (j, k) dos veces de 90 grados, o sea de 180 grados: m(j)=ij(-i)= -iji = - ki = -j y m(k)= - k.

Hemos obtenido por lo tanto una rotación alrededor de eje (O,i), pero con un ángulo doble de lo deseado. Basta con dividir los ángulos de f e h por dos para obtener la fórmula.

i= e ^iπ/2 = cos(π/2) + i.sin(π/2) = cos 90º + i·sin 90º. El número que corresponde al medio ángulo es

e^i.π/4 = (1+i)/√2 = w,

y la función que realiza la rotación pedida es r: q → wqw^~.

Este raciocinio se generaliza a cualquier eje de rotación, y no soló a los tres ejes (O,i) (O, j) y (O, k). Si se quiere girar alrededor del eje (O,u) donde u es un vector unitario, hay que considerar el plano (1, u) y otro plano perpendicular (ortogonal) en H, y emplear el número cos(θ/2) + u·sin(θ/2).

[editar] III. La fórmula

Sea q = xi + yj + zk un punto (o un vector) del espacio, u un vector unitario del mismo espacio y θ un real. La rotación alrededor del eje (0,u) de un ángulo θ envía el punto q sobre el punto q' = a'i + b'j + c'k dado por la fórmula:

Para obtener la rotación alrededor de un eje (c,u), donde c es un punto cualquiera del espacio, basta con componer la función anterior por dos translaciones:

Note que h es un cuaternión unitario, como en el caso de los complejos.

La fórmula resulta algo más complicada que en el plano complejo porque trabajamos en cuatro dimensiones con los cuaterniones pero queremos permanecer en el espacio usual de tres dimensiones. Una simple multiplicación, a la izquierda o a la derecha, daría dos rotaciones simultáneas en dos planos perpendiculares (ortogonales) en el espacio cuadridimensional.

[editar] IV. Un ejemplo

Consideremos la rotación alrededor del eje (O, i + j + k), con un ángulo de 120º o sea 2π/3 radianes. Nos proponemos calcular la imagen del vector j. Puesto que el vector i + j + k no es de unitario, lo dividemos por su norma, con lo que obtenemos el siguiente vector unitario: u = ( i + j + k)/√3.

El medio ángulo es π/3, por lo tanto h = cos(π/3) + u.sin(π/3) = (1 + i + j + k)/2, y su conjugado es h^~ = (1 - i - j - k)/2.

El vector (o punto correspondiente) j será enviado en hjh^~. hjh^~ =(1 + i + j + k)j(1 - i - j - k)/4 = (j + k - 1 - i)(1 - i - j - k)/4 y, distribuyendo los factores, hallamos:

hjh^~ = ( (j+k+1-i) + (k-j+i+1) + (-1+i+j+k) + (-i-1+k-j) )/4 = 4k/4 = k.

Del mismo modo hallaríamos que hkh^~ = i e hih^~ = j, lo que da la expresión analítica de la rotación:

r(xi + yj + zk ) = zi + xj + yk.

Cuaterniones Número complejo