scipy的curve_fit没有收到预期的结果

时间:2016-09-29 11:06:17

标签: python scipy curve-fitting

我试图用scipy curve_fit对我的数据进行多元对数回归,结果期望获得一条线,但得到一条曲线。 这是我使用的代码:

Quercetin=[23,195,6,262,272,158,79,65,136,198]
Naringenin=[11,4,8,6,6,7,6,9,7,9]
Rutin=[178,165,93,239,202,3325,4427,7607,3499,1762]
TEAC=[23,189,37,265,290,267,362,388,364,321]

import matplotlib.pyplot as plt
import scipy
from scipy.optimize import curve_fit
import numpy as np
def func(x, a, b, c, d,e):
    m=np.log(a*x[0]+b*x[1]+c*x[2])
    return d*(m)+e
x=scipy.array([Quercetin, Naringenin,Rutin])
y=scipy.array(TEAC)
popt, pcov = curve_fit(func, x ,y)
print (popt)

plt.plot(func(x,*popt),y,'ro-')
plt.show()

我得到了这个结果: enter image description here

虽然我想得到这样的东西:

enter image description here

有谁能请我暗示我做错了什么? 如果这很重要,我在Windows 10上使用Anaconda的Python 3.5。

1 个答案:

答案 0 :(得分:2)

(我假设问题实际上是关于绘图)。 您要求matplotlib绘制由直线('ro')连接的红点(-)。 Matplotlib以它们给出的顺序强制并连接它们。

如果您想绘制一条线,只需将其单独绘制:

In [58]: yres = func(x, *popt)

In [59] plt.plot(yres, y, 'ro')
Out[59]: [<matplotlib.lines.Line2D at 0x7f5c0796b828>]

In [60]: plt.plot([0, 400], [0, 400], '-')
Out[60]: [<matplotlib.lines.Line2D at 0x7f5c07a9c860>]