Math Labs

解法:塔克利用区间算术给出的严格计算机辅助证明(1999年)

第 1/8 步:Smale 第14问题:洛伦兹的蝴蝶真的是稳健的吸引子吗?
通俗地说

1963年,气象学家爱德华·洛伦兹在早期计算机上模拟了一个被极度简化的大气对流模型,发现了一件令人震惊的事:轨道从未稳定到一个不动点或一个重复的循环,而是描绘出一个无限、永不闭合的双螺旋形状,形似蝴蝶的翅膀,而起点上极其微小的变化会导致完全不同的长期路径。这正是广为人知的「蝴蝶效应」概念的由来。

但计算机模拟无论画面多么震撼,都不是证明:在混沌系统中,浮点舍入误差会不断累积,所以屏幕上的图案可能只是被打扮得像真正混沌的数值噪声,而不是一个真实的数学对象。1998年,史蒂夫·斯梅尔把这个问题列为他为21世纪列出的数学挑战清单中的第 1414 题:请严格证明洛伦兹的图案是真实存在的。

x˙=σ(y−x),y˙=x(ρ−z)−y,z˙=xy−βz,(σ,ρ,β)=(10,28,8/3)\dot{x} = \sigma(y - x), \quad \dot{y} = x(\rho - z) - y, \quad \dot{z} = xy - \beta z, \quad (\sigma, \rho, \beta) = (10, 28, 8/3)
详细分析

Lorenz 方程 x˙=σ(y−x)\dot{x} = \sigma(y - x), y˙=x(ρ−z)−y\dot{y} = x(\rho - z) - y, z˙=xy−βz\dot{z} = xy - \beta z 源自对 Rayleigh–Bénard 对流 Fourier 展开的粗暴截断——Lorenz 只保留三个模——在 (σ,ρ,β)=(10,28,8/3)(\sigma, \rho, \beta) = (10, 28, 8/3) 处一切数值模拟都显示著名的双叶吸引子。但模拟不是证明:浮点轨道并不严格地“影子”任何真实轨道,且所观测的行为在逻辑上仍可能与一条天文长度级的稳定周期轨道相容,或与只存在于舍入误差累积之内的吸引子相容。

Smale 把此问题列为他 1998 年“下世纪数学问题”的第 14 题:经典参数下的 Lorenz 系统是否具有一个鲁棒的奇怪吸引子——即在向量场的小扰动下依然存续?古典证明的障碍是结构性的。原点是位于吸引子内部的鞍型平衡点,因此系统是奇异双曲而非一致双曲:轨道可以任意接近鞍点,在那里两方向的收缩与一方向的扩张相遇,且向任何截面的回归时间在稳定流形附近无界增长。

Tucker 1999 年的证明(完整细节见 “A rigorous ODE solver and Smale's 14th problem”, Found. Comput. Math. 2002)给出了肯定答案:他把相空间分为原点邻域(用精确分析处理)与其补集(用误差界精确到机器精度的区间算术计算处理)——这是经验证数值方法首次解决一个真正未决的动力系统问题。

本步骤中的术语
奇异吸引子
相空间中一个有界区域,附近的轨道会被吸入其中并且永不离开,但其上的运动是混沌的(对初始条件敏感)而非周期性的;它的截面通常具有分形的、非整数维的结构。
鞍型平衡点
动力系统中的一个不动点,轨道沿某些方向被吸引,沿另一些方向被排斥,就像放在马鞍上的球:前后轻推是稳定的,左右轻推则不稳定。
本步骤用到的知识