行为不当的积分函数

如何解决行为不当的积分函数

我正在尝试在 100x100 点的网格上评估函数 f(x,y)。该函数是一个数值积分。但是，在这些点评估函数，然后绘制 (plt.contourf) 后，我得到了一个奇怪的结果。我怀疑集成可能是原因，但这种行为非常出乎意料，所以我不能确定。积分通过 scipy.integrate.quad(lambda s: f(s,x,y),np.Inf) 进行评估，嵌套 for 循环遍历 x 和 y 坐标（无网格）

蓝色实线也是我希望我的函数评估的部分。显然，它的形状是正确的，但有一些可疑的事情正在发生；我无法解释锯齿行为。这些点在对数空间中均匀分布。

如果有帮助，我可以显示该函数的更多详细信息。任何想法这里发生了什么？还有其他更好的数值积分技术吗？我尝试了 quadpy，但它给出了相同的结果，这可能是意料之中的。也很感激有关如何进一步调查的任何想法

我还沿着 1D 绘制了函数 f(x,y)，并且摆动似乎只沿着一个轴。 IE。

对比

感谢阅读，非常感谢您的帮助。

解决方法

也许只是使用你用数字编码的黎曼和，这样你就可以调整 dx、dy ...

很抱歉，我没有足够的声誉来写这篇评论。