用Python求解幂律分布

Question

我的数据与power law发行版非常相似。使用Python，我希望通过以以下形式求解两个方程来近似数据：

​

​

y是y轴数据。在Python中，它将是data[i]。X为i + 1。因此，在第一个数据索引处和在数据中的“随机”第二个变量处，我们得到了两个含有两个未知变量的方程：

​

​

问题归结为解决问题

​

​

因为数学上的简化。我不知道如何使用像numpy.linalg.solve这样的库来求解这样的方程。如何使用Python找到a的值？

Answer 1

好吧，我知道了。

import math

def get_power_law_variables(data):
    c = data[0]
    middle_index = len(data) / 2
    division = float(data[middle_index]) / c
    logarithm_base = middle_index + 1
    a = math.log(division, logarithm_base)
    return c, a

# Example usage
data = range(50, 150)
c, a = get_power_law_variables(data)
print c, a