罗德里格斯旋转公式证明-简洁

news/2024/10/18 14:16:14/

罗德里格斯旋转公式证明。

设旋转向量为 ( n , θ ) (n, \theta) (n,θ),设其对应的旋转矩阵为 R R R

如何证明?
R = c o s θ I + n ∧ s i n θ + ( 1 − c o s θ ) n n T R=cos\theta I + n^{\wedge}sin\theta+(1-cos\theta)nn^{T} R=cosθI+nsinθ+(1cosθ)nnT
证明过程如下:

如图所示,设旋转向量为 A ^ \hat{A} A^,记为 n n n,设三维中的点 r r r n n n旋转 θ \theta θ后得到 r ′ r^{'} r,其中 n n n为单位方向向量,向量 n n n的起点为坐标原点。

image

r 3 r_3 r3为r在 n n n上的投影,则
r 3 = ( r ⋅ n ) n (1) r_3=(r\cdot n)n \tag{1} r3=(rn)n(1)
r 1 r_1 r1为r减去r在 n n n上面的分量 r 3 r_3 r3,则
r 1 = r − r 3 (2) r_1=r-r_3 \tag{2} r1=rr3(2)
r 2 r_2 r2 n n n r 1 r_1 r1的叉乘结果向量,则
r 2 = n × r 1 (3) r_2 = n\times r_1 \tag{3} r2=n×r1(3)
因此, r 1 , r 2 , r 3 r_1,r_2,r_3 r1r2r3构成了两两垂直的坐标系,但是模长不等于1, r 1 r_1 r1 r 2 r_2 r2模长相等。

由上图所示, r ′ r^{'} r r 1 和 r 2 r_1和r_2 r1r2所在的平面上的投影为 r ′ − r 3 r^{'}-r_3 rr3,则将其用 r 1 和 r 2 r_1和r_2 r1r2表示有
r ′ − r 3 = r 1 c o s θ + r 2 s i n θ r^{'}-r_3=r_1cos\theta+r_2sin\theta rr3=r1cosθ+r2sinθ
则,
r ′ = r 1 c o s θ + r 2 s i n θ + r 3 (4) r^{'}=r_1cos\theta+r_2sin\theta+r_3 \tag{4} r=r1cosθ+r2sinθ+r3(4)
综上所述,将(1)(2)(3)代入(4)式,则
r ′ = ( r − r 3 ) c o s θ + ( n × r 1 ) s i n θ + r 3 = r c o s θ + ( n × r 1 ) s i n θ + ( 1 − c o s θ ) r 3 = r c o s θ + ( n × ( r − r 3 ) ) s i n θ + ( 1 − c o s θ ) r 3 = r c o s θ + ( n × r − n × r 3 ) s i n θ + ( 1 − c o s θ ) r 3 ( 由于 n × r 3 = 0 ) = r c o s θ + n × r s i n θ + ( 1 − c o s θ ) r 3 = r c o s θ + n ∧ s i n θ ⋅ r + ( 1 − c o s θ ) ( r ⋅ n ) n = I c o s θ ⋅ r + n ∧ s i n θ ⋅ r + ( 1 − c o s θ ) n n T ⋅ r (5) \begin{aligned} r^{'} &=(r-r_3)cos\theta+(n\times r_1) sin\theta+r_3 \\ &=rcos\theta+(n\times r_1)sin\theta+(1-cos\theta)r_3 \\ &=rcos\theta+(n\times(r-r_3))sin\theta+(1-cos\theta)r_3 \\ &=rcos\theta+(n\times r-n\times r_3)sin\theta+(1-cos\theta)r_3 \space \space(由于n\times r_3 =0)\\ &=rcos\theta+n\times r sin\theta+(1-cos\theta)r_3 \\ &=rcos\theta+n^{\wedge}sin\theta \cdot r+(1-cos\theta)(r\cdot n)n \\ &=Icos\theta\cdot r+n^{\wedge}sin\theta\cdot r+(1-cos\theta) nn^{T}\cdot r \end{aligned} \tag{5} r=(rr3)cosθ+(n×r1)sinθ+r3=rcosθ+(n×r1)sinθ+(1cosθ)r3=rcosθ+(n×(rr3))sinθ+(1cosθ)r3=rcosθ+(n×rn×r3)sinθ+(1cosθ)r3  (由于n×r3=0)=rcosθ+n×rsinθ+(1cosθ)r3=rcosθ+nsinθr+(1cosθ)(rn)n=Icosθr+nsinθr+(1cosθ)nnTr(5)
设旋转矩阵为R,则 r ′ = R ⋅ r r^{'}=R\cdot r r=Rr,由公式(5)可知
R = I c o s θ + n ∧ s i n θ + ( 1 − c o s θ ) n n T R=Icos\theta+n^{\wedge}sin\theta+(1-cos\theta)nn^{T} R=Icosθ+nsinθ+(1cosθ)nnT
证明完毕。

参考链接:

1、https://wuli.wiki/online/RotA.html

2、https://en.wikipedia.org/wiki/Rodrigues%27_rotation_formula

3、https://sites.cs.ucsb.edu/~lingqi/teaching/resources/GAMES101_Lecture_04_supp.pdf


http://www.ppmy.cn/news/1466215.html

相关文章

Linux系统维护

1. 批量安装部署 2. 初始化配置 3. 禁用Selinux 永久更改 SELinux 配置: 编辑 SELinux 配置文件:使用文本编辑器打开 /etc/selinux/config 文件: 在配置文件中,找到 SELINUX… 的行。将其值更改为以下选项之一: e…

python练习五

Title1:请实现一个装饰器,每次调用函数时,将函数名字以及调用此函数的时间点写入文件中 代码: import time time time.strftime("%Y-%m-%d %H:%M:%S", time.localtime()) # 获取当前的时间戳 # 定义一个有参装饰器来实…

C++11 在 Windows 环境下的多线程编程指南

引言 随着多核处理器的普及,利用多线程编程来提升应用程序性能变得越来越重要。C11 标准库引入了一系列用于多线程编程的API,使得多线程编程变得更加简洁和高效。 一、基本概念 在开始编写多线程代码之前,了解一些基本概念是非常重要的&am…

linux学习:进程通信 管道

目录 例子1 父进程向子进程发送一条消息,子进程读取这条消息 例子2 mkfifo 函数创建一个命名管道 例子3 mkfifo 函数创建一个命名管道处理可能出现的错误 例子4 管道文件是否已存在 例子5 除了“文件已存在”进行处理 例子6 创建一个命名管道&…

2024年最适合中小型企业的顶级12款CRM系统对比

适合中小企业的CRM系统包括:纷享销客、Zoho CRM、HubSpot CRM、Freshsales、Pipedrive、Monday Sales CRM、Salesmate、Capsule CRM、金蝶云CRM、用友CRM、有赞CRM、简信CRM。 在现代商业运营中,CRM(客户关系管理)系统已成为不可或…

小白跟做江科大32单片机之旋转编码器计次

原理部分按照下面这个链接理解即可y小白跟做江科大32单片机之对射式红外传感器计次-CSDN博客https://blog.csdn.net/weixin_58051657/article/details/139350487https://blog.csdn.net/weixin_58051657/article/details/139350487 实验过程 1.按照江科大老师给的电路图进行连接…

MySQL Doublewrite Buffer 有了解过吗?

引言:在数据库管理中,确保数据的完整性和一致性是至关重要的。然而,在持久化数据到磁盘的过程中,可能会遇到各种意外情况,如断电或系统崩溃,从而导致部分数据写入,而另一部分数据未能成功写入&a…

数据结构 实验 1

题目一:用线性表实现文具店的货品管理问题 问题描述:在文具店的日常经营过程中,存在对各种文具的管理问题。当库存文具不足或缺货时,需要进货。日常销售时需要出库。当盘点货物时,需要查询货物信息。请根据这些要求编…