三维旋转矩阵的计算

1、三维旋转矩阵的计算旋转矩阵(Rotationmatrix)是在乘以一个向量的时候有改变向量的方向但不改变大小的效果的矩阵。旋转矩阵不包括反演，它可以把右手坐标系改变成左手坐标系或反之。所有旋转加上反演形成了正交矩阵的集合。在三维空间中，旋转变换是最基本的变换类型之一，有多种描述方式，如 Euler 角、旋转矩阵、旋转轴/旋转角度、四元数等。本文将介绍各种描述方式以及它们之间的转换。1 .旋转矩阵用一个 3 阶正交矩阵来表示旋转变换，是一种最常用的表示方法。容易证明，3 阶正交阵的自由度为 3。注意，它的行列式必须等于 1,当等于-1 的时候相当于还做了一个镜像变换。2 .Euler 角根据 E

2、uler 定理，在三维空间中，任意一种旋转变换都可以归结为若干个沿着坐标轴旋转的组合，组合的个数不超过三个并且两个相邻的旋转必须沿着不同的坐标轴。因此，可以用三个沿着坐标轴旋转的角度来表示一个变换，称为 Euler 角。旋转变换是不可交换的，根据旋转顺序的不同，有 12 种表示方式，分别为：XYZ、XZY、XYX、XZX、YXZ、YZX、YXY、YZY、ZXY、ZYX、ZXZ、ZYZ,可以自由选择其中的一种。对于同一个变换，旋转顺序不同，Euler 角也不同，在指定 Euler 角时应当首先约定旋转顺序。2.1 Euler 角转化为旋转矩阵不妨设先绕 Z 轴旋转T,再绕 Y 轴旋转制最后绕 X

3、 轴旋转%即旋转顺序为 XYZ,旋转矩阵10R(a,/3,y)=R,，Ry-&=0g%3 .旋转轴/旋转角度用旋转轴的方向向量 n 和旋转角度 8 来表示一个旋转，其中=（勺,叫巩尸II=加+彳+.=1e诙示逆时针旋转。3.1 旋转轴/旋转角度转化为旋转矩阵设 v 是任意一个向量，定义Vn=(ny)n=nnvv=v=(Znn)v%=nxv=mj,=/7Xi,=川、，=-v,人A.A.A.AJ4AV|=v-v_L=(/+w2)vcy=cos(z),3=sin(2+kZ2+C02) 2-272*11*-77i=j=k=ijk=-1,i=JK=kjJ=ki=ik,k=y=-jiq|=Jx2+

4、/+z2+2q 称为单位四元数，如果|q|=1。一个单位四元数可以表示三维旋转。用单位四元数表示旋转可以保持一个光滑移动的相机的轨迹，适合动画生成。旋转轴/旋转角度转化为单位四元数根据旋转轴 n 和旋转角度 9,得到单位四元数 qv=(WZ)=sin，wg单位四元数转化为旋转轴/旋转角度.90q=(v,co)=(sinw?cos22,6=2tan-1()(o单位四元数转化为旋转矩阵T?(Z7,0)=I+sinBn+(1-cos)/72=I+2a)v+2v2-z0XyXL22-x-zyz22yz-x-yR(q)=1-2(/+z2)2(xy+za)2(xz-yeo)2(xy-za)2(xz+ya)

5、1-2(x2+z2)2(yzxco)2(yz+XG)l-2(x2+y2)4.4四元数的性质定义四元数的逆、乘法和除法，如下所示/=（匕，？），%=（%。2）qf=（匕，一:=（一匕，1q=1102=I.匕xv2+中产2+由?匕，o)xa2-%.匕)q=qjq、=q1q?=(%x匕+8一，一-%R(q)=&(%)&(%),RW)=&(%)及(1)根据该性质， 我们可以对两个旋转变换 q1 和 q2 作线性插值， 这相当于在四维空间中的超球面上对点 q1 和 q2 作球面线性插值。q(D=sin(l-t)B+sin。sintOsin。Q=cos-1(必=cosT(%+)也可