python 半曲线图中的多项式拟合

Question

我正在尝试获取适合我的数据的多项式。目前，我正在使用 numpy 中的 polyfit 来获得对数对数图中的最佳拟合。但我的目标是让数据适合半逻辑图。我的代码如下所示：

import matplotlib.pyplot as plt
import numpy as np
from scipy.optimize import curve_fit
import scipy.optimize as optimization

l = [ 0.006,   0.01,    0.014,   0.024,   0.0346,  0.049,   0.0535,  0.0736,  0.11  ]
f = [5.3375903383330048, 60.531976422513054, 89.111502526131474, 47.132498501584969, 17.447001214543118, 5.2583622688081455, 3.7779565652126865, 1.0621247249682186, 0.1922152085619766]


logx = np.log(l)
logy = np.log(f)
coeffs = np.polyfit(logx,logy,deg=3)
poly = np.poly1d(coeffs)
yfit = lambda x: np.exp(poly(np.log(x)))
plt.loglog(l,yfit(l), ':')




plt.loglog(l,f, 'o')
plt.show()

如果您建议我必须进行哪些更改才能获得 semilogy 最佳拟合曲线，我将不胜感激。另外，如果python中还有其他包，请一并提及。

Answer 1

我觉得是这样。

# log sacle    
x2 = np.linspace(np.min(l), np.max(l), 1000)
y2log = poly(np.log(x2))
plt.loglog(x2,np.exp(y2log), ':')

plt.loglog(l,f, 'o')
plt.show()

python 半曲线图中的多项式拟合

polynomial fitting in a semilogy plot in python

python

numpy

curve-fitting