洛伦兹方程

洛伦兹方程(Lorenz equation)描述空气流体运动的一个简化微分方程组.1963年，美国气象学家洛伦兹(Lorenz,E. N.)将描述大气热对流的非线性偏微分方程组通过傅里叶展开，大胆地截断而导出描述垂直速度、上下温差的展开系数x(t),y(t),z(t)的三维自治动力系统 ^[1] :

洛伦兹吸引子当 ρ = 28, σ = 10, and β = 8/3 ^[2]

[math]\displaystyle{ \begin{cases} \dot{x}=P(y-x) \\ \dot{y}=Ra x-y-xz \\ \dot{z}=xy-bz \end{cases} }[/math]

其中：

[math]\displaystyle{ P }[/math]:普朗特数

[math]\displaystyle{ Ra }[/math]:为瑞利数

他发现当Ra 不断增加时，系统就由定常态(表示空气静止)分岔出周期态(表示对流状态)，最后，当Ra>24. 74时，又分岔出非周期的混沌态(表示湍流). 如图是三维相空间的混沌态在二维平面上的投影轨线.从图上看出，轨线起初在右边从外向内绕圈子，后来随机地跳到左边从外向内绕圈子，后又再次随机地跳回右边绕圈子，……如此左右跳来跳去，每次绕的圈数，何时发生跳跃都是随机的、无规则的. 由于洛伦兹是世界上第一个从确定的方程中发现了非周期的混沌现象，所以将上述方程一般称为洛伦兹方程. ^[3]

背景及历史

基本概念

1963年，美国气象学家爱德华·洛伦兹提出的一组简化微分方程组。

[math]\displaystyle{ \begin{cases} \dot{x}=\sigma(y-x) \\ \dot{y}=rx-y-xz \\ \dot{z}=xy-bz \end{cases} }[/math]

该方程组通过傅里叶展开和截断方法，将大气热对流的非线性PDE简化为三维自治动力系统，用于描述空气流体运动^[4]。

重要地位

洛伦兹方程第一个揭示了确定性系统中的内在随机性，标志混沌理论的诞生。为非线性科学开辟了新的研究方向。

发展历史

问题起源：气象学的无穷维方程（1900-1916年）

洛伦兹方程的起源，可追溯到流体力学中的经典问题——瑞利-贝纳德对流^[5]。这是大气科学中理解环流和对流的基础，数学模型如下：

[math]\displaystyle{ \begin{cases} \nabla \cdot u = 0 \\ \displaystyle \frac{\partial u}{\partial t} + (u \cdot \nabla)u = -\frac{1}{\rho_0}\nabla p + \nu \nabla^2 u + \alpha g (T-T_0) e_z \\ \displaystyle \frac{\partial T}{\partial t} + u \cdot \nabla T = \gamma \nabla^2 T \end{cases} }[/math]

这是一组非线性PDE方程，自由度无限维。

铺路奠基：降维与截断的七阶模型（1962年）

科学家尝试求解方程，实现天气预测。但高维非线性PDE无法直接求解，需简化。Saltzman将方程展开为傅里叶级数并截断，建立了七阶模型LOM(7)，为洛伦兹方程奠定了数学基础。

横空出世：一杯咖啡与洛伦兹方程（1963年）

混沌的发现充满了戏剧性，洛伦兹计算时，将 0.506127 舍入为 0.506，喝完咖啡，发现微小初始偏差产生巨大差别，这是混沌现象的经典发现。

洛伦兹发现LOM(7)中三个是核心，将方程截断为三维ODE。洛伦兹确定性方程产生看似随机的行为，打破了“确定性方程必有确定规则解”的认知。

奇怪吸引子：蝴蝶效应及混沌的诞生（1972-1975年）

1972 年，洛伦兹在 AAAS 年会上，以“蝴蝶效应”形象阐释了初值敏感性，让混沌概念得以普及。 1975 年，李天岩与约克的论文首次给出“混沌”的严格数学定义。洛伦兹方程非周期、初值敏感、轨迹有界不重复特征，完全契合混沌定义，确认为首个经典混沌模型。

数学形式

方程的简化及形式

1）模态分析洛伦兹发现七阶系统中三个模态 X₁,X₂,X₄ 为基模及第一谐波，是对流主模态，慢模态，具有长期动力学行为。其余四个模态 X₃,X₅,X₆,X₇ 为快速衰减模态，对应高次竖直模/水平高次模。

2）方程简化设快模态 [math]\displaystyle{ \dot{x}_i=0 }[/math]（中心流形定理，快变量流形上为常数），解代数方程，将快模态表示为慢模态的函数，代入慢模态方程。 [math]\displaystyle{ \begin{cases} 0 = -a_{33}X_3 + b_{33}X_6 + c_{31}X_1X_5 \implies X_3 = K_3(X_1,X_4,X_2) \\ 0 = -b_{55}X_5 + e_{52}X_2 + f_{51}X_1X_4 \implies X_5 = K_5(X_1,X_4,X_2) \\ 0 = -b_{66}X_6 + e_{63}X_3 + f_{61}X_1X_5 \implies X_6 = K_6(X_1,X_4,X_2) \\ 0 = -a_{77}X_7 + g_{71}X_1X_5 \implies X_7 = K_7(X_1,X_4,X_2) \end{cases} }[/math] 3）方程形式令：X=X₁，Y=X₄，Z=X₂，代入并合并常数，将 Rayleigh-Bénard 对流方程简化为洛伦兹方程，即标准三维自治常微分方程组：

[math]\displaystyle{ \begin{cases} \dot{x}=\sigma(y-x) \\ \dot{y}=rx-y-xz \\ \dot{z}=xy-bz \end{cases} }[/math]

其中 [math]\displaystyle{ \dot{x}=\frac{dx}{dt},\dot{y}=\frac{dy}{dt},\dot{z}=\frac{dz}{dt} }[/math]，都是对时间 t 的导数。这是一组三维、自治、非线性、确定性常微分方程组。

方程的参数及变量

1）方程变量

[math]\displaystyle{ x(t) }[/math]：对流滚动强度，与流速场第一模态振幅成正比，对应对流环流流速大小。
[math]\displaystyle{ y(t) }[/math]：温度差，反映热传递不均匀性，正比于上升暖气流与下沉冷气流温差。
[math]\displaystyle{ z(t) }[/math]：垂直方向温度分布扭曲程度，描述温度剖面偏离线性梯度的畸变。

2）方程参数

[math]\displaystyle{ \sigma }[/math]：普朗特数，动量扩散系数与热扩散系数之比，标准取值 [math]\displaystyle{ \sigma\approx 10 }[/math]。
[math]\displaystyle{ r }[/math]：瑞利数，表征系统加热强度，为分岔参数，标准取值 [math]\displaystyle{ r=28 }[/math]。
[math]\displaystyle{ b }[/math]：几何耗散参数，与对流胞纵横比有关，标准取值 [math]\displaystyle{ b=\displaystyle\frac{8}{3}\approx 2.6667 }[/math]。

方程的拓扑及空间结构

1）拓扑传递及混合性质洛伦兹方程的混沌结构由拓扑混合性严格刻画。拓扑混合性：对任意两个非空开集 U,V⊂X，存在 T>0，当 t>T 时，[math]\displaystyle{ \phi_t(U)∩V≠∅ }[/math]，则流 [math]\displaystyle{ \phi_t }[/math] 在 X 上拓扑混合。

洛伦兹吸引子 [math]\displaystyle{ \Lambda }[/math] 上的连续流是拓扑混合的，是不可分解的单一混沌不变集。

拓扑传递性：1979 年 R.F.Williams 证明^[6]
拓扑混合性：2005 年 Luzzatto 等证明经典参数下为拓扑混合^[7]

2）空间收缩相空间体积随时间指数收缩： [math]\displaystyle{ V(t)=V(0)e^{-(\sigma+1+b)t} }[/math] 所有轨道最终被吸引到低维混沌吸引子上。

方程的数学性质

1）非线性：含 [math]\displaystyle{ xz }[/math]、[math]\displaystyle{ xy }[/math] 二次非线性项，是混沌的数学根源。 2）对称性：[math]\displaystyle{ (x,y,z)\mapsto (-x,-y,z) }[/math] 对称，系统关于 [math]\displaystyle{ z }[/math] 轴对称。 3）耗散性：散度 [math]\displaystyle{ \nabla\cdot F=-\sigma-1-b\lt 0 }[/math]，体积严格指数收缩。 4）解空间：无准周期解、无排斥不动点，所有解收敛到吸引子。

动力学性质

稳定特性

不动点分析

令 [math]\displaystyle{ \dot{x}=\dot{y}=\dot{z}=0 }[/math]，解得三组不动点：

[math]\displaystyle{ C_0(0,0,0) }[/math]：始终存在
[math]\displaystyle{ C_+ }[/math]、[math]\displaystyle{ C_- }[/math]：当 [math]\displaystyle{ r\gt 1 }[/math] 时出现

不动点存在性
序号	参数范围	[math]\displaystyle{ C_0 }[/math]	[math]\displaystyle{ C_+ }[/math]	[math]\displaystyle{ C_- }[/math]
1	[math]\displaystyle{ 0\lt r\lt 1 }[/math]	[math]\displaystyle{ (0,0,0) }[/math]	不存在	不存在
2	[math]\displaystyle{ r=1 }[/math]	[math]\displaystyle{ (0,0,0) }[/math]	不存在	不存在
3	[math]\displaystyle{ r\gt 1 }[/math]	[math]\displaystyle{ (0,0,0) }[/math]	存在	存在

局部线性稳定性（雅克比矩阵）

雅克比矩阵： [math]\displaystyle{ J= \begin{bmatrix} -\sigma & \sigma & 0 \\ r-z & -1 & -x \\ y & x & -b \end{bmatrix} }[/math]

[math]\displaystyle{ r\lt 1 }[/math]：[math]\displaystyle{ C_0 }[/math] 全局稳定
[math]\displaystyle{ r=1 }[/math]：[math]\displaystyle{ C_0 }[/math] 非双曲平衡点
[math]\displaystyle{ r\gt 1 }[/math]：[math]\displaystyle{ C_0 }[/math] 不稳定
临界值 [math]\displaystyle{ r_c=\displaystyle\frac{\sigma(\sigma+b+3)}{\sigma-b-1} }[/math]：发生 Hopf 分岔

全局稳定性

[math]\displaystyle{ r\lt 1 }[/math]：[math]\displaystyle{ C_0 }[/math] 全局稳定
[math]\displaystyle{ r\gt 1 }[/math]：[math]\displaystyle{ C_+ }[/math]、[math]\displaystyle{ C_- }[/math] 仅局部稳定，无全局稳定

稳定性总结

序号	参数范围	[math]\displaystyle{ C_0 }[/math]	[math]\displaystyle{ C_+ }[/math]	[math]\displaystyle{ C_- }[/math]	物理现象
1	[math]\displaystyle{ 0\lt r\lt 1 }[/math]	全局稳定	不存在	不存在	静止无对流
2	[math]\displaystyle{ r=1 }[/math]	非双曲平衡点	不存在	不存在	静止无对流
3	[math]\displaystyle{ 1\lt r\lt r_c }[/math]	不稳定	局部稳定	局部稳定	稳定对流
4	[math]\displaystyle{ r\gt r_c }[/math]	不稳定	不稳定	不稳定	混沌

分叉特性

固定 [math]\displaystyle{ \sigma=10 }[/math]，[math]\displaystyle{ b=\displaystyle\frac{8}{3} }[/math]，以 [math]\displaystyle{ r }[/math] 为分岔参数：

[math]\displaystyle{ r\lt 1 }[/math]：稳定不动点
[math]\displaystyle{ r=1 }[/math]：叉形分岔
[math]\displaystyle{ r\approx 13.926 }[/math]：同宿轨分岔
[math]\displaystyle{ r\approx 24.74 }[/math]：亚临界 Hopf 分岔（进入混沌）
[math]\displaystyle{ r\gt 24.74 }[/math]：混沌 + 周期窗口

分岔路径：稳定原点 → 稳定对流 → 周期对流 → 混沌

轨道特性

1）有序轨道（[math]\displaystyle{ r\lt 24.74 }[/math]）

[math]\displaystyle{ r\lt 1 }[/math]：收敛到原点
[math]\displaystyle{ 1\lt r\lt 13.926 }[/math]：收敛到 [math]\displaystyle{ C_+ }[/math]/[math]\displaystyle{ C_- }[/math]
[math]\displaystyle{ 13.926\lt r\lt 24.74 }[/math]：周期极限环

2）混沌轨道（[math]\displaystyle{ r\gt 24.74 }[/math]）

初值敏感性（蝴蝶效应）
非周期、有界
蝴蝶形双涡卷奇怪吸引子
分形维数≈2.06

解的遍历性与混合性

遍历性：轨道稠密覆盖吸引子，时间平均=空间平均
混合性：任意区域均匀扩散，初始相关指数衰减

动力性质小结

[math]\displaystyle{ r }[/math] 取值	动力学行为	稳定对象	轨道特征
[math]\displaystyle{ r\lt 1 }[/math]	纯导热稳态	原点[math]\displaystyle{ C_0 }[/math]	静止稳定
[math]\displaystyle{ 1\lt r\lt 13.926 }[/math]	稳定对流	[math]\displaystyle{ C_+ }[/math]、[math]\displaystyle{ C_- }[/math]	收敛至不动点
[math]\displaystyle{ 13.926\lt r\lt 24.74 }[/math]	周期对流	极限环	周期振荡
[math]\displaystyle{ r\approx 24.74 }[/math]	临界失稳	无	向混沌过渡
[math]\displaystyle{ r\gt 24.74 }[/math]	混沌+周期窗口	洛伦兹吸引子	蝴蝶轨道、初值敏感

混沌水车

混沌水车的由来

1972年，麻省理工学院 Willem Malkus 和 Lou Howard 建造第一台混沌水车，将洛伦兹方程转化为可观测机械系统^[8]。随注水流量增大，水车依次呈现：静止 → 单向旋转 → 无规则混沌转动。

混沌水车与动力学参数关系

物理结构：圆形转轮 + 周向水槽 + 顶部注水 + 底部漏水 + 自由转动。

关键参数：

[math]\displaystyle{ \omega(t) }[/math]：角速度
[math]\displaystyle{ m(\theta,t) }[/math]：水的质量分布
[math]\displaystyle{ Q(\theta) }[/math]：注水流量（对应洛伦兹方程 [math]\displaystyle{ r }[/math]）
[math]\displaystyle{ K }[/math]：泄漏率（对应系统耗散）
[math]\displaystyle{ I }[/math]：转动惯量
[math]\displaystyle{ v }[/math]：旋转阻尼系数

混沌水车运动到洛伦兹方程推导

1）质量守恒方程 [math]\displaystyle{ \frac{\partial m}{\partial t}=Q-Km-\omega\frac{\partial m}{\partial \theta} }[/math]

2）力矩平衡 [math]\displaystyle{ I\dot{\omega}=-v\omega + gr\int_0^{2\pi}m(\theta,t)\sin\theta d\theta }[/math]

3）傅里叶展开后，仅保留一阶模态，得到封闭系统： [math]\displaystyle{ \begin{cases} \dot{a}_1 = \omega b_1 - K a_1 \\ \dot{b}_1 = -\omega a_1 - K b_1 + q_1 \\ \dot{\omega} = \frac{1}{I}\left(-v\omega + \pi g r a_1\right) \end{cases} }[/math]

4）变量代换后，严格等价于洛伦兹方程： [math]\displaystyle{ \begin{cases} \dot{x}=\sigma(y-x) \\ \dot{y}=rx-y-xz \\ \dot{z}=xy-bz \end{cases} }[/math]

混沌水车物理演化

[math]\displaystyle{ r\lt 1 }[/math]：静止
[math]\displaystyle{ 1\lt r\lt r_0 }[/math]：单向稳定旋转
增大[math]\displaystyle{ r }[/math]：周期正反转
[math]\displaystyle{ r\approx 28 }[/math]：混沌无规则转动
更高[math]\displaystyle{ r }[/math]：混沌带、周期窗口、激变

科学意义与跨域应用