当前位置: 首页 > news >正文

[Mv]_× = M [v]_× M^T的证明

假设 M M M 是一个 3 × 3 3 \times 3 3×3 的旋转矩阵, v = ( v 1 , v 2 , v 3 ) T v = (v_1, v_2, v_3)^T v=(v1,v2,v3)T 是一个三维向量。
下面给出 [ M v ] × = M [ v ] × M T [Mv]_\times = M [v]_\times M^T [Mv]×=M[v]×MT 的推导

1. [ M v ] × [Mv]_\times [Mv]×

计算 M v Mv Mv
M v = [ m 11 m 12 m 13 m 21 m 22 m 23 m 31 m 32 m 33 ] [ v 1 v 2 v 3 ] = [ m 11 v 1 + m 12 v 2 + m 13 v 3 m 21 v 1 + m 22 v 2 + m 23 v 3 m 31 v 1 + m 32 v 2 + m 33 v 3 ] Mv = \begin{bmatrix} m_{11} & m_{12} & m_{13} \\ m_{21} & m_{22} & m_{23} \\ m_{31} & m_{32} & m_{33} \end{bmatrix} \begin{bmatrix} v_1 \\ v_2 \\ v_3 \end{bmatrix} = \begin{bmatrix} m_{11}v_1 + m_{12}v_2 + m_{13}v_3 \\ m_{21}v_1 + m_{22}v_2 + m_{23}v_3 \\ m_{31}v_1 + m_{32}v_2 + m_{33}v_3 \end{bmatrix} Mv= m11m21m31m12m22m32m13m23m33 v1v2v3 = m11v1+m12v2+m13v3m21v1+m22v2+m23v3m31v1+m32v2+m33v3
M v = ( u 1 , u 2 , u 3 ) T Mv = (u_1, u_2, u_3)^T Mv=(u1,u2,u3)T,则:

[ u ] × = [ 0 − u 3 u 2 u 3 0 − u 1 − u 2 u 1 0 ] [u]_\times = \begin{bmatrix} 0 & -u_3 & u_2 \\ u_3 & 0 & -u_1 \\ -u_2 & u_1 & 0 \end{bmatrix} [u]×= 0u3u2u30u1u2u10

具体形式为:

[ M v ] × = [ 0 − ( m 31 v 1 + m 32 v 2 + m 33 v 3 ) ( m 21 v 1 + m 22 v 2 + m 23 v 3 ) ( m 31 v 1 + m 32 v 2 + m 33 v 3 ) 0 − ( m 11 v 1 + m 12 v 2 + m 13 v 3 ) − ( m 21 v 1 + m 22 v 2 + m 23 v 3 ) ( m 11 v 1 + m 12 v 2 + m 13 v 3 ) 0 ] [Mv]_\times = \begin{bmatrix} 0 & -(m_{31}v_1 + m_{32}v_2 + m_{33}v_3) & (m_{21}v_1 + m_{22}v_2 + m_{23}v_3) \\ (m_{31}v_1 + m_{32}v_2 + m_{33}v_3) & 0 & -(m_{11}v_1 + m_{12}v_2 + m_{13}v_3) \\ -(m_{21}v_1 + m_{22}v_2 + m_{23}v_3) & (m_{11}v_1 + m_{12}v_2 + m_{13}v_3) & 0 \end{bmatrix} [Mv]×= 0(m31v1+m32v2+m33v3)(m21v1+m22v2+m23v3)(m31v1+m32v2+m33v3)0(m11v1+m12v2+m13v3)(m21v1+m22v2+m23v3)(m11v1+m12v2+m13v3)0

2. M [ v ] × M T M [v]_\times M^T M[v]×MT

再看 ( M [v]_\times M^T ):

[ v ] × = [ 0 − v 3 v 2 v 3 0 − v 1 − v 2 v 1 0 ] [v]_\times = \begin{bmatrix} 0 & -v_3 & v_2 \\ v_3 & 0 & -v_1 \\ -v_2 & v_1 & 0 \end{bmatrix} [v]×= 0v3v2v30v1v2v10

计算 ( M [v]_\times M^T ),其中 ( M^T ) 是 ( M ) 的转置:

M [ v ] × M T = [ m 11 m 12 m 13 m 21 m 22 m 23 m 31 m 32 m 33 ] [ 0 − v 3 v 2 v 3 0 − v 1 − v 2 v 1 0 ] [ m 11 m 21 m 31 m 12 m 22 m 32 m 13 m 23 m 33 ] M [v]_\times M^T = \begin{bmatrix} m_{11} & m_{12} & m_{13} \\ m_{21} & m_{22} & m_{23} \\ m_{31} & m_{32} & m_{33} \end{bmatrix} \begin{bmatrix} 0 & -v_3 & v_2 \\ v_3 & 0 & -v_1 \\ -v_2 & v_1 & 0 \end{bmatrix} \begin{bmatrix} m_{11} & m_{21} & m_{31} \\ m_{12} & m_{22} & m_{32} \\ m_{13} & m_{23} & m_{33} \end{bmatrix} M[v]×MT= m11m21m31m12m22m32m13m23m33 0v3v2v30v1v2v10 m11m12m13m21m22m23m31m32m33

具体计算每一项:

( 1 , 1 ) : ( − v 3 m 12 + v 2 m 13 ) m 11 + ( v 3 m 11 − v 1 m 13 ) m 21 + ( − v 2 m 11 + v 1 m 12 ) m 31 = 0 (1,1): (-v_3 m_{12} + v_2 m_{13}) m_{11} + (v_3 m_{11} - v_1 m_{13}) m_{21} + (-v_2 m_{11} + v_1 m_{12}) m_{31} = 0 (1,1):(v3m12+v2m13)m11+(v3m11v1m13)m21+(v2m11+v1m12)m31=0

( 1 , 2 ) : ( − v 3 m 12 + v 2 m 13 ) m 12 + ( v 3 m 11 − v 1 m 13 ) m 22 + ( − v 2 m 11 + v 1 m 12 ) m 32 = − ( m 31 v 1 + m 32 v 2 + m 33 v 3 ) (1,2): (-v_3 m_{12} + v_2 m_{13}) m_{12} + (v_3 m_{11} - v_1 m_{13}) m_{22} + (-v_2 m_{11} + v_1 m_{12}) m_{32} = -(m_{31}v_1 + m_{32}v_2 + m_{33}v_3) (1,2):(v3m12+v2m13)m12+(v3m11v1m13)m22+(v2m11+v1m12)m32=(m31v1+m32v2+m33v3)

( 1 , 3 ) : ( − v 3 m 12 + v 2 m 13 ) m 13 + ( v 3 m 11 − v 1 m 13 ) m 23 + ( − v 2 m 11 + v 1 m 12 ) m 33 = ( m 21 v 1 + m 22 v 2 + m 23 v 3 ) (1,3): (-v_3 m_{12} + v_2 m_{13}) m_{13} + (v_3 m_{11} - v_1 m_{13}) m_{23} + (-v_2 m_{11} + v_1 m_{12}) m_{33} = (m_{21}v_1 + m_{22}v_2 + m_{23}v_3) (1,3):(v3m12+v2m13)m13+(v3m11v1m13)m23+(v2m11+v1m12)m33=(m21v1+m22v2+m23v3)

( 2 , 1 ) : ( v 3 m 11 − v 1 m 13 ) m 11 + ( − v 3 m 12 + v 2 m 13 ) m 21 + ( v 3 m 11 − v 1 m 13 ) m 31 = ( m 31 v 1 + m 32 v 2 + m 33 v 3 ) (2,1): (v_3 m_{11} - v_1 m_{13}) m_{11} + (-v_3 m_{12} + v_2 m_{13}) m_{21} + (v_3 m_{11} - v_1 m_{13}) m_{31} = (m_{31}v_1 + m_{32}v_2 + m_{33}v_3) (2,1):(v3m11v1m13)m11+(v3m12+v2m13)m21+(v3m11v1m13)m31=(m31v1+m32v2+m33v3)

( 2 , 2 ) : ( v 3 m 11 − v 1 m 13 ) m 12 + ( − v 3 m 12 + v 2 m 13 ) m 22 + ( v 3 m 11 − v 1 m 13 ) m 32 = 0 (2,2): (v_3 m_{11} - v_1 m_{13}) m_{12} + (-v_3 m_{12} + v_2 m_{13}) m_{22} + (v_3 m_{11} - v_1 m_{13}) m_{32} = 0 (2,2):(v3m11v1m13)m12+(v3m12+v2m13)m22+(v3m11v1m13)m32=0

( 2 , 3 ) : ( v 3 m 11 − v 1 m 13 ) m 13 + ( − v 3 m 12 + v 2 m 13 ) m 23 + ( v 3 m 11 − v 1 m 13 ) m 33 = − ( m 11 v 1 + m 12 v 2 + m 13 v 3 ) (2,3): (v_3 m_{11} - v_1 m_{13}) m_{13} + (-v_3 m_{12} + v_2 m_{13}) m_{23} + (v_3 m_{11} - v_1 m_{13}) m_{33} = -(m_{11}v_1 + m_{12}v_2 + m_{13}v_3) (2,3):(v3m11v1m13)m13+(v3m12+v2m13)m23+(v3m11v1m13)m33=(m11v1+m12v2+m13v3)

( 3 , 1 ) : ( − v 2 m 11 + v 1 m 12 ) m 11 + ( v 3 m 11 − v 1 m 13 ) m 21 + ( − v 2 m 11 + v 1 m 12 ) m 31 = − ( m 21 v 1 + m 22 v 2 + m 23 v 3 ) (3,1): (-v_2 m_{11} + v_1 m_{12}) m_{11} + (v_3 m_{11} - v_1 m_{13}) m_{21} + (-v_2 m_{11} + v_1 m_{12}) m_{31} = -(m_{21}v_1 + m_{22}v_2 + m_{23}v_3) (3,1):(v2m11+v1m12)m11+(v3m11v1m13)m21+(v2m11+v1m12)m31=(m21v1+m22v2+m23v3)

( 3 , 2 ) : ( − v 2 m 11 + v 1 m 12 ) m 12 + ( v 3 m 11 − v 1 m 13 ) m 22 + ( − v 2 m 11 + v 1 m 12 ) m 32 = ( m 11 v 1 + m 12 v 2 + m 13 v 3 ) (3,2): (-v_2 m_{11} + v_1 m_{12}) m_{12} + (v_3 m_{11} - v_1 m_{13}) m_{22} + (-v_2 m_{11} + v_1 m_{12}) m_{32} = (m_{11}v_1 + m_{12}v_2 + m_{13}v_3) (3,2):(v2m11+v1m12)m12+(v3m11v1m13)m22+(v2m11+v1m12)m32=(m11v1+m12v2+m13v3)

( 3 , 3 ) : ( − v 2 m 11 + v 1 m 12 ) m 13 + ( v 3 m 11 − v 1 m 13 ) m 23 + ( − v 2 m 11 + v 1 m 12 ) m 33 = 0 (3,3): (-v_2 m_{11} + v_1 m_{12}) m_{13} + (v_3 m_{11} - v_1 m_{13}) m_{23} + (-v_2 m_{11} + v_1 m_{12}) m_{33} = 0 (3,3):(v2m11+v1m12)m13+(v3m11v1m13)m23+(v2m11+v1m12)m33=0

整理后的结果为:

M [ v ] × M T = [ 0 − ( m 31 v 1 + m 32 v 2 + m 33 v 3 ) ( m 21 v 1 + m 22 v 2 + m 23 v 3 ) ( m 31 v 1 + m 32 v 2 + m 33 v 3 ) 0 − ( m 11 v 1 + m 12 v 2 + m 13 v 3 ) − ( m 21 v 1 + m 22 v 2 + m 23 v 3 ) ( m 11 v 1 + m 12 v 2 + m 13 v 3 ) 0 ] M [v]_\times M^T = \begin{bmatrix} 0 & -(m_{31}v_1 + m_{32}v_2 + m_{33}v_3) & (m_{21}v_1 + m_{22}v_2 + m_{23}v_3) \\ (m_{31}v_1 + m_{32}v_2 + m_{33}v_3) & 0 & -(m_{11}v_1 + m_{12}v_2 + m_{13}v_3) \\ -(m_{21}v_1 + m_{22}v_2 + m_{23}v_3) & (m_{11}v_1 + m_{12}v_2 + m_{13}v_3) & 0 \end{bmatrix} M[v]×MT= 0(m31v1+m32v2+m33v3)(m21v1+m22v2+m23v3)(m31v1+m32v2+m33v3)0(m11v1+m12v2+m13v3)(m21v1+m22v2+m23v3)(m11v1+m12v2+m13v3)0

总结

经过计算,可以看到 [ M v ] × 和 M [ v ] × M T [Mv]_\times 和 M [v]_\times M^T [Mv]×M[v]×MT 的形式确实是一致的,因此等式 [ M v ] × = M [ v ] × M T [Mv]_\times = M [v]_\times M^T [Mv]×=M[v]×MT 是成立的。


http://www.mrgr.cn/news/60317.html

相关文章:

  • vue2+elementui日期选择器
  • Mysql之视图创建
  • autoware-传感器驱动和感知算法(sensing-perception)-笔记
  • 关于嵌入式学习的一些短浅经验
  • 分账系统适用于那些商家?
  • gateway 整合 spring security oauth2
  • flowable 去掉自带的登录权限
  • 基本法代码阅读
  • UDS诊断刷写–跳转执行
  • Sigrity Power SI Noise coupling analysis模式如何观测PDN系统的交流耦合噪声以及平面间电压波动操作指导(二)
  • 未定义项目JDK 解决办法
  • C语言数据结构学习:单链表
  • python3的基本数据类型:Bool(布尔类型)
  • Screen简介和使用
  • Notepad++如何同时检索多个关键字
  • Python自动化个人健康日志与建议:记录日常活动并提供建议
  • Javascript高级:正则表达式基础与应用
  • OOP特性 多态
  • c语言错题——#define对应的查找替换
  • 【JavaEE】【多线程】线程池
  • WebGIS开发之编辑功能(分割、融合、捕捉、追踪)
  • 2023年中国县域统计年鉴(县市卷+乡镇卷)(excel格式)
  • 声学气膜馆:品牌发布会的理想之选—轻空间
  • Labview通讯测试耗时
  • 支持向量机SVM简述
  • Spring Boot摄影工作室:构建Web版在线服务平台