We can use the 3x3 rotation matrix R to determine the effect of rotations about the principle axes:
Vafter = R.Vbefore
That's straightforward.
Does anyone know what the strategy should be for solving for unknown rotations r={rx,ry,rx} radians, given we know Vbefore and Vafter?
Then presumably you know that a rotation about the x-axis though angle rx is given by the matrix with rows [1 0 0], [0 cos(rx) -sin(rx)], [0 sin(rx) cos(rx)], that a rotation about the y-axis through angle ry is given by the matrix [cos(ry) 0 -sin(ry)], [0 1 0], [sin(ry), 0 cos(ry)], and that the rotation around the z-axis through an angle rz is given by the matrix [cos(rz) -sin(rz) 0], [sin(rz) cos(rz) 0], [0 0 1].Quote:
Originally Posted by BPW
And the rotation matrix R is the product of those matrics.
Go ahead and multiply the matrices above by the Vbefore vector, and set the result equal to the Vafter vector. That gives you three equations to solve for vx, vy, and vz. (I presume you meant {rx,ry,rz}, not {rx,ry,rx}.)Quote:
Does anyone know what the strategy should be for solving for unknown rotations r={rx,ry,rx} radians, given we know Vbefore and Vafter?
Hi Halls.....
yes I realize Rzyx = Rx.Ry.Rz and that I can multiply out the equations. You get three simultaneous equations with anglex, angley and anglez as the unknowns - the problem is these equations involve Sin and Cos terms and are rather tricky to solve. For example, and there are multiple solutions. I've also read there are problems with 'ambiguities' - by that, I think they mean the solutions are:
angle = C x Pi + function
where C is a member of the class of integers. The solutions seem to have branching. I do wonder if I have stumbled over the rationale for using quaternions!
Quote:
Originally Posted by HallsofIvy
The vectors $\displaystyle \mathbf{v}_{\text{after}}$ and $\displaystyle \mathbf{v}_{\text{before}}$ obviously have to have the same length, or you can't rotate one to coincide with the other. Let's assume that they are both unit vectors.
I'm going to start with one unit vector $\displaystyle \mathbf{v} = (a,b,c)$, with $\displaystyle a^2+b^2+c^2 = 1$. Let $\displaystyle \theta =\arccos{a}$ and let $\displaystyle d = b/\sin\theta,\ e = c/\sin\theta$. (That can't be done if $\displaystyle \sin\theta=0$, but we'll come back to that later.) Then $\displaystyle \mathbf{v} = (\cos\theta,d\sin\theta,e\sin\theta)$, and $\displaystyle d^2+e^2=1$. So we can find $\displaystyle \phi$ such that $\displaystyle \cos\phi =d,\ \sin\phi = e$, and $\displaystyle \mathbf{v} = (\cos\theta,\sin\theta\cos\phi,\sin\theta\sin\phi)$.
Let $\displaystyle \mathbf{u} = (1,0,0)$. Then it is easy to check that $\displaystyle \mathbf{v} = X_\phi Z_\theta\mathbf{u}$, where $\displaystyle X_\phi$ is a rotation of $\displaystyle \phi$ radians about the x-axis and $\displaystyle Z_\theta$ is a rotation of $\displaystyle \theta$ radians about the z-axis. The inverse transformation is given by $\displaystyle \mathbf{u} = Z_{-\theta}X_{-\phi}\mathbf{v}$.
Now if you are given two vectors $\displaystyle \mathbf{v_1}$ and $\displaystyle \mathbf{v_2}$, you can rotate $\displaystyle \mathbf{v_1}$ to $\displaystyle \mathbf{v_2}$ by rotating it first to $\displaystyle \mathbf{u}$ and then on to $\displaystyle \mathbf{v_2}$.
Thus $\displaystyle \mathbf{v_2} = X_{\phi_2}Z_{\theta_2-\theta_1}X_{-\phi_1}\mathbf{v_1}$. That is more or less what you wanted. I seem to have ended with rotations {rx,rz,rx} rather than the {rx,ry,rx} that you wanted, so you'll have to switch the axes round to get that.
Finally, what happens if $\displaystyle \sin\theta=0$? In that case, $\displaystyle \cos\theta = \pm1$ and the vector $\displaystyle \mathbf{v}$ will be $\displaystyle (\pm1,0,0)$. If the sign is +, then $\displaystyle \mathbf{v}$ is already equal to $\displaystyle \mathbf{u}$. If it is -, then the rotation $\displaystyle Z_\pi$ will swing it round to $\displaystyle \mathbf{u}$.
