方位

                       方位

 表示方位的方法有三种，分别为矩阵，欧拉角，四元数表示法。

描述方位必须要有一个参照方位，如果没有参照方位，矩阵，欧拉角和四元数表示法只可以代表一种角位移。

 

方位其实就是在原方位上的旋转，那怎么样确定旋转的正负方向呢？

在左手坐标系中，从旋转轴的正方向望向原点，顺时针为正方向，逆时针为负方向；在右手坐标系中则相反，即逆时针为正方向，顺时针为负方向。

 

本文使用左手坐标系。

 一。矩阵

表示角位移的矩阵其实就是旋转矩阵，此矩阵为绕某个旋转轴旋转。

图形API（如OPENGL，DX都是用矩阵来表示方位）。

用矩阵表示角位移需要存放9个数（如果用4X4矩阵则是16个数），占用内存大。

 

矩阵表示法的优缺点详细说明见(p142)

二．欧拉角

欧拉角用三个数表示角位移，分别为heading, pitch, bank，这三个数都是角度值，heading表示绕物体坐标系的Y轴旋转的角度，picth为X轴，bank为Z轴。旋转顺序其实是无关紧要的，我们一般就按heading->pitch->bank顺序操作，本文所提到的欧拉角都是使用这种表示方法和旋转顺序。

 

由于角度描述和旋转顺序的多样性导致了表示同一方位的欧拉角表述的不唯一性，为了排除这种不唯一性，我们将heading和bank限制在+180到-180度之间，pitch限制在+90到-90度之间。

 

欧拉角有一个有趣的问题，当pitch为+90或-90度时，将导致heading,bank的旋转轴相同（注意，heading,pitch,bank表示绕物体坐标系的轴,而不是原坐标系），称做万向锁，也叫做万向节。这将导致另外一种别名问题，由于heading,bank的旋转轴相同，它们可以相互抵消或相加，而最终结果是可以只用heading或bank一个旋转就可以了。为此，我们规定，当picth为+90或-90度时，我们用heading完成竖直轴的旋转，bank定为0。万向节还导致欧拉角描述方法丢失一个自由度，这就是著名的万向节死锁问题。

 

我们设A(φ,θ,ψ) 为heading=φ, pitch=θ, bank=ψ的方位，P（φ）为绕Y轴旋转φ个角度， N（φ,θ）为绕X轴旋转θ个角度，R（φ,θ,ψ）为绕Z轴旋转ψ个角度。有如图的计算方式：

 

A(φ,θ,ψ) = R(φ,θ,ψ)N(φ,θ)P(φ)

 

解释为A的转换等于先绕Y轴旋转φ个角度，再绕X轴旋转θ个角度，最后绕Z轴旋转ψ个角度(从右到左的运算顺序)。

 

有公式为：

 

A(δφ + φ,θ,ψ) = P(δφ)A(φ,θ,ψ)

A(φ,δθ + θ,ψ) = N(φ,δθ)A(φ,θ,ψ)

A(φ,θ,δψ + ψ) = R(φ,θ,δψ)A(φ,θ,ψ)

 

解释为物体在现有角位移（设为角位移A）上进行*转换等于在A的相应转换上加上*转换，如：物体在现有角位移A上再绕物体坐标系的Y轴旋转*个角度等于在A的Y轴转换中加上*个角度。在这里有个例外，当N（绕X轴的转换）为正或负90度时，我们丢失了一个旋转自由度，故上面三个等式均不成立，这就是万向节死锁的原因：当N为正或负90度时，我们不能再用原有方位+P或N或R转换来得出现有方位（即不能用连续的欧拉角来表示空间中连续的角位移）。

 

三．四元数

一个四元数包含一个标量和一个3D向量:

             [w, v]

              [w, ( x, y, z)]

现在我们来复习一下复数，复数包含一个实部和虚部。复数（a,b）定义为复数a+bi，a称为实数，b为虚数，满足i的平方等于-1。共轭复数使虚部变负。

 

四元数扩展了复数系统，它使用三个虚部i,j,k,它们的关系如下：

 

              i^2 = j^2 = k^2 = -1;

              ij = k;  ji = -k;

              jk = I;  kj = -I;

              ki = j;  ik = -j;

 

一个四元数[w, (x, y, z)] 定义为复数w+xi+yj+zk。

 

四元数P被解释为角位移的轴——角对方式，旋转轴n和旋转角度θ并不直接存储在四元数中，但的确是在四元数中。

              P = [cos(θ/2), sin(θ/2)n]

               = [ cos(θ/2), (sin(θ/2) _x, sin(θ/2)n_y, sin(θ/2)n_z) ]

我们试着将旋转角度θ+360度的N数倍，旋转轴不变，不会改变P所代表的角位移，但却使P变成了-P。故任意角位移都有两种四元数的表示方法，它们互相为负。

       注：-p = [-cos(θ/2), -sin(θ/2)n]

                = [-cos(θ/2), (-sin(θ/2)n_x, -sin(θ/2)n_y, -sin(θ/2)n_z) ]

 

[1, 0], [-1, 0]表示单位角位移, 旋转角度为0度或者360度的整数倍时没有角位移，旋转轴是无交紧要的。

 

四元数的叉乘的效果和矩阵乘法差不多。

                            [w₁ v₁][w₂ v₂] = [w₁w₂-v_1*v₂  w₁v₂+w₂v₁+v₂ X v₁]

 

                            (w₁  x₁i  y₁j  z₁k)(w₂  x₂i  y₂j  z₂k)

                            =    w₁w₂ + w₁x₂i + w₁y₂j +w₁z₂k +

                                  x₁w₂i + x₁x₂i² + x₁y₂ij + x₁z₂ik +

                                  y₁w₂j + y₁x₂ji + y₁y₂j² + y₁z₂jk +

                                  z₁w₂k + z₁x₂ki + z₁y₂kj + z₁z₂k²

                            =    w₁w₂ – x₁x₂ – y₁y₂ – z₁z₂ +

                                   (w₁x₂ + x₁w₂ + y₁z₂ – z₁y₂)i +

                                   (w₁y₂ + x₁z₂ + y₁w₂ – z₁x₂)j +

                                   (w₁z₂ + x₁y₂ + z₁w₂ – y₁x₂)k

                            =   [ w₁w₂ + x₁x₂ + y₁y₂ - z₁z₂  

( w₁x₂ + x₁w₂ + y₁z₂ - z₁y₂

                                    w₁y₂ + x₁z₂ + y₁w₂ - z₁x₂

                                    w₁z₂ + x₁y₂ +z₁w₂ – y₁x₂) ]

四元数叉乘满足结合律不满足交换律。

我们用[0, (x, y, z)] 代表三维空间中的q点，[cos(θ/2), (sin(θ/2)n_x, sin(θ/2)n_y, sin(θ/2)n_z) ] 代表四元素p, 进行以下运算后

q₁ = pqp^-1

你会发现得到的结果是q点绕n轴旋转了θ度。我们试着进行多次旋转,先进行a旋转再进行b旋转得到

q₁  = b(aqa^-1)b^-1, 由于四元数叉乘满足结合率不满足交换率，

 = (ba)q(a^-1b^-1), 由于a^-1b^-1 = (ba)^-1

 =(ba) q (ba)^-1

              先执行a旋转再b旋转等于ba叉乘所代表的单一旋转。

 

通过四元素的叉乘我们可以将多次旋转连接起来，不过最终得到的结果连接顺序是由后到前的。为了解决这个问题，我们可以将四元素叉乘代表旋转的标准q₁ = pqp^-1改为 q₁ = p^-1qp ，得出q₁  = b^-1(a^-1qa)b = (ab)^-1q(ab)。本文将采取这种非标准的计算方法。在本文中非特别说明，没有符号的ab代表的a和b 的叉乘。

 

||P||表示四元数的模。公式：(见P145)， 表示角位移的四元数的模为1。

四元数叉乘的模等于模的乘积：

                     ||q X p|| = ||q|| * ||p||

四元数逆的叉乘等于相反顺序的叉乘的逆：

                     q^-1 X p^-1 = (p X q)^-1

四元数共轭和复数一样使虚部变负，记为：

q* = [w, -v]。

四元数的逆等于共轭除以模，记为：

q^-1 = q*/||q||。

 

四元数的差

四元数的差可以代表从一个方位到另一个方位的角位移。已知角位移a和角位移b，求它们之间的角位移t，由于方位x先旋转a再旋转t等于方位x旋转b  à at = b。我们要求出t，如果a是标量，两边同时除以a就可得出t。但我们不可以直接除以四元数，不过我们可以乘以四元数的逆以得到相同的效果。

                     a^-1at = a^-1b

                            [1,0] t = a^-1b

                                   t= a^-1b

 

四元数点乘

在本文中，点乘用*表示，叉乘用X表示

q₁*q₂ = [w₁ v₁] * [w₂ v₂]

              = w₁*w₂ + v₁*v₂

              = [w₁ (x₁ y₁ z₁) ] * [w₂ (x₂ y₂ z₂) ]

              = w₁*w₂ + x₁*x₂ + y₁*y₂ + z₁*z₂

像向量的点乘一样，四元数的点乘也是一个标量，且点乘结果的正负并不重要，因为

q₁*q₂ = - (q₁* -q₂)

而q₂和-q₂代表相同的角位移。

对于单位四元数a,b，有 -1 <= a*b <= 1。

四元数点乘的几何意义是，点乘绝对值越大，代表这两个角位移越“相似”，四元数的点乘也代表两个四元数夹角的cos值。推导过程如下：

设两个角位移a和b，要求它们的夹角的cos值。四元数(角位移)自身只代表一种转换，它们的夹角是指某个方位x经过角位移a后的方位要和方位x经过角位移b后的方位相同必须经过的角位移t的旋转角度θ。四元数的差可以求出这个t = a^-1b，我们要求出角度θ就不难了，因为t的实部代表cos(θ/2)。我们先求出t：

                                   t= a^-1b

                                    = a^* xb

                                    = [w_a -v_a][w_b v_b]

                                    = [w_aw_b–(-v_a)*v_b   w_av_b-w_bv_a+v_b x (-v_a)]

我们已经求出t的实部为w_aw_b+v_a*v_b，这和点乘的定义是一样的，所以我们说“四元数的点乘代表两个四元数夹角的cos值”。

四元数指数，对数，标量乘运算

对数：

              log q = log([cosЄ sinЄn])

                      = [ 0, Єn ]

指数：

              exp( q )

              exp([ 0, Єn ])  = [ cosЄ sinЄn ]

 

              exp( log(q) ) = q

标量乘：

              k[w, v] = [kw, kv]

四元数幂：

四元数能作为底数，记作 q^t （不要和指数运算混淆，指数运算只接受一个四元数作为参数，而四元数求幂有两个参数 ---- 四元数和指数）。四元数求幂的意义类似于实数求幂。回忆一下，a⁰ = 1， a¹ = a，a为非零标量。当t从0变到1时，a^t从1到a。四元数求幂有类似的结论：当t从0变到1， q^t从[1, 0 ]到 q 。从四元数的叉乘可以证明这一结果。

这对四元数求幂非常有用，因为它可以从角位移中抽取"一部分"。例如，四元数 q 代表一个角位移，现在想要得到代表1/3这个角位移的四元数，可以这样计算： q^1/3。

指数超出[0, 1]范围外的几何行为和预期的一样（但有一个重要的注意事项）。例如， q²代表的角位移是 q 的两倍。假设 q 代表绕x轴顺时针旋转30度，那么 q²代表绕x轴顺时针旋转60度， q^-1/3代表绕x轴逆时针旋转10度。

上面提到的注意事项是，四元数表达角位移时使用最短圆弧，不能"绕圈"。继续上面的例子， q⁴不是预期的绕x轴顺时针旋转240度，而是逆时针80度。显然，向一个方向旋转240度等价于向相反的方向旋转80度，都能得到正确的"最终结果"。但是，在此基础上的进一步运算，产生的就可能不是预期的结果了。例如，(q4)1/2不是 q 2，尽管我们感觉应该是这样。一般来说，凡是涉及到指数运算的代数公式，如(a^s)^t = a^(st)，对四元数都不适用。

现在，我们已经理解四元数求幂可以为我们做什么了。让我们看看它的数学定义，四元数求幂定义在前一节讨论的"有用"运算上，定义如公式10.18：

注意，对于标量求幂，也有类似结论：

四元数常被用于线性插值，因为四元数的线性插值是最平滑的。最著名的就是slerp插值。具体算法见(p154).