如何在Python中集成信号

Question

我必须在我决定的范围内整合信号。

我使用过这种方法：

def integrand(x,y):
    return y[x]

result=scipy.integrate.quad(integrand,t0,t1,args=(y))*0.2E-6

y是500microsec中的2500个点的数组。例如，t0是700，t1是1500。 0.2E-6是信号的两个值之间的步骤。

这是对的吗？

Answer 1

不是真的。要集成采样功能，您可以使用scipy.integrate.simps或scipy.integrate.trapz：

result_simps = scipy.integrate.simps(y[t0:t1+1], dx=0.2E-6)
result_trapz = scipy.integrate.trapz(y[t0:t1+1], dx=0.2E-6)