Tschauner-Hempel方程：航天器椭圆轨道相对运动控制的核心技术

1. 椭圆轨道相对运动：航天器编队飞行的核心挑战

在航天器编队飞行任务中，精确描述和预测两颗航天器之间的相对运动轨迹是确保任务成功的关键技术。传统方法在处理椭圆参考轨道时往往面临计算复杂、精度不足的问题，直到Tschauner-Hempel方程的提出才从根本上解决了这一难题。

这个看似晦涩的微分方程，实际上支撑着现代卫星编队、空间站补给、深空探测等众多关键任务。想象一下，当货运飞船需要与国际空间站对接时，两者轨道高度不同（空间站约400km，货运飞船初始轨道可能低20km），且都运行在近圆椭圆轨道上。此时，Tschauner-Hempel方程就能精确计算出货运飞船应该如何调整速度和方向，才能安全、高效地完成交会对接。

2. Tschauner-Hempel方程的数学本质

2.1 从Clohessy-Wiltshire方程到TH方程的演进

1960年提出的Clohessy-Wiltshire(CW)方程是描述圆轨道相对运动的经典模型，其基本形式为：

math复制\ddot{x} - 2n\dot{y} - 3n^2x = 0 \\
\ddot{y} + 2n\dot{x} = 0 \\
\ddot{z} + n^2z = 0

其中n为平均角速度。但当参考轨道偏心率e>0.01时，CW方程的误差会急剧增大。1965年，Tschauner和Hempel通过引入真近点角f作为自变量，建立了适用于任意偏心率的统一模型：

math复制\frac{d^2x}{df^2} - 2\frac{dy}{df} - \frac{3x}{1+e\cos f} = 0 \\
\frac{d^2y}{df^2} + 2\frac{dx}{df} = 0 \\
\frac{d^2z}{df^2} + z = 0

这个转变的关键在于：用真近点角替代时间变量，使得方程不再显含时间项，从而消去了轨道偏心率带来的周期性扰动。

2.2 方程解析解的构造方法

TH方程的解析解可以通过以下步骤构建：

1. 引入辅助变量η=1+e·cosf
2. 将径向运动方程转化为二阶变系数微分方程
3. 采用Floquet理论求解周期系数系统
4. 最终得到包含4个积分常数的通解：

math复制x(f) = C_1\frac{e}{\eta}\s