旋轉變換的一些個人理解

對於座標系ww上的一點(x,y)(x,y)經過線性變換AA得到的點是(x,y)(x^\prime,y^\prime)。注意,這裏面的點還是在ww座標系下的。比如,旋轉變換L\bf{L},原基底是(1,0),(0,1)(1,0),(0,1),經過旋轉變換後得到的是(cosθ,sinθ),(sinθ,cosθ)(\cos \theta,\sin \theta),(-\sin \theta,\cos\theta),因此旋轉矩陣可以寫作爲:
(1)L=(cosθsinθsinθcosθ)\bf{L}= \left( \begin{matrix} \cos \theta&\sin \theta\\ -\sin \theta&\cos\theta \end{matrix} \right) \tag{1}
ww座標系下的任何點的座標經過旋轉變換的座標都可以由下面的式子求得:
(2)x=Lx \bf{x^{\prime}} = Lx \tag{2}
下面思考的問題是,一個向量在旋轉後的座標系(假設爲bb),則在此座標系下的座標應該怎麼求解呢?顯然,我們可以根據幾何關係很容易得到,我們這裏從另外一個角度出發:
對於一個向量mmww座標系下的表示爲:m=ai+bjm=a\vec{i}+b\vec{j},在bb座標系下的座標爲:m=ai+bjm=a^{'}\vec{i^{'}}+b\vec{j^{'}}那麼這兩個向量應該是相同的,表示爲矩陣則是:
(3)(ab)=P(ab) \left( \begin{matrix} a\\ b \end{matrix} \right) =\bf{P} \left( \begin{matrix} a^\prime\\ b^\prime \end{matrix} \right) \tag{3}
結合23式我們可以看出P=L1\bf{P}=L^{-1},對於一個二維的旋轉矩陣:
(4)L=(cosθsinθsinθcosθ) \bf{L}=\left( \begin{matrix} \cos\theta&-\sin\theta\\ \sin\theta&\cos\theta \end{matrix} \right) \tag{4}
那麼,我們可以求得旋轉後的座標系向量表示變換矩陣:
(5)P=L1=(cosθsinθsinθcosθ)\bf{P}=\bf{L}^{-1}= \left( \begin{matrix} \cos\theta&\sin\theta\\ -\sin\theta&\cos\theta \end{matrix} \right) \tag{5}
在詢問了數學系的科大學長後才得知,這是兩個不同數學角度的表達,前者是代數角度,後者是幾何角度。啊,自己還想了一個鐘頭。=-=。
加速度計測量到的值爲(注意是機身座標系)(Axb,Ayb,Azb)T(A_x^b,A_y^b,A_z^b)^T,由於在懸停過程中,飛行器在導航座標系的加速度值爲(0,0,g)T(0,0,g)^T,因此有:
RbnAb=(00g)R_b^nA^{b}=\left(\begin{matrix}0\\0\\g\end{matrix}\right)
其中,RbnR_b^n是機身到導航座標系的變換矩陣:
Rbn=(cθsϕsθcψcθsψsϕsψ+cϕsθsψcθsψcϕcψ+sϕsθsψcϕsθsψsϕcψsθcθsϕcϕcθ)R_b^n= \left( \begin{matrix} c\theta &s\phi s\theta c\psi-c\theta s\psi & s\phi s\psi +c\phi s\theta s\psi \\ c\theta s\psi & c\phi c\psi +s\phi s\theta s\psi & c\phi s\theta s\psi -s\phi c\psi\\ -s\theta & c\theta s\phi & c\phi c\theta \end{matrix} \right)
師兄說,手推出姿態角,比較懶,找到MMA求解旋轉問題的方法:

ComplexExpand[RotationMatrix[{{0, 0, -g}, {x, y, z}}],  TargetFunctions -> {Re, Im}]//Simplify

比較快速的求解出姿態角。

矩陣的線性變換的再次理解

x1, ,xn\boldsymbol{x_1},\cdots ,\boldsymbol{x_n}以及y1,yn\boldsymbol{y_1},\cdots \boldsymbol{y_n}是空間V\boldsymbol{V}的兩組基,那麼:
yi=k=1nakixi=(x1, ,xn)(a1iani) \boldsymbol{y_i} = \sum_{k = 1}^{n}a_{ki}\boldsymbol{x_i} = (\boldsymbol{x_1,\cdots,x_n}) \left( \begin{matrix} a_{1i}\\ \dots\\ a_{ni} \end{matrix} \right)
因此,可以用矩陣來表示:
y=xA \boldsymbol{y} = \boldsymbol{x}A
稱矩陣AA爲過渡矩陣。

如果基底x\boldsymbol{x}採用的是標準正交基,那麼y=A=(a1, ,an)\boldsymbol{y} = A = (\boldsymbol{a_1,\cdots,a_n}),即變換矩陣的每一列的就是另外一組基。

下面討論座標的變換

對於空間V\boldsymbol{V}中的任意一個向量l=i=1nξixi=j=1nηjyj\boldsymbol{l} = \sum_{i = 1}^{n}\xi_{i}\boldsymbol{x_i} = \sum_{j = 1}^{n}\eta_{j}\boldsymbol{y_j}

那麼:
l=(x1, ,xn)(ξ1ξn)=(y1, ,yn)A1(ξ1ξn)=(y1, ,yn)(η1ηn) \boldsymbol{l} = (\boldsymbol{x_1,\cdots,x_n})\left(\begin{matrix}\xi_{1}\\\cdots\\\xi_{n}\end{matrix}\right)= (\boldsymbol{y_1,\cdots,y_n})A^{-1}\left(\begin{matrix}\xi_{1}\\\cdots\\\xi_{n}\end{matrix}\right) = (\boldsymbol{y_1,\cdots,y_n})\left(\begin{matrix}\eta_{1}\\\cdots\\\eta_{n}\end{matrix}\right)
那麼:
(η1ηn)=A1(ξ1ξn) \left(\begin{matrix}\eta_1\\\cdots\\\eta_n\end{matrix}\right) = A^{-1}\left(\begin{matrix}\xi_1\\\cdots\\\xi_{n}\end{matrix}\right)
比如在飛控中,統一物理量在不同座標系下的表示要需要這個trick。

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