1. 课程导入与多体问题概述

当我们抬头观察夜空，行星沿各自轨道绕太阳运行，卫星环绕行星，双星彼此绕转，这些现象背后都遵循同一个物理规律——万有引力。然而，一旦系统中存在三个或更多相互施加引力的天体，问题便从"可解"跃入"复杂"。本课程《多体引力运动与轨道计算》正是围绕这一核心矛盾展开：如何描述、计算并预测多个引力天体的运动。

什么是多体问题

多体问题（N-body problem）指研究三个及以上天体在牛顿万有引力作用下的运动规律。对任意两个质点，引力大小为

F = G · m₁ · m₂ / r²

其中 G 为引力常数，m₁、m₂ 为两体质量，r 为间距。两体问题可以解析求解，轨道是圆锥曲线；但三体及以上一般不存在封闭解析解，必须借助数值积分。这正是多体问题既经典又现代的原因：它连接了天体力学、航天动力学与科学计算。

历史脉络

1687 年牛顿在《自然哲学的数学原理》中给出万有引力定律，两体问题迎刃而解。但三体问题很快成为难题：欧拉与拉格朗日找到了若干特殊解（如共线解与等边三角形解），庞加莱在 19 世纪末证明三体问题一般不可积，并由此开创混沌研究。20 世纪以来，计算机使数值求解成为主流，从阿波罗登月的轨道设计到今天的系外行星动力学，多体方法贯穿始终。

典型应用场景

行星系统：太阳系中行星之间的相互摄动会累积成轨道长期演化，例如木星对近地小行星的引力扰动直接影响其撞击风险。

卫星编队：多颗卫星在相近轨道上保持相对构型，需要精确计算彼此引力与地球非球形引力，用于分布式观测与干涉测量。

系外行星：通过恒星视向速度或凌星时刻变化反推行星质量与轨道，多体模型是解释轨道共振与稳定性的基础。

学习路径与工具准备

本课程按"物理建模—数值方法—轨道计算—工程应用"推进：先建立引力与运动方程，再引入数值积分器，随后处理轨道要素、摄动与稳定性，最后落到编队与系外行星等场景。工具上建议准备 Python（NumPy、SciPy、Matplotlib）或 MATLAB，并熟悉基本的天文单位与坐标系约定。

下图概括了本节的核心知识结构：从两体可解到多体不可积的分岔，以及课程各模块之间的依赖关系。