在下面的代码中,I_sim和t_sim的长度为501。但是,在函数方程中,t具有不同的长度(1001)。当我有两个不同的维度(t_和t_sim)时,如何在下面的代码中使用插值?这是我的错误:ValueError:操作数不能与形状(1001,)(501,)
一起广播'''
conf_sim = copy(conf)
conf_sim['D'] = 0.00042
conf_sim['k0'] = 0.0123
conf_sim['kc'] = 0.0321
V_sim, I_sim, t_sim = sim_multi_cv(conf_sim, p=2, L=1000, phi=-2.5)
t_sim = t_sim[::2]
I_sim = I_sim[::2]
plt.plot(t_sim,I_sim)
plt.show()
from scipy.interpolate import interp1d, interp2d
def fitfunction(t_, I_sim , D=2, k0=0.01,kc=0.01):
conf_residual = copy(conf)
conf_residual['D'] = D
conf_residual['k0'] = k0
conf_residual['kc'] = kc
V,I,t = sim_multi_cv(conf_residual, p=2, L=1000, phi=-2.5)
#interpolation between t_ and t_sim?
??
return I
model = Model(fit_model, independent_vars= ['t_'])
model.set_param_hint('D', value =0.0001)
model.set_param_hint('k0', value =0.01)
model.set_param_hint('kc', value =0.01)
params = model.make_params(D=2, k0=0.01,kc=0.01)
result = model.fit(I_sim, t_ = t_sim)
print(result.values)
result.plot()
plt.show()
''''''
追踪(最近通话):文件“ /data/user/0/ru.iiec.pydroid3/files/accomp_files/iiec_run/iiec_run.py”,第31行,在start(fakepyfile,mainpyfile)文件“ /data/user/0/ru.iiec.pydroid3/files/accomp_files/iiec_run/iiec_run.py”,在开始的第30行exec(open(mainpyfile).read(),main。dict)文件“”,第3行,在NameError:未定义名称“副本”
[程序完成]
缺少代码...