行为不端的积分函数

Posted 2023-03-24

【中文标题】行为不端的积分函数

【英文标题】：Misbehaving Integral Function

【发布时间】：2021-09-20 05:03:52

【问题描述】：

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

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

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

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

对

感谢阅读，非常感谢任何帮助。

【问题讨论】：

你的函数和积分看起来很小。我怀疑您在模型中看到舍入错误。也许重新调整问题，使大多数数量是订单 1？

是的，我认为这是@Nachiket 的问题。感谢您的建议

【参考方案1】：

也许只需使用您以数字方式编码的黎曼和，以便您可以调整 dx、dy ...

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