姿态矩阵/旋转矩阵/反对称阵

news/2025/2/27 0:21:09

在这里插入图片描述
物理意义,端点矢量=角速率叉乘本身向量;
负号是动系b看固定系i是相反的;

一个固定
在惯性导航解算中,旋转矢量的叉乘用于描述姿态矩阵的微分方程。你提到的公式中, ω i b b × \boldsymbol{\omega}_{ib}^b \times ωibb×表示的是一个向量叉乘的矩阵形式,通常称为叉乘矩阵反对称矩阵。以下是详细的解释和计算规则:


1. 叉乘矩阵的定义

对于任意三维向量 ω = [ ω x , ω y , ω z ] T \boldsymbol{\omega} = [\omega_x, \omega_y, \omega_z]^T ω=[ωx,ωy,ωz]T,其叉乘矩阵 ω × \boldsymbol{\omega} \times ω×定义为:
ω × = ( 0 − ω z ω y ω z 0 − ω x − ω y ω x 0 ) \boldsymbol{\omega} \times = \begin{pmatrix} 0 & -\omega_z & \omega_y \\ \omega_z & 0 & -\omega_x \\ -\omega_y & \omega_x & 0 \end{pmatrix} ω×= 0ωzωyωz0ωxωyωx0

这个矩阵的作用是将向量 ω \boldsymbol{\omega} ω的叉乘运算转化为矩阵乘法运算。对于任意向量 r \boldsymbol{r} r,有:
ω × r = ω × r \boldsymbol{\omega} \times \boldsymbol{r} = \boldsymbol{\omega} \times \boldsymbol{r} ω×r=ω×r

2. 姿态矩阵微分方程的推导

在惯性导航中,姿态矩阵 C b i \boldsymbol{C}_b^i Cbi 表示从载体坐标系( b b b 系)到惯性坐标系( i i i 系)的旋转矩阵。其微分方程为:
C ˙ b i = C b i ( ω i b b × ) \dot{\boldsymbol{C}}_b^i = \boldsymbol{C}_b^i (\boldsymbol{\omega}_{ib}^b \times) C˙bi=Cbi(ωibb×)
其中:

  • ω i b b \boldsymbol{\omega}_{ib}^b ωibb 是载体相对于惯性系的角速度,在载体坐标系( b b b 系)下表示。
  • ω i b b × \boldsymbol{\omega}_{ib}^b \times ωibb× 是角速度向量 ω i b b \boldsymbol{\omega}_{ib}^b ωibb的叉乘矩阵

3. 叉乘矩阵的作用

叉乘矩阵 ω i b b × \boldsymbol{\omega}_{ib}^b \times ωibb× 的作用是将角速度向量 ω i b b \boldsymbol{\omega}_{ib}^b ωibb的叉乘运算转化为矩阵乘法运算。具体来说:

  • 对于任意向量 r b \boldsymbol{r}^b rb,有:
    ω i b b × r b = ω i b b × r b \boldsymbol{\omega}_{ib}^b \times \boldsymbol{r}^b = \boldsymbol{\omega}_{ib}^b \times \boldsymbol{r}^b ωibb×rb=ωibb×rb
  • 对于矩阵 C b i \boldsymbol{C}_b^i Cbi,有:
    C ˙ b i = C b i ( ω i b b × ) \dot{\boldsymbol{C}}_b^i = \boldsymbol{C}_b^i (\boldsymbol{\omega}_{ib}^b \times) C˙bi=Cbi(ωibb×)

4. 具体计算规则

假设角速度向量 ω i b b = [ ω x , ω y , ω z ] T \boldsymbol{\omega}_{ib}^b = [\omega_x, \omega_y, \omega_z]^T ωibb=[ωx,ωy,ωz]T,则其叉乘矩阵为:
ω i b b × = ( 0 − ω z ω y ω z 0 − ω x − ω y ω x 0 ) \boldsymbol{\omega}_{ib}^b \times = \begin{pmatrix} 0 & -\omega_z & \omega_y \\ \omega_z & 0 & -\omega_x \\ -\omega_y & \omega_x & 0 \end{pmatrix} ωibb×= 0ωzωyωz0ωxωyωx0

姿态矩阵 C b i \boldsymbol{C}_b^i Cbi的微分方程为:
C ˙ b i = C b i ⋅ ( 0 − ω z ω y ω z 0 − ω x − ω y ω x 0 ) \dot{\boldsymbol{C}}_b^i = \boldsymbol{C}_b^i \cdot \begin{pmatrix} 0 & -\omega_z & \omega_y \\ \omega_z & 0 & -\omega_x \\ -\omega_y & \omega_x & 0 \end{pmatrix} C˙bi=Cbi 0ωzωyωz0ωxωyωx0

具体计算时, C ˙ b i \dot{\boldsymbol{C}}_b^i C˙bi的每一列是 C b i \boldsymbol{C}_b^i Cbi的对应列与叉乘矩阵的乘积。


5. 示例

假设:
ω i b b = ( ω x ω y ω z ) , C b i = ( c 11 c 12 c 13 c 21 c 22 c 23 c 31 c 32 c 33 ) \boldsymbol{\omega}_{ib}^b = \begin{pmatrix} \omega_x \\ \omega_y \\ \omega_z \end{pmatrix}, \quad \boldsymbol{C}_b^i = \begin{pmatrix} c_{11} & c_{12} & c_{13} \\ c_{21} & c_{22} & c_{23} \\ c_{31} & c_{32} & c_{33} \end{pmatrix} ωibb= ωxωyωz ,Cbi= c11c21c31c12c22c32c13c23c33

则:
C ˙ b i = C b i ⋅ ( 0 − ω z ω y ω z 0 − ω x − ω y ω x 0 ) \dot{\boldsymbol{C}}_b^i = \boldsymbol{C}_b^i \cdot \begin{pmatrix} 0 & -\omega_z & \omega_y \\ \omega_z & 0 & -\omega_x \\ -\omega_y & \omega_x & 0 \end{pmatrix} C˙bi=Cbi 0ωzωyωz0ωxωyωx0

计算结果为:
C ˙ b i = ( c 11 ⋅ 0 + c 12 ⋅ ω z + c 13 ⋅ ( − ω y ) c 11 ⋅ ( − ω z ) + c 12 ⋅ 0 + c 13 ⋅ ω x c 11 ⋅ ω y + c 12 ⋅ ( − ω x ) + c 13 ⋅ 0 c 21 ⋅ 0 + c 22 ⋅ ω z + c 23 ⋅ ( − ω y ) c 21 ⋅ ( − ω z ) + c 22 ⋅ 0 + c 23 ⋅ ω x c 21 ⋅ ω y + c 22 ⋅ ( − ω x ) + c 23 ⋅ 0 c 31 ⋅ 0 + c 32 ⋅ ω z + c 33 ⋅ ( − ω y ) c 31 ⋅ ( − ω z ) + c 32 ⋅ 0 + c 33 ⋅ ω x c 31 ⋅ ω y + c 32 ⋅ ( − ω x ) + c 33 ⋅ 0 ) \dot{\boldsymbol{C}}_b^i = \begin{pmatrix} c_{11} \cdot 0 + c_{12} \cdot \omega_z + c_{13} \cdot (-\omega_y) & c_{11} \cdot (-\omega_z) + c_{12} \cdot 0 + c_{13} \cdot \omega_x & c_{11} \cdot \omega_y + c_{12} \cdot (-\omega_x) + c_{13} \cdot 0 \\ c_{21} \cdot 0 + c_{22} \cdot \omega_z + c_{23} \cdot (-\omega_y) & c_{21} \cdot (-\omega_z) + c_{22} \cdot 0 + c_{23} \cdot \omega_x & c_{21} \cdot \omega_y + c_{22} \cdot (-\omega_x) + c_{23} \cdot 0 \\ c_{31} \cdot 0 + c_{32} \cdot \omega_z + c_{33} \cdot (-\omega_y) & c_{31} \cdot (-\omega_z) + c_{32} \cdot 0 + c_{33} \cdot \omega_x & c_{31} \cdot \omega_y + c_{32} \cdot (-\omega_x) + c_{33} \cdot 0 \end{pmatrix} C˙bi= c110+c12ωz+c13(ωy)c210+c22ωz+c23(ωy)c310+c32ωz+c33(ωy)c11(ωz)+c120+c13ωxc21(ωz)+c220+c23ωxc31(ωz)+c320+c33ωxc11ωy+c12(ωx)+c130c21ωy+c22(ωx)+c230c31ωy+c32(ωx)+c330


6. 物理意义

  • 叉乘矩阵 ω i b b × \boldsymbol{\omega}_{ib}^b \times ωibb× 描述了载体坐标系相对于惯性坐标系的旋转速率。
  • 姿态矩阵微分方程 C ˙ b i = C b i ( ω i b b × ) \dot{\boldsymbol{C}}_b^i = \boldsymbol{C}_b^i (\boldsymbol{\omega}_{ib}^b \times) C˙bi=Cbi(ωibb×)描述了姿态矩阵随时间的变化率。

总结

在惯性导航解算中,叉乘矩阵 ω i b b × \boldsymbol{\omega}_{ib}^b \times ωibb×是将角速度向量 ω i b b \boldsymbol{\omega}_{ib}^b ωibb 的叉乘运算转化为矩阵乘法运算的工具。通过姿态矩阵微分方程,可以实时更新姿态矩阵 C b i \boldsymbol{C}_b^i Cbi,从而解算载体的姿态变化。


http://www.niftyadmin.cn/n/5869308.html

相关文章

QT MD5校验文件和数据的完整性

MD5介绍: ‌MD5,全称为‌‌Message-Digest Algorithm 5(消息摘要算法5)‌,是一种广泛使用的密码散列函数。它可以将任意长度的“字节串”变换成一个128位(16字节)的散列值(hash valu…

【 C# 字符转换为 Unicode 编码】

在 C# 中,将字符转换为 Unicode 编码(即字符的 Unicode 码点)非常简单。每个字符在 C# 中都是一个 char 类型,而 char 类型本质上是一个 16 位的 Unicode 字符。你可以通过以下方法将字符转换为 Unicode 编码: 1. 使用…

JSPSSM526SpringBoot 教师人事档案管理系统功能说明

JSPSSM526SpringBoot 教师人事档案管理系统功能说明 JSPSSM526SpringBoot教师人事档案管理系统是一款专为教育机构设计的高效、安全、易用的信息化管理系统。该系统基于Spring Boot框架开发,结合JSP(JavaServer Pages)技术,旨在帮…

STM32+ESP8266局域网通信

若该文为原创文章,转载请注明原文出处。 此篇记录如何使用APP,在局域网内和STM32通信,并远程控制。 一、原理 STM32通过串口AT指令控制ESP8266(ESP-01S), ESP8266会产生一个热点,并创建了TCP服务,手机连接热点&#…

Solidity 开发环境

Solidity 开发环境 Solidity编辑器:Solidity编辑器是⼀种专⻔⽤于编写和编辑Solidity代码的编辑器。常⽤的Solidity编辑器包括 Visual Studio Code、Atom和Sublime Text。以太坊开发环境:以太坊开发环境(Ethereum Development Environment&a…

PHP Zip 文件处理指南

PHP Zip 文件处理指南 引言 ZIP文件是一种流行的压缩格式,广泛用于文件压缩和归档。PHP作为一门强大的服务器端脚本语言,提供了丰富的类和方法来处理ZIP文件。本文将详细介绍PHP中ZIP文件的处理方法,包括创建、添加文件、提取文件以及压缩和…

双重因子认证:守护数字安全的“双保险”

在数字化时代,网络安全问题日益凸显,个人信息和资产面临着诸多威胁。为了应对这些挑战,双重因子认证(Two-Factor Authentication,简称2FA)应运而生,成为一种重要的安全防护手段。 一、什么是双…

使用 Conda 管理 R 语言

1. 引言 在数据科学和统计分析领域,R 语言 是一门强大的编程语言,广泛用于数据可视化、统计建模和机器学习等任务。与 Python 相似,R 也有大量的开源包和库支持其功能。然而,在项目中使用 R 语言时,如何高效管理其环境…