剛體運動學(xué)(6):有限與無限微小轉(zhuǎn)動

\mathrm{\mathbf{I.}}有限轉(zhuǎn)動

為了尋找一種用轉(zhuǎn)動參量(轉(zhuǎn)動角度付呕、方向余弦)表示的坐標(biāo)變換表征延曙,讓我們先考慮順時針方向的主動有限旋轉(zhuǎn)芹务。

如上圖所示蝉绷,矢量\mathbf{r}順時針經(jīng)過有限角度\Phi后變成了矢量\mathbf{r}^{\prime}

\mathbf{n}定義為沿轉(zhuǎn)軸\vec{ON}方向的單位矢量枣抱,不難得出

\vec{ON} = (\rm{proj}_{\mathbf{n}}\mathbf{r})\mathbf{n} = (\mathbf{r} \boldsymbol{\cdot} \mathbf{n})\mathbf{n}\\

根據(jù)一些簡單的幾何關(guān)系熔吗,其它的矢量也可使用\mathbf{r}\mathbf{n}一并表示:

\vec{NP} = \vec{OP} - \vec{ON} = \mathbf{r} - (\mathbf{r} \boldsymbol{\cdot} \mathbf{n})\mathbf{n}\\

\vec{NV} = ||\vec{NQ}||\cos\Phi \;\hat{NP} = [ \mathbf{r} - (\mathbf{r} \boldsymbol{\cdot} \mathbf{n})\mathbf{n}]\cos\Phi \\

\vec{VQ} = ||\vec{NQ}||\sin\Phi\; \hat{(\mathbf{r} \times \mathbf{n})} = (\mathbf{r} \times \mathbf{n})\sin\Phi\\

于是,轉(zhuǎn)動后的矢量\mathbf{r}^{\prime}可表示為

\begin{align*}\mathbf{r}^{\prime} = \vec{ON} + \vec{NV} + \vec{VQ}\\\end{align*}\\

代入相關(guān)量后可以得到

\boxed{\mathbf{r}^{\prime} = \mathbf{n}(\mathbf{n} \boldsymbol{\cdot} \mathbf{r})(1 - \cos\Phi) + \mathbf{r}\cos\Phi + (\mathbf{r} \times \mathbf{n}) \sin\Phi}\\

這是有限角度的矢量轉(zhuǎn)動公式(rotation formula)佳晶,對任何旋轉(zhuǎn)皆有效桅狠,無論角度多大。

\bullet轉(zhuǎn)動角度\Phi與歐拉角之間的關(guān)系可以通過考慮轉(zhuǎn)動算子的跡在變換前后均不變得到:

 \begin{align*}1 + 2\cos\Phi + 1 &= 1 + \cos\theta + \cos\psi\cos\phi - \sin\psi\sin\phi \\ & \quad + \cos\theta\cos\Phi\cos\psi - \cos\theta\sin\phi\sin\psi\end{align*}\\

\implies \cos^2\frac{\Phi}{2} = \cos^2\frac{\theta}{2}\cos^2\frac{\phi + \psi}{2}\\

由于當(dāng)\Phi \rightarrow 0\theta,\phi,\psi \rightarrow 0中跌,去平方后等式右邊為正咨堤,所以

\boxed{\cos\frac{\Phi}{2} = \cos\frac{\theta}{2}\cos\frac{\phi + \psi}{2}}\\

\rm{\mathbf{I\!I.}}無限微小轉(zhuǎn)動

因為存在多次不同轉(zhuǎn)動的變換矩陣是由矩陣乘法連接,矩陣乘法運算不具有對易性漩符,因此一喘,除了此時的轉(zhuǎn)動為無窮小的情況,轉(zhuǎn)動變換是通常無法使用矢量來表征的嗜暴。

考慮被動轉(zhuǎn)動變換(順時針)凸克,轉(zhuǎn)動角度取無窮小,變換前后矢量的坐標(biāo)的變化僅保留到一階小量:

x_1^{\prime} = x_1 + \epsilon_{11}x_1 + \epsilon_{12}x_2 + \epsilon_{13}x_3\\...

抽象指標(biāo)形式

x_i^{\prime} = x_i + \epsilon_{ij}x_j = (\delta_{ij} + \epsilon_{ij})x_j\\

矩陣形式

\mathbf{x}^{\prime} = (\rm{I} + \boldsymbol{\epsilon})\mathbf{x}\\

其中\boldsymbol{\epsilon}是無窮小算子

\bullet微小轉(zhuǎn)動的情況不受矩陣乘法運算不對易問題的影響闷沥,因為矩陣加法具有對易性萎战,所以現(xiàn)在矩陣\rm{I} + \boldsymbol{\epsilon}對易:

(\rm{I} + \boldsymbol{\epsilon}_1)(\rm{I} + \boldsymbol{\epsilon}_2) = (\rm{I} + \boldsymbol{\epsilon}_2)(\rm{I} + \boldsymbol{\epsilon}_1) = \rm{I} + \boldsymbol{\epsilon}_2 + \boldsymbol{\epsilon}_1 =  \rm{I} + \boldsymbol{\epsilon}_2 + \boldsymbol{\epsilon}_1\\

所以微小轉(zhuǎn)動算子可以用矢量來表征。

(例)使用泰勒展開并只保留一階項舆逃,可以得到含歐拉角的無限微小轉(zhuǎn)動矩陣

\rm{R} = \begin{bmatrix}1 & (d\phi + d\psi) & 0\\-(d\phi + d\psi) & 1 & d\theta\\0 & -d\theta &1\end{bmatrix} \\

根據(jù)剛體運動學(xué)(5)撞鹉,將\lambda = 1代入本征方程,可以得到關(guān)于轉(zhuǎn)軸的方向余弦

\begin{align*}(\rm{R} - \rm{I})d\mathbf{\Omega} &= \mathbf{0}\\\boldsymbol{\epsilon} d\mathbf{\Omega} &= \mathbf{0}\end{align*}\\

 \begin{bmatrix}1 & (d\phi + d\psi) & 0\\-(d\phi + d\psi) & 1 & d\theta\\0 & -d\theta &1\end{bmatrix} \begin{pmatrix}d\Omega_1\\ d\Omega_2\\ d\Omega_3\end{pmatrix} =  \begin{pmatrix}0\\0\\0\end{pmatrix}\\

d\Omega_3 = d\phi + d\psi颖侄,解得

 d\mathbf{\Omega} = \begin{pmatrix}d\theta\\ 0\\ d\phi + d\psi \end{pmatrix} \\

\implies d\mathbf{\Omega} = \mathbf{i}d\theta + \mathbf{k}(d\phi + d\psi)\\

\bullet矩陣\rm{R} = \rm{I} + \boldsymbol{\epsilon}的逆矩陣

\rm{R}^{-1}  = \rm{I} - \boldsymbol{\epsilon}\\

證明很簡單:

\begin{align*}\rm{R}\rm{R}^{-1} &= (\rm{I} + \boldsymbol{\epsilon})(\rm{I} - \boldsymbol{\epsilon})\\&= \rm{I}^2 - I\boldsymbol{\epsilon} + \boldsymbol{\epsilon}\rm{I}\\&= \rm{I}  \end{align*}\\

根據(jù)\rm{R}的正交性鸟雏,

\rm{R}^t = \rm{R}^{-1}\\

\rm{R}^t = (\rm{I} + \boldsymbol{\epsilon})^t = \rm{I} + \epsilon^t = \rm{R}^{-1} = \rm{I} - \boldsymbol{\epsilon}\\

\implies \boldsymbol{\epsilon}^t = -\boldsymbol{\epsilon}\\

可見,無窮小算子\boldsymbol{\epsilon}是一個反對稱矩陣

\epsilon_{ij} = -\epsilon_{ji}\\

i = j览祖,反對稱矩陣的對角元必須為零孝鹊。

i \neq j,總共可以有\frac{3 \cdot 2}{2} = 3個不同的矩陣元展蒂,不妨將無窮小算子記為如下形式:

\boldsymbol{\epsilon} = \begin{bmatrix}0 & d\Omega_3 & -d\Omega_2\\ -d\Omega_3 & 0 & d\Omega_1\\ d\Omega_2 & -d\Omega_1 & 0\end{bmatrix}\\

\bullet物理量d\Omega_1又活,d\Omega_2d\Omega_3組成矢量

d\boldsymbol{\Omega} = \begin{pmatrix}d\Omega_1\\ d\Omega_2\\ d\Omega_3 \end{pmatrix}\\

對于位矢的無限微小轉(zhuǎn)動

\mathbf{r}^{\prime} = (\rm{I} + \boldsymbol{\epsilon})\mathbf{r}\\

定義位矢在參考系的微小轉(zhuǎn)動下坐標(biāo)的變化 \mathbf{r}^{\prime} - \mathbf{r} \equiv d\mathbf{r}^{\prime}

d\mathbf{r}^{\prime} = \boldsymbol{\epsilon}\mathbf{r}\\

\begin{align*}dx_1 &= x_2d\Omega_3 - x_3d\Omega_2\\dx_2 &= x_3d\Omega_1 - x_1d\Omega_3\\dx_3 &= x_1d\Omega_2 - x_2d\Omega_1 \end{align*}\\

可見锰悼,等式右側(cè)可表示為兩個矢量的叉乘

\boxed{d\mathbf{r} = \mathbf{r} \times d\boldsymbol{\Omega}}\\

\bullet位矢\mathbf{r}的正交變換為

x_{i}^{\prime} = a_{ij} x_j\\

矢量d\boldsymbol{\Omega}的正交變換與位矢\mathbf{r}的變換在形式上相似柳骄,

d\Omega_i^{\prime} = |\rm{B}|b_{ij}d\Omega_j\\

其中|\rm{B}|為變換矩陣\rm{B}的行列式。

\bullet矢量d\mathbf{\Omega}的特征可從有限轉(zhuǎn)動公式得到箕般。

\mathbf{r}^{\prime} = \mathbf{n}(\mathbf{n} \boldsymbol{\cdot} \mathbf{r})(1 - \cos\Phi) + \mathbf{r}\cos\Phi + (\mathbf{r} \times \mathbf{n}) \sin\Phi \\

在無限微小極限下耐薯,角度\Phi \rightarrow d\Phi\cos\Phi \rightarrow 1丝里,\sin\Phi \rightarrow d\Phi曲初,轉(zhuǎn)動公式變?yōu)?/p>

\mathbf{r}^{\prime} = \mathbf{r} + (\mathbf{r} \times \mathbf{n})d\Phi\\

于是

d\mathbf{r} = \mathbf{r}^{\prime} - \mathbf{r} = \mathbf{r} \times \mathbf{n} \;d\Phi\\

\implies \boxed{d\mathbf{\Omega} = \mathbf{n}\;d\Phi}\\

可見,矢量d\mathbf{\Omega}總是沿轉(zhuǎn)軸方向杯聚。

\bullet三維坐標(biāo)反演矩陣通常具有形式S_{ij} = -\delta_{ij} \quad (i,j = 1,2,3)

\rm{S} = \begin{bmatrix}-1 & 0 & 0\\0 & -1 & 0\\0& 0 & -1\end{bmatrix}\\

根據(jù)矢量在坐標(biāo)反演下表現(xiàn)出的不同性質(zhì)臼婆,可將矢量分為極矢量(polar vector)軸矢量/偽矢量(axial vector/pseudo-vector)

(1)極矢量是變換形如

d\mathbf{r} = \mathbf{r} \times d\boldsymbol{\Omega}\\

這一類的矢量。大部分矢量幌绍,如位矢及其關(guān)于時間的微商等均屬于極矢量颁褂。極矢量在坐標(biāo)反演下會改變符號故响。

設(shè)宇稱算子為\rm{P},它執(zhí)行坐標(biāo)反演

x \rightarrow -x,\quad y \rightarrow -y,\quad z \rightarrow -z\\

那么對于極矢量\mathbf{V}颁独,則有

\rm{P}\mathbf{V} = -\mathbf{V}\\

(2)最簡單的一類偽矢量由兩個極矢量的叉乘產(chǎn)生:

\mathbf{V}^{\ast} = \mathbf{D} \times \mathbf{F}\\

\mathbf{V}_i^{\ast} = D_jF_k - F_jD_k\\

可見彩届,偽矢量在反演下不改變符號

\rm{P}\mathbf{V}^{\ast} = \mathbf{V}^{\ast}\\

\bullet偽矢量的變換遵循

v_i^{\prime} = |\rm{B}|b_{ij}v_j\\

根據(jù)剛體運動學(xué)(2),正交矩陣的行列式|\rm{B}| = \pm 1奖唯。

如果是常規(guī)正交變換惨缆,|\rm{B}| = 1

v_i^{\prime} = b_{ij}v_j\\

那么極矢量與偽矢量的變換具有相同形式丰捷。

剛才提到過坯墨,矢量d\mathbf{\Omega}遵循上述變換規(guī)律,所以是一個偽矢量病往。

\bullet對于兩種矢量的點積捣染,有

\rm{P}(\mathbf{V} \boldsymbol{\cdot} \mathbf{V}^{\ast}) = -(\mathbf{V} \boldsymbol{\cdot} \mathbf{V}^{\ast})\\

所以對于標(biāo)量,我們同樣可以沿用定義停巷,將變換方式形如

\rm{P}S^{\ast} = -S^{\ast} \\

的標(biāo)量稱為偽標(biāo)量(pseudo-vector)

一個偽標(biāo)量與一般標(biāo)量的乘積仍是一個偽標(biāo)量

\rm{P}(SS^{\ast}) = -SS^{\ast}\\

一個偽標(biāo)量與偽矢量的乘積是一個極矢量

\rm{P}(S^{\ast}\mathbf{V}^{\ast}) = -(S^{\ast}\mathbf{V}^{\ast})\\

\bullet從下列關(guān)系也可以證明d\mathbf{\Omega}是一個偽矢量:

由于\mathbf{r},d\mathbf{r}均為極矢耍攘,

d\mathbf{r} = \mathbf{r} \times d\boldsymbol{\Omega}\\

\implies d\mathbf{\Omega} = d\mathbf{r} \times \mathbf{r}\\

符合偽矢量的定義。

\mathrm{\mathbf{I\!I\!I.}}記號

\bullet由于許多表達(dá)式涉及到兩個矢量間的叉乘積畔勤,為了方便表示蕾各,引入置換符號(permutation symbol)/列維齊維塔密度(Levi-Civita density)

e_{ijk} = \begin{cases}+1 & \text{偶次置換};\\-1 & \text{奇次置換}; \\0 & \text{其它}\end{cases}\\

在張量分析中,由于置換符號并非嚴(yán)格遵循張量變換(這里主要指正交變換)但具有相似形式庆揪,人們將其稱為相對張量(relative tensor)式曲,比如矢量d\mathbf{\Omega}就是一個一階相對張量。

著名的列維齊維塔符號(Levi-Civita symbol)是一個絕對張量(absolute tensor)缸榛,二者之間的關(guān)系為

\varepsilon_{ijk} = \sqrt{|\mathbf{g}|}\;e_{ijk}\\

其中 \sqrt{|\mathbf{g}|}體積元(volume element)吝羞。

使用置換符號,兩個矢量\mathbf{A},\mathbf{B}的叉乘可以表示為

C_{i} = e_{ijk}A_jB_k\\

\mathrm{\mathbf{I\!V.}}總結(jié)

接下來是逆時針方向的主動轉(zhuǎn)動變換的公式總結(jié)内颗。

\bullet轉(zhuǎn)動公式

\mathbf{r}^{\prime} = \mathbf{r}\cos\Phi + \mathbf{n}(\mathbf{n} \boldsymbol{\cdot} \mathbf{r})(1 - \cos\Phi) + (\mathbf{n} \times \mathbf{r}) \sin\Phi\\

(注意與順時針轉(zhuǎn)動公式的區(qū)別)

\bullet無限小轉(zhuǎn)動:\sin\Phi \rightarrow d\Phi钧排,\cos\Phi \rightarrow 1d\mathbf{r}^{\prime} \equiv \mathbf{r}^{\prime} - \mathbf{r}

d\mathbf{r}^{\prime} = (\mathbf{n}\;d\Phi)\times \mathbf{r} = d\mathbf{\Omega} \times \mathbf{r} \\

\bullet無限小算子\boldsymbol{\epsilon}為反對稱矩陣均澳,使用關(guān)系d\Omega_i = n_i \;d\Phi恨溜,可將其表示為

\boldsymbol{\epsilon} = \begin{bmatrix}0 & -d\Omega_3 & d\Omega_2\\d\Omega_3 & 0 & -d\Omega_1\\-d\Omega_2 & d\Omega_1 & 0\end{bmatrix} = \begin{bmatrix}0 & -n_3 & n_2\\n_3 & 0 & -n_1\\-n_2 & n_1 & 0\end{bmatrix}d\Phi \\

對于無限小轉(zhuǎn)動變換,\mathbf{r}^{\prime} = (\rm{I} + \boldsymbol{\epsilon} )\mathbf{r}负懦,

d\mathbf{r} = \boldsymbol{\epsilon}\mathbf{r} = \begin{bmatrix}0 & -n_3 & n_2\\n_3 & 0 & -n_1\\-n_2 & n_1 & 0\end{bmatrix}d\Phi \;\mathbf{r}\\

\implies \frac{d\mathbf{r}} {d\Phi} = -\rm{N}\mathbf{r}\\

其中\rm{N}是一個反對稱矩陣筒捺,N_{ij} = e_{ijk}n_k

\rm{N} = \begin{bmatrix}0 & n_3 & -n_2\\-n_3 & 0 & n_1\\n_2 & -n_1 & 0\end{bmatrix}\\

\bullet另一種比較有用的表征為

\boldsymbol{\epsilon} = n_i \mathbf{M}_i \;d\Phi \\

其中\mathbf{M}_i無限小轉(zhuǎn)動生成算子(infinitesimal rotation generator),它們分別為

\mathbf{M}_1 = \begin{bmatrix}0& 0& 0\\0&0&-1\\0&1& 0\end{bmatrix}纸厉,\mathbf{M}_2 = \begin{bmatrix}0& 0& 1\\0&0&0\\-1&0& 0\end{bmatrix}\mathbf{M}_3 = \begin{bmatrix}0& -1& 0\\1&0&0\\0&0& 0\end{bmatrix}

這些生成算子之間具有性質(zhì)

\mathbf{M}_i\mathbf{M}_j - \mathbf{M}_j\mathbf{M}_i \equiv \left[ \mathbf{M}_i, \mathbf{M}_j\right] = e_{ijk}\mathbf{M}_k\\

對易子\left[ \mathbf{M}_i, \mathbf{M}_j\right]也屬于李括號(Lie bracket)五嫂,因此用來定義含轉(zhuǎn)動參量的旋轉(zhuǎn)群群的李代數(shù)颗品。

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末肯尺,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子躯枢,更是在濱河造成了極大的恐慌则吟,老刑警劉巖,帶你破解...
    沈念sama閱讀 221,273評論 6 515
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件锄蹂,死亡現(xiàn)場離奇詭異氓仲,居然都是意外死亡,警方通過查閱死者的電腦和手機得糜,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 94,349評論 3 398
  • 文/潘曉璐 我一進(jìn)店門敬扛,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人朝抖,你說我怎么就攤上這事啥箭。” “怎么了治宣?”我有些...
    開封第一講書人閱讀 167,709評論 0 360
  • 文/不壞的土叔 我叫張陵急侥,是天一觀的道長。 經(jīng)常有香客問我侮邀,道長坏怪,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 59,520評論 1 296
  • 正文 為了忘掉前任绊茧,我火速辦了婚禮铝宵,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘按傅。我一直安慰自己捉超,他們只是感情好,可當(dāng)我...
    茶點故事閱讀 68,515評論 6 397
  • 文/花漫 我一把揭開白布唯绍。 她就那樣靜靜地躺著拼岳,像睡著了一般。 火紅的嫁衣襯著肌膚如雪况芒。 梳的紋絲不亂的頭發(fā)上惜纸,一...
    開封第一講書人閱讀 52,158評論 1 308
  • 那天,我揣著相機與錄音绝骚,去河邊找鬼耐版。 笑死,一個胖子當(dāng)著我的面吹牛压汪,可吹牛的內(nèi)容都是我干的粪牲。 我是一名探鬼主播,決...
    沈念sama閱讀 40,755評論 3 421
  • 文/蒼蘭香墨 我猛地睜開眼止剖,長吁一口氣:“原來是場噩夢啊……” “哼腺阳!你這毒婦竟也來了落君?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 39,660評論 0 276
  • 序言:老撾萬榮一對情侶失蹤亭引,失蹤者是張志新(化名)和其女友劉穎绎速,沒想到半個月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體焙蚓,經(jīng)...
    沈念sama閱讀 46,203評論 1 319
  • 正文 獨居荒郊野嶺守林人離奇死亡纹冤,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 38,287評論 3 340
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了购公。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片萌京。...
    茶點故事閱讀 40,427評論 1 352
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖君丁,靈堂內(nèi)的尸體忽然破棺而出枫夺,到底是詐尸還是另有隱情,我是刑警寧澤绘闷,帶...
    沈念sama閱讀 36,122評論 5 349
  • 正文 年R本政府宣布橡庞,位于F島的核電站,受9級特大地震影響印蔗,放射性物質(zhì)發(fā)生泄漏扒最。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 41,801評論 3 333
  • 文/蒙蒙 一华嘹、第九天 我趴在偏房一處隱蔽的房頂上張望吧趣。 院中可真熱鬧,春花似錦耙厚、人聲如沸强挫。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,272評論 0 23
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽俯渤。三九已至,卻和暖如春型宝,著一層夾襖步出監(jiān)牢的瞬間八匠,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 33,393評論 1 272
  • 我被黑心中介騙來泰國打工趴酣, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留梨树,地道東北人。 一個月前我還...
    沈念sama閱讀 48,808評論 3 376
  • 正文 我出身青樓岖寞,卻偏偏與公主長得像抡四,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子仗谆,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 45,440評論 2 359

推薦閱讀更多精彩內(nèi)容

  • 轉(zhuǎn)載自VR設(shè)計云課堂[http://www.reibang.com/u/c7ffdc4b379e]Unity S...
    水月凡閱讀 1,018評論 0 0
  • 本征值問題和久期方程 在任意時刻床嫌,剛體的方向可由正交變換來表示跨释。時間的演進(jìn)會導(dǎo)致剛體方向的變化胸私,所以它的變換矩陣該...
    有限與微小的面包閱讀 5,297評論 0 2
  • 前言 在三維空間中厌处,為了指定一個剛體的位置和姿態(tài),至少需要6個數(shù)岁疼。本章通過一個附著到物體上的參考坐標(biāo)系來描述它的位...
    play_robot閱讀 3,077評論 0 2
  • Android啟動流程 個人理解的Android啟動流程阔涉,有沒有大佬幫忙看下對不對。 1. Loader層 Boo...
    懶貓1105閱讀 693評論 0 2
  • 朋友Z在我耳邊滔滔不絕地談著她的大學(xué)生活:室友去國家大劇院表演話劇捷绒、班級里不少的狀元以及與北京其他優(yōu)秀高校的項目合...
    navy_閱讀 315評論 0 2