Jak wykreślić całkę sygnału w miarę upływu czasu?

Question

Mam sygnał zależny od czasu.

Chcę zobrazować jego integrację w czasie z czasem będącym osią X i wartością integracji będącą osią Y.

Czy jest jakiś sposób na wykonanie Python way?

Aby być bardziej konkretne:

Mam tablicę czasową, time oraz tablicę sygnału, signal. Mają ten sam wymiar.

Muszę zintegrować signal przez time z scipy.integrate.trapz().

Zamiast uzyskać ostatnią całkę, chcę zobaczyć całkę zmieniającą się wraz z upływem czasu.

Answer 1

Spróbuj pomocą scipy.integrate.cumtrapz() zamiast:

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

Oblicza tablicę zawierającą skumulowane wartości całkowite.

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

Answer 2

Nieco lepsza odpowiedź używa opcjonalnego "initial" argument. Powyżej znajduje się pełna przykład:

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() 

ten sposób unika się dziwny indeksowanie jak time[:-1]