力学

倍周期通往混沌:Lyapunov 指数的量化

导读:Strogatz 用逻辑斯蒂映射 $x_{n+1} = rx_n(1-x_n)$ 展示通往混沌最著名的路径——倍周期级联,并用 Lyapunov 指数量化混沌强度。

倍周期级联

$r$ 从 2 增到 3:稳定不动点 → 2 周期 → 4 周期 → 8 周期……分岔参数间距按 Feigenbaum 常数 $\delta = 4.669$ 几何收敛——普适常数,与具体映射无关。$r_\infty \approx 3.57$ 进入混沌。

周期窗口:混沌区内嵌 3 周期窗口($r\approx 3.83$)——按 Sharkovskii 序,「3 周期蕴含一切周期」。分岔图是必做的数值作业(本站 Duffing 文章有完整数值流程)。

Lyapunov 指数

$\lambda = \lim \frac{1}{n}\sum \ln|f'(x_i)|$:相邻轨道的指数分离率。$\lambda<0$ 收敛、$=0$ 周期分岔点、$>0$ 混沌。逻辑斯蒂映射的 $\lambda(r)$ 图像在分岔点过零、混沌区正负交替——混沌的定量身份证。

工程含义:正 Lyapunov 指数 = 初值误差指数放大 = 长期预测不可能(蝴蝶效应的严格表述)。天气预报的预测时限本质上由大气的大 Lyapunov 指数决定。

通往混沌的其他路径

准周期路径(Ruelle-Takens,三次霍普夫)与间歇性道路(Pomeau-Manneville):倍周期不是唯一入口。Duffing 振子在周期激励下即可展现多条路径——数值实验的最佳对象。

常见误区

  • 把数值伪振荡当混沌(须验证正 Lyapunov 指数 + 吸引子有界)。
  • 用固定积分步长的粗欧拉法研究混沌,把数值误差当动力性态。

吃透自检

  • 能编程绘制逻辑斯蒂映射分岔图与 Lyapunov 图
  • 能解释 Feigenbaum 常数的普适性含义
  • 能区分三条通往混沌的路径

← 力学 书目精读列表 经典研读