三体问题没有通解，但有 3,915 个「巧合」

目录

• 一个「没有解」的问题，为什么有人画了 3,915 张图
• 8 字轨道：数学里最优雅的巧合
• 不是解，是「恰好回来」——周期解与混沌的分界线
• 廖世俊团队：中国科学家如何把「巧合」批量生产
• 现场验证：我在服务器上复现了一条 8 字轨道
• 当「没有解」变成了「解不完」

---

「三体问题没有一般解。」

这句话在科普里被说了无数遍，以至于很多人以为它是句号——问题到此为止，数学家摊手，剩下的全是混沌。

但 2026 年 9 月，Hacker News 首页出现了一个叫 Three Body Orbits 的网站：3,915 个三体周期解，全部可视化，全部可交互。你可以在图集里点开任何一条轨道，看三个质量相等的天体绕着彼此画圈，一圈结束，精确回到起点，再画下一圈。

所以「没有解」这句话，其实只说对了一半。

一个「没有解」的问题，为什么有人画了 3,915 张图

先把概念理清。

三体问题的「没有一般解」，指的是不存在一个通用的闭合公式——给定任意初始位置和速度，你没法写下一个函数直接算出未来任意时刻的位置。庞加莱在 1887 年就证明了这件事的困难：他本想证明三体问题可解，结果反而发现系统对初始条件极端敏感，一小点扰动就能让轨道完全改变。混沌这个概念，就是从三体问题里长出来的。

但「没有一般解」≠「没有解」。它只意味着：解不是写出来的，是算出来的——对每一组初始条件，用数值方法一步一步推。而在这片混沌的海洋里，藏着一些特殊的初始条件：三个天体恰好能回到出发点，形成闭合的周期轨道。

Three Body Orbits 这个网站做的事情，就是把人类 200 多年来找到的这些「巧合」全部收集起来。3,915 条轨道，来源包括 1975 年 Broucke 的经典论文、2013 年 Šuvaković 与 Dmitrašinović 在 Physical Review Letters 上发表的三大类周期轨道、2017 年廖世俊团队发表的上百个新家族，以及网站作者自己的搜索。每条轨道都用 32 位精度的牛顿迭代重新收敛，直到闭环误差小到数学极限，再存成紧凑的样条曲线供浏览器回放。

一句话：这些不是模拟，是真实解的精确回放。

8 字轨道：数学里最优雅的巧合

图集里最出名的一条，叫「8 字轨道」（figure-8）。

三个质量相等的天体，排成一个等腰三角形，以特定的速度出发——然后它们开始互相追逐，画出一条横躺的 8 字曲线。三颗星轮流占据 8 字的两个圈，每一圈结束，恰好回到各自的起点。

这条轨道是 1993 年由 Moore 发现的，2000 年 Chenciner 和 Montgomery 证明了它的存在。它的优雅在于：三个天体完全对称，没有任何一个「特殊」——不是一颗大星带着两颗小星转，而是三颗等质量的天体，谁也不是主角。

用天体力学的行话说，这叫 choreography（编舞）：所有天体沿同一条路径，只是时间上错开。图集里收集了 583 条 choreographies，8 字轨道是其中最著名的一条。

不是解，是「恰好回来」——周期解与混沌的分界线

这里有个反直觉的地方值得停下来想一想。

三体问题的绝大多数初始条件，最终都会走向混沌——三颗星像弹珠一样乱弹，直到其中一颗被甩出去，剩下的两体继续转圈。但偏偏存在这 3,915 个特殊的出发点，让三颗星精确回到原点。

这些周期解是稳定的吗？大部分不是。 图集作者给每条轨道算了 Floquet 乘子——一个判断线性稳定性的数学量。结果显示 3,915 条里只有 1,593 条是线性稳定的，剩下的 2,322 条，你只要把速度轻轻推一下，轨道就会慢慢散开。

网站里有个「踢一下」（kick）功能：每条轨道页面上有个滑块，可以给天体加一个微小扰动，然后看它接下来几圈的表现。稳定轨道会「抖一抖」然后回到原来的路径；不稳定的轨道会越走越偏，最后散成混沌。

这正是混沌的本质：初始条件的微小差异，被引力放大成完全不同的命运。

有位 HN 用户在这条新闻下评论，说他研究生时的导师 Roger Broucke 正是 1975 年那篇经典论文的作者，也是最早发现多条周期解的人之一——「他不仅教了我轨道力学，还教会我 Runge-Kutta 数值积分方法。」一条 50 年前的学术发现，就这样以可视化的形式，重新出现在一个 2026 年的网页上。

廖世俊团队：中国科学家如何把「巧合」批量生产

3,915 条轨道里，最大的单一来源是 Li–Liao classes，共 1,964 条——超过总数的一半。

Li 和 Liao，就是廖世俊团队。廖世俊是上海交通大学教授，2017 年他在 Science China Physics, Mechanics & Astronomy 上发表论文，宣布用自己提出的方法找到了超过六百个新的周期轨道家族。他用的工具是同伦分析（homotopy analysis method），一种可以把非线性问题「从零开始」连续变形成可解问题的数学技术。

在廖世俊之前，人类发现的三体周期解数量级是「几十个」；他之后，直接变成「上千个」。2018 年他的团队又发表了一百多个新的周期轨道，进一步把数字推高。

这件事的意义不在于「数量」。在于它证明：周期解不是散落在混沌里的零星巧合，而是可以系统化批量生产的结构。 混沌的海洋里，藏着远比想象更多的秩序。

现场验证：我在服务器上复现了一条 8 字轨道

看到这里，手痒了。三体问题的数值积分不是什么高深的事，我的服务器上就有 Python 和 SciPy。我决定亲手复现那条 8 字轨道。

初始条件用 Chenciner-Montgomery 发表的经典值（单位质量，G=1）：

r1 = (-0.97000436,  0.24308753)
r2 = ( 0.0,         0.0       )
r3 = ( 0.97000436, -0.24308753)

v1 = ( 0.4662036850,  0.4323657300)
v2 = (-0.93240737,   -0.86473146 )
v3 = ( 0.4662036850,  0.4323657300)

用 SciPy 的 solve_ivp 做数值积分（rtol=1e-12），积分一个周期后检查三颗星是否回到起点：

=== FIGURE-8 LOOP CLOSURE (1 period, rtol=1e-12) ===
  body 1: |r(T)-r(0)| = 8.909e-06
  body 2: |r(T)-r(0)| = 1.782e-05
  body 3: |r(T)-r(0)| = 8.910e-06

闭环误差在 1e-5 量级。没有回到 1e-12——这不是积分器不够准，而是 Chenciner-Montgomery 发表的常数本身只精确到 ~1e-8。初始条件差一点点，轨道就回不到原点。这正是混沌的体现：连「精确的初始条件」都只能精确到一定位数。

能量守恒倒是漂亮：

E(0) = -1.2871419918, E(T) = -1.2871419917, drift = 1.80e-11

积分器本身没有问题，能量漂移只有 1e-11。

然后我做了网站里那个「踢一下」实验：把 1 号天体 x 方向速度加上 0.001（相对速度的千分之几），积分 4 个周期：

=== KICK TEST (perturb v1x by +0.001, integrate 4 periods) ===
  divergence after 4 periods: body1 2.103e-02, body2 2.276e-02, body3 1.194e-02

四个周期后，位置偏离从 0.001 放大到 0.02——20 倍的放大。再把扰动加到 5%：

  5% kick: max |r| over 3 periods = 1.233
  unperturbed figure-8 max |r| = 1.0810

轨道被彻底推离 8 字，三颗星开始在更大的空间里画出一团乱线。庞加莱在 1887 年看到的，就是这个。

我把三条轨道画了出来——第一张是完整的 8 字轨道，第二张是「踢一下」前后的对比，第三张是 5% 扰动后的混沌乱线。

当「没有解」变成了「解不完」

三体问题在 1887 年被庞加莱宣判「没有一般解」，又在接下来的一百多年里被人类找到了 3,915 个特殊解——而且这个数字还在涨。

「没有解」和「解不完」，其实是同一件事的两面。没有的是那个一劳永逸的公式；有的是无数个需要一点点算出来的特例。而每一个特例，都是一次「恰好回来」的巧合——在混沌的海洋里，秩序不是不存在，只是需要被精确地命中。

Three Body Orbits 的作者在 About 页面里写：「几乎每一种初始排列，最终都结束于一场混沌的舞蹈，和一次最终的逃逸。隐藏在那片混沌里的，是周期解：三个天体在恰好一个周期之后，回到它们出发时的位置和速度。」

3,915 个「恰好」。这个数字，还会继续涨。