第 10 章 · 逻辑斯蒂映射:通往混沌的餐桌
一维迭代 $x_{n+1} = rx_n(1-x_n)$ 是混沌科学的果蝇:参数从 0 调到 4,经历不动点、倍周期级联、周期窗口,最终进入混沌。所有混沌的通用语法(倍周期、普适常数、分形)都能在它身上看到。
倍周期级联的完整旅程
- $0 < r \le 1$:灭绝;$1 < r < 3$:稳定不动点(饱和种群);
- $3 < r < 3.449$:周期 2(两年 alternation);$3.449 < r < 3.544$:周期 4——每次分岔周期翻倍,参数间隔按常数比缩小:
费根鲍姆常数 $\delta$ 是普适的:任何"单峰映射"(一个极大值)都给出同一个 $\delta$——混沌不是算术巧合,是结构规律;
- $r_\infty \approx 3.5699$ 后进入混沌带,其间嵌着周期窗口(最显眼的是周期 3,$r \approx 3.83$)——"周期三蕴含混沌"(李-约克定理)。
混沌的操作性定义与李雅普诺夫指数
混沌 ≠ 复杂波形,而是对初值敏感:相邻轨道指数分离 $\delta x_n \sim e^{\lambda n}$。李雅普诺夫指数 $\lambda = \lim\frac{1}{n}\sum\ln|f'(x_i)|$;正 $\lambda$ 就是混沌的身份证(逻辑斯蒂在混沌区 $\lambda > 0$)。蝴蝶效应的定量版本:可预测时间 $\sim \frac{1}{\lambda}\ln\frac{L}{\delta_0}$——初始精度每提高十倍,预测期只延长固定一段。
不变分布与分支树
混沌轨道访问 $x$ 的频率收敛到一个不变密度 $\rho(x)$(幂律奇异,端点处 $1/\sqrt{x}$ 型);迭代图像沿对角线折叠产生分支树(cobweb)——几何法直接读出收敛、振荡、混沌。分支树 + 返回映射($x_{n+2}$ 对 $x_n$)是分析任意一维迭代的两件随身工具。
老师的提醒
- 手工实验最省的混沌课:拿计算器迭代 $rx(1-x)$,扫 $r$ 记录吸引子——一小时亲手看到倍周期与窗口,胜过读十页。
- 蝴蝶效应的正确含义:不是"蝴蝶引起风暴",而是有限精度下长期预测不可能——误差增长快于任何测量改进。
自检
- 倍周期级联的参数序列与两个普适常数?
- 李雅普诺夫指数的定义与混沌判据?
- 为什么"有限精度 ⟹ 可预测时间有限"?
更深一层:倍周期级联是"通向混沌的国道"
逻辑斯蒂映射的倍周期之路有一条"普适国道":不动点 → 周期 2 → 4 → 8 → … → 混沌,参数间隔比收敛到费根鲍姆常数 $\delta \approx 4.669$。普适的深刻在于:任何"单峰映射"(一个二次极大)——不论物理来源(流体对流、电子电路、-population-虫口)——都给出同一个 $\delta$ 与同一个 $\alpha \approx 2.503$(分支叉宽度的相空间标度)。重整化群式的解释:倍周期是"时间尺度自相似"的相变,标度常数是普适类常数——混沌门槛处的物理学与铁磁相变共享同一套重整化语法。
周期窗口的"画廊"也值得走一遍:周期 3 窗口($r \approx 3.828$,李-约克定理的明星:"周期三蕴含混沌"——一个 3 周期点强制存在一切周期的点与不可数混沌集)、周期 5、6 窗口……混沌带里嵌着无穷层窗口、窗口里又嵌混沌——分形的自相似在参数轴上展开。
推导细读:李雅普诺夫指数的"预测时钟"
李雅普诺夫指数 $\lambda = \lim_{n\to\infty}\frac1n\sum_i\ln|f'(x_i)|$ 的工程读法是预测时钟:初始误差 $\delta_0$ 增长到系统尺寸 $L$ 的时间 $T \approx \frac{1}{\lambda}\ln\frac{L}{\delta_0}$。对数依赖的残酷之处:初值精度提高 10 倍,预测期只延长 $\frac{\ln 10}{\lambda}\approx 0.87$ 个迭代步——蝴蝶效应的精确含义(有限精度下长期预测被指数增长率锁死)。逻辑斯蒂 $r = 4$ 的 $\lambda = \ln 2$:每步误差翻倍——这不是比喻而是可计算的数字。
在真实世界里:倍周期的"实验室工单"
倍周期级联在实验室的工单清单:滴水龙头(滴水间隔从恒定 → 二倍 → 四倍 → 随机)、非线性电路(Chua 电路与变容二极管电路是最方便的桌面混沌发生器)、声学(喇叭驱动的流体柱)、生理(血细胞生成的周期性粒细胞病模型)。工单的共同格式:扫一个驱动参数、记录分岔图、对费根鲍姆 $\delta$ 做标度检验——验证成本一天、收获是对"普适类"概念的终身肌肉记忆。
专题深入:倍周期级联是"通向混沌的国道"完整版
逻辑斯蒂映射的倍周期之路有一条"普适国道":不动点 → 周期 2 → 4 → 8 → … → 混沌,参数间隔比收敛到费根鲍姆常数 $\delta \approx 4.669$。普适的深刻在于:任何"单峰映射"(一个二次极大)——不论物理来源(流体对流、电子电路、虫口)——都给出同一个 $\delta$ 与同一个 $\alpha \approx 2.503$(分支叉宽度的相空间标度)。重整化群式的解释:倍周期是"时间尺度自相似"的相变,标度常数是普适类常数——混沌门槛处的物理学与铁磁相变共享同一套重整化语法。
周期窗口的"画廊"也值得走一遍:周期 3 窗口($r \approx 3.828$,李-约克定理的明星:"周期三蕴含混沌"——一个 3 周期点强制存在一切周期的点与不可数混沌集)、周期 5、6 窗口……混沌带里嵌着无穷层窗口、窗口里又嵌混沌——分形的自相似在参数轴上展开。
数学细节展开:李雅普诺夫指数的"预测时钟"完整版
李雅普诺夫指数 $\lambda = \lim_{n\to\infty}\frac1n\sum_i\ln|f'(x_i)|$ 的工程读法是预测时钟:初始误差 $\delta_0$ 增长到系统尺寸 $L$ 的时间 $T \approx \frac{1}{\lambda}\ln\frac{L}{\delta_0}$。对数依赖的残酷之处:初值精度提高 10 倍,预测期只延长 $\frac{\ln 10}{\lambda}\approx 0.87$ 个迭代步——蝴蝶效应的精确含义(有限精度下长期预测被指数增长率锁死)。逻辑斯蒂 $r = 4$ 的 $\lambda = \ln 2$:每步误差翻倍——这不是比喻而是可计算的数字。
工程案例细节:倍周期的"实验室工单"工程版
倍周期级联在实验室的工单清单:滴水龙头(滴水间隔从恒定 → 二倍 → 四倍 → 随机)、非线性电路(Chua 电路与变容二极管电路是最方便的桌面混沌发生器)、声学(喇叭驱动的流体柱)、生理(血细胞生成的周期性粒细胞病模型)。工单的共同格式:扫一个驱动参数、记录分岔图、对费根鲍姆 $\delta$ 做标度检验——验证成本一天、收获是对"普适类"概念的终身肌肉记忆。
例题
逻辑斯蒂映射 $x_{n+1} = rx_n(1-x_n)$,取 $r = 3.2$(周期 2 区)。从不动点方程出发说明不动点不稳定,并求周期 2 点。
解答
不动点:$x^* = 0$ 或 $1 - 1/r = 0.6875$。导数 $f'(x) = r(1 - 2x)$:$|f'(0.6875)| = |3.2\times(1-1.375)| = 1.2 > 1$ → 不稳定 ⟹ 轨道被踢向振荡。
周期 2:解 $f(f(x)) = x$ 的非不动点根。$f^2(x) - x = 0$ 因子分解后二次因子:
$r = 3.2$:$x_{1,2} = \frac{4.2\pm\sqrt{4.2\times0.2}}{6.4} = \frac{4.2\pm0.9165}{6.4} \Rightarrow x_1 \approx 0.513,\ x_2 \approx 0.799$
轨道在两点间轮流跳动(隔年 oscillation)——倍周期级联的第一级。