本页目录

非线性 · 混沌与奇怪吸引子

对标:Strogatz《Nonlinear Dynamics and Chaos》/ Ott / Guckenheimer & Holmes | 前置:mech-03(相空间)、fl-04 1963 年 Lorenz 把对流方程截成三个模,想做天气预报,却发现了确定性系统中的不可预测。这一页讲清混沌的判据、通向它的普适道路,以及一个深刻的界限:决定论与可预测性是两回事

Lorenz 吸引子与 Logistic 映射的倍周期分岔。

图 fl-05.1左:Lorenz 奇怪吸引子(\(\sigma=10,\ \rho=28,\ \beta=8/3\)),轨道永不重复却被限制在有限区域。右:Logistic 映射的倍周期分岔,周期 \(1\to2\to4\to\cdots\) 在 \(r_\infty\approx3.5699\) 处进入混沌。

1. Lorenz 系统:从对流到混沌

对 Rayleigh–Bénard(fl-04)作三模 Galerkin 截断:

\[\dot x = \sigma(y-x),\qquad \dot y = x(\rho-z)-y,\qquad \dot z = xy-\beta z\]

\(x\) 为对流强度、\(y,z\) 为温度分布模式;\(\rho\propto\mathrm{Ra}\)

分岔序列\(\rho<1\) 静止态稳定 → \(\rho>1\) 出现两个对流不动点 → \(\rho\approx24.74\) 两者同时失稳 → 轨道无处可去,只能在两个不稳定焦点之间永久游荡——奇怪吸引子

三个定义性特征

注意 Lorenz 系统是耗散的(相空间体积以 \(\nabla\cdot\mathbf f = -(\sigma+1+\beta)<0\) 收缩),所以吸引子是零测度的分形集(Lorenz 吸引子的分形维数约 2.06)。"体积收缩到零,却又不能收缩到点或环"——这个张力正是奇怪吸引子存在的原因。

2. 判据:Lyapunov 指数

两条邻近轨道的间距

\[|\delta(t)| \approx |\delta_0|\,e^{\lambda t}\]

最大 Lyapunov 指数 \(\lambda>0\) 即为混沌的定量判据

可预测时间:设初始误差 \(\delta_0\)、可容忍误差 \(\Delta\)

\[T_{\text{pred}} \sim \frac{1}{\lambda}\ln\frac{\Delta}{\delta_0}\]

这个对数关系是残酷的把初始测量精度提高 \(10\) 倍,预报时间只延长 \(\ln10/\lambda\)——一个常数。要把天气预报从 5 天延到 10 天,需要初始数据精度提升约 \(e^{5\lambda}\) 倍。这不是技术问题,是动力学的性质。

决定论 ≠ 可预测:Lorenz 方程完全确定,但长期行为实际不可预测。这一点在哲学上的分量,怎么强调都不过分——拉普拉斯妖的失败不需要量子力学,经典混沌就够了。

3. 通向混沌的普适道路

倍周期分岔(Feigenbaum):以 Logistic 映射 \(x_{n+1}=rx_n(1-x_n)\) 为例,周期 \(1\to2\to4\to8\to\cdots\),分岔点间距比趋于普适常数

\[\delta = 4.669201\ldots,\qquad \alpha = 2.502907\ldots\]

普适性是关键\(\delta\) 与具体系统无关,只要映射有二次极值。流体实验(Taylor–Couette、对流)、电路、化学振荡中都测到了同一个 \(\delta\)

这与 asm-03 的重整化群是同一个思想——Feigenbaum 正是用 RG(对映射作重标度变换,寻找不动点)证明了普适性。临界现象的普适类概念,在这里以完全不同的面貌重现。

另外两条道路【引用】:准周期(Ruelle–Takens,环面破裂)与阵发(intermittency)。

4. 保守系统:KAM 定理

耗散系统有吸引子;哈密顿系统(mech-03)没有——Liouville 定理保证相空间体积守恒。

KAM 定理:可积系统受小扰动时,大多数不变环面存活(那些频率比"足够无理"的),只有共振环面破裂形成混沌层。

\[\text{小扰动} \Rightarrow \text{规则与混沌区域交织共存}\]

后果:太阳系不是简单的"稳定"或"混沌",而是混合相空间——大行星轨道在数亿年尺度上规则,小天体与共振区则混沌(柯克伍德空隙的成因)。长期数值积分显示太阳系的 Lyapunov 时间约为数百万到数千万年量级——我们无法精确预测十亿年后行星的位置

5. 练习与要点

例 1(蝴蝶效应的量化) 大气 \(\lambda^{-1}\sim1\)–2 天。初始误差 \(10^{-3}\)、容忍 \(1\)\(T\sim\lambda^{-1}\ln10^3\sim7\lambda^{-1}\sim10\) 天——与业务预报的实际有效期一致"两周"不是技术极限,是动力学极限。

例 2(为什么混沌不等于随机) 混沌轨道有确定的吸引子结构、可算的分形维数、可重构的相空间(Takens 嵌入定理)——从一维时间序列可以重建吸引子。这使"是混沌还是噪声"成为可判别的问题。

例 3(混沌控制) 奇怪吸引子中嵌有稠密的不稳定周期轨道;用极小的扰动可把系统稳定到其中一条(OGY 方法)——🔗 自动化站的反馈思想在此与非线性动力学相遇。\(\blacksquare\)


下一页:给流体加上电磁场。等离子体是宇宙中占比最大的物质态,它的行为由磁流体方程支配——而"磁力线被冻结在流体中"这一条,几乎解释了它的一切奇特之处。