问题描述
如何计算该积分的值:
f_tu(t)给出为numpy.array
。该图如下所示:
我该如何实现? 我能找到的一切看起来都像这样
from scipy.integrate import quad
def f(x):
return 1/sin(x)
I = quad(f,1)
解决方法
sklearn.metrics中的auc怎么样?
import numpy as np
import numpy as np
from scipy.integrate import quad
from sklearn.metrics import auc
x = np.arange(0,100,0.001)
y = np.sin(x)
print('auc:',auc(x,y))
print('quad:',quad(np.sin,100))
auc:0.13818791291277366
方形:(0.1376811277123232,9.459751315610276e-09)
好的,所以您有一个讨厌的无穷大积分。这是我的处理方式:
import numpy as np
from scipy.integrate import quad
def f(x):
return(1/(x**2)) #put your function to integrate here
print(quad(f,np.Infinity)) #integrates from 0 to infinity
这将返回两个值。第一个是积分的估计值,第二个是积分的近似绝对误差,这很有用。
,如果要集成一个numpy数组,这是一个简单的解决方案:
import numpy as np
print(np.trapz(your numpy array here))