如何在Python中集成信号

Question

我必须在我决定的范围内积分信号。

我用过这个方法：

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

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

y 是 500 微秒内 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)