python - 時間の経過とともに信号の積分をプロットする方法は?

Question

時間依存の信号があります。

時間をx軸、積分値をy軸にして、時間の経過に伴う積分をプロットしたいと思います。

これを行うPythonの方法はありますか?

もう少し詳しく言うと：

時間配列timeと信号配列 がありsignalます。それらは同じ次元です。

と統合signalする必要があります。timescipy.integrate.trapz()

最終的な積分を取得する代わりに、時間の経過とともに変化する積分を確認したいと考えています。

Answer 1

少し良い答えは、オプションの"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]