非线性 · 混沌与奇怪吸引子
对标:Strogatz《Nonlinear Dynamics and Chaos》/ Ott / Guckenheimer & Holmes | 前置:mech-03(相空间)、fl-04 1963 年 Lorenz 把对流方程截成三个模,想做天气预报,却发现了确定性系统中的不可预测。这一页讲清混沌的判据、通向它的普适道路,以及一个深刻的界限:决定论与可预测性是两回事。
1. Lorenz 系统:从对流到混沌
对 Rayleigh–Bénard(fl-04)作三模 Galerkin 截断:
\(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 指数
两条邻近轨道的间距
最大 Lyapunov 指数 \(\lambda>0\) 即为混沌的定量判据。
可预测时间:设初始误差 \(\delta_0\)、可容忍误差 \(\Delta\),
这个对数关系是残酷的:把初始测量精度提高 \(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\) 与具体系统无关,只要映射有二次极值。流体实验(Taylor–Couette、对流)、电路、化学振荡中都测到了同一个 \(\delta\)。
这与 asm-03 的重整化群是同一个思想——Feigenbaum 正是用 RG(对映射作重标度变换,寻找不动点)证明了普适性。临界现象的普适类概念,在这里以完全不同的面貌重现。
另外两条道路【引用】:准周期(Ruelle–Takens,环面破裂)与阵发(intermittency)。
4. 保守系统:KAM 定理
耗散系统有吸引子;哈密顿系统(mech-03)没有——Liouville 定理保证相空间体积守恒。
KAM 定理:可积系统受小扰动时,大多数不变环面存活(那些频率比"足够无理"的),只有共振环面破裂形成混沌层。
后果:太阳系不是简单的"稳定"或"混沌",而是混合相空间——大行星轨道在数亿年尺度上规则,小天体与共振区则混沌(柯克伍德空隙的成因)。长期数值积分显示太阳系的 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\)
下一页:给流体加上电磁场。等离子体是宇宙中占比最大的物质态,它的行为由磁流体方程支配——而"磁力线被冻结在流体中"这一条,几乎解释了它的一切奇特之处。