卫星姿轨控中的运动学与动力学-洪萨配资

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

在航天器控制系统中，姿态轨道控制（简称“姿轨控”）是确保卫星在轨稳定运行、精确指向和轨道维持的关键技术。其理论基础主要由运动学（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−ωx0−ωzωy−ωyωz0−ωx−ωz−ωyωx0

若使用欧拉角（如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}}ω=1000cosϕ−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=−Kpeq−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 区别总结

维度	运动学	动力学
研究目标	描述“如何运动”	解释“为何运动”
是否含物理参数	否（仅几何/时间关系）	是（质量、惯量、力、力矩）
方程类型	一阶微分（如q˙=f(ω)\dot{\mathbf{q}} = f(\boldsymbol{\omega})q˙=f(ω)）	二阶微分（如Jq¨=τ\mathbf{J} \ddot{\mathbf{q}} = \boldsymbol{\tau}Jq¨=τ）
控制设计用途	状态预测、传感器融合	控制律设计、执行机构建模