Question

我正在尝试使用Python插值函数来获取给定x的值y但是我得到错误“提高ValueError（”x和y数组的长度必须沿着插值轴“，即使我的数组也是如此具有相同的大小和形状（根据我在代码中使用.shape得到的结果）。我对编程很新，所以我不知道如何检查我的数组中还有什么不同。这是我的代码：

s = []

def slowroll(y, t):

    phi, dphi, a = y
    h = np.sqrt(1/3. * (1/2. * dphi**2 + 1/2.*phi**2))
    da = h*a
    ddphi = -3.*h*dphi - phi

    return [dphi,ddphi,da]


phi_ini = 18.
dphi_ini = -0.1
init_y = [phi_ini,dphi_ini,1.]
h_ini =np.sqrt(1/3. * (1/2. * dphi_ini**2. + 1/2.*phi_ini**2.))

t=np.linspace(0.,20.,100.)

from scipy.integrate import odeint

sol = odeint(slowroll, init_y, t)

phi = sol[:,0]
dphi = sol[:,1]
a=sol[:,2]
n=np.log(a)
h = np.sqrt(1/3. * (1/2. * dphi**2 + 1/2.*phi**2))

s.extend(a*h)
x = np.asarray(s)
y = np.asarray(t)
F = interp1d(y, x, kind='cubic') 

print F(7.34858263)

Answer 1

添加所需的导入后，我无法使用2.7.12版复制您的错误。你在运行什么python版本？

import numpy as np
from scipy.interpolate import interp1d

s = []

def slowroll(y, t):

    phi, dphi, a = y
    h = np.sqrt(1/3. * (1/2. * dphi**2 + 1/2.*phi**2))
    da = h*a
    ddphi = -3.*h*dphi - phi

    return [dphi,ddphi,da]


phi_ini = 18.
dphi_ini = -0.1
init_y = [phi_ini,dphi_ini,1.]
h_ini =np.sqrt(1/3. * (1/2. * dphi_ini**2. + 1/2.*phi_ini**2.))

t=np.linspace(0.,20.,100.)

from scipy.integrate import odeint

sol = odeint(slowroll, init_y, t)

phi = sol[:,0]
dphi = sol[:,1]
a=sol[:,2]
n=np.log(a)
h = np.sqrt(1/3. * (1/2. * dphi**2 + 1/2.*phi**2))

s.extend(a*h)
x = np.asarray(s)
y = np.asarray(t)
F = interp1d(y, x, kind='cubic') 

print F(7.34858263)

输出： 2.11688518961e+20

python - 沿插值轴的不同数组长度？

1 个答案: