4. CW方程解析解与自由运动轨迹
在圆轨道参考航天器附近，当两航天器距离远小于轨道半径时，相对运动可用Clohessy-Wiltshire（CW）方程描述。设参考航天器运行于半径 rc 的圆轨道，轨道角速度 n=√(μ/rc3)，以参考航天器质心为原点建立LVLH坐标系：x轴沿径向向外，y轴沿轨道切向（速度方向），z轴沿轨道面法向。追踪器相对位置记为 r=[x, y, z]T，相对速度记为 v=[ẋ, ẏ, ż]T。CW方程在三个方向上的耦合关系为：
ẍ − 3n²x − 2nẏ = 0
ÿ + 2nẋ = 0
z̈ + n²z = 0可以看到，轨道面内（x、y）与面外（z）解耦：z方向是频率为 n 的简谐振荡，而x、y方向通过科氏项 2nẏ 与 2nẋ 相互耦合，这正是自然绕飞与漂移现象的根源。
状态转移矩阵：解析解的统一表达
CW方程是线性时不变系统，其解析解可写成状态转移矩阵形式。取状态向量 X=[x, y, z, ẋ, ẏ, ż]T，则任意时刻 t 的状态为 X(t)=Φ(t)X(0)，其中 Φ(t) 为6×6状态转移矩阵。记 s=sin(nt)，c=cos(nt)，面内部分为：
Φ_xy(t) =
[ 4−3c 0 s/n 2(1−c)/n ]
[ 6(s−nt) 1 −2(1−c)/n (4s−3nt)/n ]
[ 3ns 0 c 2s ]
[ −6n(1−c) 0 −2s 4c−3 ]面外部分为独立的2×2振荡矩阵：z(t)=z₀cos(nt)+(ż₀/n)sin(nt)，ż(t)=−z₀n sin(nt)+ż₀cos(nt)。状态转移矩阵的价值在于：给定初始相对位置与速度，无需数值积分即可直接算出任意时刻的位置和速度，是轨迹设计、机动规划与偏差传播的基础工具。
自然绕飞轨迹与周期性条件
若初始条件满足特定关系，追踪器将围绕参考航天器做周期性闭合绕飞。最典型的是面内椭圆绕飞：取初始位置 x₀=0、y₀=2ρ，初始速度 ẋ₀=−nρ、ẏ₀=0，代入解析解可得
x(t) = ρ·sin(nt)
y(t) = 2ρ·cos(nt)即轨迹为半长轴沿y方向、半短轴沿x方向的椭圆，长轴是短轴的2倍，周期等于参考轨道周期 T=2π/n。这种“2:1椭圆”是自然绕飞的经典形态，追踪器在一个轨道周期内相对参考航天器完成一圈闭合运动，无需任何推力维持。
周期性条件可从解析解中提取。面内运动出现长期项（含 nt 的项）的条件是 ẏ₀=−2nx₀，此时y方向漂移被消除，运动退化为周期椭圆。更一般地，要使面内轨迹闭合，需满足：
- 初始径向速度 ẋ₀ 与切向位置 y₀ 满足 ẋ₀ = −n·y₀/2 时，径向振荡幅度与切向漂移解耦；
- 初始切向速度满足 ẏ₀ = −2n·x₀，消除长期漂移项；
- 面外运动本身即为周期振荡，只要 z₀ 与 ż₀ 有限，就始终闭合，周期为 2π/n。
当上述条件同时满足时，追踪器相对参考航天器的运动完全周期化，形成稳定的自然绕飞。若仅满足部分条件，则轨迹呈现“椭圆叠加漂移”的形态，需要后续机动修正。
漂移与发散特性分析
CW方程解析解中最需要警惕的是长期项。当 ẏ₀ ≠ −2nx₀ 时，y方向解中出现 −(3nt)x₀ 或 (4s−3nt)/n·ẏ₀ 形式的项，这些项随时间线性增长，导致追踪器沿轨道切向持续漂移。漂移速率近似为
v_drift ≈ −(3/2)·n·x₀ （当 ẏ₀=0 时）这意味着即使初始径向位置只有几十米，经过若干轨道周期后切向漂移可达数公里。例如取 n=1.1×10⁻³ rad/s（约500 km高度圆轨道），x₀=100 m，则每圈漂移约 3πx₀ ≈ 942 m，10圈后漂移近10 km。这种发散特性对交会接近任务极为关键：若不加控制，追踪器会逐渐远离参考航天器，无法维持稳定的相对构型。
面外运动不存在长期项，z方向始终是幅值有限的简谐振荡，因此法向偏差不会发散，只会周期性往复。这一特性在轨道面选择与接近走廊设计中常被利用：将主要接近机动安排在面内，而面外偏差可通过自然振荡“消化”。
轨迹形态与参数关系图解
下面用一张参数关系图汇总CW解析解中初始条件对轨迹形态的影响。图中左侧为初始状态输入，中间为解析解分支，右侧为对应的轨迹类型与长期行为，便于快速判断“给定初始条件会得到什么轨迹”。
工程应用要点
在实际交会接近轨迹设计中，CW解析解主要用于三类工作。第一，快速构型设计：给定期望的绕飞椭圆尺寸ρ，直接反推初始条件 x₀=0、y₀=2ρ、ẋ₀=−nρ、ẏ₀=0，即可获得稳定的自然绕飞。第二，漂移评估：在接近段若无法精确满足周期性条件，可用 v_drift≈−(3/2)n·x₀ 估算漂移速率，判断多久需要一次修正机动。第三，偏差传播：利用状态转移矩阵将初始导航误差映射到终端时刻的位置误差，为交会走廊与安全区设计提供依据。
需要强调的是，CW方程假设参考轨道为圆轨道、相对距离远小于轨道半径、无摄动力。当参考轨道偏心率较大或相对距离达到数十公里以上时，解析解精度下降，需引入T-H方程或数值积分修正。但在近程交会（几十米至几公里）的初步设计与分析阶段，CW解析解仍是效率最高、物理意义最清晰的工具。