专著选读二 · NEW IN THE REVISED EDITION

第四十二章 最近邻交手:SDE数学与五个先行框架的可裁决差异

修订版新增。把本书的每一块承重构件放到最近的先行者旁边:混沌边缘、非局域耦合振子、模拟退火、Poincaré的表象空间、无量纲群家族——承认重述,指认剩余,并为每处剩余给出一条判本书输的预测。

方法论说明: 本章是修订版新增章,等级为A级(对已有文献的事实陈述)+C级(本章提出的五条对立预测)。它的功能不是宣传,而是清账:把本书的核心构件逐一放到已有文献的最近邻旁边,说明哪些是重述、哪些是剩余,并为每一处剩余给出一条可以判它输的预测。凡本章确认为「已被占位」的,第一版中相应的优先权措辞一律作废,作废清单列于§42.7。

42.1 为什么这一章必须存在

一个新框架的价值不由它自己声称的高度决定,而由它与最近的先行者之间的可裁决差异决定。这是一条比礼貌更硬的要求:若一个命题能被某个已有框架无损重述,则无论它的语言多新、跨域多广,它在发生意义上都没有产生新结构——它只是换了一套记号。

第一版在这一点上有系统性缺陷。全书四十条参考文献集中在计算数学(网格生成、有限元、非局域微积分)与哲学(Peirce、Whitehead、Lakoff-Núñez)两端,而本书赖以立论的三块中间地带——复杂性科学的临界性研究、非局域耦合振子的分岔理论、组合优化中的劣解接受策略——一条文献也没有。这不是谦虚,是漏账。漏账的直接后果是:读者无法判断「介生态」是一个新概念还是一个旧概念的新名字。

本章补上这笔账。方法是固定的四步:(一)陈述最近邻的原始主张,尽量用它自己的语言;(二)承认本书与它重合的部分,逐条列出;(三)指认剩余——本书说了而它没说的那一句;(四)为剩余给出一条对立预测,写明什么样的观测会判本书输。第四步是这一章的重心:没有第四步,前三步只是文献综述。

42.2 交手一:介生态与「混沌边缘」

最近邻。 Langton在1990年提出,元胞自动机的计算能力与一个统计参数λ(规则表中映射到非静止态的比例)相关,并主张能够进行复杂计算的规则聚集在「有序与混沌之间的相变附近」——即混沌边缘。同期与随后的工作把这一图景推得更远:Wolfram的第四类元胞自动机对应长寿命、非周期也非随机的结构;Bak、Tang与Wiesenfeld的自组织临界性给出了系统自发演化到临界点的机制;Kauffman把它用于生命起源与基因调控网络;Beggs与Plenz在皮层切片上观测到服从幂律的神经元雪崩,Chialvo等人由此提出临界脑假说。

重合部分(必须承认)。 本书第六章的三态划分——秩序态、介生态、混沌态,以及「复杂性与自组织能力在两者之间最大」的判断——与混沌边缘图景在结构上同构。第八章用N_SDE≈0.14去对应临界脑假说,实际上是在引用这一传统的结论。介生态不是本书发明的现象,它是一个已被反复观测、已被多学科命名的现象。凡第一版中暗示介生态为本书首次辨识的措辞,一律作废。

剩余(本书多说的那一句)。 混沌边缘传统的关键弱点,是它的控制参数依附于系统的具体实现。λ是元胞自动机规则表的统计量,换一个系统就没有λ。更严重的是,Mitchell、Hraber与Crutchfield在1993年的再检验中证明,λ并不可靠地预测计算能力,λ_c也不是普适常数——它依赖于在规则空间中所取的路径。因此混沌边缘至今是一个描述性图景,不是一个可外推的判据:你能在事后说某系统处在边缘,却不能在事前算出边缘在哪。

本书的剩余正在这里:N_SDE = r·ω/(v·T_c) 不是规则表的统计量,而是由四个可独立测量的物理量组成的无量纲群——纠缠半径、重叠度、传播速度、系统特征时间。它对系统的实现方式不敏感:无论载体是网格、神经组织、振子阵列还是市场,只要能测出这四个量,就能算出同一个数。若这一构造成立,混沌边缘就从「事后的描述」变成「事前的坐标」。

对立预测P1(无量纲不变性)。 取任一可控的伞场系统,在保持N_SDE数值不变的前提下,单独改变r、ω、v中的两个(例如把r加倍同时把v加倍)。本书预测:三态之间的转变位置(以功率谱形态与最大Lyapunov指数变号为判据)在N_SDE坐标上不移动,移动量小于标定误差。混沌边缘传统对此没有任何预测,因为它没有跨实现的参数。

判伪条件。 若在至少三组独立的参数组合下,转变位置在N_SDE坐标上出现系统性移动(超过标定误差的两倍),则N_SDE不是控制参数,第六章的无量纲化构造失败。这是死法F2,见附录F。

42.3 交手二:伞场动力学与非局域耦合振子

最近邻。 第十一章给出的伞场动力学方程,形式为「自身动力学+核加权的差异耦合+噪声」。这一形式在非线性动力学中有一支成熟且庞大的文献。Kuramoto与Battogtokh在2002年研究非局域耦合的相位振子环,发现同一组全同振子会自发分裂为相干区与非相干区共存的状态;Abrams与Strogatz在2004年命名了这一状态并给出精确解,称之为chimera态;Sethia、Sen与Atay在2008年进一步把时滞与非局域耦合合在一起,得到成簇的chimera。这一支的标准结论包括:非局域耦合(既非近邻也非全局)是chimera出现的关键;耦合核的宽度是核心参数;时滞与固有周期之比是标准分岔参数;Hopf分岔进入极限环、倍周期级联进入混沌是标准路径。

重合部分(必须承认)。 本书第十一章的方程是这一族的一个成员,不是一个新族。命题11.1(共振放大)的线性化—特征方程—同步/反同步模式分析是这一支的常规操作;「超临界Hopf分岔后的极限环即介生态」的刻画使用的是标准动力系统语言,本书的推荐语正确地指出了这一点。伞的来回时延t_d = 2r/v与固有频率的匹配条件,正是时滞耦合文献中的时滞—周期比。凡第一版中暗示这套动力学为本书自建的措辞,一律作废。

剩余(本书多说的那一句)。 在chimera文献中,耦合核的宽度是一个由实验者拧的旋钮:你设定它,然后观察系统。本书§31.0的算法要求的是另一件事——核宽度由场自身决定:r = r(Φ),当局部梯度增大时伞自动撑开。这是一个自指参数,它把「谁决定非局域尺度」从外部实验者手里交还给系统。据我们检索,非局域耦合振子文献中尚未系统研究核宽度自指的情形。

对立预测P2(自指核的边界自锁)。 比较两个系统:(甲)固定核宽度的非局域耦合振子环;(乙)核宽度随局部梯度自适应的同族系统,其余参数相同。本书预测:在乙中,相干区与非相干区的边界位置在长时间尺度上收敛到一个由参数决定的稳定值,其在不同初值下的方差显著小于甲;而在甲中边界位置随初值漂移,方差不随时间收敛。直观说法是:自指核会把系统锁在临界宽度上,固定核不会。

判伪条件。 若乙的边界位置方差与甲无统计显著差异(在至少50组随机初值、三组参数下),则自指参数的必要性论证失败,§31.0算法中r = r(Φ)的设计只是工程技巧,不是原理。

42.4 交手三:不合法中间态与模拟退火

最近邻。 「必须接受暂时更差的状态才能到达更好的终态」这一策略,在组合优化中不是新见解,而是1953年Metropolis接受准则与1983年Kirkpatrick、Gelatt、Vecchi模拟退火的核心内容;Wales与Doye的basin hopping、以及各种非单调线搜索方法都建立在同一直觉上。第一版把模拟退火仅作为「只接受合法中间态」的对照组处理,这个对照是不公平的:退火的全部要点正是接受劣解。这一处理方式作废。

重合部分(必须承认)。 全书第三条金句「你必须经过不合法的中间态,才能到达合法路径永远到不了的终态」,若理解为「必须接受质量更差的中间态」,则它是1983年的工程常识,本书对此没有贡献。

剩余(本书多说的那一句)。 但这两件事在数学上不是一件事,而这个区分正是本书第一编的实质。模拟退火翻越的是能量势垒:目标函数在构型空间上有高有低,温度让你以一定概率上坡,从而离开局部极小。它的前提是路径存在——你要去的地方与你所在的地方在同一个连通分支里,只是中间隔着一座山。Santos在2000年证明的是另一回事:三维点集的翻转图可以不连通。不连通意味着路径根本不存在。任何温度调度都无法把你送到另一个连通分支,因为退火只在给定的状态图上做随机游走,它不新增边。要新增边,必须离开合法流形本身——插入临时Steiner点、允许退化单元、即本书所说的隧穿。

一句话的区分:退火解决的是势垒,隧穿解决的是不连通。前者是度量问题,后者是拓扑问题。温度对拓扑无能为力。

对立预测P3(判决性实验)。 在Santos的P₂₄构造上:(甲)任意退火调度的合法flip搜索,无论初温、降温速率与运行时长;(乙)允许经过退化中间态(负体积或零体积单元)后回到合法配置的搜索。本书预测:甲永远到不了T₂所在的连通分支,乙能到;且这不是概率差异而是可达性差异——甲的到达概率恒为零,而不是很小。

判伪条件。 若存在一个纯合法flip的退火调度以非零频率从T₁到达T₂,则要么Santos定理的适用范围被误读,要么本书对「合法flip路径」的界定与定理不符——无论哪种,第一编的必要性论证作废。这是死法F1,见附录F。本书第三十一章实验一就是这个判决性实验的实施方案。

附带的诚实说明: P3的甲方结论其实是Santos定理的直接推论,因此它「几乎必然成立」。这意味着P3不是一个高风险预测,它的价值在于把两个被混为一谈的机制切开,而不在于冒险。真正冒险的是F7(工程增量),见附录F。

42.5 交手四:几何的视觉SIO偏向与Poincaré

最近邻。 §39.4称「几何是视觉SIO的数学」为「本章最锋利的一刀」,并称这是「被忽视了两千多年的事实」。这个说法必须撤回。Poincaré在1902年的《科学与假设》第四章中已经明确写下:表象空间以视觉的、触觉的、运动的三种形态出现,三者都本质地区别于几何空间;它既非均匀也非各向同性,甚至不能说它是三维的;并且他专门指出,触觉空间比视觉空间更复杂,离几何空间更远。他还给出了几何空间由位移群构造出来的路线。Nicod在1930年前后继续了这条线索,Révész在盲人空间知觉的心理学研究中给出了大量触觉空间的经验材料。

重合部分(必须承认)。 「几何不是空间本身的数学,而是某一种感知互动的数学」这一判断,Poincaré在124年前已经作出,且论证比本书第一版更细(他给出了从可补偿的感觉变化到位移群的构造)。§39.4—§39.9的诊断部分属于重述,不属于新辨识。第一版中「被忽视两千多年」「本章最锋利的一刀」两处措辞作废。

剩余(本书多说的那一句)。 Poincaré提出了问题并停在那里:他指出触觉空间不同于几何空间,但没有为触觉空间写下公理,没有给出度量,没有给出曲率。此后一百年,这个空白基本保持。本书§39.12—§39.13的工作是把它形式化:接触事件Tacton取代点作为原子对象(公理T1)、序列性取代同时性(T2)、路径依赖导致度量非对称(T3)、有效维度随压力涌现(T4),并据此定义触觉能量泛函、触觉测地线与触觉曲率κ_T = dΦ_press/ds。其中T3给出的非对称性把触觉空间推出黎曼框架、推入芬斯勒框架——这是一个有具体数学后果的断言,Poincaré没有做过。

对立预测P4(触觉度量的非对称性)。 取一段具有明确各向异性纹理的表面(如定向绒面或锯齿纹),令被试以恒定压力与速度沿同一路径分别在两个方向滑动,报告主观「触觉路程」。本书预测:两方向的判断出现系统性差异(顺纹判短、逆纹判长),差异幅度随表面各向异性程度单调增大;在各向同性表面上差异消失。

判伪条件。 若在至少两种各向异性表面、n≥30的被试上,两方向的触觉路程判断无系统差异(差异的置信区间包含零),则公理T3假,触觉几何退回黎曼框架,§39.13的芬斯勒论证与命题39.1作废。这是死法F6,见附录F。

42.6 交手五:N_SDE与无量纲群家族

最近邻。 §6.1把N_SDE类比于Reynolds数。这个类比选得不够近。更近的邻居有三个:Deborah数——材料的松弛时间与观测时间之比,用于判断同一材料表现为固体还是流体,Reiner在1964年提出;Knudsen数——分子平均自由程与系统特征长度之比,用于判断连续介质假设何时失效;Damköhler数——反应时间与输运时间之比。

重合部分(必须承认)。 把N_SDE按定义拆开:r/v是伞内单程传播时间,因此r/(v·T_c)正是一个传播时间与特征时间之比,即一个Deborah型的数;而r相对于域尺度L的比值r/L是一个Knudsen型的数——它正是「连续/局域假设何时失效」的经典判据。因此N_SDE = ω·[r/(v·T_c)]的正确定位是:一个Deborah型时间比,被重叠度ω加权。这个定位比Reynolds类比准确得多,因为Reynolds数是惯性与粘性两种力之比,与本书的构造不同源。第一版的Reynolds类比降格为修辞性说明,不再作为构造依据。

剩余(本书多说的那一句)。 Deborah数刻画时间尺度之比,Knudsen数刻画长度尺度之比,两者都是单一比值。N_SDE的构造是把长度非局域性(r)、多重覆盖(ω)与时间延迟(r/v相对T_c)合并为一个数,并主张三者对系统行为的影响只通过这一个组合起作用——这就是P1所检验的不变性主张。据检索,尚无文献提出把重叠度作为无量纲群的一个因子。

对立预测P5(ω的乘性)。 本书预测:ω在N_SDE中以一次幂进入,即在其余量固定时,把重叠度加倍与把传播时间加倍对三态边界位置的影响相同。这是一个可以直接被扫参数否定的强断言。

判伪条件。 若实验显示ω的影响呈非线性(例如需要ω^α,α显著偏离1),则N_SDE的具体函数形式错误——注意这不否定「存在一个控制参数」的主张,只否定当前形式。修订应给出α的经验值并说明其来源。

42.7 作废清单与尚未交手的名单

基于本章,第一版以下措辞正式作废,修订版已删除或改写:

(一)任何暗示介生态为本书首次辨识的表述——它是混沌边缘现象的重新命名,本书的贡献限于给出跨实现的控制参数。

(二)任何暗示伞场动力学为本书自建方程族的表述——它是非局域耦合振子族的一个成员,本书的贡献限于自指核。

(三)把模拟退火仅作为「只接受合法中间态」的对照——退火本就接受劣解,正确的区分是势垒与不连通。

(四)§39.4「被忽视了两千多年」「本章最锋利的一刀」——Poincaré 1902已作出该判断,本书的贡献限于形式化。

(五)§6.1「类比Reynolds数」作为构造依据——正确的最近邻是Deborah数。

(六)§6.7「N_SDE≈0.14恰好落在秩序—介生边界上,这不可能是巧合」——在三态阈值尚未标定之前,这句话没有内容,见§6.10的修订。

以下是本书尚未与之正面交手、留给下一版的名单,列出以免再次漏账:单位分解法与广义有限元的enrichment理论(Babuška-Melenk 1997、Belytschko 1994、Liu的再生核质点法1995)与伞场的关系;图极限(graphon)作为稠密耦合矩阵的连续极限;分数阶导数的记忆核与时间方向上的伞;几何测度论在伞场边界刻画中的位置;以及科学哲学方面,Cartwright与Hacking关于模型与现象学定律的讨论对本书「纲领」自我定位的约束。

42.8 本章的自我适用

本章要求每一个构件说出自己的最近邻与死法。这个要求同样适用于本章本身:若某位读者指出,上述五条对立预测中的任何一条其实已在某篇文献中被提出并检验,则该条的「剩余」资格立即取消,本书应在下一版中删除相应的优先权主张并致谢指出者。作废是常态,不是失败;未交手才是失败。

本页节选自《SDE数学导论(修订版)——结构·差异·纠缠:数学的重新定义》,王德生著,Demai International Press。