坐标系与坐标变换：从地面到桨毂的数学桥梁

各位同学，今天我们来聊聊倾转旋翼动力学里最基础、也最容易出错的一环——坐标系与坐标变换。说实话，我当年刚入行时，就因为坐标系搞混，在仿真里跑出了一个“直升机倒着飞”的诡异结果。嗯，从那以后，我对坐标变换就格外上心。

你想想看，一架倾转旋翼机，既有固定翼的飞行模式，又有直升机的悬停能力。它的运动描述，需要在多个坐标系之间来回切换。地面坐标系、机体坐标系、桨毂坐标系……每个坐标系都有自己的“视角”，而我们要做的，就是搭建它们之间的数学桥梁。

1. 地面坐标系（Earth Frame）

地面坐标系，也叫惯性坐标系（近似）。我习惯把它定义为：原点在地面某点，X轴指向北（或东，看你习惯），Z轴垂直向下指向地心，Y轴按右手定则确定。

为什么Z轴向下？因为飞行器设计里，高度通常用负的Z表示。这在导航里很常见。我个人建议初学者先别纠结方向，关键是记住：地面坐标系是“绝对参考系”，其他所有坐标变换最终都要回到这里。

2. 机体坐标系（Body Frame）

机体坐标系固连在飞行器上。原点在重心，X轴指向机头，Y轴指向右翼，Z轴指向下（符合右手定则）。

这里有个坑：不同厂家定义的机体坐标系可能不同。比如波音和空客的X轴方向可能相反。我在做联合仿真时，就吃过这个亏——两个团队的模型对接，结果力矩方向全反了。所以，一定要在文档里明确写出你的坐标系定义。

机体坐标系最大的用处是：描述飞行器自身的角速度和线速度。比如滚转角速度p、俯仰角速度q、偏航角速度r，都是在机体坐标系下测量的。

3. 桨毂坐标系（Hub Frame）

这个坐标系是倾转旋翼特有的。它固连在旋翼桨毂上，原点在桨毂中心。Z轴沿旋翼轴方向，X轴和Y轴在旋转平面内。

为什么要单独搞一个桨毂坐标系？因为旋翼的气动力计算，需要在旋转坐标系下进行。你想想看，桨叶在旋转，它的攻角、来流速度都是随时间变化的。在桨毂坐标系下，这些计算会简单很多。

我记得有一次，学生问我：“老师，能不能直接用机体坐标系算旋翼力？”我说：“可以，但你会被三角函数折磨死。”说白了，选择合适的坐标系，就是选择计算的便利性。

4. 欧拉角与四元数

这是坐标变换的核心工具。我们先说欧拉角。

4.1 欧拉角

欧拉角用三个角度描述刚体姿态：滚转角φ、俯仰角θ、偏航角ψ。旋转顺序我习惯用Z-Y-X（偏航→俯仰→滚转），这是航空航天领域的标准。

但欧拉角有个致命问题：万向锁。当俯仰角接近±90°时，滚转和偏航会耦合，导致自由度丢失。我在做倾转旋翼的过渡飞行仿真时，就遇到过这个问题——旋翼从垂直转到水平的过程中，俯仰角刚好经过90°，结果仿真直接发散。

4.2 四元数

四元数就是来解决这个问题的。它用四个参数表示旋转：一个标量部分和三个矢量部分。没有奇点，计算效率高，适合计算机实现。

四元数的形式是：

q = [q0, q1, q2, q3]^T
其中 q0 = cos(θ/2)
     [q1, q2, q3] = sin(θ/2) * [ex, ey, ez]

这里θ是旋转角，[ex, ey, ez]是旋转轴单位矢量。我个人习惯用四元数做内部计算，最后再转成欧拉角输出给用户看。因为人脑对欧拉角更直观，但计算机喜欢四元数。

5. 坐标变换矩阵推导

好了，重头戏来了。我们来看看坐标变换矩阵怎么推。

5.1 从地面坐标系到机体坐标系

这个变换矩阵由三个基本旋转矩阵相乘得到。按Z-Y-X顺序：

C_ground_to_body = Rx(φ) * Ry(θ) * Rz(ψ)

其中：
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φ]

注意顺序：先转的矩阵在右边。这个顺序我经常搞混，后来干脆在代码里加注释：% 先偏航，再俯仰，最后滚转。

5.2 从机体坐标系到桨毂坐标系

这个变换相对简单。因为桨毂坐标系相对于机体坐标系，通常只有一个倾转角α（倾转旋翼的倾转角）。

C_body_to_hub = Ry(α) = [cosα  0  -sinα; 0  1  0; sinα  0  cosα]

当α=0时，旋翼轴垂直向上（直升机模式）；α=90°时，旋翼轴水平向前（固定翼模式）。

5.3 完整的变换链

从地面到桨毂的完整变换是：

C_ground_to_hub = C_body_to_hub * C_ground_to_body

反过来，从桨毂到地面：

C_hub_to_ground = (C_ground_to_hub)^T

因为旋转矩阵是正交矩阵，逆等于转置。这个性质在编程里特别有用——不用求逆，直接转置就行。

6. 知识体系总览

下面这张图，是我做课程时画的。它把整个坐标系与坐标变换的逻辑串起来了：

7. 实战建议

最后，给各位几个实战建议：

• 统一符号命名：我习惯用C_a_b表示从a系到b系的变换矩阵。比如C_g_b是地面到机体。这样在代码里一目了然。
• 验证变换矩阵：写完后，用单位矩阵测试。比如当所有角度为0时，变换矩阵应该是单位阵。我每次写完都会做这个测试。
• 注意数值精度：四元数长时间积分后，可能会失去归一化。记得每隔几步做一次归一化。我在Simulink里会用norm(q)检查。
• 画图辅助理解：在MATLAB里用plot3画出坐标系，旋转一下看看。视觉反馈比纯数学推导更直观。