0
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 进行对比,但我尝试的一种方法是给我一个错误,而另一种方法是给一个奇异值。

4

1 回答 1

0

问:您为什么为此使用分析方法/库?在这种情况下

v = v0 + Sigma( da * dt ) * 常数。

常数取决于产生信号的加速度传感器。如果该传感器线性工作,那么您可以对您的 lambda变化(da) 进行采样,直到 t 具有足够小的步长 (dt),然后将结果乘以步长并累加结果。每个样本都将获得一个您可以绘制的 v 值。要找到 v 的值,请确保添加 v0(起始速度)

于 2021-04-07T14:26:09.943 回答