计算曲线内的面积

Question

我正在测量通过样品的电流，以改变其上的电压。结果是这样的电流-电压图

https://content.sciendo.com/view/journals/joeb/9/1/graphic/j_joeb-2018-0023_fig_004.jpg

我想计算前半周期曲线内的面积（第一象限曲线的一部分）。不确定在 python 中执行此操作的最佳方法是什么。我试图编写一些代码来查找具有相同 x 坐标的对，并从顶部 y 值中减去底部 y 值，然后遍历第一象限中的所有点。

def LobeAreaByPeriod(all_periods_I):
    print('Starting Lobe area by period')
    LobeArea = []
    for period in all_periods_I:
        halfperiod = round(len(period)/2)
        duration = round(halfperiod/2)
        area = 0
        for i in range(duration):
            area += (period[i] - period[halfperiod - i])
        LobeArea.append(area)
    return LobeArea

我不确定这些对会以这种方式直接位于彼此上方，而且我发现很难检查答案是否正确。有关如何执行此操作的任何提示？

Answer 1

要计算曲线内的面积，您必须进行某种近似。根据您对精度的需要，您选择的方法会有所不同。我喜欢梯形规则。 Numpy 有一个非常好的实现。

np.trapz