本页目录

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

对标：Strogatz《Nonlinear Dynamics and Chaos》/ Ott / Guckenheimer & Holmes ｜ 前置：相空间与 Hamilton 方程、流动失稳与振幅方程 1963 年 Lorenz 把对流方程截成三个模，发现了确定性系统中的实际不可预测性。本页先用 Logistic 映射做一项 先预测再观察的实验，再把有限时间诊断放回 Lorenz、Feigenbaum 普适性和 KAM 的适用语境中：决定论、混沌和不可预测不是同一个词。

Lorenz投影与Logistic数值扫描及解析周期2分支
图 fl-05.1左：DOP853 积分（相对容差 2×10⁻¹²、绝对容差 2×10⁻¹⁴），初值 (1,1,20)，舍去 t<10 后绘制到 t=40；Lorenz 三模在经典参数 \((\sigma,\rho,\beta)=(10,28,8/3)\) 下的典型奇怪吸引子；右：601 个参数逐一从 x₀=0.2 烧入 600 步，再绘出 64 步；红线叠加解析周期 2 分支。双精度扫描可能遗漏窄窗口或有长暂态。两幅图是模型语境的入口，不是由一张有限采样图推出“永不重复”或“必然混沌”的证明。

学习层：先预测，再打开 Logistic 映射

1. 固定模型与三问预测门

实验只使用一个确定的离散映射：

\[ x_{n+1}=f_r(x_n)=r x_n(1-x_n), \qquad 0\le r\le4,\quad 0\le x_0\le1. \]

这是精确递推规则，没有随机噪声。因为 \(f_r([0,1])\subseteq[0,1]\)，两条轨道的距离

\[ D_n=|x_n-y_n|\le1, \qquad y_0=x_0\pm\delta_0\in[0,1] \]

有一个硬上限。实验保留你设定的 \(x_0\)：优先向右取邻点，越过 1 时改为向左，实际方向会显示。\(\delta_0\) 是原始两初值的目标间距；共同烧入后，图中重新编号的 \(D_0\) 通常已不等于它。局部线性化只在 \(D_n\) 仍然很小时可信：

\[ D_N\approx D_0\exp\left(\sum_{n=0}^{N-1}\log|r(1-2x_n)|\right). \]

打开实验前，请先回答三件事：

  1. 规则完全确定，是否因此保证长期可预测？
  2. 邻轨的指数分离能否无限持续？
  3. 混沌参数区间是否没有周期窗口？

交互版会在提交三项预测后才显示预设、参数、时间序列、邻轨分离和

\[ \lambda_N=\frac1N\sum_{n=0}^{N-1}\log|r(1-2x_n)|. \]

揭示后可以连续调 \(r,x_0,\delta_0\)、烧入步数和有限窗 \(N\)；切换预设会保留揭示状态，重置才会重新关门。

无 JavaScript 时的静态读法：固定 \(\delta_0=10^{-8}\)。表中的“有限窗诊断”只描述给定初值、烧入和采样长度的计算，不能把一张非周期图当作定理。

预设 \(r\) \(x_0\) 烧入后读法 \(\lambda_N\) 应如何解释
稳定不动点 2.8 0.2 \(x_n\to x_*=(r-1)/r=5/7\) 长窗趋向 \(\log\lvert2-r\rvert=\log0.8<0\)
稳定周期 2 3.2 0.2 两个值交替；固定点已失稳 周期轨道的长窗平均为负
稳定周期 4 3.5 0.2 四个值轮流出现 负值只说明这条周期轨道附近的收缩
不稳定不动点 4 0.75 基准轨道恒等于 0.75，邻轨可以离开 \(\lambda_N=\log2>0\)，但这条基准轨道是周期 1
周期 3 窗口附近 3.83 0.2 周期 3 的稳定窗口附近 负值是该窗口的局部结论，不是整个混沌区间的结论
典型混沌诊断 3.9 0.2 有限样本通常不显示短周期 常见有限窗给出正值；这是轨道/测度/时间尺度相关的诊断
\(r=4\) 的异常初值 4 0.5 \(0.5\to1\to0\to0\to\cdots\) 若求和窗口含 \(x_0=0.5\)，有 \(\log0=-\infty\)

当 \(r>0\) 时，固定点方程的候选解是 \(x=0\) 和 \(x_*=(r-1)/r\)（\(r=1\) 时重合；\(r=0\) 只有 0，\(0<r<1\) 的第二个候选在区间外）；非零固定点的导数为 \(2-r\)，所以它在 \(1<r<3\) 稳定。\(r=3.2\) 的周期 2 与 \(r\approx3.83\) 的周期 3 窗口，正是提醒我们：所谓“混沌区间”不是一整块没有周期结构的平面。

2. 如何读有限时间账本

实验中的 \(\lambda_N\) 会明确处理三个容易被图形遮住的事件：

  • 若 \(x_n=1/2\)，则 \(f'_r(x_n)=0\)，对应项是 \(\log0=-\infty\)，不是 JavaScript 的“坏数值”；
  • 若两条 IEEE 双精度轨道变成同一个浮点数，后续 \(D_n=0\)。这既可能涉及舍入，也可能反映精确映射的多对一结构：\(f_r(x)=f_r(1-x)\)。仅凭数值相等不能认定原因；图上把零距离放在独立红点行，不把 \(-\infty\) 伪装成有限刻度；
  • \(D_n\le1\) 限制了误差大小，却不要求距离单调增长、停在 1 或保持某个平台。实验的 0.95 只是接近上界的显示阈值；有限距离可能反复增长、缩小或合流。

因此，\(\lambda_N>0\) 是常用的有限时间混沌诊断，却依赖轨道、初始分布或不变测度以及 \(N\)。\(\lambda_N<0\) 也可能只是落在稳定周期窗口；一次正值或一张“看起来不重复”的图都不单独构成定理证明。

3. 从一步差分到指数：被忽略的二次项

令有符号误差 \(e_n=y_n-x_n\)。直接展开同一映射，而不是先假定指数律：

\[ e_{n+1}=r(1-2x_n)e_n-r e_n^2. \]

切向量 \(v_n\) 的定义是把初始扰动先趋于零，因此满足线性变分方程

\[ v_{n+1}=r(1-2x_n)v_n,\qquad \frac1N\log\left|\frac{v_N}{v_0}\right|=\lambda_N. \]

这条乘积恒等式对切向量成立；代入有限距离时，必须控制二次项。尤其在 \(x_n=1/2\) 处，一阶项为零，但有限误差下一步是 \(-r e_n^2\)，并不一般为零。因此 \(\lambda_N=-\infty\) 不是“任意有限扰动都被完全消掉”。

烧入也改变问题：取 \(r=4,x_0=1/2\)，包含初始临界点的切向乘积为零；烧入两步后从 \(x=0\) 重新施加切向扰动，得到 \(\lambda_N=\log4\)。两者对应不同起点的导数，不矛盾。实验会分别标记烧入和采样窗中的零导数。

4. 手算一个周期，再构造一个反例

从 \(f_r(f_r(x))=x\) 除去不动点因子，得到周期 2 候选

\[ x_\pm=\frac{r+1\pm\sqrt{(r-3)(r+1)}}{2r},\qquad m_2=f'_r(x_+)f'_r(x_-)=-r^2+2r+4. \]

在 \(r>3\) 两点才分开；由 \(|m_2|<1\) 得到稳定范围 \(3<r<1+\sqrt6\)。在 \(r=3.2\)，两点约为 \(0.513045,0.799455\)，\(m_2=0.16\)，每步平均指数为 \(\tfrac12\log0.16=\log0.4<0\)。零乘子处指数取 \(-\infty\)。实验的周期检测检查末至多 80 点、至少四个完整周期块、绝对容差 \(10^{-8}\)，只是“接近周期”的有限检验。

再设 \(r=4,x_0=3/4\)。此时轨道始终是 \(3/4\)，但 \(f'_4(3/4)=-2\)，所以每个 \(N\) 都有 \(\lambda_N=\log2>0\)。这是不稳定不动点，足以否定“正指数就证明当前轨道非周期”。可以把实验的 \(r\) 调为 4、\(x_0\) 调为 0.75 亲自复算。

为何同一个 \(\log2\) 又经常代表 \(r=4\) 的典型混沌？令 \(x=\sin^2(\pi\theta)\)，映射成为 \(\theta\mapsto2\theta\pmod1\) 的投影。均匀角度推到 \(x\) 上给出概率密度

\[ p(x)=\frac1{\pi\sqrt{x(1-x)}},\qquad \int_0^1\log|4(1-2x)|p(x)\,dx =\log4+\int_0^1\log|\cos(2\pi\theta)|\,d\theta=\log2. \]

倍增映射的遍历性使这个平均描述相应测度下几乎处处的实数初值，而不是所有初值；上述不动点、临界点及其前像仍须单独处理。双精度有限状态实验也不能替代这个几乎处处结论。Fremlin 的遍历论讲义 给出测度保持与混合的严格背景。

1. Lorenz 系统：从对流到典型奇怪吸引子

对 Rayleigh–Bénard 对流作三模 Galerkin 截断，得到

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

\(x\) 可理解为对流强度，\(y,z\) 是温度分布的模式，\(\rho\) 与 Rayleigh 数成正比。平衡点为

\[ C_0=(0,0,0),\qquad C_\pm=\left(\pm\sqrt{\beta(\rho-1)},\ \pm\sqrt{\beta(\rho-1)},\ \rho-1\right)\quad(\rho>1). \]

以下取 \(\sigma,\beta,\rho>0\)。原点的一个本征值为 \(-\beta\)，另两个满足

\[ \lambda^2+(\sigma+1)\lambda+\sigma(1-\rho)=0. \]

故 \(0<\rho<1\) 时原点渐近稳定；\(\rho>1\) 时失稳并出现 \(C_\pm\)。对非零平衡的 Jacobian 取行列式，得到

\[ \lambda^3+(\sigma+\beta+1)\lambda^2 +\beta(\sigma+\rho)\lambda+2\sigma\beta(\rho-1)=0. \]

三次 Routh–Hurwitz 条件除系数为正外，还要求 \((\sigma+\beta+1)\beta(\sigma+\rho)>2\sigma\beta(\rho-1)\)。若 \(\sigma>\beta+1\)，边界为

\[ \rho_H=\frac{\sigma(\sigma+\beta+3)}{\sigma-\beta-1}. \]

经典 \((\sigma,\beta)=(10,8/3)\) 给出 \(\rho_H=470/19\approx24.7368\)。这是共轭本征值穿越虚轴的局部信息；判断 Hopf 分岔的非退化性、分支稳定性及全局吸引子还需非线性分析，不能仅凭平衡失稳推出“只能混沌”。若 \(\sigma\le\beta+1\)，这条正参数 Hopf 阈值公式不能当作一个可跨越的正阈值使用。

在 \((10,28,8/3)\) 的经典模型中，奇怪吸引子的存在已有 Tucker 的计算机辅助证明，结合严格区间计算与解析论证。普通浮点蝴蝶图只是数值示例，不是那份证明，也不是对真实大气的证明。

常见教材会把下面三项放在一起：有界、非周期、对初值敏感。它们是很有用的操作性特征组，却不是对所有“混沌”定义都普适的充分必要三条件：有限观测中的非周期可能只是周期很长，敏感性也可能是暂态；不同的拓扑、测度论或符号动力学定义会强调不同结构。

Lorenz 流的散度为

\[ \nabla\cdot\mathbf f=-(\sigma+1+\beta)<0, \]

由 Liouville 公式，沿流搬运的小体积乘以 \(e^{-(\sigma+1+\beta)t}\)。体积收缩本身不保证轨道有界：二维系统 \(\dot u=u,\dot v=-2v\) 的散度也为负，但 \(u\) 可以逃向无穷。

Lorenz 的有界性可以另外算出来。取

\[ V=\rho x^2+\sigma y^2+\sigma(z-2\rho)^2. \]

代入方程后，\(xyz\) 与 \(xy\) 交叉项抵消：

\[ \dot V=-2\rho\sigma x^2-2\sigma y^2 -2\sigma\beta(z-\rho)^2+2\sigma\beta\rho^2 \le-cV+4\sigma\beta\rho^2,\quad c=\min(2\sigma,2,\beta)>0. \]

最后一步用了 \((z-\rho)^2\ge\tfrac12(z-2\rho)^2-\rho^2\)。积分得到

\[ V(t)\le e^{-ct}V(0)+\frac{4\sigma\beta\rho^2}{c}(1-e^{-ct}). \]

因此足够大的椭球最终容纳全部轨道；这是一个宽松的吸收界，并不确定吸引子的精细形状。若要估计维数，必须说明是盒维数、相关维数还是 Lyapunov/Kaplan–Yorke 维数；本页的轨道图不输出维数，不能从散度或像素密度读出一个精确数字。

2. Lyapunov 指数与预测时间

对流的最大 Lyapunov 指数可写成切向量增长率的极限形式；对合适的光滑系统、不变遍历测度和典型轨道，最大指数为正是标准的混沌证据或判据语境。它不是脱离语境的万能测试：有限时间正值可能来自暂态或不稳定周期轨道，非遍历轨道也不能替代测度平均；Logistic 实验中的 \(\lambda_N\) 也只针对当前轨道和当前时间窗。

在局部指数增长阶段，若初始误差为 \(\delta_0\)，可容忍误差为 \(\Delta\)，则

\[ T_{\mathrm{pred}}\sim\frac1\lambda\log\frac{\Delta}{\delta_0}. \]

这个对数关系表达的是误差预算，不是无条件的长期公式：进入非线性区、碰到边界或观测误差结构改变后，单一 \(\lambda\) 不再足够。做一个纯粹的误差预算例子：假设某模型的 \(\lambda=0.5\,\mathrm{s^{-1}}\)，初始误差 \(10^{-6}\)、容限 \(10^{-2}\)（同一单位），则 \(T\approx2\log10^4=18.42\,\mathrm{s}\)。把初始误差缩小 100 倍，只多争取 \(2\log100=9.21\,\mathrm{s}\)，而不是把预测时间乘以 100。这里的秒和指数都是假设值，不能拿来标定天气或太阳系的实际预测期限。

原始 Lorenz 论文 将低维对流模型与长期预报问题联系起来。用今天的模型预测具体变量，还必须另行评估观测、误差增长阶段、模型误差和所选容限；本页不提供跨系统的统一预报年限。

3. Feigenbaum：倍周期路线的普适性边界

Logistic 映射展示了典型的倍周期序列：周期 \(1\to2\to4\to8\to\cdots\)，分岔点 \(r_k\) 的间距比趋向

\[ \delta=\lim_{k\to\infty}\frac{r_{k-1}-r_{k-2}}{r_k-r_{k-1}} \approx4.669201609\ldots, \]

而临界点附近的空间重标度给出尺度比的绝对值 \(|\alpha|\approx2.502907875\ldots\)（带方向的重标度约定可取负号）。Feigenbaum 原始论文 的重整化群把“换尺度后仍像自己”的不动点变成了普适性来源；它与 重整化群 中的换尺度思路相连。

这里的“普适”必须加限定：结论属于相应的一维、足够光滑的单峰映射类，临界点为二次极值，并满足相应的横截性、缩放和动力学正则条件。它不是“所有非线性系统都有同一个 \(\delta\)”；流体、电路、化学振荡等系统若落入相应的普适类，才会在合适观测量上显示同样的标度。

倍周期不是唯一道路。准周期环面破裂（Ruelle–Takens）和阵发（intermittency）提供了不同路线；即便在 Logistic 映射的混沌区间，也会嵌入稳定周期窗口。\(r\approx3.83\) 的周期 3 窗口和 \(r=4,x_0=1/2\) 的异常轨道，分别提醒我们区分参数结构与初值例外。

4. 保守系统：KAM 定理的测度语境

耗散 Lorenz 系统收缩体积；完整 Hamilton 系统则由 Liouville 定理保持相空间体积，因而不能把一片正体积区域压成通常意义的吸引子。KAM 讨论的是更细的结构：从可积 Hamiltonian

\[ H(I,\theta)=H_0(I)+\varepsilon H_1(I,\theta) \]

出发，在足够小扰动、足够光滑（经典版本常取解析）且非退化条件 \(\det\partial^2H_0/\partial I^2\ne0\) 下，满足 Diophantine 无理频率条件的许多不变环面继续存在。所谓“大多数”是相对于作用量区域的测度陈述：被破坏环面的测度在 \(\varepsilon\to0\) 时趋于零，而不是说每一个环面都存活。共振环面、低正则性、非小扰动和退化情形都需要另行分析。

把“足够无理”写成可以核对的条件：对 \(n\) 个自由度、\(\omega(I)=\partial_IH_0\)，要求存在 \(\gamma>0,\tau>n-1\)，使

\[ |k\cdot\omega|\ge\gamma |k|^{-\tau},\qquad 0\ne k\in\mathbb Z^n. \]

这个下界约束解同调方程时出现的小分母。在固定紧致作用量区域的内部、统一非退化与解析范数控制下，扰动上限依赖 \(\gamma\) 及这些常数。不能先任取 \(\gamma\to0\) 再沿用同一个扰动阈值。标准定量版本可把未被保留的主要 KAM 环面覆盖的相空间体积估成 \(O(\sqrt\varepsilon)\)，常数依赖区域和原系统；这仍不判定剩余集合全是混沌。Biasco–Chierchia 的定量估计 专门澄清这种依赖。

因此小扰动下常见的是规则环面、共振岛和混沌层交织的混合相空间。太阳系的长期动力学可用这个框架获得直觉，但 KAM 假设不能直接替真实多体模型盖章；“大行星在许多时间尺度上规则”和“小天体在共振区更敏感”可以并存，不能简写成无条件的“太阳系稳定”或“太阳系混沌”。

5. 具象例子：同一个“混沌”词的不同含义

6. 误区表：把诊断还原成条件句

过强说法 更准确的读法
Lorenz 非零平衡失稳后，轨道无处可去只能吸引子 失稳是局部信息；经典参数和相应吸引域中可出现奇怪吸引子，但全局结论需分析不变集与吸引域
有界、非周期、敏感依赖初值是混沌的普适充分定义 这是常用教材特征组；有限样本和暂态会误导，定义依语境而异
从蝴蝶图能直接读出精确维数 需要先选维数定义，再给数值算法、收敛和误差证据；本图不测量维数
最大 Lyapunov 指数 \(>0\) 对任何数据都证明混沌 对合适系统和不变测度是重要诊断；有限时间正值不单独等于定理证明
天气两周后必然不可预测 “两周”是依变量、形势、观测和模型而变的经验量级，不是无条件定律
太阳系几百万年后必然完全失去可预测性 Lyapunov 时间是模型与口径相关的局部尺度；不同轨道族可同时规则与混沌
KAM 说大多数环面都永远保留 “大多数”是小扰动、非退化、Diophantine 语境下的测度陈述，共振环面可能破裂
Feigenbaum 常数适用于所有非线性系统 它限于相应的一维单峰、二次临界点等普适类
Logistic 图就是 Lorenz 或真实天气的缩小版 它是无空间、无噪声的一维离散教学模型；能教局部膨胀与边界，不替代流体或天气模型

7. 模型边界与收束

Logistic 映射把高维连续流压成一个标量和一个离散步长：它没有空间输运、守恒律、观测同化、外部强迫和测量噪声，也不包含 Lorenz 的三维几何。JavaScript 版还使用 IEEE 双精度，因此轨道舍入与有限字长重复需要单独记账；而 \(-\infty\) 还可能来自精确导数为零；数学上的精确轨道与有限字长轨道必须分开说。

它仍然是很好的第一台实验机：固定规则可以产生复杂时间序列，邻轨分离有指数阶段却最终受 \(D_n\le1\) 限制，周期窗会嵌在典型混沌参数附近，有限时间 Lyapunov 读数可以帮助提出而不是替代动力学论证。用这个顺序读 Lorenz、Feigenbaum 和 KAM，才不会把一张图、一项数值或一个漂亮常数误认成跨模型的无条件定理。\(\blacksquare\)

8. 迁移练习：先写理由，再展开答案

  1. 取 \(r=4,x_0=3/4\) 和 \(r=4,x_0=1/2\)，分别计算无烧入的 \(\lambda_N\)。后一条烧入两步再计算会怎样？哪个读数可以单独证明非周期？
  2. 若 \(r=3.2\)，周期 2 的误差每两步约乘几？如果两个初值是 \(1/4,3/4\)，第一步是否必然不同？解释这与有限距离和切向量的区别。
  3. 在 Lorenz 经典参数下，算出上述吸收界的 \(c\) 与 \(V\) 上限，并说明为什么它不能证明混沌。把初值精度提高 100 倍又为何不能普遍让预测时间增加 100 倍？
展开三题推导与检查点

1. \(x=3/4\) 恒定而导数模为 2，所以 \(\lambda_N=\log2\)。\(1/2\to1\to0\) 的首项为 \(\log0\)，无烧入时 \(\lambda_N=-\infty\)；烧入两步从 0 重新定义切向扰动，\(\lambda_N=\log4\)。都不单独证明当前轨道非周期；两个例子最终都是不动点轨道。需要分清切向扰动重新施加的位置。

2. 两步乘子为 \(0.16\)，每步指数 \(\log0.4\approx-0.916291\)，这是周期轨道附近小扰动的渐近线性读数。\(f_r(1/4)=f_r(3/4)=3r/16\)，这两条精确实轨道一步合流，不需要浮点舍入。对称点相距 1/2，不能只保留差分式的一阶项。

3. \(c=\min(20,2,8/3)=2\)，\(4\sigma\beta\rho^2/c=125440/3\approx41813.33\)。对任意更大的椭球 \(V\le R\)，所有轨道最终进入；它只给有界性。原点这个不动解也满足界，故有界不能等同于混沌。在固定正 \(\lambda\) 且误差仍处于指数阶段时，精度改善 100 倍只增加 \(\log100/\lambda\)；模型误差和非线性饱和还会使这个估计失效。

下一页：等离子体与磁流体。从确定性局部模型回到带空间输运的连续介质，继续检查“磁力线冻结”需要哪些理想化条件。