插值与求值问题scipy.interpolate.RectBivariateSplin

2024-09-28 03:20:09 发布

您现在位置:Python中文网/ 问答频道 /正文

我有一个数组,它是(158575),基本上,是8575高斯数叠加在一起(所有这些都略有不同)。我需要把它转换成一个数组,它是(15127181),用来和另一个数组进行卷积。经过一番搜索,我在这里找到了一个与我的问题类似的问题,建议使用scipy.interpolate.RectBivariateSpline公司对数据进行插值,然后将函数计算为更大的数组。在

这对于插值来说似乎很有效(至少我没有收到错误消息,但不知道函数是否正确),但我似乎无法让它在不产生错误的情况下评估函数。我肯定我遗漏了一些明显的东西,但看不到。相关代码如下:

apoLSF_wavelen = np.zeros(len(apoLSF[0]))
for i in range(len(apoLSF[0])):
    apoLSF_wavelen[i] = 10**(apoLSF_l_s + apoLSF_l_d * i)
y = np.arange(15)
f_LSF = interp.RectBivariateSpline(y, apoLSF_wavelen, apoLSF)
siu_LSF = f_LSF.ev(y, siu_wavelen)

我从最后一行得到的错误如下。第一次我以为我把参数的顺序弄错了,然后把它们交换过来,结果发现错误仍然存在,但是数字交换了。在

^{pr2}$

如果上面的代码不足以分析错误,我很抱歉,但我不想复制大量不必要的东西,或者如果有一个明显的解释我忽略了这里或其他地方。在

敬上。在


Tags: 函数代码len错误npscipy数组卷积

热门问题