坐标系基础：常用坐标系与转换方法
各位同学，今天我们来聊聊导弹弹道建模里最基础、也最容易让人头晕的东西——坐标系。说实话，我刚开始做弹道仿真那会儿，也在这上面栽过跟头。坐标系选错了，后面算出来的弹道数据全是错的，还找不出原因。所以这一节，咱们把它彻底搞明白。
为什么需要多个坐标系？
你想想看，导弹在空中飞，涉及到的运动描述其实很复杂。它要绕地球转，要相对地面飞，还要自己翻滚转弯。用一个坐标系去描述所有运动，那是不可能的。就像你描述一个人走路，用地面做参考最方便；但描述他手里的手机怎么转，就得用手机本身的坐标系了。
所以，我们至少需要三个层次的坐标系：
- 描述导弹绕地球飞——用地心惯性系
- 描述导弹相对发射点的运动——用发射坐标系
- 描述导弹自身的姿态——用弹体坐标系
这三个坐标系之间，需要互相转换。嗯，这就是本节的核心内容。
地心惯性系（ECI）
地心惯性系，说白了就是一个不跟着地球转的坐标系。它的原点在地心，Z轴指向北极，X轴指向春分点（或者某个固定方向），Y轴按右手定则补齐。
为什么要用这个坐标系？因为牛顿定律只适用于惯性系。导弹飞行的动力学方程，必须建立在惯性系里才能用。我在项目中遇到过有人直接用发射坐标系写动力学方程，结果算出来的加速度全是错的——因为发射坐标系本身在跟着地球转，有惯性力。
关键点：地心惯性系是“绝对”的参考系，所有动力学方程的第一步推导，都要回到这里。
发射坐标系
发射坐标系就接地气多了。它的原点在导弹发射点，X轴指向发射方向（通常是正北或目标方向），Y轴垂直向上，Z轴按右手定则指向侧向。
这个坐标系的好处是直观。你站在发射场，看着导弹飞出去，它的位置、速度、高度，用发射坐标系描述最方便。我记得有一次做靶场测试，指挥员问“导弹现在多高”，我直接报发射坐标系里的Y坐标，大家都听得懂。
但要注意，发射坐标系是跟着地球自转的。所以它其实是一个非惯性系。如果你要在发射坐标系里写动力学方程，必须加上哥氏力和离心力项。这个坑，我踩过。
避坑指南：我曾经在发射坐标系里直接套用牛顿第二定律，忘了加惯性力项，结果弹道仿真偏差了十几公里。后来花了三天才找到原因。记住：发射坐标系不是惯性系！
弹体坐标系
弹体坐标系是固定在导弹身上的。原点在导弹质心，X轴沿弹体纵轴指向头部，Y轴在弹体对称平面内垂直X轴向上，Z轴按右手定则指向右侧。
这个坐标系用来描述导弹的姿态——俯仰、偏航、滚转。说白了，导弹怎么转，全看弹体坐标系相对发射坐标系（或惯性系）的夹角。
我建议初学者先把弹体坐标系和发射坐标系的关系搞清楚。因为导弹的控制指令（舵面偏角、推力矢量）都是在弹体坐标系里给出的，而导航数据（位置、速度）是在发射坐标系里的。这两个坐标系之间的转换，是弹道仿真的核心环节。
坐标系转换方法
坐标系转换，说白了就是旋转。把一个坐标系里的向量，通过旋转矩阵，变到另一个坐标系里去。
常用的转换有三种：
- 绕X轴旋转——滚转
- 绕Y轴旋转——偏航
- 绕Z轴旋转——俯仰
每次旋转对应一个旋转矩阵。三个矩阵乘在一起，就是完整的转换矩阵。
举个例子，从发射坐标系到弹体坐标系的转换，通常按3-2-1顺序：先绕Z轴转俯仰角，再绕新Y轴转偏航角，最后绕新X轴转滚转角。这个顺序不能乱，乱了角度定义就全错了。
个人经验：我习惯把旋转矩阵写成函数库，每次调用时传入三个角度，返回一个3x3矩阵。这样代码干净，不容易出错。千万别每次都手算矩阵元素，太容易写反正负号了。
转换矩阵的数学形式
绕Z轴旋转角度ψ（俯仰）：
Rz(ψ) = [cosψ -sinψ 0]
[sinψ cosψ 0]
[0 0 1]绕Y轴旋转角度θ（偏航）：
Ry(θ) = [cosθ 0 sinθ]
[0 1 0 ]
[-sinθ 0 cosθ]绕X轴旋转角度φ（滚转）：
Rx(φ) = [1 0 0 ]
[0 cosφ -sinφ]
[0 sinφ cosφ]完整的转换矩阵：
C = Rx(φ) * Ry(θ) * Rz(ψ)注意乘法顺序。矩阵乘法不满足交换律，顺序错了，结果就错了。我见过有人把顺序写反，仿真出来的导弹姿态一直在乱抖，查了两天没找到原因——其实就是矩阵乘反了。
知识体系总览
下面这张图，是我自己总结的坐标系知识框架。你看一眼，心里就有数了。
实际应用中的注意事项
坐标系转换看起来就是几个矩阵乘法，但实际用起来，有几个细节特别容易出错：
- 角度定义要统一：俯仰角是抬头为正还是低头为正？偏航角是左转为正还是右转为正？不同教材定义不一样。我建议你在代码开头用注释写清楚，免得后面自己都忘了。
- 注意奇点：当俯仰角接近±90度时，偏航角和滚转角会退化（万向锁问题）。如果你做的是垂直发射的导弹，这个情况很常见。我一般用四元数代替欧拉角来避免这个问题。
- 单位要一致：角度用弧度还是度？我习惯在内部计算全部用弧度，只在输入输出时转换。这样不容易乱。
总结一下：坐标系是弹道仿真的地基。地基没打好，上面盖的楼再漂亮也是危楼。我建议你花时间把这三个坐标系的关系画清楚，把转换矩阵推导一遍，再写代码验证。这一步做扎实了，后面的动力学、制导律都会顺利很多。
好了，坐标系的内容就讲到这里。下一节我们进入动力学方程，到时候你会看到，坐标系转换的功底有多重要。
公众号:蓝海资料掘金营，微信deep3321