我必须在我决定的范围内整合一个信号。
我用过这个方法:
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 是信号的两个值之间的步长。
这是对的吗?
并不真地。要集成采样函数,您可以使用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)