3. 相对运动动力学模型：C-W方程（Hill方程）的推导与假设
各位同学，今天我们来聊聊航天器交会对接里最基础、也最经典的一个模型——C-W方程，也叫Hill方程。说实话，我刚入行那会儿，第一次看到这个方程时，觉得它长得挺简单，不就是几个线性微分方程嘛。后来在实际项目中吃过亏，才明白它背后的假设有多讲究。
好，咱们一步步来。先搞清楚：这个方程到底在描述什么？
3.1 为什么要用相对运动模型？
你想啊，两个航天器在太空中飞行，如果都用绝对轨道去描述，那计算量太大了。而且交会对接时，我们关心的不是它们各自在哪儿，而是一个相对于另一个的位置和速度。
说白了，就是我们把坐标系原点放在目标航天器上，然后看追踪航天器怎么动。这样问题就简化了很多。
我个人习惯把这种思路叫做“跟着目标走”。你想想看，在工程实践中，我们不可能每秒钟都去算两个航天器的绝对轨道，那太笨了。用相对运动模型，计算量小，物理意义也直观。
3.2 坐标系定义
在推导之前，先定好坐标系。这里用的是LVLH坐标系（Local Vertical Local Horizontal），也叫轨道坐标系。
- 原点：目标航天器的质心
- x轴：沿目标航天器径向，指向地心方向（R方向）
- y轴：沿目标航天器速度方向（V方向），也就是切向
- z轴：垂直于轨道平面，按右手定则确定（W方向）
嗯，这里要注意：x轴是指向地心的，不是背向地心。很多初学者会搞反，我当年也犯过这个错。
重要概念：LVLH坐标系是一个旋转坐标系，它随着目标航天器一起在轨道上转动。所以我们在写动力学方程时，必须考虑坐标系的旋转效应。
3.3 推导前的假设条件
任何模型都有假设，C-W方程也不例外。我建议你把这些假设背下来，因为在实际工程中，一旦偏离了这些假设，结果就可能出大问题。
- 目标航天器在圆轨道上运行——偏心率e=0。这是最关键的假设。
- 仅考虑地球中心引力——忽略J2摄动、大气阻力、太阳光压等。
- 追踪航天器与目标航天器距离很近——相对于轨道半径来说是小量。
- 追踪航天器不施加控制力——推导的是自由运动方程。
避坑指南：我曾经在一个项目中，直接用C-W方程去设计椭圆轨道交会的控制律，结果仿真时发现误差越来越大。后来才意识到，目标轨道偏心率0.05，已经超出了方程的适用范围。所以，用之前一定要检查轨道偏心率！
3.4 推导过程
好，现在开始推导。我会尽量写得简洁，但关键步骤不会省。
首先，在惯性坐标系中，追踪航天器的绝对运动方程为：
r̈_t = -μ * r_t / |r_t|³目标航天器的绝对运动方程为：
r̈_c = -μ * r_c / |r_c|³其中，r_t和r_c分别是追踪航天器和目标航天器的位置矢量。μ是地球引力常数。
定义相对位置矢量：
ρ = r_t - r_c我们的目标就是求出ρ的动力学方程。
这里有个技巧：因为LVLH坐标系在旋转，所以ρ在惯性系中的二阶导数，和它在旋转系中的二阶导数之间，有一个转换关系。这个关系涉及到角速度和角加速度。
经过一系列推导（嗯，中间过程我就不一步步写了，教科书上都有），最终得到：
ẍ - 2nẏ - 3n²x = 0
ÿ + 2nẋ = 0
z̈ + n²z = 0这就是大名鼎鼎的C-W方程。其中n是目标航天器的轨道角速度，n = √(μ/a³)，a是轨道半长轴。
个人经验：我建议你记住这个方程的形式。x方向有一个-3n²x项，这是重力梯度项；y方向和x方向有耦合项2nẋ和2nẏ，这是科里奥利力。z方向是独立的简谐振动。这些物理意义搞清楚了，后面设计控制律就顺手多了。
3.5 方程的解析解
C-W方程是线性常系数微分方程组，有解析解。这个解在工程中非常有用，比如可以用来做飞行的预测。
给定初始条件ρ₀ = [x₀, y₀, z₀]ᵀ，v₀ = [ẋ₀, ẏ₀, ż₀]ᵀ，解为：
x(t) = (ẋ₀/n) sin(nt) - (3x₀ + 2ẏ₀/n) cos(nt) + (4x₀ + 2ẏ₀/n)
y(t) = (2ẋ₀/n) cos(nt) + (6x₀ + 4ẏ₀/n) sin(nt) - (6nx₀ + 3ẏ₀)t + (y₀ - 2ẋ₀/n)
z(t) = z₀ cos(nt) + (ż₀/n) sin(nt)你看，z方向就是简谐振动，周期等于轨道周期。x和y方向耦合在一起，还有一个长期项——那个带t的项。这意味着，如果初始条件不满足特定关系，追踪航天器会慢慢漂走。
关键点：y方向的那个长期项，是C-W方程最重要的特征之一。它告诉我们，在圆轨道交会中，如果不加控制，两个航天器会沿着轨道方向漂移。这也是为什么交会对接时，我们通常从后方接近目标——因为沿轨道方向的漂移是不可避免的。
3.6 知识体系结构图
下面我用一张SVG图来总结本章的知识结构，方便你整体把握：
3.7 工程应用中的注意事项
最后，我想聊聊在实际工程中怎么用这个方程。毕竟，理论是理论，工程是工程，中间差着十万八千里。
- 适用范围：C-W方程只适用于目标轨道偏心率小于0.01的情况。如果偏心率大一点，就要用T-H方程（针对椭圆轨道）。
- 距离限制：相对距离一般不超过轨道半径的1%。比如在400km的轨道上，相对距离最好在4km以内。
- 摄动影响：J2摄动会引起轨道面进动，长时间交会时需要考虑。我做过一个项目，交会时间超过2小时，C-W方程的预测误差就明显变大了。
小技巧：在实际工程中，我通常会用C-W方程做初步设计和快速预测，然后用高精度数值仿真做验证。这样既保证了效率，又保证了精度。记住，任何模型都是工具，关键是要知道它的边界在哪里。
好了，关于C-W方程的推导和假设，就讲到这里。这个方程虽然简单，但它是整个交会对接相对运动控制的基石。后面的章节，我们会在这个基础上，加入控制力，设计各种交会策略。