力学 · 从零学起 第 5 / 6 节

第 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——每次分岔周期翻倍,参数间隔按常数比缩小:
$$\frac{r_n - r_{n-1}}{r_{n+1} - r_n} \to \delta \approx 4.669$$

费根鲍姆常数 $\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$)是分析任意一维迭代的两件随身工具。

原文要点:第 10 章强调三个"普适":费根鲍姆 $\delta$(参数标度)、$\alpha \approx 2.5029$(相空间标度:分支叉的宽度比)、以及倍周期路径本身的普适性(二次极大 → 同一路径);书中把李-约克定理"周期三 ⟹ 周期一切 + 不可数混沌集"作为混沌存在性的严格陈述。

老师的提醒

  • 手工实验最省的混沌课:拿计算器迭代 $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$ 因子分解后二次因子:

$$x^2 - \frac{r+1}{r}x + \frac{r+1}{r^2} = 0\ \Rightarrow\ x_{1,2} = \frac{r+1\pm\sqrt{(r+1)(r-3)}}{2r}$$

$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)——倍周期级联的第一级。

本页是《非线性动力学与混沌》第 5 节的学习笔记:为便于初学,讲解、例题与要点概括均为本站原创整理,公式与章节定位供对照原书使用。

📑 本书目录 经典研读