Apollo卫星地月系统运动模型与轨道仿真

资 源 简 介

本项目主要用于构建和描述Apollo卫星在地球-月球系统中的运动数学模型，核心基于限制性三体问题（Restricted Three-Body Problem）的动力学方程。该项目提供了一个标准的MATLAB函数文件，用于计算定义的微分方程组，描述卫星在两个主要天体（地球和月球）引力作用下的平面运动轨迹。具体而言，函数接受当前的时间标量和包含位置及速度的状态向量作为输入，根据牛顿运动定律和旋转坐标系下的动力学关系（包含科里奥利力和离心力项），计算并返回状态变量的时间导数。该模型是进行航天动力学数值模拟的基础，能够直接与MATLAB的ode45、ode113等常微分方程数值求解器无缝对接，用于精确积分Apollo卫星的飞行轨迹，分析其在拉格朗日点附近的稳定性，或进行地月借力飞行等复杂轨道机动的仿真实验。

详 情 说 明

Apollo卫星数学模型与轨道仿真

项目简介

本项目旨在构建和模拟Apollo卫星在地球-月球系统中的运动轨迹。核心数学模型基于限制性三体问题 (Restricted Three-Body Problem, RTBP) 的动力学方程。通过MATLAB数值计算平台，项目实现了对卫星在两个主要天体（地球和月球）引力及旋转参考系下科里奥利力和离心力共同作用下的平面运动进行精确积分。

该仿真程序不仅能够根据给定的初始状态生成复杂的轨道几何，还包含完整的分析工具，用于验证数值积分的能量守恒特性（雅可比常数）以及计算系统中关键的力学平衡点（拉格朗日点）。

功能特性

• 高精度动力学模型：基于旋转坐标系下的牛顿运动方程，包含地球引力、月球引力、离心力及科里奥利力项。
• 数值积分求解：利用MATLAB的 ode45 求解器，配合高精度误差控制（相对误差和绝对误差均为 1e-12），进行长时间跨度的轨道推演。
• 多维度可视化：

* 轨迹图：在旋转坐标系中绘制卫星相对于地球和月球的平面运动轨迹。 * 状态演化图：提供位置 ($x, y$) 和速度 ($v_x, v_y$) 随无量纲时间变化的详细曲线。 * 守恒量验证：实时计算并绘制雅可比常数（Jacobi Constant），用于评估数值解的稳定性和精度。

• 拉格朗日点分析：自动计算并标记地月系中的全部五个拉格朗日点（L1-L5），其中共线点（L1, L2, L3）采用数值迭代求解，三角点（L4, L5）采用解析解。

系统要求

• MATLAB R2016a 或更高版本
• Optimization Toolbox（用于 fzero 函数求解拉格朗日点）

使用方法

1. 将项目代码保存为MATLAB脚本文件。
2. 直接运行主函数。
3. 程序将自动执行以下操作：

* 初始化地月系统参数（质量比 $mu$）。 * 设置特定的初始状态向量（模拟周期性轨道或特定Apollo任务特征）。 * 调用数值积分器求解微分方程。 * 弹出三个图形窗口展示仿真结果。

1. 观察控制台输出，了解积分进度及雅可比常数的统计均值。

代码实现与算法详解

本项目的主要逻辑封装在一个主流程中，通过定义参数、求解ODE、后处理分析三个阶段实现。

1. 参数初始化与ODE求解

程序首先定义了地月系统的无量纲质量参数 $mu approx 0.012277$，这决定了两个主天体的引力强度分布。初始状态向量 $y_0 = [x_0, y_0, v_{x0}, v_{y0}]$ 被硬编码为一组特定的数值，这组初值对应于三体问题中的一种典型轨道（可能是周期解或借力飞行轨道）。

时间跨度设定为约17个无量纲时间单位。为了保证轨道积分在长时间运行下的闭合性或准确性，求解器 ode45 的相对误差 (RelTol) 和绝对误差 (AbsTol) 被严格设定为 $10^{-12}$。

2. 核心动力学方程 (Satellite Dynamics)

这是系统的物理核心，描述了卫星在旋转坐标系中的运动。

• 坐标系定义：地球位于 $(-mu, 0)$，月球位于 $(1-mu, 0)$。
• 距离计算：实时计算卫星到地球 ($r_1$) 和到月球 ($r_2$) 的距离。
• 加速度计算：根据限制性三体问题的二阶微分方程计算加速度：

* $a_x = 2v_y + x - frac{(1-mu)(x+mu)}{r_1^3} - frac{mu(x-1+mu)}{r_2^3}$ * $a_y = -2v_x + y - frac{(1-mu)y}{r_1^3} - frac{mu y}{r_2^3}$ * 方程中清晰地体现了科里奥利力项 ($2v_y, -2v_x$) 和包含引力及离心力的势能梯度项。

3. 拉格朗日点计算与绘制

程序通过混合算法确定五个拉格朗日点的位置：

• L4, L5：利用几何解析解直接计算，位于地球和月球构成的等边三角形顶点。
• L1, L2, L3：这三个点位于两个天体的连线（x轴）上。程序通过构建关于位置 $x$ 的非线性方程（有效势能导数为零），并利用 fzero 函数在特定的区间内（如地月之间、月球外侧、地球外侧）进行数值搜索，精确求解其坐标。

4. 雅可比常数验证

在限制性三体问题中，雅可比常数 $C$ 是唯一的运动积分（能量守恒量）。程序包含一个专门的函数用于事后验证：

• 计算公式：$C = 2Omega(x,y) - (v_x^2 + v_y^2)$
• 其中伪势能 $Omega = frac{1}{2}(x^2+y^2) + frac{1-mu}{r_1} + frac{mu}{r_2}$。
• 通过绘制 $C$ 随时间的变化曲线，用户可以直观地判断数值积分过程中是否存在严重的能量耗散或计算误差。

5. 可视化输出

• Figure 1 (轨迹)：在 $x-y$ 平面上绘制蓝色的卫星轨迹，并用不同大小和颜色的标记点画出地球和月球，同时叠加显示计算出的所有拉格朗日点位置。
• Figure 2 (状态)：使用四个子图分别展示 $x, y, v_x, v_y$ 随时间的演变，便于分析轨道在特定方向上的振荡或逃逸趋势。
• Figure 3 (验证)：展示雅可比常数的时间历程，理想情况下应为一条水平直线。