繁星无法超越——三体问题溯源

2024年4月11日发(作者：六年级上册数学试卷沁阳)

第

40

卷第

1

期

2021

年

1

月

大学物理

COLLEGE

PHYSICS

Vol.40

No.1

Jan.

2021

繁星无法超越

一一

三体问题溯源

邵瀚雍

(北京师范大学物理学系

，

北京

100875

)

摘要

：

三体问题是一个困扰学术界多年的难题

.

自牛顿在

《

自然哲学的数学原理

》

中首次定义三体问题开始，

人们已经在

这个问题上努力了

300

年之久

.

欧拉和拉格朗日先后研究了一些较为简单的三体系统问题

，

而庞加莱首先考虑了限制性的三

体问题

，

但发现其轨迹十分复杂而无法描绘

.

他转而证明了三体问题不存在解析解

，

并着手创立了一门崭新的学科

一

非线

性动力学

.

充满生命力的非线性动力学或能为人们解决这个问题提供不同的思路

.

关键词

：三体问题;庞加莱

；

非线性动力学

中图分类号

：

O4-1

文献标识码

：

A

文章编号

：

1000-

0712(

2021

)01-0060-

06

【

DOI

】

10.16854/.1000-0712.200077

在科幻小说

《

三体

》

中

，

作者构建了一个特殊的

恒星系统:某行星系中有三颗质量不等的恒星

，

在彼

问题

，

也就是上述

N

=

3

时的三体问题⑷

-

颇为可惜的是

，

即使聪明如牛顿

，

至临终前也没

有给出任何进一步的解答

.

从物理上讲

，

这个问题

此的引力下相互作用

，

而该行星系的文明则位于

3

颗恒星附近的一颗行星上

.

因为恒星无规则的运

动

，

这个行星上的文明饱受摧残

•

一开始他们试图

预测恒星运动的规律

，

但每次的尝试都以失败告终

.

并不复杂

，

但是它对数学的要求很高

.

在牛顿生活

的

17

世纪

，

微积分学才刚刚被创立不久

，

当时人们

的数学水平远远不能解决这个复杂的物理问题

.

于

是,这个悬而未决的问题又被拖了整整

100

年之久

，

等待着下一位数学家来解决

.

1772

年

，

拉格朗日

(

Lagrange

)

和他的导师欧拉

他们转而决定逃离行星系

，

奔向地球

•

这个难题在

1900

年的第一次世界数学家大会

上就已经被提出了

.

著名数学家希尔伯特

(Hilbert

)

在会上提出了

23

个困难的数学问题

；

与

(Euler

)

接下了研究三体问题的火炬

.

他们首先开始

此同时

，

他特别谈到了两个经典的完美数学问题

，

其一是我们熟知的费马猜想

；

其二便是

N

体问题

.

考虑研究限制性的三体问题

，

在这个问题中

，

一个天

体的质量被忽略不计

，

而只受到其它两个大天体的

引力作用而运动

.

作为数学分析创始人之一的拉格

这个问题要求人们给出在已知质量和初始时刻坐

标

、

速度的

N

个质点在彼此引力作用下的后续运

朗日

，

他通过一系列严密的计算

，

求出了这个限制性

动规律

.

问题的后两个特解(前

3

个为欧拉所求得)

.

这个

5

个解正是我们熟悉的拉格朗日点

[2

]

-

1

三体溯源

：

数学家的尝试

事实上

，

三体问题有着深刻的数学和物理背景

.

作为经典力学的主要创立者

，

牛顿在研究万有引力

定律时

，

就已经基本解决了最简单的两体问题

，

并预

拉格朗日点

，

又称平动点

，

是指当一个小物体在

两个大物体的引力作用下

，

可以在空间中找到一点

，

在该点处

，

小物体相对于两大物体基本保持静止

.

拉格朗日点的发现

，

对于三体问题的求解有着巨大

的意义:人们至少看到了希望的曙光

.

相应的

，

人们

测了双星系统运动规律

.

合理外推是科学探索的核心理念之一

：

当人们

在天文发现中也证实了拉格朗日点的确存在

.

在图

1

的

5

个拉格朗日点中

，

有两个拥有良好

知道一元二次方程的通解

，

就会自然而然地思考一

元三次方程的通解

；

当人们知道两体问题的解时

，也

会思考三体乃至

N

体问题的通解

，

而日

-

地

-

月系统

的稳定性,人们可以发射航天器稳定在拉格朗日点

，

这样航天器就可以利用两个天体的引力达到相对静

止的状态

，

比如詹姆斯

•

韦伯空间望远镜

(JWST

)就

正是我们最熟悉的三体系统之一

.

作为外延

，

牛顿

首先探讨了

3

个天体在引力作用下的运动系统演化

处在日

-

地系统的第二拉格朗日点

.

除了宏观三体

收稿日期

：

2020-03-12

；

修回日期

：

2020-04-21

作者简介

：

邵瀚雍

(

2000

—

)

，

男，

四川德阳人

，

北京师范大学物理学系

2018

级本科生

.

第

1

期

邵瀚雍:繁星无法超越

—

—

三体问题溯源

61

系统的引力场内存在拉格朗日点外

，

也有人发现类

似的点电荷电场中也存在类拉格朗日点捋

I

图

1

三体系统中的拉格朗日点

自拉格朗日点发现以后

，

人们对于三体问题的

研究又停滞了近

100

年之久

，

直到某次机缘巧合

，

才

让这个被遗忘已久的问题又重新回到了人们的视

野中

.

1889

年

，

瑞典与挪威的国王奥斯卡二世

(

Oscar

n

)

要准备自己的

60

岁大寿

.

而奥斯卡二世在大学

期间对数学很感兴趣

，

对科学研究也很支持

，

于是资

助创办了

《

数学学报》

杂志

.

为了这次贺寿

，

身为

《

数

学学报

》

主编的米塔格

・

列夫勒建议

，

设立一个竞

赛

，

让人研究太阳系是否稳定

•

而这正是一个

N

体

问题的变式

•

此前的天体力学家虽然有过繁杂的计

算

，

但是仍然无法解决小分母问题

•

当年拉普拉斯

通过计算证明了太阳系的稳定性

，

但是他只能后推

900

年

，

对于几十亿年龄的太阳系来说

，

这不过是弹

指一挥间

，

难以真正回答这个问题

•

竞赛公告的第一段是

：

“

具有任意多个质点的

系统

，

其中任何两点间的作用力满足牛顿定律

，

在任

意两个点不发生碰撞的情况下

，

试给出每个点的坐

标

，

这个坐标可以以时间的某个已知函数作为变量

的级数表示

，并且对于所有的取值，

该级数是一致收

敛的

.

”

在这个时候

，

一位在数学领域造诣极高的数学

家接下了火炬

，

他正是庞加莱

(

Poincare

)

•

他的研究

领域涉及数学

、

物理和天体力学

，

并被公认为是

19-

20

世纪之交的领袖数学家

.

在这次公开悬赏之前

，

庞加莱已经凭借自己在

自守函数论和微分方程等方面的工作而被人们公认

为世界一流数学家了

.1883

年

，

29

岁的庞加莱发表

了第一篇天体力学方面的论文;一年之后

，

他又发表

了第二篇关于天体力学的论文

.

奥斯卡二世的这个

问题进一步激发了庞加莱对这个天体力学问题的

图

2

庞加莱像

兴趣

.

1885

年

，

庞加莱开始集中精力研究三体问题,

此时距离竞赛截至日期还有大约

3

年时间

.

庞加莱

先采用了希尔研究月球时所采用的模型

，

但经过不

断坚持与搁置的循环尝试后

，

他终于理清了方向

•

1887

年

，

庞加莱发现这个问题不能用牛顿

、

拉格朗

日与拉普拉斯等人所采用的小参数展开等传统方法

解决

，

而必须改换一个新的思路

•

庞加莱首先考虑了较为理想的情况

，

也就是我

们熟知的限制性三体问题

：

一个尘埃漂浮在两个巨

大天体形成的体系里

，

应该如何预测尘埃的轨迹

•

在计算机尚未出现的时代里

，

他依靠惊人的计算能

力

，

绘出了尘埃前期的运动轨迹

•

当他试图在纸上

画出那些轨迹的时候

，

他写道

：

这难以画出来的图

形的复杂性令我震惊

！

”

.

图

3

复杂的尘埃轨迹示意图

随后

，

庞加莱另辟蹊径

，

独自发展出一套方法

：

微分方程定性理论

.

他开始尝试运用渐近展开与积

62

大学物理

第

40

卷

分不变性的方法

，

定性研究小尘埃的轨迹⑸

，

并且

深入研究了小尘埃在同宿轨道和异宿轨道

（

相当于

奇点

）

附近的行为

•

但是在研究过程中

，

庞加莱发

现

，

尘埃的轨迹几乎不可能被预测

•

经过整整三年

的努力

，

他认为这个问题无法完全解决

，

并随后把自

己的研究成果寄到论文评审委员会

，

还在论文开头

写了

一句\":

u

Nunquam

praescriptis

transibunt

sidera

fines

”

（

繁星无法超越

）

.

尽管庞加莱并没有完全满足奥斯卡二世给出的

要求

，

也没有完全解决

N

体问题

，

但他长达

160

页

的论文仍然令评审团的三位数学大师兴奋无比

•

他

们认为庞加莱对三体问题的研究取得了重大突破

：

太阳系的相对稳定得到确认

•

魏尔斯特拉斯

（

Weierstrass

）

在给杂志主编的信

中写道

：

请告诉您的国王

，

这个工作虽然不能被视

为对所求问题的完整解答

，

但是它的重要性使得它

的出现将标志着天体力学的一个新时代的诞生

.

因

此

，

陛下预期的公开竞赛的目的

，

可以认为已经达到

了

”

.

于是国王高兴地把奥斯卡奖

------

枚金质奖

章和

2500

瑞典克朗

，

授予了庞加莱

•

3

夜幕降临

：

三体问题无解

对于三体问题

，

人们希望能找到一个长周期

（

哪怕这个周期非常长

）

或者是一个解析解

（

收敛的

级数解

）

来描绘三体运动的规律

•

但是

，

庞加莱证明

了三体问题的解析解是不可能存在的

•

这个由简单的代数表达式和积分给出的三体运

动微分方程没有通常意义上的解析解

，

并且

，

除特殊

情况外

，

3

个物体的运动是不重复的

•

在一般三体问

题中

，

每一个天体在其他两个天体的万有引力作用

下的运动方程都可以表示成

3

个二阶的常微分方程

（

如图

4

所示

）

，

或

6

个一阶的常微分方程

•

因此

，

一

般三体问题的运动方程为

18

阶方程

，

必须通过

18

个积分才能得到完全解

•

然而

，

目前还只能得到三

体问题的

16

个初始积分,还不能解决三体问题

\'

图

4

3

个二阶常微分方程

其实

，

早在

1887

年

，

勃隆斯

（

Bruns

）

就证明了

三体问题中质点的速度和位置无法用函数定量表

达

•

几年后

，

庞加莱进一步证明了位置

、

速度和质量

比的关系无法用定量的解析式表达

•

除了解析解不存在之外

，

数值解的情况也不容

乐观

•

为了得到具有意义的数值解

，

计算机往往需

要进行长时间的工作

•

早在

20

年前

，

就已经有对于

N

体问题进行数值解的研究了

，

但是效果不佳

*

］

•

并且

，

即使可以通过现代的量子计算方式极大

加快计算进程

，

也难以避免三体系统的一个核心问

题:误差

•

在进行数值分析的过程中

，

计算机往往会

近似估计相对的高阶小量

，

或因为计算机的数字尺

度限制

，

无法进行完整的数据迭代

•

短时间内也许

看不出什么区别

，

但是随着时间的积累

，

数值计算所

得到的结论可能与真实情况差别极大

•

除了宏观世界中存在的天体引力系统外

，

人们

也研究过微观世界中氮原子核与两个电子根据库仑

定律相互作用的三体系统

，

不过依然不能精确求解

出氮原子核的运动规律

［

9

\'

不过在最近的科学研究中

，

借助于计算机的高

速发展

，

已经有许多文献找到了不少

N

体问题运动

的特解

，

如图

5

和图

6

，

这些解十分优美

［

⑼.

第

1

期

邵瀚雍:繁星无法超越

—

—

三体问题溯源

63

4

夜空星光

：

混沌理论创立

在庞加莱证明了三体问题的无解性后

，

他并没

有停下这个工作

.

在研究这个问题的时候

，

他第一

个敏锐的发现了这种运动规律不循环的

、

对初始条

件敏感的行为

，

并着手研究这类型的问题

.

这类问

题被他统称为混沌

（

Chaos

）

.

庞加莱首先考虑了二维平面内的动力学问题,

并提出了一个重要的定理

，

即庞加莱一本迪克松定

理

[

1

]

•

该定理由庞加莱提出

，

并由本迪克松给出了

完整的证明过程

•

如果一个二维动力学系统内运动连续可微

，

且

动力学流在一个区域内封闭

，

而区域内也没有定点

，

则该区域的动力学流将收敛于一条闭合轨道

•

这个定理指出

，

在有限的二维系统内

，

物体的运

动形式只有两种

，

要么归于平衡态

，

要么满足周期运

动

.

换言之

，

庞加莱指出了二维系统中并不存在混

沌行为

•

图

7

是在进行节拍器同步过程中自由度为

1.5

的相平面轨迹

，

可以明显看出有一个定点

（

平衡

态

）

和一个较弱的闭合轨道

•

XY

Graph

Plot

°

400

300

200

100

亠

0

-100

-200

-300

-400

-400

-200

0

200

400

x

图

7

节拍器同步过程中的相平面轨迹图

不过

，

由于庞加莱所处时代的数学限制

，

他的研

究只是指出了运动的非线性方程所造成的动力系统

的非线性

.

尽管如此

，

他对于三体问题的研究也为

后者研究混沌系统打下来理论基础

•

在庞加莱之后的很长一段时间

，

非线性动力学

还是停留在非线性振子等比较基础的问题上

•

直到

20

世纪

60

年代初

，

美国的气象学家洛伦兹

（

Lorenz

）

在使用计算机进行气象问题研究时

，

发现了

“确定

性非周期流

”

•

他的研究指出

，

气候对初始条件的变

化极为敏感

，

且不同的条件会导致不同的演化轨迹

，

从而大大增加系统的不确定性

.

他将之通俗的比喻

为

“

蝴蝶效应

”

12

]

.

与此同时

，

1971

年

，

人们开始研究湍流问题的

纳维

-

斯托克斯方程

P

d

V

2

d

t

=

pg

一

Vp

+

“V

V

此时发现湍流解的形式也存在奇异吸引子的情

况⑴

]

，

这与之前庞加莱的研究结果很相近

，

而湍流

问题也正是一种困扰人们多年的混沌系统

•

随着混沌

、

分形和分岔理论的不断完善

，

多数学

者认为

，

混沌系统内有

3

个重要的特殊情况

[

4

]

:

1

）

系统对初始条件极其敏感

，

哪怕是最初状态

的轻微变化或是微小的误差都能大幅度改变结果的

形态

，

这直接导致混沌系统无法被预测

•

如图

8

，

这

是一个两面弧形的球桌

，

我们把球的初始方向从虚

线改变为实线

，

可以发现其运动轨迹完全不同

•

图

8

小球运动轨迹示意图

2

）

系统内部存在分形的情况

.

分形意味着系统

在不同的标度下有着自相似的性质

，

并且系统可以

产生具有结构和规则的难以发现的有序模式

.

如图

9

所示

，

这是限制性三体问题中尘埃在

10

个不同的初

始条件下

，

所演化的

10

条不同轨道的分形结果

，

但

这个看似无规则的轨迹有着一定的自相似性

.

图

9

尘埃运动轨迹的不同分形结果

64

大学物理

第

40

卷

3

）

系统内存在奇异吸引子

•

吸引子

，

是系统被

吸引并且最终固定为某一状态的性质

•

—

般来说,

吸引子分成三种

：

点吸引子

、

极限环吸引子和混沌吸

引子

•

前两者一般是收敛的

，

而奇异吸引子则不同

，

它可以使系统偏离收敛性的吸引子

，

而改变系统的

状态

，

增加系统的不可预测性

•

图

10

是由

Python

绘

出的洛伦兹

（

奇异

）

吸引子的示意图

•

非线性还表现在许多别的方面

，

除了经典的混

沌情况之外

，

系统也可能会出现多稳态的情况

，

这表

示系统有可能在两个或多个不同的互斥状态之间不

断切换

，

这也是非线性的表现之一

.

除此之外

，

非周

期振荡

、

孤立波等现象也不断被发现

•

随着人们对非线性方程研究的不断加深

，

非线

性系统正在作为一门崭新的学科指引人们前进的道

路

•

有的学者认为

，

混沌理论是继相对论和量子力

学之后本世纪的第三次科学革命

.

由此而生的非线

性物理也有着许多的分支学科

•

这其中包括了非线

性动力学

、

非线性声学

、

非线性光学和量子混沌学等

许多新兴学科

.

6

结语

在人类为了三体问题而奋斗的几百年里

，

虽然

结果是让人有些失望的

，

但是中间的过程是精彩而

又值得铭记的

•

费马猜想就被叫做

“

会生金蛋的

鸡

”

，

其原因正在于人们在证明的过程中诞生的种

种珍贵的思想和方法

•

而作为另一只下金蛋的

“

鸡

”

，

三体问题已经为

人类生下了不少珍贵的方法和理论

•

时至今日

，

人

们对于三体问题的兴趣仍然没有衰减

•

2013

年

，

贝尔格莱德物理研究所的物理学家发

现了等质量

、

零角动量的三体问题的

13

族特解

，

正

式将特解扩展为

16

族

口

5

1

2017

年

，

上海交通大学

的廖世俊研究小组利用超级计算机

，

找到了三体问

题的

600

多个新的周期解家族

•

如图

11

所示的几个

周期解

口

6

\'

1.01.0

亠

0.0

亠

0.0

丄?

1.5

-1.0

-0.5

0.0

0.5

1.0

1.5

0.8

------------

些

--------

°-

4

\'

0.5

0.0

-1.0

-0.5

0.0

0.5

1.0

1.5

©

0.5

0.5

0.0

0.0

图

11

新发现的一些周期解家族

除了对周期解的研究之外

，

在探寻数值解的道

路上

，

为了达到高精度和长时间的预测

，

有学者已经

尝试改用数值积分的方式

，

通过神经网络来进行机

器学习

，

以达到预测轨迹的目的⑴

：

•

参考文献

[1

]

Newton

I.

Mathematical

principles

of

natural

philosophy

[M

]

.

A.

Strahan

，

1802

:

205-303.

[2

]

Sosnitskii

S

P.

On

the

stability

of

triangular

Lagrangian

points

in

the

restricted

three-body

problem

[

J

]

.

The

As­

tronomical

Journal

，

2007

，

135

(

1

)

:

187.

[3

]

傅慎明

.

经典电磁学范畴内的圆型限制性三体问题

[

J

]

.

大学物理

，

2013

，

32

(5

)

:

52-52.

[4

]

Diacu

F.

The

solution

of

the

n-body

problem[R

]

.

1995.

[

5

]

Poincare

H.

New

methods

of

celestial

mechanics

[

M

]

.

Springer

Science

&

Business

Media

，

1992.

[6

]

Poincare

H.

Sur

le

probleme

des

trois

corps

et

les

equations

de

la

dynamique

[

J

]

.

Acta

mathematica

，

1890

，

13

(1

)

:

A3-A270.

[7

]

卢圣治

.

理论力学基本教程

[M

]

.

北京师范大学出版

社

，

2004

:

280-283.

[8

]

Greengard

L.

The

numerical

solution

of

the

n

-

body

problem

[

J

]

.

Computers

in

physics

，

1990

，

4

(

2

)

:

142-152.

第

1

期

邵瀚雍:繁星无法超越

—

—

三体问题溯源

65

[9

]

Griffiths

D

J

，

Schroeter

D

F.

Introduction

to

quantum

me­

chanics

M

]

.

Cambridge

University

Press

，

2018

：

311.

biology

，

chemistry

，

and

engineering

[

M

]

.

CRC

press

，

2018.

[

15

]

Suvakov

M

，

Dmitrasinovic

V.

Three

classes

of

Newtonian

[

10

]

Montgomery

R.

A

new

solution

to

the

three

-

body

problem

[

J

]

.

differential

equations

，

2001

，

1001

:

d2xi.

[

11

]

Hirsch

M

W

，

Smale

S

，

Devaney

R

L.

Differential

equa-

three-body

planar

periodic

orbits

[

J

]

.

Physical

review

letters

，

2013

，

110

(

11

):

114301.

[

16

]

Li

X

M

，

Liao

S

J.

More

than

six

hundred

new

families

of

tions

，

dynamical

systems

，

and

an

introduction

to

chaos

[

M

]

.

Academic

press

，

2012.

[

12

]

Lorenz

E

N.

Deterministic

nonperiodic

flow

[

J

]

.

Journal

Newtonian

periodic

planar

collisionless

three

-

body

orbits

[J

]

.

SCIENCE

CHINA

Physics

，

Mechanics

&

Astronomy

，

2017

，

60

(

12

):

129511.

[

17

]

Breen

P

G

，

Foley

C

N

，

Boekholt

T

，

et

al.

Newton

vs

the

of

the

atmospheric

sciences

，

1963

，

20

(

2

):

130-141.

[

13

]

Ruelle

D

，

Takens

F.

On

the

nature

of

turbulence

[

J

]

.

Les

rencontres

physiciens

-

mathematiciens

de

machine

:

solving

the

chaotic

three-body

problem

using

deep

neural

networks

[

J

]

.

arXiv

preprint

:

1910.

07291

，

2019.

Strasbourg-RCP25

，

1971

，

12

:

1-44.

[

14

]

Strogatz

S

H.

Nonlinear

dynamics

and

chaos

with

student

solutions

manual

:

With

applications

to

physics

，

Stars

cannot

be

surpassed

------

The

origin

of

three

-

body

problem

SHAO

Han-yong

(

Department

of

Physics

，

Beijing

Normal

University

，

Beijing

100875

，

China

)

Abstract

:

Since

Newton

in

the

Mathematical

Principles

of

Nature

Philosophy

for

the

first

time

to

define

three

body

problem

，

people

have

been

working

on

this

problem

for

300

years.

At

first

，

Euler

and

Lagrange

studied

some

relatively

simple

three

-

body

system

problems.

Poincare

first

considered

the

restricted

three

-

body

problem

，

but

found

that

its

trajectory

is

too

complex

to

be

described.

He

turned

to

prove

that

there

is

no

analytical

solution

for

the

three-body

problem

，

and

created

a

new

subject

，

nonlinear

dynamics.

Key

words

:

three-body

problem

；

Poincare

；

nonlinear

dynamics

(上接

50

页)

［

8

］

洪延姬，

窦志国

，

王明东

，

等

.

“

大学物理网络教学系

统

”

网络课程内容简介

［

〕］

•

大学物理

，

2001

，

20

(3

)

:

［

9

］

徐晓飞

，

李廉

，

战德臣

，

等

.

未来工程教育形态

：

面向可

持续竞争力的服务型教育

［

J

］

.

中国大学教学，

2019

,

39-42.

41

(

5

):

4-9.

Practice

and

introspection

on

the

online

teaching

of

college

physics

LI

Rui-shan

，

ZHANG

Pei-zeng

，

FENG

You-cai

(

1.

School

of

Science

，

Lanzhou

University

of

Technology

，

Lanzhou

，

Gansu

730050

，

China

；

2.

Department

of

Basic

Courses

，

Lanzhou

Institute

of

Technology

，

Lanzhou

，

Gansu

730050

，

China

)

Abstract

:

Recently

due

to

the

impact

of

the

pandemic

，

online

teaching

was

adopted

on

a

massive

scale

in

col

­

leges

and

universities

across

the

country.

In

this

paper

，

some

experiences

，

including

communication

forms

，

advan

­

tages

of

online

teaching

，

and

ideological

teaching

etc

，

are

concluded

from

the

teaching

practices

of

online

college

physics

，

and

difference

and

connection

between

online

teaching

and

classroom

teaching

are

explored.

Online

teach

­

ing

has

become

a

parallel

method

or

an

auxiliary

means

of

classroom

teaching.

In

the

long

term

，

the

in-depth

think

­

ing

of

the

online

teaching

practice

can

provide

valuable

ideas

for

further

improving

the

teaching

effect

of

college

physics.

Key

words

:

college

physics

；

online

teaching

；

classroom

teaching

；

ideological

teaching