Scipy 辛普森积分以获取值数组

Question

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 package 将 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 积分和 plot 它与 t 得到一个速度值列表或数组，但我尝试的一种方法是给我一个错误，而另一种方法是给一个奇异值。

Answer 1

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

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

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