Как построить интеграл сигнала с течением времени?

Question

У меня есть зависящий от времени сигнал.

Я хочу построить интеграцию с течением времени, когда ось x и значение интегрирования являются осью y.

Есть ли Python way?

Чтобы быть более конкретным:

У меня есть время массив, time и массив сигнала, signal. Они имеют одинаковую размерность.

Мне нужно интегрировать signal над time с scipy.integrate.trapz().

Вместо того, чтобы получить окончательный интеграл, я хочу видеть, как интеграл меняется со временем.

Answer 1

Попробуйте использовать scipy.integrate.cumtrapz() вместо:

plt.plot(time[:-1], scipy.integrate.cumtrapz(signal, x=time)) 
plt.show() 

Он вычисляет массив, содержащий кумулятивные целые значения.

http://docs.scipy.org/doc/scipy-0.10.1/reference/generated/scipy.integrate.trapz.html

Answer 2

Немного лучше ответ использует дополнительный аргумент "initial". Вот полный пример:

import scipy.integrate as it 
import numpy as np 
import matplotlib.pyplot as plt 
t=np.linspace(0,1, 100) 
y=t**2 
y_int = it.cumtrapz( y , t, initial=0.0) # y_int is same size as t 
plt.plot(t, y_int) 
plt.show() 

Это позволяет избежать странной индексации как time[:-1]