1. 课程导入与轨道动力学概览
航天器轨道动力学研究的是航天器在引力及其他外力作用下的运动规律。它要回答三个层层递进的问题：航天器在哪里、将要去哪里、如何让它按人的意图到达并保持在目标位置。第一个问题属于轨道确定，第二个属于轨道预报，第三个属于轨道控制与任务设计。三者共同构成轨道动力学的完整闭环，也是本课程后续所有章节的骨架。
从建模角度看，轨道动力学的核心是建立运动方程。最基础的模型是二体问题：把地球视为质点，航天器视为质点，只考虑中心引力。此时轨道是严格的圆锥曲线，能量与角动量守恒，六个轨道根数（半长轴 a、偏心率 e、倾角 i、升交点赤经 Ω、近地点幅角 ω、真近点角 f）唯一确定轨道状态。二体模型给出了解析解，是理解轨道形状、周期、速度与位置关系的基础。但真实航天器还受到地球非球形引力（主要是 J2 项）、大气阻力、太阳光压、第三体引力等摄动。这些摄动力虽小，却会随时间累积，导致轨道根数缓慢漂移。因此工程上通常采用"二体解析 + 摄动修正"的建模思路：用二体模型快速估算，用数值积分或摄动方程精确预报。
仿真的作用在于把模型变成可计算、可验证、可复现的过程。给定初始状态和力模型，通过数值积分器（如 RK4、RKF45、Dormand-Prince）逐步推进状态向量，就能得到任意时刻的位置和速度。仿真的价值不仅在于"算出来"，更在于支撑设计决策：例如评估某轨道在 5 年内的漂移量、比较不同控制策略的燃料消耗、验证星座构型是否满足覆盖要求。下面这张图给出了从任务需求到轨道方案再到仿真验证的整体知识地图，帮助建立全局认识。
图中自上而下展示了任务需求、建模手段与仿真验证的对应关系。以对地观测任务为例，若要求卫星在 700 km 高度、降交点地方时 10:30 成像，就需要选择太阳同步轨道。太阳同步条件由 J2 摄动引起的升交点赤经进动率决定：当轨道倾角约 98.2°、半长轴约 7078 km 时，进动率约为 0.9856°/天，恰好与地球绕日公转的平均角速率匹配，从而保证光照条件稳定。这一设计过程正是"任务需求→轨道根数→摄动分析→仿真验证"的典型链条。
再看星座维持。以 Walker 24/6/2 构型为例，24 颗卫星分布在 6 个轨道面，每面 4 颗，相位因子为 2。理想情况下各卫星保持固定相对几何，但 J2 摄动会使不同轨道面的升交点赤经以不同速率漂移，导致面间夹角变化，覆盖性能下降。仿真时需要输出各卫星的星下点轨迹、面间相位差随时间的变化曲线，并据此设计保持策略：当相位偏差超过阈值（例如 ±2°）时施加小脉冲修正。这类问题的关键不是单颗卫星的轨道，而是多星之间的相对运动与长期演化。
深空探测则带来另一类挑战。地月转移轨道需要同时考虑地球、月球引力，属于限制性三体问题。仿真中常用圆型限制性三体模型（CR3BP）寻找平动点与转移轨道，再通过微分修正将初值收敛到目标轨道。与近地轨道相比，深空轨道的力模型更复杂、积分时间更长、对初值更敏感，因此数值积分器的精度和稳定性尤为关键。
本课程后续将按"基础模型→摄动分析→数值方法→任务应用"的顺序展开。学习时建议把握三条主线：一是坐标系与时间系统的统一，所有位置、速度、力都必须在同一参考系下描述；二是量纲与单位的一致性，长度用 km、时间用 s、速度用 km/s、角度用 rad 或 deg，混用会导致数量级错误；三是"先解析估算、再数值精算"的工作习惯，先用二体模型判断可行性，再用摄动模型评估长期演化。下面给出一个最小可运行的仿真参数示例，用于建立直观感受。
# 示例：近地圆轨道二体仿真参数（合理假设，非实测）
mu_earth = 398600.4418 # km^3/s^2，地球引力常数
r_earth = 6378.137 # km，地球赤道半径
altitude = 700.0 # km，轨道高度
a = r_earth + altitude # km，半长轴 ≈ 7078.137
e = 0.0 # 圆轨道
i = 98.2 # deg，太阳同步倾角
T = 2*pi*sqrt(a**3/mu_earth) # 轨道周期 ≈ 5927 s ≈ 98.8 min
v_circ = sqrt(mu_earth/a) # 圆轨道速度 ≈ 7.50 km/s
运行上述参数可得到：轨道周期约 98.8 分钟，每天绕地球约 14.6 圈，星下点轨迹在赤道附近西移约 24.6°（由地球自转引起）。若加入 J2 摄动，升交点赤经每天进动约 0.9856°，与太阳同步条件吻合。这些数字是后续章节反复使用的基准，建议读者在首次仿真时逐一核对，确认模型、单位与坐标系设置正确。
轨道动力学不是孤立的数学推导，而是连接任务需求与工程实现的桥梁。理解研究对象、掌握建模手段、善用仿真验证，是后续学习轨道确定、轨道控制、星座设计等内容的共同基础。