晚上好!
考虑下图和显示实验数据的图表:
这是一个有两个参数的图; x(在 x 轴上)和 y(我选择值 1 到 5 示例性地保持可读性,从左到右)。值 z - 显示在 y 轴上 - 应该类似于函数
1-exp(-(x^a/y^b))
的输出,其中 a
和 b
是未知参数。
我想找到那些参数的值。我会使用 Scipy 的
curve_fit
,例如与:
import numpy as np
from scipy.optimize import curve_fit
data = np.loadtxt('data.csv',delimiter=',')
x = data[:,0]
y = np.linspace(1,10,10)
z = data[???]
#Fitfunktion
def func(X, a, b):
x,y = X
return 1-np.exp(-(np.power(x, a)/np.power(y, b)))
print(curve_fit(func, (x,y), z))
但是正如您所看到的,我不知道如何将实验中的值获取到我的代码中。 数据保存在 CSV 文件中(第一列是 x 值,第一行是 y 值),但我不知道如何在我的代码中使用它。
如果有任何建议,我将不胜感激。
致以诚挚的问候, 吉斯康纳