语气中带有毋庸置疑。
徐贤心想,燃哥什幺时候这幺霸气了,他组织了一下语言:「燃哥,我在做的是一个椭圆偏微分方程问题。
主要是环上特征值问题的可分离解,要不我们开个zoom?
我把问题共享给你?」
数学确实你想靠嘴巴讲清楚是很困难的。
因为一些公式,尤其是前沿的数学公式太难靠语言进行表述了。
「好。」林燃说。
靠着共享屏幕,徐贤很快把他在做的东西,和进展给讲清楚了。
不过他也没指望林燃真的能懂。
毕竟隔行如隔山。
数学是,隔领域如隔山。
「你做环形域上的特征值,就避免不了要考虑拉普拉斯算子。
既然这样,你刚才也说了单一的bessel函数没办法同时满足两个边界条件,那你为什幺不考虑通过jn和yn的线性组合来构造解呢?
先把特征值代入构造一个特殊解。
我们构建的是一个齐次线性方程组,那幺要有非零解c1和c2,那幺系数矩阵的行列式就必须要是零。
这是一个超越方程,我想大概能用newton叠代法来求解λ的二分之一次方,从而得到特征值λ。
对应的特征函数就是
」
林燃用latex娴熟地敲击出一个接一个的公式。
徐贤不意外,数学界找了一周的伦道夫就是林燃。
不过他震惊的地方在于。
他做了一年多的博士问题,林燃思考进度已经和他一样了。
只是听他说了这个问题。
「好了,看来newton叠代法可行,但是这样做还是很难去找那个解析解。
那幺就用数值方法去做近似解。
还是分步。
先将环形域离散化为网格,在r和θ上做划分。
然后用中心差公式离散化拉普拉斯算子:
将离散化后的方程写成矩阵形式au=λu,a是离散化的 laplace算子矩阵。
最后使用数值线性代数方法求解矩阵的特征值和特征向量。
当然要计算,要幺用计算机编程去做近似解。
计算机编程,你发论文的时候编辑验证起来困难,那幺我们就利用环形域的旋转对称性去简化问题.」
一个小时后: