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

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

To～ 73.4天。

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

对应的就是太阳巨大的磁场势阱底下，单位锺摆周期内，外面的高能粒子撞击势垒一次的频率。它没法直接告诉你主峰何时爆发。

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

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

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

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

【量级：卡灵顿级】

【後验置信度：78.6%】

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

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

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

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

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

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

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

三笔同向事件叠入MCMC！

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

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

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

（本章未完，请点击下一页继续阅读）

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

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

李东失笑，敲下键盘。

十一笔。

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

如果说昨晚是还是弱证据的话，那现在，用悟空号八年全量数据跑出来的结果，就是强证据了。深吸一口气，李东知道，这第一关算是过了。

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

这层东西就是……

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

不过需要改上一点……

核内势阱改“太阳磁场势阱”、α粒子改“高能耀斑外溢通量”……李东在草稿纸上魔改着微分方程。十一笔事件的振幅递增曲线被套入方程，MCMC（马尔可夫链蒙特卡洛）算法瞬间跑完五万步采样。後验分布的均值，最後钉死在了一个坐标上。

【真信号多少？】

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

李东点点头。

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

【在十几个TeV能段，找形状符合AD 774/5外推预测的事件，x/.f.卡在1.1以内。】小黑没有废话，又开始蹦蹦跳跳起来。

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

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

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

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

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

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

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

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

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

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

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

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

【Lomb-ScarglePer iodogram】

Lomb-Scargle周期图。

第304章 顺手的事 (第1/3页)

第二天清晨。

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

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

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

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

李东嘴角抽了一下。

阅读我的学习群里全是真大佬最新章节 请关注舞文小说网(www.wushuxs.net)