我不明白interp1d报告的结果。我在应该接收号码的地方收到NAN。
In [131]: bb
Out[131]:
array([ 0. , 1.80286595, 1.87443683, 2.70410611, 3.02764722,
3.11305985, 3.11534355, 3.18695351, 3.20693444])
In [132]: alphas1
Out[134]:
array([ 3.80918778e+00, 2.06547222e+00, 1.99234191e+00,
7.55942418e-01, 2.56971574e-01, 1.05144676e-01,
9.30852046e-02, 1.52574183e-02, 1.23664407e-07])
In [135]: bb.shape
Out[135]: (9,)
In [136]: alphas1.shape
Out[140]: (9,)
In [141]: pol = interp1d(alphas1, bb, bounds_error=False)
In [149]: pol(pol.x)
Out[149]: array([ nan, nan, nan, nan, nan, nan, nan, nan, nan]) # I was expecting to receive nan only at the borders.
我认为,如果您检查
interp1d
类的source code,即_check_bounds
方法,就可以看到这个问题:该方法检查您试图输入的x值是否小于
self.x[0]
,这是x
的第一个元素(在您的例子中是alphas1
)。由于alphas1[0]
是x
列表中最大的元素,因此之后的每个元素都将“超出范围”,即小于第一个元素。解决这个问题的一种方法是颠倒您的
x
和y
列表:现在
alphas1
将如scipy所期望的那样增加,并且pol(pol.x)
将如预期那样返回bb
(现在反转)。相关问题 更多 >
编程相关推荐