3. 坐标系与运动变量
在无人机气动参数辨识与标定中，坐标系不是数学装饰，而是所有测量、建模与参数估计的“共同语言”。同一组加速度、角速率和空速数据，若没有明确其表达坐标系，辨识结果可能完全错误。本章围绕三套核心坐标系与六个关键运动变量展开，建立从传感器到气动力模型的统一描述框架。
一、三套核心坐标系及其物理意义
无人机运动描述通常涉及三套右手直角坐标系，它们各自承担不同职责：
- 地面坐标系（惯性系，OeXeYeZe）：原点固定于地面某点，Xe指向北，Ye指向东，Ze垂直向下。用于描述航迹、位置、地速和重力方向，是导航与轨迹分析的基准。
- 机体坐标系（ObXbYbZb）：原点位于无人机重心，Xb指向机头，Yb指向右翼，Zb指向机腹。用于描述机体角速率、加速度计和陀螺仪的安装轴向，是飞控与惯导的直接表达系。
- 气流坐标系（OaXaYaZa）：原点同样位于重心，Xa沿空速矢量方向，Za在机体对称面内垂直于Xa向下，Ya由右手定则确定。用于定义升力、阻力和侧力，是气动参数辨识的核心表达系。
三套坐标系之间通过欧拉角与气流角相互转换。地面系到机体系由偏航角ψ、俯仰角θ、滚转角φ三个欧拉角描述；机体系到气流系则由迎角α和侧滑角β描述。理解这一转换链，是后续建立六自由度运动方程和观测方程的前提。
二、欧拉角与角速率的运动学关系
欧拉角描述的是机体系相对地面系的姿态，而角速率p、q、r是机体系下的旋转速度。两者并非简单导数关系，而是通过欧拉角运动学方程耦合：
φ̇ = p + q·sinφ·tanθ + r·cosφ·tanθ
θ̇ = q·cosφ − r·sinφ
ψ̇ = q·sinφ/cosθ + r·cosφ/cosθ当俯仰角θ接近±90°时，tanθ和1/cosθ发散，出现万向节死锁。这也是为什么在气动辨识中，常直接使用四元数或旋转矩阵进行姿态更新，而仅在输出和解释时转换为欧拉角。对于大多数固定翼和小型多旋翼，θ通常远离±90°，欧拉角仍可安全使用。
标定提示：陀螺仪输出的p、q、r是机体系角速率，而飞控姿态解算输出的φ、θ、ψ是欧拉角。在参数辨识中，若用欧拉角差分求角速率，必须经过上述运动学方程转换，否则在滚转或偏航机动中会引入显著误差。
三、迎角与侧滑角的几何定义
迎角α和侧滑角β是气流坐标系相对机体坐标系的两个角度，直接决定气动力和力矩：
- 迎角α：空速矢量在机体对称面（XbZb平面）内的投影与机体Xb轴的夹角。α增大通常使升力增大，但超过临界迎角后发生失速。
- 侧滑角β：空速矢量与机体对称面的夹角。β不为零时产生侧力和偏航力矩，是协调转弯和风标稳定性的关键变量。
由机体系速度分量(u, v, w)可计算：
α = arctan(w / u)
β = arcsin(v / V)
V = sqrt(u² + v² + w²)在气动参数辨识中，α和β通常由空速管、迎角传感器或惯性/气压组合估计得到。若直接使用GPS地速代替空速，在风场较大时α和β估计会严重偏离真实值，导致气动导数辨识失真。
四、坐标系选择与变量对照清单
下表汇总了三套坐标系下常用变量的表达方式，便于在建模和数据处理时快速对照：
| 变量类型 | 地面坐标系 | 机体坐标系 | 气流坐标系 |
|---|---|---|---|
| 速度 | Ve = (ẋe, ẏe, że) | Vb = (u, v, w) | Va = (V, 0, 0) |
| 角速率 | ψ̇, θ̇, φ̇ | p, q, r | 不直接定义 |
| 力 | 重力、推力分量 | 机体轴力 X, Y, Z | 升力 L、阻力 D、侧力 Y |
| 力矩 | 不直接使用 | 滚转 L、俯仰 M、偏航 N | 不直接使用 |
| 角度 | 航迹倾角、航向角 | 安装角、舵面偏角 | 迎角 α、侧滑角 β |
五、典型误区与复盘
误区一：混用坐标系不声明。将机体角速率p、q、r直接当作欧拉角导数使用，在滚转机动中误差可达数倍。复盘时应检查所有方程是否明确标注了参考系。
误区二：用GPS地速计算迎角。地速包含风速，而迎角定义基于空速。有风时两者差异显著，尤其在低速盘旋和起降阶段。正确做法是使用空速管或大气数据系统。
误区三：忽略安装角与重心偏移。传感器安装位置与重心不重合时，加速度计会敏感到角加速度引起的附加分量。标定前应测量安装矩阵和杆臂，或在辨识中将其作为待估参数。
误区二：用GPS地速计算迎角。地速包含风速，而迎角定义基于空速。有风时两者差异显著，尤其在低速盘旋和起降阶段。正确做法是使用空速管或大气数据系统。
误区三：忽略安装角与重心偏移。传感器安装位置与重心不重合时，加速度计会敏感到角加速度引起的附加分量。标定前应测量安装矩阵和杆臂，或在辨识中将其作为待估参数。
掌握三套坐标系与六个运动变量的定义、转换和物理意义，是后续建立六自由度模型、设计激励机动和解释辨识结果的基础。在实际工程中，建议在数据预处理阶段就统一坐标系约定，并在代码和文档中显式标注每个变量的参考系，避免因“隐性约定”导致返工。