我想使用scipy.integrate.quad评估函数的积分。这是被积数的样子:
我们可以注意到,该被积物的大部分贡献将来自.1到4或5 ish。当x = 10或更大时,尽管该函数具有振荡性(很难在图片上说出),但它非常小并且会越来越小。
这是从0到某个上限的积分结果。积分的上限在x轴上。
在这里,虽然结果对于前一百个x来说似乎是稳定的,但即使我期望一条直线,情况也不再如此。
我是python的新手,我不知道什么是最好的操作方法。现在,我最好的猜测是将小于100的某个值作为我的积分的上限,并舍弃其他值,因为这只是来自tegral.quad的不良收敛。
编辑: 要绘制第二张图,我使用了scipy.integrate.quad函数。但是,如果仅用生成的点来绘制被积物(第一个图),然后在scipy.integrate.simps中使用该点,并改变最大积分倍数,我将得到一致的结果。