Làm thế nào để vẽ tích phân của một tín hiệu khi thời gian trôi qua?

Question

Tôi có tín hiệu phụ thuộc vào thời gian.

Tôi muốn âm mưu tích hợp theo thời gian với thời gian là trục x và giá trị tích hợp là trục y.

Có cách nào Cách Python làm việc này không?

Để cụ thể hơn:

Tôi có một mảng thời gian, time, và một mảng tín hiệu, signal. Chúng có cùng kích thước.

Tôi cần tích hợp signal qua time với scipy.integrate.trapz().

Thay vì nhận được tích phân cuối cùng, tôi muốn thấy sự khác biệt tích phân như thời gian trôi qua.

Answer 1

Hãy thử sử dụng scipy.integrate.cumtrapz() thay vì:

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

Nó tính toán một mảng chứa các giá trị không thể thiếu tích lũy.

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

Answer 2

Một câu trả lời tốt hơn một chút sử dụng tùy chọn "initial" tranh cãi. Dưới đây là một ví dụ hoàn chỉnh:

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

Điều này tránh việc lập chỉ mục lạ như time[:-1]