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

第 9 章选读 · Lorenz 系统与奇怪吸引子

1963 年 Lorenz 把大气对流的三维截断模型丢进计算机,看到了永不相交又永不停息的轨道——奇怪吸引子与现代混沌观自此诞生。这一节是三维混沌的入门地图。

Lorenz 方程与耗散流

$$\dot x = \sigma(y - x), \quad \dot y = rx - y - xz, \quad \dot z = xy - bz$$

参数 $\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 原文标题就是"确定性非周期流"):集合预报、概率天气预报由此发端。
原文要点:第 9 章强调 Lorenz 系统的"确定性"与"非周期"并存:方程完全确定、初值完全已知,轨道仍不可长期预测——敏感性取代了外来的随机性;书中以原点与 $C^\pm$ 的失稳序列、以及"轨道绕翅数的不规则性"作为混沌不可预测性的最小演示。

老师的提醒

  • 判断吸引子"怪不怪"的三查:体积收缩吗(耗散)?李雅普诺夫指数正吗(混沌)?维度整数吗(分形)?三查齐备即奇怪吸引子。
  • 二维自治系统不可能混沌(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)$:

$$x_{\pm} = \pm\sqrt{\frac83\times27} = \pm\sqrt{72} \approx \pm8.49, \qquad z = 27$$

原点鞍、$C^\pm$ 在 $r = 28$ 已失稳(Hopf 于 $r\approx24.7$)——轨道无家可归,永不重复地绕两翼。

  • 三个指纹自查:体积收缩(耗散)✓、$\lambda_{max} > 0$(混沌)✓、分形维(奇怪)✓。

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

📑 本书目录 经典研读