3.2 分岔图实操:参数扫描与费根鲍姆常数测量


3.2 分岔图实操:参数扫描与费根鲍姆常数测量

本节摘要:把上一节的"逐个参数试"升级为"整条参数轴一次扫完",就得到混沌研究辨识度最高的图像——分岔图。本节给出完整的参数扫描实验流程,逐段判读图像上的每一种结构:倍周期级联、周期窗口、带合并;然后测量级联分裂点的收敛比率,重现费根鲍姆常数的发现过程,理解"普适性"为什么让这个数字超越了 logistic 映射本身。

上一节的复盘结束在一个悬念上:周期二的诞生在三,周期四在约三点四四九,周期八在约三点五四四……这一串分裂点越排越密,最后在约三点五六九九处彻底挤在一起。中间的排列规律是什么?混沌区内部是否还有秩序?这些问题都无法靠单点实验回答——需要一张全景图。本节就来做这张图,并在图上完成一个诺贝尔奖级别的业余发现的重演:一位物理学家用一台可编程计算器,从分裂间距里算出了一个属于所有平方映射的公共常数。

学习目标

阅读完本节,你应当能够:

  1. 独立完成分岔图的参数扫描实验,说出暂态抛弃量与采样策略为什么决定图的成色;
  2. 判读分岔图上的五类结构:主干、分岔、级联、窗口、带合并;
  3. 从分裂点序列测量费根鲍姆常数,说明其收敛方式与测量误差来源;
  4. 解释"普适性"的确切含义,指出它为什么对湍流等真实系统有启发价值。

一、实验设计:一张图怎么扫出来

分岔图的做法一句话就能说完:横轴放参数 r,纵轴放长期虫口值;对每个 r,迭代若干步把暂态扔掉,再记录一段轨迹的所有取值,作为纵坐标的散点竖着打上去。稳定区每个 r 只有一个纵值,图上是一根细线;周期二区每个 r 是两个交替值,图上分出两根线;混沌区每个 r 散成一团有结构的竖云。上一节的"去重计数"实验等于在竖着切片看这张图,本节则把几千个切片横着连起来。

实验设计里有两个决定成败的细节。其一,暂态抛弃要足够长。参数靠近分岔点时收敛极慢(临界慢化),抛弃步数不足会让"假的暂态值"混进散点,图像发虚;规范做法是抛弃数百步起步,靠近级联尖端还要加码。其二,参数步长与采样量要匹配目的。看结构,步长千分之一量级、每参数采样几百点足够;测分裂点(下一小节),则要在分裂点附近局部加密、加大抛弃步数、提高数值精度。工程上这叫"先粗扫定位、再细扫测量",和实验物理里先粗调再微调是一个道理。

# 分岔图扫描实验(伪代码) 对 r 从 2.5 到 4.0, 步长 0.001: x = 0.4 抛弃: 迭代 400 步 采样: 再迭代 400 步, 每步输出一行 (r, x) 绘图: 每行一个像素点, r 为横轴, x 为纵轴, 点重叠处加深 成色检查清单: - 周期二区应是两根干净的线 (抛弃不足会出现第三根模糊线) - 混沌区竖云应有清晰的疏密带 (采样过少会糊成均匀噪点) - 周期三窗口应在 r 约 3.83 处可见三根白线间的黑带

二、判读:图上的五类结构

粗扫得到的图,从左往右读过去有五类地标。主干:r 从二点五到三的单线,稳定不动点的领地。第一次分岔:r 等于三处单线劈成两根,周期二开场;此后三根分岔点逐级出现,两根变四根、四根变八根,这就是倍周期级联——在约三点四四九、约三点五四四、约三点五六四……一路排向尖端。级联尖端:约三点五六九九,分裂在有限参数范围内累积无穷次,线结构让位于带结构的混沌云周期窗口:混沌云中嵌着的"秩序岛",最显眼的是周期三窗口,其中还能看到窗口内部自己的小级联——分岔图有自相似的分形结构,放大任何窗口都像缩小整张图。带合并:混沌区的竖云并非均匀一团,而是分成若干条横带;参数推移时带数按二、四、八反向翻倍地合并,这是级联在混沌区的镜像续篇。

图 3-2 分岔图全貌示意:从稳定到混沌的完整路线图

图 3-2 分岔图全貌示意:从稳定到混沌的完整路线图

三、测量复盘:费根鲍姆常数的诞生

现在做本节的高潮实验。把级联里前几次分裂的参数位置尽量精确地定出来(局部细扫:在猜测位置附近以小步长扫描,观察去重计数何时从一个值变两个),记作 r1、r2、r3、r4……然后计算相邻间距之比:r2 减 r1 比 r3 减 r2、r3 减 r2 比 r4 减 r3……一组典型的细扫结果如下:

局部细扫记录(每点抛弃 5000 步,采样 2000 步,精度 1e-6): r1 = 3.000000 (周期二诞生) r2 = 3.449490 (周期四诞生) r3 = 3.544090 (周期八诞生) r4 = 3.564407 (周期十六诞生) r5 = 3.568759 (周期三十二诞生) 间距: d1 = 0.449490 d2 = 0.094600 d3 = 0.020317 d4 = 0.004352 比率: d1/d2 = 4.7514 d2/d3 = 4.6564 d3/d4 = 4.6682 -> 收敛于 4.6692...

比率没有停在某个"漂亮的"有理数上,而是收敛到一个约四点六六九的无理数——费根鲍姆常数 δ。发现它的过程值得一提:二十世纪七十年代中期,费根鲍姆用一台可编程台式计算器反复迭代平方型映射,注意到不同映射的分裂点间距之比都滑向同一个数。普适性的含义就在这里:这个常数不属于 logistic 映射,而属于所有"峰顶光滑的抛物线型"一维映射——把迭代函数换成正弦映射、换成二次项系数不同的其他映射,只要峰值处二阶导数非零,级联间距之比依然收敛到同一个数。这不是巧合,后来被重整化群方法从数学上论证:倍周期之路在尺度变换下存在普适的不动点结构。这也是混沌研究里罕见的"可以定量预言"的规律——第 6 章会提到,湍流 onset 实验与电路实验都曾按这个常数的预言对上过号。

测量复盘里还藏着一条方法论:用低精度计算器算出普适常数,靠的不是算力而是外推——级联越深,间距越小,单次测量的相对误差越大,但相邻比率却越来越稳。做数值实验时,"找一个随深度稳定的量"往往比"追求单点精度"更接近物理。第 8.1 节的工作流规范会把这条经验写进清单。

四、变式实验:窗口里还有一整棵小树

最后一个变式留给周期三窗口。把扫描区间缩到三点八二附近、步长压到万分之一量级重扫,你会看到窗口内部自成一体:三根线先稳定,随后其中一根再次分岔,长出一棵完整的倍周期小树,直到窗口中心也变成混沌,再以逆向级联退出窗口回到混沌带。此外,"周期三意味着混沌"这句定理级的话——李与约克在一九七五年证明:一维映射只要出现周期三轨道,就必然存在任意周期的轨道与不可数的混沌轨道集。你可以动手验证前半句:在窗口左侧边界细扫,恰好是从"最大两个值"直接跳到"三个交替值",没有任何中间周期——数学定理在分岔图上看得见。

把本节与上一节连起来看,logistic 的复盘就完整了:四种命运是定性地图,分岔图是全景路线,费根鲍姆常数是路线图上的公理级坐标。下一节换到连续时间的世界——那组三行方程将把同样的戏剧演成蝴蝶形状。

本节要点回顾

  • 扫描设计:横轴参数、纵轴长期值,暂态抛弃数百步起步;先粗扫看结构、再细扫做测量。
  • 五类地标:主干、分岔、级联、周期窗口、带合并;窗口内部藏有自相似的迷你分岔树。
  • 常数测量:分裂点间距之比收敛于约四点六六九;靠的是级联深处的稳定性而非单点精度。
  • 普适性:常数属于所有峰顶光滑的平方型映射,重整化群给其提供了理论根基,是混沌领域可定量预言的规律。
  • 周期三定理:出现周期三即蕴涵全部周期与混沌轨道集,李-约克定理让分岔图上的窗口有了定理身份。

下一节进入连续时间与三维相空间——洛伦兹系统将把拉伸折叠演成一场肉眼可见的双翼之舞。


作者与出处
原作者: 灏天文库
来源:灏天文库
整理: 灏天文库整理
由灏天文库平台收录,内容或由平台用户上传,仅供学习交流
发布者: 作者: 灏天文库 转发
评论区 (0)
U