Wie zeichnet man das Integral eines Signals mit der Zeit auf?

Question

Ich habe ein zeitabhängiges Signal.

Ich möchte seine Integration über die Zeit mit Zeit als x-Achse und Integrationswert als y-Achse darstellen.

Gibt es Python Weg dies zu tun?

Um genauer zu sein:

Ich habe eine Zeit Array, time, und ein Signal Array, signal. Sie sind von gleicher Dimension.

Ich muss signal über time mit scipy.integrate.trapz() integrieren.

Anstatt das letzte Integral zu erhalten, möchte ich sehen, dass sich das Integral mit der Zeit ändert.

Answer 1

Versuchen Sie es mit scipy.integrate.cumtrapz() statt:

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

Es berechnet ein Array mit den kumulativen Integralwerte enthält.

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

Answer 2

Eine etwas bessere Antwort wird das optionale "initial" Argument. Hier ist ein komplettes Beispiel:

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

Dies vermeidet die seltsame Indizierung wie time[:-1]