企业官网建设流程全解析

卫星姿轨控中的运动学与动力学

在航天器控制系统中，姿态轨道控制（简称“姿轨控”）是确保卫星在轨稳定运行、精确指向和轨道维持的关键技术。其理论基础主要由运动学（Kinematics）与动力学（Dynamics）两大部分构成。虽然两者紧密关联，但在建模目的、研究对象和数学表达上存在显著差异。本文将系统介绍姿轨控中运动学与动力学的基本概念、核心公式、研究对象、实际应用场景，并通过具体例子阐明二者的区别与联系。

一、基本概念

1.1 运动学（Kinematics）

运动学描述的是物体运动的几何特性，不涉及引起运动的力或力矩。在卫星姿轨控中，运动学主要关注：

• 姿态角随时间的变化（如欧拉角、四元数）
• 角速度与姿态之间的微分关系
• 位置与速度随时间的演化（轨道运动学）

其本质是“如何运动”，而非“为何如此运动”。

1.2 动力学（Dynamics）

动力学研究的是力（或力矩）与运动之间的因果关系，即由外力/力矩引起的状态变化。在姿轨控中，动力学描述：

• 角动量变化与外力矩的关系（欧拉方程）
• 轨道摄动与外力（如重力、大气阻力、太阳光压）的关系

其核心问题是“为何这样运动”，依赖牛顿力学或拉格朗日/哈密顿力学框架。

二、研究对象

三、核心公式

3.1 姿态运动学

常用四元数q=[q0,q1,q2,q3]T\mathbf{q} = [q_0, q_1, q_2, q_3]^Tq=[q0​,q1​,q2​,q3​]T表示卫星本体坐标系相对于惯性系的旋转。四元数与角速度ω=[ωx,ωy,ωz]T\boldsymbol{\omega} = [\omega_x, \omega_y, \omega_z]^Tω=[ωx​,ωy​,ωz​]T的关系为：

q˙=12Ω(ω)q \dot{\mathbf{q}} = \frac{1}{2} \mathbf{\Omega}(\boldsymbol{\omega}) \mathbf{q}q˙​=21​Ω(ω)q

其中，

Ω(ω)=[0−ωx−ωy−ωzωx0ωz−ωyωy−ωz0ωxωzωy−ωx0] \mathbf{\Omega}(\boldsymbol{\omega}) = \begin{bmatrix} 0 & -\omega_x & -\omega_y & -\omega_z \\ \omega_x & 0 & \omega_z & -\omega_y \\ \omega_y & -\omega_z & 0 & \omega_x \\ \omega_z & \omega_y & -\omega_x & 0 \end{bmatrix}Ω(ω)=​0ωx​ωy​ωz​​−ωx​0−ωz​ωy​​−ωy​ωz​0−ωx​​−ωz​−ωy​ωx​0​​

若使用欧拉角（如3-2-1滚动-俯仰-偏航顺序），角速度与欧拉角速率θ˙=[ϕ˙,θ˙,ψ˙]T\dot{\boldsymbol{\theta}} = [\dot{\phi}, \dot{\theta}, \dot{\psi}]^Tθ˙=[ϕ˙​,θ˙,ψ˙​]T的关系为：

ω=[10−sin⁡θ0cos⁡ϕsin⁡ϕcos⁡θ0−sin⁡ϕcos⁡ϕcos⁡θ]θ˙ \boldsymbol{\omega} = \begin{bmatrix} 1 & 0 & -\sin\theta \\ 0 & \cos\phi & \sin\phi \cos\theta \\ 0 & -\sin\phi & \cos\phi \cos\theta \end{bmatrix} \dot{\boldsymbol{\theta}}ω=​100​0cosϕ−sinϕ​−sinθsinϕcosθcosϕcosθ​​θ˙

3.2 姿态动力学（刚体欧拉方程）

对于刚体卫星，忽略柔性模态，其角动量H=Jω\mathbf{H} = \mathbf{J} \boldsymbol{\omega}H=Jω，其中J\mathbf{J}J为惯量张量。由角动量定理H˙=τ\dot{\mathbf{H}} = \boldsymbol{\tau}H˙=τ，可得：

Jω˙+ω×(Jω)=τ \mathbf{J} \dot{\boldsymbol{\omega}} + \boldsymbol{\omega} \times (\mathbf{J} \boldsymbol{\omega}) = \boldsymbol{\tau}Jω˙+ω×(Jω)=τ

其中τ\boldsymbol{\tau}τ为作用在卫星上的总力矩（包括控制力矩τc\boldsymbol{\tau}_cτc​和干扰力矩τd\boldsymbol{\tau}_dτd​）。

该方程是非线性的，体现了陀螺耦合效应（gyroscopic coupling）。

3.3 轨道运动学

在惯性系中，卫星位置r\mathbf{r}r与速度v\mathbf{v}v满足：

r˙=v,v˙=a \dot{\mathbf{r}} = \mathbf{v}, \quad \dot{\mathbf{v}} = \mathbf{a}r˙=v,v˙=a

其中a\mathbf{a}a为加速度——这本身是运动学关系。但a\mathbf{a}a的来源由动力学决定。

3.4 轨道动力学（二体问题）

在理想二体引力场中（忽略摄动），动力学方程为：

r¨=−μr3r \ddot{\mathbf{r}} = -\frac{\mu}{r^3} \mathbf{r}r¨=−r3μ​r

其中μ=GM\mu = G Mμ=GM为中心天体引力常数，r=∥r∥r = \|\mathbf{r}\|r=∥r∥。

若考虑摄动力ap\mathbf{a}_pap​（如大气阻力、太阳光压、地球非球形引力），则：

r¨=−μr3r+ap+Fthrustm \ddot{\mathbf{r}} = -\frac{\mu}{r^3} \mathbf{r} + \mathbf{a}_p + \frac{\mathbf{F}_{\text{thrust}}}{m}r¨=−r3μ​r+ap​+mFthrust​​

四、实际应用举例

4.1 姿态运动学应用：星敏感器数据融合

星敏感器输出的是以四元数或方向余弦矩阵表示的姿态。控制系统需通过运动学方程将角速度陀螺仪（测量ω\boldsymbol{\omega}ω）的数据与星敏姿态进行融合（如使用扩展卡尔曼滤波 EKF）。此处仅需运动学模型即可预测姿态演化：

q^k+1=qk⊗(cos⁡∥ω∥Δt2,ω∥ω∥sin⁡∥ω∥Δt2) \hat{\mathbf{q}}_{k+1} = \mathbf{q}_{k} \otimes \left( \cos\frac{\|\boldsymbol{\omega}\| \Delta t}{2}, \frac{\boldsymbol{\omega}}{\|\boldsymbol{\omega}\|} \sin\frac{\|\boldsymbol{\omega}\| \Delta t}{2} \right)q^​k+1​=qk​⊗(cos2∥ω∥Δt​,∥ω∥ω​sin2∥ω∥Δt​)

无需求解动力学方程。

4.2 姿态动力学应用：飞轮控制律设计

当设计反作用飞轮（Reaction Wheel）控制律时，必须考虑卫星本体的转动惯量和陀螺耦合。例如，采用 PD 控制：

τc=−Kpeq−Kdω \boldsymbol{\tau}_c = -K_p \boldsymbol{e}_q - K_d \boldsymbol{\omega}τc​=−Kp​eq​−Kd​ω

但执行该力矩时，需通过动力学方程反推飞轮所需的角加速度：

Jwω˙w=−τc \mathbf{J}_w \dot{\boldsymbol{\omega}}_w = -\boldsymbol{\tau}_cJw​ω˙w​=−τc​

其中Jw\mathbf{J}_wJw​为飞轮惯量。若忽略动力学，仅用运动学设计控制器，将导致实际控制性能下降甚至不稳定。

4.3 轨道运动学 vs 动力学：GPS 定轨

• 运动学定轨：仅利用 GPS 接收机测量的位置/速度数据，通过插值或平滑得到轨道（不依赖物理模型），适用于高精度事后处理。
• 动力学定轨：结合动力学方程与观测数据，通过最小二乘或滤波方法估计轨道状态和摄动参数，适用于实时导航和长期预报。

五、运动学与动力学的区别与联系

5.1 区别总结

5.2 联系与协同

在完整的姿轨控系统中，二者缺一不可：

• 运动学提供状态转移模型：用于状态估计（如 EKF 中的预测步）。
• 动力学提供控制输入-输出映射：用于设计鲁棒控制器。

例如，在模型预测控制（MPC）中，预测模型通常同时包含运动学与动力学方程，以高精度预测未来状态并优化控制输入。

更进一步，二者可通过微分几何或李群/李代数统一描述（如SO(3)SO(3)SO(3)上的姿态运动），但工程实现中仍常分开建模以简化计算。

六、结语

卫星姿轨控中的运动学与动力学，如同“描述语言”与“物理法则”的关系。运动学关注姿态与轨道的几何演化规律，而动力学揭示其背后的力学因果机制。在高精度任务（如对地观测、深空探测、编队飞行）中，正确建模并耦合二者，是实现亚角秒级指向精度和厘米级轨道控制的前提。

理解二者的区别与联系，不仅有助于构建更准确的仿真模型，也为先进控制算法（如自适应控制、滑模控制、非线性观测器）的设计奠定理论基础。