91书院(91shuyuan.com)更新快,无弹窗!
徐贤心想,燃哥什麽时候这麽霸气了,他组织了一下语言:「燃哥,我在做的是一个椭圆偏微分方程问题。
主要是环上特徵值问题的可分离解,要不我们开个zoom?
我把问题共享给你?」
数学确实你想靠嘴巴讲清楚是很困难的。
因为一些公式,尤其是前沿的数学公式太难靠语言进行表述了。
「好。」林燃说。
靠着共享屏幕,徐贤很快把他在做的东西,和进展给讲清楚了。
不过他也没指望林燃真的能懂。
毕竟隔行如隔山。
数学是,隔领域如隔山。
「你做环形域上的特徵值,就避免不了要考虑拉普拉斯算子。
既然这样,你刚才也说了单一的Bessel函数没办法同时满足两个边界条件,那你为什麽不考虑通过Jn和Yn的线性组合来构造解呢?
先把特徵值代入构造一个特殊解。
我们构建的是一个齐次线性方程组,那麽要有非零解c1和c2,那麽系数矩阵的行列式就必须要是零。
这是一个超越方程,我想大概能用NewTon叠代法来求解λ的二分之一次方,从而得到特徵值λ。
对应的特徵函数就是
」
林燃用Latex娴熟地敲击出一个接一个的公式。
徐贤不意外,数学界找了一周的伦道夫就是林燃。
不过他震惊的地方在于。
他做了一年多的博士问题,林燃思考进度已经和他一样了。
只是听他说了这个问题。
「好了,看来Newton叠代法可行,但是这样做还是很难去找那个解析解。
那麽就用数值方法去做近似解。
还是分步。
先将环形域离散化为网格,在r和θ上做划分。
然后用中心差公式离散化拉普拉斯算子:
将离散化后的方程写成矩阵形式Au=λu,A是离散化的Laplace算子矩阵。
最后使用数值线性代数方法求解矩阵的特徵值和特徵向量。
当然要计算,要麽用计算机编程去做近似解。
计算机编程,你发论文的时候编辑验证起来困难,那麽我们就利用环形域的旋转对称性去简化问题.」
一个小时后:
「总之环形域上的特徵值问题由于边界条件的复杂性,解析解难以直接获得。
使用Bessel函数的线性组合并结合数值方法求解超越方程是一种可行的解析-数值混合策略。
而我们再结合了有限差分法,这样就提供了通用的数值解法。
后续你还可以根据具体需求,例如精度丶计算资源或理论洞察,选择适合的方法进一步探索。」
徐贤是真麻了。
人已经彻底麻了。
属于是那种,不知道自己是谁,自己在什麽地方,自己要干什麽的麻。
从来没有如此麻过。
「燃哥,我们之间已经隔了一层可悲的厚壁障。」林燃最后的总结说完后,徐贤说道。
他旁边床位的室友扭头看了一眼,觉得徐贤真是莫名其妙。
「怎麽?你闰土了是吧。」林燃一下就知道徐贤在玩什麽梗。
徐贤这才想起来寝室里还有室友在呢。
当下一些高校没开学,开学了的高校也号召大家别离开校园。
燕大好点,但好的不多,大部分人都挤在图书馆。
如果你早上没能占到座位,那就只能在寝室了。
今天他们寝室两个人。
很不幸,都没有抢到座位,只能呆在寝室。
徐贤低声道:「不是,燃哥,你这也太夸张了。
都说富人靠科技,穷人靠变异。
尽管我印象里一直判断不了你是富还是穷,但你这到底是用了科技还是用了变异?
你知道这个问题,是我导师给我留的博士五年做出来的大问题。
你一个半小时就咔咔咔给我算完了,甚至这过程中要用到的数值计算都是直接手撕。
这未免有点太变态了吧?
虽然我现在还判断不了你说的对不对,但这也已经很变态了。
至少我这个问题想了半年多,挑不出你说的方法有什麽毛病。
这个问题算是我导珍藏的,说是给年轻博士打怪升级的问题,说我五年能做出来就不错了,可以在燕大顺利毕业。
你之前应该不知道才对。
如果你做的是对的,相当于我博士五年只要再做一个小问题就能从我导那毕业了?
相当于三年不到就毕业了,不是燃哥,我前脚才说带我飞,后脚你就直接帮我跳过练气直接到金丹期了,你到底是变异还是科技啊?」
一般来说,你去念数学博士,好的导师会给你两个问题,一个大问题,