4. 刚体运动学基础：位置、速度、加速度在惯性系下的表达

各位同学，欢迎来到这门课的第一章。我是你们这门课的主讲，一个在飞行器设计、数学建模和控制工程里摸爬滚打了十几年的老工程师。今天咱们要聊的，是所有飞行器动力学方程的基石——刚体运动学。

说白了，就是搞清楚一个硬邦邦的物体（比如飞机、导弹、火箭）在空间里是怎么“动”的。你想想看，要控制它，首先得知道它在哪、跑多快、加速度多大，对吧？

我个人习惯，在讲任何复杂公式之前，先把这个“坐标系”的概念理清楚。坐标系选不对，后面全是白搭。我在项目里见过太多人，就是因为惯性系和体坐标系搞混了，仿真结果飞到了外太空去。

4.1 为什么非得是“惯性系”？

咱们先解决第一个问题：什么是惯性系？

牛顿他老人家说过，力等于质量乘以加速度。这个公式成立的前提，就是你得在一个“不动的”或者“匀速直线运动”的参考系里看问题。这个参考系，就叫惯性参考系，简称惯性系。

在地球表面做飞行器设计，我们通常把“地心惯性坐标系”当作惯性系。它的原点在地心，Z轴指向北极，X轴指向春分点。嗯，听起来有点绕，但你只要记住：在这个坐标系里，牛顿定律是成立的。

我曾经接手过一个项目，新同事直接用“地球固连坐标系”（就是跟着地球自转的那个）来推导控制律。结果飞机在高速巡航时，侧向控制总是有静差。查了三天，才发现是忘了考虑地球自转带来的牵连加速度。嗯，从那以后，我要求团队所有人在推导基础方程时，必须从惯性系开始。

4.2 位置：你在哪？

在惯性系里描述一个刚体的位置，其实很简单。我们用一个三维向量 r 来表示。

假设惯性系是 O-XYZ，刚体的质心在点 P。那么位置向量就是：

r = [x, y, z]^T

这里的 x, y, z 就是 P 点在三个坐标轴上的投影。没什么花头，就是初中几何。

但这里有个坑：刚体不仅有位置，还有姿态。位置告诉你质心在哪，姿态告诉你它“脑袋”朝哪。不过这一章咱们先只聊平动，姿态（欧拉角、四元数）后面章节会专门讲。

符号	含义	单位
r	位置向量（惯性系下）	米 (m)
x, y, z	位置分量	米 (m)

4.3 速度：你跑多快？

位置对时间求导，就是速度。在惯性系下，速度向量 v 就是：

v = dr/dt = [dx/dt, dy/dt, dz/dt]^T

记作：v = [u, v, w]^T。注意，这里的 u, v, w 是惯性系下的速度分量，不是机体坐标系下的。

我个人习惯，在推导时会把速度写成这样：

v_I = [u_I, v_I, w_I]^T

加个下标 I，时刻提醒自己这是惯性系下的量。这个小习惯，能帮你避免后面90%的坐标系混淆错误。

4.4 加速度：你受力多大？

速度对时间再求一次导，就是加速度 a：

a = dv/dt = d²r/dt² = [d²x/dt², d²y/dt², d²z/dt²]^T

记作：a = [a_x, a_y, a_z]^T。

这个加速度，就是牛顿第二定律里的那个 a。你只要把作用在飞行器上的所有外力（推力、气动力、重力）加起来，除以质量，就能得到这个加速度。

你想想看，是不是很简单？

F = m * a
=> a = F / m

但这里有个关键点：这个 F 必须是惯性系下的合力。如果你在体坐标系下算出了气动力，必须通过坐标变换矩阵转到惯性系下，才能代入这个公式。

4.5 知识体系总览

说了这么多，咱们用一张图来总结一下这一章的核心逻辑。我特意画了张SVG图，帮你把位置、速度、加速度的关系理清楚。

这张图你看懂了吗？从上往下看，位置求导得速度，速度求导得加速度。而加速度通过牛顿第二定律，和外力直接挂钩。这就是整个飞行器动力学方程的起点。

4.6 一个小例子

咱们来走一个最简单的例子。假设一架飞机在惯性系下做直线飞行，位置随时间变化为：

x(t) = 100 * t
y(t) = 0
z(t) = 10000

那么它的速度就是：

u = dx/dt = 100 m/s
v = dy/dt = 0 m/s
w = dz/dt = 0 m/s

加速度呢？

a_x = du/dt = 0 m/s²
a_y = 0
a_z = 0

匀速直线运动，加速度为零。如果这时候飞机受到的推力正好等于气动阻力，重力被升力平衡，那这个运动状态就是对的。

嗯，这个例子虽然简单，但它揭示了运动学方程的本质：你只要知道了位置随时间的变化规律，速度和加速度就是纯粹的数学运算。真正的难点，在于怎么通过力和力矩来求出这个位置变化规律——那是动力学的事，咱们后面再聊。

好了，这一章的内容就到这里。记住我今天说的：坐标系是根，求导是本。把这两个搞扎实了，后面的欧拉角、四元数、动力学方程，你学起来会轻松很多。

---