坐标系与参考点：全局坐标系与局部坐标系
做飞行器设计这么多年，我见过太多新手在坐标系上栽跟头。说实话，坐标系选错了，后面所有计算都是白搭。今天咱们就好好聊聊这个看似基础、实则关键的话题。
一、全局坐标系与局部坐标系
先说说全局坐标系。你可以把它想象成「上帝视角」——一个固定在地球上的坐标系。我习惯把原点放在发射场或者机场跑道起点，Z轴指向天，X轴指向北或者指向飞行方向。这个坐标系不会动，所有飞行器的绝对位置都用它来描述。
局部坐标系就不一样了。它跟着飞行器走，原点通常在飞行器的质心上。X轴指向机头，Y轴指向右翼，Z轴指向下——这就是经典的机体坐标系。你想想看，飞行员说的「向左滚转」，其实就是绕局部坐标系的X轴转动。
关键区别：
- 全局坐标系：固定不动，描述绝对位置
- 局部坐标系：随飞行器运动，描述相对姿态
我在项目中遇到过一件事。有个同事把惯性张量算错了，查了两天发现——他把局部坐标系下的数据直接当成全局坐标系用了。嗯，这种错误其实很常见。
二、参考点的选择原则
参考点选在哪？我个人习惯遵循三条原则：
- 物理意义明确：首选质心、几何中心或关键结构点
- 计算方便：尽量让坐标轴与对称轴重合
- 便于测量：实际装配时能轻松找到这个点
举个例子。我做某型无人机时，把参考点选在了机头顶点。为什么？因为机头顶点在CAD模型里好定位，装配时也好测量。但后来发现质心计算要反复转换坐标，麻烦得很。后来我学乖了，直接以质心为参考点。
我的建议：如果条件允许，直接把参考点设在飞行器质心上。这样惯性张量计算最简洁，后续动力学分析也省事。
三、坐标变换基础
坐标变换说白了就是「换个角度看问题」。从全局坐标系到局部坐标系，或者反过来，都需要做变换。
最基本的变换是平移加旋转。平移好理解，就是原点挪个位置。旋转呢？用旋转矩阵。
import numpy as np
# 绕Z轴旋转角度theta
def rotation_z(theta):
c = np.cos(theta)
s = np.sin(theta)
return np.array([
[c, -s, 0],
[s, c, 0],
[0, 0, 1]
])
# 从全局到局部：先平移，再旋转
def global_to_local(pos_global, origin_local, theta):
# 平移：减去局部原点
translated = pos_global - origin_local
# 旋转：用旋转矩阵
R = rotation_z(theta)
return R @ translated
这段代码看着简单，但实际用起来坑不少。我曾经因为旋转顺序搞反了，算出来的姿态角差了180度。你想想看，飞机头朝东算成朝西，这要是真飞起来还得了？
注意：旋转矩阵的乘法不满足交换律。先绕X轴转再绕Y轴转，和先绕Y轴转再绕X轴转，结果完全不同。一定要搞清楚你的旋转顺序！
四、知识体系总览
下面这张图把本章的核心逻辑串起来了。我画图时习惯用这种结构化的方式，一目了然。
五、实际应用中的坑
讲几个我踩过的坑，希望对你有帮助。
坑一：坐标系方向搞反。 有一次我算某型导弹的惯性张量，结果怎么都对不上实验数据。查了三天，发现是Y轴方向定义反了。我习惯Y轴朝右，但那个项目的标准是Y轴朝左。就这么一个符号，所有交叉项全错了。
坑二：参考点漂移。 飞行器在飞行中燃料消耗，质心会移动。如果你一直用初始质心做参考点，那惯性张量就得实时更新。我曾经偷懒没更新，结果仿真结果跟实际飞行差了10%。
坑三：单位搞混。 这个说起来有点丢人。有次我把毫米当成米用了，算出来的惯性矩大了10^12倍。还好在评审阶段被同事发现了，不然图纸都发出去了。
避坑指南：
- 开工前先确认坐标系定义，最好写在文档里
- 每个坐标变换都做一次正向和反向验证
- 单位统一用国际单位制（米、千克、秒）
- 关键数据用两个独立方法交叉验证
六、小结
坐标系和参考点，说白了就是给飞行器建立一个「身份系统」。全局坐标系告诉别人你在哪，局部坐标系告诉别人你朝哪，参考点告诉别人从哪开始算。这三者配合好了，后面的质心和惯性张量计算才能顺风顺水。
我个人习惯在项目一开始就把坐标系定义写进技术规范里，每个参与的人都要签字确认。别嫌麻烦，这步省了，后面可能要花十倍的时间来填坑。
公众号:蓝海资料掘金营，微信deep3321