运行。
笔记本的风扇开始加速,屏幕上不再是简单的曲线,而是一组复杂的、由多个对称“山峰”构成的“角图”(cornerplot)。
各参数的一维边缘分布多为近似对称的单峰;
两两关联的子图则呈现出明显的倾斜椭圆,显示出在强噪声与错核拉长下的参数相关性。
“我靠,寧神,”
宋子阳的脑袋凑了过来,表情像是看见了外星生物,“你这……电脑屏幕上又是什么神图?”
林允寧没理他,刚摘下耳机,准备將结果截图,一个清脆又带著一丝不確定的声音从身旁传来。
“林允寧……你这个是角图吧?这是……在用蒙特卡洛方法吗?”
他转过头,发现秦雅不知何时已经站在了他的座位旁。
她一只手拿著本物理习题册,另一只手指著屏幕,语气有些不確定,但眼睛里闪烁著求知的光芒,“我之前看一篇关於分子动力学模擬的论文,里面提到过用这种隨机採样的方法,来处理高维参数空间。你……你这是在用马尔可夫链蒙特卡洛吗?”
她补充道,指尖下意识地攥紧了校服的衣角。
林允寧有些意外地看了她一眼,没想到她不仅认识角图,还能准確地说出mcmc方法。
这小丫头,知识面可以啊。
“对,是马尔可夫链蒙特卡洛。”
他言简意賅地解释道,“我这里有个反演问题,没法直接求出精確解,所以用它来框定参数的后验概率范围。”
秦雅若有所思地点了点头,没有再追问细节,只是將“后验概率”这几个字牢牢记在心里,然后才想起了本来的目的,举起手里的习题册:
“这道关於电容器的题……我有点不明白。”
林允寧瞥了一眼,是道有点复杂的电磁学题。
他隨手拿过一支笔,在草稿纸上画了个草图,三两句就点明了关键。
最前排,赵宇轩的眼角余光偷偷將这一切都瞧在心里,咬了咬牙,將头埋得更深。
联考在即,分数才是一切。
搞这些里胡哨的东西,不过是譁眾取宠。
在秦雅道谢离开后,林允寧將结果截图,附上一段简短的说明,通过qq发给了韩至渊。
这一次,qq却弹出了自动回覆:
【正在开会,有事请留言,稍后联繫。】
林允寧鬆了口气,正准备起身活动活动,qq上陈正平的头像却闪动起来。
【陈正平】:林师弟,韩老师刚才给我看了你的结果,第一关过得漂亮!他说如果你觉得第二关也完成了,让我將最后一个考验的数据发给你。他开会回来再和你聊。
林允寧笑了笑,打字回覆:【多谢师兄,第二关的代码我已经发给韩老师了,你能把后面的数据传过来吗?】
【陈正平】:(惊讶表情)这么快?第二组是『脏数据,信噪比特別低,irf也不对称,你怎么改进算法的?
【林允寧】:我用了mcmc和贝叶斯框架,我想著与其去找那个不可能存在的精確解,不如用蒙特卡洛模擬搜索解的概率分布。
【陈正平】:mcmc?我靠,这个思路……绝了!放弃精確解,转而去勾勒整个后验概率分布……林师弟,你的思维太灵活了!我这就把最后一组数据发给你。
很快,两个文件stress_test_data。dat和instrument_response_b。dat便传了过来。