坐标系与坐标变换：飞行器建模的“语言”
做飞行器动力学仿真，说白了就是跟各种坐标系打交道。我刚开始接触这个领域时，总觉得坐标系这东西太抽象了——不就是个参考系嘛，能有多大区别？直到有一次，我把地心惯性系和地固系搞混了，算出来的轨道参数差了十几公里。嗯，从那以后我再也不敢小看坐标系了。
今天咱们就来聊聊飞行器建模中最核心的四个坐标系，以及它们之间的变换方法。你想想看，没有统一的“语言”，你怎么描述飞行器在空中的姿态和位置？
1. 地心惯性坐标系（ECI）
地心惯性坐标系，简称ECI。它是我们描述飞行器轨道运动的基础。
定义：
- 原点：地球质心
- X轴：指向春分点方向（J2000.0历元）
- Z轴：指向地球自转轴（北极方向）
- Y轴：由右手定则确定
这个坐标系最大的特点就是——它不跟着地球转。说白了，它是一个“静止”的参考系，适合用来描述卫星、导弹等飞行器的惯性运动。
核心要点：ECI坐标系是惯性系，牛顿定律在这里直接成立。你算轨道、算引力，都得用它。
我个人习惯在仿真开始时，把所有初始条件都转换到ECI系下。这样后续的积分计算会干净很多，不用操心地球自转带来的科里奥利力。
2. 地心地固坐标系（ECEF）
地心地固坐标系，简称ECEF。这个坐标系是跟着地球一起转的。
定义：
- 原点：地球质心
- X轴：指向本初子午线与赤道的交点
- Z轴：指向地球自转轴
- Y轴：由右手定则确定
ECEF和ECI最大的区别是什么？ECEF随地球自转，而ECI不转。这就意味着，同一个空间点，在两个坐标系下的坐标是不一样的，而且这个差异会随时间变化。
实用技巧：GPS输出的位置信息就是ECEF坐标。如果你要做地面站跟踪，用ECEF最方便。但如果你要算轨道动力学，记得转回ECI。
我曾经在做一个卫星仿真项目时，直接用GPS的ECEF坐标去算轨道积分，结果发现卫星“飞”得越来越偏。后来才意识到——ECEF下的速度包含了地球自转分量，直接积分会出大问题。
3. 机体坐标系（Body Frame）
机体坐标系是固定在飞行器上的。它随飞行器一起运动、一起旋转。
定义：
- 原点：飞行器质心
- X轴：指向机头方向（通常沿机体纵轴）
- Y轴：指向右翼方向
- Z轴：由右手定则确定（通常指向机腹）
这个坐标系用来描述飞行器自身的姿态、角速度、以及作用在机体上的力和力矩。说白了，传感器（IMU、陀螺仪）测量的数据，默认都是在机体坐标系下的。
注意：不同飞行器厂商对机体坐标系的定义可能略有差异。有的把Z轴指向上方，有的指向下方。做联合仿真时一定要确认清楚，否则姿态解算会反着来。
4. 欧拉角与四元数
好了，坐标系定义清楚了，接下来就是怎么在它们之间来回切换。这里有两个主流工具：欧拉角和四元数。
4.1 欧拉角
欧拉角用三个角度来描述旋转：偏航角（ψ）、俯仰角（θ）、滚转角（φ）。
旋转顺序：通常采用Z-Y-X顺序（偏航→俯仰→滚转）
// 欧拉角转旋转矩阵（Z-Y-X顺序）
R = Rz(ψ) * Ry(θ) * Rx(φ)
其中：
Rz(ψ) = [cosψ -sinψ 0; sinψ cosψ 0; 0 0 1]
Ry(θ) = [cosθ 0 sinθ; 0 1 0; -sinθ 0 cosθ]
Rx(φ) = [1 0 0; 0 cosφ -sinφ; 0 sinφ cosφ]
欧拉角直观易懂，但有个致命问题——万向锁。当俯仰角接近±90°时，偏航和滚转的旋转轴会重合，导致自由度丢失。
避坑指南：我曾经在无人机仿真中只用欧拉角，结果做大机动飞行时姿态解算直接炸了。后来换成四元数，问题迎刃而解。
4.2 四元数
四元数用四个参数表示旋转：q = [q0, q1, q2, q3]ᵀ，其中q0是标量部分，[q1, q2, q3]是矢量部分。
优势：
- 无万向锁问题
- 插值平滑
- 计算效率高（避免了三角函数）
// 四元数转旋转矩阵
R = [
1-2(q2²+q3²) 2(q1q2-q0q3) 2(q1q3+q0q2)
2(q1q2+q0q3) 1-2(q1²+q3²) 2(q2q3-q0q1)
2(q1q3-q0q2) 2(q2q3+q0q1) 1-2(q1²+q2²)
]
我的建议：做飞行器动力学仿真，姿态更新用四元数，人机交互显示用欧拉角。两者结合，既保证了计算稳定性，又方便调试观察。
5. 坐标系变换流程
在实际工程中，我们经常需要在ECI、ECEF、机体坐标系之间来回切换。下面这张图展示了典型的变换流程：
从图中可以看出，最常用的路径是：ECI → ECEF → 机体。ECI到ECEF的变换由地球自转决定（格林尼治恒星时角GMST），ECEF到机体的变换则由飞行器的姿态决定。
6. 工程实践中的坐标系选择
| 应用场景 | 推荐坐标系 | 原因 |
|---|---|---|
| 轨道动力学积分 | ECI | 惯性系，牛顿定律直接适用 |
| 地面站跟踪 | ECEF | 与地球固定，方便计算方位角、仰角 |
| 姿态控制 | 机体坐标系 | 传感器测量值在此系下 |
| 导航解算 | ECEF + 机体 | 需要来回变换 |
重要提醒：坐标系变换时，一定要注意变换矩阵的正交性。数值误差累积会导致矩阵不再正交，进而引起姿态发散。我建议每1000步做一次正交化修正。
好了，坐标系和坐标变换这部分就聊到这儿。记住一句话：坐标系选对了，仿真就成功了一半。下次咱们接着聊飞行器动力学方程的具体构建方法。
公众号:蓝海资料掘金营，微信deep3321