本科航空航天工程专业三年级专业核心课教案：轨道力学中的微分方程建模与数值优化

  本科航空航天工程专业三年级专业核心课教案：轨道力学中的微分方程建模与数值优化

一、教学理念与背景深度分析

本教学设计立足于新时代工程教育范式——“新工科”建设的核心要求，旨在打破传统《轨道力学》或《航天器动力学》课程中数学工具讲授与工程应用相对割裂的现状。我们秉持“数学是工程的灵魂语言，计算是设计的实践之手”这一核心理念，将课程定位为一门深度融合数学理论、计算实践与工程设计的跨学科、项目式高阶专业课程。

学生群体为本科三年级航空航天工程专业学生。其前置知识体系应包括：工科数学分析（重点在微分方程、多元微积分）、线性代数、理论力学（刚体动力学、拉格朗日方程）、程序设计基础（如Python或MATLAB）。然而，传统教学反馈表明，学生普遍存在“知识孤岛”现象：虽能解算标准微分方程，却不理解其物理背景的简化与抽象过程；虽知晓数值算法原理，却无法将其有效迁移至非标准的工程微分方程组求解；虽了解轨道参数，却难以从第一性原理出发，自主完成一个完整轨道的设计与分析流程。

因此，本教案的设计逻辑是：以复杂的真实工程问题（如地火转移轨道优化）为驱动，以建立并求解控制方程为核心数学任务，以数值仿真与可视化作为验证与探索的主要手段。教学重心从“传授轨道知识”转变为“锻造通过数学与计算解决轨道设计问题的能力”。我们强调“设计即建模，建模即计算”，引导学生体验从物理世界抽象出数学模型、选择或构造数值方法、实现代码求解、分析结果并迭代优化的完整科研与工程闭环。这不仅是对学生已有知识的系统整合与升华，更是对其未来从事航天器系统设计、深空探测任务规划等高端工作的核心胜任力的关键培养。

二、教学目标

本课程为连续8学时（两次4学时连堂）的专题模块，教学目标分为三个维度：

1.知识与技能维度：

1.核心概念深化：能准确阐述二体问题、限制性三体问题的基本假设及其适用边界；深入理解轨道能量、角动量、轨道根数与状态向量（位置、速度）之间的数学转换关系。

2.数学建模能力：能独立从牛顿万有引力定律出发，推导出在惯性坐标系和旋转坐标系下的质点运动微分方程（常微分方程组，ODEs）；能将摄动力（如地球非球形引力、大气阻力、太阳光压）项作为微分方程的右函数项进行数学表达。

3.数值算法应用：掌握常微分方程初值问题数值解法的基本原理，重点对比分析龙格-库塔法（特别是四阶经典Runge-Kutta，RK4）与线性多步法（如Adams-Bashforth）的优缺点、稳定性与适用场景；理解轨道预报中“轨道积分”的数值实质。

4.计算实现技能：熟练运用科学计算语言（以Python为例），编程实现上述微分方程组的数值积分模块；实现轨道根数与直角坐标的相互转换函数；实现轨道三维可视化与二维投影绘图。

2.过程与方法维度：

1.系统性设计流程：经历完整的轨道初步设计过程：任务目标分析→动力学模型选择与简化→建立数学模型（微分方程组）→初始条件确定与猜测→数值积分求解→结果可视化分析与判读→基于误差或性能指标的设计迭代。

2.问题分解与算法选择：面对复杂的轨道动力学问题，能将其分解为多个可数值求解的子问题；能根据问题的精度要求、计算效率需求和受力模型的复杂性，合理选择或组合数值积分算法。

3.计算思维培养：形成通过编写程序进行“计算实验”以探索参数空间、验证设计假设的思维习惯。学会利用数值工具进行敏感性分析，理解初始条件误差、模型误差对长期轨道预报的影响。

3.情感、态度与价值观维度：

1.工程严谨性：深刻认识轨道设计的高度精确性要求，理解数学模型的近似性及其带来的风险，树立“失之毫厘，谬以千里”的严谨工程态度。

2.创新与探索精神：鼓励在遵守物理定律的前提下，通过调整数值模型和算法参数进行创新性轨道构型探索（如弱稳定边界转移、混沌轨道利用），激发对深空探索前沿领域的兴趣。

3.跨学科协作意识：理解轨道设计作为航天任务大系统中的一个环节，其输入（如发射窗口、推进系统参数）与输出（如到达时间、燃料消耗）如何与其它分系统（测控、载荷、热控）紧密耦合，初步建立系统工程观念。

三、教学重点与难点

1.教学重点：

1.2.从物理模型到数学方程的严谨推导：重点讲解如何从矢量形式的牛顿第二定律和万有引力定律，导出在具体坐标系（地心惯性系）下的标量形式常微分方程组。这是后续一切数值计算的基石。

2.3.数值积分算法的工程化实现：不仅讲解RK4等算法的数学公式，更重点讲授其编程实现中的关键技巧，如右函数的编写、步长控制策略、数值稳定性判断，以及如何将高阶微分方程化为一阶方程组的标准形式。

3.4.轨道设计问题的完整求解流程：强调将数学工具嵌入到工程问题解决的标准流程中，使学生掌握从需求到代码、从代码到结果的分析方法。

5.教学难点：

1.6.模型简化与真实性的权衡：学生难以判断在何种任务背景下可以忽略哪些摄动力。例如，近地轨道需考虑大气阻力和J2项，而地月转移则需考虑三体效应。需要引导学生建立“模型复杂度与任务需求相匹配”的思维。

2.7.初始条件的逆向求解：给定目标轨道（如火星捕获轨道），反向求解所需的地球出发轨道（霍曼转移的初始速度），涉及超越方程的数值求解（如兰伯特问题），对学生的数值计算和编程能力提出高要求。

3.8.数值误差的累积与识别：长时间数值积分会导致误差累积，可能产生物理上不真实的结果（如能量不守恒）。教会学生通过监测轨道能量、角动量等积分常数来检验数值解的可靠性，是一项关键的、需要经验积累的技能。

四、教学资源与课前准备

1.硬件与环境：计算机网络教室，确保每位学生可使用安装有Python科学计算环境（Anaconda发行版，含NumPy,SciPy,Matplotlib库）的计算机。教师端配备多媒体教学系统，可进行代码演示和结果投屏。

2.软件与工具：

1.3.核心计算平台：JupyterNotebook。它将代码、公式、图表和文字说明集成在一个文档中，非常适合分步教学和过程重现。

2.4.扩展库：除基础科学计算库外，可引入Astropy或Skyfield库用于天文常数和星历获取，引入Plotly或Mayavi用于更生动的三维交互式可视化（若硬件允许）。

5.学习材料：

1.6.自编讲义《轨道设计的数学与计算实践》，包含理论概要、推导过程、算法伪代码和关键代码片段。

2.7.在线知识库（课前发布链接）：包含经典轨道力学教材章节（如《FundamentalsofAstrodynamics》）、数值分析算法百科（如SciPy文档）、NASA/JPL发布的真实轨道数据案例。

3.8.课前预习任务单：要求学生复习二体问题微分方程推导，自学Python中NumPy数组操作和Matplotlib基本绘图，并尝试编写一个简单的欧拉法求解一维运动方程的程序。

9.项目驱动案例：本模块的核心驱动项目为“地球-火星能量最优转移轨道初步设计与分析”。教师需准备该问题的背景资料，包括地球和火星的轨道平均半径、轨道周期、当前历元下的近似轨道根数，以及任务约束（如转移时间在200-300天之间）。

五、教学过程详细设计与实施（8学时）

第一篇章：从物理定律到可计算的模型（第1-2学时）

1.环节一：情境导入与问题锚定（20分钟）

1.2.教师活动：展示“天问一号”火星探测器发射、轨道修正、火星捕获的动画视频。提出核心驱动问题：“如果我们接到一个任务，要设计一条从地球到火星的转移轨道，并确保探测器在预定时间以预定速度到达火星附近，我们该如何用数学和计算机来完成这项设计？”引导学生认识到，这远非一个简单的圆锥曲线拼接问题，而是一个受多种力作用、需要高精度预报和可能优化的动态过程。

2.3.学生活动：观看视频，参与讨论。基于前置知识，提出初步设想（如霍曼转移），并思考其忽略的因素（行星引力摄动、发射窗口等）。

3.4.设计意图：以国家重大工程实例激发兴趣，明确本模块要解决的真实、复杂问题，建立学习的目标感和价值感。

5.环节二：动力学模型的重构与数学表述（40分钟）

1.6.教师活动：

1.2.7.回顾与批判：快速回顾二体问题的解析解（圆锥曲线），明确指出其“仅考虑中心天体引力”的强假设在真实深空任务中的局限性。提问：“如果考虑太阳、地球、火星同时的引力，运动方程会变得怎样？还能解析求解吗？”

2.3.8.第一性原理推导：带领学生在黑板上/电子白板上，从牛顿万有引力定律开始，逐步推导探测器的运动方程。首先在地心惯性系中，写出地球引力项。然后，逐步“添加”物理项：太阳引力项（引入日心位置矢量）、火星引力项（引入火心位置矢量）。最终得到一个形式简洁但内涵复杂的矢量微分方程：r

⃗

¨

=

−

μ

⊕

r

3

r

⃗

−

μ

⊙

(

r

⃗

−

r

⃗

⊙

∣

r

⃗

−

r

⃗

⊙

∣

3

+

r

⃗

⊙

r

⊙

3

)

−

μ

M

a

r

s

(

r

⃗

−

r

⃗

M

a

r

s

∣

r

⃗

−

r

⃗

M

a

r

s

∣

3

+

r

⃗

M

a

r

s

r

M

a

r

s

3

)

\ddot{\vec{r}}=-\frac{\mu_{\oplus}}{r^3}\vec{r}-\mu_{\odot}(\frac{\vec{r}-\vec{r}_{\odot}}{|\vec{r}-\vec{r}_{\odot}|^3}+\frac{\vec{r}_{\odot}}{r_{\odot}^3})-\mu_{Mars}(\frac{\vec{r}-\vec{r}_{Mars}}{|\vec{r}-\vec{r}_{Mars}|^3}+\frac{\vec{r}_{Mars}}{r_{Mars}^3})

r<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

¨=−r3μ⊕​​r<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

−μ⊙​(∣r<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

−r<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

⊙​∣3r<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

−r<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

⊙​​+r⊙3​r<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

⊙​​)−μMars​(∣r<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

−r<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

Mars​∣3r<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

−r<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

Mars​​+rMars3​r<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

Mars​​)。解释最后两项的推导源于在非惯性系中引入的间接项。

3.4.9.模型简化讨论：引导学生讨论：在从地球到火星的不同飞行阶段，哪些力占主导？是否可以分段采用不同的简化模型？例如，在地球影响球内以地心惯性系为主，在日心转移段以日心惯性系下的限制性三体问题（太阳-探测器-目标行星）为主。

4.5.10.化为标准ODE形式：强调数值求解只处理一阶方程组。定义状态向量Y

⃗

=

[

x

,

y

,

z

,

v

x

,

v

y

,

v

z

]

T

\vec{Y}=[x,y,z,v_x,v_y,v_z]^T

Y<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

=[x,y,z,vx​,vy​,vz​]T，将二阶微分方程化为6个一阶微分方程组成的系统：d

Y

⃗

d

t

=

f

(

t

,

Y

⃗

)

\frac{d\vec{Y}}{dt}=f(t,\vec{Y})

dtdY<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

​=f(t,Y<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

)。这里，f

f

f就是需要编程实现的“右函数”。

6.11.学生活动：跟随推导，记录关键步骤。理解从物理模型到数学方程每一步的假设和意义。尝试在JupyterNotebook中，用注释和Markdown单元格，将上述推导过程整理成文档，并初步定义状态向量和右函数的框架。

7.12.设计意图：这是本课的基石。通过“白板推导”，将隐藏在前置课程中的数学细节暴露出来，让学生亲手触摸到模型的复杂性，理解“数值求解”的必然性。将方程化为标准形式，是为后续的编程实现做直接准备。

13.环节三：数值积分算法原理精讲（30分钟）

1.14.教师活动：

1.2.15.从欧拉法切入：回顾最简单的显式欧拉法：Y

n

+

1

=

Y

n

+

h

f

(

t

n

,

Y

n

)

Y_{n+1}=Y_n+hf(t_n,Y_n)

Yn+1​=Yn​+hf(tn​,Yn​)。通过几何图示解释其“用当前斜率预测下一步”的思想。同时，通过一个简单例子（如指数衰减方程）演示其精度低、稳定性差的缺点。

2.3.16.引入龙格-库塔法的哲学：解释RK法的核心思想：在步长h

h

h内，通过计算多个“试探斜率”的加权平均，来获得比欧拉法更高阶的精度。详细推导经典四阶RK（RK4）的公式，并解释其四个斜率k

1

,

k

2

,

k

3

,

k

4

k_1,k_2,k_3,k_4

k1​,k2​,k3​,k4​的物理/几何意义：它们代表了在步长开始、中间点和结束点对微分方程右函数的采样。

3.4.17.算法对比与工程选择：对比RK4与多步法（如Adams）。用图表展示在相同精度下，RK4需要更多次函数调用，但启动简单、步长易变；多步法效率高，但需要历史信息，启动复杂。引导学生得出结论：在轨道动力学这种右函数计算较复杂、且可能需要频繁变步长（如近天体时）的场合，单步的RK家族方法是更稳健和方便的选择。介绍SciPy中solve_ivp

函数提供的多种RK和变步长算法。

5.18.学生活动：理解RK4公式的由来。在笔记本上，手动计算一个简单微分方程（如y

˙

=

−

y

\dot{y}=-y

y˙​=−y）的1-2步RK4解，与解析解对比，直观感受其精度。

6.19.设计意图：算法是连接数学方程和计算机代码的桥梁。不仅要学生“会用”RK4，更要理解其“为何有效”以及“为何适合轨道问题”。为后续自主编写和调用积分器打下坚实的理论基础。

第二篇章：计算实现与基础仿真（第3-4学时）

1.环节四：编程实践——构建轨道积分器（70分钟）

1.2.教师活动：

1.2.3.示范右函数编写：在Jupyter中现场编码。首先，定义天文常数（G,M

⊕

M_{\oplus}

M⊕​,M

⊙

M_{\odot}

M⊙​等）。然后，编写一个函数orbital_odes(t,state,masses,positions)

。该函数输入当前时间t和状态向量state，根据给定的中心天体质量和位置（对于日心系，positions是行星位置，需通过简单历书模型或调用库函数计算），计算引力加速度，并返回6维的状态导数Y

˙

\dot{Y}

Y˙。这是本节课最关键的代码。要详细讲解矢量运算如何通过NumPy的数组操作高效实现，避免循环。

2.3.4.实现RK4积分器：编写一个通用的rk4_integrate(func,y0,t_span,dt)

函数。强调其通用性，使其可用于任何ODE系统。讲解如何通过循环，利用RK4公式更新状态。同时，引入结果存储（将每一步的t和Y存入列表或数组）。

3.4.5.首个集成测试——地球圆轨道验证：设置初始条件为近地圆轨道（速度v

=

μ

⊕

/

r

v=\sqrt{\mu_{\oplus}/r}

v=μ⊕​/r<pathd="M263,681c0.7,0,18,39.7,52,119

c34,79.3,68.167,158.7,102.5,238c34.3,79.3,51.8,119.3,52.5,120

c340,-704.7,510.7,-1060.3,512,-1067

l0-0

c4.7,-7.3,11,-11,19,-11

H40000v40H1012.3

s-271.3,567,-271.3,567c-38.7,80.7,-84,175,-136,283c-52,108,-89.167,185.3,-111.5,232

c-22.3,46.7,-33.8,70.3,-34.5,71c-4.7,4.7,-12.3,7,-23,7s-12,-1,-12,-1

s-109,-253,-109,-253c-72.7,-168,-109.3,-252,-110,-252c-10.7,8,-22,16.7,-34,26

c-22,17.3,-33.3,26,-34,26s-26,-26,-26,-26s76,-59,76,-59s76,-60,76,-60z

M100180h400000v40h-400000z">

​）。在仅考虑地球中心引力（即二体问题）的右函数下，用自编的RK4积分器进行约10圈的积分。演示如何绘制轨迹图（2D和3D）。引导学生观察：理想的圆轨道积分后是否闭合？能量是否守恒？通过此测试，验证积分器和右函数的正确性。

4.5.6.步长影响探究：设置不同的积分步长（如100秒，1000秒），运行相同的仿真。展示步长过大导致的轨迹发散（数值不稳定）和能量漂移。引出“步长选择需远小于系统特征时间（轨道周期）”的经验法则。

6.7.学生活动：跟随编码。这是“手脑并用”的关键环节。学生在自己的电脑上，逐行复现或借鉴教师的代码，完成右函数和RK4积分器的编写。运行地球圆轨道测试，并尝试改变步长，观察结果，记录现象。遇到错误时，鼓励调试和同伴讨论。

7.8.设计意图：将理论转化为可运行的代码。通过“从零构建”一个积分器，学生能最深切地理解数值积分的每个细节。首个测试案例（二体圆轨道）简单且预期明确，能快速给予学生正向反馈，建立信心。步长实验则培养了其参数敏感性意识。

9.环节五：从二体到多体——引入摄动（50分钟）

1.10.教师活动：

1.2.11.升级右函数：在之前的右函数基础上，增加地球非球形引力J2项的近似加速度表达式a

⃗

J

2

=

.

.

.

\vec{a}_{J2}=...

a<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-4.6670-9.167-1.833-13.5-5.5S337184337178c0-12.66715.667-32.33347-59

H213l-171-1c-8.667-6-13-12.333-13-190-4.6674.333-11.33313-20h359

c-16-25.333-24-45-24-59z">

J2​=...。解释J2项如何导致轨道进动和交点线旋转。

2.3.12.近地轨道摄动仿真：对一个高度约400km的近地圆轨道，先进行纯二体积分，再进行包含J2项的积分。将两者结果叠加绘图，直观展示J2项引起的轨道面旋转（轨道法向量方向变化）和近地点幅角进动。

3.4.13.引入太阳引力摄动：演示如何修改右函数，加入第三体（太阳）的引力摄动项。采用简化模型：假设地球和太阳在惯性空间中的位置已知（可用近似圆周运动模型计算）。对一个地球同步轨道进行长时间（如30天）积分，比较仅二体、二体+J2、二体+J2+太阳引力三种模型下的轨道差异。

4.5.14.可视化与分析工具：教授学生绘制轨道根数（半长轴a、偏心率e、轨道倾角i等）随时间变化的曲线。这需要实时将积分得到的直角坐标状态向量转换为轨道根数。提供转换函数或引导学生调用现有库（如astropy.coordinates

）。

6.15.学生活动：修改自己的右函数，加入J2项。运行近地轨道摄动仿真，观察并记录J2效应的可视化结果。尝试计算轨道根数的变化，定量分析摄动影响。思考：对于低轨遥感卫星，哪个摄动力是必须考虑的？对于高轨通信卫星呢？

7.16.设计意图：引导学生将模型从理想走向现实。通过对比仿真，让学生亲眼看到不同物理因素对轨道的具体影响，理解模型复杂化的必要性。轨道根数时间序列图是一种更专业、更深入的分析工具，引导学生从“看轨迹”深入到“看参数演化”。

第三篇章：项目实践——地火转移轨道设计探索（第5-7学时）

1.环节六：项目启动与初始猜测生成（40分钟）

1.2.教师活动：

1.2.3.问题重述与约束明确：正式发布项目任务书：“使用数值方法，设计一条从地球停泊轨道（高度约400km圆轨道）到火星影响球边界的转移轨道，目标是最小化所需的速度增量（Δ

V

\DeltaV

ΔV），转移时间约束在200-300天。”提供地球和火星在当前历元下的平均轨道半径和周期。

2.3.4.引入概念与简化：讲解霍曼转移作为理论最优解的前提（共面圆轨道、脉冲推力）。引导学生计算日心系下的理论霍曼转移参数（转移轨道半长轴、所需日心速度增量）。

3.4.5.从日心到地心：讲解如何将日心转移轨道的出发条件，转换为地心惯性系中探测器的逃逸超速（v

∞

v_{\infty}

v∞​）。提供一个简化的初始猜测：在地球停泊轨道上，施加一个脉冲，使探测器获得一个相对于地球的逃逸速度矢量，其大小和方向由简化模型估算得出。这个矢量与地球公转速度矢量的和，应近似等于日心霍曼转移所需的出发速度。

4.5.6.提供初始猜测代码框架：给出一个函数，输入期望的出发日期（简化处理），输出探测器在地心惯性系中的初始位置和速度矢量（猜测值）。

6.7.学生活动：理解任务要求。记录霍曼转移的理论计算结果。运行教师提供的初始猜测生成代码，获得自己项目的第一个“设计点”。讨论这个猜测的潜在问题（如未考虑行星引力辅助、未精确对齐火星位置）。

7.8.设计意图：将大项目分解为可操作的步骤。提供初始猜测是降低项目门槛的关键，避免学生陷入“无从下手”的困境。同时，通过理论计算与数值探索的对比，埋下后续优化的伏笔。

9.环节七：轨道积分与“命中”判断（80分钟）

1.10.教师活动：

1.2.11.构建仿真系统：指导学生建立项目的主仿真循环。动力学模型采用日心惯性系下的限制性三体问题（太阳-探测器-目标行星），或更精确的简化多体模型。右函数需包含太阳、地球、火星的引力（当探测器靠近相应天体时，该天体引力切换为主项）。

2.3.12.实现积分与终止条件：使用更高效的积分器（如SciPy的solve_ivp

，方法选用DOP853

等高阶RK），从初始猜测开始积分。设置终止条件：当探测器与火星的距离小于火星影响球半径（约0.5百万公里）时，或总积分时间超过400天时，停止积分。

3.4.13.“脱靶”分析与可视化：首次积分很可能无法“命中”火星。教师演示如何可视化结果：在日心黄道面绘制地球轨道、火星轨道、探测器轨迹以及积分终止时火星的位置。引导学生分析“脱靶”情况：是到达时间不对（探测器跑到火星轨道时，火星不在那里），还是轨道能量不对（转移轨道半长轴误差）？

4.5.14.引入误差度量：定义“脱靶”的量化指标。最直接的是在目标到达时间（由霍曼转移周期估算），计算探测器与火星的位置矢量差Δ

r

⃗

\Delta\vec{r}

Δr<pathd="M37720c0-5.3331.833-105.5-14S39103970c4.66708.6671.667125

3.3332.6676.667910196.66724.66720.33343.66741577.3334.66711

10.667111806-110-312s-6.6675-149c-28.66714.667-53.66735.667-7563

-1.3331.333-3.1673.5-5.56.5s-44.833-55.5c-1.667-2.51.333-4.52s-4.3331

-71c-