同步:自己对上拍子
同步是涌现里最干净的一个例子:每个振子只知道邻居(或者整体的平均相位)在哪儿,没有谁在指挥,可是耦合一旦超过某个门槛,几千个各走各的东西会突然走到一起。而且这件事有一个精确的临界点,可以算出来,也可以量出来。
本页实验台
LAB 01
Kuramoto 相变
把耦合强度 K 从 0 拉到 5,看序参量 r 在临界耦合处从 0 抬起来
问题
东南亚的萤火虫能整片同时闪,心脏里上亿个细胞能同时收缩,电网上所有发电机保持同一个频率。这些振子各有各的固有频率,为什么会突然对上拍子?而且这个"突然"是一个真正的相变。
模型
Kuramoto 模型:每个振子被所有其他振子拉向平均相位。
dθi/dt = ωi + (K/N) Σj sin(θj − θi); r eiψ = (1/N) Σj eiθj
ωi 是固有频率,K 是耦合强度,序参量 r ∈ [0,1] 量同步程度。固有频率服从标准差 σ 的高斯时,临界耦合 Kc = σ√(8/π)——这是一个可以当场核对的解析结果。当前的序参量 r—
理论临界耦合 σ√(8/π)—
K / Kc—
现在是什么状态—
每个振子都有自己的固有频率,有快有慢。耦合项做的事很简单:如果整体的平均相位跑在你前面,就加速;落在你后面,就减速。仅此而已。可是把 K 从 0 慢慢拖大,右边那个圆上的点会从均匀散开突然收成一束——这中间有一个明确的相变点。Kuramoto 在 1975 年把这个临界值算了出来:固有频率服从标准差为 σ 的正态分布时,Kc = σ·√(8/π) ≈ 1.596σ。实测和它对得极好:σ=0.5 时 Kc=0.798,K/Kc = 0.5 处 r 只有 0.087(就是有限规模的噪声底),K/Kc = 1 处 r 跳到 0.276,K/Kc = 2 处已经是 0.92。注意这个转变的形状:临界点之前 r 基本贴着 0,之后按 √(K−Kc) 上升——这是一个连续相变,和磁铁在居里点的磁化曲线是同一类。圆上那些颜色也值得看。橙色是固有频率接近平均值的振子,它们最先被锁住;蓝色是跑得太快或太慢的,它们在部分同步区里始终锁不住,绕着圆一直转。这就是"部分同步"的含义:不是所有人都加入,而是频率足够接近的那一部分先形成核心,然后核心越来越大。(现实例子:拍手会自发变整齐,但总有几个人跟不上;电网里发电机要锁在 50Hz,锁不住的那一台会被切出去。)
LAB 02
萤火虫与脉冲耦合
每次闪光把邻居往前推一点,看一屏乱闪几十秒后同时亮起
问题
上一台的振子是连续相互作用的,可萤火虫只在"闪"的那一瞬间影响别人。这种脉冲式的耦合还能同步吗?1990 年 Mirollo 和 Strogatz 证明:只要响应函数是凸的,几乎所有初值都会同步。
模型
积分-发放振子:相位线性增长,到阈值发放并把邻居往前推一点。
dφ/dt = 1; φ = 1 时发放并归零,同时所有邻居 φj ← min(1, f−1(f(φj) + ε))
f 是把相位映到"状态变量"的响应函数,ε 是耦合强度。关键条件是 f 凹(边际响应递减)——这保证任意两个振子之间的相位差单调收缩,从而全局同步。已经跑了多少步—
最近的同闪峰值—
相位的离散程度—
同步了吗—
这个模型和上一个不一样:振子之间不是连续耦合,而是靠脉冲。每只萤火虫的"充能"匀速上升,到顶就闪一下然后归零;而每当有别人闪光,你的充能就往前跳一点点。规则就这一条。按播放,看着一屏乱闪的点在几十步之内聚成一团团,最后整齐地一起闪。Mirollo 和 Strogatz 在 1990 年证明了一个很强的结论:只要耦合是正的(脉冲让别人提前而不是延后),而且充能曲线是凸的,几乎所有初始条件最终都会完全同步。不需要中心节点,不需要全局信息,每只萤火虫只知道"刚才有人闪了"。把推动量 ε 拖到 0:永远不会同步(各闪各的);拖大:同步得更快,但如果太大会出现"跳过头",反而出现来回摆动。这个机制的应用面比听起来广得多:无线传感器网络的时钟同步用的就是它(每个节点广播一个脉冲,收到就调表),优点是完全分布式、没有单点故障、节点加入退出都不影响。心脏的窦房结细胞、肠道的起搏细胞、甚至一排节拍器放在同一块活动板上,用的都是同一个数学。
LAB 03
耦合映射格子
把一排 Logistic 映射耦合起来,看同步簇、行波和时空混沌怎么切换
问题
把一排各自混沌的系统耦合起来,结果是更混沌,还是反而同步?答案两种都有,取决于耦合强度——而中间还会出现同步簇、行波这些两者都不是的状态。
模型
耦合映射格子:局部混沌映射 + 扩散耦合。
xit+1 = (1−ε) f(xit) + (ε/2)[ f(xi−1t) + f(xi+1t) ], f(x) = a x(1−x)
ε 是耦合强度,a 是 Logistic 参数。ε 小是时空混沌,ε 大是完全同步,中间有同步簇与行波。这是最简单的能同时展示"空间自由度 + 时间混沌"的模型,常被当作湍流的玩具模型。相邻格子的平均差异—
最大值与最小值之差—
有多少条"分界线"—
处在什么状态—
把一排混沌的逻辑斯蒂映射用最弱的方式连起来:每个格子的下一步 = (1−ε)×自己的迭代 + ε/2×左右邻居的迭代。单独一个格子在 a=3.9 时是完全混沌的,连起来以后画面会随 ε 走过三个完全不同的阶段。ε 很小(0.05 附近):时空混沌。每个格子基本各走各的,时空图上是细密的雪花点,看不出任何结构。ε 中等(0.2–0.3):出现同步簇。相邻的一片格子锁在一起,时空图上出现纵向的条带和清晰的分界线,而分界线本身还会缓慢移动、合并。ε 很大(0.4 以上):整体完全同步。所有格子走同一条混沌轨道——注意这时系统仍然是混沌的(时间上不可预测),但空间上完全一致。这两件事是分开的:混沌说的是时间,同步说的是空间。耦合映射格子是研究时空混沌最便宜的模型。它在物理上对应耦合的激光阵列、化学上对应催化剂表面的反应波、生态上对应通过迁移相连的若干个种群斑块——最后这个有一个很实际的结论:适度的迁移会让各斑块的种群数量同步,而同步意味着它们会同时崩溃,整体灭绝的风险反而上升。保护生物学里因此有"不要让栖息地连得太紧"这样的建议。
LAB 04
同步的代价与风险
电网、人群走桥、金融市场:同一套模型里,同步既是功能也是灾难
问题
同步通常被当成好事——可 2000 年伦敦千禧桥开放当天就因为行人同步而剧烈摇晃,2008 年金融危机里各家银行的风险模型同步失效。什么时候同步是功能,什么时候是灾难?
模型
同一套耦合振子方程,把"耦合"换成不同的物理量就得到不同的系统。
桥:行人步频 ↔ 桥面横向位移(正反馈)
电网:发电机相角 ↔ 潮流(需要同步)
市场:风险模型 ↔ 抛售行为(同步 = 流动性消失)
共同结构是:个体为了自身利益去跟随群体,而跟随行为本身放大了群体信号。区别只在于这个放大是被设计进去的(电网)还是意外的(桥、市场)。千禧桥最后靠加阻尼器把耦合压下去解决。电网:发电机相角 ↔ 潮流(需要同步)
市场:风险模型 ↔ 抛售行为(同步 = 流动性消失)
临界人数—
20 人时的振幅—
320 人时的振幅—
放大了多少倍—
2000 年 6 月 10 日,伦敦千禧桥开放,当天就因为剧烈横向摇晃而关闭。调查结论是一个漂亮的正反馈:桥轻微晃动 → 行人为了保持平衡不自觉地同步调整步伐 → 同步的步伐给桥更大的推力 → 桥晃得更厉害。这一页把这个回路做成模型:N 个行人的步伐相位耦合到桥的一个横向模态上。关键在那条振幅曲线:它不是随人数平滑上升的,而是在某个临界人数处突然起飞。下面那张同步度曲线在同一个位置跳起来——两件事是同一个相变。在临界人数以下,行人的步伐互不相关,他们的推力互相抵消,桥几乎不动;越过临界点,一点点同步就会带来更大的晃动,而更大的晃动又带来更多同步。这就是为什么千禧桥的问题在试运行时没被发现:人少的时候一切正常,而开放日涌上来的人数正好越过了那条线。把阻尼滑杆拖大再看:临界人数明显右移,甚至在这个范围里不再出现。这正是最后的解决方案——桥上装了 37 个阻尼器,不是让桥更硬,而是让它更"吸能"。同样的机制出现在很多地方:金融市场里"大家都看同一个指标"会放大波动;电网里发电机的同步失稳会级联;人群在密集场所的步调一致会产生远超个体的合力。同步在这些场合不是功能,是故障模式。
这一页想说的
同步是一个有精确临界点的涌现现象:Kuramoto 模型给出 Kc = σ√(8/π),实测与理论几乎重合;脉冲耦合即使没有任何中心节点也能保证收敛;耦合映射格子说明"混沌"和"同步"是两个独立的维度。最后一点最值得带走:同步既能是功能(心跳、时钟、发电机),也能是灾难(桥、电网、市场),而且它出现得非常突然。