01 它问的是什么问题(The Question It Poses)
1960 年代初的共识大致是这样:一个有耗散的系统——也就是会把能量磨掉的、现实里的任何东西——跑久了总会安定下来。要么停在某个状态不动(不动点),要么落进一个周而复始的循环(周期轨道)。除此之外没有第三种下场。
洛伦兹当时在做的事更朴素:他想知道气象学界流行的"用统计方法做长期预报"到底靠不靠谱。为此他需要一个小到能在桌上机器里跑、又保留了大气非线性本质的例子。他从 Barry Saltzman 一组描述热对流 (thermal convection,下面热、上面冷,流体自己转起来的那件事)的方程里砍掉了绝大部分,只留三个。
于是问题变成:这样一个完全确定、只有三个变量、还带耗散的系统,能不能既不停下、也不重复?
02 规则本身(The Rules)
dx/dt = σ (y − x) dy/dt = x (ρ − z) − y dz/dt = x y − β z
x ≈ 对流翻滚的强度:正负号表示往哪个方向转。
y ≈ 上升的那股流和下沉的那股流之间的温差。
z ≈ 垂直方向的温度分布偏离"从下到上均匀变化"的程度。
σ 是流体的黏性与导热之比(普朗特数 Prandtl number),ρ 衡量下面加热得多厉害(瑞利数 Rayleigh number 的相对值),β 是容器形状带来的一个几何常数。
洛伦兹用的那组数:σ = 10,ρ = 28,β = 8/3 。整套系统的所有名声都来自这一组数——换一组,故事完全不同(见最后一节)。
非线性只藏在两个地方:xz 和 xy 这两个乘积项。整套方程里没有任何随机成分,你给同一个初值,它永远给同一条轨迹。
下面热、上面冷,流体自己转起来
冷(上板)
热(下板)
一对对流卷
加热越强,翻滚越猛 — 这就是 ρ 在调的东西
x
翻滚的强度与方向
y
上升流与下沉流的温差
z
温度分布偏离均匀的程度
三个数,就是这套系统的全部状态
方程是把一层受热流体的行为压缩成三个数之后的样子。它保留了非线性,扔掉了几乎所有别的东西。
03 跑起来会看到什么(What You See When It Runs)
轨迹绕着一个中心转几圈,忽然跳到另一边,绕几圈,再跳回来。跳与不跳没有规律:跳之前的那几圈,和不跳之前的那几圈,看不出区别。
把这条轨迹画在以 x、y、z 为坐标的空间里 (相空间 → 参考页) ,就是那只著名的蝴蝶。
一条轨迹跑 70 个时间单位(x-z 投影)
x →
z →
两叶
= 对流的两个转向
线永不自交、
永不重复、
也永不离开
这不是一条线画满了一个面:放大看,它是无穷多层薄片叠起来的分形集合,维数约 2.06。
三件值得记住的事:
一、它同时是收缩的和发散的。 系统有耗散,所以相空间里的一小团体积会不断被压扁,最后体积等于零;可轨迹又对初值敏感,一小团里的点会被拉开。压扁 + 拉开 + 折回来,反复无穷次,得到的就是那个分形结构——这也是"奇怪吸引子 "(strange attractor)里"奇怪"二字的来历。
二、误差每约 0.77 个时间单位翻一倍。 这组参数下最大李雅普诺夫指数(Lyapunov exponent)λ ≈ 0.906,正数就意味着混沌。
三、两叶之间的切换序列不重复。 把"在左叶"记成 0、"在右叶"记成 1,得到的 01 串没有周期。
同一条轨迹,只看 x 随时间
右叶
左叶
有时绕一圈就翻回去,有时连绕好几圈 — 24 个时间单位里翻了 13 次,看不出周期
时间 →
横过中线一次 = 对流翻了个向。什么时候翻,事先看不出来。
04 它解释了现实中的什么(What It Explains)
严格地说,它没有"解释"任何一个具体现象——它推翻了一个信念。1963 年之前,"确定的规则 + 有耗散 ⇒ 最终归于安定"被当成常识;洛伦兹给出了三行方程作为反例。非周期不需要外部的随机来源,系统自己就能一直不重复下去。
它留下的两个概念比方程本身走得更远:奇怪吸引子 (乱是有形状的),以及可预测地平线 (能预测多久是可以算出来的一个量,不是态度问题)。今天集合预报之所以那样做,源头就在这里。
另外,它不只活在纸上。麻省理工的 Willem Malkus 做过一个混沌水轮 :一个转轮,边上挂着一圈漏水的小桶,上方往下浇水。桶接满了就往下坠、水又漏掉,轮子于是时而顺转、时而逆转,切换没有规律。这个装置的运动方程和洛伦兹方程是同一组。单模激光的方程(Haken,1975)也与它同构。
它不能解释什么(What It Cannot Explain)
它不是天气模型。 三个变量,砍掉了水汽、地形、辐射、地球自转——几乎所有东西。它能说明"天气为什么难预报",不能用来算任何一次真实的对流。洛伦兹本人对这一点讲得比后来的科普都清楚。
离开那组参数,它就不混沌了。 ρ < 1 时系统安安静静停在不动点(不对流);1 < ρ < 24.74 时它落到两个稳定的定常对流状态之一——完全可预测。混沌只出现在 ρ 更大的一段里,而 28 是洛伦兹挑的。"洛伦兹系统是混沌的"这句话,省略了"在这组参数下"。 任何拿它讲道理的场合,这个省略都要补回来。
它不能预测任何一次切换。 能给的是统计:在左叶的时间比例、驻留时长的分布。哪一圈会跳,属于这个模型明确拒绝回答的问题。
方程同构不等于系统同构。 水轮和激光确实服从同一组方程,但这是因为它们被抽象到了同一个骨架上。看到某个现实系统"有点像蝴蝶"就说它是洛伦兹系统,是把类比当推导——要主张这件事,得先把变量和参数一一对上,并给出可证伪的预言。
连"蝴蝶效应"本身,在真实大气里都还没定论。 模型里任意小的扰动都会被放大;但真实大气存在有限尺度效应,一个真正微小的扰动能否一路上传到天气尺度,学界仍有争论(Palmer 等人所说的"真正的蝴蝶效应"与这幅经典图像并不相同)。模型给的是一个机制,不是一张已经兑现的支票。