我希望y=mx+c
直线适合我的数据点,但是以日志形式。为此我使用curve_fit
模块。我的简单代码是
def func(x,m,c):
return (x*m + c)
x=log10(xdata)
y=log10(ydata)
err=log10(error)
coeff, var = curve_fit(func,x,y,sigma=err)
yfit = func(x,coeff[0],coeff[1])
pl.plot(x,y,'r0')
pl.plot(x,yfit,'k-')
pl.show()
该图给出了y轴上的负数,因为我的y值是mV。有没有办法在原木拟合的情节上使用原始xdata和ydata(以mV为单位)?
答案 0 :(得分:1)
改为绘制变换后的变量。
plot(10**x, 10**yfit, 'k-')
并且可能以对数刻度显示图表
set_yscale('log')