python - 在 Python 中拟合麦克斯韦-玻尔兹曼分布

Question

是否可以像 matplotlib 中的数据或 python 中的类似模块那样对麦克斯韦-玻尔兹曼进行拟合？

Answer 1

scipy.stats 支持 maxwell distribution .

---

import scipy.stats as stats
import matplotlib.pyplot as plt
import numpy as np

maxwell = stats.maxwell
data = maxwell.rvs(loc=0, scale=5, size=10000)

params = maxwell.fit(data, floc=0)
print(params)
# (0, 4.9808603062591041)

plt.hist(data, bins=20, normed=True)
x = np.linspace(0, 25, 100)
plt.plot(x, maxwell.pdf(x, *params), lw=3)
plt.show()

第一个参数是从零开始的位置或位移。第二个参数是缩放参数，在the wikipedia page上用a表示。 .

要使用此分布生成随机变量(随机数据)，请使用其 rvs 方法:

newdata = maxwell.rvs(*params, size=100)