我有一个函数,它接受字典并返回加速度
a
。
我将某个参数的数据列表(我们称之为
alpha
)放入该字典中,该函数返回一个包含 100 个结果的数组(相对于 alpha)。
现在我想将这些值一一积分,以获得相对于参数
v
的数组速度 alpha
(以便我可以在 matplotlib 中绘制这些值)。
编辑:
实际上我有机械问题。梯子靠墙放置,由于其
friction coefficient
太小,所以开始滑动。
所以我写了一个函数,它返回任何 acceleration
的梯子 angle between 30° and 90°
(30 是起始角度,90 是梯子在地面上时的角度)。
无论如何,我已经根据角度
acceleration
计算了alpha
。我将数据存储在数组a
中。
我尝试了 scipy.integrate.simps()
但该方法返回最终速度 v
。我也尝试过 scipy.integrate.quad()
但返回相同的结果。
我想生成一个速度数组
v
,其中包含每个角度的值,就像加速度数组一样 - a
scipy.integrate.cumtrapz
应该是你想要的。这将使用梯形规则对提供的数组执行累积积分,产生与输入大小相同的输出数组。如果初始值非零,则可以通过单独的参数控制初始值。