Интеграл данных со временем (Python)

Question

У меня есть временной ряд x(t), который является массивом NumPy. В моем задании мне сообщается, что мне нужно найти интеграл этих данных со временем.

Как я должен это делать? Это не функция, которую мне нужно интегрировать, это список данных.

Answer 1

Это зависит от постановки задачи. Грубый подход будет что-то вроде этого

import numpy as np 
import scipy as sp 

t = np.linspace(-1, 1, 100) 
x = t*t 

delta = t[1] - t[0] 
I = sum(delta*x) 

Answer 2

Вы можете использовать правило Симпсона. Обычная процедура, которая делает это для вас, - simps в spicy.integrate.

>>> help(scipy.integrate.simps) 
Help on function simps in module scipy.integrate.quadrature: 

simps(y, x=None, dx=1, axis=-1, even='avg') 
    Integrate y(x) using samples along the given axis and the composite 
    Simpson's rule. If x is None, spacing of dx is assumed. 

    If there are an even number of samples, N, then there are an odd 
    number of intervals (N-1), but Simpson's rule requires an even number 
    of intervals. The parameter 'even' controls how this is handled.