import numpy as np
from scipy.integrate import simpson
我正在尝试整合下面给出的加速度信号:
f2 = lambda t: 0.1*np.sin(np.pi*t) + 0.1*np.sin(20*np.pi*t) + 0.01
对于时间点:
t = np.arange(0, 20, 0.1)
我正在尝试使用 scipy 中的 simps/simpson 包将 f2 集成到 t 上,但出现此错误:
v2 = simpson(f2, t)
IndexError: tuple index out of range
如果我使用如下列表理解创建加速数据点列表:
a2 = [f2(i) for i in t]
并使用辛普森规则整合上面的列表,我只得到一个值:
v2 = simps(a2, t)
这是
v2 = 0.2005702922522627
我想在 t 期间使用 Simpson 积分获取速度值的列表或数组,并将其与 t 进行对比,但我尝试的一种方法是给我一个错误,而另一种方法是给一个奇异值。