天才一秒记住【热天中文网】地址：https://www.rtzw.net
他看向孙婧，露出了一个迷死人不偿命的笑容，“我们不去消除误差，我们理解误差。
“第三张王牌，就是正向拟合。
“我们用求解到的谱函数，去和我们已知的IRF做卷积，然后用这个卷积后的理论结果，去拟合我们测到的、最原始的、失真的实验数据。
“这样，我们拟合的，不再是信号，而是‘信号通过这台有缺陷的仪器后，应该是什么样子’。”
他摊了摊手，一脸的理所当然。
这番话，如同平地惊雷，让孙婧彻底呆住了。
她感觉自己多年建立起来的实验物理三观，正在被对方这种天马行空的逻辑剧烈地冲击着。
在她身边，陈正平的眼睛却猛地亮了起来，呼吸都变得有些急促。
他死死盯着白板，只见上面没有一个长公式，只有三条逻辑闭环的箭头。
如同三枚精准的图钉，将那个原本飘忽不定的反演问题，死死地钉在了一个极小的的解空间里！
他终于明白了林允宁想做什么！
“你的意思是……用一套带物理约束的迭代算法，直接从我们测到的，包含了所有仪器误差和噪声的光谱里，把最底层的‘声子自能’给硬生生构造出来？！”
“不是‘构造’，”
林允宁放下笔，纠正道，“自能函数作为物理实在，它本来就在那儿，不会因为我们算不出来，它就不存在。
“所以，问题不是‘构造’，而是用‘反演’来‘观测’??用一种更聪明的办法，透过噪声的迷雾，直接看到它。”
他的目光，最后落在韩至渊身上，进行了最后的总结。
“想要得到同时描述‘线宽’和‘频移’的自能函数，常规方法是正向推导，从原因到结果。
现在咱们反过来，给他来个‘定果求因’??直接把声子自能Σ这个‘黑箱’的盖子掀开，把它抓在手上。
“到那时候，我们就可以用这个统一的理论框架，来讨论它里面到底藏着非谐效应、缺陷散射还是电声子耦合。
那才是在用同一种语言说话。”
话音落下，他打了个哈欠，长长吐出一口浊气。
好像说了这么长的一段话，让他有些精疲力尽似的。
屋里安静了几秒，只有窗外传来几声模糊的鸟鸣。
“物理图像……很完美，不过算法呢？”
这次，轮到了韩至渊发问。
他一针见血地指出了从图像到现实的最后一道障碍，“想要解决这个问题，要有合适的算法和工具。
“但据我所知，现有的软件包，都解不了你这个带复杂物理约束的反演问题。
这个世界上，还没有人写过这样的‘反演翻译器’……”
“那就从我们这里开始写！”
林允宁的回答云淡风轻，仿佛在说一件微不足道的小事，“求解器用用投影梯度算法（POCS），保证每次迭代都在‘因果一致’这个圈里，外层套一个带界约束的拟牛顿法（L-BFGS-B）处理非负约束。
“至于Σ本身，我们也不需要知道它的精确解析式，更不需要用复杂的函数去猜。
我们可以用非负样条函数或者有理函数来参数化，把求解一个无限维的函数问题，简化成求解几十个参数的优化问题。
这样可以最大限度地减少数值病态。”
他看了一眼已经彻底呆滞的陈正平和孙婧，干脆直接又扔下一颗“重磅炸弹”
：
“我已经用Python写了个原型，在您借我那台X40笔记本电脑上试过了。
用这种方法，即使是老奔腾处理器，也能在半个小时以内收敛。”
“你……你说啥？”
……
本章未完，请点击下一章继续阅读！若浏览器显示没有新章节了，请尝试点击右上角↗️或右下角↘️的菜单，退出阅读模式即可，谢谢！