La traducción de esta página aún no se ha actualizado a la versión más reciente. Haga clic aquí para ver la última versión en inglés.

rotmat

Convertir un cuaternión en matriz de rotación

Sintaxis

rotationMatrix = rotmat(quat,rotationType)

Descripción

rotationMatrix = rotmat(quat,rotationType) convierte el cuaternión, quat, en una representación de la matriz de rotación equivalente.

Ejemplos

Convertir un cuaternión en matriz de rotación para rotación de puntos

Defina un cuaternión para usarlo en la rotación de puntos.

theta = 45;
gamma = 30;
quat = quaternion([0,theta,gamma],"eulerd","ZYX","point")

quat = quaternion
       0.8924 +  0.23912i +  0.36964j + 0.099046k

Convierta el cuaternión en una matriz de rotación.

rotationMatrix = rotmat(quat,"point")

rotationMatrix = 3×3

    0.7071   -0.0000    0.7071
    0.3536    0.8660   -0.3536
   -0.6124    0.5000    0.6124

Para verificar la matriz de rotación, cree directamente dos matrices de rotación correspondientes a las rotaciones sobre los ejes x e y. Multiplique las matrices de rotación y compare con la salida de rotmat.

theta = 45;
gamma = 30;

ry = [cosd(theta)   0           sind(theta) ; ...
      0             1           0           ; ...
     -sind(theta)   0           cosd(theta)];
 
rx = [1             0           0           ;      ...
      0             cosd(gamma) -sind(gamma) ;     ...
      0             sind(gamma) cosd(gamma)];

rotationMatrixVerification = rx*ry

rotationMatrixVerification = 3×3

    0.7071         0    0.7071
    0.3536    0.8660   -0.3536
   -0.6124    0.5000    0.6124

Convertir un cuaternión en matriz de rotación para rotación de marcos

Defina un cuaternión para usar en la rotación de marcos.

theta = 45;
gamma = 30;
quat = quaternion([0,theta,gamma],"eulerd","ZYX","frame")

quat = quaternion
       0.8924 +  0.23912i +  0.36964j - 0.099046k

Convierta el cuaternión en una matriz de rotación.

rotationMatrix = rotmat(quat,"frame")

rotationMatrix = 3×3

    0.7071   -0.0000   -0.7071
    0.3536    0.8660    0.3536
    0.6124   -0.5000    0.6124

Para verificar la matriz de rotación, cree directamente dos matrices de rotación correspondientes a las rotaciones sobre los ejes x e y. Multiplique las matrices de rotación y compare con la salida de rotmat.

theta = 45;
gamma = 30;

ry = [cosd(theta)   0           -sind(theta) ; ...
      0             1           0           ; ...
     sind(theta)   0           cosd(theta)];
 
rx = [1             0           0           ;      ...
      0             cosd(gamma) sind(gamma) ;     ...
      0             -sind(gamma) cosd(gamma)];

rotationMatrixVerification = rx*ry

rotationMatrixVerification = 3×3

    0.7071         0   -0.7071
    0.3536    0.8660    0.3536
    0.6124   -0.5000    0.6124

Convertir un vector de cuaternión en matrices de rotación

Cree un vector de cuaternión normalizado de 3 por 1.

qVec = normalize(quaternion(randn(3,4)));

Convierta el arreglo de cuaterniones en matrices de rotación. Las páginas de rotmatArray corresponden al índice lineal de qVec.

rotmatArray = rotmat(qVec,"frame");

Suponga que qVec y rotmatArray corresponden a una secuencia de rotaciones. Combine las rotaciones de cuaterniones en una sola representación, luego aplique la rotación de cuaterniones a puntos cartesianos inicializados arbitrariamente.

loc = normalize(randn(1,3));
quat = prod(qVec);
rotateframe(quat,loc)

ans = 1×3

    0.9524    0.5297    0.9013

Combine las matrices de rotación en una sola representación, luego aplique la matriz de rotación a los mismos puntos cartesianos iniciales. Verifique que la rotación del cuaternión y la matriz de rotación den como resultado la misma orientación.

totalRotMat = eye(3);
for i = 1:size(rotmatArray,3)
    totalRotMat = rotmatArray(:,:,i)*totalRotMat;
end
totalRotMat*loc'

ans = 3×1

    0.9524
    0.5297
    0.9013

Argumentos de entrada

quat — Cuaternión para convertir

objeto quaternion | arreglo de objetos quaternion

Cuaternión que se va a convertir, especificado como un objeto quaternion o un arreglo de objetos quaternion de cualquier dimensionalidad.

rotationType — Tipo de rotación

"frame" | "point"

Tipo de rotación representado por la salida rotationMatrix, especificado como "frame" o "point".

Tipos de datos: char | string

Argumentos de salida

rotationMatrix — Representación de matriz de rotación

matriz numérica de 3 por 3 | arreglo numérico de 3 por 3 por N

Representación de matriz de rotación, devuelta como una matriz numérica de 3 por 3 o un arreglo numérico de 3 por 3 por N.

• Si quat es un escalar, rotationMatrix se devuelve como una matriz de 3 por 3.

• Si quat es no escalar, rotationMatrix se devuelve como un arreglo de 3 por 3 por N, donde rotationMatrix(:,:,i) es la matriz de rotación correspondiente a quat(i).

El tipo de datos de la matriz de rotación es el mismo que el tipo de datos subyacentes de quat.

Tipos de datos: single | double

Algoritmos

Dado un cuaternión con el formato

$$q=a+bi+cj+dk,$$

la matriz de rotación equivalente para rotación de marcos se define como

$$\left[\begin{array}{ccc}2a^2-1+2b^2& 2bc+2ad& 2bd-2ac\\ 2bc-2ad& 2a^2-1+2c^2& 2cd+2ab\\ 2bd+2ac& 2cd-2ab& 2a^2-1+2d^2\end{array}\right].$$

La matriz de rotación equivalente para rotación de puntos es la trasposición de la matriz de rotación de marcos:

$$\left[\begin{array}{ccc}2a^2-1+2b^2& 2bc-2ad& 2bd+2ac\\ 2bc+2ad& 2a^2-1+2c^2& 2cd-2ab\\ 2bd-2ac& 2cd+2ab& 2a^2-1+2d^2\end{array}\right].$$

Referencias

[1] Kuipers, Jack B. Quaternions and Rotation Sequences: A Primer with Applications to Orbits, Aerospace, and Virtual Reality. Princeton, NJ: Princeton University Press, 2007.

Capacidades ampliadas

Generación de código C/C++
Genere código C y C++ mediante MATLAB® Coder™.

Historial de versiones

Introducido en R2018a

Consulte también

Funciones

• rotvec | rotvecd | euler | eulerd

Objetos

• quaternion