Combinar vectores de rotación de ejes

Question

Estoy experimentando con el uso de vectores de ejes y rotaciones en mi motor de juegos de hobby. Este es un vector de 3 componentes a lo largo del eje de rotación con una longitud de la rotación en radianes. Me gustan porque:

• A diferencia de las matrices quats o de rotación, de hecho me puedo ver los números y visualizar la rotación en mi mente
• Son un poco menos memoria que los cuaterniones o matrices.
• puedo representar valores fuera del rango de -Pi a Pi (Esto es importante si almaceno una velocidad angular)

Sin embargo, tengo un bucle estrecho que actualiza la rotación de todos mis objetos (decenas de miles) en función de su velocidad angular. Actualmente, la única forma que conozco de combinar dos vectores de ejes de rotación es convertirlos en cuaterniones, multiplicarlos y luego convertir el resultado a un eje/ángulo. A través del perfil, he identificado esto como un cuello de botella. ¿Alguien sabe un enfoque más directo?

Answer 1

Su representación es equivalente a quaternion rotation, siempre que sus vectores de rotación sean unidad de longitud. Si no desea utilizar alguna estructura de datos cuaternarios enlatados, simplemente debe asegurarse de que los vectores de rotación sean de longitud unitaria, y luego calcular el equivalente quaternion multiplications/reciprocal computation para determinar la rotación agregada. Es posible que pueda reducir el número de multiplicaciones o adiciones.

Si su ángulo es lo único que está cambiando (es decir, el eje de rotación es constante), entonces puede simplemente usar una escala lineal del ángulo y, si lo desea, modificarlo para que esté en el rango [0, 2 π). Por lo tanto, si usted tiene una velocidad de rotación de α raidans por segundo, a partir de un ángulo inicial de θ en el tiempo t , entonces el ángulo de rotación final en el tiempo t está dada por:

θ (t) = 0 θ + α (tt 0 ) mod 2 π

a continuación, sólo solicita que la rotación a su colección de vectores.

Si nada de esto mejora su rendimiento, debe considerar el uso de una biblioteca de cuaterniones enlatados, ya que las cosas ya están optimizadas para los tipos de aplicaciones que está ejecutando.

Answer 2

Debe usar unidades cuaternarias en lugar de vectores escalados para representar sus rotaciones. Se puede demostrar (no yo) que cualquier representación de rotaciones que utilice tres parámetros tendrá problemas (es decir, es singular) en algún momento. En su caso, ocurre cuando su vector tiene una longitud de 0 (es decir, la identidad) y en longitudes de 2pi, 4pi, etc. En estos casos, la representación se vuelve singular. Los cuaterniones unitarios y las matrices de rotación no tienen este problema.

Según su descripción, parece que está actualizando su estado de rotación como resultado de la integración numérica. En este caso, puede actualizar su estado de rotación convirtiendo su tasa de rotación (\ omega) en una tasa de cuaternión (q_dot).Si representamos el cuaternión como q = [q0 T1 T2 T3], donde q0 es la parte escalar a continuación:

q_dot = E*\omega 

donde

[ -q1 -q2 -q3 ] 
E = [ q0 -q3 q2 ] 
    [ q3 q0 -q1 ] 
    [ -q2 q1 q0 ] 

Luego, su actualización se convierte en

q (k + 1) = q (k) + q_dot * dt

para una integración simple. Puede elegir un integrador diferente si lo desea.

Answer 3

Puede mantenerlos como valores de eje de ángulo.

Construya una matriz de producto cruzado (anti-symmetric) utilizando los valores del eje angular (x,y,z) y pondere los elementos de esta matriz multiplicándolos por el valor del ángulo. Ahora resuma todas estas matrices de productos cruzados (one for each angle axis value) y encuentre la matriz de rotación final usando la matriz exponencial.

Si la matriz A representa esta matriz de productos cruzados (construido a partir de ángulo valor Axis) entonces,

exp(A) es equivalente a la matriz de rotación R(i.e., equivalent to your quaternion in matrix form).

Por lo tanto,

exp (A1 + A2) = R1 * R2 

probablemente un calucation más caro al final ...