I want to rotate a point using Quaternion function

6 次查看（过去 30 天）
Hassan Bosha
Hassan Bosha 2020-4-6
编辑： James Tursa 2020-4-6
i want to rotate point (1,2,3) about x-axis by 45 using quaternion function
i tried to put it in identity matrix and rotate it but it's not working also i want to plot is q.plot()
here is the function of the matlab
http://www.petercorke.com/RTB/r9/html/Quaternion.html
P = [1 2 3];
R = transl(P)*eul2tr([0 0 0])
q = Quaternion(R)
R1 = R * trotx(45);
q1 = Quaternion(R1)
  2 个评论
James Tursa
James Tursa 2020-4-6
编辑：James Tursa 2020-4-6
I read through the link you provided, and unfortunately like 99% of all quaternion stuff I read online this one also lacks in specifying the quaternion convention used. There is no description of the Quaternion(R) algorithm used or the Quaternion.R method algorithm. Can you generate a random direction cosine matrix, perform both these operations, and post the results? Then maybe I can figure it out.
EDIT Nevermind ... I found a different link that has the direction cosine matrix to quaternion conversion code. I can figure it out from that.
James Tursa
James Tursa 2020-4-6
Can you give a numeric example of the inputs and desired output? I.e., give us the starting vectors and angles, and tell us what you would expect to get as output.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

James Tursa
James Tursa 2020-4-6
编辑：James Tursa 2020-4-6
I took a look at this related link:
https://github.com/petercorke/robotics-toolbox-matlab/blob/master/Octave/%40Quaternion/Quaternion.m
In there is the code for converting a direction cosine matrix to quaternion called tr2q.m. Based on that code I get the following comparison:
>> q = rand(1,4)*2-1 % start with a random quaternion
q =
0.6294 0.8116 -0.7460 0.8268
>> q = q/norm(q)
q =
0.4155 0.5357 -0.4925 0.5457
>> dcm = quat2dcm(q) % Aerospace Toolbox conversion to dcm
dcm =
-0.0807 -0.0741 0.9940
-0.9812 -0.1697 -0.0923
0.1755 -0.9827 -0.0590
>> tr2q(dcm) % Peter Corke published code conversion dcm to quaternion
ans =
0.4155 -0.5357 0.4925 -0.5457
>> dcm2quat(dcm) % Aerospace Toolbox conversion dcm to quaternion
ans =
0.4155 0.5357 -0.4925 0.5457
And a comparison with the Robotics Toolbox functions:
>> quat2dcm(q) % Aerospace Toolbox
ans =
-0.0807 -0.0741 0.9940
-0.9812 -0.1697 -0.0923
0.1755 -0.9827 -0.0590
>> quat2rotm(q) % Robotics Toolbox
ans =
-0.0807 -0.9812 0.1755
-0.0741 -0.1697 -0.9827
0.9940 -0.0923 -0.0590
>> dcm2quat(dcm) % Aerospace Toolbox
ans =
0.4155 0.5357 -0.4925 0.5457
>> rotm2quat(dcm) % Robotics Toolbox
ans =
0.4155 -0.5357 0.4925 -0.5457
So you can see that the quaternion conversions are conjugates of each other. So you need to be very careful how you use the Peter Corke code or the Robotics Toolbox code when comparing with or using MATLAB quaternion functions from the Aerospace Toolbox.
If I were to guess, maybe the Peter Corke and Robotics Toolbox code is intended for active vector rotations (rotating a vector within the same coordinate frame). But that is just a guess on my part since I have never used this code.
See this link for a discussion of the MATLAB toolbox quaternion conventions:
https://www.mathworks.com/matlabcentral/answers/465053-rotation-order-of-quatrotate

更多回答（2 个）

Matt J
Matt J 2020-4-6
q1*[1;2;3]
  4 个评论
显示 2更早的评论
Hassan Bosha
Hassan Bosha 2020-4-6
it works just fine on 2014 but not my version 2017
q = Quaternion( rpy2tr(0.1, 0.2, 0.3) )
Matt J
Matt J 2020-4-6
Well, it's 3rd party code, so you will probably have to consult with the person who wrote it.

请先登录，再进行评论。

Matt J
Matt J 2020-4-6
You could also try this FEX submission instead of the Quaternion class
https://www.mathworks.com/matlabcentral/fileexchange/30864-3d-rotation-about-shifted-axis?focused=5191309&tab=function
>> AxelRot([1;2;3], 45, [1 0 0], [])
ans =
1.0000
-0.7071
3.5355

类别

在 Help Center 和 File Exchange 中查找有关 Coordinate Transformations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by