How to find rotation matrix from vector to another?

Question

ha ha 2017-9-20

2
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/357494-how-to-find-rotation-matrix-from-vector-to-another

编辑： ha ha 2018-11-24

I have object with 3 vector. Let's say:

e1=[a1 b1 c1];
e2=[a2 b2 c2];
e3=[a3 b3 c3];

The global coordinator system (Ox, Oy, Oz)

n1=[1 0 0];    % Ox
n2=[0 1 0];    % Oy
n3=[0 0 1];    % Oz

How to find the rotation matrix R to rotate the object to match with xyz (e1 // n1, e2 // n2, e3 // n3) ? (// :parallel)

How to find R?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Jan 2017-9-20

6
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/357494-how-to-find-rotation-matrix-from-vector-to-another#answer_282373

编辑：Jan 2017-9-20

在 MATLAB Online 中打开

Exactly as in your former question:

R = [e1; e2; e3]

is the rotation matrix already, when we assume, that these are the normalized orthogonal vectors of the local coordinate system. To convert between the two reference systems all you need is R and R.' (as long as the translation is ignored).

A vector v=[x;y;z] in the global reference system is

R * v

in the local system. Then you can convert it back to the global system by:

R.' * R * v

and because this must reply v again you get the identity:

R.' * R == eye(3)

such that

R.' == inv(R) % except for rounding errors

You can split this rotation matrix to 3 matrices around specific axes to get a set of Euler or Cardan angles.

See also https://en.wikipedia.org/wiki/Rotation_matrix.

You will find many discussions in the net, which end in flamewars, because the users cannot decide if R or R' is the actual rotation, because sometimes the rotation of the vector is meant, and sometimes the rotation of the reference system. Each field of science has its own preferences in this point. See https://www.mathworks.com/matlabcentral/answers/313150-why-is-the-result-of-quaternion-rotation-an-matrix-multiplication-not-the-same

Maybe you are looking for the "helical axes", which can define the relative attitude between two coordinate systems by a vector of translation and a rotation around this vector. If so, please clarify this.

6 个评论
显示 4更早的评论隐藏 4更早的评论

ha ha 2018-5-10

编辑：ha ha 2018-5-11

在 MATLAB Online 中打开

Corrected code:

n1=[1 0 0];n2=[0 1 0];n3=[0 0 1];e1=[1 3 0];e2=[3 -1 0];e3=[0 0 1];%input data
e1 = e1/norm(e1);e2 = e2/norm(e2);e3 = e3/norm(e3);%e1, e2, e3 must be normalized firstly
R=[e1;e2;e3];  % your general formulas R
old=[2/3;2;0];
new=R*old;

or

clear;clc;P=[2/3 2 0;1 3 0];%input data in format kx3 matrix
n1=[1 0 0];n2=[0 1 0];n3=[0 0 1];e1=[1 3 0];e2=[3 -1 0];e3=[0 0 1];%input data
e1 = e1/norm(e1);e2 = e2/norm(e2);e3 = e3/norm(e3);%e1, e2, e3 must be normalized firstly
R=[e1;e2;e3];  % rotation matrix R: 3x3 
new=(R*P')';%output data in format kx3 matrix

John Razzano 2018-5-15

Does this account for the fact that the new coordinate system is not only rotated but translated from the original? For example, if you use [0 0 0] as the "old" point you still get [0 0 0] after it has been transformed to the new frame using this method - I don't think this is correct unless I am misunderstanding something.

请先登录，再进行评论。