python - 曲线拟合 zipf 分布 matplotlib python

Question

我尝试将下面的图(红点)与 Python 中的 Zipf 分布 PDF 拟合，F~x^(-a)。我只是选择了a=0.56并绘制了y = x^(-0.56)，我得到了如下所示的曲线。

曲线显然是错误的。我不知道如何进行曲线拟合。

Answer 1

不确定您到底在寻找什么，但如果您想让模型(函数)适合数据，请使用scipy.optimize.curve_fit:

from scipy.optimize import curve_fit
from scipy.special import zetac


def f(x, a):
    return (x**-a)/zetac(a)


result = curve_fit(f, x, y, p0=[0.56])
p = result[0]

print p

如果您不信任标准化，请添加第二个参数 b 并对其进行拟合。