CurvG
登录

第 3 个 · 微分方程 · 混沌

方程完全确定,为什么结果却无法预测?

三条描述空气对流的方程,让一个点画出一只永不重复的"蝴蝶"。两个只差十万分之一的起点,十几秒后就分道扬镳。拖动换个角度,调一调参数,看混沌怎样出现。

为什么?

讲解

混沌不是随机:每一步都由方程精确决定,但起点上再小的误差都会被指数放大,所以只能预测形状,预测不了位置。

01三条方程

气象学家洛伦兹把空气受热翻滚的过程简化成三个量:xx 是翻滚的快慢,yy 和 zz 描述温度分布偏离均匀的程度。它们随时间的变化是

x˙=σ(y−x),y˙=x(ρ−z)−y,z˙=xy−βz\dot x=\sigma(y-x),\quad \dot y=x(\rho-z)-y,\quad \dot z=xy-\beta z

这里取他当年用的参数 σ=10\sigma=10、β=8/3\beta=8/3、ρ=28\rho=28。方程里没有任何随机项:给定起点,以后的每一刻都完全确定。

画面上的轨迹用四阶龙格–库塔法算出,时间步长 0.005。

02蝴蝶效应

在这个系统里,两个相邻起点的距离平均每过约 0.77 个时间单位就翻一倍(最大李雅普诺夫指数约 0.906)。指数增长非常快:故事里两个点起初只差 10−510^{-5},大约 16 个时间单位后就相距超过 1,跟着便是完全不同的两条路。

洛伦兹 1961 年在计算机上重算一段天气模拟,把一个数从 6 位小数四舍五入成 3 位,结果几周后的"天气"面目全非。他在 1972 年的一次演讲题目里问:巴西一只蝴蝶扇动翅膀,会不会在得克萨斯引起一场龙卷风?"蝴蝶效应"由此得名。

这也是天气预报很难准确超过两周的原因之一:初始观测再精确,误差也会被这样放大。

03确定,却不可预测

混沌有三个特征同时成立:

  • 确定:没有随机性,同样的起点永远给出同样的结果;
  • 敏感:起点的微小差别会指数放大;
  • 有界且不重复:轨迹永远困在有限的区域里,却从不回到走过的路上。

所以单个点的位置无法长期预测,但三千个点最后铺出的形状每次都一样。这个形状叫奇异吸引子,它是一个分形,维数约为 2.06:比曲面"多"一点点。

04调大 ρ,蝴蝶才出现

ρ\rho 代表加热的强弱(σ、β 保持不变):

  • ρ<1\rho<1:加热太弱,所有点都回到原点,空气不动;
  • 1<ρ<24.061<\rho<24.06:点最终停在两只翅膀的中心之一,空气稳定地朝一个方向翻滚(ρ\rho 在 13.93 到 24.06 之间时,可能要先乱跳很久才停下);
  • 24.06<ρ<24.7424.06<\rho<24.74:蝴蝶和两个静止点同时存在,停在哪里取决于起点;
  • ρ>24.74\rho>24.74:两个中心都不稳定了,只剩蝴蝶。

两只翅膀的中心在 (±β(ρ−1), ±β(ρ−1), ρ−1)\left(\pm\sqrt{\beta(\rho-1)},\ \pm\sqrt{\beta(\rho-1)},\ \rho-1\right),ρ=28\rho=28 时约为 (±8.49,±8.49,27)(\pm8.49,\pm8.49,27)。

常见问题

蝴蝶真的能引起龙卷风吗?

这是个比喻。它说的是:在混沌系统里,任何微小的差别都可能在很久以后变成巨大的差别,所以无法追溯"哪一只蝴蝶"负责,也无法长期预测。

既然误差会放大,屏幕上的轨迹还可信吗?

单条轨迹算久了,舍入误差也会被放大,和"真实"轨迹分开。但吸引子的形状是可靠的:2002 年 Warwick Tucker 用计算机辅助证明,确认这个奇异吸引子确实存在(这是数学家斯梅尔列出的 18 个问题中的第 14 个)。

混沌和随机有什么区别?

随机是规则里本来就有不确定性;混沌的规则完全确定,只是对起点太敏感。同样的起点重算一遍,混沌系统给出一模一样的结果,随机过程不会。

为什么形状像一只蝴蝶?

两只"翅膀"各绕着一个中心,它们是方程的两个平衡点,但都不稳定:点被推着绕中心转圈、越转越远,远到一定程度就被甩到另一边。

课本里的这部分

资料: Lorenz system — Wikipedia · Butterfly effect — Wikipedia