第 4 章 数值实验台：把三体亲手跑起来
第 4 章 · 数值实验台：把三体亲手跑起来
章节摘要：本章跟着一次完整的数值实验走——从把方程翻译成积分器，到给积分器做误差对账，再到用庞加莱截面给混沌拍照，最后组装出第 3 章同款的周期解扫描流水线。四节连起来就是一间最小但五脏俱全的三体实验室，所有器材都只有百余行代码。
一条主线：一次数值实验的完整生命周期
设想你要回答一个具体问题："这个初值出发的三体轨道，是规则的还是混沌的？"第 2 章的方程给不了答案，第 3 章的猎物图谱也对不上号——你只能亲手把它跑出来。于是实验开始：先把牛顿方程装进一个积分器，选择格式与步长；跑之前先做误差对账，确认机器没有编造轨道；跑完的原始轨迹是一条时间序列，肉眼看不出规则与否，需要庞加莱截面把它压缩成一张可判读的照片；如果照片显示这附近有周期解的影子，再上扫描流水线精确定位。
这条主线串起的四件事，对应实验的四个环节：造仪器、验仪器、读数据、采样本。环节之间有严格的依赖——仪器没验过，照片不可信；照片没读懂，扫描就是大海捞针。本章的顺序就是实验的正确开工顺序。
图 4-1 数值实验的协作回路
沿途站点
第一站「最小可用的三体积分器」把第 2 章造好的受力零件装进四阶龙格-库塔框架，讨论步长的两难——大步长省时失真、小步长精确烧机时——并用二体解析解当标尺给整台仪器刻度。这一站结束时，你拥有一台能跑任何三体初值的通用机器。
第二站「辛积分器与误差对账」处理仪器最阴险的故障模式：RK4 跑长程时能量单向漂移，轨道被机器慢慢"加热"，最后呈现的可能是算法制造的伪混沌。蛙跳与维雷特这类辛格式用几何性质锁住能量误差的有界性，本站用长程对账实验把两种格式的差异量出来，并给出格式选择的工程准则。
第三站「庞加莱截面」解决"轨迹看不懂"的问题：把四维轨迹流切过一张二维剖面的点集留下来，规则轨道在照片上画成闭合曲线，混沌轨道铺成雾状散点。本站给出截面的构造代码与判读手册，它是第 5 章一切混沌讨论的视觉基础。
第四站「周期解搜索大实验」把三站器材组装成第 3 章猎人们的流水线：网格扫描负责广度、周期判据负责淘汰、牛顿精修负责精度。本站复盘搜索预算的分配逻辑——哪一步最烧算力、哪一步最容易造假阳性，让"上千族周期解"从新闻数字变成可复算的工程流程。
拐点与结论
本章的认知转折在第二站：数值误差不是噪音，而是一种会被几何放大的"伪物理"。非辛格式的能量漂移在混沌区被指数放大，机器的舍入误差顺着李雅普诺夫通道涌入轨道——不做对账的数值实验，产出的一半结论可能属于算法而非引力。三体数值实验的全部严谨性，都悬在这道对账工序上。
章末落点是一个比例感：一台百余行的积分器加一套两小时的流程，就足以复现 2013 年级别的大部分周期解发现。三体数值实验的门槛从来不在算力规模，而在流程纪律——仪器要验、照片要判读、判据要校准。
读完你应该
- 能从牛顿方程出发写出三体积分器的完整接口，说清状态向量、导数函数、步进循环各自的职责；
- 能用二体解析解给积分器做刻度，解释为什么收敛阶要用"步长折半、误差缩几倍"来量；
- 能陈述辛格式的几何含义，用能量误差曲线区分真混沌与伪混沌；
- 能构造限制性三体的庞加莱截面，按"闭合曲线与散点雾"判读规则与混沌；
- 能组装扫描加精修的周期解流水线，并为它分配合理的算力预算；
- 能列出三体数值实验的三道质检工序及各自的淘汰标准。
下一章的接力
仪器造好了、照片拍出来了，接下来是把照片背后的现象命名与定量化：什么算混沌、混沌有多快、结局有几种、规则与混沌的分界长什么样。下一章把本章产出的实验事实升级成概念框架——李雅普诺夫时间、结局分类学与分形边界。