python - 用 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