pyplot 等高线图的索引错误

Question

我正在尝试使用来自电路的实验数据创建带有 matplotlib.pyplot 的等高线图。我正在使用带有 3 个数据参数的 scipy curve_fit 函数（由大小为 251 的数组给出），但我不断收到索引错误。这是我的情节代码的一部分：

Npts = 10000
v1scan = np.zeros(Npts) #voltage1
v2scan =np.zeros(Npts) #voltage2
tauscan = np.zeros(Npts) #time constant
chi_dof = np.zeros(Npts)
udata = np.zeros(Npts) #uncertainty
i = 0

dof = len(time) - len(popt)

for v1par in np.linspace(0, 8, 100, True):
    for v2par in np.linspace(0, 8, 100, True): 
        for taupar in np.linspace(5*(10**-5), 7*(10**-5), 100, True):
            v1scan[i] = v1par
            v2scan[i] = v2par
            tauscan[i] = taupar
            dymin = (vcap - capvoltage(time, v1par, v2par, taupar))/(0.02*vcap) #vectorized
            chi_dof[i] = sum(dymin*dymin)/dof
            i = i + 1
plt.figure() #creates new figure
ncols = 10
plt.tricontourf(v1par, v2par, taupar, chi_dof, ncols)
plt.colorbar()
plt.show()

我得到的错误信息是：

IndexError                                Traceback (most recent call last)
<ipython-input-175-7562791ba165> in <module>()
     12     for v2par in np.linspace(0, 8, 100, True):
     13         for taupar in np.linspace(5*(10**-5), 7*(10**-5), 100, True):
---> 14             v1scan[i] = v1par
     15             v2scan[i] = v2par
     16             tauscan[i] = taupar

IndexError: index 10000 is out of bounds for axis 0 with size 10000

我该如何解决这个问题？

Answer 1

100 乘以 100 乘以 100 = 1000000 而不是 10000

一旦 i 达到 10000（其最大值应为 9999），您就会尝试访问内存不足区域。

编辑：修复：

Npts = 1000000