adorerdu adorerdu
关注数: 4 粉丝数: 0 发帖数: 17 关注贴吧数: 17
数值求解逆变换和解析解不符 我的问题需要用到mellin逆变换，但没有解析解，所以想用个简单算例比较下数值积分与解析解的误差。代码如下： r = 2; f[s_] := Gamma[s]*Gamma[1 - s] InverseMellinTransform[f[s], s, r] NIntegrate[1/(2 \[Pi] I) f[s] r^-s, {s, 1 - 5 I, 1 + 5 I} (*Exclusions->{0}*), WorkingPrecision -> 20, MaxRecursion -> 20] // Timing NIntegrate[ Re[1/(2 \[Pi]) f[1 + p I]/ r (Cos[p Log[r]] - I Sin[p Log[r]])], {p, -5, 5} (*Exclusions->{0}*), WorkingPrecision -> 20, AccuracyGoal -> 8, MaxRecursion -> 100, PrecisionGoal -> 100, Method -> {"GlobalAdaptive", "MaxPoints" -> 1000}] // Timing 第一个输出的是解析解的结果，第二个是直接按照mellin逆变换的积分定义写的数值积分，第三个是按照定义写出积分后令积分变量s=1+i*p（以使得积分上下限为实数）。 我画出了这个函数实部的图像，虚部的图像，是关于y轴对称且数值集中在[-3,3]，所以数值积分区间取了[-5,5]。但是第一个数值解给出的结果是复数，且实部也不等于解析解结果。第二钟数值解因为虚部积分是0，为了避免奇点，我直接在实部积分，但是结果与前一种一样，依旧对不上解析解。 想知道有没有大佬做过类似的逆变换问题，有什么解决的思路吗？
有没有大哥知道p3 re莉莉姆怎么得到 如题，特殊合成的时候需要莉莉姆，去下层刷了会也没有，看网上要荒魂，但是荒魂也没有。崩溃了
这样只爆了一张，以后拿个保底得了，寻也白寻，没的意思。
1 下一页