4. C-W方程的解析解：自由运动解、脉冲机动解、周期性条件
各位同学，今天我们来聊聊C-W方程的解析解。说实话，这玩意儿在交会对接里太重要了。我当年刚入行时，师傅就告诉我：搞不懂C-W方程，就别想碰交会对接。嗯，这话虽然有点绝对，但确实有道理。
C-W方程，说白了就是描述两个航天器在近圆轨道上相对运动的线性化方程。它把复杂的非线性问题简化了，让我们能用解析方法去分析相对运动。我个人习惯把它的解分成三类：自由运动解、脉冲机动解，还有周期性条件。咱们一个一个来看。
核心要点：C-W方程适用于目标航天器在近圆轨道、追踪器与目标距离远小于轨道半径的情况。说白了就是“小偏差、近圆轨”这两个前提。
4.1 自由运动解
自由运动，就是没有外力作用下的相对运动。你想想看，两个航天器在太空中，如果不施加控制力，它们会怎么动？
C-W方程的标准形式是这样的：
ẍ - 2nẏ - 3n²x = 0
ÿ + 2nẋ = 0
z̈ + n²z = 0其中n是目标轨道角速度。解这个方程组，我一般会先看z方向——它独立于x和y，就是个简谐振动。x和y方向是耦合的，解起来稍微麻烦点。
自由运动的解析解为：
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)这里x₀、y₀、z₀是初始位置，ẋ₀、ẏ₀、ż₀是初始速度。看着挺复杂是吧？其实拆开来看就清楚了：
- z方向：就是个简谐振动，周期等于轨道周期
- x-y平面：包含周期项和长期项。那个带t的项就是长期漂移项
我的经验：我在做某型号交会对接仿真时，发现自由运动解里那个长期项特别关键。它决定了追踪器会不会飘走。如果初始条件选得不好，几圈下来两个航天器就分道扬镳了。
4.2 脉冲机动解
实际工程中，我们很少连续施加推力，更多是用脉冲机动——说白了就是“打一下，然后自由飞一会儿”。
脉冲机动解的核心思想是：在某个时刻施加一个速度增量Δv，然后系统就按新的初始条件自由运动。我习惯用状态转移矩阵来表达：
[x(t), y(t), z(t), ẋ(t), ẏ(t), ż(t)]ᵀ = Φ(t, t₀) · [x₀, y₀, z₀, ẋ₀, ẏ₀, ż₀]ᵀΦ(t, t₀)就是状态转移矩阵，它把初始状态映射到任意时刻的状态。对于脉冲机动，我们只需要在机动时刻更新速度项：
ẋ⁺ = ẋ⁻ + Δvₓ
ẏ⁺ = ẏ⁻ + Δvᵧ
ż⁺ = ż⁻ + Δv_z然后继续用自由运动解传播就行了。简单吧？
注意：我曾经犯过一个错误——把脉冲机动当成连续推力来处理，结果仿真结果跟实际飞行数据对不上。后来才意识到，脉冲机动是瞬时改变速度，位置不变。这个细节千万别搞错。
4.3 周期性条件
为什么要讨论周期性条件？因为交会对接中，我们希望追踪器能“乖乖地”绕着目标转，而不是越飘越远。周期性条件就是保证相对运动不出现长期漂移的条件。
从自由运动解可以看出，y方向那个带t的项要想消失，必须满足：
6nx₀ + 3ẏ₀ = 0化简一下：
ẏ₀ = -2nx₀这就是周期性条件！满足这个条件，相对运动就是周期性的，轨迹会形成一个闭合的椭圆或圆。
我整理了一下周期性运动的特点：
| 参数 | 条件 | 物理意义 |
|---|---|---|
| 初始径向速度 | ẋ₀任意 | 不影响周期性 |
| 初始法向速度 | ẏ₀ = -2nx₀ | 消除长期漂移 |
| 初始横向位置 | y₀任意 | 决定椭圆中心位置 |
| z方向 | 自动周期 | 独立简谐振动 |
实用技巧：如果你想让追踪器在目标前方100米做周期性绕飞，设x₀ = -100m，然后ẏ₀ = -2n·(-100) = 200n m/s。这样它就会在x-y平面画出一个椭圆。我在某次在轨试验中就用这个办法，效果非常好。
周期性条件在实际工程中太有用了。比如空间站附近的货运飞船，经常需要保持在一个相对稳定的位置，等指令来了再靠近。这时候就要用周期性条件来设计初始速度。
嗯，这里还要提一句：周期性条件只适用于x-y平面内的运动。z方向本来就是周期的，不用额外操心。
最后总结一下我的体会：C-W方程的解析解看着公式多，但核心就三条——自由运动解告诉你无控状态下的运动规律，脉冲机动解告诉你怎么用最少的燃料去改变运动状态，周期性条件则告诉你如何设计一个稳定的相对轨道。这三条搞明白了，交会对接的轨道控制问题就解决了一大半。