我的学习群里全是真大佬 第452节

　　【看完后告诉我，能不能把刚才的事件分成两堆：一堆幻峰，一堆真信号。】

　　【嗯！小黑听主人的。】

　　压缩包刚喂过去的瞬间，小黑没有卖萌打嗝。

　　它在原地缩了一下，然后头上的进度条开始疯长。

　　【0.013】

　　【0.027】

　　【0.063】

　　【0.094】

　　直到【0.1】，它才堪堪停住。

　　【主人，这份东西很简单嘛，小黑可以做到的呀。】

　　李东咽了一口口水，立刻下达指令。

　　【把你刚才吃下去的几万条事件，按这个规矩全部过一遍，挑出真信号，按时间排好！】

　　【小黑学习中……】

　　仅仅过了三秒钟。

　　小黑头上就出现了一行字。

　　【主人，挑好啦，一共三万八千四百二十七条事件。】

　　【真信号：九条，其余都是幻峰。】

　　李东微微张着嘴，发出阿巴阿巴的声音。

　　三秒，我勒个亲娘呢！

　　这份活儿如果丢在旁边这台笔记本上，满载跑也得大半个晚上，即便明天上了金陵大学的超算，估计也不会比小黑快。

　　而这个黑黢黢的小东西……用了不到三秒。

第304章 顺手的事

　　第二天清晨。

　　李东合上电脑时，窗外的天色已经微微亮了。

　　草稿纸散落一地，他揉了揉有些发胀打的太阳穴，然后看向了屏幕右下角。

　　小黑头上的进度条已经来到了【0.18】。

　　【主人，一百九十六万条，全吃完啦。嗝……】

　　李东失笑，敲下键盘。

　　【真信号多少？】

　　【八万六千三百九十二条，其余全是幻峰。】

　　李东点点头。

　　【小黑，挂载昨天的脚本，叠第二层筛子。】

　　【在十几个TeV能段，找形状符合AD 774/5外推预测的事件，χ?/d.o.f.卡在1.1以内。】

　　小黑没有废话，又开始蹦蹦跳跳起来。

　　到二十秒，一行字弹了出来。

　　【主人，挑好啦，一共十一笔。】

　　李东嘴角抽了一下。

　　十一笔。

　　跟昨晚他用那三十七笔“待复核”事件包外推出来的数字，分毫不差！

　　如果说昨晚是还是弱证据的话，那现在，用悟空号八年全量数据跑出来的结果，就是强证据了。

　　深吸一口气，李东知道，这第一关算是过了。

　　接下来才是真正的难搞的事。

　　他要确认十八个月、卡灵顿级的窗口期到底是不是正确的，甚至他还想将三宅级的时间也给推出来。

　　他抽出桌底那沓空白草稿纸。

　　李东压根没去等金陵大学超算的批文。

　　他要走的这条路，跟算力没什么关系。

　　十一笔真信号的时间节点被他重新铺开。

　　但是样本其实还是太少了点，传统的傅里叶变换拉不开频率分辨率。

　　于是他在纸上飞快写下了一行字。

　　【Lomb-Scargle Periodogram】

　　Lomb-Scargle周期图。

　　这是天文学界专门拿来对付“不均匀采样时间序列”的老把式，五十年前由两位天文学家先后独立提出。

　　它不要求数据点之间间隔均匀，而是在频率空间里直接做最小二乘拟合。

　　低样本量、不均匀间隔、信号叠加在噪声里，这三件事，是天文学家最难受的事。

　　而Lomb-Scargle就是为这三件事量身定制的。

　　李东拿起笔，开始倒推频率。

　　T?≈ 73.4天。

　　也就是说，如果这十一笔事件背后藏着一只钟，那这只钟，每七十三天敲一下。

　　对应的就是太阳巨大的磁场势阱底下，单位钟摆周期内，外面的高能粒子撞击势垒一次的频率。

　　它没法直接告诉你主峰何时爆发。

　　但只要再叠上一层东西，主峰自身的“激发概率”就能反推出来。

　　这层东西就是……

　　伽莫夫的量子隧穿衰变框架。

　　不过需要改上一点……

　　核内势阱改“太阳磁场势阱”、α粒子改“高能耀斑外溢通量”……李东在草稿纸上魔改着微分方程。

　　十一笔事件的振幅递增曲线被套入方程，MCMC（马尔可夫链蒙特卡洛）算法瞬间跑完五万步采样。

　　后验分布的均值，最后钉死在了一个坐标上。

　　【主峰抵达窗口期：T?+ 18.3个月（± 2.4个月）】

　　【量级：卡灵顿级】

　　【后验置信度：78.6%】

　　和昨晚的结果几乎一样，但这一次，它的可信度更高了。

　　就在这时，桌上的手机突然震动了起来，是个陌生号码。

　　“李东老师吗？我是金陵大学超算中心的，您申请的大规模回测节点批下来了，A100满规格，请问您什么时候来开账号？”

　　李东愣了一下，看了一眼屏幕右下角打饱嗝的小黑。

　　“不用了，老师。”李东笑了笑。

　　“计划有变，这节点让给别人排吧，我已经跑完了。”

　　电话那头一时间安静的有些诡异。

　　两百万条全量数据，你一天不到跑完了？拿什么跑的？算盘吗？

　　没等对方回过神，李东直接挂断了电话。

　　草稿纸开始散落在房间的各个地方，书桌、地板……

　　李东要把±2.4个月的窗口期继续往里压！

　　他开始引入羊八井观测站三年来的“待复核”异常邮件。

　　三笔同向事件叠入MCMC！

　　置信度飙升至83.4%！窗口期压至±1.9个月！

　　然后再引入LHAASO（拉索）观测站的TeV级宇宙线月报。

　　置信度推到86.7%！窗口期死锁在±1.5个月！

　　李东所有基础属性直接拉满，彻底进入了那种“洞穿真理”的专注状态。

　　天黑天亮，门外的风吹草动，全被他屏蔽了。

　　他的世界里只剩下那只在地球上空敲了八年的“钟”。

　　他正一笔一笔地，把死神敲门的倒计时，从宇宙的噪声里给硬生生抠出来！

　　……

　　两天后。

　　金陵某高层公寓。

　　“咔哒”一声，备用钥匙拧开了防盗门。

　　沈澈和张文平眉头紧锁地走了进来。

　　屋子里的气味有些发霉。

　　书桌前，一个人影正埋在几乎要将他淹没的草稿纸堆里。

　　“李东？”

　　李东没有回应。