2. Lorenz系统数学模型建立

Lorenz系统最初来自气象学家爱德华·洛伦兹对二维热对流（Rayleigh–Bénard对流）的简化研究。他把流体层中的温度场与速度场用傅里叶模态截断，只保留少数几个模态，再经过无量纲化与变量重标定，最终得到一组仅含三个状态变量、三个参数的一阶常微分方程。这组方程之所以成为混沌研究的经典模型，不在于它精确复现了真实大气，而在于它以极低的维度保留了“确定性系统出现不可预测行为”的本质机制。

一、从物理模型到三个状态变量

在推导之前先明确变量约定，全文与后续仿真保持一致：x 表示对流环流的强度（与速度场第一模态成正比），y 表示上升流与下降流之间的温度差（水平温度模态），z 表示垂直温度剖面对线性分布的偏离量。三者都是无量纲量，时间 t 也是无量纲时间。参数方面：σ 为普朗特数，反映流体黏性扩散与热扩散之比；ρ（文献中常记作 r）为瑞利数相对临界瑞利数的比值，衡量加热强度；β 为几何因子，与对流胞的水平/垂直尺度比相关。

对原始偏微分方程做模态截断后，代入无量纲变量并消去中间量，可整理出如下标准形式：

dx/dt = σ · (y - x)
dy/dt = x · (ρ - z) - y
dz/dt = x · y - β · z

这三个方程必须联立求解，任何一个方程的输出都是另外两个方程的输入，这正是“耦合”的含义。下面这张图把变量、参数与方程之间的依赖关系画清楚，便于在写代码前先建立结构认识。

二、逐项解读方程结构

第一式 ẋ = σ(y − x)。 它描述环流强度如何变化：y 与 x 之差构成驱动力，σ 决定响应快慢。若 y > x，x 增大；若 y < x，x 减小。这一式本身是线性的，但它把 y 的信息注入 x，使 x 成为后续两个非线性项的“燃料”。

第二式 ẏ = x(ρ − z) − y。 展开为 xρ − xz − y。其中 xρ 是加热驱动项，ρ 越大，温度差被拉得越强；−xz 是关键的非线性负反馈：当垂直温度偏离 z 增大时，它会抑制 y 的增长；−y 是线性耗散。正是 xz 这一乘积项，让系统在 ρ 超过临界值后无法稳定在静止解。

第三式 ż = xy − βz。 xy 表示环流强度与水平温度差的乘积对垂直温度偏离的“加热”贡献，βz 表示几何尺度决定的耗散。当 x、y 同号且较大时，z 被推高；z 升高后又通过第二式的 −xz 反过来压制 y，形成闭合的反馈回路。

三、非线性来源与参数作用

三个方程中只有两处乘积：xz 与 xy。它们使方程组成为二次非线性系统，也是全部复杂行为的根源。线性系统不可能产生混沌，因为线性系统的解只能是指数增长、衰减或振荡的组合；而二次项允许轨迹在相空间中“折叠”，使相邻初值以指数速度分离，这就是对初始条件敏感性的几何基础。

参数的作用可以这样理解：σ 只影响 x 的响应速度，不改变平衡点的位置；β 决定 z 的衰减快慢，影响轨迹在 z 方向的收缩程度；ρ 是最关键的分岔参数。当 0 < ρ < 1 时，原点 (0,0,0) 是全局稳定的，对应无对流状态；当 ρ > 1 时原点失稳，出现一对对称平衡点 C± = (±√(β(ρ−1)), ±√(β(ρ−1)), ρ−1)，对应顺时针与逆时针两种稳定对流；继续增大 ρ，C± 也失稳，系统进入混沌。经典取值 σ=10、ρ=28、β=8/3 时，轨迹在两个不稳定平衡点附近交替盘旋，形成著名的“蝴蝶”吸引子。

四、为数值仿真做的准备

建立数学模型之后，下一步是把它变成可计算的初值问题。需要确定三件事：初值 (x₀, y₀, z₀)、参数 (σ, ρ, β)、积分方法与步长。以经典参数为例，取初值 (1, 1, 1)，用四阶龙格–库塔法、步长 h=0.01 积分 50 个时间单位，即可观察到轨迹先短暂过渡、随后进入两个“翼”之间不规则切换的过程。判断是否进入混沌区间，可以观察轨迹是否既不收敛到 C±、也不发散，同时相邻初值（如把 x₀ 改为 1.001）的轨迹在若干时间后明显分离。

需要强调的是，上述参数与初值属于教学示例设定，具体数值结果依赖积分步长与实现细节，实际运行时应当以自己程序的输出为准。建模阶段真正要固化下来的是方程形式、变量含义、参数约定与符号方向，这些一旦确定，后续所有仿真、绘图与结果分析都必须与之保持一致，避免因符号或坐标约定不同而误判轨迹形态。