3.5 倍周期瀑布：逻辑斯蒂映射全解剖

3.5 倍周期瀑布：逻辑斯蒂映射全解剖

本节摘要：定义先行——倍周期瀑布是"参数缓慢滑过临界序列、稳定周期逐级翻倍、最终坠入混沌"的标准路径，其教科书载体是逻辑斯蒂映射 x_{n+1} = r·x_n·(1−x_n)。本节报出全部关键里程碑的精确数值（r=3、3.4495、3.5441、3.5644，累积于 3.5699），拆解费根鲍姆普适常数，并解释为什么混沌带里反而嵌着"周期三窗口"。这是 2.2 节"连续对离散"对照的总兑现。

1976 年，生态学家梅在《自然》杂志发表一篇后来被引用数万次的论文，标题直白得像求救信号——"简单的数学模型具有非常复杂的动力学"。他研究的对象简单到令人发笑：一条抛物线的反复迭代。本节沿他的路走一遍，你会发现通往混沌的路线图上每个里程碑都有精确到小数点后几位的坐标。

一条虫子引发的数学

逻辑斯蒂映射原本是虫口模型：x_n 是第 n 代种群占承载量的比例，r 是增长率。抛物线迭代 xn+1 = r·xn·(1−xn) 的一句话推导：世代不重叠，下一代数量与今代成正比、与剩余资源成正比。2.2 节已经证明它的连续版兄弟永守秩序，而离散版这里将经历完整的"升级之路"。分析从不动点开始：不动点 x* 满足 x* = r·x*(1−x*)，得 x* = 1 − 1/r（x=0 之外）。稳定性判据（2.2 节表格）是 |f'(x*)|<1，而 f'(x) = r(1−2x)，代入驻点得 |2−r|<1，即 1<r<3。**r=3 是第一道门**：越过它，不动点失稳，两点循环接棒——倍周期开始了。

瀑布全图与精确里程碑

图：逻辑斯蒂映射倍周期瀑布图

把里程碑钉成一张表，这些数字值得背：

δ 的来历值得多看一眼：各级瀑布的参数间隔 0.449、0.0946、0.0204……相邻相除收敛到 4.6692。惊人的是，把抛物线换成任何"单峰"迭代（正弦、三角），这个数一成不变——普适常数意味着"通往混沌的路线图"跨系统通用。1979 年前后 Libchaber 在液氦对流实验里测到了同样的收敛比，数值预言与实验在第四位小数上握手，混沌从数学玩具变成物理事实。

周期三窗口：混沌带里的秩序孤岛

分岔图在混沌带中间有一块清晰的"白带"：r 从 3.8284 起，系统突然只循环三个点（近似值 0.160、0.514、0.956），随后这个周期三自己也逐级倍周期、再次坠入混沌。秩序不是被混沌"挤走"后再也不回来，而是嵌在混沌带里反复重现。李-约克定理那句名言"周期三意味着混沌"讲的就是这个结构的另一面：一条轨道若周期三存在，则同系统内任意周期的轨道与不可数条混沌轨道共存——窗口与混沌带从来不是邻居，而是同一套迭代结构的两张面孔。秩序与混沌的交错，正是离散引擎最深的性格：你手里那个"增长率"旋钮每拧动一点，系统可能在秩序与混沌之间来回穿越好几次。

展开案例：用代码钉下每个里程碑

背景：分岔图上的里程碑应当亲手复现，而不是背表。

操作：写一个"周期探测器"——迭代去掉暂态、对末段相近值聚类，数值即周期：

import numpy as np

def period_of(r, drop=800, tail=400, tol=1e-3):
    x = 0.5
    for _ in range(drop):                    # 丢弃暂态
        x = r * x * (1 - x)
    vals = []
    for _ in range(tail):                    # 收集长期样本
        x = r * x * (1 - x)
        vals.append(x)
    uniq = []
    for v in sorted(vals):                   # 相近值归为同一点
        if not uniq or v - uniq[-1] > tol:
            uniq.append(v)
    return len(uniq)

for r in (2.9, 3.2, 3.5, 3.55, 3.84, 3.9):
    print(f"r={r}: 数出的周期 = {period_of(r)}")
​

结果（典型输出）：2.9→1、3.2→2、3.5→4、3.55→8、3.84→3、3.9→混沌（几十个互异值）。

解读：r=3.84 数出 3，正是周期三窗口的数值指纹；r=3.9 数出大量互异值且每次重跑散布不同，是混沌的直接体验。**变式**：把 tol 调小重跑 r=3.9 会数出更多"伪周期"——混沌轨道在有限精度下永远接近却从不等于旧值，这个挫败感本身就是测度意义上的遍历在敲门。

⚠️ 常见坑：倍周期是离散映射与二维以上连续系统的共同现象，但一维连续流没有它（它连振荡都没有，2.2 节的老门槛）。另外瀑布只在"单峰"迭代上普适——形状有两个以上峰的映射另有路线，别把普适性用过头。

💡 关键直觉：δ≈4.669 意味着"越靠近瀑布尽头，参数越不值钱"——想用微调 r 的办法逃离混沌，步长必须指数式缩小；工程上调参避开临界区，要预留的是指数级余量。

本节要点回顾

• 里程碑成串背：3 → 3.4495 → 3.5441 → 3.5644 → 3.5699；每级间隔按 δ≈4.669 收敛。
• 推导练手：不动点 x*=1−1/r，稳定性 |2−r|<1，第一道门 r=3 是两行代数的事。
• 周期三窗口（1+√8≈3.828）：秩序嵌在混沌带里，李-约克定理把它与混沌绑成一体。
• 普适性有实验背书：费根鲍姆常数跨映射一致、液氦对流实测命中，这是非线性动力学的高光时刻。

第三章图鉴到此收齐。下一章把镜头拉远：分形、噪声、耦合与时延——真实系统给吸引子图鉴补的四个注脚。