“剩余速度”制导法入门教程

一、教程目标

本教程要帮助读者完成一件明确的任务：理解并实现一套基于椭圆轨道理论的“剩余速度”制导方法，使飞行器从当前状态出发，通过不断消除与目标轨道之间的速度差，最终转移到一条经过期望落点的轨道上。

适合阅读本文的读者包括：从事航天任务仿真、轨道机动算法研究的开发者；使用轨道仿真工具进行任务设计的工程人员；以及具备一定数学与物理基础、希望理解制导算法内在逻辑的技术爱好者。即使你此前没有接触过轨道力学，只要具备中学阶段的解析几何与牛顿力学基础，也能跟上本文的思路。

需要提前说明的是，本文讨论的方法建立在若干理想化假设之上：将飞行器视为质点、忽略空气阻力与摄动力、基于二体问题下的椭圆轨道理论。这些假设在真实工程中需要进一步修正，但作为理解制导逻辑的起点，它们能让推导过程保持清晰。

二、适用场景

“剩余速度”制导法适用于以下典型场景：

• 轨道转移与落点制导：已知期望落点，需要反推当前应具有的速度方向与大小。
• 仿真环境中的制导验证：在轨道仿真环境中快速评估不同弹道倾角下的转移方案。
• 教学与算法原型开发：作为理解“正问题”（由速度增量推轨道）与“逆问题”（由目标轨道推速度增量）对偶关系的入门案例。

如果你熟悉轨道机动节点规划的概念，可以把本方法理解为它的逆过程：轨道规划是给定速度增量、计算转移后的轨道；而“剩余速度”法是给定目标轨道、反推出所需的速度增量。

三、前置准备

在动手实现之前，建议先确认以下基础条件：

1. 数学基础：理解椭圆的标准方程、焦点、离心率、半长轴、真近点角等基本概念；了解向量运算与三角函数。
2. 物理基础：掌握牛顿万有引力定律与二体问题下的轨道运动规律，理解“速度方向与水平面夹角”这一几何量的含义。
3. 计算环境：任意支持浮点运算与向量运算的编程语言即可，无需特定框架。建议准备一个可绘制轨道曲线的绘图库，用于可视化验证。
4. 坐标系约定：明确你使用的惯性坐标系（如以星球中心为原点的二维平面坐标系），并约定角度正方向。后续所有推导都依赖这一约定的一致性。

四、核心思路概述

整个方法的逻辑链条可以概括为四步：

1. 将飞行器视为质点，忽略阻力等干扰，问题退化为二体问题下的椭圆轨道几何问题。
2. 根据期望落点与当前位置，确定一条经过这两点的椭圆轨道。
3. 计算该轨道在当前位置处应具有的速度（称为“目标速度”）。
4. 将目标速度与当前实际速度作差，得到“剩余速度”，通过控制发动机推力方向不断消除该差值，飞行器便逐步转移到目标轨道上。

这里有一个容易被忽略的细节：仅凭两个点的坐标无法唯一确定一条椭圆轨道。过两点的椭圆有无数条，因此必须引入额外的约束参数，才能把轨道“锁定”下来。这正是下一节要解决的问题。

五、实施步骤

步骤一：引入弹道倾角作为约束参数

做什么：在已知当前位置与目标位置的基础上，额外指定当前位置处的弹道倾角（即速度方向与当地水平面的夹角），以此唯一确定一条椭圆轨道。

为什么这么做：两个点只提供两个几何约束，而椭圆轨道需要更多参数才能确定。可供选择的附加参数包括离心率、半长轴、真近点角等，选择弹道倾角的好处是它具有明确的物理意义，且与发动机推力方向直接相关，便于后续制导控制。

注意什么：弹道倾角的定义必须与坐标系约定一致。通常规定水平飞行时倾角为零，竖直向上为正方向。不同约定会导致公式符号差异，实现时务必统一。

步骤二：将目标坐标转换到惯性空间

做什么：把目标点的经纬度坐标转换为惯性坐标系下的位置矢量，并考虑星球自转带来的位置漂移。

为什么这么做：星球（包括地球及多数天体）具有自转角速度，目标点在星球表面上的位置在惯性空间中并非静止，而是随时间缓慢移动。而我们要确定的轨道是在惯性空间中静止的，因此必须先把目标位置换算到同一参考系下，否则计算出的轨道会存在系统性偏差。

注意什么：换算时需要知道目标点相对星球本体的固定坐标、星球自转角速度以及当前时刻。若飞行时间较长，还需考虑飞行过程中目标点的移动，必要时采用迭代方式修正落点。

步骤三：求解椭圆轨道并计算目标速度

做什么：根据当前位置、惯性空间中的目标位置以及指定的弹道倾角，求解出唯一的椭圆轨道，并计算该轨道在当前位置处的速度矢量。

为什么这么做：这一步是整个方法的“逆问题”核心——由几何约束反推运动状态。求解过程通常涉及椭圆轨道方程的联立，得到一个关于轨道参数的方程，解出后即可代入速度公式。

示意性伪代码：

输入：当前位置 r0，目标位置 rf，弹道倾角 gamma
1. 建立以星球中心为原点的平面坐标系
2. 由 r0、rf 与 gamma 联立椭圆轨道方程
3. 解出轨道参数（离心率、半长轴等）
4. 代入活力公式与方向约束，得到目标速度矢量 v_target
5. 计算剩余速度 dv = v_target - v_current
输出：dv

注意什么：方程可能存在多解或无解的情况（例如给定倾角下无法命中目标），实现时应加入解的合法性判断，并对无效输入给出明确提示。

步骤四：消除剩余速度

做什么：将剩余速度矢量作为推力方向的参考，控制发动机沿该方向施加加速度，使当前速度逐步逼近目标速度。

为什么这么做：剩余速度本质上就是当前状态与目标轨道之间的“速度缺口”。只要持续减小这个缺口，飞行器的轨道就会收敛到经过目标点的那条椭圆轨道上。

注意什么：实际推力有限，消除过程需要时间，期间飞行器位置会变化，因此剩余速度需要实时重算，形成闭环控制。仿真中可采用小步长迭代，每步更新当前位置与速度。

步骤五：结果验证

做什么：检查飞行器最终是否经过期望落点。

验证方法：在仿真中记录飞行轨迹，绘制轨道曲线，确认其与目标点相交；同时检查末端剩余速度是否收敛到接近零。若存在偏差，回溯弹道倾角选取与目标坐标换算环节。

六、常见问题与排查

• 轨道无法经过目标点：多因弹道倾角选取超出可行范围，尝试调整倾角或改用其他约束参数。
• 落点存在系统性偏移：检查是否遗漏星球自转修正，确认目标坐标已转换到惯性空间。
• 剩余速度反复震荡不收敛：可能是推力方向更新过于激进，减小步长或加入阻尼项。
• 多解导致轨道异常：对方程解加入物理合理性筛选，如排除负半长轴或超出合理范围的离心率。

七、优化建议

从稳定性角度，建议对剩余速度做限幅处理，避免初始偏差过大导致推力方向剧烈跳变。从可维护性角度，将轨道求解、坐标换算、速度计算拆分为独立函数，便于单元测试与替换约束参数。从精度角度，可在长飞行时间场景中引入迭代落点预测，逐步修正自转带来的目标漂移。

八、总结

本文围绕“剩余速度”制导法，从约束参数引入、惯性空间坐标换算、椭圆轨道反解到闭环消除速度差，给出了完整的实施路径。关键要点有三：用弹道倾角锁定唯一轨道、把目标位置换算到惯性空间、以剩余速度作为闭环控制量。掌握这套逻辑后，你可以进一步尝试替换约束参数、引入摄动修正，或将其与轨道规划方法对照学习，加深对正逆问题对偶关系的理解。