第二讲：坐标系那些事儿
各位同学好，我是你们的老张。
今天咱们聊聊坐标系。说实话，我刚入行那会儿，最头疼的就是这玩意儿。地面坐标系、弹体坐标系、速度坐标系……三个坐标系来回倒腾，稍不留神就搞混了。有一次仿真结果怎么都对不上，查了三天，最后发现是坐标系转换矩阵写反了——嗯，从那以后我再也不敢小看坐标系了。
1. 地面坐标系：我们的“绝对参考”
地面坐标系，说白了就是站在地面上看导弹怎么飞。它是个惯性坐标系，不跟着导弹转。
定义很简单：
- 原点O：通常选在导弹发射点
- Ox轴：指向目标方向（水平面内）
- Oy轴：垂直向上（指向天）
- Oz轴：按右手定则确定
你想想看，这个坐标系就像我们平时看地图一样。东、北、天三个方向，清清楚楚。我在项目中经常用地面坐标系来记录导弹的飞行轨迹，因为它最直观——导弹飞了多远、多高、偏了多少，一目了然。
3. 弹体坐标系：跟着导弹一起转
弹体坐标系就更有意思了。它固定在导弹上，导弹怎么转，它就怎么转。
定义是这样的：
- 原点O：导弹质心
- Ox轴：沿弹体纵轴，指向弹头
- Oy轴：在弹体对称面内，垂直于Ox轴，指向上方
- Oz轴：按右手定则确定
为什么需要这个坐标系？因为导弹上的传感器（比如陀螺仪、加速度计）都是装在弹体上的。它们测量的数据，天然就是弹体坐标系下的。我做过一个项目，需要把弹体坐标系下的角速度转换成地面坐标系下的姿态角——那转换矩阵写得我头都大了。
3. 速度坐标系：顺着气流的方向
速度坐标系，也叫“气流坐标系”。它跟导弹飞行的方向有关。
定义：
- 原点O：导弹质心
- Ox轴：沿速度矢量方向
- Oy轴：在弹体对称面内，垂直于Ox轴，指向上方
- Oz轴：按右手定则确定
这个坐标系有什么用？说白了，导弹受到的空气动力（升力、阻力、侧向力）都是相对于气流方向来定义的。你想想看，导弹迎风飞行时，阻力方向永远和速度方向相反——这就是速度坐标系的意义。
我个人习惯在计算气动参数时，先把所有数据转到速度坐标系下处理，算完了再转回弹体坐标系。这样逻辑清晰，不容易出错。
4. 坐标系之间的转换：核心中的核心
好了，三个坐标系都认识了。现在的问题是：它们之间怎么互相转换？
转换的核心就是——旋转矩阵。说白了，就是把一个坐标系下的向量，通过三次旋转，转到另一个坐标系下。
4.1 地面坐标系 → 弹体坐标系
这个转换需要三个旋转：
- 绕Oy轴转偏航角ψ
- 绕新的Oz轴转俯仰角θ
- 绕新的Ox轴转滚转角γ
转换矩阵长这样：
C_g2b = [cosθ·cosψ, sinθ, -cosθ·sinψ;
-sinθ·cosψ·cosγ + sinψ·sinγ, cosθ·cosγ, sinθ·sinψ·cosγ + cosψ·sinγ;
sinθ·cosψ·sinγ + sinψ·cosγ, -cosθ·sinγ, -sinθ·sinψ·sinγ + cosψ·cosγ]看着复杂吧？别怕。我在实际项目中从来不手算这个，都是写个函数库，直接调用。但你要理解它的物理意义——每个元素代表一个方向余弦。
4.2 地面坐标系 → 速度坐标系
这个转换需要两个角度：
- 弹道偏角ψv：速度矢量在水平面内的投影与Ox轴的夹角
- 弹道倾角θv：速度矢量与水平面的夹角
转换矩阵：
C_g2v = [cosθ_v·cosψ_v, sinθ_v, -cosθ_v·sinψ_v;
-sinθ_v·cosψ_v, cosθ_v, sinθ_v·sinψ_v;
sinψ_v, 0, cosψ_v]4.3 弹体坐标系 → 速度坐标系
这个转换最关键，因为它涉及攻角α和侧滑角β：
- 攻角α：速度矢量在弹体对称面内的投影与弹体纵轴的夹角
- 侧滑角β：速度矢量与弹体对称面的夹角
转换矩阵：
C_b2v = [cosα·cosβ, sinα, -cosα·sinβ;
-sinα·cosβ, cosα, sinα·sinβ;
sinβ, 0, cosβ]5. 知识体系总览
说了这么多，咱们用一张图来总结一下：
6. 避坑指南
最后，分享几个我踩过的坑：
- 旋转顺序不能乱：我曾经在代码里把偏航和俯仰的顺序搞反了，结果导弹在天上画了个8字。记住：地面到弹体是偏航→俯仰→滚转，顺序不能变。
- 角度正负号要统一：不同教材对角度正负的定义可能不同。我建议统一用右手定则——绕轴逆时针为正。这样不容易乱。
- 数值精度问题：角度接近90度时，三角函数值变化很快。这时候用四元数代替欧拉角会更稳定。嗯，这个后面会专门讲。
好了，坐标系的内容就讲到这里。记住一句话：坐标系是制导控制的“语言”，语言不通，什么都白搭。下一讲咱们聊聊导弹飞行的动力学方程——那才是真正有意思的东西。