Efficient Vector Rotation Matrix Routine
Suppose you have a unit vector described by f = [ai bj ck] and you would like to rotate f such that its result is t = [di ej fk]. This routine will find R(f,t) such that R(f,t)*f’ = t’. Where R(f,t) is known as the 3 x 3 transformation matrix needed to rotate f into t.
This MATLAB routine was based on a published article titled “Efficiently Building a Matrix to Rotate One Vector to Another” written by Tomas Moller and John Hughes in 1999.
This method features no square roots or trigonometric function calls and is reported to be faster than any other vector rotation matrix method tested by Moller and Hughes. In fact, the Goldman method (fastest method tested) was 50% slower than this routine in conversion speed tests.
Working Example:
Unit Vector f:
>> f = rand(1,3);
>> f = f./norm(f);
Unit Vector t:
>> t = rand(1,3);
>> t = t./norm(t);
The transformation matrix R(f,t) is found:
>> R = vecRotMat(f,t);
Now, test that R*f’= t’
>> R*f’
Citar como
Darin Koblick (2024). Efficient Vector Rotation Matrix Routine (https://www.mathworks.com/matlabcentral/fileexchange/37597-efficient-vector-rotation-matrix-routine), MATLAB Central File Exchange. Recuperado .
Compatibilidad con la versión de MATLAB
Compatibilidad con las plataformas
Windows macOS LinuxCategorías
Etiquetas
Agradecimientos
Inspiración para: rotVecAroundArbAxis(unitVec2Rotate,rotationAxisUnitVec,theta)
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!Descubra Live Editor
Cree scripts con código, salida y texto formateado en un documento ejecutable.