歐拉角奇異性產生的原因

1 歐拉角奇異性的原因.

1.1 奇異性的定義

奇異性,英文Singularity, wiki中的解釋爲

  In mathematics, a singularity is in general a point at which a given mathematical object is not defined, or a point of an exceptional set where it fails to be well-behaved in some particular way, such as differentiability.

翻譯成中文爲

  數學中,奇異性通常指的是在一個點處,某個數學對象未被定義,或者一個集合中的某個點處的某個數學特性不好,比如不可微.

1.2 歐拉角的奇異性問題

歐拉角有不同的描述方法,詳見wiki,這裏用常見的來描述飛機姿態角爲例來說明,即繞着內部座標系(intrinsic)或者機體座標系,旋轉順序爲3-2-1的姿態角描述,如圖1,其中x-y-z爲地面座標系, \psi,\theta,\psi分別爲偏航角,俯仰角,滾轉角.

圖1 歐拉角定義[2]

首先,我們來考慮這樣一個問題

給定一個姿態,是否只有一種旋轉方法從地面座標系到到給定姿態,也即是否只有一種\psi,\theta,\psi組合到給定姿態嗎?

答案是否定的,考慮一種特殊情況,當機頭向上時,即飛機呈現豎直向上飛時的姿態時,此時,\theta\pi/2, 但是偏航角和滾裝角有多種組合方式,比如,\psi=\pi/2, \phi=0\psi=\pi, \phi=-\pi/2, 如下圖

圖2 兩種不同的旋轉方法得到同樣的姿態

下面給出數學上的嚴格推導,導出歐拉角在某些特殊位置存在不唯一的解.這裏將藉助於旋轉矩陣,即

R_x(\phi)=\begin{bmatrix} 1& 0 & 0\\ 0& cos(\phi) & sin(\phi)\\ 0& -sin(\phi)& cos(\phi) \end{bmatrix}

R_y(\theta)=\begin{bmatrix} cos(\theta) & 0 & -sin(\theta)\\ 0 & 1 & 0\\ sin(\theta) & & cos(\theta) \end{bmatrix}

R_z(\psi)=\begin{bmatrix} cos(\psi) & sin(\psi) & 0\\ -sin(\psi) & cos(\psi) & 0\\ 0 & 0 & 1 \end{bmatrix}

R=R_x(\phi)\cdot R_y(\theta)\cdot R_z(\psi)=\begin{bmatrix} cos(\theta)cos(\psi) & cos(\psi)sin(\theta)sin(\phi)+sin(\psi)cos(\phi) &-cos(\psi)sin(\theta)cos(\phi)+sin(\psi)sin(\phi) \\ -cos(\theta)sin(\psi)& -sin(\psi)sin(\theta)sin(\phi)+cos(\psi)cos(\phi)&sin(\psi)sin(\theta)cos(\phi)+cos(\psi)sin(\phi) \\sin(\theta) & -sin(\phi)cos(\theta)&cos(\phi)cos(\theta) \end{bmatrix} (1)

由地面座標系到給定的飛機姿態可以得到唯一的旋轉矩陣R,詳見[3].但從R是否可以推出存在唯一的 \psi,\theta,\psi呢,如果是,則從地面座標系到給定的飛機姿態可以推出唯一的歐拉角,但是根據上面的舉例可以發現答案是否定的.那麼我們現在嘗試從R來推出\psi,\theta,\psi,記

R=\begin{bmatrix} r_{11} & r_{12} &r_{13} \\ r_{21} & r_{22} &r_{23} \\ r_{31} & r_{32} &r_{33} \end{bmatrix} (2)

則由(1) 可得

\theta=arcsin(r_{31})

\phi=arctan2(-r_{32},r_{33}) \quad (3)

\psi=arctan2(-r_{21},r_{11})

\theta=\pm \pi/2時,

R=\begin{bmatrix} 0 &sin(\psi \pm \phi) & -cos(\psi \pm \phi)\\ 0& cos(\psi \pm \phi)& sin(\psi \pm \phi)\\ \pm 1 & 0& 0\end{bmatrix}

角度\psi,\phi,不再能從(3)求出, 而從R的右上角的矩陣只可以推出\psi \pm \phi,即此時存在多個解.即從R求出唯一\psi,\theta,\psi的性質不能保證,此時把這種情況成爲奇異

參考文獻

[1] https://en.wikipedia.org/wiki/Singularity_(mathematics)

[2] https://en.wikipedia.org/wiki/Euler_angles

[3] https://en.wikipedia.org/wiki/Rotation_matrix#Rotation_matrix_from_axis_and_angle

[4] https://en.wikipedia.org/wiki/Rotation_formalisms_in_three_dimensions#Conversion_formulae_between_formalisms

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章