第343节
第343节
他清晨想的那一套物理切口,必须要有数据支撑。 不是说他物理那一头脑补得再漂亮,就一定是对的。 他得拿数学这边能跑出多远来做参照。 万一数学方向真的只是小卡一下,后面一路畅通呢? 那他这一晚上的物理切口,就成了无用功了。 他不能拿一个直觉去推翻整个组的工程进度。 所以他得先去把数学这边推到尽头。 李东把包放在桌上,也不开电脑。 他先把草稿纸和笔拿了出来。 李东是搞纯数的。 反问题,严格的说他不算专业。 可是这些天,鄂伟南院士给的那一份书单,他已经啃下了一大半。 Hadamard的良态性条件、Tikhonov正则化、谱截断、Backus-Gilbert方法…… 这些名字加上对应的核心定理,全都扎扎实实地挂在他记忆宫殿的【应用数学/反演】那一格。 再加上他理解东西的速度,本来就在常人之上。 因此他也有信心不比一般走应用数学的学者差。 他坐下来,开始推导。 第一步:测量算子K:T→S的形式定义,T是配位场张量空间,S是隧穿谱空间。 第二步:核函数展开,把K拆成本征函数序列。 第三步:上吉洪诺夫正则化,调整参数α。 第四步:把K的奇异值谱分布,对应到T的低秩近似。 第五步…… 他一步一步地推。 每一步,他都对照着昨晚已经搭好的那一套物理参数,看数学这一头能不能合上。 基本上前面五六步,都没有出大问题。 Tikhonov那一段他卡了不到十分钟,就用一个加权范数把它处理掉了。 基函数的选择那一段,他换了一组Whittaker型的局部基底,规避了一部分相位混淆。 …… 直到第七步。 李东看着自己写出来那一行,皱起了眉头。 【在(L?(R?,dμ),〈・,・〉_K)上对相位算子做谱分解,要求三组基底的过渡矩阵收敛于一个无相位歧义的酉变换……】 “是这样吗?” “不对呀。” 这一段,本来是要把不同基底的展开“硬合”在一起。 按列正常思路,是用一个隐性约束量去反卡这个过渡矩阵的上界。 可是李东推到这里,就是过不去。 他换了三种正则化,结果一种比一种糟糕。 他试着把鄂伟南院士点出来的那一根“隐性约束”硬塞进去。 卡得住,又卡不死。 第三峰从糊里头能拎出来一点点边,但拎到一半,又掉回去了。 李东盯着那一行字,叹了一口气。 “……果然啊。” “全国十几个组,全都卡在这。” …… 就在他第三次放下笔的时候,背后突然冒出一个声音。 “卡住了?” 李东猛地回过头。 吴开教授不知什么时候已经站在了他的身后。 李东赶紧站起来。 “吴老师。” “您来了。” 吴开摆摆手,走过来在他对面那张椅子上坐下。 “我刚从化院开完会过来。” 他看了一眼李东桌上的草稿纸。 看见最上头那一张推到一半的过渡矩阵谱分解。 吴开点了点头。 这一块,他也是看得懂的。 不是因为他数学多好。 是因为这一块,全圈儿的人都卡在这儿。 卡得太多了,他不熟也熟了。 吴开笑了笑。 “李东,你别急啊。” “这个东西,已经卡好几年了。” “国内国外大家都卡在你这一步呢。” “你今天上午就推到这儿,已经比绝大多数人快了。” “再说,你也别给自己压力。” 他温和地看着李东。 “我也没想过你一来就能解决问题啥的。” “咱们组的进度计划,是按‘数学这一块磨一两年’排的。” “前面那一堆样品和针尖,咱们可以先把后续的衔接做了。” “等你这边反推回来,咱们再把它接上,也是可以的。” 吴开这话说得是真心的。 他怕李东受挫。 李东毕竟还年轻。 基础数学那一头,他跑得是真的太快了。 这种年少成名的孩子,最怕的就是在另一个领域里被一面墙撞得头破血流。 吴开过来的不光是问候。 他是来给李东托个底的。 …… 李东点了点头。 他听出了吴开的好意。 心里挺感激。 他正要开口说“老师我会好好磨的”。 结果话还没出口,他想起来一件事。 他刚才推到第七步卡住的时候,没再继续往第七步那一堵墙上撞…… 他用了另一种“绕”的方法,从纯数那一边借了一招过来,把第七步先放到一边,跳到了第八步。 第八步往下,他一口气推了好几页。 也不知道这样行不行,毕竟大家都没用纯数的方法来,万一在这里不行呢? 所以他问道。 “老师,这后头那一段我用的自己的方法,您给看一下这样行不行。” 李东说完,就低头去翻草稿纸最下面那几张。 然而听到这话的吴开整个人就在旁边僵住了。 ??? 啥意思啊?你推过去了? 吴开数学不算精通。 他这辈子混的是化学合成。 但数学这一头,因为这个项目卡了好几年,他这一两年也开始“陪着读”。 反演这里头,每一步在哪儿、为什么会卡、卡在哪个具体的点上…… 他门儿清。 刚才李东桌面上那一张过渡矩阵谱分解。