叉乘的矩阵形式

                     

贡献者: addis

预备知识 1 矢量叉乘

   对任意矢量 $ \boldsymbol{\mathbf{a}} $ 和 $ \boldsymbol{\mathbf{b}} $,令

\begin{equation} \boldsymbol{\mathbf{c}} = \boldsymbol{\mathbf{a}} \boldsymbol\times \boldsymbol{\mathbf{b}} ~. \end{equation}
该运算可以看作列矢量 $ \boldsymbol{\mathbf{b}} $ 到列矢量 $ \boldsymbol{\mathbf{c}} $ 的线性变换。我们知道线性变换可以用矩阵表示,所以必存在矩阵 $ \boldsymbol{\mathbf{A}} $,满足
\begin{equation} \boldsymbol{\mathbf{c}} = \boldsymbol{\mathbf{A}} \boldsymbol{\mathbf{b}} ~. \end{equation}
令 $ \boldsymbol{\mathbf{a}} $ 的坐标为 $(a_x, a_y, a_z)$,根据叉乘的分量表达式(式 12 ),易得变换矩阵为
\begin{equation} \boldsymbol{\mathbf{A}} = \begin{pmatrix} 0 & -a_z & a_y\\ a_z & 0 & -a_x\\ -a_y & a_x & 0 \end{pmatrix} ~\end{equation}
这是一个反对称矩阵,即
\begin{equation} \boldsymbol{\mathbf{A}} ^{\mathrm{T}} = - \boldsymbol{\mathbf{A}} ~. \end{equation}

   类似地,也有

\begin{equation} \boldsymbol{\mathbf{c}} ^{\mathrm{T}} = \boldsymbol{\mathbf{a}} ^{\mathrm{T}} \boldsymbol{\mathbf{B}} ~, \end{equation}
其中 $ \boldsymbol{\mathbf{B}} $ 的定义和 $ \boldsymbol{\mathbf{A}} $ 一致。

   同理,式 1 也可以看作是 $ \boldsymbol{\mathbf{a}} $ 到 $ \boldsymbol{\mathbf{c}} $ 的线性变换

\begin{equation} \boldsymbol{\mathbf{c}} = \boldsymbol{\mathbf{B}} \boldsymbol{\mathbf{a}} ~. \end{equation}
其中
\begin{equation} \boldsymbol{\mathbf{B}} = \begin{pmatrix} 0 & b_z & -b_y\\ -b_z & 0 & b_x\\ b_y & -b_x & 0 \end{pmatrix} ~ \end{equation}
这恰好与式 3 符号相反。

1. 叉乘矩阵的旋转变换

预备知识 2 三维旋转矩阵

   令 $S'$ 坐标系与 $S$ 坐标系的原点重合,且 $S'$ 坐标系中的矢量 $ \boldsymbol{\mathbf{r}} $ 到 $S$ 系中的矢量 $ \boldsymbol{\mathbf{r}} '$ 的旋转变换矩阵为 $ \boldsymbol{\mathbf{R}} $,满足 $ \boldsymbol{\mathbf{r}} ' = \boldsymbol{\mathbf{R}} \boldsymbol{\mathbf{r}} $。令两个几何矢量在 $S$ 系中的坐标为 $ \boldsymbol{\mathbf{u}} , \boldsymbol{\mathbf{v}} $,在 $S'$ 中坐标为 $ \boldsymbol{\mathbf{u}} ', \boldsymbol{\mathbf{v}} '$。那么有

\begin{equation} \boldsymbol{\mathbf{u}} ' = \boldsymbol{\mathbf{R}} \boldsymbol{\mathbf{u}} , \qquad \boldsymbol{\mathbf{v}} ' = \boldsymbol{\mathbf{R}} \boldsymbol{\mathbf{v}} ~. \end{equation}
由于叉乘的几何定义不依赖于坐标系,那么必定有
\begin{equation} \boldsymbol{\mathbf{R}} ( \boldsymbol{\mathbf{u}} \boldsymbol\times \boldsymbol{\mathbf{v}} ) = \boldsymbol{\mathbf{u}} ' \boldsymbol\times \boldsymbol{\mathbf{v}} '~. \end{equation}

   如果用 $ \boldsymbol{\mathbf{U}} $ 来代替 $ \boldsymbol{\mathbf{u}} \boldsymbol\times $,用 $ \boldsymbol{\mathbf{U}} '$ 来代替 $ \boldsymbol{\mathbf{u}} ' \boldsymbol\times $,那么上式变为

\begin{equation} \boldsymbol{\mathbf{R}} \boldsymbol{\mathbf{U}} \boldsymbol{\mathbf{v}} = \boldsymbol{\mathbf{U}} ' \boldsymbol{\mathbf{v}} ' = \boldsymbol{\mathbf{U}} ' \boldsymbol{\mathbf{R}} \boldsymbol{\mathbf{v}} ~. \end{equation}
这对所有 $ \boldsymbol{\mathbf{v}} $ 都成立,所以 $ \boldsymbol{\mathbf{R}} \boldsymbol{\mathbf{U}} = \boldsymbol{\mathbf{U}} ' \boldsymbol{\mathbf{v}} ' = \boldsymbol{\mathbf{U}} ' \boldsymbol{\mathbf{R}} $,或者
\begin{equation} \boldsymbol{\mathbf{U}} ' = \boldsymbol{\mathbf{R}} \boldsymbol{\mathbf{U}} \boldsymbol{\mathbf{R}} ^{\mathrm{T}} ~. \end{equation}


致读者: 小时百科一直以来坚持所有内容免费,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此,我们请求广大读者热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 10 元,我们一个星期内就能脱离亏损, 并保证在接下来的一整年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。

                     

友情链接: 超理论坛 | ©小时科技 保留一切权利