U盘插进接口,笔记本风扇随之嗡嗡作响。闷头转了好几秒,屏幕上才迟迟跳出文件目录。整整四十个G,两百多个子文件夹,按患者编号依次排开。林允宁没急着去找孟兰的数据,而是先理了一遍目录结构,以摸清程新竹的分类习惯。每个患者对应一个主文件夹,点进去按采集日期分层,最底端躺着逐通道的原始时间序列文件。标准的欧洲数据格式(.edf)。他顺着列表往下划,点开了编号“Ad-02-mXL-001”的文件夹。mXL正是孟筱兰拼音的首字母缩写。里面并排躺着十七个按日期命名的子文件夹,最新的一份生成于上周。程新竹做事很细,每个目录下都额外塞了个纯文本标注,用来记录当天采集时的特殊事件。林允宁从中过滤出含有“幽灵吸引子”字样的三次记录,挑了最近的一次点开。屏幕上刷新出六十四个通道的原始时间序列,采样率为一千赫兹。他将标注里名为“高相干窗口”的时间段单独截取出来。从第1247秒到第1264秒,一共十七秒的切片。一千赫兹的采样率乘上六十四个通道,意味着这短短十七秒内挤着一万七千个时间点,而每个点上又叠加着六十四个维度的数值。林允宁抽过一张草稿纸,随手写下泛函C[∮]的定义式,笔尖悬在了半空。怎么把这堆庞杂的数据塞进公式里?C[∮]的逻辑其实很清晰:输入场构型∮,输出一个标量,用以衡量该构型在局部区域内维持凝聚态的能力。难点在于,这里的∮到底该指代什么?在NS方程里它是速度场;在杨-米尔斯理论中它是规范场的联络形式。无论哪种,都是定格在某一时刻的静态对象,据此,C[p]才能给出“当前构型下凝聚态是否稳定”的判定。然而,杂乱的活体脑电数据根本给不出如此干净的切片。六十四个电极通道,勉强能类比成物理空间中的离散点。若把每个时间点上这六十四个读数揉成一个高维向量,硬生生套作某时刻的“场构型”∮去跑C[∮]的计算,倒也能吐出一个结果。林允宁顺着这个思路试了一把。他用第1247秒的读数算了一次,得到的C[p]值是0.73。用第1248秒的读数算了一次,得到的是1.42。两个相邻时间点,一毫秒之隔,C[p]的值翻了将近一倍。他又接连算了几个点:0.31,1.87,0.58.这五个数字落在草稿纸上,毫无规律地剧烈震荡着。假设临界值为1.0,那按照这组数据的逻辑,孟筱兰的大脑简直是以毫秒为单位,在“凝聚”与“崩溃”的边缘疯狂闪烁。这根本说不通。幽灵吸引子在临床上会展现出十五到二十秒的高相干窗口期,期间患者的认知能力会得到显著提升,直至窗口期猝然崩塌。倘若微观上的凝聚态真的一毫秒翻转一次,宏观层面绝不可能支撑起长达十七秒的连续相干行为。出现问题的原因也很简单。脑电信号的毫秒级波动里裹挟了太多高频噪声,单点读数根本没资格充当所谓的“场构型”。必须引入时间窗口做平均处理。他先切了个一百毫秒的滑动窗口重新跑了一遍。经过平滑,C[∮]的曲线温顺了不少,在整个高相干期内稳稳停留在1.2到1.5之间,似乎确立了临界值之上的优势。可当他把窗口缩窄到五十毫秒时,均值瞬间跌到了0.8附近,大半截身子淹没在及格线之下;要是把窗口拉长到二百毫秒,均值又一路高歌猛进到了1.9。他不信邪,又凑了二十、一百五、三百毫秒几个梯度挨个试水。结果,六种不同的窗口宽度,硬是跑出了六条南辕北辙的曲线。面对同样的数据和泛函,一百毫秒的尺度信誓旦旦地昭示“凝聚态稳定”,五十毫秒的尺度矢口否认它的存在,而到了二百毫秒那里,凝聚态又变得牢不可破。仅仅调整了平滑参数,得出的结论便如同儿戏般推倒重来。显然,这种高度依赖人为设定尺度所榨出的结果,毫无说服力。林允宁将这六条杂乱的曲线叠画在同一张纸上,蹙眉端详。渐渐地,他察觉到了一丝异样。尽管这些曲线在绝对数值上各唱各的调,但它们共享着一个极其隐蔽的特征——无论窗口设定得多宽多窄,在第1261秒至1262秒的区间内,所有曲线都遭遇了一次断崖式的暴跌。暴跌前后的走势或许千奇百怪,但这一处塌方却如铁律般横跨了所有的尺度。而且坡度极其陡峭。在任何一种窗口下,C[p]的值几乎都在一两个时间步内,从高位瞬间被清零。这是一次跃变。整整十七秒的高相干窗口,其崩解过程居然全部压缩在极短的瞬间完成,连半点缓冲的余地都没留。林允宁放下笔,回去看原始数据。他把第1261秒到第1262秒之间的毫秒级原始时间序列单独拉出来,不做任何平均和滤波,直接看六十四个通道的相干性指标。结果一目了然。在第秒之前,多通道的相位同步指数稳稳咬在0.7以上,维持着典型的高相干态;但在跃过秒的界限后,仅仅三个毫秒,该指数就从0.72垂直砸到了0.09。短短三毫秒,整个系统从高度协同滑向彻底失联,仿佛被人一刀切断了引线。这种形态太眼熟了。它像极了NS方程中的涡量爆破,也与杨-米尔斯场在特定构型下的相变如出一辙。事实上,当初设计C[p]这个泛函,正是为了去捕捉这种极端的物理行为。但在这里,C[p]却彻底宕机了。根本原因在
本章未完，请点击下一页继续阅读》》