第 9 章选读 · Lorenz 系统与奇怪吸引子
1963 年 Lorenz 把大气对流的三维截断模型丢进计算机,看到了永不相交又永不停息的轨道——奇怪吸引子与现代混沌观自此诞生。这一节是三维混沌的入门地图。
Lorenz 方程与耗散流
参数 $\sigma = 10, b = 8/3, r = 28$ 时:三个不动点(原点 + 一对 $C^\pm$)先后失稳(原点在 $r = 1$ 跨临界、$C^\pm$ 在 $r \approx 24.7$ 亚临界 Hopf),系统落入一个有界但永不重复的吸引集:蝴蝶吸引子(两个"翅膀",轨道绕几圈后随机换边)。
奇怪吸引子的三个指纹
- 有界 + 耗散:体积收缩($\nabla\cdot\mathbf{F} = -(\sigma + 1 + b) < 0$),吸引子体积零;
- 拉伸折叠:局部指数分离(最大李雅普诺夫指数 $\lambda \approx 0.9 > 0$)+ 全局折叠回归——面团的揉制动作;
- 分形结构:维度非整数(关联维 $D \approx 2.06$:介于面与体之间)。混沌与分形是一体两面:时间上的不可预测 = 空间上的分形几何。
从 Lorenz 到混沌学的世界观
- Poincaré 截面:三维流降为二维映射,截面上的结构(折叠一次的曲线)直接显示拉伸折叠——分析高维混沌的标准手术;
- 分岔视角:$r$ 扫过时 Lorenz 经历"不动点 → 混沌 → 周期窗口 → 混沌"的非单调旅程——混沌不是终点而是领地;
- 可预测性的极限成为大气科学的方法论转折(Lorenz 原文标题就是"确定性非周期流"):集合预报、概率天气预报由此发端。
老师的提醒
- 判断吸引子"怪不怪"的三查:体积收缩吗(耗散)?李雅普诺夫指数正吗(混沌)?维度整数吗(分形)?三查齐备即奇怪吸引子。
- 二维自治系统不可能混沌(P-B 定理),至少三维流或一维半以上映射才有混沌——维度门槛是基本的。
自检
- Lorenz 方程中三项各自的物理来源?
- 奇怪吸引子的三个指纹与相互关系?
- Poincaré 截面如何把流降为映射?
更深一层:Lorenz 系统的"耗散三件套"
Lorenz 系统 $\dot x = \sigma(y - x)$、$\dot y = rx - y - xz$、$\dot z = xy - bz$ 的分析三件套:耗散($\nabla\cdot\mathbf{F} = -(\sigma + 1 + b) < 0$——体积以 $e^{-13.67t}$ 收缩,吸引子体积为零);对称($(x, y) \to (-x, -y)$ 不变——蝴蝶双翼的对称来源);有界(能量函数论证轨线最终进入一个球——永不逃逸)。三件套合体:轨线被压缩到一个零体积的集合上、仍指数分离($\lambda_{max} \approx 0.9$)——"压扁 + 拉伸折叠"的奇怪吸引子,分形维数 $\approx 2.06$(面与体之间)。
Poincaré 截面是分析三维混沌的标准手术:取 $z = r - 1$ 平面截轨线,三维流降为二维映射——截面上轨道点的"折叠带"结构直接展示拉伸折叠机制。截面法是高维混沌的显微镜。
推导细读:不动点失稳序列的"分岔编年史"
Lorenz 的参数编年史(固定 $\sigma = 10, b = 8/3$,扫 $r$):$r < 1$ 原点全球吸引(静止对流);$r = 1$ 原点跨临界失稳(对流开始);$1 < r < 24.06$ 一对 $C^\pm$ 稳定(定常对流);$r \approx 24.06$ 亚临界:$C^\pm$ 失稳前后混沌带已出现($r \approx 24.7$ 处 $C^\pm$ 的 Hopf 只决定"局部还回不回得去");经典混沌区 $r \approx 28$(Lorenz 1963 的原参数);更高处穿插周期窗口与倍周期——混沌不是"终点"而是参数轴上的"领地",与周期窗口交错如群岛。
在真实世界里:从 Lorenz 到"可预测性科学"
Lorenz 的遗产超出了吸引子图形:集合预报(既然初值不确定性指数放大,预报就从"一个未来"变成"未来的概率云"——现代天气预报每天跑几十个扰动成员)、混沌控制(OGY 方法:利用吸引子的无穷不稳定周期轨道,小扰动把系统"钉"在选中的周期轨道上——混沌不全是敌人,是"可编程的丰富性")、保密通信(混沌同步:两个耦合 Lorenz 系统可同步——混沌载波通信的历史提案)。Lorenz 1963 年那 10 页论文,是"确定性"与"不可预测"和解的出生证明。
专题深入:Lorenz 系统的"耗散三件套"完整版
Lorenz 系统 $\dot x = \sigma(y - x)$、$\dot y = rx - y - xz$、$\dot z = xy - bz$ 的分析三件套:耗散($\nabla\cdot\mathbf{F} = -(\sigma + 1 + b) < 0$——体积以 $e^{-13.67t}$ 收缩,吸引子体积为零);对称($(x, y) \to (-x, -y)$ 不变——蝴蝶双翼的对称来源);有界(能量函数论证轨线最终进入一个球——永不逃逸)。三件套合体:轨线被压缩到一个零体积的集合上、仍指数分离($\lambda_{max} \approx 0.9$)——"压扁 + 拉伸折叠"的奇怪吸引子,分形维数 $\approx 2.06$(面与体之间)。
Poincaré 截面是分析三维混沌的标准手术:取 $z = r - 1$ 平面截轨线,三维流降为二维映射——截面上轨道点的"折叠带"结构直接展示拉伸折叠机制。截面法是高维混沌的显微镜。
数学细节展开:不动点失稳序列的"分岔编年史"完整版
Lorenz 的参数编年史(固定 $\sigma = 10, b = 8/3$,扫 $r$):$r < 1$ 原点全球吸引(静止对流);$r = 1$ 原点跨临界失稳(对流开始);$1 < r < 24.06$ 一对 $C^\pm$ 稳定(定常对流);$r \approx 24.06$ 亚临界:$C^\pm$ 失稳前后混沌带已出现($r \approx 24.7$ 处 $C^\pm$ 的 Hopf 只决定"局部还回不回得去");经典混沌区 $r \approx 28$(Lorenz 1963 的原参数);更高处穿插周期窗口与倍周期——混沌不是"终点"而是参数轴上的"领地",与周期窗口交错如群岛。
工程案例细节:从 Lorenz 到"可预测性科学"(工程版)
Lorenz 的遗产超出了吸引子图形:集合预报(既然初值不确定性指数放大,预报就从"一个未来"变成"未来的概率云"——现代天气预报每天跑几十个扰动成员)、混沌控制(OGY 方法:利用吸引子的无穷不稳定周期轨道,小扰动把系统"钉"在选中的周期轨道上——混沌不全是敌人,是"可编程的丰富性")、保密通信(混沌同步:两个耦合 Lorenz 系统可同步——混沌载波通信的历史提案)。Lorenz 1963 年那 10 页论文,是"确定性"与"不可预测"和解的出生证明。
例题
Lorenz 系统:$\dot x = \sigma(y-x)$、$\dot y = rx - y - xz$、$\dot z = xy - bz$($\sigma=10, b=8/3$)。(1)计算相体积收缩率;(2)求 $r = 28$ 时三个不动点。
解答
(1) $\nabla\cdot\mathbf{F} = -\sigma - 1 - b = -(10 + 1 + 8/3) = -13.67 < 0$——任意初始体积元按 $e^{-13.67t}$ 收缩 → 耗散系统,吸引子体积为零(但维数非整数:$D\approx2.06$,蝴蝶吸引子)。
(2) 不动点:$(0,0,0)$ 与 $C^\pm = (\pm\sqrt{b(r-1)}, \pm\sqrt{b(r-1)}, r-1)$:
原点鞍、$C^\pm$ 在 $r = 28$ 已失稳(Hopf 于 $r\approx24.7$)——轨道无家可归,永不重复地绕两翼。
- 三个指纹自查:体积收缩(耗散)✓、$\lambda_{max} > 0$(混沌)✓、分形维(奇怪)✓。