证明旋转矩阵R的指数映射

旋转矩阵如 R = [ r 11 r 12 r 13 r 21 r 22 r 23 r 31 r 32 r 33 ] \boldsymbol{R}=\left[\begin{array}{lll}r_{11} & r_{12} & r_{13} \\ r_{21} & r_{22} & r_{23} \\ r_{31} & r_{32} & r_{33}\end{array}\right] R= r11r21r31r12r22r32r13r23r33 ,对于给定的 R ∈ S O ( 3 ) R \in SO(3) RSO(3),存在 ω ∈ ℜ 3 , ∥ ω ∥ = 1 \omega \in \Re^{3}, \quad\|\omega\|=1 ω3,ω=1 θ ∈ ℜ \theta \in \Re θ,使得 R = exp ⁡ ( ω ^ θ ) ∘ \boldsymbol{R}=\exp (\hat{\boldsymbol{\omega}} \theta)_{\circ} R=exp(ω^θ)

证明: log ⁡ ( R ) = { θ 2 sin ⁡ ( θ ) ( R − R T ) ,  if  θ ≠ 0 0 ,  if  θ = 0 \log (\mathbf{R})=\left\{\begin{array}{ll}\frac{\theta}{2 \sin (\theta)}\left(\mathbf{R}-\mathbf{R}^{T}\right), & \text { if } \theta \neq 0 \\ 0, & \text { if } \theta=0\end{array}\right. log(R)={2sin(θ)θ(RRT),0, if θ=0 if θ=0

证明:
根据定义 R = [ r 11 r 12 r 13 r 21 r 22 r 23 r 31 r 32 r 33 ] \boldsymbol{R}=\left[\begin{array}{lll}r_{11} & r_{12} & r_{13} \\ r_{21} & r_{22} & r_{23} \\ r_{31} & r_{32} & r_{33}\end{array}\right] R= r11r21r31r12r22r32r13r23r33 ,令 v θ = 1 − cos ⁡ θ , c 9 = cos ⁡ θ , s θ = sin ⁡ θ v_{\theta}=1-\cos \theta, \quad c_{9}=\cos \theta, \quad s_{\theta}=\sin \theta vθ=1cosθ,c9=cosθ,sθ=sinθ,则罗德里格斯公式可以改写如下形式:

e ω θ ˙ = I + ω ^ sin ⁡ θ + ω ^ 2 ( 1 − cos ⁡ θ ) = [ 1 − v θ ( ω 2 2 + ω 3 2 ) ω 1 ω 2 v θ − ω 3 s θ ω 1 ω 3 v θ + ω 2 s θ ω 1 ω 2 v θ + ω 3 s θ 1 − v θ ( ω 1 2 + ω 3 2 ) ω 2 ω 3 v θ − ω 1 s θ ω 1 ω 3 v θ − ω 2 s θ ω 2 ω 3 v θ + ω 1 s θ 1 − v θ ( ω 1 2 + ω 2 2 ) ] = [ ω 1 2 v θ + c θ ω 1 ω 2 v θ − ω 3 s θ ω 1 ω 3 v θ + ω 2 s θ ω 1 ω 2 v θ + ω 3 s θ ω 2 2 v θ + c θ ω 2 ω 3 v θ − ω 1 s θ ω 1 ω 3 v θ − ω 2 s θ ω 2 ω 3 v θ + ω 1 s θ ω 3 2 v θ + c θ ] \begin{aligned} \mathrm{e}^{\dot{\omega \theta}} & =\boldsymbol{I}+\hat{\boldsymbol{\omega}} \sin \theta+\hat{\boldsymbol{\omega}}^{2}(1-\cos \theta)= \\ & {\left[\begin{array}{ccc}1-v_{\theta}\left(\omega_{2}^{2}+\omega_{3}^{2}\right) & \omega_{1} \omega_{2} v_{\theta}-\omega_{3} s_{\theta} & \omega_{1} \omega_{3} v_{\theta}+\omega_{2} s_{\theta} \\ \omega_{1} \omega_{2} v_{\theta}+\omega_{3} s_{\theta} & 1-v_{\theta}\left(\omega_{1}^{2}+\omega_{3}^{2}\right) & \omega_{2} \omega_{3} v_{\theta}-\omega_{1} s_{\theta} \\ \omega_{1} \omega_{3} v_{\theta}-\omega_{2} s_{\theta} & \omega_{2} \omega_{3} v_{\theta}+\omega_{1} s_{\theta} & 1-v_{\theta}\left(\omega_{1}^{2}+\omega_{2}^{2}\right)\end{array}\right]=} \\ & {\left[\begin{array}{ccc}\omega_{1}^{2} v_{\theta}+c_{\theta} & \omega_{1} \omega_{2} v_{\theta}-\omega_{3} s_{\theta} & \omega_{1} \omega_{3} v_{\theta}+\omega_{2} s_{\theta} \\ \omega_{1} \omega_{2} v_{\theta}+\omega_{3} s_{\theta} & \omega_{2}^{2} v_{\theta}+c_{\theta} & \omega_{2} \omega_{3} v_{\theta}-\omega_{1} s_{\theta} \\ \omega_{1} \omega_{3} v_{\theta}-\omega_{2} s_{\theta} & \omega_{2} \omega_{3} v_{\theta}+\omega_{1} s_{\theta} & \omega_{3}^{2} v_{\theta}+c_{\theta}\end{array}\right] }\end{aligned} eωθ˙=I+ω^sinθ+ω^2(1cosθ)= 1vθ(ω22+ω32)ω1ω2vθ+ω3sθω1ω3vθω2sθω1ω2vθω3sθ1vθ(ω12+ω32)ω2ω3vθ+ω1sθω1ω3vθ+ω2sθω2ω3vθω1sθ1vθ(ω12+ω22) = ω12vθ+cθω1ω2vθ+ω3sθω1ω3vθω2sθω1ω2vθω3sθω22vθ+cθω2ω3vθ+ω1sθω1ω3vθ+ω2sθω2ω3vθω1sθω32vθ+cθ

根据R的定义与改写的罗德里格斯公式,可得: trace ⁡ ( R ) = r 11 + r 22 + r 33 = 1 + 2 cos ⁡ θ \operatorname{trace}(\boldsymbol{R})=r_{11}+r_{22}+r_{33}=1+2 \cos \theta trace(R)=r11+r22+r33=1+2cosθ

为证明该式有解,可利用R的迹等于其特征值之和的关系。由于R能保持长度不变且det®=+1,其特征值为1,并且为复共轭。由 − 1 ⩽ t r a c e ( R ) ⩽ 3 -1 \leqslant trace (\boldsymbol{R}) \leqslant 3 1trace(R)3,得:

θ = arccos ⁡ ( trace ⁡ ( R ) − 1 2 ) \theta=\arccos \left(\frac{\operatorname{trace}(\boldsymbol{R})-1}{2}\right) θ=arccos(2trace(R)1),其中 θ \theta θ值不确定,它可在 θ ± 2 π n 或 − θ ± 2 π n \theta \pm 2 \pi n 或 -\theta \pm 2 \pi n θ±2πnθ±2πn

再令R和 exp ⁡ ( ω ^ θ ) \exp (\hat{\boldsymbol{\omega}} \theta) exp(ω^θ)的非对角线元素相等,有

r 32 − r 23 = 2 ω 1 s θ r 13 − r 31 = 2 ω 2 s θ r 21 − r 12 = 2 ω 3 s θ \begin{array}{l}r_{32}-r_{23}=2 \omega_{1} s_{\theta} \\ r_{13}-r_{31}=2 \omega_{2} s_{\theta} \\ r_{21}-r_{12}=2 \omega_{3} s_{\theta}\end{array} r32r23=2ω1sθr13r31=2ω2sθr21r12=2ω3sθ

θ ≠ 0 \theta \neq 0 θ=0,则取

ω ^ = 1 2 s θ [ r 32 − r 23 r 13 − r 31 r 21 − r 12 ] \hat{\boldsymbol{\omega}} =\frac{1}{2 s_{\theta}}\left[\begin{array}{l}r_{32}-r_{23} \\ r_{13}-r_{31} \\ r_{21}-r_{12}\end{array}\right] ω^=2sθ1 r32r23r13r31r21r12

1 2 sin ⁡ ( θ ) ( R − R T ) = [ 0 − w 3 w 2 w 3 0 − w 1 − w 2 w 1 0 ] = ω ^ ∧ \frac{1}{2 \sin (\theta)}\left(\mathbf{R}-\mathbf{R}^{T}\right)=\left[\begin{array}{ccc}0 & -w_{3} & w_{2} \\ w_{3} & 0 & -w_{1} \\ -w_{2} & w_{1} & 0\end{array}\right]=\hat{\boldsymbol{\omega}} ^{\wedge} 2sin(θ)1(RRT)= 0w3w2w30w1w2w10 =ω^

所以,旋转矩阵R的对数映射 log ⁡ ( R ) = { θ 2 sin ⁡ ( θ ) ( R − R T ) ,  if  θ ≠ 0 0 ,  if  θ = 0 \log (\mathbf{R})=\left\{\begin{array}{ll}\frac{\theta}{2 \sin (\theta)}\left(\mathbf{R}-\mathbf{R}^{T}\right), & \text { if } \theta \neq 0 \\ 0, & \text { if } \theta=0\end{array}\right. log(R)={2sin(θ)θ(RRT),0, if θ=0 if θ=0

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值