Lorenz Chaotic Dynamics and Attractors
洛伦茨系统（Lorenz System）是由气象学家爱德华·洛伦茨（Edward Lorenz）于 1963 年提出的三维非线性自治微分方程组。该系统最初用于模拟大气对流，其简化版本却展现出了极其复杂的混沌行为 (Chaotic Dynamics)与经典的蝴蝶效应 (Butterfly Effect)。
洛伦茨系统的常微分方程定义如下：
$$\begin{aligned} \frac{dx}{dt} &= a(y - x), \\[0.3em] \frac{dy}{dt} &= rx - y - xz, \\[0.3em] \frac{dz}{dt} &= xy - bz, \end{aligned}$$ 其中，$x, y, z$ 代表系统的状态变量（分别对应流体翻转速率、水平温差和垂直温度梯度），而参数 $a$ (Prandtl 数)、$r$ (Rayleigh 数) 和 $b$ (几何比例因子) 则决定了系统的动力学相图结构。
当参数设定为经典的 $a = 10, b = 8/3, r = 28$ 时，系统展现出极度敏感的初值依赖性。下表总结了在不同 Rayleigh 数 $r$ 下系统的动力学演化特征：
| Rayleigh 数 ($r$) | 相空间形态 | 系统稳定性分析 |
| $0 < r < 1$ | 原点不动点 $C_0(0,0,0)$ | 全局渐近稳定（不发生对流） |
| $1 < r < 24.74$ | 一对对称不动点 $C_1, C_2$ | 稳定焦点（形成稳态对流） |
| $r > 24.74$ | 奇异吸引子 (Strange Attractor) | 不稳定且存在混沌双翼轨道 |
数值积分通常采用四阶 Runge-Kutta 算法（RK4），其离散迭代表达如下：
$\mathbf{x}_{n+1} = \mathbf{x}_n + \frac{\Delta t}{6} (k_1 + 2k_2 + 2k_3 + k_4)$
即便积分步长取值较小（例如 $\Delta t = 0.001$），相轨迹也会在两个蝴蝶翅膀（吸引子双叶）之间进行不可预测的随机跳跃。
洛伦茨吸引子不仅揭示了确定性系统中内建的随机性，更奠定了现代非线性动力学与混沌理论的基础。更有趣的是，可以通过以下过程探索系统的复杂相变现象：
- 通过调整控制台中的 r (Rayleigh) 参数，观察系统从稳定焦点向混沌状态分岔（Bifurcation）的过渡过程。
- 检查在不同 Zoom 缩放比例下吸引子的自相似分形结构：${\color{#2f54eb}D_H \approx 2.06 \pm 0.01}$。